Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1377 1434 (N30W24 L=179)  1137    364        2 (1255/1319-1472/1556) -47.0 W15-W30 N22-N38
      2 1611 1517 (N37W48 L=203)  1396    189        1 (1532/1477-1678/1545)   0.0 W38-W55 N34-N39
      3  386 1440 (N25E52 L=104)   441     63        2 ( 357/1416- 415/1461)  45.0 E55-E49 N23-N26
      4 1817 1328 (N24W65 L=221)  1363    217        0 (1731/1257-1876/1391) -28.1 W55-W73 N19-N28
      5 1749  856 (S05W53 L=208)   527     33        1 (1725/ 832-1774/ 883) -11.3 W50-W55 S06-S03
      6  344  640 (S26E51 L=105)  2636    411        2 ( 257/ 492- 434/ 760) -63.4 E62-E45 S36-S18
      7 1675  681 (S17W49 L=205)   404     85        0 (1658/ 644-1692/ 724)  90.0 W48-W50 S19-S13
      8 1788  611 (S21W66 L=222)  1350    180        2 (1753/ 543-1817/ 695) -65.6 W62-W73 S26-S15
      9  614  507 (S34E29 L=127)   480     94        0 ( 573/ 491- 662/ 525) -15.9 E31-E24 S34-S32
  

    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