Matlab计算两点（经纬度坐标）距离（大弧法和Haversine法）

A-Chin

已于 2022-03-18 15:36:08 修改

阅读量4.8k

点赞数 2

分类专栏： matlab 文章标签： matlab 经纬度坐标距离大弧法 Haversine法、

于 2021-12-31 00:14:37 首次发布

本文链接：https://blog.csdn.net/l_j_kin/article/details/122247893

版权

matlab 专栏收录该内容

35 篇文章 18 订阅

订阅专栏

利用大弧法和Haversine法计算两个经纬度之间的距离。Python版本请参考 Python实现Haversine公式计算两点（经纬度坐标）距离。

R = 6371;

% great cricle arc
arcLAT1 = LAT * pi / 180.0;
arcLAT2 = LAT_minus * pi / 180.0;
arcLON1 = LON * pi / 180.0;
arcLON2 = LON_minus * pi / 180.0;

cos_sita = sin(arcLAT1) .* sin(arcLAT2) + ...
    cos(arcLAT1) .* cos(arcLAT2) .* cos(arcLON1 - arcLON2);
sita = acos(cos_sita);
tyTransDist_gca = R * sita;

% haversine
arcLAT1 = LAT * pi / 180.0;
arcLAT2 = LAT_minus * pi / 180.0;
arcLON1 = LON * pi / 180.0;
arcLON2 = LON_minus * pi / 180.0;
dLon = arcLAT1 - arcLAT2;
dLat = arcLON1 - arcLON2;

a = sin(dLat / 2) .^ 2 + ...
    + (1 - sin(dLat / 2) .^ 2 - sin((arcLAT1 + arcLAT2) / 2)) .* sin(dLon / 2) .^ 2;
c = 2 * asin(sqrt(a));
tyTransDist_hav = R .* c;