13
Вот пример.
library(raster)
# example data
x <- raster(system.file("external/test.grd", package="raster"))
Чтобы получить прямоугольную степень
e <- extent(x)
# coerce to a SpatialPolygons object
p <- as(e, 'SpatialPolygons')
Чтобы получить многоугольник, который окружает ячейки, которые не являются NA
# make all values the same. Either do
r <- x > -Inf
# or alternatively
# r <- reclassify(x, cbind(-Inf, Inf, 1))
# convert to polygons (you need to have package 'rgeos' installed for this to work)
pp <- rasterToPolygons(r, dissolve=TRUE)
# look at the results
plot(x)
plot(p, lwd=5, border='red', add=TRUE)
plot(pp, lwd=3, border='blue', add=TRUE)
r <- r > -Inf
части. Что именно делает? И насколько он отличается отvalues(r)[!is.na(values(r))] <- 1
(который устанавливает все местоположения, которые не являются NA в r в 1).r > -Inf
является основнымR
. Пойдуc(1,NA,3,NA) > -Inf
посмотреть, как это работает. Я добавил альтернативу (reclassify
). Ваша альтернатива работает, но она не подходит для больших объектов.rasterToPolygons()
деталь, и тогда я оценю ваш ответ. Спасибо за-Inf
бит, это будет очень полезно!aggregate
поскольку вы все равно не увидите разницу.aggregate()
он не обрабатывается (компьютер оснащен 16 ГБ ОЗУ). Проблема в том, что мне нужно сохранить границу как можно более родной, поскольку матрицу высот необходимо извлечь из результирующего многоугольника. Если я продолжу понижать частоту моего растра, я потеряю эту границу. Любая работа вокруг?