在R中编写多个GeoTiff文件

Question

我是新用户并在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和here之前询问了问题，here提供了答案，bylayer提供了writeRaster选项here。

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

类似

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