Chcę mieć obraz Geotiff w R. Jak mogę to zrobić? Oto prosty przykład. Chcę mieć wyprodukowaną fabułę w formacie Geotiff.jak uzyskać obraz geograficzny
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")
Konwersja SpatialPixelsDataFrame do rastra i napisać go:
> require(raster)
> r = raster(x["var1.pred"])
> plot(r)
> writeRaster(r,"r.tiff","GTiff")
plot jest opcjonalne.
Napisz SpatialPixelsDataFrame bezpośrednio do formatu rastrowego przy użyciu rgdal, który zapewnia powiązania z GDAL. Istnieje wiele rzeczy można zrobić z gdal, takie jak czytanie i/lub pisaniu szeroką gamę formatów rastrowych, typów pikselowych, zrobić transformacji współrzędnych, obsługiwać wiele zespołów, itp
require(rgdal)
writeGDAL(x["var1.pred"], "var1_pred.tif")
To nie ma wiele zmysłu. Geotiff to raster, twoje punkty to dane wektorowe. Co planujesz zrobić z Geotiffem, kiedy go masz? – Spacedman
Edytowałem kod. to, czego chcę, to w końcu mapa przewidywań. –