Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  499 1538 (N22E42 L=211)  1555    693        1 ( 350/1445- 699/1655) -18.0 E61-E24 N17-N30
      2  420 1137 (S03E40 L=213)   438    174        1 ( 362/1081- 482/1189) -26.6 E44-E34 S06-N00
      3  201  986 (S13E58 L=194)   523    110        3 ( 166/ 958- 235/1016) -14.0 E62-E54 S14-S10
      4  390  950 (S14E41 L=212)   484    180        4 ( 330/ 916- 458/ 982)  11.3 E45-E35 S16-S12
      5  211  919 (S17E57 L=195)   364    109        3 ( 186/ 900- 246/ 939)  18.4 E59-E53 S17-S15
      6  972  851 (S16E02 L=251)   982    315        6 ( 887/ 737-1095/ 938)  57.3 E06-W05 S23-S09
      7  676  882 (S17E20 L=232)   199     65        1 ( 660/ 849- 697/ 914)  63.4 E21-E19 S18-S14
      8  466  729 (S27E36 L=217)   418    151        0 ( 423/ 675- 518/ 775)  41.2 E39-E31 S31-S24
      9 1249  721 (S21W17 L=270)   237     87        0 (1219/ 683-1269/ 772)  80.5 W15-W18 S23-S17
     10 1623  601 (S23W47 L=299)  1182    176        3 (1554/ 552-1678/ 669)  53.1 W42-W50 S27-S18
  

    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