[GRASS-user] spgrass6 in GRASS script

Martin Wegmann wegmann at biozentrum.uni-wuerzburg.de
Thu Feb 22 08:44:28 EST 2007


Hello,

I am trying to insert R commands in a GRASS script.

The indirect way works:

...GRASS script ...
echo "
x <- read.table('$points.xy', sep=\"|\")
...
" > R_temp
echo 'source("R_temp")' | R --vanilla --slave
rm R_temp
... GRASS script ... end

but I would like to skip the export/import routines in the script and do it 
instead with spgrass6 like:

... GRASS script ...
R
require(spgrass6)
mydata <- readFLOAT6sp(c("myraster","myraster2"))
... R analysis ...
... GRASS script ...

but this approach fails with: 
mean_extract_R: line 293: syntax error near unexpected token `spgrass6'
mean_extract_R: line 293: `require(spgrass6)'

and:
echo "
require(spgrass6)
mydata <- readFLOAT6sp(c("myraster"))
" > R_temp
echo 'source("R_temp")' | R --vanilla --slave

fails with:
Loading required package: spgrass6
Loading required package: sp
Loading required package: maptools
Loading required package: foreign
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Error in seq(along = vname) : object "myraster" not found
Execution halted

So far I did not find a threat concerning R/spgrass6-GRASS communication in 
scripts.
I appreciate any help/links/documents where R spgrass6 commands are included 
in a GRASS script.

thanks in advance, Martin




More information about the grass-user mailing list