我可以读取图像,光栅,将值限制在10-100之间。我不能做的就是将限制转换为可以对所有值求和的矩阵。
library(raster)
DEM <- raster("img.JPG")
image(DEM, zlim=c(10,100))
我想转换
的结果image(DEM, zlim=c(10,60))
放入一个我可以进行计算的矩阵中。
目标是仅对红色圆圈求和。
答案 0 :(得分: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