Я пытаюсь объединить (наложить) ~ 1000 (всего 12 МБ) прозрачных изображений * .png в один PNG или другой растровый файл.Какой пакет или функцию в R необходимо использовать?
Изображения содержат данные облаков радара.Каждый день радар составляет ~ 350 png.Образец: 001 , 002 , 003 , 004
Вот легенда цветов легенда
После этого я сравню результат с данными 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)