Teresa,
This is a quick and dirty script. There is no error checking!!
The script expects the site file to be in your current mapset.
Two arguments should be entered on the command line:
proj.site infile outfile
where:
infile is the name of the site_lists file to project
outfile is the name for the new projected site file.
The output file is placed under your current directory or the "full path"
Contents of file proj.site------------------
:
file=$1
outfile=$2
OUTPARMS="+proj=aea +ellps=clark66"
INPARMS=`$GISBASE/mapgen/interface/get_proj_inf`
sed -e "s/|/TAB/g" $LOCATION/site_lists/$file |proj $INPARMS -I - >tmp.$$
cat tmp.$$|proj $OUTPARMS -|sed -e "s/TAB/|/g" >$outfile
rm tmp.$$
------------------------end
note!!! be sure to use the TAB key in the sed statement.
Additionally, you could setup OUTPARMS as an argument (OUTPARMS=$3).
Then you would enter the command as:
proj.site infile outfile "+proj=aea +ellps=clark66"
Hope this helps!
Jill
__________________________________________
____ __ ______ _____ _______ .
/ \ // / // / // // // / \
/ \ // / // / // / \
/ \ // /____// / //____/// /.....\/\ /\
/ \// / \\ / // / \_/\_/ \
_/ \/__/ \\_/_____//______ // /\/ _/\ ---------
________________________________________/
Jill R. Schuler
USDA
Natural Resource Conservation Service
2627 Redwing, Suite 120
Ft. Collins, Co 80526 FAX: (970) 282-1930
VCOMM:9010-2418
e-mail: jschuler@itc.nrcs.usda.gov Phone: (970) 282-2418