On 07/25/2015 08:08 PM, patrick s.
wrote:
Dear all
Is there a function in GRASS70 to match multiple rasters as values
to a point-vector. Some of the rasters are of type integer, some
of type float and the columns should ideally respect the datatype.
It might be possible as a bash.loop on v.what.rast, but I do not
find a way to query the datatype of a raster and get warnings:
You can use r.info -g to get information about each rast, including
datatype.
So:
eval `r.info -g $i`
will put all raster parameters into shell variables, among them
datatype.
if ( $datatype == 'CELL' ) ...
for i in `g.list type=rast ` #ideally only type integer
OR float
do
echo "adding layer '$i'"
v.db.addcolumn map=pt_grid column="$i double" #adapt
datatype of column according to raster type
v.what.rast map=pt_grid rast=$i col=$i
done;
Happy about any help and hint,
Patrick
This mail was received via Mail-SeCure System.
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
This mail was received via Mail-SeCure System.
|
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
