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

Question

在adehabitatHR包中，有一些使用kernelbb函数计算野猪的布朗桥运动模型的示例代码。我想计算得到的物体的面积。我还想知道这个地区的单位。我认为NTF (巴黎)/ 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)

Answer 1

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

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

我假设单位是m^2。要检查坐标参考系(在本例中没有设置)，您可以这样做：proj4string(sp)。