Регрессия движущегося окна

Я хочу выполнить регрессию движущегося окна для каждого пикселя двух растровых стеков, представляющих данные Band3 и Band4 Landsat. Результатом должно быть два дополнительных стека, один представляющий Перехват, а другой представляющий наклон регрессии. Таким образом, уровень 1 стека "B3" и стека "B4" приводит к уровню 1 стека "intercept" и стека "slope". Уровень 2 стека B3 и стека B4 приводит к уровню 2, .... и так далее.

Я уже пришел вместеgwrфункции, но хотите остаться в растровом пакете. Я как-то знаю чтоfocal должен быть включен для установки моего движущегося окна (которое должно быть 3х3 пикселя) и как-то линейной модели, например:lm(as.matrix(b3)~as.matrix(b4)) хотя я не думаю, что это дает мне значения по пикселям ...

Вместо растрового стека также возможен послойный подход. (Так что это не обязательно должен быть растерстак Band3.

У кого-нибудь есть клей, как программировать это на R?

Ответы на вопрос(1)

Ваш ответ на вопрос