Маскирповка облаков

[_taras_]

Доброго времени!
При расчете среднего снимка за период (r.series) возникает проблема с пикселями "нет значений", которые образовались после удаления облаков. Подскажите чем эти пустые места лучше заполнить дабы получить целый без пропусков растр?

[gamm]

Их лучше пропустить при подсчете среднего. Соответствующие функции есть и в R и в Python

[_taras_]

gamm, спасибо за ответ!
Поскольку мне надобно получить цельное, без особых пропусков, усредненную картинку (данные) за месяц не будет ли ересью следующий подход к маскировке проплешин от облаков.
1 Маской облаков условно 15 дня вырезаем участки от более ранних и поздних снимков.
2 Для них считаем среднее и объединяем с 15 днем
3 Так для всех снимков месяца.

[gamm]

чего-то мудрено ... обычно просто грузим очередной растр, считаем маску 0/1 (облака/не облака), и прибавляем снимок, умноженный на маску, в одно место, а маску прибавляем в другое. Потом делим сумму снимков на сумму масок, там, где сумма масок равна нулю, ставим NoData (не было ни одного безоблачного дня для этого пикселя). Место под сумму снимков и сумму масок заготавливаем при загрузке первого снимка (предполагается, что все снимки одинаковой геометрии). Если снимки разной геометрии, то после загрузки трансформируем на нужны растр (приводим к единой геометрии)..