I have lat/long data of two animals tracked in Western Australia and I'd like to find their home ranges using adehabitatHR.
library(sp)
library(rgdal)
library(raster)
library(adehabitatHR)
library(sf)
quolls<-read.csv("quolls.csv")
head(quolls)
Latitude Longitude animal_ID 1 -22.62271 117.1247 1 2 -22.62286 117.1246 1 3 -22.62192 117.1223 1 4 -22.62021 117.1224 1 5 -22.61989 117.1244 1 6 -22.62022 117.1260 1
But the home range estimates of each animal are obviously too small. I think the EPSG must be wrong but after a very long time looking I still can't find the right one.
Can anyone point me in the right direction please?
# make a SpatialPoints dataframe without a CRS
quolls2 <- quolls
quoll.latlong<-data.frame(x=quolls2$Longitude,y=quolls2$Latitude)
coordinates(quolls2) <- quoll.latlong
# add crs
proj4string(quolls2) <- CRS(SRS_string = "EPSG:4283")
mcp<-mcp(quolls2[,7],percent=95,unout = c("ha"))
mcp
Home range for animal 1 is 1.217428e-08 and animal 2 is 6.253689e-08.
And likewise with kernel density estimation;
quoll_ud <- adehabitatHR::kernelUD(quolls2[7],grid = 450)
quoll_hr <- adehabitatHR::getverticeshr(quoll_ud, 99)
print(quoll_hr)
which estimates animal 1 at 2.36917592701502e-08 and animal 2 at 1.16018636413173e-07.