In [97]:
import re
import struct
import datetime

# This gives bullshit
def weeksecondstoutc(gpsweek,gpsseconds,leapseconds):
    import datetime, calendar
    datetimeformat = "%Y-%m-%d %H:%M:%S"
    epoch = datetime.datetime.strptime("1980-01-06 00:00:00",datetimeformat)
    elapsed = datetime.timedelta(days=(gpsweek*7),seconds=(gpsseconds+leapseconds))
    return datetime.datetime.strftime(epoch + elapsed,datetimeformat)

gpsdata = open('data/data.gps','rb').read()
#print(gpsdata, 'x')
for msg in re.findall('\xb5b\r.{33}', gpsdata):
    (ch, flags, count, wnR, wnF, towMsR, towSubMsR, towMsF, towSubMsF, accEst) = \
        struct.unpack('<BBHHHIIIII', msg[5:33])
    print(count,wnR, towMsR)

    #print 1e-6*towMsR,1e-9*towSubMsR

In [ ]:


In [ ]: