1

我正在尝试计算许多动物利用率分布的交集量。我已经估计了他们move::brownian.bridge.dyn()对每只动物使用的 UD，然后使用getVolumeUD()这给了我 UD 作为RasterLayer.

mod1 <- brownian.bridge.dyn(indiv, raster = xy, location.error = 10) #using default margin and window size
mod1.ud <- getVolumeUD(mod1)

为了计算我试图使用的重叠，adehabitatHR::kerneloverlapHR它要求 UD 是一个estUDm类型对象。我可以将单个 UD 变成一个estUD对象

mod1.px <- as(mod1.ud, "SpatialPixelsDataFrame")
mod1.estud <- new("estUD",mod1.px)

但无法弄清楚如何将它们合二为一estUDm。我试过了

all.ud <- c(mod1.estud, mod2.estud, ...)
class(all.ud) <- "estUDm"
vi <- kerneloverlapHR(all.ud, meth="VI", percent=95, conditional=TRUE)
Error in kerneloverlaphr(all.ud) : x should not be a volume under UD

我还尝试通过堆叠 UD 栅格来制作 estUDm 对象。

有任何想法吗？我也愿意用另一种方法来计算体积交点，只是还没有找到另一种方法。

谢谢！

4

0 回答 0