adehabitathr - 使用 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 <- mcp((data)[,1],percent=95, unin = c("m"), unout = c("m2")) 计算了我的数据的主范围

它绘制的情节（cp，轴=真，边界=彩虹（12））

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

任何帮助，将不胜感激！！

谢谢

胡安

Answer 1

这是一个使用包中的示例数据的示例adehabitatHR，因为您并没有真正提供可重现的示例。

library(adehabitatHR)
data("puechabonsp")

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

一种方法是使用ggplot2and sf：

library(sf)
library(tidyverse)

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