Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  721 1813 (N42E43 L= 92)  1372    126        2 ( 669/1766- 771/1851) -10.0 E53-E36 N37-N45
      2  980 1683 (N35E15 L=120)   252     38        1 ( 961/1664-1004/1704) -14.0 E16-E12 N33-N36
      3 1138 1662 (N36E03 L=132)   288     50        0 (1121/1628-1152/1690) -78.7 E04-E02 N33-N38
      4  473 1641 (N25E53 L= 83)   651     79        1 ( 438/1604- 505/1672) -20.0 E58-E47 N23-N27
      5 1512 1637 (N43W28 L=163)   555     49        1 (1491/1606-1539/1664)  50.2 W25-W31 N39-N45
      6 1630 1335 (N23W36 L=171)  1038    306        0 (1543/1224-1701/1458) -42.7 W28-W41 N17-N29
      7 1865 1326 (N29W63 L=198)   771     82        0 (1823/1301-1901/1349)   8.1 W55-W69 N26-N30
      8  578 1062 (S10E28 L=107)   344     74        1 ( 538/1038- 615/1086)  21.8 E30-E25 S11-S08
      9  108  995 (S17E75 L= 60)  4341    190        2 (  78/ 885- 154/1072)  61.4 E85-E66 S23-S12
     10  359  899 (S22E44 L= 92)   477    115        0 ( 305/ 861- 399/ 927) -14.0 E48-E40 S23-S20
     11  845  884 (S17E09 L=127)   286     65        1 ( 807/ 863- 873/ 911)  45.0 E11-E07 S18-S14
     12 1547  888 (S05W35 L=170)   216     40        0 (1521/ 865-1567/ 910) -26.6 W32-W36 S06-S04
     13 1488  687 (S18W36 L=171)   364     72        1 (1454/ 655-1517/ 723)  78.7 W33-W37 S20-S15
     14 1528  584 (S23W43 L=178)   494     88        1 (1493/ 561-1575/ 611)  -7.1 W39-W47 S24-S22
     15  579  438 (S51E29 L=107)   592     87        0 ( 556/ 392- 601/ 483)  81.9 E30-E26 S54-S47
     16 1053  307 (S52W20 L=156)   369     70        0 (1015/ 293-1094/ 326)  15.9 W16-W23 S53-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