[R-sig-Geo] kriging: export to grid

Nicolas Meurisse nmeuriss at ulb.ac.be
Sun Dec 14 10:24:37 CET 2008


Hello,

After kriging with the use of the  krige.conv function, I would like to 
export my result under the grid format. In order to view it into a GIS.

It was suggested to me to use the writeGDAL function into the rgdal 
package. However It looks like I have a problem of supported formats
(dataset object of SpatialGridDataFrame-class or 
SpatialPixelsDataFrame-class ) that is not the same as the one resulting 
from my krige.conv function (kriging -class apparently)

I tried several function to achieve this conversion (getRasterData, 
SpatialGridDataFrame, SpatialPixelsDataFrame, create2GDAL, 
points2grid...). But without any success.

#my actual kriging script:
pred.grid = expand.grid(seq(170000, 130000, l = 50), seq(190000, 160000, 
l = 50))
myPred = krige.conv(myGData, loc = pred.grid, krige = 
krige.control(obj.m = myOLSMod))
image(myPred, loc = pred.grid, col = gray(seq(1, 0.1, l = 30)))

#my conversion attempts
tf <- tempfile()
writeGDAL(myPred, tf, drivername="GTiff", type="Byte", options=NULL)
[Erreur dans function (classes, fdef, mtable) : unable to find an 
inherited method for function "gridded", for signature "kriging"]

Thanks a lot,

Nicolas




More information about the R-sig-Geo mailing list