#!/bin/sh -x

export RESOLUTION=1/256in

cd "`dirname $0`" &&
mkdir -p tmp &&

dataset=data/NED_55723387/ned_55723387/w001001.adf
gdal_translate -of PNM -ot UInt16 -scale 0 400 0 65535 "$dataset" tmp/sf.pgm &&
echo 'P2 1 1 1 1' > tmp/flat.pgm &&
pnmflip -r90 tmp/sf.pgm > tmp/sfrot.pgm &&
out/arrange \
  tmp/sfrot.pgm@xsize=2.5in,ysize=3in,zsize=3/4in,z=1/4in \
  tmp/flat.pgm@xsize=2.5in,ysize=3+1/8in,z=1in \
  tmp/flat.pgm@x=2+1/2in,xsize=1/4in,ysize=3+1/8in,z=7/4in \
  tmp/flat.pgm@x=2+3/4in,xsize=1/4in,ysize=3+1/8in \
  > tmp/sfgoal.pgm &&
out/arrange \
  tmp/sfgoal.pgm \
  tmp/flat.pgm@xsize=2+3/4in,ysize=3+1/8in,z=0 \
  tmp/flat.pgm@x=2+3/4in,xsize=1/4in,ysize=3+1/8in \
  > tmp/sfstate0.pgm &&
rm -rf tmp/debug* &&
CPUPROFILE=tmp/toolpath.prof \
out/toolpath -d tmp/debug1 -s 1/16in -t 1/8in \
  tmp/sfstate0.pgm tmp/sfgoal.pgm -o tmp/sfstate1.pgm > tmp/sfpath1.rml &&
out/toolpath -d tmp/debug2 -s 1/64in -t 1/8in \
  tmp/sfstate1.pgm tmp/sfgoal.pgm -o tmp/sfstate2.pgm > tmp/sfpath2.rml