Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  703 1786 (N39E51 L=140)  1408    164        0 ( 634/1735- 762/1839)  -8.7 E68-E41 N36-N42
      2  731 1640 (N29E37 L=155)   951    238        1 ( 614/1598- 813/1690)  47.7 E43-E33 N23-N34
      3 1561 1518 (N42W29 L=220)  1661    248        2 (1484/1415-1617/1609) -50.2 W21-W34 N35-N48
      4 1493 1348 (N28W23 L=215)  1871    575        0 (1323/1130-1575/1542) -40.0 W07-W32 N16-N37
      5  413 1497 (N13E54 L=137)   580     93        2 ( 384/1451- 447/1537) -81.9 E55-E52 N09-N16
      6 1154 1343 (N19W01 L=192)   240     76        1 (1124/1321-1179/1384)  76.0 E01-W02 N17-N22
      7 1806 1181 (N25W54 L=245)  1052    119        2 (1778/1122-1828/1237)  76.0 W50-W56 N21-N29
      8  260 1036 (S17E54 L=137)  2112    400        3 ( 170/ 890- 335/1144)  50.2 E66-E47 S27-S09
      9  495  900 (S21E31 L=160)   388    146        2 ( 439/ 863- 537/ 959)  45.0 E35-E28 S23-S16
     10 1266  851 (S07W18 L=209)   228     99        0 (1222/ 832-1311/ 872)  38.7 W15-W20 S08-S04
     11 1488  803 (S05W33 L=225)   226     93        1 (1447/ 782-1532/ 819)  31.0 W31-W36 S06-S02
     12 1435  733 (S10W32 L=223)   255     79        1 (1413/ 703-1456/ 770) -53.1 W30-W33 S11-S07
     13  882  713 (S24E02 L=190)   310    122        0 ( 824/ 681- 921/ 756)  56.3 E04-W00 S26-S20
     14  348  621 (S43E46 L=145)  1151    142        4 ( 280/ 592- 415/ 656)  14.0 E57-E37 S45-S39
     15 1316  630 (S19W28 L=219)   324     88        1 (1268/ 618-1355/ 653)  26.6 W24-W30 S20-S16
     16 1454  474 (S25W45 L=236)   644    236        0 (1363/ 400-1548/ 534)  63.4 W41-W49 S31-S18
  

    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