Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  797 1803 (N35E46 L=257)  2599    379        4 ( 664/1705- 878/1866)  25.0 E57-E33 N28-N40
      2  492 1682 (N19E60 L=243)   737    123        0 ( 467/1629- 517/1740) -48.8 E64-E56 N15-N23
      3  596 1599 (N16E45 L=258)   414    158        0 ( 559/1548- 629/1663) -71.6 E47-E43 N11-N20
      4 1288 1326 (N19W06 L=309)   366    153        0 (1251/1268-1331/1396) -32.0 W02-W10 N16-N21
      5 1562 1067 (N12W30 L=333)   557    283        0 (1472/ 975-1641/1178) -19.4 W21-W38 N09-N16
      6  263 1063 (S23E53 L=250)   377     66        0 ( 242/1025- 287/1097)  68.2 E55-E51 S25-S20
      7 1819  860 (N10W58 L=  0)   905    151        1 (1792/ 796-1847/ 942)  90.0 W55-W59 N06-N15
      8 1698  905 (N08W45 L=347)   362     65        2 (1673/ 874-1726/ 935)  68.2 W42-W46 N05-N10
      9  160  871 (S36E68 L=235)  1020    110        0 ( 125/ 820- 201/ 920)  22.4 E79-E60 S39-S32
     10  954  884 (S15E00 L=303)   310     84        0 ( 914/ 853- 994/ 913)  68.2 E01-W01 S17-S12
     11  313  874 (S33E46 L=257)   339     47        1 ( 286/ 855- 340/ 896)  -9.5 E48-E42 S33-S31
     12  537  728 (S36E23 L=280)   319     63        0 ( 502/ 716- 577/ 744)  15.9 E26-E19 S37-S34
     13 1801  667 (S00W65 L=  8)  1558    139        1 (1762/ 631-1854/ 733)  -3.2 W57-W76 S02-N01
     14  976  532 (S35W13 L=316)   832    224        1 ( 877/ 511-1082/ 561)  30.3 W07-W19 S38-S30
     15 1509  284 (S30W70 L= 13)  1318     99        2 (1464/ 259-1565/ 315)  22.6 W64-W79 S32-S27
     16 1347  238 (S39W62 L=  4)   571     72        0 (1323/ 221-1382/ 269)   0.0 W55-W66 S40-S37
  

    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