[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