Регрессия движущегося окна
Я хочу выполнить регрессию движущегося окна для каждого пикселя двух растровых стеков, представляющих данные 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?