Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1140 1813 (N48W44 L= 31)   384     54        0 (1125/1778-1154/1844)  32.0 W40-W48 N45-N50
      2 1536 1686 (N26W61 L= 48)   792     89        2 (1477/1665-1579/1704) -26.6 W55-W65 N24-N28
      3 1557 1504 (N16W49 L= 36)  1847    532        3 (1321/1411-1744/1567) -41.2 W35-W62 N05-N26
      4 1108 1388 (N23W15 L=  2)  1222    427        5 ( 929/1320-1299/1451) -46.6 W06-W23 N13-N31
      5 1520 1307 (N07W37 L= 24)   243     47        0 (1501/1285-1539/1331)  71.6 W35-W38 N05-N08
      6 1599 1267 (N02W41 L= 28)   508    114        2 (1566/1221-1635/1312)  20.6 W37-W45 N00-N05
      7  845 1209 (N20E05 L=342)   226     66        0 ( 812/1190- 877/1225) -14.0 E07-E03 N18-N20
      8 1924 1146 (S15W71 L= 58)   928     58        1 (1910/1116-1943/1178) -78.7 W68-W76 S17-S12
      9  152  998 (N23E67 L=280)  1281    115        2 ( 104/ 968- 216/1017)  -6.0 E76-E57 N21-N24
     10 1355  985 (S06W17 L=  5)   237     73        0 (1327/ 958-1377/1015)  26.6 W15-W19 S07-S04
     11 1277  833 (S12W09 L=357)   208    109        2 (1252/ 804-1300/ 872)  53.1 W07-W10 S13-S09
     12 1853  841 (S30W55 L= 43)   410     36        1 (1837/ 815-1872/ 865) -71.6 W53-W58 S31-S28
     13 1839  738 (S37W54 L= 42)   558     76        2 (1809/ 700-1868/ 771) -37.9 W49-W61 S40-S33
     14 1462  609 (S31W16 L=  3)  2042    300        6 (1404/ 486-1512/ 730)  51.8 W09-W21 S39-S23
     15 1702  585 (S42W38 L= 25)  2896    339        6 (1555/ 532-1839/ 638) -15.9 W22-W66 S51-S37
     16  907  317 (S34E31 L=316)   417     87        0 ( 892/ 279- 928/ 357)  56.3 E33-E27 S37-S31
  

    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