Как рассчитать площадь модели движения броуновского моста?

В пакете adehabitatHR есть пример кода для расчета модели движения броуновского моста для кабана с использованием функции kernelbb. Я хотел бы рассчитать площадь объекта, который получается. Я также хотел бы знать единицы площади. Я думаю, что NTF (Париж)/зона Ламберта III", код EPSG 27573, является CRS.

Вот код:

require(adehabitatHR)

data(puechcirc)
x <- puechcirc[1]
x

lik <- liker(x, sig2 = 58, rangesig1 = c(1, 10))

tata <- kernelbb(x, sig1 = 6.23, sig2 = 58, grid = 50)
tata

image(tata)
plot(getverticeshr(tata, 95), add=TRUE, lwd=2)

r gps kernel-density adehabitathr
person adkane    schedule 03.09.2019    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Функция getverticeshr возвращает объект класса SpatialPolygonsDataFrame. Вы можете использовать, например, gArea из пакета rgeos для расчета площади.

library(rgeos)
sp <- getverticeshr(tata, 95 )
gArea(sp)
## [1] 1576647

Я предполагаю, что единицы m^2. Чтобы проверить систему отсчета координат (которая не была установлена ​​в этом примере), вы должны сделать: proj4string(sp).

person johannes    schedule 03.12.2019