Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  465 1665 (N26E62 L=227)  5071    257        7 ( 424/1556- 518/1763) -32.0 E79-E52 N18-N33
      2  713 1621 (N28E37 L=252)   299     78        0 ( 681/1596- 744/1640)  11.3 E38-E33 N26-N29
      3 1676 1524 (N46W44 L=333)   610     94        0 (1628/1495-1716/1546)   0.0 W37-W49 N44-N47
      4 1502 1396 (N31W24 L=313)   577    182        0 (1449/1320-1561/1459) -35.0 W19-W29 N27-N34
      5 1324 1413 (N28W11 L=300)   816    272        1 (1207/1346-1433/1449)  -3.4 W02-W19 N26-N30
      6 1386 1177 (N15W19 L=308)   210     65        0 (1363/1145-1412/1207) -45.0 W16-W20 N13-N16
      7  591 1011 (S11E26 L=263)   333     85        0 ( 548/ 992- 631/1031)  31.0 E28-E23 S12-S09
      8 1865  960 (N13W64 L=352)   821    113        0 (1829/ 905-1889/ 998) -39.8 W58-W67 N09-N14
      9  289  766 (S33E51 L=238)   524    131        0 ( 241/ 718- 343/ 818)  29.7 E58-E44 S36-S28
     10  783  735 (S24E09 L=280)   232     80        1 ( 747/ 718- 828/ 756)   9.5 E11-E05 S24-S22
     11 1535  662 (S12W41 L=330)   456    110        0 (1485/ 638-1585/ 677)  14.0 W37-W45 S13-S10
     12 1701  647 (S09W56 L=345)   395     44        0 (1685/ 616-1720/ 674) -76.0 W54-W57 S11-S07
     13 1539  478 (S23W51 L=340)  1574    209        1 (1491/ 392-1609/ 562) -67.4 W48-W56 S28-S16
     14  569  505 (S45E23 L=266)   476     64        0 ( 544/ 473- 596/ 534)  45.0 E25-E20 S47-S42
     15 1033  476 (S35W14 L=303)   591    176        0 ( 958/ 456-1113/ 506)  26.6 W08-W18 S37-S31
     16  347  465 (S55E58 L=231)   667    105        1 ( 311/ 433- 388/ 500)  12.9 E80-E45 S58-S50
  

    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