Quiero realizar una regresión de ventana móvil en cada píxel de dos pilas ráster que representan los datos de Band3 y Band4 de Landsat. El resultado debe ser dos pilas adicionales, una que representa la intersección y la otra que representa la pendiente de la regresión. Entonces, la capa 1 de la pila "B3" y la pila "B4" dan como resultado la capa 1 de la pila "interceptar" y la pila "pendiente". La capa 2 de la pila B3 y la pila B4 dan como resultado la capa 2, .... y así sucesivamente.

Ya vine por elgwrfunción, pero desea permanecer en el paquete ráster. De alguna manera sé quefocal debe incluirse para configurar mi ventana móvil (que debe ser de 3x3 píxeles) y de alguna manera un modelo lineal como:lm(as.matrix(b3)~as.matrix(b4)) aunque no creo que esto me consiga los valores en píxeles ...

En lugar de un rasterstack, también es posible un enfoque de capa por capa. (Por lo tanto, no necesariamente debe ser un rasterstack de Band3.

¿Alguien sabe cómo programar esto en R?