Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1606 1514 (N18W54 L=234)  1967    430        0 (1432/1385-1764/1604) -55.7 W46-W62 N05-N27
      2 1583 1394 (N12W46 L=226)   309     89        0 (1534/1383-1625/1403) -23.2 W42-W49 N10-N13
      3 1152 1385 (N24W18 L=198)   224     73        0 (1117/1371-1193/1399) -36.9 W15-W20 N22-N25
      4  211 1014 (N23E60 L=120)  1244    236        1 ( 173/ 939- 270/1123) -56.3 E65-E53 N18-N30
      5 1483  927 (S10W24 L=204)   439    149        0 (1454/ 854-1509/ 990) -90.0 W23-W25 S15-S05
      6  166  747 (N07E75 L=105)  1010     98        2 ( 142/ 696- 195/ 785) -59.0 E79-E70 N04-N09
      7 1128  617 (S19E04 L=176)   334    109        0 (1082/ 568-1158/ 663)  29.7 E08-E01 S20-S16
      8 1781  611 (S42W49 L=229)  1220    120        1 (1740/ 568-1827/ 635)   0.0 W43-W57 S45-S38
      9 1675  533 (S44W36 L=216)   406     37        0 (1655/ 515-1700/ 554) -33.7 W33-W39 S45-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