r - 合并不同范围的栅格，对R中重叠的单元格值求和

Question

我正在尝试合并具有不同范围的光栅化折线，以创建指示单元重叠次数的单个表面。

由于计算限制(考虑到我的研究区域的大小)，我无法使用 extend然后 stack对于每个栅格(总计数 = 67)。

我遇到过merge R 中的函数，这允许我将栅格合并到一个表面中。然而，它似乎不喜欢我插入一个函数来计算重叠单元格的总和。

也许我遗漏了一些明显的东西，或者这是 merge 的限制。功能。关于如何生成此输出的任何建议，避免 extend & stack将不胜感激!

代码:

# read in specific route rasters
raster_list <- list.files('Data/Raw/tracks/rasterized/', full.names = TRUE)

for(i in 1:length(raster_list)){

   # get file name
   file_name <- raster_list[i]

   # read raster in
   road_rast_i <- raster(file_name)

   if(i == 1){

   combined_raster <- road_rast_i

   } else {

   # merge rasters and calc overlap
   combined_raster <- merge(combined_raster, road_rast_i, 
                            fun = function(x, y){sum(x@data@values, y@data@values)})
  }
 }

当前输出图像:



单个路线的图像(示例):



修复图像:

Answer 1

解决了。有一个马赛克功能，它允许以下内容:

combined_raster <- mosaic(combined_raster, road_rast_i, fun = sum)