[R-sig-Geo] Global variance estimation for evaluating resulting krigging quality

2011-05-01 Thread Sébastien Durand
Hello to all, I have been trying to find the following since quite a while but I guess I was simply using the wrong terms to express myself. The following shall describe what I am really looking for in R. I wish to know if there is an easy way to obtain that value using already preexisting

[R-sig-Geo] Minor tick marks in spplot

2011-03-27 Thread Sébastien Durand
Hello to all, I wish to add minor tick mark to the already existing major tick marks on the following plot I wish those minor tick marks to be shorter and grey Here is an example... data(meuse.grid) m = SpatialPixelsDataFrame(points = meuse.grid[c(x, y)], data = meuse.grid) spplot(m,

[R-sig-Geo] R2 values from SSErr fit.variogram attribute

2011-03-15 Thread Sébastien Durand
Thanks a lot, I understand what you meant now. I was like you said the general the fact that the computed means used to produce de SSErr value was obviously weighted and there to get a true R2 the same weighted means had to be used in the computation for the SSTot... Thanks a lot for your

Re: [R-sig-Geo] R2 values from SSErr fit.variogram attribute

2011-03-14 Thread Sébastien Durand
Hello Dr. Pebesma and all other 1989 readers... Believe me, I would have been much happier, if I could comprehend this tiny little detail that you refer to as : It looks quite OK; I suspect however that instead of using mean(s.v$gamma) for SSTot, for this mean the same weighting scheme

Re: [R-sig-Geo] Polygon dissolution or union issue

2010-12-05 Thread Sébastien Durand
]]@polygons[[...@coords - round(p...@polygons[[i]]@polygons[[...@coords, digits) } ### union - unionSpatialPolygons(pol, ID=rep(1, times=length(p...@polygons))) plot(r) plot(union, add=T) Robert On Sat, Dec 4, 2010 at 10:00 PM, Sébastien Durand v8ex...@gmail.com wrote: Hello

[R-sig-Geo] Polygon edge smoothing

2010-12-05 Thread Sébastien Durand
Hello to all, Is there an elegent way or an already exisiting function that would allow to smooth out the edge of a polygons (an object of class SpatialPolygons SP) In my case I have a polygon edge made out of stair cases... so if would appear much better smoothed... Thanks to all. S.

Re: [R-sig-Geo] Creating polygons out of an irregular data set to represent specific classes

2010-12-02 Thread Sébastien Durand
Ok here is my latest development on the topic! But I still have a question! At the end I still get segmented polygons! library(akima) library(maptools) library(raster) data(akima) # define de dimension of grid dimgrid=50 # Interpolate to regular grid ak.li=interp(akima$x, akima$y, akima$z,