làm thế nào để có được hình ảnh địa lý

Question

Tôi muốn có một hình ảnh địa lý trong R.How Tôi có thể làm điều đó? Sau đây là một ví dụ đơn giản. Tôi muốn có cốt truyện được sản xuất ở định dạng 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

Chuyển đổi SpatialPixelsDataFrame của bạn đến một raster và viết nó:

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

Các plot là không bắt buộc.

Answer 2

Viết SpatialPixelsDataFrame trực tiếp vào định dạng raster sử dụng rgdal, cung cấp các liên kết đến GDAL. Có nhiều điều bạn có thể làm với GDAL, chẳng hạn như đọc và/hoặc viết một loạt các định dạng raster, loại pixel, thực hiện chuyển đổi tọa độ, hỗ trợ nhiều dải, v.v.

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