Google地球(GPS)坐标之地图坐标偏移

转自：http://www.xkshow.cn/?p=33

这个问题涉及国家的法律问题，很重要。根据我国相关法律，在国内发行的任何民用中国地图产品不准使用真实的坐标，必须加入一定的偏移量(大约600米)，据说这一做法是“出于国家安全方面的考虑”。

偏移数据是不公开的，五六百米是实测结果。我国法律条文规定在国内必须使用国家测绘局的坐标体系(西安80)，而国际通行的坐标体系是WGS84。另外，民航使用的坐标系是国际标准WGS84，否则国际航班没法正常飞行。

Google Maps就属于在国内发行的民用中国地图产品，自然也不例外地加入了偏移量。所有在中国市场上能买到的GPS中国地图(比如灵图，凯立德，Route66和Garmin等)也都带有600米左右的偏差，其显示的坐标与真实坐标相差600米，没有任何实用价值.

例如，在软件里查看GPS坐标，同一地点的数据如下：

官方地图坐标 N22.540938666666665，E113.95289433333335

偏移地图坐标 N22.5383656877113157，E113.9581618420152

GPS数据查看器显示的坐标与纠偏地图坐标相同，官方地图与纠偏地图误差大约有500米。（可以简单的使用偏移差值计算，但并非所有的确均相同）

Google Map的中国地图，坐标做过偏移，没法直接和GPS得来的经纬度数据匹配使用，得要在相关行业管理部门取得加密插件，植入软件里，才能正确显示位置。

国家保密插件，也叫做加密插件或者加偏或者SM模组，其实就是对真实坐标系统进行人为地加偏处理，按照几行代码的算法，将真实的坐标加密成虚假的坐标，而这个加偏并不是线性的加偏，所以各地的偏移情况有所不同，而加密后的坐标也常被人称为火星坐标系统，让国外的政府和企业等，都看不懂我国坐标系统。

现在，所有电子地图所有导航设备，都需要加入国家保密插件。第一步，地图公司测绘地图，测绘完成后，送到国家测绘局，将真实坐标的电子地图，加密成“火星坐标“，这样的地图才是可以出版和发布的，然后才可以让GPS公司处理。第二步，所有GPS公司，只要需要汽车导航的，需要用到导航电子地图的，统统需要在软件中加入国家保密算法，将Com口读出来的真实的坐标信号，加密转换成国家要求的保密的坐标，这样，GPS导航仪和导航电子地图就可以完全匹配，GPS也就可以正常工作。

所以，所有导航电子地图公司，像四维图新、瑞图、易图通、高德等，都需要将自己的电子地图拿到国家测绘局进行加密处理，而所有导航软件公司，都需要将自己的导航软件中加入国家测绘局提供的加密算法的代码，而这一段代码，就是国家的机密，不是公开的，每次去国家测绘局加密处理都需要预约并在封闭的环境中进行，编译完成后需要是在主程序的exe中，而不能编译在外部dll等文件中。

当然，并不是说这样就没有使用GPS的可能了。通过自定义GmapType(ditu)和自定义投影坐标Gprojection(dituProjection)来实现偏差纠正，原理是修改ditu层的投影换算函数:

.fromLatlngToPixel()+offset

.fromPixelToLatlng()-offset

这样就无须修改上层Overlay的坐标参数，直接使用正确的Lat,Lng即可。即通过重新调整地图的投影，从而实现地图的正确显示。

附：

自定义投影坐标(GProjection)来实现Google Map中国地图和卫星图的偏差纠正

原理:

通过自定义GMapType(ditu)和自定义投影坐标GProjection(dituProjection) 来实现偏差纠正。 

原理是修改ditu层的投影换算函数，
在.fromLatLngToPixel()里面 +Offset
在.fromPixelToLatLng()里面 -Offset

这样子就无需修改上层Overlay的坐标参数，直接使用正确的Lat,Lng即可

代码:(以上海为例)

function load_gmap()
{

    //shanghai
    var xOffset= -0.001889737;
    var yOffset= 0.004844069;

    var map = null;


    if(GBrowserIsCompatible()) {
	

	var Cditu = new GCopyright(1, 
				   new GLatLngBounds(
						     new GLatLng(-90,-180),
						     new GLatLng(90,180) ),
				   0,
				   "Mapabc.com"); 
	
	var copyright = new GCopyrightCollection("");
	
	copyright.addCopyright(Cditu);
	var dituTileLayer = new GTileLayer(copyright, 1, 17);
	
	dituTileLayer.getTileUrl = function(tile, zoomlevel,x) {
	    var url = G_NORMAL_MAP.getTileLayers()[0].getTileUrl(tile,zoomlevel,x);
	    var bits = url.split("&");
	    bits[1] = "http://servicetile.mapabc.com/googlechina/maptile?" + bits[1];
	    /**bits[1] = "http://mapgoogle.mapabc.com/googlechina/maptile?" + bits[1];**/
	    bits.shift();
	    url = bits.join("&");
	    return url; 
	};

	
	function dituProjection(xOffset,yOffset){
	    this.xOffset = xOffset;
	    this.yOffset = yOffset;
	}

	dituProjection.prototype = new GProjection();

	dituProjection.prototype.fromLatLngToPixel = function(latlng,  zoom){
	    return (G_NORMAL_MAP.getProjection()).fromLatLngToPixel(new GLatLng(latlng.lat()+this.xOffset,latlng.lng()+this.yOffset),zoom);
	};

	dituProjection.prototype.fromPixelToLatLng = function(pixel,zoom,unbounded)
	{
	    var latlng = (G_NORMAL_MAP.getProjection()).fromPixelToLatLng(pixel,zoom,unbounded);
	    return new GLatLng(latlng.lat()-this.xOffset,latlng.lng()-this.yOffset);
	}

	dituProjection.prototype.tileCheckRange = function(tile, zoom, tilesize)
	{
	    return (G_NORMAL_MAP.getProjection()).tileCheckRange(tile,zoom,tilesize);
	}
	dituProjection.prototype.getWrapWidth = function(zoom)
	{
	    return (G_NORMAL_MAP.getProjection()).getWrapWidth(zoom);
	}
	

	var ditu = new GMapType([dituTileLayer],
				new dituProjection(xOffset,yOffset), 
				"Ditu", 
				{ shortName: "ditu", alt: "layer from ditu.google.com" }
				);


	map = new GMap2(document.getElementById('gmap'),{mapTypes:[ditu,/*G_NORMAL_MAP,*/G_SATELLITE_MAP]});
	map.addControl(new GLargeMapControl());
	/**map.addControl(new GMapTypeControl());**/
	map.addControl(new GOverviewMapControl());
	map.addControl(new GScaleControl());

	map.enableDoubleClickZoom();
	/**map.enableScrollWheelZoom();**/
	
	map.setCenter(new GLatLng(31.231628368031693, 121.47645235061645 ), 16);

	ditu.getName = function(){ return '中国地图';};
	ditu._getName = function(){ return 'ditu';};
	/*
	G_NORMAL_MAP.getName = function(){ return '世界地图';};
	G_NORMAL_MAP._getName = function(){ return 'Map';};
	*/
	G_SATELLITE_MAP.getName = function(){ return '卫星图';};
	G_SATELLITE_MAP._getName = function(){ return 'Satellite';};
	
	/** the map type **/
	var _mt = ditu;

	map.setMapType(_mt);

	map.addControl(new GMapTypeControl());

	var point = new GLatLng(31.231628368031693, 121.47645235061645);

	var marker = new GMarker(point);
	GEvent.addListener(marker,"click",function(){marker.openInfoWindowHtml('Text Info Here!');});

	map.addOverlay(marker);

    }



}

囧。。。貌似这个方法也解决不了非线性的问题。