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

[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 model.frame.default(terms(formula), as(data, "data.frame"), na.action = na.fail) :
  object is not a matrix


> coKrige.fit <- fit.lmc(coKrige1, coKrige2, vgm(1, "Gau", 350, 1), TRUE, 1.0)Error in fit.lmc(coKrige1, coKrige2, vgm(1, "Gau", 350, 1), TRUE, 1) :
  gstatVariogram var1 not present


> coKrige.predict <- predict(coKrige.fit, rGrid)Error in UseMethod("predict") :
  no applicable method for 'predict' applied to an object of class
"c('variogramModel', 'data.frame')"




On Fri, Aug 21, 2015 at 3:33 PM, Rajabhushanam. C <rajcheruk at gmail.com>
wrote:

> 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 use *fit.lmc()*, I get an error message *"gstatVariogram
> var1 not present"*. I have done some preliminary findings, and there are
> no pertinent messages in this message board that have Raster as an input to
> cokriging.
>
> Really appreciate if anybody can point me some relevant code or pseudo
> code to perform cokriging with Landsat one band imagery.
>
> Thank You
> C. Rajabhushanam.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150830/89cb88a4/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Sample1_Band3.tif
Type: image/tiff
Size: 20384 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150830/89cb88a4/attachment.tif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Image_Fit_coKriging.R
Type: text/x-r-source
Size: 1008 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150830/89cb88a4/attachment.bin>