经纬度的转换，百度、高德、腾讯、GPS互转

/// <summary>
/// 经纬度转换
/// </summary>
public static class LocationUtil
{
    private const double Pi = 3.1415926535897932384626;
    private const double A = 6378245.0;
    private const double EE = 0.00669342162296594323;
    private const double X_PI = Pi * 3000.0 / 180.0;

    /// <summary>
    /// 是否在中国区域
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static bool IsInChina(double lat, double lng)
    {
        return lng >= 72.004 && lng <= 137.8347 && lat >= 0.8293 && lat <= 55.8271;
    }

    #region GPS <=> GCJ02

    /// <summary>
    /// 将WGS84坐标转换为GCJ02坐标（中国国测局坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Wgs84ToGcj02(double lat, double lng)
    {
        if (!IsInChina(lat, lng)) return (lat, lng);
        var (dLat, dLng) = Delta(lat, lng);
        return (lat + dLat, lng + dLng);
    }

    /// <summary>
    /// 将GCJ02坐标转换为WGS84坐标（GPS坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Gcj02ToWgs84(double lat, double lng)
    {
        if (!IsInChina(lat, lng)) return (lat, lng);
        var (dLat, dLng) = Delta(lat, lng);
        return (lat - dLat, lng - dLng);
    }

    private static (double, double) Delta(double lat, double lng)
    {
        double dLat = TransformLat(lng - 105.0, lat - 35.0);
        double dLng = TransformLng(lng - 105.0, lat - 35.0);
        double radLat = lat / 180.0 * Pi;
        double magic = Math.Sin(radLat);
        magic = 1 - EE * magic * magic;
        double sqrtMagic = Math.Sqrt(magic);
        dLat = (dLat * 180.0) / ((A * (1 - EE)) / (magic * sqrtMagic) * Pi);
        dLng = (dLng * 180.0) / (A / sqrtMagic * Math.Cos(radLat) * Pi);
        return (dLat, dLng);
    }

    private static double TransformLat(double x, double y)
    {
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y +
                     0.1 * x * y + 0.2 * Math.Sqrt(Math.Abs(x));
        ret += (20.0 * Math.Sin(6.0 * x * Pi) +
                20.0 * Math.Sin(2.0 * x * Pi)) * 2.0 / 3.0;
        ret += (20.0 * Math.Sin(y * Pi) +
                40.0 * Math.Sin(y / 3.0 * Pi)) * 2.0 / 3.0;
        ret += (160.0 * Math.Sin(y / 12.0 * Pi) +
                320.0 * Math.Sin(y * Pi / 30.0)) * 2.0 / 3.0;
        return ret;
    }

    private static double TransformLng(double x, double y)
    {
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x +
                     0.1 * x * y + 0.1 * Math.Sqrt(Math.Abs(x));
        ret += (20.0 * Math.Sin(6.0 * x * Pi) +
                20.0 * Math.Sin(2.0 * x * Pi)) * 2.0 / 3.0;
        ret += (20.0 * Math.Sin(x * Pi) +
                40.0 * Math.Sin(x / 3.0 * Pi)) * 2.0 / 3.0;
        ret += (150.0 * Math.Sin(x / 12.0 * Pi) +
                300.0 * Math.Sin(x / 30.0 * Pi)) * 2.0 / 3.0;
        return ret;
    }

    #endregion

    #region GCJ02 <=> BD09

    /// <summary>
    /// 将GCJ02坐标转换为BD09坐标（百度坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Gcj02ToBd09(double lat, double lng)
    {
        double x = lng, y = lat;
        double z = Math.Sqrt(x * x + y * y) + 0.00002 * Math.Sin(y * X_PI);
        double theta = Math.Atan2(y, x) + 0.000003 * Math.Cos(x * X_PI);
        double bdLng = z * Math.Cos(theta) + 0.0065;
        double bdLat = z * Math.Sin(theta) + 0.006;
        return (bdLat, bdLng);
    }

    /// <summary>
    /// 将BD09坐标转换为GCJ02坐标（百度坐标转国测局坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Bd09ToGcj02(double lat, double lng)
    {
        double x = lng - 0.0065, y = lat - 0.006;
        double z = Math.Sqrt(x * x + y * y) - 0.00002 * Math.Sin(y * X_PI);
        double theta = Math.Atan2(y, x) - 0.000003 * Math.Cos(x * X_PI);
        double ggLng = z * Math.Cos(theta);
        double ggLat = z * Math.Sin(theta);
        return (ggLat, ggLng);
    }

    #endregion

    #region WGS84 <=> BD09 (组合)

    /// <summary>
    /// 将WGS84坐标转换为BD09坐标（百度坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Wgs84ToBd09(double lat, double lng)
    {
        var (gLat, gLng) = Wgs84ToGcj02(lat, lng);
        return Gcj02ToBd09(gLat, gLng);
    }

    /// <summary>
    /// 将BD09坐标转换为WGS84坐标（GPS坐标）
    /// </summary>
    /// <param name="lat">纬度(X，100多那个)</param>
    /// <param name="lng">经度(Y，30多那个)</param>
    /// <returns></returns>
    public static (double lat, double lng) Bd09ToWgs84(double lat, double lng)
    {
        var (gLat, gLng) = Bd09ToGcj02(lat, lng);
        return Gcj02ToWgs84(gLat, gLng);
    }

    #endregion
}