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 _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
