Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1397 (N15E24 L=285)  2851    918       10 ( 587/1066- 971/1693) -34.6 E48-E03 N00-N31
      2 1115 1652 (N40E09 L=301)   300     64        1 (1083/1632-1151/1670)   0.0 E11-E05 N38-N41
      3 1543 1227 (N21W29 L=339)   267     95        0 (1498/1184-1580/1261) -23.2 W25-W32 N19-N22
      4 1700 1091 (N16W43 L=353)   550    143        1 (1653/1032-1746/1153) -26.6 W38-W48 N13-N19
      5  545  935 (S16E28 L=281)  1582    498        5 ( 357/ 884- 738/1008)  24.8 E42-E16 S21-S08
      6 1727  778 (S02W53 L=  2)   477    129        1 (1679/ 751-1774/ 816)  35.5 W49-W56 S04-N01
      7 1765  600 (S11W65 L= 14)  1018    237        2 (1715/ 567-1808/ 622)   0.0 W58-W73 S12-S09
      8 1045  479 (S35W14 L=324)  2012    442        9 ( 821/ 440-1247/ 551)  27.4 E01-W26 S42-S25
      9  715  469 (S43E10 L=300)   266     53        1 ( 683/ 455- 743/ 488)  11.3 E12-E07 S44-S41
  

    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