CoordTransform.js 3.64 KB
Newer Older
xinzhedeai's avatar
xinzhedeai committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121
var CoordTransform = new function(){
	//定义一些常量
	const BD_FACTOR = (3.14159265358979324 * 3000.0) / 180.0;
	const PI = 3.1415926535897932384626;
	const RADIUS = 6378245.0;
	const EE = 0.00669342162296594323;

	/**
	 * 判断是否在国内,不在国内则不做偏移
	 * @param lng
	 * @param lat
	 * @returns {boolean}
	 */
	function out_of_china(lng, lat) {
		return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false);
	};

	return {
		/**
		   * BD-09(百度坐标系) To GCJ-02(火星坐标系)
		   * @param lng
		   * @param lat
		   * @returns {number[]}
		   */
		BD09ToGCJ02: function(){
			let x = +lng - 0.0065;
			let y = +lat - 0.006;
			let z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * BD_FACTOR);
			let theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * BD_FACTOR);
			let gg_lng = z * Math.cos(theta);
			let gg_lat = z * Math.sin(theta);
			return [gg_lng, gg_lat];
		},

		/**
		   * GCJ-02(火星坐标系) To BD-09(百度坐标系)
		   * @param lng
		   * @param lat
		   * @returns {number[]}
		   * @constructor
		   */
		GCJ02ToBD09: function(){
			let z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * BD_FACTOR);
			let theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * BD_FACTOR);
			let bd_lng = z * Math.cos(theta) + 0.0065;
			let bd_lat = z * Math.sin(theta) + 0.006;
			return [bd_lng, bd_lat];
		},
		/**
		   * GCJ-02(火星坐标系) To WGS-84(世界大地坐标系)
		   * @param lng
		   * @param lat
		   * @returns {number[]}
		   * @constructor
		   */
		GCJ02ToWGS84: function(lng, lat){
			if (out_of_china(lng, lat)) {
			  return [lng, lat];
			} else {
			  let d = this.delta(lng, lat);
			  let mgLng = lng + d[0];
			  let mgLat = lat + d[1];
			  return [lng * 2 - mgLng, lat * 2 - mgLat];
			}
		},

		/**
		   *
		   * @param lng
		   * @param lat
		   * @returns {number[]}
		   */
		delta: function(lng, lat){
			let dLng = this.transformLng(lng - 105, lat - 35);
			let dLat = this.transformLat(lng - 105, lat - 35);
			const radLat = (lat / 180) * PI;
			let magic = Math.sin(radLat);
			magic = 1 - EE * magic * magic;
			const sqrtMagic = Math.sqrt(magic);
			dLng = (dLng * 180) / ((RADIUS / sqrtMagic) * Math.cos(radLat) * PI);
			dLat = (dLat * 180) / (((RADIUS * (1 - EE)) / (magic * sqrtMagic)) * PI);
			return [dLng, dLat];
		},
		
		transformLng: function(lng, lat){
			lat = +lat;
			lng = +lng;
			let ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
			ret += ((20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0) / 3.0;
			ret += ((20.0 * Math.sin(lng * PI) + 40.0 * Math.sin((lng / 3.0) * PI)) * 2.0) / 3.0;
			ret += ((150.0 * Math.sin((lng / 12.0) * PI) + 300.0 * Math.sin((lng / 30.0) * PI)) * 2.0) / 3.0;
			return ret;
		},
		transformLat: function(lng, lat){
			let ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
			ret += ((20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0) / 3.0;
			ret += ((20.0 * Math.sin(lat * PI) + 40.0 * Math.sin((lat / 3.0) * PI)) * 2.0) / 3.0;
			ret += ((160.0 * Math.sin((lat / 12.0) * PI) + 320 * Math.sin((lat * PI) / 30.0)) * 2.0) / 3.0;
			return ret;
		},
		/**
		   * WGS-84(世界大地坐标系) To GCJ-02(火星坐标系)
		   * @param lng
		   * @param lat
		   * @returns {number[]}
		   */
		WGS84ToGCJ02: function(lng, lat){
			if (this.out_of_china(lng, lat)) {
			  return [lng, lat];
			} else {
			  let d = this.delta(lng, lat);
			  return [lng + d[0], lat + d[1]];
			}
		},
		out_of_china: function(lng, lat){
			return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
		},
	};
};

export var CoordTransform = CoordTransform;