Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  423 1542 (N30E54 L=145)  1575    161        3 ( 357/1495- 499/1586) -10.3 E66-E44 N27-N32
      2  300 1029 (S03E51 L=148)   291     51        1 ( 283/1004- 325/1050)  33.7 E52-E48 S05-S02
      3  273  885 (S13E54 L=145)   592     76        0 ( 243/ 846- 296/ 927)  50.2 E57-E50 S15-S09
      4  210  837 (S16E62 L=137)   516     73        0 ( 194/ 805- 227/ 873) -63.4 E63-E59 S18-S14
      5  468  752 (S19E36 L=163)   483    133        1 ( 408/ 718- 502/ 796)  26.6 E41-E33 S21-S16
      6  771  768 (S16E14 L=185)   217     24        0 ( 749/ 752- 797/ 785)  18.4 E15-E12 S16-S14
      7 1205  741 (S14W14 L=212)   640    222        0 (1099/ 711-1297/ 774)  -4.4 W06-W19 S15-S12
      8 1677  629 (S19W52 L=251)  1687    289        6 (1593/ 561-1747/ 679)  30.3 W45-W59 S23-S15
      9  777  374 (S44E15 L=184)   451     78        0 ( 746/ 349- 807/ 410)  39.8 E18-E12 S46-S41
     10  534  345 (S50E44 L=155)   532     89        0 ( 491/ 324- 569/ 366) -14.0 E48-E39 S51-S48
  

    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