Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1496 (N32E02 L=287)   374    100        0 (1061/1443-1092/1540) -71.6 E03-W00 N28-N34
      2 1333 1194 (N15W18 L=307)   249     80        1 (1312/1164-1351/1231) -53.1 W16-W19 N12-N16
      3 1679 1107 (N13W44 L=333)   661    259        0 (1618/ 998-1725/1205) -53.1 W38-W48 N06-N18
      4 1878 1063 (N11W66 L=355)  1694    221        2 (1859/ 959-1891/1172)  69.0 W63-W69 N05-N18
      5  147 1006 (S10E70 L=220)   903    118        1 ( 102/ 975- 197/1041)  -5.7 E81-E61 S10-S07
      6  832  882 (S09E10 L=279)   624    165        1 ( 794/ 806- 874/ 955)  78.7 E12-E07 S14-S04
      7  340  856 (S17E46 L=243)  1055    213        4 ( 271/ 779- 383/ 909) -69.4 E53-E42 S21-S13
      8  216  792 (S22E62 L=228)  1889    188        6 ( 174/ 714- 267/ 880) -74.7 E68-E55 S27-S16
      9  305  742 (S24E51 L=238)   894    137        2 ( 276/ 672- 327/ 800) -82.9 E54-E49 S28-S20
     10 1551  572 (S22W43 L=333)   490    166        0 (1490/ 544-1600/ 624)  26.6 W37-W46 S23-S17
     11 1274  399 (S37W28 L=317)   299     36        0 (1263/ 370-1285/ 432) -68.2 W26-W29 S39-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