使用amt包时步长sl_未以米显示，是否需转换为UTM坐标？

问题原因 #

你猜测的完全正确，问题根源在坐标系（CRS）上。你使用的WGS84（EPSG:4326）是地理坐标系，单位为度，所以steps()计算出的sl_是经纬度的差值（度数），而非实际距离的米数。

解决方案：转换为UTM投影 #

UTM（通用横轴墨卡托投影）是平面坐标系，单位为米，是计算实际移动距离的理想选择。只需先将轨迹转换到数据所在区域的UTM坐标系，再计算步长即可。

具体操作步骤 #

1. 确定UTM分区：你的数据位于东经7°、北纬50°附近，对应UTM分区为32N（EPSG:32632）。如果数据在其他区域，可使用amt::utm_zone()自动匹配对应UTM代码。
2. 转换坐标系后重新计算步长。

修改后的代码示例 #

library(amt)
library(tidyverse)
library(lubridate)

# 加载你提供的示例数据
data <- structure(list(ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), 
    E = c(7.3900622244, 7.3919191098, 7.3993466517, 7.4051295235, 
    7.4166952672, 7.4190826913, 7.418817422, 7.4220537081, 7.4201437687, 
    7.4050234158), N = c(50.3376469514, 50.3356839582, 50.3288400089, 
    50.323905999, 50.3209349823, 50.322261329, 50.3302724634, 
    50.3436420387, 50.3467191632, 50.3597173613), datetime = c("2022-11-11 07:00:05", 
    "2022-11-11 08:00:12", "2022-11-11 11:02:05", "2022-11-11 12:00:35", 
    "2022-11-11 13:01:07", "2022-11-11 14:00:24", "2022-11-11 16:00:04", 
    "2022-11-11 17:00:12", "2022-11-11 22:01:01", "2022-11-11 23:00:44"
    )), class = "data.frame", row.names = c(NA, -10L))

# 标准化日期时间格式（用lubridate更稳定）
data$datetime <- ymd_hms(data$datetime, tz = "Europe/Berlin")

# 创建轨迹并转换为UTM投影
trk.test <- make_track(data, E, N, datetime, crs = 4326) %>% 
  st_transform(crs = 32632)

# 自动匹配UTM分区的方法（替代手动输入EPSG）
# trk.test <- make_track(data, E, N, datetime, crs = 4326) %>% 
#   st_transform(crs = utm_zone(.))

# 计算步长
stp.test <- steps(trk.test)

# 查看结果，sl_已转为米单位
view(stp.test)

验证结果 #

转换坐标系后，sl_列的数值会变成实际的米数，比如示例数据中的第一个步长会从约7.9e-4度变为约217米，符合真实移动距离。

内容的提问来源于stack exchange，提问作者fx0003