0

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

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

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

纬度 经度 动物_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.12424 1 6 -212.620 117.

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

4

1 回答 1

0

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

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

于 2021-08-06T08:54:31.663 回答