Here is a solution for a non-interactive conversion of sites data to raster
that was recommended by Michael Shapiro [I believe Michael has also placed
source code on moon which provides non-interactive versions of the individual
s.menu options]:
:
# received from shapiro 911202 as sites_to_cell solution
# then modified
USAGE="$0 [-q quadsize] [-t title] sites_list [raster]"
if [ $# -eq 0 -o ${1-""} = "help" ]
then
echo $USAGE
echo " sites_list = name of the input sites list"
echo " raster = name of the output raster file (default: same as sites)"
echo ' quadsize = size of each site in cells'
echo ' (default: "0" == 1 cell, see s.menu)'
echo ' title = title for raster map'
echo ' (default: "created from sites file $sites")'
exit -1
fi
# continue only if user is executing GRASS
if [ "$GISRC" = "" ]
then
echo "Sorry, you are not running GRASS" >&2
exit 4
else
USERGRASS="$GISDBASE/$LOCATION_NAME/$MAPSET"
fi
# set defaults, then test for user-specified settings
while [ `echo $1 | cut -c1` = "-" ]
do
case $1 in
-q) quadsize=$2
;;
-t) title=$2
;;
*) echo "$1 is not a valid flag"
exit -1
;;
esac
shift 2
done
sites=$1
raster=$2
# remove the raster map to be safe
g.remove ${raster:=$sites} >/dev/null
# output sites file to a file in /tmp
tmp=/tmp/$$
rm -f $tmp
s.out.ascii -ad $sites fs='|' > $tmp
# run the program which converts sites to cell
# the command line has unix-path name of sites file
# the input is
# quadsize
# name of output raster map
# title for raster map
echo "converting sites list $sites to raster $raster"
$GISBASE/etc/sites/sites_to_cell $tmp << EOF
${quadsize:=0}
$raster
${title:="created from sites file $sites"}
EOF
status=$?
rm -f $tmp
exit $status
******************************************************************************
Greg Koerper Internet: greg@towhee.cor2.epa.gov
ManTech Environmental Technology, Inc. UUCP: hplabs!hp-pcd!orstcs!koerper
200 SW 35th St., JSB
Corvallis, OR 97333 "The 90s will make the 60s look like the 50s.
(503) 754-4490 Just ask your kids."
fax: (503) 754-4338
******************************************************************************