【发布时间】:2015-11-05 16:16:07
【问题描述】:
如何使用 MATLAB R2015a(以米为单位)计算两个世界地图坐标(纬度和经度)之间的距离?
【问题讨论】:
标签: matlab coordinates distance
如何使用 MATLAB R2015a(以米为单位)计算两个世界地图坐标(纬度和经度)之间的距离?
【问题讨论】:
标签: matlab coordinates distance
如果您有权访问映射工具箱,那么this 页面上描述的功能可能会对您有所帮助。
在这种情况下,您可以使用函数distance(LAT1,LON1,LAT2,LON2) 来获取连接两个点的大圆弧的长度(以度为单位)。然后您可以使用deg2km、deg2nm 或deg2sm 将其转换为公里或英里。
一个以公里为单位的班轮将是:
deg2km(distance(lat1, lon1, lat2, lon2))
【讨论】:
如果您无法访问 MATLAB 映射工具箱,那么一个简单的近似方法是使用 Haversine 公式。以下是链接的摘录:
haversine 公式是一个在导航中很重要的方程,它给出了球体上两点与它们的经度和纬度之间的大圆距离。它是球面三角学中一个更一般的公式的特例,即半正弦定律,将球面三角形的边和角联系起来。
Here 是一个 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
【讨论】: