Как я могу написать функцию для гибкой многодневной комбинации спутниковых снимков? - PullRequest
Новая
       37

Как я могу написать функцию для гибкой многодневной комбинации спутниковых снимков?

0 голосов
/ 20 января 2020

Я пытаюсь написать функцию, которая будет выполнять роль временного фильтра для списков растровых данных. Цель этого - заполнить пиксели NA значениями за день до или после этого дня. Я хочу объединить 3-5 дней изображений в одно с меньшим облачным покровом. Я просто не уверен, как отформатировать функцию в l oop и установить порог для остановки дальнейших комбинаций, когда пиксели NA составляют менее 30% изображения.

Формат растров - это списки с 1:30 или 1:31: 

[[30]]
class      : RasterLayer 
dimensions : 931, 1541, 1434671  (nrow, ncol, ncell)
resolution : 0.004166667, 0.004166667  (x, y)
extent     : -122.6792, -116.2583, 49, 52.87917  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : memory
names      : layer 
values     : 0, 239  (min, max)

Функция, которую я написал ранее для объединения трех спутниковых изображений в один и тот же день, была выполнена с использованием функция покрытия также. Я не знаю, поможет ли это: 

Replace <- function(VIIRS, AQUA, TERRA){
  month_combined <- list()
  for(i in 1:length(VIIRS)){
      month_combined[[i]] <- cover(TERRA[[i]], VIIRS[[i]], AQUA[[i]])
  }
  return(month_combined)
}

, который охватывает пиксели NA (облако, отсутствующие данные и т. Д. c.) С пикселем другого спутника со значением пикселя в тот же день , Мне нужно написать функцию, которая будет делать это, но со значениями за день до или после.

Заранее спасибо!

...