python 计算GCJ-02坐标两点直线距离

在Python中计算GCJ-02坐标系下两点之间的直线距离，通常需要先将坐标从高斯投影坐标转换为直角坐标，然后利用勾股定理计算两点间的距离。由于GCJ-02坐标并非常见的经纬度系统，而是中国自主研发的一种加偏算法，所以在处理时需要注意其坐标变换规则。

以下是基本步骤：

1. **坐标转换**：首先你需要了解GCJ-02坐标到笛卡尔坐标的转换公式。这种转换涉及到一些复杂的数学运算，包括对数、指数和乘法。通常在网上能找到相关的转换函数，比如百度地图API提供的`gcj02_to_wgs84`函数。

2. **计算距离**：在直角坐标系（如WGS84或EPSG:4326）下，你可以直接使用Haversine公式或者Vincenty's Inverse Formula计算两点间的大圆距离。这两个公式考虑了球面地球的曲率。

以下是一个简单的示例，假设你已经有了两个GCJ-02坐标点`(lat1, lon1)`和`(lat2, lon2)`：

from math import radians, sin, cos, sqrt, atan2

def haversine(lat1, lon1, lat2, lon2):
    R = 6371  # 地球平均半径，单位公里
    dLat = radians(lat2 - lat1)
    dLon = radians(lon2 - lon1)
    a = sin(dLat / 2) * sin(dLat / 2) + cos(radians(lat1)) * cos(radians(lat2)) * sin(dLon / 2) * sin(dLon / 2)
    c = 2 * atan2(sqrt(a), sqrt(1 - a))
    distance = R * c
    return distance

# 使用之前转换得到的直角坐标
distance_km = haversine(lat1, lon1, lat2, lon2)