r - 在 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 上找到您的问题的答案。之前在这里和这里被问过，这里提供了答案，但这里也有提供的选项。bylayerwriteRaster

bylayer = T如果您不需要将它们取消堆叠到列表中，请尝试。

就像是

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