我有GPS数据,我正在尝试计算行程最远部分的方位。我的数据如下所示:
Time Lat Long 18:11:25 -38.621449 142.931061 18:13:25 -38.62212 142.929871 18:15:33 -38.623333 142.92868 ... 23:16:26 -38.793068 142.827881 23:18:54 -38.795414 142.838516 23:23:02 -38.795803 142.839066 ... 09:59:49 -38.621532 142.930786
我正在尝试计算起点 (-38.621449 142.931061) 和最远点 (-38.793068,142.827881) 之间的方位角。
我已经能够通过从电子表格中手动提取 2 个位置并使用我的脚本计算方位角来做到这一点:
` bearing <- function(lat1,lon1,lat2,lon2){
lat1 <- angledim(lat1,"degrees","radians")
lon1 <- angledim(lon1,"degrees","radians")
lat2 <- angledim(lat2,"degrees","radians")
lon2 <- angledim(lon2,"degrees","radians")
dlon <- lon1 - lon2
theta <- atan2(cos(lat2)*sin(dlon),cos(lat1)*sin(lat2)-sin(lat1)*cos(lat2)*cos(dlon))
theta <- angledim(theta,"radians","degrees")
theta <- (theta+360) %% 360
theta
}`
有没有更简单的方法来计算这一点,而无需通过无休止的电子表格?
(如果格式不正确,请提前道歉,这是我的第一个问题!
谢谢!
我在这里使用的是化石包中的 earth.bear 函数,但您可以使用自己的自定义编写函数。
earth.bear(mydata$Long[1],mydata$Lat[1],
tail(mydata$Long, n=1),tail(mydata$Lat, n=1))