7

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")
4

2 に答える 2

8

SpatialPixelsDataFrame をラスターに変換して書き込みます。

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

はplotオプションです。

于 2012-12-17T10:49:58.653 に答える
0

GDALrgdalへのバインドを提供するを使用して、SpatialPixelsDataFrame をラスター形式に直接書き込みます。GDAL でできることはたくさんあります。たとえば、さまざまなラスター形式、ピクセル タイプの読み取りと書き込み、座標変換の実行、複数バンドのサポートなどです。

require(rgdal)
writeGDAL(x["var1.pred"], "var1_pred.tif")
于 2013-10-03T05:12:33.960 に答える