如何獲取geotiff圖像

我想在R中擁有Geotiff圖像。我可以這樣做嗎？以下是一個簡單的例子。我想要以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")

來源

2012-12-17 Mona Eshraghi

這並不會產生很多的意義。 Geotiff是一個柵格，你的點是矢量數據。你有什麼打算與Geotiff做？ – Spacedman

我編輯了代碼。我最後想要的是一張預測圖。 –

將您SpatialPixelsDataFrame爲柵格，寫它：

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

的plot是可選的。

來源

2012-12-17 10:49:58 Spacedman

使用rgdal直接將SpatialPixelsDataFrame寫入光柵格式，該格式提供了與GDAL的綁定。使用GDAL可以做很多事情，例如讀取和/或寫入各種光柵格式，像素類型，進行座標變換，支持多個波段等。

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

來源

2013-10-03 05:12:33

如何獲取geotiff圖像

回答

相關問題