Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  722 1747 (N41E44 L=320)   541     90        1 ( 691/1708- 748/1789) -66.8 E47-E41 N37-N44
      2 1164 1510 (N32W00 L=  4)   736    303        1 (1084/1387-1248/1616) -42.9 E07-W08 N25-N38
      3  494 1410 (N14E43 L=321)   631    136        2 ( 470/1341- 516/1472) -45.0 E47-E39 N09-N17
      4 1512 1261 (N22W28 L= 32)   705    123        4 (1468/1216-1547/1321)  54.5 W25-W30 N19-N26
      5 1403 1247 (N19W20 L= 24)   251    130        0 (1358/1213-1450/1277) -15.9 W16-W23 N17-N21
      6  714  913 (S12E17 L=347)   204    104        0 ( 704/ 866- 733/ 966) -59.0 E18-E15 S14-S09
      7  620  879 (S16E23 L=341)   743    223        2 ( 553/ 811- 669/ 938)  90.0 E27-E19 S20-S12
      8  217  786 (S29E61 L=303)  1687    231        6 ( 159/ 714- 280/ 864)  24.8 E79-E51 S34-S22
      9 1119  618 (S24W13 L= 17)  3258    853        8 ( 802/ 493-1458/ 837)  40.7 E06-W30 S36-S04
     10  566  518 (S41E26 L=338)  1179    195        4 ( 491/ 480- 653/ 560)  26.6 E34-E18 S45-S36
     11 1716  519 (S19W66 L= 70)  1241    141        2 (1648/ 492-1781/ 540)  17.4 W58-W76 S22-S16
     12 1450  412 (S31W46 L= 50)  1002    142        2 (1380/ 380-1522/ 453)  33.7 W40-W51 S34-S27
     13 1225  307 (S44W35 L= 39)   862    103        2 (1186/ 275-1271/ 344)  -5.7 W29-W39 S46-S40
  

    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