Ich möchte für jedes Pixel von zwei Raster-Stapeln, die Band3- und Band4-Landsat-Daten darstellen, eine Regression mit sich bewegenden Fenstern ausführen. Das Ergebnis sollten zwei zusätzliche Stapel sein, von denen einer den Achsenabschnitt und der andere die Steigung der Regression darstellt. Die Schicht 1 des Stapels "B3" und des Stapels "B4" führt zu der Schicht 1 des Stapels "Intercept" und des Stapels "Slope". Schicht 2 von Stapel B3 und Stapel B4 führt zu Schicht 2, ... und so weiter.

Ich bin schon am @ vorbei gekommgwrFunktion, möchte aber im Rasterpaket bleiben. Ich weiß irgendwie, dassfocal muss enthalten sein, um mein bewegliches Fenster (das 3x3 Pixel groß sein sollte) und irgendwie ein lineares Modell wie das folgende festzulegen:lm(as.matrix(b3)~as.matrix(b4)) obwohl ich nicht denke, dass dies mir die Werte pixelweise bringt ...

nstelle eines Rasterstacks ist auch ein Layer-by-Layer-Ansatz möglich. (Es muss also nicht unbedingt ein Rasterstack von Band3 sein.

Hat jemand eine Ahnung, wie man das in R programmiert?