How to use coKriging with LANDSAT ETM Imagery
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>