Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1132 1597 (N42W14 L=136)   262     51        2 (1110/1574-1149/1624)  53.1 W11-W15 N39-N43
      2  719 1566 (N42E21 L=102)   647    131        0 ( 663/1546- 797/1582)   5.2 E25-E14 N40-N43
      3  190 1235 (N20E64 L= 58)  1441    335        3 ( 131/1186- 247/1273) -14.0 E72-E56 N17-N22
      4 1888 1225 (N08W71 L=194)   813    115        1 (1862/1167-1904/1264)  80.5 W67-W73 N04-N10
      5 1751 1137 (N05W52 L=174)   584     82        2 (1730/1094-1769/1184)  71.6 W50-W54 N02-N08
      6 1672 1098 (N04W44 L=167)   600    145        3 (1639/1047-1699/1166)  45.0 W41-W47 N01-N07
      7 1534  919 (S05W32 L=155)   233     52        1 (1516/ 891-1549/ 950)  56.3 W30-W33 S07-S03
      8 1617  686 (S21W39 L=162)   630    151        1 (1564/ 656-1683/ 716)   6.3 W34-W44 S22-S20
      9  772  638 (S18E19 L=103)   277     97        0 ( 723/ 605- 821/ 665)  23.2 E22-E15 S19-S16
  

    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