Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  547 1602 (N15E48 L=263)   339     77        1 ( 527/1569- 570/1626) -36.9 E50-E46 N13-N16
      2  329 1536 (N07E65 L=246)  1625    148        1 ( 279/1467- 374/1613)  90.0 E67-E62 N01-N12
      3  780 1565 (N19E31 L=280)   263    102        1 ( 739/1534- 817/1602)  -8.1 E34-E27 N17-N20
      4  961 1414 (N15E15 L=296)   336    106        0 ( 930/1368- 982/1467) -80.5 E16-E13 N11-N17
      5 1537 1072 (N12W28 L=339)   384    106        0 (1512/1028-1574/1121) -26.6 W25-W31 N09-N13
      6 1663 1087 (N17W37 L=347)   493    145        1 (1623/1052-1714/1116)  33.7 W34-W40 N13-N18
      7 1693  884 (N06W45 L=355)  2077    325        2 (1613/ 777-1741/ 964) -28.3 W36-W51 N01-N10
      8 1848  828 (N09W62 L= 12)   580    106        0 (1821/ 787-1870/ 861) -20.6 W57-W65 N07-N10
      9  552  700 (S37E22 L=289)   534    146        0 ( 496/ 666- 622/ 734)  41.6 E25-E16 S40-S32
     10 1683  437 (S16W67 L= 18)  1179    256        1 (1658/ 382-1698/ 505) -22.6 W61-W73 S18-S12
  

    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