简单介绍下anttna.com提供的两个纠偏/反纠偏接口

(1) 纠偏接口

http://www.anttna.com/goffset/goffset1.php?lat=xxx&lon=xxx

其中，lat为点分十进制纬度值，lon为点分十进制经度值。

返回数据格式为:纠偏后的纬度,纠偏后的经度

通过该接口，可以将GPS设备获取到的实际经纬度值准确地标记在电子地图上。

示例:实际经纬度为22.559151,113.943581 ,通过查询 http://www.anttna.com/goffset/goffset1.php?lat=22.559151&lon=113.943581 可得到纠偏后的经纬度 22.556129054934,113.94844652715。

(2)反纠偏接口

http://www.anttna.com/goffset/goffset2.php?lat_offset=xxx&lon_offset=xxx

其中，lat_offset为点分十进制纠偏后纬度值，lon_offset为点分十进制纠偏后的经度值

返回数据格式为:实际纬度,实际经度

通过该接口，可以将纠偏后的经纬度转换为实际经纬度。例如，将google maps 地图中通过google map api取到的经纬度值转换为实际经纬度。

示例:纠偏后经纬度为22.556129054934,113.94844652715,通过查询 http://www.anttna.com/goffset/goffset2.php?lat_offset=22.556129054934&lon_offset=113.94844652715 可得到实际经纬度22.559151,113.943581

ant做了张google”基站”和移动基站位置的对比图，这张图可以很清晰的反映出两者之间的关系：
注意:本文涉及的所有数据均指移动GSM基站数据。CDMA,3G不在本文讨论范围内。

点击图片查看大图

图中蓝色点是移动基站实际的位置。这个点是怎么得来的呢？hoho，很简单，目视！
基站有3个扇区，每个扇区120度，分别指向正北，东南，西南三个方向，各扇区的cellid值不相同。
红色点是通过google查询lac,cellid得到的经纬度。google认为一个cellid就是一个基站，所以3个扇区就有3个经纬度。把三个红色点连成三角形，你会发现实际的基站位置大概就在三角形的中心附近。

从这张图上也可以看出，Google把一个基站分成了3个”基站”，”基站”的区域范围和移动的单基站比缩小了很多。所以，ant认为Google的基站数据更适合于LBS应用。

新版谷歌地图的偏移原理
以前的谷歌地图(ditu.google.cn)，电子地图和卫星地图图片都和mas.google.com一样。用户在使用浏览器访问谷歌地图时，客户端浏览器把要显示的电子地图和卫星地图下载到本地，再调用google的偏移接口获取偏移值，用javascript进行地图纠偏。一切操作都是在客户端浏览器进行的。现在改进后的谷歌地图就不一样了，谷歌地图在服务器上预先把标准的卫星地图按照偏移量切割成新的带偏移的图片，用户在使用浏览器访问谷歌地图时，客户端浏览器直接下载偏移的地图图片，在和谷歌地图的交互过程中不传输任何偏移值。

新版谷歌地图的纠偏方法
谷歌地图的这次升级很完美，以前大家普遍采用的纠偏方法都失效了。但是，任何系统都有漏洞，关键就看怎么去挖掘。现在的切入点就是，谷歌地图升级后，同一经纬度的卫星图片 ditu.google.cn和maps.google.com不一样。
我们看看对于经纬度(22.535076,114.021229)，通过ditu.google.cn和maps.google.com访问得到zoom 14级别卫星图片的差别。关于经纬度和google maps的图片编号x，y转换关系这里就不叙述了。
用公式计算出(22.535076,114.021229)在zoom 14下的图片编号为x=13381,y=7138
获取ditu.google.cn图片:
http://mt1.google.cn/vt/lyrs=s@58&gl=cn&x=13381&y=7138&z=14&s=Gali
获取maps.google.com图片:
http://khm1.google.com/kh/v=58&x=13381&y=7138&z=14&s=Gali

同一个点在cn地图和com地图上的位置是不一样的，两个位置坐标的差值就是偏移量!

下面就重点研究如何获取这个差值。单纯取一个像素点来算会有很大误差，而且图片中的一个正确点经过偏移后，很可能已经“跑”出这张256×256图片的范围。ant对以前通过接口获取的偏移值统计分析发现，zoom 18级别地图的X方向的偏移<2000,Y方向<1000，也就是16×8 张图片的范围。
为了计算方便和提高精度，我们用zoom 14下整张图片(256×256)做对比的方法来计算偏移，因为zoom 14和zoom 18大小相差16倍，通过计算可知，zoom 14下的偏移范围是X<125,Y<63。所以在zoom 14级别下可以保证一张标准卫星图片经过偏移后不会“跑”出以该图片为中心的9张 (3×3)图片范围。
(1)首先读取ditu.google.cn下(22.535076,114.021229)经过偏移后的卫星地图(图片a) ，x=13381,y=7138

(2)读取maps.google.com下以坐标点(13381,7138)为中心的9张图片，分别是
(13380,7137),(13381,7137),(13382,7137),
(13380,7138),(13381,7138),(13382,7138),
(13380,7139),(13381,7139),(13382,7139)

合并成一张图片(图片b)

(3)剩下的事就是找出有偏移的图片a在图b中的位置(x1,y1),标准图片13381-7138在图b中的实际位置是(256,256),所以偏移值就是X=x1-256,Y=y1-256

纠偏算法验证
ant用vc写了个example来验证算法的正确性，感谢虫牙童鞋在算法上的帮助。给定经纬度后，程序计算该经纬度对应的zoom 18级的偏移值。

(22.535076,114.021229)实际偏移值是:938,568，测试程序计算出的结果是942,566。偏差x方向4个像素，y方向2个像素。

ant又找了10个不同的点，这些点都保存了以前通过google接口查询的标准值，用来验证算法的正确性。

从结果看，偏移值误差大部分在4个像素内，也就是说在google maps卫星地图最高的清晰度级别(zoom 18)时，算法的误差大概是1～2米。

如何快速获取谷歌电子地图的偏移

上面提到的方法挺复杂，其实最简单的方法就是使用anttna.com的偏移接口: http://www.anttna.com/goffset/goffset1.php?lat=xxx&lon=xxx

虽然google接口关闭了，但是ant的goffset接口仍然有效。
goffset的具体使用方法见ant的上一篇blog 完美解决Google电子地图偏移问题 。

注意:2010/03/31 因Google更换了API，本文所提的方法已经失效。

请使用本站的偏移修正接口： http://www.anttna.com/goffset/goffset1.php?lat=纠偏前的纬度值&lon=纠偏前的经度值

Google电子地图偏移的原理

google 地图在显示指定经纬度(lat,lon)的坐标点(x,y)点时，根据偏移函数算出该点的偏移值(a,b)，把该点的坐标减去偏移值后的结果(x-a,y-b)显示在电子地图上。

Google电子地图偏移的修正方法

知道了原理后，修正的方法就很简单了：用同样的方法获取该点的偏移值(a,b)，预先给定位点加上这个偏移值，把结果(x+a,y+b)提交给google map,google map再经过上面的算法处理后(x+a-a,y+b-b)，显示在地图上的就是(x,y)了。

如何得到该点的偏移值(a,b)?

网上已经有人研究出了google获取偏移值的方法，接口地址为 http://ditu.google.cn/maps/vp?spn=0.0,0.0&z=18&vp=xxxx,xxxx ，只要提交同样的经纬度到该接口，返回的结果里就包含了该经纬度的偏移值。

为了简化操作，ant封装了该接口，可以提交经纬度到接口 http://www.anttna.com/goffset/goffset.php?lat=xxx&lon=xxx 直接获得偏移值。返回格式为: zoom,offset_x,offset_y

zoom: 这个值在做坐标运算时需要用到。

offset_x: 指定经纬度X方向的偏移值

offset_y:指定经纬度Y方向的偏移值

有了这三个值就可以在google map里修正坐标了

function convertCoord(vertex,zoom,offsetX,offsetY)
 {
 var projection = gMap.getCurrentMapType().getProjection();
var point = projection.fromLatLngToPixel(vertex, zoom);
 point = new GPoint(point.x + offsetX, point.y + offsetY);
 vertex = projection.fromPixelToLatLng(point, zoom);
 return vertex;
 }
var point = new GLatLng(lat,lng);//原经纬度点
//通过接口http://www.anttna.com/goffset/goffset.php?lat=xxx&lon=xxx获得的数据
var zoom=18;
var offset_x=100;
var offset_y=200;
//修正经纬度
point=convertCoord(point,zoom,offset_x,offset_y);
//

效果测试

之前和zhjtan讨论过关于google map的偏移问题，这里仍然用该经纬度(22.535076,114.021229)做测试。

未修正的地图（电子地图和卫星地图不一致）

下面是修正后的地图（电子地图和卫星地图一致）

从测试效果来看，使用该修正方法可以完美的解决google地图偏移的问题。

电子地图

卫星地图

怀着极大的好奇心，对google maps的数据进行了简单的分析。为了减少GSM，GPRS网络对分析的影响，打开了手机的飞行模式。重设google maps的数据，接上usb通过active sync上网。打开wireshark抓包。。

打开飞行模式


wireshark抓到的数据包

64.233.167.147是google的服务器，google maps在打开的时候和/glm/mmap做了几次数据的交互，分析之前先baidu,google一下“/glm/mmap”，免得做无用功

http://bbs.godeyes.cn/Announce/Announce.asp?BoardID=100&ID=274997 帖子里提到了google maps会把基站id作为参数提交给google服务器，通过基站id得出当前所在的经纬度。
我手机当前获得的基站id是myl:x:x:16647:42553，转换成16进制
16647->0×4107 , 42553->0xA639
在post到/glm/mmap的数据里找找，果然有41 07 和 a6 39 .. 的字样

推断：
1.google最初并没有基站位置信息（中国移动不可能提供），
2.通过全球数以万计使用google maps GPS定位功能的用户搜集到大部分国家的GSM基站id和GPS经纬度，构建成一个庞大的对应数据库，实时更新。
3.有了google的数据库，在没有gps功能的手机设备上也可以根据基站信息精确的定位手机使用者。