#! /bin/sh

DB=${3:-raster}

echo "Computing path from $1 to $2 in $DB"

cd $GEOHOME/demos/raster

/bin/rm oids_out

(
monitor -TN -c "retrieve (the_oid= roads.oid, frompnt= PointSelector(To2Pln(roads.geo_pline), 1),topnt= PointSelector(To2Pln(roads.geo_pline), -1), length= Length2Pln(To2Pln(roads.geo_pline)) / 1500.0, dir="1"::int4) where roads.category = 1" $DB

monitor -TN -c "retrieve (the_oid= roads.oid, frompnt= PointSelector(To2Pln(roads.geo_pline), 1),topnt= PointSelector(To2Pln(roads.geo_pline), -1), length= Length2Pln(To2Pln(roads.geo_pline)) / 1000.0, dir="1"::int4) where roads.category != 1" $DB
) |
$GEOHOME/analyses/cpath -d $1 $2

monitor $DB << EOF

delete PathOids \g
copy PathOids from "$GEOHOME/demos/raster/oids_out" \g

EOF