gis经纬度坐标转换多格式兼容:支持字符串/数组/GeoJSON
格式
1 2 3 4 5 6 | let coordinatesStrReg = /((-*[1][0-9]{0,2}|0)(\.[0-9]{1,6})*),\s{0,2}((-*[1-9][0-9]{0,1}|0)(\.[0-9]{1,6})*)/g str.replace(coordinatesStrReg, (str, $1, $2, $3, $4, $5) => { // lat=$1 lng lat=$4 console.log($1, $4) }) |
代码,https://www.npmjs.com/package/coordtransform2
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 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 | /* eslint-disable */ /** * Created by Wandergis on 2015/7/8. * 提供了百度坐标(BD09)、国测局坐标(火星坐标,GCJ02)、和WGS84坐标系之间的转换 * https://www.npmjs.com/package/coordinate-convert WGS84 GPS全球定位系统使用的坐标系统,GPS设备直接返回的坐标即为WGS84,另外Google地图使用的也是WGS84坐标 GCJ02 由中国国家测绘局制定的地理信息系统的坐标,国内出版的各种地图坐标系统(包括电子地图),必须至少采用GCJ02对WGS84进行首次加密。 BD09 百度在GCJ02的基础上进行的二次加密格式 */ // 定义一些常量 const x_PI = 3.14159265358979324 * 3000.0 / 180.0 const PI = 3.1415926535897932384626 const a = 6378245.0 const ee = 0.00669342162296594323 class CoordsTransformer { constructor (decimal = 6) { this .decimal = decimal } setDecimal (num) { this .decimal = num } toFixed (arr) { let lng = (arr[0] - 0).toFixed( this .decimal) let lat = (arr[1] - 0).toFixed( this .decimal) return [+lng, +lat] } /** * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换 * 即 百度 转 谷歌、高德 * @param bd_lon * @param bd_lat * @returns {*[]} */ bd09togcj02 (bd_lon, bd_lat) { let x_pi = 3.14159265358979324 * 3000.0 / 180.0 let x = bd_lon - 0.0065 let y = bd_lat - 0.006 let z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi) let theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi) let gg_lng = z * Math.cos(theta) let gg_lat = z * Math.sin(theta) return this .toFixed([gg_lng, gg_lat]) } /** * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换 * 即谷歌、高德 转 百度 * @param lng * @param lat * @returns {*[]} */ gcj02tobd09 (lng, lat) { let z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_PI) let theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_PI) let bd_lng = z * Math.cos(theta) + 0.0065 let bd_lat = z * Math.sin(theta) + 0.006 return this .toFixed([bd_lng, bd_lat]) } /** * WGS84转GCj02 * @param lng * @param lat * @returns {*[]} */ wgs84togcj02 (lng, lat) { if (out_of_china(lng, lat)) { return this .toFixed([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 = lat + dlat let mglng = lng + dlng return this .toFixed([mglng, mglat]) } } /** * GCJ02 转换为 WGS84 * @param lng * @param lat * @returns {*[]} */ gcj02towgs84 (lng, lat) { if (out_of_china(lng, lat)) { return this .toFixed([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 = lat + dlat let mglng = lng + dlng return this .toFixed([lng * 2 - mglng, lat * 2 - mglat]) } } } function transformlat (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 } function transformlng (lng, lat) { 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 } /** * 判断是否在国内,不在国内则不做偏移 * @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 ) } export default CoordsTransformer |
代码
转载本站文章《gis经纬度坐标转换多格式兼容:支持字符串/数组/GeoJSON》,
请注明出处:https://www.zhoulujun.cn/html/GIS/WebGIS/8227.html
标签:
GeoJSON
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了