1
/// <summary>
2
/// 经纬度转换
3
/// </summary>
4
public static class LocationUtil
5
{
6
    private const double Pi = 3.1415926535897932384626;
7
    private const double A = 6378245.0;
8
    private const double EE = 0.00669342162296594323;
9
    private const double X_PI = Pi * 3000.0 / 180.0;
10

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

22
    #region GPS <=> GCJ02
23

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

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

50
    private static (double, double) Delta(double lat, double lng)
51
    {
52
        double dLat = TransformLat(lng - 105.0, lat - 35.0);
53
        double dLng = TransformLng(lng - 105.0, lat - 35.0);
54
        double radLat = lat / 180.0 * Pi;
55
        double magic = Math.Sin(radLat);
56
        magic = 1 - EE * magic * magic;
57
        double sqrtMagic = Math.Sqrt(magic);
58
        dLat = (dLat * 180.0) / ((A * (1 - EE)) / (magic * sqrtMagic) * Pi);
59
        dLng = (dLng * 180.0) / (A / sqrtMagic * Math.Cos(radLat) * Pi);
60
        return (dLat, dLng);
61
    }
62

63
    private static double TransformLat(double x, double y)
64
    {
65
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y +
66
                     0.1 * x * y + 0.2 * Math.Sqrt(Math.Abs(x));
67
        ret += (20.0 * Math.Sin(6.0 * x * Pi) +
68
                20.0 * Math.Sin(2.0 * x * Pi)) * 2.0 / 3.0;
69
        ret += (20.0 * Math.Sin(y * Pi) +
70
                40.0 * Math.Sin(y / 3.0 * Pi)) * 2.0 / 3.0;
71
        ret += (160.0 * Math.Sin(y / 12.0 * Pi) +
72
                320.0 * Math.Sin(y * Pi / 30.0)) * 2.0 / 3.0;
73
        return ret;
74
    }
75

76
    private static double TransformLng(double x, double y)
77
    {
78
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x +
79
                     0.1 * x * y + 0.1 * Math.Sqrt(Math.Abs(x));
80
        ret += (20.0 * Math.Sin(6.0 * x * Pi) +
81
                20.0 * Math.Sin(2.0 * x * Pi)) * 2.0 / 3.0;
82
        ret += (20.0 * Math.Sin(x * Pi) +
83
                40.0 * Math.Sin(x / 3.0 * Pi)) * 2.0 / 3.0;
84
        ret += (150.0 * Math.Sin(x / 12.0 * Pi) +
85
                300.0 * Math.Sin(x / 30.0 * Pi)) * 2.0 / 3.0;
86
        return ret;
87
    }
88

89
    #endregion
90

91
    #region GCJ02 <=> BD09
92

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

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

125
    #endregion
126

127
    #region WGS84 <=> BD09 (组合)
128

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

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

153
    #endregion
154
}