[GRASS-user] GRASS6.4 Problem with v.db.addcol

[BLANDENIER Lucien]

Dear all,

As I explained in my previous post, I have some troubles using v.db.addcol (see below). After making different tests, I realized there is a problem using script with v.db.addcol. I tried the following script:

#!/bin/bash
v.db.addcol map="extracting_point_ETo_RFE_test" columns="RF_2001001 double precision" --overwrite

And I received the following error message
GRASS 6.4.2 (Merti_WGS84_dd)> bash list_extract_RFE.sh
list_extract_RFE.sh: line 2: v.db.addcol: command not found

The command is working when I write it in the shell... but I don't want to do it 2000 times!!!
Im using GRASS 6.4 on windows.

do someone has an idea where the problem comes from?

thanks


Lucien
________________________________________
De : grass-user-bounces at lists.osgeo.org [grass-user-bounces at lists.osgeo.org] de la part de BLANDENIER Lucien [lucien.blandenier at unine.ch]
Date d'envoi : mercredi 10 avril 2013 16:09
À : Marco Sciaini; grass-user at lists.osgeo.org
Objet : [GRASS-user] RE : Fwd: GRASS 6.4 How to extract point-values    from    multiple rasters

Dear Marco,

Thank you for your answer and your script. I tried it but I meet some problem.

I used the script :

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

and I receive the following error message for each raster of my mapset :

extract_PointValue_from_RFE_raster.sh: line 13: v.db.addcol: command not found
ERREUR :Colonne <RFE_2001001> non trouvÚe
Spatial query of RFE_2001001 finished at Wed Apr 10 13:54:36 GMT 2013
extract_PointValue_from_RFE_raster.sh: line 13: v.db.addcol: command not found
ERREUR :Colonne <RFE_2001002> non trouvÚe
Spatial query of RFE_2001002 finished at Wed Apr 10 13:54:36 GMT 2013
...


I separated commands (v.db.addcol extracting_point_ETo_RFE_test columns="test double precision") and it worked

Do someone has an idea of the problem?

thanks


Lucien

________________________________________
De : grass-user-bounces at lists.osgeo.org [grass-user-bounces at lists.osgeo.org] de la part de Marco Sciaini [sciaini.marco at gmail.com]
Date d'envoi : jeudi 4 avril 2013 16:58
À : grass-user at lists.osgeo.org
Objet : [GRASS-user] Fwd: GRASS 6.4 How to extract point-values from    multiple rasters

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