有谁知道,为什么我的步长sl_不显示在米?
我正在使用amt包,并尝试执行基于gps运动数据的集成步骤选择分析。
我对这个很陌生,在谷歌上搜索没有发现太多。我的猜测是,由于CRS,它混乱了,我必须使用UTM坐标吗?
不幸的是,我不能发布我的数据集,如果你需要更多的信息,我想我可以创建一些随机的运动数据。
感谢您抽出时间来帮忙。
导入数据后,我将其转换为曲目。
trk <- make_track(data, E, N, datetime, crs = 4326)
> class(trk)
[1] "track_xyt" "track_xy" "tbl_df" "tbl" "data.frame"
然后执行steps()
函数。
stp <- steps(trk)
它工作得很好,只是我的sl_是这样显示的:
7.914206e-04
5.160599e-04
...
library(amt)
library(tidyverse)
library(lubridate)
setwd("path")
data <- read.csv2("test.csv", skip = 1, header = TRUE, sep = ",")
glimpse(data)
# format chr to dttm
data$datetime <- strptime(data$datetime, "%Y-%m-%d %H:%M:%S")
# set timezone
tz(data$datetime) <- "Europe/Berlin"
# make track out of data frame
trk.test <- make_track(data, E, N, datetime, crs = 4326)
class(trk.test)
# default plotting methods
plot(trk.test)
# turning tracks into steps
stp.test <- steps(trk.test)
view(stp.test)
WGS 84, EPSG 4326
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))
我明白了。它实际上只是加上了参数lonlat = TRUE
。steps()
-function默认为lonlat = FALSE
。