当前位置: 代码迷 >> Web前端 >> googlemap经纬度偏移修正算法完美解决方案
  详细解决方案

googlemap经纬度偏移修正算法完美解决方案

热度:330   发布时间:2012-10-08 19:54:56.0
google地图经纬度偏移修正算法完美解决方案

原帖:http://blog.sina.com.cn/s/blog_65bd7eef0100hu5o.html

?

由于受到国家一些法律法规限制,所有的电子地图服务提供商都需要给地图数据加上偏移和加密,因此广受大家关注的google地图也不例外。但是作为一些用惯了免费服务的网友们,始终对google MAP上的标注偏移很头疼,我收集了一些资料,今天在这里和大家一起分享。

所谓的地图数据加密偏移,其实就是用一个偏移算法对地图的经纬度做一个加减偏移量,从而达到与实际地图不一致。这个偏移算法本身是没有什么规律可言的,每家地图服务商都有一套自己的加密偏移算法,既然算法上没有什么规律可言,但是对于算法中生成的偏移量是否有规律呢?这个是可以肯定的,但是偏移算法中生成的偏移量是有规律而言的。偏移量的规律很难得到,要是能拿到这个偏移量,就可以说是破解了某一个地图服务商的地图加密。


那我们怎么解决这个偏差呢!我在网上找到了一个强人破解google map偏移的方法

毕竟做地图开发的都存在卫星图和地图的切换,谁不想让切换以后的地图对上呢!后来我发现google的地图服务,ditu.gogle开始的都没有偏差,maps.google开头的服务就有偏差,我就开始查找google的取偏移量算法,事前我还是图破解google手机版本的数据,没有成功,估计是使用了自己的压缩或加密算法,最后也没有找到规律,后来才尝试破解web版本的不过web版本的接口我对于js不是特别熟悉,所以本次破解放弃了分析源代码的步骤,而是直接采用排除法那就是把地图部分访问的全部地址,一个接一个封杀掉,查找那个起偏移做用的网址 http://ditu.google.cn/maps/vp?spn=0.0,0.0&z=18& vp=39.111195,117.148067 最后找到了就是这个,记住每次测试用清理浏览器缓存哦,使用fixfox的fildder和adblock就够了 然后就是分析每次返回内容的规律,皇天不负有心人啊,我总算是搞定了,下面是我整理的接口

Google 中国地图偏移接口

?

1.????? 接口地址http://ditu.google.cn/maps/vp?spn=0.0,0.0&z=18&vp=39.111195,117.148067

?????? (注:.cn和.com都可以,我用国内服务器就会选择.cn用美国服务器就会选择.com)

2.????? 返回内容中的有效部分:

3.????? (39.111195, 117.148067, 18, [9, -2, 18, -4, 37, -8, 74, -16, 149, -33, 298, -67, 596, -135, 1193, -270])

4.????? Spn参数暂时未知实际意义,但是需要上发spn参数,任意两个小数用逗号分开Vp参数纬经度值,用逗号分隔,z参数为地图缩放级别,无实际意义

5.????? 取回的部分中有效数字为[9, -2, 18, -4, 37, -8, 74, -16, 149, -33, 298, -67, 596, -134, 1192, -268]这个数组总共有8组数字,每两个为一组,分为别从11级到18级的地图和卫星图的偏移像素数量,我们前一组数字精确的等于后一组数字除二,我们为了得到最精确的偏移,故选择第18级的偏移量1193-2701193x方向上精度的偏移像素,-270y方向上维度偏移像素

6.????? 经纬度的偏移转换我们需要江经纬度39.111195117.148067转化为18级像素值25620917 55392414,然后分别加上偏移量-270,1193,然后再转化为经纬度39.11231854918217 117.15446412563324,即位偏移后的经纬度

?

google的地图采用将地球圆表面投影成平面的方式进行贴图

假设zoom=15;

横坐标从左至右像素为0-256*2的15次幂,也就是每增加一级,地图大小横纵坐标加倍,256为一个标准图片的大小

显示-180度到+180的范围,经度越大x越大

纵坐标从上到下像素为0-256*2的15次幂,显示+85到-85度的范围,纬度越小y越大

我们先看一下经度的转换

经度的转换,我就不多说了,均匀分布,大家简单的看一下下面的公式就明白了

经度到像素X值

public static double lngToPixel(double lng, int zoom) {

return (lng + 180) * (256L << zoom) / 360;

}

像素X到经度

public static double pixelToLng(double pixelX, int zoom) {

return pixelX * 360 / (256L << zoom) - 180;

}

纬度到像素Y

public static double latToPixel(double lat, int zoom) {

double siny = Math.sin(lat * Math.PI / 180);

double y = Math.log((1 + siny) / (1 - siny));

return (128 << zoom) * (1 - y / (2 * Math.PI));

}

像素Y到纬度

public static double pixelToLat(double pixelY, int zoom) {

double y = 2 * Math.PI * (1 - pixelY / (128 << zoom));

double z = Math.pow(Math.E, y);

double siny = (z - 1) / (z + 1);

return Math.asin(siny) * 180 / Math.PI;

}

?

维度的这个转换,单纯去理解不是很好理解,我也没有太深入的了解,从公式来看,采用了一种非线性变化,也就是靠近赤道的地方单位像素表现的纬度间距大,越靠近两极越小,可能是因为 经度在靠近两极的方向均匀变化,导致纬度也要进行拉伸,否则 ?靠近两极的地方,地理形状就该发生变化了,总之上面的转换公式大家还是可以研究一下的,google维度的表示范围是-85到+85,这个可以求出来!

上面得到的像素XY是像素的坐标,并非是google地图取地图的那个XY那个XY是把像素所标除以256得到的商,也就是每张图片的大小

?

1 楼 cookieanson 2011-03-23  
google于2010.04把坐标偏移解析做到内部服务器里了,你这个方法已经out了。
2 楼 snow_fish 2012-08-14  
楼主,你好,那个链接,我获取来的是这样的,并没有你说的那个数组啊?
window.GAppFeatures && window.GAppFeatures({cb:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
,traffic:{bounds:[{s:29750000,w:121375000,n:30000000,e:121750000,ix:63}]}
,transit_l:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
,transit:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2},{s:29840644,w:121289062,n:30145127,e:121640625,ix:-1}]}
,ob:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
,bike:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
,rmi:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
,mapmaker:{bounds:[{s:29917407,w:121634174,n:29917407,e:121634174,ix:-2}]}
},16);
  相关解决方案