用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)

我想要18个Geotiff文件蒙面从我的shapefile，但输出只有一个tif文件，它不打开。我在网上使用搜索选项，但找不到合适的答案。

谢谢

Answer 1

你一定会找到你问题的答案的。它在here和here之前被问到，答案提供了here，但也提供了由writeRaster提供的bylayer选项的here。

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

就像这样

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