计算两个经纬度之间的距离

计算两个经纬度之间的距离，通常可以使用Haversine公式或者类似的球面三角学方法。这些公式假设地球是一个完美的球体，虽然实际上地球是一个椭球体，但这种假设在大多数情况下已经足够精确。以下是计算两个经纬度之间距离的详细步骤和公式：

1. 假设和变量定义

• 假设地球的平均半径为R（通常取6371.004千米）。
• 设第一点A的经纬度为(LonA, LatA)，第二点B的经纬度为(LonB, LatB)。
• 经度(Longitude)和纬度(Latitude)需要转换成弧度进行计算，可以使用公式：弧度 = 度数 * π / 180。

2. 转换经纬度为弧度

# 假设有以下经纬度（这里仅作为示例，实际计算时应替换为具体值）
lonA, latA = 116.4074, 39.9042  # 北京的经纬度
lonB, latB = 121.4737, 31.2304  # 上海的经纬度

# 转换为弧度
rad_lonA, rad_latA = lonA * math.pi / 180, latA * math.pi / 180
rad_lonB, rad_latB = lonB, latB * math.pi / 180

3. 使用Haversine公式计算距离

Haversine公式是一种计算两点间大圆距离（即球面最短距离）的公式。假设我们有两个点 $A(la{t}_1,lo{n}_1)$ 和 $B(la{t}_2,lo{n}_2)$，其中 $lat$ 是纬度，$lon$ 是经度，$R$ 是地球半径（通常取6371公里）。

1. 计算两点间经纬度差

   • $\Delta lat=la{t}_2-la{t}_1$
   • $\Delta lon=lo{n}_2-lo{n}_1$

2. 将角度转换为弧度

   • $la{t}_1,lo{n}_1,la{t}_2,lo{n}_2$ 都需要转换为弧度（如果输入是度的话）

3. 使用Haversine公式

   • $a={\sin }^2\left(\frac{\Delta lat}{2}\right)+\cos (la{t}_1)\cdot \cos (la{t}_2)\cdot {\sin }^2\left(\frac{\Delta lon}{2}\right)$
   • $c=2\cdot \text{atan2}\left(\sqrt{a},\sqrt{1-a}\right)$
   • $d=R\cdot c$

其中，$d$ 是两点间的距离（单位：公里）。

在Excel中计算两个经纬度之间的距离，可以使用多种公式。这里提供一个常用的、基于哈弗辛公式（Haversine formula）的简化版本，该公式考虑了地球的曲率，因此在较长距离的计算中也能保持较高的准确性。请注意，由于Excel中的三角函数默认以弧度为单位，因此需要将经纬度从度转换为弧度。

4. Excel公式计算

1. 准备数据：

   • 在Excel的A列输入第一个地点的经度。
   • 在B列输入第一个地点的纬度。
   • 在C列输入第二个地点的经度。
   • 在D列输入第二个地点的纬度。

2. 输入公式：

   • 在E列（或任意空白列）的第一个单元格中输入以下公式来计算两个地点之间的距离（以千米为单位）：

     =6371*ACOS(COS(RADIANS(B2))*COS(RADIANS(D2))*COS(RADIANS(C2-A2))+SIN(RADIANS(B2))*SIN(RADIANS(D2)))
     

     其中，6371是地球的平均半径（以千米为单位），RADIANS()函数将度转换为弧度，ACOS()函数计算反余弦值，COS()和SIN()函数分别计算余弦和正弦值。

3. 应用公式：

   • 将上述公式拖拽或复制到E列的其他单元格中，以计算其他经纬度对之间的距离。

5. JavaScript代码实现

function haversineGreatCircleDistance(latitudeFrom, longitudeFrom, latitudeTo, longitudeTo, earthRadius = 6371) {
  // 将十进制度数转化为弧度
  const latFrom = toRadians(latitudeFrom);
  const lonFrom = toRadians(longitudeFrom);
  const latTo = toRadians(latitudeTo);
  const lonTo = toRadians(longitudeTo);

  // 应用haversine公式
  const deltaLat = latTo - latFrom;
  const deltaLon = lonTo - lonFrom;
  
  const a = Math.sin(deltaLat / 2) * Math.sin(deltaLat / 2) +
            Math.cos(latFrom) * Math.cos(latTo) *
            Math.sin(deltaLon / 2) * Math.sin(deltaLon / 2);
  const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
  
  return earthRadius * c; // 返回距离，单位千米
}

function toRadians(angle) {
  return angle * (Math.PI / 180);
}

// 使用示例
const distance = haversineGreatCircleDistance(34.0522, -118.2437, 36.1699, -115.1398);
console.log(distance); // 输出两点间距离，单位公里

这段代码提供了一个使用Haversine公式计算两点间球面距离的JavaScript函数，并给出了一个使用示例。你可以根据自己的需要调整经纬度输入和地球半径（如果需要的话）。

6. Python代码实现

import math

# 地球平均半径
R = 6371.004

# 转换经纬度为弧度（已在第2步中完成）

# 计算a
a = math.sin((rad_latA - rad_latB) / 2)**2 + math.cos(rad_latA) * math.cos(rad_latB) * math.sin((rad_lonA - rad_lonB) / 2)**2

# 计算c
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))

# 计算距离
distance = R * c

# 输出结果
print(f"两点之间的距离是：{distance:.2f}千米")

7. SQL代码实现

ALTER FUNCTION [dbo].[F_GEO_GetDistance]
(
	@lon1 FLOAT = 116.407526,
	@lat1 FLOAT = 39.904030,
	@lon2 FLOAT = 121.473701,
	@lat2 FLOAT = 31.230416
)
RETURNS VARCHAR(100)
AS
BEGIN
  
-- 地球平均半径（公里）  
DECLARE @earthRadius FLOAT = 6371.137;  
  
-- SQL查询，计算两个点之间的距离  
RETURN   
    (@earthRadius *   
    ACOS(  
        COS(RADIANS(@lat1)) * COS(RADIANS(@lat2)) *   
        COS(RADIANS(@lon2) - RADIANS(@lon1)) +   
        SIN(RADIANS(@lat1)) * SIN(RADIANS(@lat2))  
    ))
END
/******************************************
PRINT dbo.F_GEO_GetDistance(116.407526,39.904030,121.473701,31.230416)
*******************************************/

8. 注意事项

• 上述公式和代码是基于地球是一个完美球体的假设，这在大多数情况下是足够的。
• 在实际应用中，地球的形状对距离计算的影响可以忽略不计，除非在极高精度的应用中（如导航系统）。
• 如果需要更高的精度，可以考虑使用WGS-84椭球体模型等更复杂的地球模型。
• 如果需要以其他单位（如英里或米）表示距离，可以相应地调整地球的平均半径值。例如，以米为单位时，可以使用6371000（因为1千米=1000米）。
• 纬度相差一度，在地球表面的实际距离大约是111千米（在赤道附近），但这个值会随着纬度的增加而减少。经度相差一度的实际距离则取决于所在的纬度，计算公式为$111\times \cos (\text{纬度})$千米。