[STATSGRASS] kriging, gstat, GRASS

Thomas Adams Thomas.Adams at noaa.gov
Thu Dec 15 13:00:16 EST 2005


Edzer,

I'm re-visiting an old issue. Previously, using GRASS 5.3 & 5.4, I had 
no problems with using gstat & GRASS together. Now, I've gone to GRASS 
6.1 and gstat 2.4.5 (compiled with GRASS support) and I get errors. I 
apologize that I may have contacted you before about this, but I can not 
find any old emails. Attached is my gstat input file and a file with the 
gstat error I get from within my GRASS6.1 session. I'm getting:

WARNING: Vector used instead of sites.

I know that the R sp package will read GRASS vector location data and I 
could go that route if need be, but then there is the issue of getting 
the grids generated through kriging back into GRASS, which (I think) is 
more tedious.

Regards,
Tom



Edzer J. Pebesma wrote:

> Thomas Adams wrote:
>
>> Horacio,
>>
>> The best way to do this is to use gstat. Go to: 
>> http://grass.itc.it/statsgrass/index.html. Scroll to near the bottom 
>> of the page and follow the "gstat offers GRASS support" link 
>> (http://www.gstat.org/grass.html). I have used gstat with GRASS for a 
>> number of projects and it works wonderfully.
>>
>> Regards,
>> Tom
>
>
> Tom, thanks for the good news. There's two ways to get gstat's
> functionality in GRASS: either through the gstat binary reading
> and writing GRASS data, or through Grass-R, using the
> gstat R package. The latter will be easier if you're used to R.
> -- 
> Edzer
>


-- 
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL:	thomas.adams at noaa.gov

VOICE:	937-383-0528
FAX:	937-383-0033

-------------- next part --------------
# file: gstat.cmd.temps
#
# gstat command file
#
data(temp_20050801_08): 'temp_20050801_08', x=1, y=2, v=1;
variogram(temp_20050801_08): 0.2772 Nug(0) + 0.85 Sph(987521);
-------------- next part --------------
GRASS-GRID > gstat gstat.cmd.temps
gstat: Linux version 2.4.5 (17 June 2005)
Copyright (C) 1992, 2004 Edzer J. Pebesma
using Marsaglia's random number generator
data(temp_20050801_08): WARNING: Vector used instead of sites.
WARNING: Invalid absolute datetime format
Memory fault 



More information about the grass-stats mailing list