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

【问题标题】：How do I calculate the area of a Brownian bridge movement model?如何计算布朗桥运动模型的面积？
【发布时间】：2023-03-14 14:15:01
【问题描述】：

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

【解决方案1】：

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

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

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

【讨论】：