Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  495 1733 (N56E46 L= 80)   745     59        0 ( 472/1703- 524/1763) -22.4 E55-E38 N52-N59
      2 1347 1613 (N32W32 L=157)   362     84        1 (1307/1585-1387/1635) -36.9 W29-W34 N30-N34
      3  227 1409 (N33E63 L= 63)   984    113        2 ( 193/1347- 259/1452) -25.2 E73-E55 N28-N36
      4  660 1322 (N22E20 L=106)  1309    430        4 ( 463/1272- 825/1390)   2.2 E34-E08 N18-N23
      5 1861  942 (S15W61 L=186)  1875    291        0 (1839/ 806-1885/1085)  86.4 W58-W64 S22-S05
      6  302  715 (S10E54 L= 71)  1022    170        0 ( 260/ 652- 327/ 799)  45.0 E61-E51 S13-S05
      7  568  635 (S18E35 L= 91)   225     50        0 ( 550/ 607- 585/ 661) -90.0 E35-E33 S20-S16
      8  370  533 (S22E56 L= 70)  1117    207        1 ( 307/ 470- 436/ 600)  15.3 E68-E46 S25-S18
      9  533  550 (S23E40 L= 85)   354    124        0 ( 488/ 523- 577/ 592) -45.0 E42-E37 S25-S20
     10  729  348 (S40E35 L= 91)   401     65        1 ( 705/ 326- 770/ 368)  15.9 E37-E30 S41-S38
     11  615  273 (S44E52 L= 74)   461     62        0 ( 583/ 254- 645/ 290)   4.8 E58-E46 S45-S42
  

    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