GCJ-02转WGS-84

  • WGS-84: 地球坐标系, 国际通用坐标系
  • GCJ-02: 火星坐标系, WGS84坐标系加密后的坐标系; Google国内地图、高德、QQ地图 使用
  • BD-09: 百度坐标系, GCJ02坐标系加密后的坐标系

javascript 方法

  1. (function (root, factory) {
  2.     if (typeof define === 'function' && define.amd) {
  3.                 // amd,注册为匿名模块
  4.         define([], factory);
  5.     } else if (typeof module === 'object' && module.exports) {
  6.                 //不能与严格的CommonJS一起工作，而只能在类似CommonJS的环境中工作，这些环境支持mode.exports，比如Node。
  7.         module.exports = factory();
  8.     } else {
  9.                 // 浏览器全局参数(根目录为窗口)
  10.         root.coordtransform = factory();
  11.     }
  12. }(this, function () {
  13.             // 定义一些常量
  14.     var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
  15.     var PI = 3.1415926535897932384626;
  16.     var a = 6378245.0;
  17.     var ee = 0.00669342162296594323;
  18.     /**
  19.     * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02) 的转换
  20.     * 即 百度 转 谷歌、高德
  21.     * @param bd_lng
  22.     * @param bd_lat
  23.     * @returns {*[]}
  24.     */
  25.     var bd09togcj02 = function bd09togcj02(bd_lng, bd_lat) {
  26.         var bd_lng = +bd_lng;
  27.         var bd_lat = +bd_lat;
  28.         var x = bd_lng - 0.0065;
  29.         var y = bd_lat - 0.006;
  30.         var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_PI);
  31.         var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_PI);
  32.         var gg_lng = z * Math.cos(theta);
  33.         var gg_lat = z * Math.sin(theta);
  34.         return [gg_lng, gg_lat]
  35.     };
  37.     /**
  38.     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
  39.     * 即 谷歌、高德 转 百度
  40.     * @param lng
  41.     * @param lat
  42.     * @returns {*[]}
  43.     */
  44.     var gcj02tobd09 = function gcj02tobd09(lng, lat) {
  45.         var lat = +lat;
  46.         var lng = +lng;
  47.         var z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_PI);
  48.         var theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_PI);
  49.         var bd_lng = z * Math.cos(theta) + 0.0065;
  50.         var bd_lat = z * Math.sin(theta) + 0.006;
  51.         return [bd_lng, bd_lat]
  52.     };
  54.     /**
  55.     * WGS-84 转 GCJ-02
  56.     * @param lng
  57.     * @param lat
  58.     * @returns {*[]}
  59.     */
  60.     var wgs84togcj02 = function wgs84togcj02(lng, lat) {
  61.         var lat = +lat;
  62.         var lng = +lng;
  63.         if (out_of_china(lng, lat)) {
  64.             return [lng, lat]
  65.         } else {
  66.             var dlat = transformlat(lng - 105.0, lat - 35.0);
  67.             var dlng = transformlng(lng - 105.0, lat - 35.0);
  68.             var radlat = lat / 180.0 * PI;
  69.             var magic = Math.sin(radlat);
  70.             magic = 1 - ee * magic * magic;
  71.             var sqrtmagic = Math.sqrt(magic);
  72.             dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
  73.             dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
  74.             var mglat = lat + dlat;
  75.             var mglng = lng + dlng;
  76.             return [mglng, mglat]
  77.         }
  78.     };
  80.     /**
  81.     * GCJ-02 转换为 WGS-84
  82.     * @param lng
  83.     * @param lat
  84.     * @returns {*[]}
  85.     */
  86.     var gcj02towgs84 = function gcj02towgs84(lng, lat) {
  87.         var lat = +lat;
  88.         var lng = +lng;
  89.         if (out_of_china(lng, lat)) {
  90.             return [lng, lat]
  91.         } else {
  92.             var dlat = transformlat(lng - 105.0, lat - 35.0);
  93.             var dlng = transformlng(lng - 105.0, lat - 35.0);
  94.             var radlat = lat / 180.0 * PI;
  95.             var magic = Math.sin(radlat);
  96.             magic = 1 - ee * magic * magic;
  97.             var sqrtmagic = Math.sqrt(magic);
  98.             dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
  99.             dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
  100.             var mglat = lat + dlat;
  101.             var mglng = lng + dlng;
  102.             return [lng * 2 - mglng, lat * 2 - mglat]
  103.         }
  104.     };
  106.     var transformlat = function transformlat(lng, lat) {
  107.         var lat = +lat;
  108.         var lng = +lng;
  109.         var ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
  110.         ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
  111.         ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;
  112.         ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
  113.         return ret
  114.     };
  116.     var transformlng = function transformlng(lng, lat) {
  117.         var lat = +lat;
  118.         var lng = +lng;
  119.         var ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
  120.         ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
  121.         ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;
  122.         ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;
  123.         return ret
  124.     };
  126.     /**
  127.     * 判断是否在国内, 不在国内则不做偏移
  128.     * @param lng
  129.     * @param lat
  130.     * @returns {boolean}
  131.     */
  132.     var out_of_china = function out_of_china(lng, lat) {
  133.         var lat = +lat;
  134.         var lng = +lng;
  135.                 // 纬度 3.86~53.55, 经度 73.66~135.05
  136.         return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
  137.     };
  139.     return {
  140.         bd09togcj02: bd09togcj02,
  141.         gcj02tobd09: gcj02tobd09,
  142.         wgs84togcj02: wgs84togcj02,
  143.         gcj02towgs84: gcj02towgs84
  144.     }
  145. }));
  147. // 使用方法
  148. //国测局坐标(火星坐标,比如高德地图在用),百度坐标,wgs84坐标(谷歌国外以及绝大部分国外在线地图使用的坐标)
  149. //百度经纬度坐标转国测局坐标
  150. var bd09togcj02 = coordtransform.bd09togcj02(116.404, 39.915);
  151. //国测局坐标转百度经纬度坐标
  152. var gcj02tobd09 = coordtransform.gcj02tobd09(116.404, 39.915);
  153. //wgs84转国测局坐标
  154. var wgs84togcj02 = coordtransform.wgs84togcj02(116.404, 39.915);
  155. //国测局坐标转wgs84坐标
  156. var gcj02towgs84 = coordtransform.gcj02towgs84(116.404, 39.915);
  157. console.log(bd09togcj02);
  158. console.log(gcj02tobd09);
  159. console.log(wgs84togcj02);
  160. console.log(gcj02towgs84);
  161. //result
  162. //bd09togcj02:   [ 116.39762729119315, 39.90865673957631 ]
  163. //gcj02tobd09:   [ 116.41036949371029, 39.92133699351021 ]
  164. //wgs84togcj02:  [ 116.41024449916938, 39.91640428150164 ]
  165. //gcj02towgs84:  [ 116.39775550083061, 39.91359571849836 ]