|
# S ^6 B; N# F7 Y: l l 为什么要进行坐标系转换?
6 z4 O3 [* Q1 T- V: m# f 因为国内对地理坐标有特殊的政策,GPS获取的原始坐标不能在互联网上对外提供,所有电子地图必须对位置做偏移处理,这就导致从底层接口得到的坐标展示在地图上会有偏移。这种偏移不是线性的,不能简单地通过加减来校正,并且不同的地图提供商采用的算法也不一样。不同的坐标系之间可能有几十到几百米的偏移,所以在地图上绘制矢量图形、做地理数据可视化时,需要修正不同坐标系之间的偏差。 / D& n7 @+ u: e$ ]5 x8 ~4 e
l 关于坐标系
9 X2 v) n2 `% b4 d( f" ? 通常GPS获取的是WGS84坐标体系,包括很多国外地图厂商提供的地图数据也是WGS84坐标体系,但是在我国,由于政策原因,在地图发布和出版的时候,坐标至少需要经过国家测绘局加密,也就是对WGS84坐标进行一次非线性加偏,得到GCJ02坐标系,俗称火星坐标。
8 i1 s! m3 Y# v5 x 另外,国内许多地图公司,由于各种原因,会在火星坐标的基础上进行二次非线性加偏,得到自己的坐标系统,各地图厂商使用坐标系情况如下 。
! A0 u1 `) T. _- C ! O$ I9 i, }4 Y
l 坐标转换方法& Y0 v6 o2 M/ o# J
整理了一下坐标系转换的方法,仅限单个数组转换
' U1 V2 e% ?' K9 o7 p 火星坐标转wgs84
: K- m! a0 f' F6 A. W // 火星坐标转wgs84$ ^( {# Z+ `. z `$ M( H
function gcj02towgs84(lng, lat) {- m, b8 w0 M; ?' z7 y* L
var x_PI = 3.14159265358979324 * 3000.0 / 180.0;5 j' P8 x T( N9 w/ \
var PI = 3.1415926535897932384626;
9 r2 `( \0 ^7 `: X, q3 y var a = 6378245.0;7 v& k, E) V0 X9 P8 _7 C! A' E
var ee = 0.00669342162296594323;6 [) s- W' I5 T0 T# t }$ |
if ((lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)) {3 f/ C( c, I( ~6 U
return [lng, lat]3 |- y5 o6 a8 G, y/ _0 ]2 w
}0 b9 X; Q" N: A6 H
else {
8 [- L/ A s) P1 m% P var dlat = ((lng, lat) => {
3 G/ M. D$ Z! d. T% e, o 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));
% M) a9 `, R0 s ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;! U; n, t- H/ w; j8 U$ `/ x) @
ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;% D8 N; ?& e) ^" b# Q. y1 Z: q
ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
( a, X- p/ _$ o0 o. l. p- j1 U( O return ret+ a( i- {9 _- ?* h( C) {
})(lng - 105.0, lat - 35.0);! y5 C$ U0 V {: ]
var dlng = ((lng, lat) => {
, M7 D' S" w5 } var ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
" {; s1 {) B1 n6 D ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;4 C; B6 A: Q* j
ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;
' D2 _, K4 d. y- z ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;( T/ s' I# p/ L- t8 m: v
return ret! ]0 ~ }* K% U6 g9 [2 {2 }% N; ]' e
})(lng - 105.0, lat - 35.0);) o' O; ?& m4 L7 s
var radlat = lat / 180.0 * PI;4 ~, s: \$ L7 ^, X
var magic = Math.sin(radlat);
* e4 D' h1 @9 t) x7 G0 D( X magic = 1 - ee * magic * magic;
" ~: S, d' N: ^/ s7 S- g var sqrtmagic = Math.sqrt(magic);# x- A, b3 J" X5 o7 v
dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);; G& M7 V: H+ [. b% j7 X! |1 _# `
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI); |# j/ W& `0 L( o: w: z `
mglat = lat + dlat;) t( ]1 D7 o& A7 q$ W
mglng = lng + dlng;8 H4 C2 Z" L$ o( C6 N; L
return [lng * 2 - mglng, lat * 2 - mglat]
7 D5 t+ o% W: t' K( D$ c }
2 e( C1 a5 `3 W1 K6 f$ W5 T } % A( g* I- m1 W9 v2 g! _9 |* @
wgs84转火星坐标 " Q; G; [. _1 P$ E/ A( c+ W
// wgs84转火星坐标% r, H5 Q- R' ]/ k Q& d
function wgs84togcj02(lng, lat) {1 o }0 q" V$ T7 Y4 C3 o
var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
& W' A; n& D; H' ^ I F+ p var PI = 3.1415926535897932384626;
1 `# L+ F9 j- l! _9 c w- W var a = 6378245.0;
# K: D6 b/ L: R1 h7 k4 I E var ee = 0.00669342162296594323;
7 u1 T' s6 R2 b, U( S if ((lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)) {
& h% S" f S' ^) O return [lng, lat]
+ |( l1 ?: O- _* t0 G' O } else {
( \+ s5 n+ Z: ]% S5 N5 x; } var dlat = ((lng, lat) => {
! v4 p, D% T, M( ]: u 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));
0 k* P1 S! l- y) H6 A& V% t: v ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;/ z5 n; i, ^' f
ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;
( r: X" o1 Q7 c ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
( B/ L) { r3 k6 j return ret
! C% J) ^: [, z6 c7 c4 ^ })(lng - 105.0, lat - 35.0);" v" v: u7 u" I$ l' o
var dlng = ((lng, lat) => {3 D6 t+ G7 u% k5 g
var ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));4 S& |8 N* a* @, k) I9 R
ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;' B- t, d$ Q8 _) c
ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;: r: j u3 r) {9 D, I
ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;
~* s& Q1 s4 D0 u$ N1 g0 H8 m) d return ret
8 X0 ?& u9 d/ J C, w }" P5 X7 R1 J! \/ t: p) a& k, b
)(lng - 105.0, lat - 35.0);" `+ b/ \2 ~2 d) J# \$ R
var radlat = lat / 180.0 * PI;7 b4 c: M* [; l, U9 P4 k$ }$ k0 w
var magic = Math.sin(radlat);
) s% m4 Q2 i0 \6 `* v" K magic = 1 - ee * magic * magic;
' C* l6 N+ w7 F, ?; `7 n4 ~ var sqrtmagic = Math.sqrt(magic);
7 f2 B0 B; P' O1 j dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);; o% W/ ~6 b& U' K7 [" f6 q' s! B
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
1 m/ i6 F4 ]+ ~ var mglat = lat + dlat;
- `# {! \3 m! v, K8 \' P9 M var mglng = lng + dlng;: o$ ^. R% X/ Y9 j' [- C
return [mglng, mglat]
3 Y% i, N4 t+ w }
3 J' W5 c4 c7 z i } - Y8 `% z6 C/ ?4 b; Q
百度转火星坐标
# P9 J' Y0 b4 Y // 百度转火星坐标
5 l+ i' \6 @* i& y7 O# I const bd2xh = (bd_lng, bd_lat) => {" y* n b0 L, L/ h& n5 ]
const x_pi = Math.PI * 3000.0 / 180.0;
: W4 R7 b) o) X const x = bd_lng - 0.0065;
! t0 w4 E& {; s: F1 A const y = bd_lat - 0.006;. b: o$ [0 j0 M9 n$ J$ a! Q
let z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
$ a/ K" c1 _# y0 [- L4 u let theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);$ d4 ~! r# N* l1 l, e
return [z * Math.cos(theta), z * Math.sin(theta)];
) X. U2 m1 C# R" T' E( M! U }
$ K$ B$ L8 f7 o$ F) Q5 V: ~& { l 坐标转换工具 8 P2 z. N9 _% [, g
有的静态数据不需要动态转换也有现成的在线坐标系转换工具可以直接拿来用 + d/ Y5 d+ Y1 R) d5 \9 N
( M! p+ y8 z; n, i: X* S+ ^ 支持常见百度地图(BD09)、火星坐标系(GCJ02)、WGS84坐标系之间的相互转换 ( k+ `' A( S+ V" b: Y+ a3 j9 C
具体支持的坐标系转换如下
" c* a ?9 Y6 g- X. C6 F' b 1)WGS-84 → GCJ-02
6 ~0 n E2 ~) V; s( q 2)WGS-84 → BD-09
9 C8 R" I8 p$ A# _0 m 2)GCJ-02 → WGS-84
! ^ N1 P0 ]* Q' U$ Q" v' C! |' ~+ n d 3)GCJ-02 → BD-09
5 H7 ^5 n" M3 n# s7 g 4)BD-09 → GCJ-02 c" n7 j& P% @& M4 z0 K: X, g
5)BD-09 → WGS-84
- V! O" \, B$ \ 支持单个数组、不包含嵌套关系的对象、标准GeoJSON
* X! k; r9 M N6 ~2 C" v5 K; x 这个工具也还可以获取行政边界数据、地图数据模拟
* w: n/ I/ g* h. a, b& W; a+ F2 E# n0 j 7 m \: a. S- n/ \* A+ A! @
关注【官方微信公众号】,查看更多可视化资讯~等你呦1、EasyV数字孪生可视化2、易知微数字孪生世界 EasyV数字孪生可视化平台免费试用易知微官网,申请EasyV免费试用,个人版已上线,立即去免费体验! " _9 c& a& d b4 ?3 W1 P4 C/ O
免费试用期间:运营、产品、设计、技术人员专业一对一拉群指导,解答疑惑。
# F3 K: C- E5 ]4 H, g- |# i$ J" }( `$ h8 @1 P. d& W9 T+ r
" V7 \% {; I' {, Z; z* V8 n, H4 R# Y% e4 C
( i6 ]8 b% b: S2 v
|