Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread Edzer J. Pebesma
temiz wrote: > remove.duplicates(el3, zero = 0.0, remove.second = TRUE) Error: cannot allocate vector of size 2259572 Kb Is there another way rather than choosing smaller region ? Yes, wait for sp 0.9-12, which is "on its way to CRAN", where zerodist should work without a large memory footpr

Re: [R-sig-Geo] Good projection for N/S America?

2007-04-06 Thread Roger Bivand
Since this topic is of general interest, I've made an exception and allowed (this once!) a posting of more than 200K. In general, if graphics are big, please consider either an alternative device (png is often OK), or posting just a URL to the real file. With apologies to list members on dial-u

Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread temiz
Edzer J. Pebesma wrote: > Most probably you have duplicate points in your observations. > > You may try to use zerodist() to find them, or remove.duplicates() to > remove them but this may currently run out of memory as the current > implementation forms the full distance matrix. In that case you

Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread Edzer J. Pebesma
Most probably you have duplicate points in your observations. You may try to use zerodist() to find them, or remove.duplicates() to remove them but this may currently run out of memory as the current implementation forms the full distance matrix. In that case you'd have to go through it in part

Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread temiz
Edzer J. Pebesma wrote: > temiz wrote: >> ... >> I regarded some modifications in >> http://grass.gdf-hannover.de/wiki/GRASS_6_Tutorial#GRASS_and_R_kriging_interpolation >> >> >> >> I encountered a new problem : >> > OK_pred <- krige(z ~ 1, el3, newdata=mask_SG, model=efitted) >> [using ordinar

Re: [R-sig-Geo] [R] how to image.plot a XY grid file into a lat-lon map

2007-04-06 Thread Edzer J. Pebesma
May I suggest you take a look at the Spatial Task View? r-project.org -> CRAN -> Select Mirror -> Task Views -> Spatial. -- Edzer Linda Smith wrote: > Could you please let me know more details about rgval package? I have never > used it before. I only tried image.plot() where if you have regular

Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread Edzer J. Pebesma
temiz wrote: > ... > I regarded some modifications in > http://grass.gdf-hannover.de/wiki/GRASS_6_Tutorial#GRASS_and_R_kriging_interpolation > > I encountered a new problem : > > OK_pred <- krige(z ~ 1, el3, newdata=mask_SG, model=efitted) > [using ordinary kriging] > > "memory.c", line 57: can't

Re: [R-sig-Geo] problem in krige function

2007-04-06 Thread temiz
Roger Bivand wrote: > On Thu, 5 Apr 2007, temiz wrote: > > >> following the topic GRASS and R kriging interpolation from grass6 tutorial >> >> I came here without problems >> >> But here: >> >> zz <- krige(z ~ 1, locations = el1 , newdata = mask_SG , model = efitted) >> Error in predict.gstat(g,