Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1600 (N18E40 L=184)   426    101        1 ( 618/1561- 664/1663) -40.6 E44-E36 N15-N21
      2  353 1554 (N10E61 L=162)   758     76        2 ( 322/1517- 379/1593) -24.0 E66-E57 N08-N12
      3 1640 1526 (N41W31 L=255)   340     29        0 (1624/1508-1660/1554) -26.6 W29-W33 N39-N42
      4  529 1482 (N09E42 L=182)   571    152        1 ( 497/1413- 575/1553) -35.0 E46-E36 N05-N12
      5  293 1500 (N06E65 L=159)   512     50        1 ( 271/1475- 312/1524)  90.0 E67-E62 N04-N07
      6  407 1259 (S07E44 L=180)   317     79        0 ( 365/1236- 444/1280)   0.0 E47-E40 S07-S05
      7  335 1205 (S11E49 L=175)   360     80        0 ( 297/1191- 379/1221)   7.1 E52-E44 S12-S09
      8  936 1147 (S02E08 L=216)   205     45        0 ( 905/1130- 961/1163)   0.0 E09-E06 S03-S01
      9  396 1045 (S19E41 L=183)   253     31        2 ( 377/1019- 415/1068) -56.3 E42-E38 S20-S17
     10  492  961 (S22E32 L=192)   461    147        0 ( 466/ 907- 534/1032) -40.6 E34-E27 S24-S18
     11  292  983 (S25E49 L=174)   279     32        0 ( 273/ 964- 311/1009) -33.7 E51-E47 S25-S23
     12  236  878 (S32E56 L=168)   335     15        0 ( 219/ 857- 256/ 896)  26.6 E57-E53 S33-S30
     13  219  785 (S38E60 L=164)   610     71        0 ( 179/ 771- 253/ 806)   9.5 E66-E54 S38-S36
  

    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