Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1515 (N31E58 L= 37)   557     71        1 ( 331/1497- 408/1529)   5.7 E64-E53 N29-N31
      2  701 1378 (N24E24 L= 72)   850    304        3 ( 573/1350- 850/1407)  -2.7 E33-E12 N22-N25
      3  433 1238 (N14E42 L= 54)   512    134        0 ( 404/1175- 474/1289)  66.8 E43-E38 N09-N16
      4  829 1251 (N16E13 L= 82)   198     57        0 ( 791/1235- 859/1264) -14.0 E15-E11 N15-N16
      5 1239 1197 (N14W13 L=108)   251     86        0 (1197/1169-1276/1223) -21.8 W10-W15 N12-N15
      6 1805  917 (S03W58 L=153)   465     77        1 (1762/ 887-1840/ 943)  20.6 W53-W61 S04-S01
      7 1452  761 (S12W29 L=124)   267    127        0 (1404/ 725-1494/ 810)  39.8 W25-W31 S14-S09
      8 1778  719 (S15W59 L=154)   538     72        1 (1756/ 679-1798/ 760)  78.7 W56-W60 S17-S12
      9  863  476 (S34E10 L= 85)  1044    283        4 ( 806/ 356- 939/ 581) -70.6 E13-E04 S43-S26
     10 1452  528 (S28W33 L=129)   463    112        1 (1409/ 497-1492/ 562)  45.0 W30-W35 S30-S25
     11 1456  274 (S50W51 L=146)   593     72        3 (1425/ 247-1484/ 303) -16.4 W44-W61 S52-S46
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y