Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1858 (N72E35 L=151)   774     92        0 ( 740/1820- 783/1885) -29.1 E49-E25 N66-N76
      2  410 1488 (N38E50 L=137)   332     46        1 ( 388/1461- 428/1510) -56.3 E52-E46 N36-N39
      3 1791 1479 (N24W74 L=260)   771     63        2 (1770/1456-1812/1510)  76.0 W69-W79 N22-N26
      4  649 1465 (N36E25 L=161)   267     70        1 ( 613/1449- 683/1481)  -9.5 E28-E22 N34-N36
      5 1149 1471 (N32W13 L=199)   338    123        1 (1087/1458-1209/1481)   0.0 W08-W17 N31-N33
      6 1054 1369 (N26W05 L=191)   544    179        1 (1022/1286-1084/1446) -90.0 W03-W06 N20-N31
      7  346 1367 (N30E52 L=135)   372     80        1 ( 337/1323- 363/1409) -68.2 E53-E49 N27-N32
      8 1456 1371 (N23W33 L=219)   310    107        0 (1411/1339-1501/1398) -38.7 W30-W36 N20-N24
      9 1754 1220 (N10W55 L=241)   818    133        2 (1732/1164-1779/1285) -90.0 W52-W56 N06-N13
     10  706 1166 (N16E19 L=167)   708    226        2 ( 647/1088- 760/1224) -76.0 E23-E15 N11-N19
     11 1751  858 (S12W50 L=237)   646     59        3 (1717/ 830-1781/ 886)   9.5 W47-W53 S14-S10
     12 1725  710 (S22W50 L=236)   458     88        0 (1685/ 675-1757/ 747)  33.7 W45-W52 S23-S19
     13 1538  638 (S24W33 L=219)  1360    300        0 (1463/ 506-1592/ 745)  61.9 W28-W36 S33-S17
     14 1769  668 (S26W56 L=243)  1138    191        0 (1702/ 604-1834/ 736)  23.2 W49-W63 S29-S21
     15  233  630 (S17E70 L=117)  1334    113        2 ( 198/ 598- 280/ 676) -15.3 E76-E62 S19-S13
     16  885  595 (S21E13 L=174)   464    154        0 ( 815/ 554- 952/ 635)  26.6 E17-E07 S23-S18
     17  973  443 (S33E09 L=177)   273     75        0 ( 955/ 403- 989/ 476)  59.0 E10-E07 S35-S30
  

    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