与Adehabitat一起策划个人家庭

Question

我试图将我的研究人员的名字放在多边形家庭范围内，但经过多次尝试后我仍然无法实现。

这里和我的数据示例：X和Y是坐标，id是个人

X   Y   id          
29  29  4           
44  28  7           
57  57  5           
60  81  11          
32  41  4           
43  29  7           
57  57  5           
46  83  11          
32  41  4           
43  29  7           
57  56  5           
60  82  11          
35  40  4           
43  28  7           
62  55  5           
54  73  11          
27  40  4           
43  28  7           
61  54  5

首先，我用MPC计算了我的数据的原始范围 cp＆lt; - mcp（（data）[，1]，percent = 95，unin = c（“m”），unout = c（“m2”））

情节吧 plot（cp，axes = TRUE，border = rainbow（12））

但是我不知道哪些多边形对应于每个人，如果可能的话，我需要在每个多边形中包含我个体的id

任何帮助都将不胜感激!!

由于

涓

Answer 1

以下是使用adehabitatHR包中的示例数据的示例，因为您并未真正提供可重现的示例。

library(adehabitatHR)
data("puechabonsp")

cp <- mcp(puechabonsp$relocs[, 1], percent=95, unin = c("m"), unout = c("m2"))

一种方法是使用ggplot2和sf：

library(sf)
library(tidyverse)

st_as_sf(cp) %>% ggplot(., aes(fill = id)) + geom_sf(alpha = 0.5) +
  scale_fill_discrete(name = "Animal id")