用 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((数据)[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"))

一种方法是使用 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")