[GRASS-user] Fwd: GRASS 6.4 How to extract point-values from multiple rasters

Marco Sciaini sciaini.marco at gmail.com
Thu Apr 4 07:58:52 PDT 2013


Dear Lucien,

here is a little bash script, which updates the attribut-table from
your vector file and
add all extracted values from multiple rasters in a mapset.

for A_Raster in $( eval "g.mlist type=rast mapset="..." )
do
v.db.addcol map="..." columns="$A_Raster double precision" --overwrite
v.what.rast vector="..." raster=$A_Raster column=$A_Raster --overwrite
echo "Spatial query of $A_Raster finished at $(date)"
done

Kind regards
Marco


More information about the grass-user mailing list