How to use coKriging with LANDSAT ETM Imagery
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 08/21/2015 12:03 PM, Rajabhushanam. C 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. [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster, Heisenbergstra?e 2, 48149 M?nster, Germany; +49 251 83 33081 Journal of Statistical Software: http://www.jstatsoft.org/ Computers & Geosciences: http://elsevier.com/locate/cageo/ Spatial Statistics Society http://www.spatialstatistics.info -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150821/d96a52db/attachment.bin>