Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  726 1833 (N35E51 L=106)  1749    332        4 ( 687/1773- 759/1907) -42.7 E59-E46 N30-N42
      2  578 1714 (N23E53 L=105)   481     96        0 ( 540/1667- 614/1746)  90.0 E54-E51 N18-N25
      3  350 1607 (N11E68 L= 90)  3874    501        7 ( 267/1536- 444/1672)  38.2 E76-E60 N06-N17
      4  926 1537 (N20E21 L=137)   605    109        1 ( 878/1490- 973/1567)  -7.1 E24-E16 N18-N22
      5  291  783 (S39E46 L=112)   324     54        0 ( 270/ 763- 311/ 805)  36.9 E47-E43 S40-S37
      6 1826  787 (N06W61 L=219)   452    157        0 (1781/ 771-1863/ 797)  15.3 W55-W66 N04-N07
      7 1230  729 (S16W20 L=178)   663    185        0 (1137/ 684-1327/ 788)  55.0 W16-W24 S20-S10
      8 1832  723 (N03W65 L=223)   456     89        0 (1798/ 712-1860/ 740)  20.6 W61-W69 N01-N04
      9 1324  571 (S22W32 L=190)   237     84        0 (1302/ 549-1352/ 611) -11.3 W29-W34 S22-S19
     10  856  447 (S45W09 L=167)  1050    248        3 ( 729/ 412- 978/ 486)  17.5 E01-W18 S48-S41
     11 1201  361 (S38W36 L=194)   372     89        0 (1152/ 345-1250/ 383)  63.4 W33-W38 S40-S34
     12 1473  295 (S31W59 L=216)   865    370        0 (1410/ 259-1558/ 338)  10.0 W51-W68 S32-S28
  

    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