如何用R的leaflet包从多CSV生成热力图？栅格合并/重采样问题

问题背景 #

• 基于GPS位置及对应属性生成热力图，现有约400个包含GPS位置的CSV文件，单文件可正常生成栅格热力图
• 多文件批量处理时，栅格合并遇到两类问题：
  • 直接合并：因分辨率/原点不匹配报错
  • 按原栅格（KernelDensityRaster）的res/ext/crs参数重采样后，栅格值全为NA

现有代码 #

合并代码 #

raster_layers<-lapply(processed_data_list, function(item) item$raster_data)

merged_raster_stack<- do.call(merge,raster_layers)

重采样尝试代码 #

resampled_raster<- KernelDensityRaster # 导致合并错误
# or
# resampled_raster<- raster(res = res(KernelDensityRaster), ext = extent(KernelDensityRaster), crs = crs(KernelDensityRaster)) # 生成全NA值的栅格

栅格数值统计 #

原栅格（KernelDensityRaster）统计 #

print(summary(values(KernelDensityRaster)))
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00000 0.00000 0.04892 6.14610 9.73696 56.01518

重采样后栅格统计 #

print(summary(values(resampled_raster)))
Mode NA's
logical 10000

需求 #

寻求解决多CSV栅格合并与重采样问题的可行方法。

内容的提问来源于stack exchange，提问作者andyroberts1234