根据经纬度计算地球上两点之间的距离——Haversine公式介绍及计算步骤

最新推荐文章于 2024-10-11 11:59:05 发布
最新推荐文章于 2024-10-11 11:59:05 发布
阅读量2.3w 收藏 189
点赞数 42
文章标签： 几何学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/gaocuisheng/article/details/126060795
版权

目录

摘要

1.半正矢公式（Haversine Formula）介绍

2.半正矢公式应用

3.半正矢公式计算

3.1 主要思路

3.2 计算步骤

3.2.1 平面向量计算方法

3.2.2 空间向量计算方法

摘要

写本文的出发点是需要在Qlik中根据经纬度计算地球上两点间的距离。我在社区上搜到了相关公式的分享，这个公式叫做“半正矢公式”。对于“半正矢”知之甚少的我，决定从头到尾将公式计算一遍，并通过这篇文章记录所有的步骤，希望它能对需要的人有所帮助。这也算是温习高中所学的三角函数了。

本文的部分图片源自其他网站，剩下的图都是通过GeoGebra制作的。这个工具非常方便，推荐大家使用。

(GeoGebra - the world’s favorite, free math tools used by over 100 million students and teachers)

1.半正矢公式（Haversine Formula）介绍

公式： Haversin(\theta )=sin^2{\left (\theta/2 \right ) }

(i)

如图 (i), OAC是一个圆，O是其圆心，A和C是圆上的点，r 是⊙OAC的半径。

设：r = 1, 那么 OA = OC = r =1;

过点A做一条垂线，交线段OC于点B，∠OBA=90°。

设：∠AOB=θ;

因此，AB = sinθ, OB = cosθ, BC = 反正弦 = 1 – cosθ.

所谓的半正矢就是“一半的反正弦”（Haversine = “half reversed sin”），所以：

 

最后得出： Haversin(\theta )=sin^2{\left (\theta/2 \right ) }

2.半正矢公式应用

回到应用场景，已知地球上的两个点A和B，现在我们想知道这两点之间的最短距离。地球可以近似看作为一个球体，那么AB间的最短距离就是过AB两点的大圆圆弧，这个大圆是地球的内切圆，半径和圆心点都于地球相同。

(ii)

设：A点坐标 = (Ø1, λ1); B点坐标 = (Ø2, λ2);

Ø = 纬度, λ = 经度

d = A 和 B 的最短距离

R = 地球半径 = 6371 km

把 sin^{2}\left ( \frac{\theta }{2} \right ) ​​​替换为 Haversin(θ)  可得：

 

过去的水手们可以查余弦函数、半正矢函数和半正矢函数的逆函数的数值表，从而求出距离d。只要有数值表，还是挺方便的。[1]

[1]https://plus.maths.org/content/lost-lovely-haversine

3.半正矢公式计算

3.1 主要思路

(iii)

如图 (iii), ⊙OAB 是地球的内切圆，圆心为O。

距离计算步骤：

  1. 计算AB，即点A和B之间的直线距离；
  2. 通过AB和地球半径R，计算∠AOB 的角度；
  3. 通过∠AOB 和R计算 \widehat{AB} ，即⊙ OAB上AB的弧长，这个弧长就是最短距离。

3.2 计算步骤

3.2.1 平面向量计算方法

 (iv)

与第2节相同，设：

A点坐标 = (Ø1, λ1)； B点坐标 = (Ø2, λ2)；

Ø = 纬度，λ = 经度；

d = AB间的最短距离；

R = 地球半径 = 6371 km；

O = 地球球心；

N = 北极点；

点C和E与点A在同一条经度线上；

点D和F与点B在同一条经度线上;

点D 与点A在同一条纬度线上；

点C 与点B在同一条纬度线上；

点E和F又是在赤道上，点A所在经度线与点B所在经度线相较于北极点。

可得：

C点坐标 = (Ø2, λ1)； D点坐标 = (Ø1, λ2)；

E点坐标 = (0, λ1)； F点坐标 = (0, λ2)；

引入 sin(\frac{\theta }{2}) ：

(v)

如图(v)，OXY是一个圆，设r 为⊙OXY的半径，且r = 1；

OX = OY = r = 1, OZ是线段XY的垂线；

设∠XOY=θ，那么∠XOZ=∠YOZ=θ/2；

所以 XZ=sin(\frac{\theta }{2})  ， 且 XY=2XZ=2sin(\frac{\theta }{2}) 。

回到图 (iv), ∠AOC = Ø2 – Ø1;

因为 OA = OC =R，所以 AC =  ;

同理，BD = AC；EF = ；

点O’ 是点A和点D所在圆的圆心，O’A 长度等于圆的半径。

过点A做一条垂线，相交OE于G，∠AOG = ∠AOE = Ø1；

所以 :  ；

由于 ∠GOO’ = ∠AGO = 90°，OG = O’A，所以 O’A = ；

∠AO’D = λ2 – λ1, 可得 AD = ；

同理， BC =  ；

知道AD和BC的长度后，我们在等腰梯形ADBC中就可以计算AB的长度了。

 (vi)

在等腰梯形ADBC中，过点A做一条垂线相交BC于点E。

那么：

；

；  

；

 ；

现在，AB线段的长度已经计算出来，接下来就是最后一步计算 \widehat{AB}的弧长。

(vii)

 

如图(vii)，⊙OAB 是地球的内接圆，圆心为O，OA = R.

设：AC = √a, ∠AOB = b, ∠OCA = 90°；

那么 AC = AB/2;

 ；

 ；

 ； 

   ；

 

​​​​​​​

   ；

 

把 sin^{2}\left ( \frac{\theta }{2} \right ) ​​​替换为 Haversin(θ)  ：  

与第二节的公式一致，计算完毕。

或者我们可以用正切函数tan:

  ​​​​​​​

​​​​​​​

​​​​​​​

​​​​​​​

这与Qlik 社区分享的计算公式一致[2]。

[2] https://community.qlik.com/t5/QlikView-App-Dev/Haversine-formula-to-find-distance-between-two-lat-long-points/m-p/234630

3.2.2 空间向量计算方法

如通过空间向量计算距离，我们需要把经纬度的坐标转换为三维空间坐标。

如图 (viii)，我们建立了一个坐标系，Z轴就是地轴所在的轴，平面XOY是赤道所在的平面，平面XOZ是本初子午线所在的平面。A点是地球上任意一点。

设：A点坐标 = (Ø1, λ1)，R是地球半径；

A 点的空间坐标 = (X1, Y1, Z1)

 (viii)

过A点做一条垂线相交平面XOY于A’点，则：

 ；

X1 = OC, Y1 = CA’, Z1 = A’A

因此， ；

同理，我们有一点 B (Ø2, λ2)，则  ；

那么：

​​​​​​​

 设 ，那么我们可通过余弦定理计算θ 的值：

  ​​​​​​​

​​​​​​​

 把cos(θ) 替换为 Haversin(θ)，结果与上文等价。

参考文章

Haversine formula to find distance between two lat... - Qlik Community - 234630

Lost but lovely: The haversine | plus.maths.org

地理空间距离计算及优化（根据两个点经纬度计算距离）_ArthurKingYs的博客-CSDN博客_两个经纬度算距离公式及方法

半正矢公式（Haversine公式）

如何计算地球上两点的距离（附公式推导）_呆呆papa的博客-CSDN博客_地球两点距离

Gcs4real
关注
计算球面上两点间最短距离多语言版本
Alvin's Blog
11-23 2498
Haversine球面半正矢公式，用来计算球面上连点之间的最短距离， 这里有几种语言版本的实现方式，就摘了几种，还有其他例如Smalltalk, Tcl, NCL的就不贴了。//C# using System; namespace HaversineFormula { /// <summary> /// 距离类型，英里/公里 /// </summary> publ
Haversine公式：用于计算两个球面位置之间的球面距离
weixin_38256051的博客
08-28 5185
然而，对于较长距离，特别是跨越大片陆地或海洋的距离，Haversine 公式可能会引入一定的误差。需要注意的是，Haversine 公式是一种近似计算方法，适用于许多实际应用中的球面距离计算。然而，对于更高精度的计算，特别是在大范围内或需要考虑地球形状的其他因素时，可能需要使用更复杂的算法或地理空间扩展。Haversine 是一种常用的球面三角函数，用于计算两个球面位置之间的球面距离。它是根据地球的球体模型进行计算的，可以用于近似计算两个经纬度坐标之间的球面距离。
Haversine 算法计算地球上两个点之间的最短距离
YOUSEEYA的博客
09-12 453
Haversine 算法用于计算地球上两个点之间的最短距离（假设地球是一个理想的球体），这种距离被称为大圆距离。它是航海和航空领域常用的一种方法，用来估算两个地理坐标点（经度和纬度）之间的距离。该算法的计算原理基于球面三角学。
经纬度计算距离
最新发布
zheng_ruiguo的专栏
10-11 1253
经纬度可以使用球面三角法来计算距离。一种常用的方法是使用Haversine公式。计算地球上两点之间的距离时，一般采用大圆距离近似公式，即球面上两点之间的距离等于两点之间的弧长。如果要考虑地球的椭圆形状和半径的不均匀误差，可以使用Vincenty公式来计算。这个公式考虑了地球的椭圆形状和半径的不均匀性，可以得到更精确的结果。
Haversine公式：精确计算地球上任意两点之间的距离
qq_44103359的博客
05-11 2010
在地理信息系统（GIS）和相关应用中，根据地球上的两点的经纬度计算它们之间的距离是一个常见的需求。Haversine公式是一种用于计算地球上两点之间距离的精确方法，它考虑了地球的曲率。Haversine公式提供了一种准确计算地球上两点之间距离的方法。通过将经纬度转换为弧度，计算经纬度差，并应用Haversine公式，我们可以得到两点之间的最短距离。公式考虑了地球的半径和两点之间的经纬度差，来计算两点之间的最短距离（大圆距离）。将转换后的经纬度和经纬度差代入Haversine公式，计算两点之间的距离。
haversine-formula
CY_TEC的博客
05-27 4298
haversine formula 给定一个点的经纬度和一个范围长度，计算以这个点为中心的正方形边界，初步实现代码： https://github.com/SunnyCat2013/haversine-distance/blob/master/haversine.py 球面上任意两点的中心角求解公式： hav(dr)=hav(ϕ2−ϕ1)+cos(ϕ1)cos(ϕ2)hav(λ2−λ1)...
Haversine公式
anzheng6118的博客
08-30 2551
1 Number.prototype.toRadians = function() { 2 return this * Math.PI / 180; 3 } 4 5 function distance(latitude1, longitude1, latitude2, longitude2) { 6 //R是地球的半径，以KM为单位 7 ...
如何计算地球上两点的距离（附公式推导）
热门推荐
呆呆papa的博客
01-16 1万+
前段时间，看了一些电子围栏的算法，对其中一段计算球面上两点距离的代码有些不解，然后找了一下相关算法，在维基百科的大圆距离词条中记录了相关的计算公式，大致思路就是求出这两点间的弧长对应的圆心角的余弦或正弦，然后利用反三角函数计算出圆心角的弧度，最后求出：弧长=弧度值 × 地球半径。 注：上图使用的是百度地图测距功能，测量湖北省襄阳市火车站出站口和位于吉林省长春市的地铁1号线所途径的长春站北地铁站的距离 一、具体实现 假设球面上有两点A(λ1 , φ1)、B(λ2, φ2)，λ 和 φ 分别表示它们在地图中的
经纬度计算距离,经纬度计算距离公式,matlab
09-10
本篇文章将深入探讨如何利用经纬度来计算地球上两点之间的距离，并介绍一个基于MATLAB的GUI工具，方便用户进行此类计算。 首先，让我们理解经纬度的基本概念。经度是指从本初子午线向东或向西测量的角度，0°经线...
使用Haversine公式计算距离：使用Haversine公式计算两个位置之间的距离-matlab开发
05-30
Haversine 公式给出了地球表面上两点之间的最短距离，忽略海拔、山丘等...... 该程序接受两个输入，它们可以是指定为度、分和秒或十进制度的位置。 输出是以公里、海里和英里为单位的距离。 例子： 半正弦('53 08...
C#通过经纬度计算2个点之间距离的实现代码
09-04
这个公式基于球面三角学，用于计算地球上两点之间的大圆距离。大圆距离是最短的路径，也就是飞机飞行的路线。哈弗辛公式假设地球是一个完美的球体，实际地球更接近椭球体，但这种简化在大多数情况下足够精确。 以下...
haversine:Haversine是一个实现Haversine公式的Go库
04-27
Haversine 软件包Haversine是一个实现Haversine公式的Go库。 Haversine公式根据其经度和纬度给出球体上两个点之间的大圆距离。 在这种情况下，球体是地球的表面。 黄点线是一个大圆弧。 它给出了两个黄色点之间的最短距离。 图片由USGS提供。 安装 去获取github.com/umahmood/haversine cd $ GOPATH / src / github.com / umahmood / haversine / 去测试-v。/ ... 用法 以下示例显示了如何计算地球表面上两个坐标之间的最短路径。 package main import ( "fmt" "github.com/umahmood/haversine" ) func main() { oxford := haversine.Coord{Lat: 51.
Haversine：Haversine公式的Crystal实现，可根据给定的纬度和经度来计算两个点之间的距离
02-05
Haversine：Haversine公式的Crystal实现，可根据给定的纬度和经度来计算两个点之间的距离
haversine:计算地球上连续点之间的半正弦（大圆）距离（以米为单位）。-matlab开发
05-30
计算地球表面上连续点之间的半正弦（大圆）距离（以米为单位）。 这些点被指定为纬度和经度的向量。
java实现经纬度计算距离
03-23
哈弗辛公式基于球面三角学，用于计算地球上任意两点之间的最短距离。 哈弗辛公式的计算过程如下： 1. 将经纬度转换为弧度：地球的平均半径约为6371公里，1度经线对应的弧长约为111.320公里，1度纬线弧长大约在111....
半正矢公式（Haversine公式）——根据经度纬度计算两点间距离
Williamcsj的博客
08-17 8663
半正矢公式（Haversine公式）——根据经度纬度计算两点间距离一、说明二、参考 一、说明 该计算得到的距离是两点之间的直线距离。 二、参考 坐标测算直线距离与坐标转换 半正矢公式：根据经度纬度计算两点间距离 半正矢公式（Haversine公式） ...
利用经纬度求解两点球面距离(Haversine formula)
Hardict的算法日记
04-02 5909
利用经纬度求解两点球面距离 Haversine formula h(θ)=sin2(θ2)=1−cos(θ)2h(\theta)=sin^2(\frac{\theta}{2})=\frac{1-cos(\theta)}{2}h(θ)=sin2(2θ​)=21−cos(θ)​ 则h(θ)=h(dR)=h(Δβ)+cos(β1)cos(β2)h(Δα)h(\theta)=h(\frac{d}{R})...
haversine公式，计算地图上两点之间的距离
唯有美食与爱不可辜负
12-14 2509
测试环境 python 3.6 django 3.1 代码部分 haversine 公式 用于计算地图上两点之间的直线距离 A点经纬度(lon1, lat1) B点经纬度(lon2, lat2) def haversine(lon1, lat1, lon2, lat2): """ 计算地图两点之间的的距离 在地球上（以十进制度数指定） lon1: 经度1 lat1: 纬度1 lon2: 经度2 lat2: 纬度2 """ # 将十进制度
Swift代码中的Haversine公式
weixin_45843292的博客
11-06 226
Here is an example of the implementation of the Haversine formula in Swift code. The Haversine formula is used to measure two gps point codes: swiftswiftimportUIKitfunchaversineDinstance(la1:Double, l...
Haversine公式通过经纬度计算两点距离
07-11
Haversine公式是一种常用的方法，通过经纬度计算地球上两点之间的距离。该公式基于球面三角形的概念，假设地球是一个完美的球体。 公式如下： d = 2r * arcsin(sqrt(sin^2((lat2-lat1)/2) + cos(lat1) * cos(lat2) * sin^2((lon2-lon1)/2))) 其中，d表示两点之间的距离，r表示地球的半径，lat1和lat2分别表示两点的纬度，lon1和lon2分别表示两点的经度。这里的纬度和经度需要使用弧度表示。 需要注意的是，Haversine公式对于小距离的计算结果比较准确，但在计算大距离时可能存在一定误差。此外，由于地球并非完全球形，其他更精确的公式也可以用于计算两点之间的距离。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值