Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1928 (N60W58 L=263)   823     76        1 (1107/1900-1166/1957)   6.3 W47-W75 N57-N61
      2  924 1842 (N60W21 L=226)  1340    233        2 ( 838/1781-1027/1911)   6.9 W09-W42 N55-N63
      3 1682 1066 (S09W42 L=247)   662    132        2 (1648/1002-1716/1121)  39.8 W38-W45 S11-S06
      4 1793  693 (S35W52 L=256)   808    139        0 (1757/ 660-1848/ 745) -23.2 W46-W60 S37-S30
      5 1708  653 (S35W42 L=246)   309     32        1 (1692/ 630-1725/ 676)  45.0 W40-W43 S36-S33
      6 1001  648 (S19E09 L=195)   241     51        0 ( 976/ 631-1035/ 668) -18.4 E10-E07 S20-S17
      7  901  567 (S22E18 L=187)   253     95        0 ( 863/ 538- 950/ 583)   0.0 E20-E14 S22-S20
      8  710  308 (S35E43 L=162)   921    212        0 ( 618/ 288- 813/ 343) -30.3 E48-E35 S38-S30
      9 1414  284 (S58W14 L=218)   535     79        2 (1365/ 268-1458/ 305) -10.3 W08-W19 S59-S55
  

    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