Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  375 1086 (N24E42 L=317)   564    148        1 ( 314/1051- 439/1116) -26.6 E46-E36 N21-N26
      2  310  969 (N18E49 L=310)   614    130        0 ( 263/ 930- 361/1011) -40.6 E52-E44 N15-N21
      3 1379  915 (S08W18 L= 17)   366     84        1 (1334/ 883-1425/ 936)   0.0 W14-W21 S08-S06
      4  663  830 (N04E26 L=333)   425     82        2 ( 636/ 790- 693/ 875) -68.2 E26-E24 N01-N06
      5 1198  753 (S12W04 L=  2)   649    216        1 (1161/ 656-1227/ 861)  61.4 W00-W06 S17-S06
      6  661  729 (S02E28 L=330)   220     31        0 ( 639/ 710- 683/ 753) -90.0 E29-E26 S03-N00
      7 1193  521 (S27E02 L=357)   713    252        0 (1076/ 435-1301/ 616)  18.4 E12-W06 S30-S23
      8  843  435 (S24E27 L=332)   257     71        1 ( 815/ 409- 872/ 466)  26.6 E29-E23 S25-S22
      9  837  238 (S39E38 L=320)   777    248        1 ( 812/ 172- 867/ 299)  40.2 E45-E32 S45-S34
  

    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