Haversine 距离算法详解

Haversine 距离算法是计算地球表面两点球面直线距离的经典算法，日常用的地图测距、打车软件预估里程，背后都有它的身影。

一、算法的核心用途

我们生活的地球是一个近似球体，如果要计算两个地点（比如北京到上海）的'直线距离'，不能直接用平面几何的勾股定理（因为地球表面是曲面）。

Haversine 算法的作用，就是基于两点的经纬度坐标，计算它们在地球球面上的最短距离（这个最短距离也叫大圆距离，即穿过球心的平面切割球面形成的圆弧长度）。

二、必须掌握的前置知识

在理解公式前，先记住 3 个关键概念：

1. 经纬度的定义
   • 纬度 (latitude)：衡量地点南北位置，范围是 [-90°, 90°]，赤道是 0°，北极是 90°N，南极是 90°S。
   • 经度 (longitude)：衡量地点东西位置，范围是 [-180°, 180°]，本初子午线是 0°，向东为东经，向西为西经。
2. 角度与弧度的转换数学中三角函数（sin、cos）的计算需要弧度值，而我们日常用的经纬度是角度值，因此必须先转换：
   • 弧度 = 角度 × π / 180°
   • 角度 = 弧度 × 180° / π
3. 地球半径的取值地球不是完美球体，赤道半径约 6378km，极半径约 6357km。日常计算取平均半径 R = 6371km 即可满足精度需求。

三、Haversine 公式拆解

1. 公式的数学表达式

假设地球表面有两点：

• 点 A：纬度 lat1，经度 lon1
• 点 B：纬度 lat2，经度 lon2

Haversine 公式的最终形式为：

d = 2R · arcsin(√(sin²(Δφ/2) + cos(φ1)·cos(φ2)·sin²(Δλ/2)))

其中：

• Δφ = lat2 − lat1：两点的纬度差
• Δλ = lon2 − lon1：两点的经度差
• R：地球平均半径（6371km）
• d：两点的球面直线距离

2. 公式的通俗理解（分 4 步计算）

我们不用死记公式，而是把计算拆成 4 个简单步骤：

1. 角度转弧度将 lat1, lon1, lat2, lon2 全部转换成弧度值（记为 φ1, λ1, φ2, λ2）。
2. 计算差值计算纬度差 Δφ=φ2−φ1，经度差 Δλ=λ2−λ1。
3. 计算核心根式代入公式计算根号内的部分：a=sin²(Δφ/2)+cos(φ1)·cos(φ2)·sin²(Δλ/2)
4. 计算最终距离代入公式计算距离：d=2R·arcsin(√a)

四、零基础能看懂的 Python 代码实现

下面给出完整的 Python 代码，逐行添加中文注释，你可以直接复制运行。我们以北京（39.9042°N, 116.4074°E）到上海（31.2304°N, 121.4737°E）为例计算距离。

# 导入数学库，用于三角函数和弧度转换
import math

def haversine_distance(lat1, lon1, lat2, lon2):
    
    
    R = 
    
    
    phi1 = math.radians(lat1)  
    lambda1 = math.radians(lon1)  
    phi2 = math.radians(lat2)  
    lambda2 = math.radians(lon2)  
    
    
    delta_phi = phi2 - phi1
    delta_lambda = lambda2 - lambda1
    
    
    a = math.sin(delta_phi / ) **  + \
        math.cos(phi1) * math.cos(phi2) * \
        math.sin(delta_lambda / ) ** 
    
    
    c =  * math.atan2(math.sqrt(a), math.sqrt( - a))
    
    distance = R * c
     distance


beijing_lat, beijing_lon = , 
shanghai_lat, shanghai_lon = , 


distance = haversine_distance(beijing_lat, beijing_lon, shanghai_lat, shanghai_lon)
()