Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1685 (N28W56 L=238)   871    168        1 (1400/1661-1563/1708) -33.7 W51-W63 N24-N32
      2 1147 1678 (N40W32 L=214)   299     68        0 (1117/1654-1183/1698) -11.3 W29-W35 N38-N40
      3 1248 1448 (N23W26 L=208)   377    102        1 (1191/1422-1287/1466) -38.7 W22-W28 N20-N24
      4  656 1082 (N18E20 L=162)   501    191        0 ( 584/1050- 767/1109)  -8.7 E25-E12 N16-N19
      5 1289  884 (S09W11 L=193)   252     65        0 (1272/ 849-1305/ 918)  76.0 W09-W11 S10-S06
      6  652  810 (N03E27 L=155)   440    139        0 ( 618/ 746- 679/ 875)  81.9 E27-E25 N00-N07
      7  188  780 (N11E68 L=114)  1401    134        4 ( 154/ 717- 214/ 837)  53.1 E72-E64 N06-N14
      8  472  816 (N08E39 L=143)   305     65        2 ( 442/ 796- 502/ 833) -26.6 E40-E36 N06-N09
      9  620  331 (S25E51 L=131)  1344    220        1 ( 578/ 230- 660/ 435)  53.1 E58-E45 S33-S17
     10 1492  367 (S51W14 L=196)  1165    146        2 (1424/ 339-1565/ 405) -14.0 W06-W22 S53-S47
  

    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