#!/bin/sh -x

export RESOLUTION=1/256in

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

for file in N40E028 N40E029 N41E028 N41E029; do
  (echo P5 1201 1201 65535; cat data/$file.hgt) |
  pnmnorm -bvalue 0 -wvalue 300 |
  ppmchange white black | ppmtopgm > tmp/$file.pgm
done

out/arrange \
  tmp/N40E028.pgm@x=0,y=0 \
  tmp/N40E029.pgm@x=1200,y=0 \
  tmp/N41E028.pgm@x=0,y=1200 \
  tmp/N41E029.pgm@x=1200,y=1200 |
pnmcut 1018 1086 327 152 > tmp/istanbul.pgm &&

echo 'P2 1 1 1 1' > tmp/flat.pgm &&

out/arrange \
  tmp/istanbul.pgm@xsize=5in,ysize=3in,zsize=1+1/4in \
  tmp/flat.pgm@xsize=5+1/8in,ysize=3+1/8in,z=1+1/4in \
  > tmp/istgoal.pgm &&

out/arrange \
  tmp/istgoal.pgm \
  tmp/flat.pgm@xsize=5+1/8in,ysize=3+1/8in,z=0 \
  > tmp/iststate0.pgm &&

rm -rf tmp/debug* &&
out/toolpath -d tmp/debug1 -s 1/16in -t 1/8in \
  tmp/iststate0.pgm tmp/istgoal.pgm -o tmp/iststate1.pgm > tmp/istpath1.rml &&
out/toolpath -d tmp/debug2 -s 1/256in -t 1/8in \
  tmp/iststate1.pgm tmp/istgoal.pgm -o tmp/iststate2.pgm > tmp/istpath2.rml