我有在西澳大利亚追踪到的两种动物的后期/长期数据，我想用adehabitatHR找到它们的活动范围。
library(sp)
library(rgdal)
library(raster)
library(adehabitatHR)
library(sf)
quolls<-read.csv("quolls.csv")
head(quolls)
经纬度animal_ID1 -22.62271 117.12472 -22.62286 117.12463 -22.62192 117.12234 -22.62021 117.12245 -22.61989 117.12446 -22.62022 117.1260 1
但是每种动物的家园范围估计显然太小了。我想EPSG一定是错的，但经过很长时间的寻找，我仍然找不到正确的。
谁能给我指个正确的方向?# 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
动物1的起始范围为1.217428e-08，动物2的起始范围为6.253689e-08。
核密度估计也是如此;
quoll_ud <- adehabitatHR::kernelUD(quolls2[7],grid = 450)
quoll_hr <- adehabitatHR::getverticeshr(quoll_ud, 99)
print(quoll_hr)
估算动物1为2.36917592701502e-08，动物2为1.16018636413173e-07。
偶然发现答案…是EPSG 28350。
我最终放弃了原始的lats和long，而是使用st_read导入我拥有的动物数据的shapefile，从而使其工作。然后st_transform到28350。然后，由于mcp只接受SpatialPoints，我将对象转换为as(obj, "Spatial").