Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  749 1712 (N59E09 L=359)  1562    207        4 ( 681/1643- 831/1752)  10.6 E17-W01 N53-N62
      2 1307 1696 (N43W42 L= 50)  2097    435        2 (1131/1656-1502/1728) -20.1 W26-W58 N35-N48
      3 1529 1504 (N25W47 L= 55)   416    184        0 (1481/1470-1566/1533) -26.6 W42-W49 N23-N27
      4 1336 1266 (N15W24 L= 32)  1229    266        5 (1261/1156-1376/1357) -85.9 W20-W25 N08-N22
      5  610 1337 (N32E23 L=345)   311     89        1 ( 573/1312- 655/1356)   0.0 E26-E19 N31-N33
      6  765 1223 (N23E13 L=355)   283    113        1 ( 751/1159- 780/1268)  81.9 E14-E11 N18-N25
      7  595 1081 (N17E26 L=342)   508    161        0 ( 527/1047- 661/1113)  11.3 E31-E21 N15-N17
      8 1757  714 (S28W49 L= 57)  2402    219        6 (1712/ 614-1805/ 802) -66.8 W43-W54 S35-S21
      9  939  220 (S48E29 L=339)   531     30        1 ( 915/ 202- 969/ 242) -26.6 E31-E25 S49-S45
  

    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