有谁知道,为什么我的步长sl_没有以米显示?
我正在使用amt软件包,并尝试执行基于gps移动数据的综合步骤选择分析。
我对此非常陌生,也没有在google中发现太多的搜索。我的猜测是,因为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))发布于 2022-11-15 15:33:50
这样啊,原来是这么回事。实际上,它只是在添加参数lonlat = TRUE。steps()-function的缺省值是lonlat = FALSE。
https://stackoverflow.com/questions/74388857
复制相似问题