GPS基本解析思路

时间:2025-01-25 07:24:22
世界时区的划分图 全球时区的划分: 每个时区跨15°经度。以0°经线为界向东向西各划出7.5°经度,作为0时区。即0时区的经度范围是7.5°W——7.5°E。从7.5°E与7.5°W分别向东、向西每15°经度划分为一个时区,直到东11区和西11区。东11区最东部的经度是172.5°E,由172.5°E——180°之间就是东12区。西11区最西部的经度是172.5°W,由172.5°W——180°之间就是西12区。东、西12区各占经度7.5°,合成一个完整的时区,即全球总共划分为24个时区。东、西12区钟点相同,日期相差1天,因此180°称为理论上的国际日期变更线。 由于地球的自转运动,不同地区有不同的地方时间,为了解决时间混乱的问题,采取了划分时区的办法。每个时区*经线所在地的地方时间就是这个时区共用的时间,称为区时。在实际应用中各国不完全按照区时来定时间,许多国家制定一个法定时,作为该国统一使用的时间,例如我国使用120°E的地方时间,称为北京时间。 GPS 上电后,每隔一定的时间就会返回一定格式的数据,数据格式为: $信息类型,x,x,x,x,x,x,x,x,x,x,x,x,x每行开头的字符都是'$',接着是信息类型,后面是数据,以逗号分隔开。 一行完整的数据如下: $GPRMC,080655.00,A,4546.40891,N,12639.65641,E,1.045,328.42,170809,,,A*60 GPS信息类型: GPGSV:可见卫星信息 GPGLL:地理定位信息 GPRMC:推荐最小定位信息 GPVTG:地面速度信息 GPGGA:GPS定位信息 GPGSA:当前卫星信息 GPRMC 最小定位信息: 数据详解:$GPRMC,<1>,<2>,<3>,<4>,<5>,<6>,<7>,<8>,<9>,<10>,<11>,<12>*hh <1> UTC 时间,hhmmss(时分秒)格式 <2> 定位状态,A=有效定位,V=无效定位 <3>纬度(度分)格式(前面的0也将被传输) <4> 纬度半球N(北半球)或S(南半球) <5>经度(度分)格式(前面的0也将被传输) <6> 经度半球E(东经)或W(西经) <7>地面速率(000.0~999.9节,前面的0也将被传输) <8>地面航向(000.0~359.9度,以真北为参考基准,前面的0也将被传输) <9> UTC 日期,ddmmyy(日月年)格式 <10>磁偏角(000.0~180.0度,前面的0也将被传输) <11> 磁偏角方向,E(东)或W(西) <12>模式指示(仅NMEA01833.00版本输出,A=自主定位,D=差分,E=估算,N=数据无效) 解析内容: 1.时间,这个是格林威治时间,是世界时间(UTC),我们需要把它转换成北京时间(BTC),BTC和UTC差了8个小时,要在这个时间基础上加8个小时。 2. 定位状态,在接收到有效数据前,这个位是‘V’,后面的数据都为空,接到有效数据后,这个位是‘A’,后面才开始有数据。 3. 纬度,我们需要把它转换成度分秒的格式,计算方法:如接收到的纬度是:4546.40891 4546.40891/100=45.4640891可以直接读出45度, 4546.4089145*100=46.40891, 可以直接读出4646.4089146 =0.40891*60=24.5346读出24秒, 所以纬度是:454624秒。 4. 南北纬,这个位有两种值‘N’(北纬)和‘S’(南纬) 5. 经度的计算方法和纬度的计算方法一样 6. 东西经,这个位有两种值‘E’(东经)和‘W’(西经) 7.速率,这个速率值是海里/时,单位是节,要把它转换成千米/时,根据:1海里=1.85公里,把得到的速率乘以1.858. 航向,指的是偏离正北的角度 9. 日期,这个日期是准确的,不需要转换 GPGGA GPS定位数据 数据详解:$GPGGA,<1>,<2>,<3>,<4>,<5>,<6>,<7>,<8>,<9>,M,<10>,M,<11>,<12>*xx<CR><LF> $GPGGA:起始引导符及语句格式说明(本句为GPS定位数据); <1> UTC 时间,格式为; <2> 纬度,格式为(第一位是零也将传送); <3> 纬度半球,N 或S(北纬或南纬) <4> 经度,格式为(第一位零也将传送); <5> 经度半球,E 或W(东经或西经) <6> 定位质量指示,0=定位无效,1=定位有效; <7>使用卫星数量,从0012(第一个零也将传送) <8>水平精确度,0.599.9 <9>天线离海平面的高度,-9999.99999.9米M指单位米 <10>大地水准面高度,-9999.99999.9米M指单位米 <11>差分GPS数据期限(RTCMSC-104),最后设立RTCM传送的秒数量 <12>差分参考基站标号,从00001023(首位0也将传送)。 解析内容: 第910 个字段,海平面高度和大地水准面高度,单位是米 GPVTG 地面速度信息 $GPVTG,<1>,T,<2>,M,<3>,N,<4>,K,<5>*hh <1> 以正北为参考基准的地面航向(000~359度,前面的0也将被传输) <2> 以磁北为参考基准的地面航向(000~359度,前面的0也将被传输) <3> 地面速率(000.0~999.9节,前面的0也将被传输) <4> 地面速率(0000.0~1851.8公里/小时,前面的0也将被传输) <5> 模式指示(仅NMEA0183 3.00版本输出,A=自主定位,D=差分,E=估算,N=数据无效 GPGSV 可视卫星状态 例:$GPGSV,2108063324045103607447162107844173631342*78 标准格式: $GPGSV,(1),(2),(3),(4),(5),(6),(7),…(4),(5),(6),(7)*hh(CR)(LF)  各部分含义为: (1)总的GSV语句电文数;2; (2)当前GSV语句号:1; (3)可视卫星总数:08; (4)PRN码(伪随机噪声码) 也可以认为是卫星编号 (5)仰角(00~90度):33; (6)方位角(000~359度):240; (7)信噪比(00~99dB):45dB(后面依次为第101617号卫星的信息); *总和校验域;  hh 总和校验数:78; (CR)(LF)回车,换行。 注:每条语句最多包括四颗卫星的信息,每颗卫星的信息有四个数据项,即:   (4)-卫星号,(5)-仰角,(6)-方位角,(7)-信噪比。 例: $GPGSV,3,1,10,24,82,023,40,05,62,285,32,01,62,123,00,17,59,229,28*70 每条语句包含四部分内容,例如:第一部分是“24,82,023,40”,第二部分是“05,62,285,32”等等。 每部分的第一个词为PRC,第二个词为卫星高程,跟着为方位角和信号强度。 这个语句里最重要的指标应该算是“信号躁声比(signal-to-noise ratio)”(以下简称为SNR)。 这个数值标示卫星信号的接收率。我们知道,卫星是以相同的强度发射信号,但是传播过程中难免会遇到诸如树和墙之类的 障碍物,这样就影响了信号的识别。 典型的SNR值在050之间,其中50表示非常好的信号。(SNR可以达到99)。 GPGSA 当前卫星信息 例:$GPGSA,A,3,01,20,19,13,,,,,,,,,40.4,24.4,32.2*0A   字段0:$GPGSA,语句ID,表明该语句为GPS DOP and Active Satellites(GSA)当前卫星信息 字段1:定位模式,A=自动手动2D/3D,M=手动2D/3D 字段2:定位类型,1=未定位,2=2D定位,3=3D定位 字段3:PRN码(伪随机噪声码),第1信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段4:PRN码(伪随机噪声码),第2信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段5:PRN码(伪随机噪声码),第3信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段6:PRN码(伪随机噪声码),第4信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段7:PRN码(伪随机噪声码),第5信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段8:PRN码(伪随机噪声码),第6信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段9:PRN码(伪随机噪声码),第7信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段10:PRN码(伪随机噪声码),第8信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段11:PRN码(伪随机噪声码),第9信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段12:PRN码(伪随机噪声码),第10信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段13:PRN码(伪随机噪声码),第11信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段14:PRN码(伪随机噪声码),第12信道正在使用的卫星PRN码编号(00)(前导位数不足则补0) 字段15:PDOP综合位置精度因子(0.5 - 99.9) 字段16:HDOP水平精度因子(0.5 - 99.9) 字段17:VDOP垂直精度因子(0.5 - 99.9) 字段18:校验值 以MTK平台主机获取到的数据为例: $GPRMC,092927.000,A,2235.9058,N,11400.0518,E,0.000,74.11,151216,,D*49 $GPVTG,74.11,T,,M,0.000,N,0.000,K,D*0B $GPGGA,092927.000,2235.9058,N,11400.0518,E,2,9,1.03,53.1,M,-2.4,M,0.0,0*6B $GPGSA,A,3,29,18,12,25,10,193,32,14,31,,,,1.34,1.03,0.85*31 $GPGSV,3,1,12,10,77,192,17,25,59,077,42,32,51,359,39,193,49,157,36*48 $GPGSV,3,2,12,31,47,274,25,50,46,122,37,18,45,158,37,14,36,326,18*70 $GPGSV,3,3,12,12,24,045,45,26,17,200,18,29,07,128,38,21,02,174,*79 GPRMC的信息以上格式可以解析为: 09h29m27s,有效定位,维度2235.9058(海里/时,单位是节),北半球,经度11400.0518,东经,地面速率0.000,地面航向74.11(偏离正北的角度),UTC日期151216年,磁偏角,磁偏角方向, 部分解析代码 以下为网友提供的部分解析代码可供参考: 创建一个GPS数据结构体: typedef data struct{ double latitude; //经度 double longitude; //纬度 int latitude_Degree; //度 int latitude_Cent; //分 int latitude_Second; //秒 int longitude_Degree; //度 int longitude_Cent; //分 int longitude_Second; //秒 float speed; //速度 float direction; //航向 float height; //海拔高度 int satellite; U8 NS; U8 EW; DATE_TIME D; }GPS_INFO; 时间结构体: typedef struct{ int year; int month; int day; int hour; int minute; int second; }DATE_TIME; 核心算法就是解析GPRMC数据,得到经纬度,日期时间,速度,航向: int GPS_RMC_Parse(char *line, GPS_INFO *GPS) { U8 ch, status, tmp; float lati_cent_tmp, lati_second_tmp; float long_cent_tmp, long_second_tmp; float speed_tmp; char *buf = line; ch = buf[5]; status = buf[GetComma(2, buf)]; if (ch == 'C') //如果第五个字符是C,($GPRMC) { if (status == 'A') //如果数据有效,则分析 { GPS->NS = buf[GetComma(4, buf)]; GPS->EW = buf[GetComma(6, buf)]; GPS->latitude = Get_Double_Number(&buf[GetComma(3, buf)]); GPS->longitude = Get_Double_Number(&buf[GetComma(5, buf)]); GPS->latitude_Degree = (int)GPS->latitude / 100; //分离纬度 lati_cent_tmp = (GPS->latitude - GPS->latitude_Degree * 100); GPS->latitude_Cent = (int)lati_cent_tmp; lati_second_tmp = (lati_cent_tmp - GPS->latitude_Cent) * 60; GPS->latitude_Second = (int)lati_second_tmp; GPS->longitude_Degree = (int)GPS->longitude / 100; //分离经度 long_cent_tmp = (GPS->longitude - GPS->longitude_Degree * 100); GPS->longitude_Cent = (int)long_cent_tmp; long_second_tmp = (long_cent_tmp - GPS->longitude_Cent) * 60; GPS->longitude_Second = (int)long_second_tmp; speed_tmp = Get_Float_Number(&buf[GetComma(7, buf)]); //速度(单位:海里/时) GPS->speed = speed_tmp * 1.85; //1海里=1.85公里 GPS->direction = Get_Float_Number(&buf[GetComma(8, buf)]); //角度 GPS-> = (buf[7] - '0') * 10 + (buf[8] - '0'); //时间 GPS-> = (buf[9] - '0') * 10 + (buf[10] - '0'); GPS-> = (buf[11] - '0') * 10 + (buf[12] - '0'); tmp = GetComma(9, buf); GPS-> = (buf[tmp + 0] - '0') * 10 + (buf[tmp + 1] - '0'); //日期 GPS-> = (buf[tmp + 2] - '0') * 10 + (buf[tmp + 3] - '0'); GPS-> = (buf[tmp + 4] - '0') * 10 + (buf[tmp + 5] - '0') + 2000; UTC2BTC(&GPS->D); return 1; } } return 0; } line是串口接收的一行数据buf GetComma函数作用是一行数据中第几个逗号后面那个字符在这行数据中的位置 Get_Double_Number函数作用是把给定字符串第一个逗号之前的字符转化成双精度型,在这里就是把代表经度和纬度的字符串转换成数字,同样的函数还有Get_Float_Number UTC2BTC函数是将世界时间转换成北京时间(相差8小时) 在LCD显示程序中把GPS_INFO结构体的已经被赋值的变量显示到屏上相应的位置即可 还有一个GPGGA信息段可以提供海拔高度和卫星数量信息 int GPS_GGA_Parse(char *line, GPS_INFO *GPS) { U8 ch, status; char *buf = line; ch = buf[4]; status = buf[GetComma(2, buf)]; if (ch == 'G') //$GPGGA { if (status != ',') { GPS->height = Get_Float_Number(&buf[GetComma(9, buf)]); GPS->satellite = Get_Int_Number(&buf[GetComma(7, buf)]); return 1; } } return 0; }