r语言 - 如何计算布朗桥运动模型的面积?

Question

在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)

Answer 1

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

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

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