Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1769 (N41E20 L=261)   420     53        1 (1044/1748-1092/1794)  71.6 E22-E18 N38-N42
      2 1157 1702 (N38E11 L=271)   734     99        1 (1108/1665-1205/1738)   0.0 E14-E06 N36-N40
      3  818 1668 (N26E31 L=251)   356     98        1 ( 806/1625- 833/1712) -50.2 E33-E28 N23-N29
      4  585 1509 (N11E39 L=243)  2204    531        3 ( 517/1361- 721/1672) -29.1 E52-E25 N04-N20
      5 1125 1301 (N11E01 L=281)  1332    414        4 (1086/1125-1158/1523) -71.6 E06-W04 N00-N24
      6 1842 1030 (N16W54 L=336)  1381    183        2 (1802/ 946-1876/1112) -35.0 W48-W58 N11-N19
      7  341  971 (S25E43 L=238)  1034    262        0 ( 247/ 910- 432/1035)  33.7 E51-E36 S29-S19
      8  705 1012 (S15E18 L=264)   189     74        0 ( 666/ 998- 741/1026)  11.3 E20-E15 S15-S14
      9  139  828 (S36E63 L=219)  1052    284        0 ( 106/ 795- 176/ 875)  32.0 E66-E58 S37-S32
     10 1805  814 (N03W57 L=339)   658    129        2 (1774/ 787-1841/ 869) -26.6 W54-W61 N01-N05
     11 1470  582 (S20W40 L=322)   371     71        0 (1447/ 543-1488/ 623) -80.5 W38-W41 S22-S16
  

    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