I am looking for suggestions on a query of worldclim data
(http://worldclim.org/).
I have used something similar to below for previous data extraction where I am
looking for data associated with a particular site, or list of sites, where a
single lon lat coordinate will return a single value. However, in this case, I
need to pull an entire longitudinal band for a particular latitude.
For example instead of grabbing data from a single coordinates with "east_north
= -92.50, 46.51”, I want to grab all the data in the longitudinal band at
latitude 46.51.
If it were possible something like “east_north = *, 46.51”
Any suggestions on this point would be much appreciated. Thanks in advance -
Kirk
==============================================
if test "$GISBASE" = ""; then
echo "You must be in GRASS to run this program"
exit
fi
# Use input text file for coords in the format: lon lat (easting northing)
single
# space between coords. Example: -98.42072 55.91481. Must use real coordinates,
no blanks.
# Script will create file for each $MAP in the list. g.mlist can also take the
"pattern"
# argument which takes regular expressions. For example: pattern="*" returns
all
# maps in the database. If script fails check for hidden characters in text
file, best to
# use vi or emacs to eliminate them
for MAP in `g.mlist type=rast pattern=“HIST_tmin*"`; do
r.what --verbose input="$MAP" < ~/Desktop/ai/coords.txt >
~/Desktop/ai/"$MAP".txt;
echo "$MAP"
done
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
