Блог

У меня есть 200 растров с температурами воздуха на каждый день, хранящихся в каталоге, и я хотел бы получить совокупную сумму температур за каждый день, например, day1 = day1; day2 = day1 day2; day3 = day1 day2 day3 и т.д. Я пытался сделать это в raster пакете, но cumsum функция вычисляет совокупную сумму всех ячеек в каждом растре, а не совокупную сумму отдельных растров (по крайней мере, мне так кажется из моих результатов). Я попытался сделать это так:

 library(raster)

setwd("C:/Air_temperatures/AT") # set working directory

AT.files <- list.files(pattern="AT") # list files with name AT
All.AT <- unique(AT.files) # unique files
  for (i in 1:200) {
    AT.i <- All.AT[i] # make i
    AT.raster.i<-raster(AT.i) # make rasters from files
    AT.sum.i <- cumsum(AT.raster.i) # ???calculate cumulative sum???
    writeRaster(AT.sum.i,
               paste(AT.i, "_cumsum.tif", sep = ""),
              datatype='FLT4S', overwrite=TRUE) # write new files and paste new name
  } 
 

Этот цикл работал, когда я пытался, например, добавить константу к каждому растру и так далее, Но я понятия не имею, как рассчитать совокупную сумму.

Reduce Здесь может помочь функция. Эта функция последовательно объединяет объекты на основе указанной функции. Поскольку существует метод для Raster* объектов, вы можете использовать его как функцию и получить совокупную сумму значений растра.

Сначала создайте 200 растровых объектов в списке:

 theATs <- lapply(All.AT, raster)
 

Затем используйте Reduce функцию с accumulate аргументом, равным TRUE

 res <- Reduce(" ", theATs, accumulate = TRUE)
 

Затем запишите результаты в файлы:

 lapply(seq_along(res), function(x) {
    writeRaster(res[[x]],
# Apparently, using .tif files threw an error
#     paste(All.AT[x], "_cumsum.tif", sep = ""),
# Alternatively, omit file extension to write default type (usually .grd)
      paste(All.AT[x], "_cumsum", sep = ""),
      datatype = 'FLT4S', overwrite = TRUE)
    })
 

Растры, которые вы создаете здесь, хранятся в памяти, поэтому, если исходные файлы очень большие, у вас могут возникнуть проблемы.