MATLAB函数计算两个坐标之间的距离（纬度和经度）

Question

如何使用 MATLAB R2015a（以米为单位）计算两个世界地图坐标（纬度和经度）之间的距离？

Answer 1

如果您无法访问 MATLAB 映射工具箱，则一个简单的近似方法是使用半正弦公式。以下是链接的摘录：

半正矢公式是一个在导航中很重要的方程，根据经度和纬度给出球体上两点之间的大圆距离。它是球面三角学中更通用的公式（半正矢定律）的一个特例，涉及球面三角形的边和角。

这是一个 MATLAB 实现：

function rad = radians(degree) 
% degrees to radians
    rad = degree .* pi / 180;
end; 

function [a,c,dlat,dlon]=haversine(lat1,lon1,lat2,lon2)
% HAVERSINE_FORMULA.AWK - converted from AWK 
    dlat = radians(lat2-lat1);
    dlon = radians(lon2-lon1);
    lat1 = radians(lat1);
    lat2 = radians(lat2);
    a = (sin(dlat./2)).^2 + cos(lat1) .* cos(lat2) .* (sin(dlon./2)).^2;
    c = 2 .* asin(sqrt(a));
    arrayfun(@(x) printf("distance: %.4f km\n",6372.8 * x), c);
end;

[a,c,dlat,dlon] = haversine(36.12,-86.67,33.94,-118.40); % BNA to LAX