Столкнулся с задачей, которую я решал в R
Есть площадка 100 х 100 м. В ней на равном расстоянии друг от друга (10 м) стоят ловушки, образуя сетку 10 на 10 ловушек (итого 100 ловушек) - получается сетка как в морском бое. На это площадке отлавливали меченных зверьков.
Мне нужно было посчитать площадь индивидуального участка каждого зверька, это я сделал с помощью пакета adehabitat. Рассчитал MCP (минимальный выпуклый полигон). Я только начинающий пользователь, поэтому мой скрипт очень неловкий:
Код: Выделить всё
#ID=17
x17 <-c(40,40,70,10,50,20,40,40)
y17 <-c(80,100,100,90,100,90,10,80)
xy17 <-cbind(x17,y17)
xysp17 <- SpatialPoints(xy17)
clu17 <- clusthr(xysp17)
class(clu17)
plot(clu17)
sp17<-mcp(xysp17, percent=95, unin = c("m"),
unout = c("m2"))
plot(sp17)
sp17a<-mcp(xysp17, percent=100, unin = c("m"),
unout = c("m2"))
plot(sp17a)
print (sp17)
print (sp17a)
Этот скрипт я повторял вручную для каждого зверька. Наверное это выглядит также, как если бы я делал таблицу в фотошопе – сделать можно, но очень неудобно))
Теперь передо мной стоит задача визуализации результатов – мне нужно нанести на сетку ловушек полигоны индивидуальных участков зверьков. Нескольких зверьков одновременно (например всех зверьков одного пола, или всех зверьков одного вида).
У меня такой вопрос:
1) как построить свою сетку координат – она без привязки к географическим координатам, просто квадраты, как в морском бое.
2) как на этой сетке построить полигоны нескольких участков обитания одновременно
3) какой пакет лучше для этого использовать
Заранее благодарен
