Mathieu Basille
Thu, 02 Sep 2010 07:09:47 -0700
Tom,
Did you have a look at 'area2spol'?
From the example of 'kernelUD':
data(puechabon)
loc <- puechabon$locs[, c("X", "Y")]
id <- puechabon$locs[, "Name"]
(ud <- kernelUD(loc, id))
ver <- getverticeshr(ud, 95)
Then, you can do something like that:
spver <- lapply(ver, area2spol)
Which returns a list of SpatialPolygons (sp):
class(spver[[1]])
[1] "SpatialPolygons"
attr(,"package")
[1] "sp"
I don't really know spatstat but there might be some conversion tools
between sp and spatsat.
Hope this helps, Mathieu. Le 02/09/2010 05:58, Tom_R a écrit :
Hi List, I am trying to calculate the area of overlap amongst a list of animal home ranges; ten ants each, from five colonies. The help page for kernelUD {adehabitat} states, " getverticeshr stores the home range contour as objects of class area in a list of class kver, with one component per animal." [the 'components' are of class 'area'] . However, I am having difficulties converting the list (class 'kver') of home ranges (class 'area') to a list of polygons that can be manipulated with spatstat. xy<- data.frame(cbind(x[colony==i& ant==j], y[colony==i& ant==j])) who<- as.factor(ant[colony==i& ant==j]) ud<- kernelUD(xy, who, h = "href", same4all = TRUE) ver<- getverticeshr(ud, 50) class(ver)[1] "kver"After converting the 'ver' list to polygons, I would convert to owin objects (in the spatstat package) and then use the functions intersect.owin and area.owin to find the area of the overlap. In summary: 1) How to convert a list of class kver, with one area object per animal? Many thanks in advance!! Tom Richardson
-- ~$ whoami Mathieu Basille, Post-Doc ~$ locate Laboratoire d'Écologie Comportementale et de Conservation de la Faune + Centre d'Étude de la Forêt Département de Biologie Université Laval, Québec ~$ info http://ase-research.org/basille ~$ fortune ``If you can't win by reason, go for volume.'' Calvin, by Bill Watterson. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo