我想在R中有一个Geotiff图像,该怎么办?以下是一个简单的示例。我想以Geotiff格式制作剧情。

require(gstat)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m)
spplot(x["var1.pred"], main = "ordinary kriging predictions")

最佳答案

将您的SpatialPixelsDataFrame转换为栅格并编写:

> require(raster)
> r = raster(x["var1.pred"])
> plot(r)
> writeRaster(r,"r.tiff","GTiff")


plot是可选的。

关于r - 如何获取geotiff图像,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/13910161/

10-12 23:29