Quero executar uma regressão da janela em movimento em cada pixel de duas pilhas de varredura que representam os dados da Banda3 e da Banda4 dos Landsat. O resultado deve ser duas pilhas adicionais, uma representando a interceptação e a outra representando a inclinação da regressão. Portanto, a camada 1 da pilha "B3" e a pilha "B4" resultam na camada 1 da pilha "interceptar" e "inclinação" da pilha. A camada 2 da pilha B3 e a pilha B4 resultam na camada 2, .... e assim por diante.

Eu já vim juntogwrfunção, mas deseja permanecer no pacote raster. De alguma forma eu sei quefocal deve ser incluído para definir minha janela móvel (que deve ter 3x3 pixels) e, de alguma forma, um modelo linear como:lm(as.matrix(b3)~as.matrix(b4)) embora eu não ache que isso me dê os valores em pixels ...

Em vez de um rasterstack, também é possível uma abordagem camada por camada. (Portanto, não deve ser necessariamente um rasterstack do Band3.

Alguém tem uma cola como programar isso em R?