Re: [R-sig-Geo] Alternate methods: Polygon Algebra / Polygon Overlay with R Spatial objects..

2010-01-15 Thread Roger Bivand
On Fri, 15 Jan 2010, rick reeves wrote: Hello List: I have been reviewing techniques within R for performing the GIS-like operation 'polygon algebra' (computing the union and intersection of polygons within two SpatialPolygonDataFrame objects). The goal is to combine two polygon data frames,

[R-sig-Geo] Alternate methods: Polygon Algebra / Polygon Overlay with R Spatial objects..

2010-01-15 Thread rick reeves
Hello List: I have been reviewing techniques within R for performing the GIS-like operation 'polygon algebra' (computing the union and intersection of polygons within two SpatialPolygonDataFrame objects). The goal is to combine two polygon data frames, to produce a new polygon data frames that

[R-sig-Geo] ordinary kriging and elevation with GSTAT

2010-01-15 Thread Tobin Cara
Hello, If I generate a variogram with an outside program on the x-y plane, would this variogram be valid for ordinary kriging in gstat? I am giving gstat the x,y,z coordinates and I am trying to predict at the elevations of a digital elevation model. Thank you in advance for taking the time to

Re: [R-sig-Geo] RE : Inverse Distance Weighting and Treating kriged data before cross validation

2010-01-15 Thread Edzer Pebesma
mykrige.cv = gstat:::krige.cv.spatial mykrige.cv = edit(mykrige.cv) then, you will find out that you will have to do the same with gstat.cv, and that gstat.cv is rather cryptic as it deals with cokriging of which ordinary kriging is a special case. Maybe it wasn't such a good idea after all i

[R-sig-Geo] RE : Inverse Distance Weighting and Treating kriged data before cross validation

2010-01-15 Thread Tobin Cara
Hello Prof. Pebesma, I would like to use leave-one-out cross validation to have a fair comparison with my other methods. Could you explain further how to copy krige.cv to my own function and modify? Can I open up the krige.cv code somewhere? I am still new to R. Thank you in advance, Cara __

Re: [R-sig-Geo] Inverse Distance Weighting and Treating kriged data before cross validation

2010-01-15 Thread Edzer Pebesma
Dear Cara, one way would be to copy krige.cv to your own function and modify this; I would recommend this if you want leave-one-out cross validation. If n-fold cross validation is good enough for a few folds, e.g. because you have plenty of data, then I would do the whole procedure "manually"