使用Python(3.1或2.6),我试图从GPS接收器生成的二进制数据文件中读取数据.每小时的数据存储在一个单独的文件中,每个文件大约18 MiB.数据文件有多个可变长度记录,但是现在我需要从其中一个记录中提取数据.
我已经有了能够解码头部的程度.我有点说,因为有些数字没有意义,但大多数都没有.花了几天时间(我开始学习使用Python编程),我没有取得进展,所以是时候寻求帮助了.
参考指南为我提供了消息头结构和记录结构.标头长度可变,但通常为28个字节.
HeaderField # Field Name Field Type Desc Bytes Offset1 Sync char Hex 0xAA 1 02 Sync char Hex 0x44 1 13 Sync char Hex 0x12 1 24 Header Lgth uchar Length of header 1 35 Message ID ushort Message ID of log 2 48 Message Lgth ushort length of message 2 811 Time Status enum Quality of GPS time 1 1312 Week ushort GPS week number 2 1413 Milliseconds GPSec Time in ms 4 16RecordField # Data Bytes Format Units Offset1 Header 02 Number of SV Observations 4 integer n/a H *For first SV Observation* 3 PRN 4 integer n/a H 44 SV Azimuth angle 4 float degrees H 85 SV Elevation angle 4 float degrees H 126 C/N0 8 double db-Hz H 167 Total S4 8 double n/a H 24...27 L2 C/N0 8 double db-Hz H 14828 *For next SV Observation* SV Observation is satellite - there could be anywhere from 8 to 13 in view.
这是我试图理解标题的代码:
import structfilename = "100301_110000.nvd"f = open(filename, "rb")s = f.read(28)x, y, z, lgth, msg_id, mtype, port, mlgth, seq, idletime, timestatus, week, millis, recstatus, reserved, version = struct.unpack("<cccBHcBHHBcHLLHH", s)print(x, y, z, lgth, msg_id, mtype, port, mlgth, seq, idletime, timestatus, week, millis, recstatus, reserved, version)
它输出:
b'\xaa' b'D' b'\x12' 28 274 b'\x02' 32 1524 0 78 b'\xa0' 1573 126HeaderField # Field Name Field Type Desc Bytes Offset1 Sync char Hex 0xAA 1 02 Sync char Hex 0x44 1 13 Sync char Hex 0x12 1 24 Header Lgth uchar Length of header 1 35 Message ID ushort Message ID of log 2 48 Message Lgth ushort length of message 2 811 Time Status enum Quality of GPS time 1 1312 Week ushort GPS week number 2 1413 Milliseconds GPSec Time in ms 4 16RecordField # Data Bytes Format Units Offset1 Header 02 Number of SV Observations 4 integer n/a H *For first SV Observation* 3 PRN 4 integer n/a H 44 SV Azimuth angle 4 float degrees H 85 SV Elevation angle 4 float degrees H 126 C/N0 8 double db-Hz H 167 Total S4 8 double n/a H 24...27 L2 C/N0 8 double db-Hz H 14828 *For next SV Observation* SV Observation is satellite - there could be anywhere from 8 to 13 in view.
0 10485760 3545 35358
解决方法:
18 MB应该可以很好地适应内存,所以我只需将整个内容放入一个大的字节串中,其中一个打开(文件,’rb’)为f:data = f.read(),然后执行所有操作在片上“解析”以按记录推进记录.它更方便,并且可能比在文件中从此处和那里进行许多小读取更快(尽管它不会影响下面的逻辑,因为在任何一种情况下,“数据中的当前关注点”始终在移动[ [总是转发,因为它发生]]基于一次几个字节的struct-unpack计算的数量,以找到标题和记录的长度).
给定“记录开始”偏移量,您可以通过查看一个字节(“字段四”,从头部开头偏移3,与记录开头相同)来确定其标头的长度,并查看消息ID(下一个字段) ,2个字节)看看它是否是你关心的记录(所以只需要那3个字节的结构解包就足够了).
无论它是否是您想要的记录,您接下来需要计算记录的长度(要么跳过它要么全部得到);为此,您计算实际记录数据的开始(记录的开始加上标题的长度加上记录的下一个字段(标题后面的4个字节)乘以观察的长度(如果我正确读取,则为32个字节) ).
这样你就可以将要给予struct.unpack的子字符串隔离开来(当你最终到达你想要的记录时),或者只是将标题记录的总长度添加到“记录开始”偏移量,以获得偏移量下一条记录的开头.
来源:https://www.icode9.com/content-1-343101.html联系客服