swift - 从 DMS 计算方位

Question

我很难让它发挥作用。我尝试了许多不同的公式，但总是得到相同的结果。 基于相同的数据(坐标)我可以计算距离。但不是轴承。 我已经没有想法和关键字来寻找可行的公式了。

感谢大家的帮助!

let fLat = d2r(degrees: lat)
let fLng = d2r(degrees: long)
let tLat = d2r(degrees: lat2)
let tLng = d2r(degrees: long2)

let a = CGFloat(sin(fLng - tLng) * cos(tLat))
let b = CGFloat(cos(fLat) * sin(tLat) - sin(fLat) * cos(tLat) * cos(fLng-tLng))

let bearing = atan2(a,b)

func d2r(degrees : Double) -> Double {
    return degrees * Double.pi / 180.0
}

Answer 1

鉴于评论中的示例，您似乎只是遇到了符号错误:

CYWG = N 49 54 6 W097 14 4, CYYQ = N58 44 4 W094 03 9, I should have an heading of 11 degrees True.

let lat = 49.0 + 54.0 / 60.0 + 6.0 / 3600.0
let long = -(97.0 + 14.0 / 60.0 + 4.0 / 3600.0)
let lat2 = 58.0 + 44.0 / 60.0 + 4.0 / 3600.0
let long2 = -(94.0 + 3.0 / 60.0 + 9.0 / 3600.0)

let fLat = d2r(degrees: lat)
let fLng = d2r(degrees: long)
let tLat = d2r(degrees: lat2)
let tLng = d2r(degrees: long2)

let a = sin(fLng - tLng) * cos(tLat)
let b = cos(fLat) * sin(tLat) - sin(fLat) * cos(tLat) * cos(fLng - tLng)

let bearing = atan2(a,b)
bearing * 180.0 / Double.pi // -10.58473709161037

所以你的代码的答案是〜-11。检查我自己的 Vincenty's formula 的实现我得到的准确方位角是

let (d, az, raz) = vincentyInverse(p1, to: p2, ellipsoid: .wgs84)
az.toDegrees // 10.61078843895496

所以看来您的公式的所有问题都在于您从 fLng 中减去了 tLng，而不是反之亦然。将两个地方的 fLng - tLng 替换为 tLng - fLng 会得到 10.58 作为答案。

当然，如果您想要精确的方位角，您应该使用 Vincenty 或 Robbins 或其他一些考虑地球非球形形状的公式 - 您使用的公式仅在完美球体上准确......