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

Question

在包中有一些示例代码，用于使用该函数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)

Answer 1

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

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

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