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

sisma66 at tiscali.it sisma66 at tiscali.it
Fri Aug 25 08:38:46 EDT 2006


Thanks for all information you told me, so I undestood:

1) this isn't a bug but is an R problem that he can't manage a vector 
so big
2)both the sizes of the data and newdata are the same of grass region 
(I hope this can help more)
3) reducing the resolution It works and also using maxdist
4) sorry for my stupid question, next time I will try to solve alone 
my issue

Thanks to all

Pierluigi

Il giorno gio, 24/08/2006 alle 18.36 +0200, Roger Bivand ha scritto:
On Thu, 24 Aug 2006, Edzer J. Pebesma wrote:
> 
> > Pierluigi De Rosa wrote:
> > > Dear all users of gstat
> 
> I think it is maxdist= again, or some other way of reducing the 
size. The 
> GRASS region is roughly 800 by 900, but the sizes of both the data 
and 
> newdata are unknown. Pierluigi also cross-posted, which does 
indicate a 
> propensity to appeal for help rather than to look at the archives. 
Using 
> the search facility on the GRASS site with "gstat memory", there 
were lots 
> of promising hits from this list just in the last few months.
> 
> > >
> > > 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!! :-)
> > >   
> > Perhaps he can, but remember he isn't a daily grass user. I doubt 
it is 
> > a bug because you do get a clear error message: you (through how 
you use 
> > gstat) ask the computer more memory than your computer has 
available, 
> > and that's what you are told. The question is why.
> > 
> > Without the commands how you retrieved your data and of what size 
> > (number of rows and columns) they are there is little I can do, 
from 
> > here, however.
> > --
> > Edzer
> > > 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)
> > >
> > > _______________________________________________
> > > statsgrass mailing list
> > > statsgrass at grass.itc.it
> > > http://grass.itc.it/mailman/listinfo/statsgrass
> > >   
> > 
> > _______________________________________________
> > statsgrass mailing list
> > statsgrass at grass.itc.it
> > http://grass.itc.it/mailman/listinfo/statsgrass
> > 
> 


Naviga e telefona senza limiti con Tiscali     
Scopri le promozioni Tiscali adsl: navighi e telefoni senza canone Telecom

http://abbonati.tiscali.it/adsl/




More information about the grass-stats mailing list