C#实现两点经纬度距离计算

基本原理与准备工作

地球并非完美的球体，而是一个近似椭球体。在实际开发中，为了简化计算并保证足够的精度（通常在几米误差范围内），我们通常假设地球为一个平均半径的球体。标准的大圆距离计算主要依赖Haversine公式，该公式能有效避免在小距离计算时的浮点数精度丢失问题。

• 地球平均半径 (R)：通常取 6371.004 km 或近似值 6378.137 km。
• 角度转换：经纬度输入通常为十进制度数，而三角函数库需要弧度值。转换关系为：弧度 = 度数 × π / 180。

C# 距离计算代码实现

以下是完整的C#实现代码，包含辅助方法和主计算逻辑：

private const double EARTH_RADIUS = 6378.137; // 地球半径，单位千米

/// 

/// 将角度转换为弧度
/// 
private static double rad(double d)
{
    return d * Math.PI / 180.0;
}

/// 

/// 根据两点经纬度计算距离
/// 
/// 起点纬度
/// 起点经度
/// 终点纬度
/// 终点经度
/// 两点间距离，单位为千米
public static double GetDistance(double lat1, double lng1, double lat2, double lng2)
{
    double radLat1 = rad(lat1);
    double radLat2 = rad(lat2);
    
    // 纬度差
    double a = radLat1 - radLat2;
    // 经度差
    double b = rad(lng1) - rad(lng2);

    // Haversine 公式核心计算
    double s = 2 * Math.Asin(Math.Sqrt(Math.Pow(Math.Sin(a/2),2) +
        Math.Cos(radLat1)*Math.Cos(radLat2)*Math.Pow(Math.Sin(b/2),2)));
    
    s = s * EARTH_RADIUS; // 乘以地球半径得到距离
    s = Math.Round(s * 10000) / 10000; // 保留四位小数以提高精度感观
    
    return s;
}

关键步骤解析

第1步：理解参数含义

方法接收四个 double 类型的参数，分别代表起点的纬度(lat1)、经度(lng1)和终点的纬度(lat2)、经度(lng2)。注意：在北半球和赤道附近，这些数值为正；南半球和高纬度地区的数值需根据具体坐标系规范处理，但Haversine公式本身对符号不敏感，因为它使用的是正弦和余弦的平方或绝对差值。

第2步：执行几何计算

代码首先通过 rad() 方法将所有角度转换为弧度。接着计算纬度和经度的差值。核心公式 2 * Asin(Sqrt(...)) 是Haversine距离的标准形式，它直接给出了两点在球面上的圆心角（弧度制）。

第3步：单位换算与结果输出

将圆心角乘以地球半径 EARTH_RADIUS 即可得到以千米为单位的弧长距离。最后通过 Math.Round 对结果进行微调，避免因浮点数运算产生的微小误差干扰业务逻辑判断。

注意事项与验证

• 单位一致性：返回的距离单位与定义的 EARTH_RADIUS 单位一致。若需英里(mile)，可将半径设为 3958.76 或在返回前乘以换算系数 0.621371。
• 精度范围：该公式适用于大多数常规GIS应用场景。对于极高精度要求的科学测量，应考虑使用更复杂的椭球体模型（如Vincenty公式）。
• 数据源：经纬度数据通常来自GPS模块（WGS84坐标系）或地图服务商API，确保来源坐标系统一是关键。

简要总结

上述代码提供了一个轻量级且高效的C#解决方案，用于快速估算地理空间距离。通过封装 GetDistance 静态方法，开发者可以轻松集成到地图标记、路线规划或本地生活服务等业务场景中。以上就是C#根据两点经纬度计算距离的详细 implementation 内容，更多关于C# GIS开发技巧的资料请关注本站其它相关文章！