在 R 中编写多个 GeoTiff 文件

Question

我是 R 新用户，在 R 中读取了 18 个 Geotiff 文件，并根据我感兴趣的区域屏蔽了 tiff。我现在无法编写 18 掩码的 Geotiff 文件。

require(raster)
raster_data <- list.files(path=getwd())
s <- stack(raster_data)
spf<-readShapePoly("basin.shp")
rc<-crop(s, extent(spf))
rm<-mask(rc, spf)
rf <- writeRaster(rm, filename=outputFile, overwrite=TRUE)

我想从我的 shapefile 中屏蔽 18 个 Geotiff 文件，但输出只有一个 tif 文件，而且它打不开。我在互联网上使用搜索选项，但找不到适合我查询的答案。

谢谢

Answer 1

您一定会在 SO 上找到问题的答案。在 here and here and answers provided here but also here 之前使用 writeRaster 提供的 bylayer 选项被询问过。

如果您不需要将它们拆分到列表中，请尝试 bylayer = T。

类似

library(raster)
r1 <- raster(ncol=10, nrow=10)
r1[] <- 1:100
s <- stack(r1, r1)
writeRaster(s, '~:/r.tif', bylayer = T)

在 R 中编写多个 GeoTiff 文件

Writing Multiple GeoTiff File in R

r

raster

geotiff