如何计算布朗桥运动模型的面积？

Question

在 adehabitatHR 包中有一些示例代码使用 kernelbb 函数计算野猪的布朗桥运动模型。我想计算结果对象的面积。我还想知道该地区的单位。我认为 NTF (Paris) / Lambert zone III", EPSG code 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)

Answer 1

函数 getverticeshr returns class SpatialPolygonsDataFrame 的一个对象。例如，您可以使用 rgeos 包中的 gArea 来计算面积。

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

我假设单位是m^2。要检查坐标参考系统（本例中未设置），您可以执行以下操作：proj4string(sp).

如何计算布朗桥运动模型的面积？

How do I calculate the area of a Brownian bridge movement model?

r

gps

kernel-density

adehabitathr