使用geosphere::bearing
我可以计算两条线的轴承,但是是否可以计算两个轴承之间的角度?
当然,您可以尝试减去或总结轴承,但在特定情况下,一个为负,另一个为正,这是行不通的。
例如,如果ber1
= - 175 和ber2
= 175,则两者之间的角度应为 10。
有什么建议吗?
我不确定现成的软件包,但如果您对解决方案感兴趣,那么您可以尝试
angle_diff <- function(theta1, theta2){
theta <- abs(theta1 - theta2) %% 360
return(ifelse(theta > 180, 360 - theta, theta))
}
它给出了示例轴承之间的角度 -175 和 175 作为
angle_diff(-175, 175)
#[1] 10