Конвертировать растр в матрицу

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