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

[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 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>