Растр в точки с пропущенными значениями

Question

Я использую функцию rasterToPoints() из пакета raster.По умолчанию эта функция пропускает значения NA.

raster <- raster(matrix(c(1,2,NA, 1,2,3, 1,2,3), nrow = 3, ncol = 3 ))
raster_points <- rasterToPoints(raster, na.rm = FALSE)

Поэтому я помещаю растр с 200 ячейками и получаю data.frame только с 150 строками.Как я могу получить точку с пустыми ячейками в моем конечном data.frame?

Answer 1

Вы можете объединить xyFromCell, values и cbind. Это возвращает центральные координаты каждой ячейки (когда вы инициализировали растр, он принимает экстент 1 x 1, если не определено иначе) и его значение, включая NAs, в матричной форме:

cbind(xyFromCell(raster, 1:ncell(raster)), values(raster))