geotiff resmi nasıl edinilir

Question

R'de bir Geotiff görüntüsüne sahip olmak istiyorum Bunu nasıl yapabilirim? Aşağıda basit bir örnek. Geotiff formatında üretilen arsa sahip olmak istiyorum.

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

bir raster için SpatialPixelsDataFrame dönüştürme ve yazmak:

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

plot isteğe bağlıdır.

Answer 2

SpatialPixelsDataFrame ürününü GDAL numaralı bağlamalara bağlayan rgdal kullanarak bir tarama biçimine doğrudan yazın. Böyle okuma ve/veya raster formatları, piksel türleri geniş bir yelpazede yazma gibi GDAL ile yapabileceğiniz birçok şey, vb, koordinat dönüşümleri yapmak çoklu bantları destek vardır çok yapmaz

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