我正在尝试在一个png或其他光栅文件中合并(叠加)〜1000张(共12MB)透明图像* .png。必须使用R中的哪个程序包或函数?
图像包含雷达云数据。埃弗里·戴夫(Evrey)那天,雷达使〜350 png。 样本:001,002,003,004
这里是legend颜色的传说
然后,我将结果与120个雨量计的数据进行比较。
现在我使用此代码
library(base)
library(png)
homeserverDir = "S:/R/radarTEST"
dDir_new = paste(homeserverDir, "dd_stack", sep = "/")
rr_dir = paste(homeserverDir, "dd_stackD", sep = "/")
setwd(dDir_new)
filenames <- list.files(path = dDir_new, pattern="*.png")
theATs <- lapply(filenames, raster)
STACK1 <- stack(theATs)
res <- Reduce("+", theATs, accumulate = TRUE)
setwd(rr_dir)
lapply(seq_along(res), function(x) {
writeRaster(res[[x]], paste(filenames[x], "_cumsum", sep = ""), datatype = 'FLT4S', overwrite = TRUE)
})
Grids_list = list.files(rr_dir, "grd")
Radar_sum = raster(Grids_list[length(Grids_list)])
plot(Radar_sum)