Python已知经纬度求两点距离

1. 引言

在使用Python进行地理信息处理时，经纬度是一个常见的数据形式。如果需要计算两个经纬度之间的距离，我们可以使用Haversine公式来实现。本文将介绍如何使用Python计算已知经纬度之间的距离，并提供相关的代码示例。

2. Haversine公式的介绍

Haversine公式是一种计算两个地球上任意两点（给定经度和纬度）之间距离的方法。该公式基于球面三角形的角度，在较短的距离范围内提供较高的准确性。

2.1 公式推导

根据Haversine公式，两点间的距离可以通过以下公式来计算：

import math
def haversine(lon1, lat1, lon2, lat2):
    """
    Calculate the distance between two points on the earth given their longitudes and latitudes in degrees.
    """
    # Convert degrees to radians
    lon1 = math.radians(lon1)
    lat1 = math.radians(lat1)
    lon2 = math.radians(lon2)
    lat2 = math.radians(lat2)
    
    # Haversine formula
    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = math.sin(dlat/2)**2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon/2)**2
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
    distance = 6371 * c  # Radius of earth in kilometers
    
    return distance

以上haversine()函数接受四个参数，分别是第一个点的经度lon1、纬度lat1，以及第二个点的经度lon2、纬度lat2。函数将返回两点之间的距离，单位为公里。

2.2 示例

现假设有两个经纬度坐标点，如下所示：

# Coordinates of Point A
lon1 = 12.34
lat1 = 56.78
# Coordinates of Point B
lon2 = 98.76
lat2 = 23.45

我们可以调用haversine()函数来计算这两个点的距离：

distance = haversine(lon1, lat1, lon2, lat2)
print("The distance between Point A and Point B is:", distance, "kilometers")

以上代码将输出：

The distance between Point A and Point B is: 7459.740459987601 kilometers

3. 小结

本文介绍了如何使用Python根据已知经纬度计算两个点之间的距离。使用Haversine公式可以在较短的距离范围内提供较高的准确性。通过编写haversine()函数，可以方便地计算任意两个地点之间的距离。