[STATSGRASS] Grass link to Gstat

[Edzer J. Pebesma]

Agustin Lobo wrote:
> 
> I've been succesful to get a variogram from
> a sites file, but not from a raster file.
> I do:
> GRASS:~/grass/gavarres_gis/h21/cell > gstat -i
> 
> and then:
> 1. Select
> 2. Choose identifier or enter new one []: a
> (where a is just a dummy name)
> 3.Enter data file name [] : test
> (where test is a raster file in $LOCATION/cell
> 4.Enter variable column number [0]:
> (I let it to be 0 as I don't know what to input, this
> is a raster file, not an x,y,(z),value file)
> 5. Hit return for Log transf. and for Average dups.
> 
> and the program quits with no comment or clue.
> 
> I'm sure I'm doing something wrong,
> but what? how can I actually tell gstats
> that my file is a raster file and not a sites file?
> 
I see what you mean -- gstat usually reads (point) data
from both lists or from maps, but the grass version does
not; the actual error message is hidden when the interactive
user interfaces jumps back to normal screen.

One solution is to convert your raster data to a sites
list. The other solution is downloading gstat-2.4.2
sources (now on gstat.org) which tries to do so; but 
please configure with:

./configure --with-grass --without-gdal

as current, highly experimental, gdal support still seems 
to mess things up.
--
Edzer