将光栅转换为矩阵

Question

我可以读取图像，光栅，限制数值从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))会导致

​

目标是只与红圈相加。

​

​

Answer 1

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