首先需要知道GPS的坐标系。 GPS坐标系符合WGS-84标准,在此标准下GPS芯片可发行不同的分组格式。 根据该数据帧报头,GPS数据可以分类为GPGGA、GPGSA、GPGSV、GPRMC等。 这些标头标识后续帧中数据的结构。 通常,诸如纬度、经度、速度和时间之类的感兴趣的定位数据能够从GPRMC帧获取。
这次不说明具体的帧的格式,但是可以在网上简单地检索。 您也可以在总芯片组的接口文档中找到数据帧格式。 通过解密GPRMC软件包,可以获得“3040.8639,n,10405.7573,e”等位置数据。 该字段的含义在文档中。 在此分析的结果为北纬3040.8639 ‘、东经1045.7573 ‘。
得到这两个数据后,试着在地图上确定位置吧。 我发现很多地图的位置都使用小数形式的坐标。 在此需要说明如何将GPS的原始数据转换为小数形式的数据。
示例: 10614’15 ‘的转换
因为度分分秒秒都是60进制
可以这样转换。
15/60=0.25分
(14 0.25 )/60=0.2375度
106 0.2375=106.2375度
所以最后的结果是106.2375
通过这个变换后,我的位置坐标【北纬3040.8639 ‘,东经1045.7573 ‘】变成了【30.681065N,104.095955E】(这就是WGS-84的坐标)。 此坐标是GPS的物理定位,根据国际标准,此坐标需要进行GCJ-02偏移转换。 转换后的坐标可以在谷歌地图、黄金地图和腾讯地图上定位。 (上述三家公司均符合GCJ-02加密标准。 为什么要这么做,是因为要加密坐标。 因为GCJ-02是不可逆的变换。
这里记载了从标准导WGS-84标准转换为GCJ-02的C++源代码。
http://www.Sina.com/viewplaincopyconstdoublepi=3.14159265358979324; //krasovsky 1940///a=6378245.0,1/f=298.3//b=a* (1- f )/ee=(a^2-b^ ) )/a^2; constdoublea=6378245.0; constdoubleee=0.00669342162296594323; staticbooloutofchina(doublelat,doublelon ) if ) lon72.0
04 || lon > 137.8347) return true; if (lat < 0.8293 || lat > 55.8271) return true; return false; } static double transformLat(double x, double y) { double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * sqrt(abs(x)); ret += (20.0 * sin(6.0 * x * pi) + 20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0; ret += (20.0 * sin(y * pi) + 40.0 * sin(y / 3.0 * pi)) * 2.0 / 3.0; ret += (160.0 * sin(y / 12.0 * pi) + 320 * sin(y * pi / 30.0)) * 2.0 / 3.0; return ret; } static double transformLon(double x, double y) { double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * sqrt(abs(x)); ret += (20.0 * sin(6.0 * x * pi) + 20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0; ret += (20.0 * sin(x * pi) + 40.0 * sin(x / 3.0 * pi)) * 2.0 / 3.0; ret += (150.0 * sin(x / 12.0 * pi) + 300.0 * sin(x / 30.0 * pi)) * 2.0 / 3.0; return ret; } /* 参数 wgLat:WGS-84纬度wgLon:WGS-84经度 返回值: mgLat:GCJ-02纬度mgLon:GCJ-02经度 */ void gps_transform( double wgLat, double wgLon, double& mgLat, double& mgLon) { if (outOfChina(wgLat, wgLon)) { mgLat = wgLat; mgLon = wgLon; return; } double dLat = transformLat(wgLon – 105.0, wgLat – 35.0); double dLon = transformLon(wgLon – 105.0, wgLat – 35.0); double radLat = wgLat / 180.0 * pi; double magic = sin(radLat); magic = 1 – ee * magic * magic; double sqrtMagic = sqrt(magic); dLat = (dLat * 180.0) / ((a * (1 – ee)) / (magic * sqrtMagic) * pi); dLon = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * pi); mgLat = wgLat + dLat; mgLon = wgLon + dLon; }
经过以上的转换,我们的坐标就可以在遵循GCJ-02标准的地图删定位了,这里给出大家一个比较好用的测试定位的网址:http://www.gpsspg.com/maps.htm
百度地图有些特殊,它在GCJ-02基础上又进行了一次加密,百度把这个加密标准叫做BD-09,但是这个加密并没有对外公开。但是万能的网络以及无所不能的网友总会有办法解决这类问题。这里我就贴出一个可行的转换办法,经过测试,定位很准确。
[cpp] view plain copy #include <math.h> const double x_pi = 3.14159265358979324 * 3000.0 / 180.0; //将 GCJ-02 坐标转换成 BD-09 坐标 void bd_encrypt(double gg_lat, double gg_lon, double &bd_lat, double &bd_lon) { double x = gg_lon, y = gg_lat; double z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi); double theta = atan2(y, x) + 0.000003 * cos(x * x_pi); bd_lon = z * cos(theta) + 0.0065; bd_lat = z * sin(theta) + 0.006; } void bd_decrypt(double bd_lat, double bd_lon, double &gg_lat, double &gg_lon) { double x = bd_lon – 0.0065, y = bd_lat – 0.006; double z = sqrt(x * x + y * y) – 0.00002 * sin(y * x_pi); double theta = atan2(y, x) – 0.000003 * cos(x * x_pi); gg_lon = z * cos(theta); gg_lat = z * sin(theta); }
这样,我们的坐标也可以在百度地图上准确定位了,满心欢喜~~
另外再贴上一个比较容易解释坐标为何需要转换的图例,以方便大家理解:
国际通用的转换做法:
百度的做法:
这也解释了,为什么GPS坐标需要经过转换的原因。感谢贡献这些资料到网上的人们!
参考:
http://m.blog.csdn.net/blog/wildboy2001/39497681
http://www.haodaima.net/art/2441684
http://www.gpsspg.com/maps.htm
http://wenda.haosou.com/q/1365472754066079
http://www.geekcome.com/content-10-1464-1.html