在栅格上计算R中的像素级线性回归

时间:2020-09-08 13:20:16

标签: linear-regression

ı想为每个网格单元计算线性回归。首先,ı试图将光栅文件转换为数组以使其更容易

例如,有三个栅格数据,并且ı尝试将这些栅格转换为1个数组,并为每个网格单元计算线性回归,该怎么做?

x <- raster(ncol=94, nrow=192, xmn=-180, xmx=180, ymn=-90, ymx=90)
x1 <- raster(ncol=94, nrow=192, xmn=-180, xmx=180, ymn=-90, ymx=90)
x2 <- raster(ncol=94, nrow=192, xmn=-180, xmx=180, ymn=-90, ymx=90)
x3 <- raster(ncol=94, nrow=192, xmn=-180, xmx=180, ymn=-90, ymx=90)
values(x)  <- 1:ncell(x)
values(x1) <- 1:ncell(x)
values(x2) <- 1:ncell(x)
values(x3) <- 1:ncell(x)
set.seed(0)
values(x)  <- runif(ncell(x1))
values(x1) <- runif(ncell(x1))
values(x2) <- runif(ncell(x2))
values(x3) <- runif(ncell(x3))

0 个答案:

没有答案