Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  866 1584 (N21E26 L=280)   626    119        1 ( 814/1530- 902/1623) -20.6 E30-E22 N18-N23
      2 1675 1277 (N28W34 L=340)  1784    389        0 (1635/1142-1745/1429) -51.3 W28-W42 N19-N37
      3 1866 1246 (N34W56 L=  3)   922    147        1 (1819/1176-1899/1304) -28.6 W49-W61 N30-N36
      4  173  996 (S29E64 L=242)  1696    194        1 (  86/ 969- 272/1020)   4.8 E87-E51 S30-S25
      5  831  832 (S23E06 L=301)   266     90        0 ( 812/ 789- 857/ 872) -53.1 E07-E02 S24-S20
      6 1707  849 (N05W47 L=353)   402     70        1 (1668/ 834-1751/ 864)  23.2 W43-W50 N02-N06
      7  661  835 (S27E17 L=290)   234     79        0 ( 632/ 795- 691/ 862)  71.6 E17-E14 S30-S24
      8 1125  652 (S24W17 L=324)   590    218        0 (1027/ 624-1233/ 681)  38.7 W12-W22 S27-S19
      9 1463  554 (S19W42 L=349)   417     67        1 (1438/ 522-1484/ 587) -76.0 W40-W45 S20-S16
     10 1555  502 (S18W52 L=358)   340     51        2 (1528/ 477-1581/ 520)  71.6 W49-W53 S19-S16
     11 1546  429 (S22W57 L=  3)  1564    183        2 (1442/ 396-1610/ 473)  42.0 W50-W62 S27-S18
     12  666  354 (S60E02 L=304)   372     65        0 ( 625/ 342- 701/ 369)  26.6 E07-W01 S62-S57
  

    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