r - 如何获取geotiff图像

Question

我想在 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")

Answer 1

将您的 SpatialPixelsDataFrame 转换为栅格并写入：

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

plot是可选的。

Answer 2

使用 将 SpatialPixelsDataFrame 直接写入光栅格式rgdal，它提供与GDAL的绑定。您可以使用 GDAL 做很多事情，例如读取和/或写入各种栅格格式、像素类型、进行坐标转换、支持多波段等。

require(rgdal)
writeGDAL(x["var1.pred"], "var1_pred.tif")