在涉及到经纬度的地方,我们在编码的时候,经常需要进行转换。而且转换代码似乎非常复杂,真正理解,其实需要做一些工作,这里我将了解的经纬度坐标相关的知识做一个梳理。
首先了解经纬度坐标系统。
因为地球并不是一个规则的球体,而是一个椭球体。如何对椭球体进行坐标划定,这里有以下几种:
这里面提到的WGS84,也叫大地坐标系,它是原始坐标系统,为了数据安全和保密,通过地形图非线性保密处理算法(俗称火星加密)加密过的WGS84坐标系,俗称国测局坐标系,或火星坐标系就是我们今天所要提到的GCJ02,目前谷歌地图(中国cn)、腾讯地图、高德地图,使用的都是GCJ02,只有百度地图没有使用这种加密算法,而是使用的是BD09,从名字上可以看出,GCJ02是2002年提出来的算法,BD09则是2009年提出来的,虽然百度地图没有使用GCJ02加密算法,但是他却是在GCJ02基础上做了一个二次加密,所以说,从WGS84坐标系不能直接转BD09,中间需要跨越一个GCJ02,反过来,需要将GCJ02或者BD09转为WGS84就是纠偏算法,相当于逆向解密,同样的DB09直接到不了WGS84,中间还需要转为GCJ02,所以现在的很多算法,如果你看到有百度坐标转大地坐标,基本上需要借助火星坐标来计算。
上面说了坐标转换算法之间的关系,下面来说说具体的算法:
一般,你的算法里面,可能会有如下几个变量:
public static double pi = 3.1415926535897932384626;
public static double a = 6378245.0;
public static double ee = 0.00669342162296594323;
结合上面的表格,你就知道,它采用了何种椭球系数。显然,这是采用了克拉索索夫斯基椭球系数。有的地方称ee是扁率,其实不对,其实是我们上面提到的第一偏心率e的平方。
WGS84转GCJ02的理论公式:
上面公式里面有个经纬度偏移值:{105,35}。这个其实是中华人民共和国大地原点坐标。位置在中国陕西省咸阳市下的泾阳县。
这个公式看着很复杂,最麻烦的在于对经纬度做一个多项式转换。即使很麻烦,但是根据公式,我们通过代码也能一步步