Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  312 1534 (N06E67 L=246)  1803    163        2 ( 262/1468- 362/1614)  79.7 E71-E64 N01-N12
      2  768 1565 (N18E31 L=281)   249     81        0 ( 738/1541- 800/1597)  -9.5 E34-E28 N17-N19
      3  932 1421 (N14E16 L=296)   461    135        0 ( 891/1363- 960/1479) -82.9 E18-E15 N09-N17
      4 1518 1084 (N12W27 L=339)   317     97        1 (1491/1037-1552/1131) -26.6 W23-W29 N10-N13
      5 1637 1102 (N17W34 L=347)   357     93        0 (1612/1072-1679/1126)  51.3 W32-W37 N14-N19
      6 1674  898 (N07W43 L=355)  1844    477        4 (1585/ 794-1720/ 979) -19.7 W34-W49 N02-N10
      7 1833  841 (N09W59 L= 12)   541    115        0 (1811/ 801-1854/ 886) -33.7 W55-W63 N07-N11
      8  544  713 (S36E22 L=290)   575    161        0 ( 487/ 679- 623/ 746)  37.9 E26-E16 S39-S32
      9 1672  440 (S17W65 L= 18)  1082    235        1 (1648/ 391-1690/ 506) -24.4 W60-W71 S18-S13
  

    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