racerxdl / channelmuxer.py

#!/usr/bin/env python import sys, struct, os, packetmanager FRAMESIZE = 892 M_PDUSIZE = FRAMESIZE - 6 EXPORTCORRUPT = False SEQUENCE_FLAG_MAP = { 0: "Continued Segment", 1: "First Segment", 2: "Last Segment", 3: "Single Data" } def ParseMSDU(data): o = struct.unpack(">H", data[:2])[0] version = (o & 0xE000) >> 13 type = (o & 0x1000) >> 12 shf = (o & 0x800) >> 11 apid = (o & 0x7FF) o = struct.unpack(">H", data[2:4])[0] sequenceflag = (o & 0xC000) >> 14 packetnumber = (o & 0x3FFF) packetlength = struct.unpack(">H", data[4:6])[0] -1 data = data[6:] return version, type, shf, apid, sequenceflag, packetnumber, packetlength, data def CalcCRC(data): lsb = 0xFF msb = 0xFF for c in data: x = ord(c) ^ msb x ^= (x >> 4) msb = (lsb ^ (x >> 3) ^ (x << 4)) & 255 lsb = (x ^ (x << 5)) & 255 return (msb << 8) + lsb def CheckCRC(data, crc): c = CalcCRC(data) if not c == crc: print " Expected: %s Found %s" %(hex(crc), hex(c)) return c == crc def SavePacket(channelid, packet): global totalCRCErrors global totalSavedPackets try: os.mkdir("channels/%s" %channelid) except: pass if packet["apid"] == 2047: print " Fill Packet. Skipping" return filename = "channels/%s/%s_%s.lrit" % (channelid, packet["apid"], packet["version"]) print "- Saving packet to %s" %filename if packet["framesdropped"]: print " WARNING: Some frames has been droped for this packet." datasize = len(packet["data"]) if not datasize - 2 == packet["size"]: # CRC is the latest 2 bytes of the payload print " WARNING: Packet Size does not match! Expected %s Found: %s" %(packet["size"], len(packet["data"])) if datasize - 2 > packet["size"]: datasize = packet["size"] + 2 print " WARNING: Trimming data to %s" % datasize data = packet["data"][:datasize-2] crc = packet["data"][datasize-2:datasize] crc = struct.unpack(">H", crc)[0] crc = CheckCRC(data, crc) if not crc: print " WARNING: CRC does not match!" totalCRCErrors += 1 if crc or (EXPORTCORRUPT and not crc): firstorsinglepacket = packet["sequenceflag_int"] == 1 or packet["sequenceflag_int"] == 3 f = open(filename, "w" if firstorsinglepacket else "a") f.write(data[10:] if firstorsinglepacket else data) # Remove transport layer size f.close() if packet["sequenceflag_int"] == 2 or packet["sequenceflag_int"] == 3: packetmanager.manageFile(filename) totalSavedPackets += 1 else: print " Corrupted frame, skipping..." def CreatePacket(data): while True: if len(data) < 6: return -1, data version, type, shf, apid, sequenceflag, packetnumber, packetlength, data = ParseMSDU(data) pdata = data[:packetlength+2] if apid != 2047: pendingpackets[apid] = { "data": pdata, "version": version, "type": type, "apid": apid, "sequenceflag": SEQUENCE_FLAG_MAP[sequenceflag], "sequenceflag_int": sequenceflag, "packetnumber": packetnumber, "framesdropped": False, "size": packetlength } print "- Creating packet %s Size: %s - %s" % (apid, packetlength, SEQUENCE_FLAG_MAP[sequenceflag]) else: apid = -1 if not packetlength+2 == len(data) and packetlength+2 < len(data): # Multiple packets in buffer SavePacket(sys.argv[1], pendingpackets[apid]) del pendingpackets[apid] data = data[packetlength+2:] #print " Multiple packets in same buffer. Repeating." else: break return apid, "" if len(sys.argv) < 2: print "Usage: ./channeldecode.py CHANNELID" print "This will open channels/channel_CHANNELID.bin" exit() filename = "channels/channel_%s.bin" % sys.argv[1] f = open(filename, "r") fsize = os.path.getsize(filename) readbytes = 0 pendingpackets = {} lastFrameNumber = -1 totalFrameDrops = 0 totalCRCErrors = 0 totalSavedPackets = 0 lastAPID = -1 buff = "" while readbytes < fsize: if fsize - readbytes < FRAMESIZE: print " Some bytes at end of file was not enough for filling a frame. Remaining Bytes: %s - Frame Size: %s" % (fsize-readsize, FRAMESIZE) break # Read Data data = f.read(FRAMESIZE) versionNumber = (ord(data[0]) & 0xC0) >> 6 scid = (ord(data[0]) & 0x3F) << 2 | (ord(data[1]) & 0xC0) >> 6 vcid = (ord(data[1]) & 0x3F) counter = struct.unpack(">I", data[2:6])[0] counter &= 0xFFFFFF00 counter >>= 8 # Check for dropped Frames if not lastFrameNumber == -1 and not lastFrameNumber+1 == counter: print " Frames dropped: %s" % (counter-lastFrameNumber-1); totalFrameDrops += counter-lastFrameNumber-1; if not lastAPID == -1: # Fill pendingpackets[lastAPID]["data"] += "\x00" * 878 pendingpackets[lastAPID]["framesdropped"] = True #print "SC: %s ID: %s Frame Number: %s" % (scid, vcid, counter) # Demux M_PDU data = data[6:] # Strip channel header fhp = struct.unpack(">H", data[:2])[0] & 0x7FF data = data[2:] # Strip M_PDU Header #print " First Packet Header: %s" %fhp #data is now TP_PDU if not fhp == 2047: # Frame Contains a new Packet # Data was incomplete on last FHP and another packet starts here. if lastAPID == -1 and len(buff) > 0: #print " Data was incomplete from last FHP. Parsing packet now" if fhp > 0: buff += data[:fhp] lastAPID, data = CreatePacket(buff) if lastAPID == -1: buff = data else: buff = "" if not lastAPID == -1: # We are finishing another packet if fhp > 0: pendingpackets[lastAPID]["data"] += data[:fhp] SavePacket(sys.argv[1], pendingpackets[lastAPID]) del pendingpackets[lastAPID] lastAPID = -1 # Try to create a new packet buff += data[fhp:] lastAPID, data = CreatePacket(buff) if lastAPID == -1: buff = data else: buff = "" else: if len(buff) > 0 and lastAPID == -1: #print " Data was incomplete from last FHP. Parsing packet now" buff += data lastAPID, data = CreatePacket(buff) if lastAPID == -1: buff = data else: buff = "" elif len(buff) > 0: print " PROBLEM!" elif lastAPID == -1: buff += data lastAPID, data = CreatePacket(buff) if lastAPID == -1: buff = data else: buff = "" else: #print " Appending %s bytes to %s" % (lastAPID, len(data)) pendingpackets[lastAPID]["data"] += data lastFrameNumber = counter readbytes += FRAMESIZE # One packet can be still in pending packets for i in pendingpackets.keys(): SavePacket(sys.argv[1], pendingpackets[lastAPID]) print "\n\nReport:" print "\tTotal Frames Dropped: %s" %totalFrameDrops print "\tTotal Saved Packets: %s" %totalSavedPackets print "\tTotal Packet CRC Fails: %s" %totalCRCErrors f.close()

#!/usr/bin/env python import sys, struct, os from packetmanager import * if len(sys.argv) < 2: print "Usage: packetinfo.py filename.lrit" exit(1) filename = sys.argv[1] f = open(filename, "r") fsize = os.path.getsize(filename) readbytes = 0 t = 0 data = f.read() headers = getHeaderData(data) hsize = 0 ftype = 0 for i in headers: if i["type"] == 0: hsize = i["headerlength"] ftype = i["filetypecode"] printHeaders(headers) data = data[hsize:] if ftype == 2: print data f.close()

#!/usr/bin/env python import os, struct ''' File Type Codes: 0 - Image 2 - Text 130 - DCS Data ''' def manageFile(filename): f = open(filename, "r") try: k = readHeader(f) print k type, filetypecode, headerlength, datalength = k except: print " Header 0 is corrupted for file %s" %filename return newfilename = filename while f.tell() < headerlength: data = readHeader(f) if data[0] == 4: print " Filename is %s" % data[1] newfilename = data[1] break f.close() if filename != newfilename: print " Renaming %s to %s/%s" %(filename, os.path.dirname(filename), newfilename) os.rename(filename, "%s/%s" %(os.path.dirname(filename), newfilename)) #os.unlink(filename) else: print " Couldn't find name in %s" %filename def getHeaderData(data): headers = [] while len(data) > 0: type = ord(data[0]) size = struct.unpack(">H", data[1:3])[0] o = data[3:size] data = data[size:] td = parseHeader(type, o) headers.append(td) if td["type"] == 0: print "Header Size: %s" % td["headerlength"] data = data[:td["headerlength"]-size] return headers def parseHeader(type, data): if type == 0: filetypecode, headerlength, datalength = struct.unpack(">BIQ", data) return {"type":type, "filetypecode":filetypecode, "headerlength":headerlength, "datalength":datalength} elif type == 1: bitsperpixel, columns, lines, compression = struct.unpack(">BHHB", data) return {"type":type, "bitsperpixel":bitsperpixel, "columns":columns, "lines":lines, "compression":compression} elif type == 2: projname, cfac, lfac, coff, loff = struct.unpack(">32sIIII", data) return {"type":type, "projname":projname, "cfac":cfac, "lfac":lfac, "coff":coff, "loff":loff} elif type == 3: return {"type":type, "data":data} elif type == 4: return {"type":type, "filename":data} elif type == 5: days, ms = struct.unpack(">HI", data[1:]) return {"type":type, "days":days, "ms":ms} elif type == 6: return {"type":type, "data":data} elif type == 7: return {"type":type, "data":data} elif type == 128: imageid, sequence, startcol, startline, maxseg, maxcol, maxrow = struct.unpack(">7H", data) return {"type":type, "imageid":imageid, "sequence":sequence, "startcol":startcol, "startline":startline, "maxseg":maxseg, "maxcol":maxcol, "maxrow":maxrow} elif type == 129: signature, productId, productSubId, parameter, compression = struct.unpack(">4sHHHB", data) return {"type":type, "signature":signature, "productId":productId, "productSubId":productSubId, "parameter":parameter, "compression":compression} elif type == 130: return {"type":type, "data":data} elif type == 131: flags, pixel, line = struct.unpack(">HBB", data) return {"type":type, "flags":flags, "pixel":pixel, "line":line} elif type == 132: return {"type":type, "data": data} else: return {"type":type} def readHeader(f): global t type = ord(f.read(1)) size = f.read(2) size = struct.unpack(">H", size)[0] data = f.read(size-3) if type == 0: filetypecode, headerlength, datalength = struct.unpack(">BIQ", data) return type, filetypecode, headerlength, datalength elif type == 1: bitsperpixel, columns, lines, compression = struct.unpack(">BHHB", data) return type, bitsperpixel, columns, lines, compression elif type == 2: projname, cfac, lfac, coff, loff = struct.unpack(">32sIIII", data) return type, projname, cfac, lfac, coff, loff elif type == 3: return type, data elif type == 4: return type, data elif type == 5: days, ms = struct.unpack(">HI", data[1:]) return type, days, ms elif type == 6: return type, data elif type == 7: return type, data elif type == 128: imageid, sequence, startcol, startline, maxseg, maxcol, maxrow = struct.unpack(">7H", data) return type, imageid, sequence, startcol, startline, maxseg, maxcol, maxrow elif type == 129: signature, productId, productSubId, parameter, compression = struct.unpack(">4sHHHB", data) return type, signature, productId, productSubId, parameter, compression elif type == 130: return type, data elif type == 131: flags, pixel, line = struct.unpack(">HBB", data) return type, flags, pixel, line elif type == 132: return type, data else: return type def printHeaders(headers, showStructuredHeader=False, showImageDataRecord=False): for head in headers: type = head["type"] if type == 0: print "Header type: %s File Type Code: %s Header Length %s Data Field Length: %s" %(type, head["filetypecode"], head["headerlength"], head["datalength"]) elif type == 1: print "Image Structure Header: " print " Bits Per Pixel: %s" %head["bitsperpixel"] print " Columns: %s" %head["columns"] print " Lines: %s" %head["lines"] print " Compression: %s" %head["compression"] elif type == 2: print "Image Navigation Record" print " Projection Name: %s" %head["projname"] print " Column Scaling Factor: %s" %head["cfac"] print " Line Scaling Factor: %s" %head["lfac"] print " Column Offset: %s" %head["coff"] print " Line Offset: %s" %head["loff"] elif type == 3: print "Image Data Function Record" if showImageDataRecord: print " Data: %s" %head["data"] else: print " Data: {HIDDEN}" elif type == 4: print "Annotation Record" print " Filename: %s" %head["filename"] elif type == 5: print "Timestamp Record" print " Delta from 1 January 1958" print " Days: %s" %head["days"] print " Miliseconds: %s" %head["ms"] elif type == 6: print "Ancillary Text" print " Data: " t = head["data"].split(";") for i in t: print " %s" %i elif type == 7: print "Key Header" print " Data: %s" %head["data"] elif type == 128: print "Segment Identification Header" print " Image Id: %s" %head["imageid"] print " Sequence: %s" %head["sequence"] print " Start Column: %s" %head["startcol"] print " Start Line: %s" %head["startline"] print " Number of Segments: %s" %head["maxseg"] print " Width: %s" %head["maxcol"] print " Height: %s" %head["maxrow"] elif type == 129: print "NOAA Specific Header" print " Signature: %s" %head["signature"] print " Product ID: %s" %head["productId"] print " Product SubId: %s" %head["productSubId"] print " Parameter: %s" %head["parameter"] print " Compression: %s" %head["compression"] elif type == 130: print "Header Structured Record" if showImageDataRecord: t = head["data"].split("UI") print " Data: " for i in t: print " %s" %i else: print " Data: {HIDDEN}" elif type == 131: print "Rice Compression Record" print " Flags: %s" %head["flags"] print " Pixel: %s" %head["pixel"] print " Line: %s" %head["line"] elif type == 132: # Got in DCS Data print "DCS Data: " print " Data: %s" %head["data"] else: print "Type not mapped: %s" % type print ""