在adehabitatHR中，我的数据应该采用何种格式进行归属范围分析？

Question

我有南非跟踪动物的latlong数据，我使用adehabitatHR进行分析。这是我的数据样本：

Latitude        Longitude
-25.870265      27.947412
-25.816235      28.022442
-25.751107      28.1113
-25.670537      28.185403
-25.619823      28.290013

我需要转换数据，以便确定家庭范围。我一直在使用本教程来帮助我，他的数据看起来也像十进制度http://www.mikemeredith.net/blog/1212_Data_for_home_range_analysis_in_R.htm 这是我的代码：

library(raster) 
library(rgdal)
library(maptools)
library(adehabitatHR)
library(sp)

data <- read.csv("stackoverflowEg.csv", sep = ",", header = T)
head(data)
coordinates(data) <- c("X", "Y")

proj4string(data) <- CRS("+init=epsg:4326")

## I got the projected coordinate system from here [http://epsg.io/22235][1]

track <- spTransform(data, CRS("+proj=longlat +init=epsg:22235"))

summary(track)

cp <- mcp(track, percent=95)

cp

问题是mcp的结果值太小（2.230023e-07），所以我觉得我的预测做错了。任何帮助将非常感激。

Answer 1

试试这个

coordinates(df) <- ~Longitude+Latitude
proj4string(df) <- CRS('+init=epsg:4326')
dfutm <- spTransform(df, CRS('+init=epsg:32735')) # utm 35S for SAfrica
dfutm
SpatialPoints:
     Longitude Latitude
[1,]  594921.4  7138341
[2,]  602485.7  7144268
[3,]  611454.0  7151409
[4,]  618966.7  7160268
[5,]  629521.1  7165787
Coordinate Reference System (CRS) arguments: +init=epsg:32735 +proj=utm +zone=35 +south
+datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0