How to get hypsometry and zonal stats

1. make swaths in ArcEdit
% you can use copy parallel to baselines to get better positions
(in grid)
create zz1 def_poly
generate zz1
 fishnet
 origin coords: x,y
 y-axis coords: x,y
 length,width
 nrows,ncols:   25,1
 (done)

build
clean
save

2. clip and convert to grid in ARC
clip zz1 def_poly zz1_poly poly #
polygrid zz1_poly zz1_grid zz1_poly-id # #
 % where poly-id is the attribute to assign to the grid

ZONAL STATS
3. get and output stats in GRID and TABLES
(in grid)
zz1_stat = zonalstats(zz1_grid,elv,all,data)
(in tables)
 sel zz1_stat
 % need to expand area,sum fields
 alter
  Enter item name: area
  Item name: area
  Output width: 12
  blah
  Item decimal places: 0
  blah

  Enter item name: sum
  Item name: sum
  Output width: 12
  blah
  Item decimal places: 0
  blah

(back in tables)
% to unload all of the info in columnar format
unload zz1_stat.dat columnar init
 
HYPSOMETRY
4. run aml script to i) select elevations in each swath, ii) unload 
value, count to file 

Use the following script for this:


  &do i &list [listunique GRIDNAME.vat -info value]
    /* get the grid
    tmpgrid = con(GRIDNAME == %i%,elv)

    /* in tables
    arc tables
    sel tmpgrid.vat
    unload tmp.dat columnar init
    &sys mv tmp.dat hyps/GRIDNAME/hyps_%i%.dat
    quit

    /* back in grid
    kill tmpgrid

    &ty %i%
  &end

5. use Matlab scripts to load, average, clip, and plot data