Hello,
Is there a straightforward method for exporting a kernel home range
utilization distribution produced by kernelUD to an ESRI shapefile? I
have the following code:
xy<-Turtles[,c("x","y")]
id<-Turtles$Beach
kud<-kernelUD(xy,id)
I get a nice picture of the kernel home range utilization distribution
with contours. I want to create my maps in ArcMap (I'm not that good
with maps in R yet), so I'd like to export each homerange UD as an
ESRI shapefile, keeping the contour information so I can color ramp it
appropriately and do some analysis in ArcMap. This is what I think
I'm supposed to do, from reading
https://wiki.faunalia.it/dokuwiki/doku.php?id=public:animove_howto:
kver<-getverticeshr(kud)
spol <- kver2spol(kver)
df <- data.frame(Beach="Baja",row.names=("Agua Blanca, Baja California Sur"))
spdfBaja <- SpatialPolygonsDataFrame(spol, df, match.ID = TRUE)
writeOGR(spdfBaja,paste(td,"R-Output/Shapefiles",sep=""),"bajaUD","ESRI
Shapefile")
But I think this only gets me one contour level (95%), right? Once I
get that into ArcMap it's just a single polygon. It doesn't look like
getverticieshr takes a list of levels, just one. Do I need to repeat
this and glue them together for all the levels? Is there a different /
better way to do this?
Also, somewhat related, I'd like to get my projection information into
the SpatialPolygon object, but it's created in the kver2spol
function... is there a way to set the CRS still?
Thanks for any help or pointers.
Corrie
--
Corrie Curtice
Research Associate
Marine Geospatial Ecology Lab
Nicholas School of the Environment, Duke University
http://mgel.env.duke.edu
em: [email protected]
_______________________________________________
AniMov mailing list
[email protected]
http://lists.faunalia.it/cgi-bin/mailman/listinfo/animov
