[STATSGRASS] The gstat R library's: memory problem

Pierluigi De Rosa sisma66 at tiscali.it
Thu Aug 24 10:29:33 EDT 2006


Dear all users of gstat

I'm using the gstat library for univariate geostatistics analysis,
I have to interpolate with kriging value of water level from wells to evaluate 
the water table.
I use the univarsal kriging and the indipendent variable is the digital 
elevation model (DEM)
If I use a DEM with resolution of 20 mt I have an error memory

This is the string command I use for do that:

OK_res20=krige(LEV_MAG91~DEM,loc=pozzi91_meno,newdata=DEM20,model=vgm_fit_mag91_DEM)

This is the error I have
Error: cannot allocate vector of size 3995610 Kb

If I use a DEM with resolution of 30 mt (using r.resample.rst)

(This is the string) 
OK_res30=krige(LEV_MAG91~DEM,loc=pozzi91_meno,newdata=DEM30,model=vgm_fit_mag91_DEM)

every thinks works ok, is this a bug?? Maybe  Dr. Pebesma can help me!! :-)

Thanks, Pierluigi

P.S.
this is my region setting

g.region -p
projection: 99 (Transverse Mercator)
zone:       0
datum:      rome40
ellipsoid:  international
north:      4829500
south:      4813540
west:       2280740
east:       2298420
nsres:      20
ewres:      20
rows:       798
cols:       884

--
Pierluigi De Rosa

Department of Enviromental & Civil engineering 
Faculty of Engineering 
Perugia (Italy)




More information about the grass-stats mailing list