Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  716 1642 (N39E21 L=217)   552     97        0 ( 666/1627- 776/1654)   0.0 E25-E16 N38-N40
      2  188 1079 (N05E61 L=177)  1257    135        1 ( 146/1008- 227/1131) -35.0 E66-E56 N00-N08
      3  238  758 (S15E60 L=178)   553     92        0 ( 221/ 715- 265/ 799) -90.0 E61-E57 S17-S11
      4 1083  750 (S19W03 L=240)   258     65        0 (1062/ 714-1111/ 780)  53.1 W01-W04 S21-S17
      5 1499  689 (S25W31 L=269)   579    183        0 (1431/ 637-1563/ 760)  33.7 W26-W35 S27-S20
      6  224  687 (S19E65 L=173)  1385    166        4 ( 190/ 627- 260/ 754) -76.0 E69-E60 S22-S14
      7  553  506 (S33E39 L=199)   674    191        0 ( 536/ 415- 573/ 598)  52.4 E44-E34 S39-S26
      8  369  465 (S34E60 L=177)   664     95        2 ( 345/ 420- 393/ 509)  19.7 E68-E54 S36-S31
      9 1129  315 (S51W06 L=244)   267     32        0 (1115/ 292-1143/ 341) -76.0 W04-W07 S52-S48
     10  876  267 (S54E20 L=218)   336     44        0 ( 861/ 238- 891/ 295)  59.0 E22-E18 S56-S51
  

    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