关于度分秒到数字度和几大坐标转换的工具类
public class MapLocationHelper { #region 构造函数 /// <summary> /// 构造函数 /// </summary> public MapLocationHelper() { } #endregion #region 方法 /// <summary> /// 数字经纬度和度分秒经纬度转换 (Digital degree of latitude and longitude and vehicle to latitude and longitude conversion) /// </summary> /// <param name="digitalLati_Longi">数字经纬度</param> /// <return>度分秒经纬度</return> static public string ConvertDigitalToDegrees(string digitalLati_Longi) { double digitalDegree = Convert.ToDouble(digitalLati_Longi); return ConvertDigitalToDegrees(digitalDegree); } /// <summary> /// 数字经纬度和度分秒经纬度转换 (Digital degree of latitude and longitude and vehicle to latitude and longitude conversion) /// </summary> /// <param name="digitalDegree">数字经纬度</param> /// <return>度分秒经纬度</return> static public string ConvertDigitalToDegrees(double digitalDegree) { const double num = 60; int degree = (int)digitalDegree; double tmp = (digitalDegree - degree) * num; int minute = (int)tmp; double second = (tmp - minute) * num; string degrees = "" + degree + "°" + minute + "′" + second + "″"; return degrees; } /// <summary> /// 度分秒经纬度(必须含有'°')和数字经纬度转换 /// </summary> /// <param name="digitalDegree">度分秒经纬度</param> /// <return>数字经纬度</return> static public double ConvertDegreesToDigital(string degrees) { const double num = 60; double digitalDegree = 0.0; int d = degrees.IndexOf('°'); //度的符号对应的 Unicode 代码为：00B0[1]（六十进制），显示为°。 if (d < 0) { return digitalDegree; } string degree = degrees.Substring(0, d); digitalDegree += Convert.ToDouble(degree); int m = degrees.IndexOf('′'); //分的符号对应的 Unicode 代码为：2032[1]（六十进制），显示为′。 if (m < 0) { return digitalDegree; } string minute = degrees.Substring(d + 1, m - d - 1); digitalDegree += ((Convert.ToDouble(minute)) / num); int s = degrees.IndexOf('″'); //秒的符号对应的 Unicode 代码为：2033[1]（六十进制），显示为″。 if (s < 0) { return digitalDegree; } string second = degrees.Substring(m + 1, s - m - 1); digitalDegree += (Convert.ToDouble(second) / (num * num)); return digitalDegree; } /// <summary> /// 度分秒经纬度(必须含有'/')和数字经纬度转换 /// </summary> /// <param name="digitalDegree">度分秒经纬度</param> /// <param name="cflag">分隔符</param> /// <return>数字经纬度</return> static public double ConvertDegreesToDigital_default(string degrees) { char ch = '/'; return ConvertDegreesToDigital(degrees, ch); } /// <summary> /// 度分秒经纬度和数字经纬度转换 /// </summary> /// <param name="digitalDegree">度分秒经纬度</param> /// <param name="cflag">分隔符</param> /// <return>数字经纬度</return> static public double ConvertDegreesToDigital(string degrees, char cflag) { const double num = 60; double digitalDegree = 0.0; int d = degrees.IndexOf(cflag); if (d < 0) { return digitalDegree; } string degree = degrees.Substring(0, d); digitalDegree += Convert.ToDouble(degree); int m = degrees.IndexOf(cflag, d + 1); if (m < 0) { return digitalDegree; } string minute = degrees.Substring(d + 1, m - d - 1); digitalDegree += ((Convert.ToDouble(minute)) / num); int s = degrees.Length; if (s < 0) { return digitalDegree; } string second = degrees.Substring(m + 1, s - m - 1); digitalDegree += (Convert.ToDouble(second) / (num * num)); return digitalDegree; } #endregion }
上面是度分秒转度
下面是几个坐标互转
1 /* ---------------------------------------------------------- 2 * 文件名称：MapConverter.cs 3 * 4 * 开发环境： 5 * Visual Studio V2017 6 * 7 * 版本历史： 8 * V1.0 2017年05月16日 9 * 坐标转换器 10 * 11 * 说明： 12 * WGS84：为一种大地坐标系，也是目前广泛使用的GPS全球卫星定位系统使用的坐标系。 13 * GCJ02：又称火星坐标系，是由中国国家测绘局制定的地理坐标系统，是由WGS84加密后得到的坐标系。 14 * BD09：为百度坐标系，在GCJ02坐标系基础上再次加密。其中bd09ll表示百度经纬度坐标，bd09mc表示百度墨卡托米制坐标。 15 * 16 * 参考资料： 17 ------------------------------------------------------------ */ 18 /// <summary> 19 /// WGS-84、GCJ-02（火星坐标系）、BD-09（百度坐标系）之间的坐标转换器 20 /// </summary> 21 public static class MapConverter 22 { 23 /// <summary> 24 /// 圆周率 25 /// </summary> 26 private const double PI = 3.1415926535897932384626; 27 private const double X_PI = PI * 3000.0 / 180.0; 28 29 /// <summary> 30 /// 地理位置是否位于中国以外 31 /// </summary> 32 /// <param name="wgLat">WGS-84坐标纬度</param> 33 /// <param name="wgLon">WGS-84坐标经度</param> 34 /// <returns> 35 /// true：国外 36 /// false：国内 37 /// </returns> 38 public static bool OutOfChina(double wgLat, double wgLon) 39 { 40 if (wgLon < 72.004 || wgLon > 137.8347) return true; 41 if (wgLat < 0.8293 || wgLat > 55.8271) return true; 42 43 return false; 44 } 45 46 /// <summary> 47 /// WGS-84坐标系转火星坐标系 (GCJ-02) 48 /// </summary> 49 /// <param name="wgLat">WGS-84坐标纬度</param> 50 /// <param name="wgLon">WGS-84坐标经度</param> 51 /// <param name="mgLat">输出：GCJ-02坐标纬度</param> 52 /// <param name="mgLon">输出：GCJ-02坐标经度</param> 53 public static void WGS84ToGCJ02(double wgLat, double wgLon, out double mgLat, out double mgLon) 54 { 55 if (OutOfChina(wgLat, wgLon)) 56 { 57 mgLat = wgLat; 58 mgLon = wgLon; 59 } 60 else 61 { 62 double dLat; 63 double dLon; 64 Delta(wgLat, wgLon, out dLat, out dLon); 65 mgLat = wgLat + dLat; 66 mgLon = wgLon + dLon; 67 } 68 } 69 70 /// <summary> 71 /// 火星坐标系 (GCJ-02)转WGS-84坐标系 72 /// </summary> 73 /// <param name="mgLat">GCJ-02坐标纬度</param> 74 /// <param name="mgLon">GCJ-02坐标经度</param> 75 /// <param name="wgLat">输出：WGS-84坐标纬度</param> 76 /// <param name="wgLon">输出：WGS-84坐标经度</param> 77 public static void GCJ02ToWGS84(double mgLat, double mgLon, out double wgLat, out double wgLon) 78 { 79 if (OutOfChina(mgLat, mgLon)) 80 { 81 wgLat = mgLat; 82 wgLon = mgLon; 83 } 84 else 85 { 86 double dLat; 87 double dLon; 88 Delta(mgLat, mgLon, out dLat, out dLon); 89 wgLat = mgLat - dLat; 90 wgLon = mgLon - dLon; 91 } 92 } 93 94 /// <summary> 95 /// 火星坐标系 (GCJ-02)转WGS-84坐标系 96 /// </summary> 97 /// <param name="mgLat">GCJ-02坐标纬度</param> 98 /// <param name="mgLon">GCJ-02坐标经度</param> 99 /// <param name="wgLat">输出：WGS-84坐标纬度</param> 100 /// <param name="wgLon">输出：WGS-84坐标经度</param> 101 public static void GCJ02ToWGS84Exact(double mgLat, double mgLon, out double wgLat, out double wgLon) 102 { 103 const double InitDelta = 0.01; 104 const double Threshold = 0.000001; 105 106 double dLat = InitDelta; 107 double dLon = InitDelta; 108 double mLat = mgLat - dLat; 109 double mLon = mgLon - dLon; 110 double pLat = mgLat + dLat; 111 double pLon = mgLon + dLon; 112 113 double nLat; 114 double nLon; 115 116 int i = 0; 117 do 118 { 119 wgLat = (mLat + pLat) / 2; 120 wgLon = (mLon + pLon) / 2; 121 122 WGS84ToGCJ02(wgLat, wgLon, out nLat, out nLon); 123 124 dLat = nLat - mgLat; 125 dLon = nLon - mgLon; 126 if ((Math.Abs(dLat) < Threshold) && (Math.Abs(dLon) < Threshold)) break; 127 128 if (dLat > 0) pLat = wgLat; else mLat = wgLat; 129 if (dLon > 0) pLon = wgLon; else mLon = wgLon; 130 } while (++i <= 30); 131 } 132 133 /// <summary> 134 /// 百度坐标系 (BD-09)转火星坐标系 (GCJ-02) 135 /// </summary> 136 /// <param name="bdLat">百度坐标系纬度</param> 137 /// <param name="bdLon">百度坐标系经度</param> 138 /// <param name="mgLat">输出：GCJ-02坐标纬度</param> 139 /// <param name="mgLon">输出：GCJ-02坐标经度</param> 140 public static void BD09ToGCJ02(double bdLat, double bdLon, out double mgLat, out double mgLon) 141 { 142 double x = bdLon - 0.0065; 143 double y = bdLat - 0.006; 144 double z = Math.Sqrt(x * x + y * y) - 0.00002 * Math.Sin(y * X_PI); 145 double theta = Math.Atan2(y, x) - 0.000003 * Math.Cos(x * X_PI); 146 mgLat = z * Math.Sin(theta); 147 mgLon = z * Math.Cos(theta); 148 } 149 150 /// <summary> 151 /// 火星坐标系 (GCJ-02)转百度坐标系 (BD-09) 152 /// </summary> 153 /// <param name="mgLat">GCJ-02坐标纬度</param> 154 /// <param name="mgLon">GCJ-02坐标经度</param> 155 /// <param name="bdLat">输出：百度坐标系纬度</param> 156 /// <param name="bdLon">输出：百度坐标系经度</param> 157 public static void GCJ02ToBD09(double mgLat, double mgLon, out double bdLat, out double bdLon) 158 { 159 double x = mgLon; 160 double y = mgLat; 161 double z = Math.Sqrt(x * x + y * y) + 0.00002 * Math.Sin(y * X_PI); 162 double theta = Math.Atan2(y, x) + 0.000003 * Math.Cos(x * X_PI); 163 bdLat = z * Math.Sin(theta) + 0.006; 164 bdLon = z * Math.Cos(theta) + 0.0065; 165 } 166 167 /// <summary> 168 /// WGS-84坐标系转百度坐标系 (BD-09) 169 /// </summary> 170 /// <param name="wgLat">WGS-84坐标纬度</param> 171 /// <param name="wgLon">WGS-84坐标经度</param> 172 /// <param name="bdLat">输出：百度坐标系纬度</param> 173 /// <param name="bdLon">输出：百度坐标系经度</param> 174 public static void WGS84ToBD09(double wgLat, double wgLon, out double bdLat, out double bdLon) 175 { 176 double mgLat; 177 double mgLon; 178 179 WGS84ToGCJ02(wgLat, wgLon, out mgLat, out mgLon); 180 GCJ02ToBD09(mgLat, mgLon, out bdLat, out bdLon); 181 } 182 183 /// <summary> 184 /// 百度坐标系 (BD-09)转WGS-84坐标系 185 /// </summary> 186 /// <param name="bdLat">百度坐标系纬度</param> 187 /// <param name="bdLon">百度坐标系经度</param> 188 /// <param name="wgLat">输出：WGS-84坐标纬度</param> 189 /// <param name="wgLon">输出：WGS-84坐标经度</param> 190 public static void BD09ToWGS84(double bdLat, double bdLon, out double wgLat, out double wgLon) 191 { 192 double mgLat; 193 double mgLon; 194 195 BD09ToGCJ02(bdLat, bdLon, out mgLat, out mgLon); 196 GCJ02ToWGS84(mgLat, mgLon, out wgLat, out wgLon); 197 } 198 199 public static double Distance(double LatA, double LonA, double LatB, double LonB) 200 { 201 const double EarthR = 6371000.0; 202 203 double x = Math.Cos(LatA * PI / 180.0) * Math.Cos(LatB * PI / 180.0) * Math.Cos((LonA - LonB) * PI / 180); 204 double y = Math.Sin(LatA * PI / 180.0) * Math.Sin(LatB * PI / 180.0); 205 double s = x + y; 206 if (s > 1) s = 1; 207 if (s < -1) s = -1; 208 209 return Math.Acos(s) * EarthR; 210 } 211 212 private static void Delta(double Lat, double Lon, out double dLat, out double dLon) 213 { 214 const double AXIS = 6378245.0; 215 const double EE = 0.00669342162296594323; 216 217 dLat = TransformLat(Lon - 105.0, Lat - 35.0); 218 dLon = TransformLon(Lon - 105.0, Lat - 35.0); 219 double radLat = Lat / 180.0 * PI; 220 double magic = Math.Sin(radLat); 221 magic = 1 - EE * magic * magic; 222 double sqrtMagic = Math.Sqrt(magic); 223 dLat = (dLat * 180.0) / ((AXIS * (1 - EE)) / (magic * sqrtMagic) * PI); 224 dLon = (dLon * 180.0) / (AXIS / sqrtMagic * Math.Cos(radLat) * PI); 225 } 226 227 private static double TransformLat(double x, double y) 228 { 229 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)); 230 ret += (20.0 * Math.Sin(6.0 * x * PI) + 20.0 * Math.Sin(2.0 * x * PI)) * 2.0 / 3.0; 231 ret += (20.0 * Math.Sin(y * PI) + 40.0 * Math.Sin(y / 3.0 * PI)) * 2.0 / 3.0; 232 ret += (160.0 * Math.Sin(y / 12.0 * PI) + 320 * Math.Sin(y * PI / 30.0)) * 2.0 / 3.0; 233 return ret; 234 } 235 236 private static double TransformLon(double x, double y) 237 { 238 double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.Sqrt(Math.Abs(x)); 239 ret += (20.0 * Math.Sin(6.0 * x * PI) + 20.0 * Math.Sin(2.0 * x * PI)) * 2.0 / 3.0; 240 ret += (20.0 * Math.Sin(x * PI) + 40.0 * Math.Sin(x / 3.0 * PI)) * 2.0 / 3.0; 241 ret += (150.0 * Math.Sin(x / 12.0 * PI) + 300.0 * Math.Sin(x / 30.0 * PI)) * 2.0 / 3.0; 242 return ret; 243 } 244 }
经纬度计算完了以后就是距离计算代码
public class LocationUtil { ///// <summary> /////计算两点GPS坐标的距离（单位：米） ///// </summary> ///// <param name="n1">第一点的纬度坐标</param> ///// <param name="e1">第一点的经度坐标</param> ///// <param name="n2">第二点的纬度坐标</param> ///// <param name="e2">第二点的经度坐标</param> //public static double Distance(double n1, double e1, double n2, double e2) //{ // double jl_jd = 102834.74258026089786013677476285;//每经度单位米; // double jl_wd = 111712.69150641055729984301412873;//每纬度单位米; // double b = Math.Abs((e1 - e2) * jl_jd); // double a = Math.Abs((n1 - n2) * jl_wd); // return Math.Sqrt((a * a + b * b)); //} ///// <summary> ///// 获取维度差 ///// </summary> ///// <param name="km">千米</param> //public static double GetLatitudeDifference(double km) //{ // return km * (1d / 111d); //} ///// <summary> ///// 获取经度差 ///// </summary> ///// <param name="km">千米</param> //public static double GetLongitudeDifference(double km, double latitude) //{ // return km * (1d / (111d * Math.Cos(latitude))); //} private const double EARTH_RADIUS = 6378.137; //地球半径 private static double Rad(double d) { return d * Math.PI / 180.0; } public static double Distance(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); 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; } }