Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1357 1677 (N47W34 L=340)   375     71        1 (1321/1656-1396/1695)   5.7 W29-W39 N45-N48
      2 1474 1626 (N42W43 L=349)   411     79        1 (1437/1609-1510/1647) -18.4 W39-W46 N40-N43
      3 1472 1221 (N15W30 L=336)   396    143        0 (1408/1183-1524/1260) -35.5 W26-W33 N12-N17
      4 1760 1174 (N10W54 L=  0)   315     50        1 (1741/1144-1775/1202)  56.3 W51-W55 N08-N11
      5 1117 1136 (N11W06 L=312)   454     93        0 (1094/1104-1157/1181) -59.0 W04-W08 N08-N13
      6 1800  801 (S14W58 L=  4)   613     80        1 (1773/ 754-1828/ 833) -39.8 W54-W62 S16-S11
      7 1732  782 (S14W50 L=357)   397    114        0 (1712/ 717-1743/ 827)  81.9 W49-W51 S18-S11
      8  362  647 (S20E50 L=256)   544    185        0 ( 306/ 568- 394/ 704) -73.3 E54-E47 S25-S15
      9  636  592 (S23E28 L=278)   300     59        1 ( 615/ 563- 659/ 624)  53.1 E30-E26 S24-S20
     10 1514  482 (S34W37 L=344)   360     99        0 (1484/ 464-1549/ 512) -23.2 W33-W40 S35-S31
  

    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