如何计算不跨越国际日期变更线的恒向线距离与方位角？

解决方案 #

核心逻辑为判断两点经度差的绝对值，若超过180°则调整其中一点的经度偏移360°，将两点经度差压缩到±180°范围内，再直接调用geosphere原生函数计算即可，无需复杂三角函数推导。

实现代码

library(geosphere)

# 计算不跨国际日期变更线的恒向线方位角与距离
get_rhumb_no_dateline <- function(p1, p2) {
  lon1 <- p1[["lon"]]
  lon2 <- p2[["lon"]]
  lat1 <- p1[["lat"]]
  lat2 <- p2[["lat"]]
  
  # 计算原始经度差
  lon_diff <- lon2 - lon1
  # 经度差绝对值超过180°时调整目标点经度
  if (abs(lon_diff) > 180) {
    lon2 <- ifelse(lon_diff > 0, lon2 - 360, lon2 + 360)
  }
  p2_adj <- c(lon = lon2, lat = lat2)
  
  # 计算调整后的恒向线参数
  bearing <- (bearingRhumb(p1, p2_adj) + 360) %% 360
  dist <- distRhumb(p1, p2_adj)
  
  # 若需生成绘制用的路径采样点，可打开以下注释，自动修正经度到正常范围
  # line_points <- destPointRhumb(p1, bearing, seq(0, dist, length.out = 100))
  # line_points[, "lon"] <- ifelse(line_points[, "lon"] > 180, line_points[, "lon"] - 360, line_points[, "lon"])
  
  return(list(bearing = bearing, dist = dist))
}

# 测试用例：中国到加拿大
p1 <- c("lon" = 104, "lat" = 36.6)
p2 <- c("lon" = -102, "lat" = 57.7)
res <- get_rhumb_no_dateline(p1, p2)

# 输出结果
print(paste0("不跨日界线恒向线方位角：", round(res$bearing, 2), "°"))
print(paste0("不跨日界线恒向线距离：", round(res$dist/1000, 2), "km"))

说明

• 调整经度的逻辑仅用于计算恒向线参数，最终绘制路径时，只需将采样点中大于180°的经度减去360°，即可得到适配常规地图的坐标，不会出现跨日界线的异常长线。
• 该方案完全复用geosphere包的原生计算逻辑，准确率和兼容性远高于手动三角函数推导，适配所有点对点的恒向线计算场景。

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