将栅格转换为矩阵

Posted 2021-03-28

我可以读取图像，栅格，将值限制在10-100。我不能做的是将限制转换为矩阵，我可以将所有值相加。

library(raster)

DEM <- raster("img.JPG")

image(DEM, zlim=c(10,100))

我想转换结果

image(DEM, zlim=c(10,60))

到我可以执行计算的矩阵中。

image(DEM, zlim=c(10,60))将导致

目标是仅对红色圆圈求和。

答案

library(raster)
d <- raster("img.JPG")
dd <- reclassify(d, rbind(c(-Inf, 10, NA), c(60, Inf, NA)))
# or:  d[d< 10 | d > 60] <- NA
plot(dd)
hist(dd)

或者，如果你真的想要一个矩阵

m <- matrix(d)
m[m<10] <- NA