-1

我正在使用 Haversine 公式来计算两个经纬度对的距离。

function getDistanceFromLatLonInKm(lat1,lon1,lat2,lon2) {
    var R = 6371; // Radius of the earth in km
    var dLat = deg2rad(lat2-lat1);
    var dLon = deg2rad(lon2-lon1); 
    var lat1 = deg2rad(lat1);
    var lat2 = deg2rad(lat2);

    var a = 
        Math.sin(dLat/2) * Math.sin(dLat/2) + 
        Math.sin(dLon/2) * Math.sin(dLon/2) * 
        Math.cos(lat1) * Math.cos(lat2);

    var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
    var d = R * c;
    return d;
}

给定一个起点（lat1，lat2），在直线上移动所需的距离和角度，我需要确定终点（如 lat2 和 lon2）。

请参阅下面的尝试：

function getFinalLatLon(lat1, lon1, distance, angle) {
    var R = 6371; // Radius of the earth in km
    var c = distance/R;
    // Math.atan2(Math.sqrt(a), Math.sqrt(1-a)) = c/2
    var a = // stuck here

    // looking for this part of the code

    return [lat2, lon2];
}
4

2 回答 2

4

如果您是水平移动，您可以将经度增加distance / (R * cos(lat)). 不需要atan。

编辑：由于您想要一个适用于一般情况的公式，请考虑以下几何推导：

  • 正视图：

    在此处输入图像描述

  • 侧面图：

    在此处输入图像描述

  • 整体设置：

    在此处输入图像描述

笔记：

  • r是你的起始位置的单位向量，s是端点。

    在此处输入图像描述

  • a, b, c是帮助计算的中间向量。

  • (θ, φ)是（纬度，经度）坐标。
  • γ是您将要行驶的方向的方位。
  • δ是经过的角度（距离/半径R = 6400000m）。

我们需要a, b垂直于北r并a与北对齐。这给出了：

在此处输入图像描述

c由（简单三角函数）给出：

在此处输入图像描述

因此我们得到s（通过一些非常乏味的代数）：

在此处输入图像描述

现在我们可以计算使用的最终（纬度，经度）坐标s：

在此处输入图像描述

代码：

function deg2rad(deg) { return deg * (Math.PI / 180.0) }
function rad2deg(rad) { return rad * (180.0 / Math.PI) }

function getFinalLatLong(lat1, long1, distance, angle, radius) {
    // calculate angles
    var delta = distance / radius,
        theta = deg2rad(lat1),
        phi   = deg2rad(long1),
        gamma = deg2rad(angle);

    // calculate sines and cosines
    var c_theta = Math.cos(theta), s_theta = Math.sin(theta);
    var c_phi   = Math.cos(phi)  , s_phi   = Math.sin(phi)  ;
    var c_delta = Math.cos(delta), s_delta = Math.sin(delta);
    var c_gamma = Math.cos(gamma), s_gamma = Math.sin(gamma);

    // calculate end vector
    var x = c_delta * c_theta * c_phi - s_delta * (s_theta * c_phi * c_gamma + s_phi * s_gamma);
    var y = c_delta * c_theta * s_phi - s_delta * (s_theta * s_phi * c_gamma - c_phi * s_gamma);
    var z = s_delta * c_theta * c_gamma + c_delta * s_theta;

    // calculate end lat long
    var theta2 = Math.asin(z), phi2 = Math.atan2(y, x);

    return [rad2deg(theta2), rad2deg(phi2)];
}

测试用例：

  1. 输入(lat, long) = (45, 0), angle = 0, distance = radius * deg2rad(90)=> (45, 180)（如我之前所说）

  2. 输入(lat, long) = (0, 0), angle = 90, distance = radius * deg2rad(90)=> (0, 90)（如预期 - 从赤道开始，向东行经 90 经度）

  3. 输入(lat, long) = (54, 29), angle = 36, distance = radius * deg2rad(360)=> (54, 29)（如预期的那样 - 从任意随机位置开始并沿任意方向绕一圈）

  4. 有趣的案例： input (lat, long) = (30, 0)，其他一切都一样。=> (0, 90)（我们预计(30, 90)？ - 不是从赤道开始，向北行进 90 度）

    原因是向北 90 度不是东方（如果你不在赤道）！该图应说明原因：

    在此处输入图像描述

    如您所见，向北 90 度的运动路径不是向东方向。

于 2017-07-21T10:20:05.093 回答
0

我刚刚在这里发现了一个类似的问题，然后我按照解决方案提出了一个适用于我的情况的函数。

希望它可以帮助别人：

function getFinalLatLon(lat1, lon1, distance, angle){ 
    function deg2rad(deg) {
        return deg * (Math.PI/180)
    }
    // dy = R*sin(theta) 
    var dy = distance * Math.sin(deg2rad(angle)) 
    var delta_latitude = dy/110574
    // One degree of latitude on the Earth's surface equals (110574 meters
    delta_latitude = parseFloat(delta_latitude.toFixed(6));

    // final latitude = start_latitude + delta_latitude
    var lat2 = lat1 + delta_latitude

    // dx = R*cos(theta) 
    var dx = distance * Math.cos(deg2rad(angle)) 
    // One degree of longitude equals 111321 meters (at the equator)
    var delta_longitude = dx/(111321*Math.cos(deg2rad(lat1))) 
    delta_longitude = parseFloat(delta_longitude.toFixed(6));

    // final longitude = start_longitude + delta_longitude
    var lon2 = lon1 + delta_longitude

    return [lat2, lon2];
}

水平移动的角度为 0 度。你可以随意切换。如果有人向北移动，那将是 90 度。西北135度等等...

于 2017-07-21T10:40:57.590 回答