#! /bin/sh

DB=${3:-raster}

echo "Computing district from $1 with max distance $2 in $DB"

cd $GEOHOME/demos/raster

/bin/rm -f 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) 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) where roads.category != 1" $DB
) |
$GEOHOME/analyses/district $1 $2


monitor $DB << EOF

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

EOF

exit 0

# Future extension when indices work OK...

monitor -c "destroy PathOids" $DB
$GEOHOME/analyses/makePathOids $DB

monitor $DB << EOF

copy PathOids from "$GEOHOME/demos/raster/oids_out" \g
define index PathOidsIndex on PathOids using btree (the_oid oid_ops) \g

EOF