adehabitatHR 家庭范围估计太小答案

【问题标题】：adehabitatHR home range estimation is too smalladehabitatHR 家庭范围估计太小
【发布时间】：2021-10-11 00:36:34
【问题描述】：

我有在西澳大利亚追踪的两只动物的经纬度数据，我想使用 adehabitatHR 找到它们的栖息地。

library(sp)
library(rgdal)
library(raster)
library(adehabitatHR)
library(sf)

quolls<-read.csv("quolls.csv")
head(quolls)

纬度经度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

但每只动物的家园范围估计显然太小了。我认为 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 adehabitathr

【解决方案1】：

刚刚偶然发现答案.. 它是 EPSG 28350。

最终我放弃了原始的经纬度，而是使用 st_read 导入了我拥有的动物数据的 shapefile。然后 st_transform 到 28350。然后由于 mcp 只接受 SpatialPoints，我将对象转换为 as(obj, "Spatial").

【讨论】：