рассчитать объем пересечения распределений использования - PullRequest
Новая
       46

рассчитать объем пересечения распределений использования

0 голосов
/ 19 июня 2020

Я пытаюсь рассчитать объем пересечения для распределения использования многих животных. Я рассчитал их UD, используя move::brownian.bridge.dyn() для каждого животного, затем использовал 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

Я также пытался создать объект estUDm, складывая растры UD.

Есть идеи? Я открыт для другого метода расчета пересечения объемов, просто не нашел другого метода.

Спасибо!

...