Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  986 1880 (N60E36 L=118)  1026     92        2 ( 941/1847-1030/1919) -11.3 E50-E27 N56-N64
      2  675 1571 (N27E36 L=118)   478    183        0 ( 631/1544- 738/1625)   0.0 E41-E31 N25-N30
      3  335 1546 (N20E67 L= 87)  1114    116        2 ( 283/1520- 389/1567)   6.7 E76-E58 N17-N20
      4  578 1460 (N18E38 L=115)   494     55        2 ( 554/1427- 604/1492) -31.0 E40-E35 N16-N19
      5 1592 1219 (N21W35 L=188)   531    146        0 (1575/1140-1606/1286) -77.5 W32-W35 N16-N25
      6  253  836 (S25E55 L= 98)   829    140        0 ( 234/ 767- 283/ 894)  52.1 E58-E51 S29-S20
      7  399  793 (S25E40 L=113)   330     89        0 ( 366/ 778- 446/ 818)   0.0 E42-E35 S26-S24
      8 1105  621 (S24W13 L=166)  2323    602        8 ( 866/ 502-1313/ 750)  46.3 W00-W22 S35-S12
      9 1559  605 (S17W45 L=198)   544    120        1 (1519/ 565-1606/ 637)  51.3 W42-W48 S19-S14
     10  618  246 (S68E31 L=122)   593     89        0 ( 575/ 236- 663/ 258)  13.5 E46-E20 S71-S64
  

    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