Hello Everyone...
I am fairly new to R. I would like to use the adehabitatHR package to create
kernel density and isopleths from my sea turtle GPS data. I?m running into some
issues?
Basically I am having trouble assigning IDs and XY fields within R to create KDs
and MCPs.
Thank you in advance for any help!
Gwen Lockhart
GIS Research Specialist
Virginia Aquarium and Marine Science Center Foundation
Division of Research and Conservation
glockhar at virginiaaquarium.com<mailto:glockhar at virginiaaquarium.com>
Desk: 757-385-6486
Cell: 757-748-6757
#loading packages
`library(adehabitatHR)
library(raster)
library(rgdal)
library(maptools)`
#read CSV with UTM xy and ids
track<-
read.table("T:/GIS/Data/Tracking/state_space_model/20150210_AbHabiatatHRmodel/All_prelim_model_6hr_utm_forage_10A.csv",
header=TRUE, sep=",", na.strings="NA", dec=".",
strip.white=TRUE)
# Turn track into a SpatialPointsDataFrame sby specifying that the "X"
and "Y" columns are the coordinates:
coordinates(track) <- c("X", "Y")
class(track)
plot(track)
#Project into utm:
proj4string(track) <- CRS("+init=epsg:32618")
#read shore;ine file
shore <-
readShapeSpatial("C:/Users/gemme001/Desktop/R_state_space/STATES_VA_COAST_UTM.shp",
delete_null_obj=TRUE)
plot(shore)
proj4string(shore) <- CRS("+init=epsg:32618")
#Add to list the list with the names "map" and "relocs"
my.homerange.data <- list(map = shore, relocs = track)
#####THIS IS WHERE I AM RUNNING INTO ISSUES#####
#Assign IDs and XY ? the IDs work but coordinates don?t work. I get the
following error: Error in `[.data.frame`(x at data, i, j, ..., drop = FALSE) :
undefined columns selected
id<-my.homerange.data$relocs$Name
xy<-(my.homerange.data$relocs["X","Y"]
#Create CP ? this works
cp <- mcp(my.homerange.data$relocs[,1], percent=95)
class(cp)
plot(cp)
#Create KUD. This doesn?t work. I get the following: Error in xy.coords(x, y,
xlabel, ylabel, log) : 'x' is a list, but does not have components
'x' and 'y'
kud <- kernelUD(track[,1], h="href")
[[alternative HTML version deleted]]