MATLAB函数计算两个坐标(经纬度)之间的距离

Question

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

Answer 1

如果您没有访问MATLAB映射工具箱，那么一个简单的近似是使用Haversine公式。以下是链接的节选:

哈弗斯公式是一个在导航中很重要的方程，它给出了球面上两点之间从经度和纬度的大圆距离。它是球面三角学中一个更一般的公式——哈弗斯定律的一个特例，它与球面三角形的边和角有关。

这里是一个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 kmn",6372.8 * x), c);
end;
[a,c,dlat,dlon] = haversine(36.12,-86.67,33.94,-118.40); % BNA to LAX

Answer 2

如果您有访问映射工具箱，那么可能在此页上描述的功能可以帮助您。

在这种情况下，您可以使用函数distance(LAT1,LON1,LAT2,LON2)来获得连接两个点的大圆弧的长度(以度为单位)。然后您可以使用deg2km, deg2nm或deg2sm将其转换为公里或英里。

以km为单位计算距离的一条直线是:

deg2km(distance(lat1, lon1, lat2, lon2))