2012-12-17 2 views
6

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") 
+0

의미의. 지오 티프는 래스터이고 점은 벡터 데이터입니다. Geotiff를 가지고있을 때 무엇을 할 계획입니까? – Spacedman

+0

코드를 편집했습니다. 내가 원하는 것은 예측지도입니다. –

답변

7

는 래스터에 SpatialPixelsDataFrame 변환하고 쓰기 :

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

plot는 선택 사항입니다.

0

GDAL에 대한 바인딩을 제공하는 rgdal을 사용하여 SpatialPixelsDataFrame을 직접 래스터 형식으로 씁니다. 이러한 읽기 및/또는 래스터 형식, 픽셀 다양한 유형의 쓰기와 같이 GDAL과 함께 할 수있는 많은 일들이, 등, 좌표 변환을 할 다중 대역을 지원하고 있습니다 많이하지 않습니다

require(rgdal) 
writeGDAL(x["var1.pred"], "var1_pred.tif") 
관련 문제