Re: [R-sig-Geo] How to use coKriging with LANDSAT ETM Imagery

2015-08-31 Thread Edzer Pebesma
Please find a working script below. The error message was given because the variable "value" was not present in new.data; the variable was called "Sample1_Band3". It is quite unclear what you try to do: interpolating an image at the image pixels themselves (which will, by definition, return the

Re: [R-sig-Geo] How to use coKriging with LANDSAT ETM Imagery

2015-08-29 Thread Rajabhushanam. C
Hi, I am using a Landsat ETM Scene imagery for calculating the Linear Corregionalization Method in CoKriging. (gstat). Can you kindly debug my snippet (See enclosed). Really appreciate your help. Thanking in advance. coKrige1 - variogram(coKrige2)Error in

Re: [R-sig-Geo] How to use coKriging with LANDSAT ETM Imagery

2015-08-21 Thread Edzer Pebesma
With the amount of detail you give (no script) I can only guess. In case you try to feed the gstat routines a Raster or RasterStack object, it doesn't work. In that case you should first convert it; say if your raster is called x, pass it to gstat() as as(x, SpatialPointsDataFrame). On

[R-sig-Geo] How to use coKriging with LANDSAT ETM Imagery

2015-08-21 Thread Rajabhushanam. C
Hi, I am using a Landsat ETM Scene imagery for calculating the Linear Corregionalization Method in CoKriging. (gstat). Specifically, there is one raster band imagery and then using readGDAL and gstat function calls, I am making a gstat object that is input to the Variogram call. Now, while I