Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  496 1578 (N20E47 L= 46)   923    174        2 ( 451/1506- 529/1661) -76.0 E50-E44 N14-N26
      2  629 1601 (N23E37 L= 56)   402    120        1 ( 604/1557- 655/1654) -82.9 E38-E35 N19-N27
      3  338 1606 (N21E66 L= 27)   893    100        2 ( 318/1558- 363/1651) -37.9 E72-E61 N17-N24
      4  977 1579 (N27E13 L= 80)   597    151        2 ( 912/1521-1031/1623) -26.6 E18-E08 N23-N29
      5  368 1517 (N15E56 L= 37)   811    114        2 ( 319/1493- 429/1544)  -4.4 E62-E49 N13-N16
      6  315 1442 (N10E57 L= 35)   368     69        1 ( 297/1407- 329/1483) -51.3 E60-E54 N07-N12
      7  333 1143 (S08E48 L= 45)  1368    259        1 ( 282/1043- 394/1251)  65.2 E51-E43 S14-S01
      8 1683 1155 (N13W40 L=133)   301     22        0 (1660/1134-1704/1177)   0.0 W38-W42 N12-N14
      9  400 1087 (S11E41 L= 52)   354    102        0 ( 374/1055- 425/1120) -36.9 E43-E39 S12-S09
     10  334 1045 (S14E46 L= 47)   418     71        1 ( 310/1015- 355/1083)  68.2 E48-E44 S16-S11
     11 1832  940 (N04W57 L=150)  1136    214        3 (1803/ 831-1849/1035) -74.7 W53-W60 S01-N10
     12  492  972 (S17E33 L= 60)   241     42        0 ( 469/ 954- 515/ 997) -33.7 E34-E31 S17-S15
     13  169  906 (S24E65 L= 28)  1129    185        5 ( 128/ 812- 210/ 987) -61.4 E70-E60 S29-S18
     14  571  811 (S26E26 L= 67)  1091    366        1 ( 532/ 647- 637/ 976)  79.2 E29-E22 S36-S15
     15  394  841 (S26E41 L= 52)   802    178        1 ( 356/ 760- 425/ 927)  70.0 E43-E37 S31-S20
     16  716  676 (S32E15 L= 78)   264     86        0 ( 672/ 658- 754/ 703)  38.7 E17-E12 S33-S29
     17  849  596 (S35E04 L= 89)   256     83        0 ( 831/ 552- 872/ 633)  80.5 E05-E02 S38-S32
     18 1353  402 (S38W39 L=132)   266     47        1 (1329/ 373-1370/ 430) -63.4 W37-W41 S39-S35
  

    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