| | |
| | | return path; |
| | | } |
| | | |
| | | //定义一些常量 |
| | | const PI = 3.1415926535897932384626; |
| | | const a = 6378245.0; //长半轴 |
| | | const ee = 0.00669342162296594323; //扁率/*** GCJ02 转换为 WGS84* @param lng* @param lat* @returns {*[]}*/ |
| | | |
| | | /** |
| | | * 判断是否在国内,不在国内则不做偏移 |
| | | * @param lng |
| | | * @param lat |
| | | * @returns {boolean} |
| | | */ |
| | | function out_of_china(lng, lat) { |
| | | // 纬度3.86~53.55,经度73.66~135.05 |
| | | return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55); |
| | | } |
| | | |
| | | function transformlat(lng, lat) { |
| | | // lat = +lat lng = +lng |
| | | 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; |
| | | } |
| | | |
| | | function transformlng(lng, lat) { |
| | | // lat = +latlng = +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; |
| | | } |
| | | |
| | | /** |
| | | * 将gps经纬度转换为高德地图经纬度 |
| | | * @param {*} lnglats |
| | | * @param {*} callback |
| | | */ |
| | | function _convertLatlng(index, coor, lnglats, callback) { |
| | | function _convertLatlng(index, coor, lnglats, type = 'gps', callback) { |
| | | if (index < coor.length) { |
| | | var path = parse2LngLat(coor[index]); |
| | | // eslint-disable-next-line no-undef |
| | | AMap.convertFrom(path, 'gps', function (status, result) { |
| | | AMap.convertFrom(path, type, function (status, result) { |
| | | if (result.info === 'ok') { |
| | | lnglats.push.apply(lnglats, result.locations); |
| | | _convertLatlng(index + 1, coor, lnglats, callback); |
| | | _convertLatlng(index + 1, coor, lnglats, type, callback); |
| | | } |
| | | }); |
| | | } else { |
| | |
| | | var lng = ((dx / ed + (pos[0] * Math.PI) / 180.0) * 180.0) / Math.PI; |
| | | var lat = ((dy / ec + (pos[1] * Math.PI) / 180.0) * 180.0) / Math.PI; |
| | | return [lng, lat]; |
| | | }, |
| | | |
| | | /** |
| | | * 获取两个经纬度之间的角度(0度-360度) |
| | | */ |
| | | getAngle(lng_a, lat_a, lng_b, lat_b) { |
| | | var a = ((90 - lat_b) * Math.PI) / 180; |
| | | var b = ((90 - lat_a) * Math.PI) / 180; |
| | | var AOC_BOC = ((lng_b - lng_a) * Math.PI) / 180; |
| | | var cosc = |
| | | Math.cos(a) * Math.cos(b) + Math.sin(a) * Math.sin(b) * Math.cos(AOC_BOC); |
| | | var sinc = Math.sqrt(1 - cosc * cosc); |
| | | var sinA = (Math.sin(a) * Math.sin(AOC_BOC)) / sinc; |
| | | var A = (Math.asin(sinA) * 180) / Math.PI; |
| | | var res = 0; |
| | | if (lng_b > lng_a && lat_b > lat_a) res = A; |
| | | else if (lng_b > lng_a && lat_b < lat_a) res = 180 - A; |
| | | else if (lng_b < lng_a && lat_b < lat_a) res = 180 - A; |
| | | else if (lng_b < lng_a && lat_b > lat_a) res = 360 + A; |
| | | else if (lng_b > lng_a && lat_b == lat_a) res = 90; |
| | | else if (lng_b < lng_a && lat_b == lat_a) res = 270; |
| | | else if (lng_b == lng_a && lat_b > lat_a) res = 0; |
| | | else if (lng_b == lng_a && lat_b < lat_a) res = 180; |
| | | return res; |
| | | }, |
| | | |
| | | /** |
| | | * 获取两经纬度间的距离 |
| | | */ |
| | | getDistance(lng1, lat1, lng2, lat2) { |
| | | lat1 = lat1 || 0; |
| | | lng1 = lng1 || 0; |
| | | lat2 = lat2 || 0; |
| | | lng2 = lng2 || 0; |
| | | |
| | | var rad1 = (lat1 * Math.PI) / 180.0; |
| | | var rad2 = (lat2 * Math.PI) / 180.0; |
| | | var a = rad1 - rad2; |
| | | var b = (lng1 * Math.PI) / 180.0 - (lng2 * Math.PI) / 180.0; |
| | | var r = 6378137; |
| | | var distance = |
| | | r * |
| | | 2 * |
| | | Math.asin( |
| | | Math.sqrt( |
| | | Math.pow(Math.sin(a / 2), 2) + |
| | | Math.cos(rad1) * Math.cos(rad2) * Math.pow(Math.sin(b / 2), 2) |
| | | ) |
| | | ); |
| | | |
| | | return distance; |
| | | }, |
| | | |
| | | /** |
| | |
| | | return coors_GD; |
| | | }, |
| | | |
| | | convertFromGPS: function (gps, callback) { |
| | | convertFromGPS: function (gps, callback, type = 'gps') { |
| | | var coor = _prepare4convert(gps); |
| | | _convertLatlng(0, coor, [], function (result) { |
| | | _convertLatlng(0, coor, [], type, function (result) { |
| | | var gd = []; |
| | | result.forEach((r) => { |
| | | gd.push([r.lng, r.lat]); |
| | | }); |
| | | callback(gd); |
| | | }); |
| | | }, |
| | | |
| | | /** |
| | | * 高德地图坐标转GPS坐标算法 |
| | | */ |
| | | gcj02towgs84(lng, lat) { |
| | | // lat = +latlng = +lng |
| | | if (out_of_china(lng, lat)) { |
| | | return [lng, lat]; |
| | | } else { |
| | | let dlat = transformlat(lng - 105.0, lat - 35.0); |
| | | let dlng = transformlng(lng - 105.0, lat - 35.0); |
| | | let radlat = (lat / 180.0) * PI; |
| | | let magic = Math.sin(radlat); |
| | | magic = 1 - ee * magic * magic; |
| | | let sqrtmagic = Math.sqrt(magic); |
| | | dlat = (dlat * 180.0) / (((a * (1 - ee)) / (magic * sqrtmagic)) * PI); |
| | | dlng = (dlng * 180.0) / ((a / sqrtmagic) * Math.cos(radlat) * PI); |
| | | let mglat = Math.round((lat * 2 - lat - dlat) * 1000000) / 1000000; |
| | | let mglng = Math.round((lng * 2 - lng - dlng) * 1000000) / 1000000; |
| | | return [mglng, mglat]; |
| | | } |
| | | }, |
| | | |
| | | //从GPS转高德 |
| | | wgs84_To_Gcj02(lon, lat) { |
| | | if (out_of_china(lon, lat)) { |
| | | return [lon, lat]; |
| | | } else { |
| | | let dLat = transformlat(lon - 105.0, lat - 35.0); |
| | | let dLon = transformlng(lon - 105.0, lat - 35.0); |
| | | let radLat = (lat / 180.0) * PI; |
| | | let magic = Math.sin(radLat); |
| | | magic = 1 - ee * magic * magic; |
| | | let sqrtMagic = Math.sqrt(magic); |
| | | dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * PI); |
| | | dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * PI); |
| | | let mgLat = lat + dLat; |
| | | let mgLon = lon + dLon; |
| | | return [mgLon, mgLat]; |
| | | } |
| | | } |
| | | }; |