#!/usr/bin/env python # # A Python AIVDM/AIVDO decoder # # This file is Copyright (c) 2010 by the GPSD project # BSD terms apply: see the file COPYING in the distribution root for details. # # This code runs compatibly under Python 2 and 3.x for x >= 2. # Preserve this property! from __future__ import absolute_import, print_function, division from array import array import exceptions import re import sys try: BaseException.with_traceback def reraise_with_traceback(exc_type, exc_value, exc_traceback): raise exc_type(exc_value).with_traceback(exc_traceback) except AttributeError: def reraise_with_traceback(exc_type, exc_value, exc_traceback): raise exc_type, exc_value, exc_traceback # This decoder works by defining a declarative pseudolanguage in which # to describe the process of extracting packed bitfields from an AIS # message, a set of tables which contain instructions in the pseudolanguage, # and a small amount of code for interpreting it. # # Known bugs: # * Doesn't join parts A and B of Type 24 together yet. # * Only handles the broadcast case of type 22. The problem is that the # addressed field is located *after* the variant parts. Grrrr... # * Message type 26 is presently unsupported. It hasn't been observed # in the wild yet as of Jan 2010; not a lot of point in trying util # we have test data. We'd need new machinery to constrain how many # bits the data spec eats in order to recover the radio bits after it. # * No support for IMO236 and IMO289 special messages in types 6 and 8 yet. # # Decoding for 1-15, 18-21, and 24 have been tested against live data. # Decoding for 16-17, 22-23, and 25-27 have not. # Here are the pseudoinstructions in the pseudolanguage. class bitfield(object): "Object defining the interpretation of an AIS bitfield." # The only un-obvious detail is the use of the oob (out-of-band) # member. This isn't used in data extraction, but rather to cut # down on the number of custom formatting hooks. With this we # handle the case where the field should be reported as an integer # or "n/a". def __init__(self, name, width, dtype, oob, legend, validator=None, formatter=None, conditional=None): self.name = name # Fieldname, for internal use and JSON self.width = width # Bit width self.type = dtype # Data type: signed/unsigned/string/raw self.oob = oob # Out-of-band value to be shown as n/a self.legend = legend # Human-friendly description of field self.validator = validator # Validation checker self.formatter = formatter # Custom reporting hook. self.conditional = conditional # Evaluation guard for this field class spare(object): "Describes spare bits,, not to be interpreted." def __init__(self, width, conditional=None): self.width = width self.conditional = conditional # Evaluation guard for this field class dispatch(object): "Describes how to dispatch to a message type variant on a subfield value." def __init__(self, fieldname, subtypes, compute=lambda x: x, conditional=None): self.fieldname = fieldname # Value of view to dispatch on self.subtypes = subtypes # Possible subtypes to dispatch to self.compute = compute # Pass value through this pre-dispatch self.conditional = conditional # Evaluation guard for this field # Message-type-specific information begins here. There are four # different kinds of things in it: (1) string tables for expanding # enumerated-type codes, (2) hook functions, (3) instruction tables, # and (4) field group declarations. This is the part that could, in # theory, be generated from a portable higher-level specification in # XML; only the hook functions are actually language-specific, and # your XML definition could in theory embed several different ones for # code generation in Python, Java, Perl, etc. cnb_status_legends = ( "Under way using engine", "At anchor", "Not under command", "Restricted manoeuverability", "Constrained by her draught", "Moored", "Aground", "Engaged in fishing", "Under way sailing", "Reserved for HSC", "Reserved for WIG", "Reserved", "Reserved", "Reserved", "Reserved", "Not defined", ) def cnb_rot_format(n): if n == -128: return "n/a" elif n == -127: return "fastleft" elif n == 127: return "fastright" else: return str((n / 4.733) ** 2) def cnb_latlon_format(n): return str(n / 600000.0) def cnb_speed_format(n): if n == 1023: return "n/a" elif n == 1022: return "fast" else: return str(n / 10.0) def cnb_course_format(n): return str(n / 10.0) def cnb_second_format(n): if n == 60: return "n/a" elif n == 61: return "manual input" elif n == 62: return "dead reckoning" elif n == 63: return "inoperative" else: return str(n) # Common Navigation Block is the format for AIS types 1, 2, and 3 cnb = ( bitfield("status", 4, 'unsigned', 0, "Navigation Status", formatter=cnb_status_legends), bitfield("turn", 8, 'signed', -128, "Rate of Turn", formatter=cnb_rot_format), bitfield("speed", 10, 'unsigned', 1023, "Speed Over Ground", formatter=cnb_speed_format), bitfield("accuracy", 1, 'unsigned', None, "Position Accuracy"), bitfield("lon", 28, 'signed', 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, 'signed', 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("course", 12, 'unsigned', 0xe10, "Course Over Ground", formatter=cnb_course_format), bitfield("heading", 9, 'unsigned', 511, "True Heading"), bitfield("second", 6, 'unsigned', None, "Time Stamp", formatter=cnb_second_format), bitfield("maneuver", 2, 'unsigned', None, "Maneuver Indicator"), spare(3), bitfield("raim", 1, 'unsigned', None, "RAIM flag"), bitfield("radio", 19, 'unsigned', None, "Radio status"), ) epfd_type_legends = ( "Undefined", "GPS", "GLONASS", "Combined GPS/GLONASS", "Loran-C", "Chayka", "Integrated navigation system", "Surveyed", "Galileo", ) type4 = ( bitfield("year", 14, "unsigned", 0, "Year"), bitfield("month", 4, "unsigned", 0, "Month"), bitfield("day", 5, "unsigned", 0, "Day"), bitfield("hour", 5, "unsigned", 24, "Hour"), bitfield("minute", 6, "unsigned", 60, "Minute"), bitfield("second", 6, "unsigned", 60, "Second"), bitfield("accuracy", 1, "unsigned", None, "Fix quality"), bitfield("lon", 28, "signed", 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, "signed", 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("epfd", 4, "unsigned", None, "Type of EPFD", validator=lambda n: n >= 0 and n <= 8 or n == 15, formatter=epfd_type_legends), spare(10), bitfield("raim", 1, "unsigned", None, "RAIM flag "), bitfield("radio", 19, "unsigned", None, "SOTDMA state"), ) ship_type_legends = ( "Not available", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Wing in ground (WIG) - all ships of this type", "Wing in ground (WIG) - Hazardous category A", "Wing in ground (WIG) - Hazardous category B", "Wing in ground (WIG) - Hazardous category C", "Wing in ground (WIG) - Hazardous category D", "Wing in ground (WIG) - Reserved for future use", "Wing in ground (WIG) - Reserved for future use", "Wing in ground (WIG) - Reserved for future use", "Wing in ground (WIG) - Reserved for future use", "Wing in ground (WIG) - Reserved for future use", "Fishing", "Towing", "Towing: length exceeds 200m or breadth exceeds 25m", "Dredging or underwater ops", "Diving ops", "Military ops", "Sailing", "Pleasure Craft", "Reserved", "Reserved", "High speed craft (HSC) - all ships of this type", "High speed craft (HSC) - Hazardous category A", "High speed craft (HSC) - Hazardous category B", "High speed craft (HSC) - Hazardous category C", "High speed craft (HSC) - Hazardous category D", "High speed craft (HSC) - Reserved for future use", "High speed craft (HSC) - Reserved for future use", "High speed craft (HSC) - Reserved for future use", "High speed craft (HSC) - Reserved for future use", "High speed craft (HSC) - No additional information", "Pilot Vessel", "Search and Rescue vessel", "Tug", "Port Tender", "Anti-pollution equipment", "Law Enforcement", "Spare - Local Vessel", "Spare - Local Vessel", "Medical Transport", "Ship according to RR Resolution No. 18", "Passenger - all ships of this type", "Passenger - Hazardous category A", "Passenger - Hazardous category B", "Passenger - Hazardous category C", "Passenger - Hazardous category D", "Passenger - Reserved for future use", "Passenger - Reserved for future use", "Passenger - Reserved for future use", "Passenger - Reserved for future use", "Passenger - No additional information", "Cargo - all ships of this type", "Cargo - Hazardous category A", "Cargo - Hazardous category B", "Cargo - Hazardous category C", "Cargo - Hazardous category D", "Cargo - Reserved for future use", "Cargo - Reserved for future use", "Cargo - Reserved for future use", "Cargo - Reserved for future use", "Cargo - No additional information", "Tanker - all ships of this type", "Tanker - Hazardous category A", "Tanker - Hazardous category B", "Tanker - Hazardous category C", "Tanker - Hazardous category D", "Tanker - Reserved for future use", "Tanker - Reserved for future use", "Tanker - Reserved for future use", "Tanker - Reserved for future use", "Tanker - No additional information", "Other Type - all ships of this type", "Other Type - Hazardous category A", "Other Type - Hazardous category B", "Other Type - Hazardous category C", "Other Type - Hazardous category D", "Other Type - Reserved for future use", "Other Type - Reserved for future use", "Other Type - Reserved for future use", "Other Type - Reserved for future use", "Other Type - no additional information", ) type5 = ( bitfield("ais_version", 2, 'unsigned', None, "AIS Version"), bitfield("imo_id", 30, 'unsigned', 0, "IMO Identification Number"), bitfield("callsign", 42, 'string', None, "Call Sign"), bitfield("shipname", 120, 'string', None, "Vessel Name"), bitfield("shiptype", 8, 'unsigned', None, "Ship Type", # validator=lambda n: n >= 0 and n <= 99, formatter=ship_type_legends), bitfield("to_bow", 9, 'unsigned', 0, "Dimension to Bow"), bitfield("to_stern", 9, 'unsigned', 0, "Dimension to Stern"), bitfield("to_port", 6, 'unsigned', 0, "Dimension to Port"), bitfield("to_starbord", 6, 'unsigned', 0, "Dimension to Starboard"), bitfield("epfd", 4, 'unsigned', 0, "Position Fix Type", validator=lambda n: n >= 0 and n <= 8 or n == 15, formatter=epfd_type_legends), bitfield("month", 4, 'unsigned', 0, "ETA month"), bitfield("day", 5, 'unsigned', 0, "ETA day"), bitfield("hour", 5, 'unsigned', 24, "ETA hour"), bitfield("minute", 6, 'unsigned', 60, "ETA minute"), bitfield("draught", 8, 'unsigned', 0, "Draught", formatter=lambda n: n/10.0), bitfield("destination", 120, 'string', None, "Destination"), bitfield("dte", 1, 'unsigned', None, "DTE"), spare(1), ) type6_dac_or_fid_unknown = ( bitfield("data", 920, 'raw', None, "Data"), ) type6_dispatch = {} type6_dispatch[0] = type6_dac_or_fid_unknown # DAC 235 and 250 (UK, Rep. of Ireland) type6_dac235_dispatch = {} type6_dac235_dispatch[0] = type6_dac_or_fid_unknown type6_dac235_fid10 = ( bitfield("ana_int", 10, 'unsigned', None, "Supply voltage"), bitfield("ana_ext1", 10, 'unsigned', None, "Analogue (Ext#1)"), bitfield("ana_ext2", 10, 'unsigned', None, "Analogue (Ext#2)"), bitfield("racon", 2, 'unsigned', None, "RACON status"), bitfield("light", 2, 'unsigned', None, "Light status"), bitfield("health", 1, 'unsigned', None, "Health"), bitfield("stat_ext", 8, 'unsigned', None, "Status (ext)"), bitfield("off_pos", 1, 'unsigned', None, "Position status"), ) type6_dac235_dispatch[10] = type6_dac235_fid10 type6_dac235 = ( dispatch("fid", type6_dac235_dispatch, lambda m: m if m in type6_dac235_dispatch else 0), ) type6_dispatch[235] = type6_dac235 type6_dispatch[250] = type6_dac235 type6 = ( bitfield("seqno", 2, 'unsigned', None, "Sequence Number"), bitfield("dest_mmsi", 30, 'unsigned', None, "Destination MMSI"), bitfield("retransmit", 1, 'unsigned', None, "Retransmit flag"), spare(1), bitfield("dac", 10, 'unsigned', 0, "DAC"), bitfield("fid", 6, 'unsigned', 0, "Functional ID"), dispatch("dac", type6_dispatch, lambda m: m if m in type6_dispatch else 0), ) type7 = ( spare(2), bitfield("mmsi1", 30, 'unsigned', 0, "MMSI number 1"), spare(2), bitfield("mmsi2", 30, 'unsigned', 0, "MMSI number 2"), spare(2), bitfield("mmsi3", 30, 'unsigned', 0, "MMSI number 3"), spare(2), bitfield("mmsi1", 30, 'unsigned', 0, "MMSI number 4"), spare(2), ) # # Type 8 have subtypes identified by DAC (Designated Area Code) and # FID (Functional ID) # def type8_latlon_format(n): return str(n / 60000.0) type8_dac_or_fid_unknown = ( bitfield("data", 952, 'raw', None, "Data"), ) type8_dispatch = {} type8_dispatch[0] = type8_dac_or_fid_unknown # DAC 1 (international) type8_dac1_dispatch = {} type8_dac1_dispatch[0] = type8_dac_or_fid_unknown # DAC 1, FID 11: IMO236 Met/Hydro message def type8_dac1_fid11_airtemp_format(n): return str(n * 0.1 - 60) def type8_dac1_fid11_dewpoint_format(n): return str(n * 0.1 - 20) def type8_dac1_fid11_pressure_format(n): return str(n + 800) def type8_dac1_fid11_visibility_format(n): return str(n * 0.1) def type8_dac1_fid11_waterlevel_format(n): return str(n * 0.1 - 10) def type8_dac1_fid11_cspeed_format(n): return str(n * 0.1) def type8_dac1_fid11_waveheight_format(n): return str(n * 0.1) type8_dac1_fid11_seastate_legend = ( "Calm", "Light air", "Light breeze" "Gentle breeze", "Moderate breeze", "Fresh breeze", "Strong breeze", "High wind", "Gale", "Strong gale", "Storm", "Violent storm", "Hurricane force", "Reserved", "Reserved", "Reserved" ) def type8_dac1_fid11_watertemp_format(n): return str(n * 0.1 - 10) type8_dac1_fid11_preciptype_legend = ( "Reserved", "Rain", "Thunderstorm", "Freezing rain", "Mixed/ice", "Snow", "Reserved", "Reserved" ) def type8_dac1_fid11_salinity_format(n): return str(n * 0.1) type8_dac1_fid11_ice_legend = ( "Yes", "No" ) type8_dac1_fid11 = ( bitfield("lat", 24, "signed", 2**24-1, "Latitude", formatter=type8_latlon_format), bitfield("lon", 25, "signed", 2**25-1, "Longitude", formatter=type8_latlon_format), bitfield("day", 5, 'unsigned', 0, "ETA day"), bitfield("hour", 5, 'unsigned', 24, "ETA hour"), bitfield("minute", 6, 'unsigned', 60, "ETA minute"), bitfield("wspeed", 7, 'unsigned', 127, "Wind speed"), bitfield("wgust", 7, 'unsigned', 127, "Wind gust"), bitfield("wdir", 9, 'unsigned', 511, "Wind direction"), bitfield("wgustdir", 9, 'unsigned', 511, "Wind gust direction"), bitfield("airtemp", 11, 'unsigned', 2047, "Air temperature", formatter=type8_dac1_fid11_airtemp_format), bitfield("humidity", 7, 'unsigned', 127, "Relative humidity"), bitfield("dewpoint", 10, 'unsigned', 1023, "Dew point", formatter=type8_dac1_fid11_dewpoint_format), bitfield("pressure", 9, 'unsigned', 511, "Atmospheric pressure", formatter=type8_dac1_fid11_pressure_format), bitfield("pressuretend", 2, 'unsigned', 3, "Atmospheric pressure tendency"), bitfield("visibility", 8, 'unsigned', 255, "Horizontal visibility", formatter=type8_dac1_fid11_visibility_format), bitfield("waterlevel", 9, 'unsigned', 511, "Water level", formatter=type8_dac1_fid11_waterlevel_format), bitfield("leveltrend", 2, 'unsigned', 3, "Water level trend"), bitfield("cspeed", 8, 'unsigned', 255, "Surface current speed", formatter=type8_dac1_fid11_cspeed_format), bitfield("cdir", 9, 'unsigned', 511, "Surface current direction"), bitfield("cspeed2", 8, 'unsigned', 255, "Current speed #2", formatter=type8_dac1_fid11_cspeed_format), bitfield("cdir2", 9, 'unsigned', 511, "Current direction #2"), bitfield("cdepth2", 5, 'unsigned', 31, "Current measuring level #2"), bitfield("cspeed3", 8, 'unsigned', 255, "Current speed #3", formatter=type8_dac1_fid11_cspeed_format), bitfield("cdir3", 9, 'unsigned', 511, "Current direction #3"), bitfield("cdepth3", 5, 'unsigned', 31, "Current measuring level #3"), bitfield("waveheight", 8, 'unsigned', 255, "Significant wave height", formatter=type8_dac1_fid11_waveheight_format), bitfield("waveperiod", 6, 'unsigned', 63, "Significant wave period"), bitfield("wavedir", 9, 'unsigned', 511, "Significant wave direction"), bitfield("swellheight", 8, 'unsigned', 255, "Swell height", formatter=type8_dac1_fid11_waveheight_format), bitfield("swellperiod", 6, 'unsigned', 63, "Swell period"), bitfield("swelldir", 9, 'unsigned', 511, "Swell direction"), bitfield("seastate", 4, 'unsigned', 15, "Sea state", formatter=type8_dac1_fid11_seastate_legend), bitfield("watertemp", 10, 'unsigned', 1023, "Water temperature", formatter=type8_dac1_fid11_watertemp_format), bitfield("preciptype", 3, 'unsigned', 7, "Precipitation type", formatter=type8_dac1_fid11_preciptype_legend), bitfield("salinity", 9, 'unsigned', 511, "Salinity", formatter=type8_dac1_fid11_salinity_format), bitfield("ice", 2, 'unsigned', 3, "Ice?", formatter=type8_dac1_fid11_ice_legend), spare(6) ) type8_dac1_dispatch[11] = type8_dac1_fid11 type8_dac1 = ( dispatch("fid", type8_dac1_dispatch, lambda m: m if m in type8_dac1_dispatch else 0), ) type8_dispatch[1] = type8_dac1 type8 = ( spare(2), bitfield("dac", 10, 'unsigned', 0, "DAC"), bitfield("fid", 6, 'unsigned', 0, "Functional ID"), dispatch("dac", type8_dispatch, lambda m: m if m in type8_dispatch else 0), ) def type9_alt_format(n): if n == 4094: return ">=4094" else: return str(n) def type9_speed_format(n): if n == 1023: return "n/a" elif n == 1022: return "fast" else: return str(n) type9 = ( bitfield("alt", 12, 'unsigned', 4095, "Altitude", formatter=type9_alt_format), bitfield("speed", 10, 'unsigned', 1023, "SOG", formatter=type9_speed_format), bitfield("accuracy", 1, 'unsigned', None, "Position Accuracy"), bitfield("lon", 28, 'signed', 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, 'signed', 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("course", 12, 'unsigned', 0xe10, "Course Over Ground", formatter=cnb_course_format), bitfield("second", 6, 'unsigned', 60, "Time Stamp", formatter=cnb_second_format), bitfield("regional", 8, 'unsigned', None, "Regional reserved"), bitfield("dte", 1, 'unsigned', None, "DTE"), spare(3), bitfield("assigned", 1, 'unsigned', None, "Assigned"), bitfield("raim", 1, 'unsigned', None, "RAIM flag"), bitfield("radio", 19, 'unsigned', None, "Radio status"), ) type10 = ( spare(2), bitfield("dest_mmsi", 30, 'unsigned', None, "Destination MMSI"), spare(2), ) type12 = ( bitfield("seqno", 2, 'unsigned', None, "Sequence Number"), bitfield("dest_mmsi", 30, 'unsigned', None, "Destination MMSI"), bitfield("retransmit", 1, 'unsigned', None, "Retransmit flag"), spare(1), bitfield("text", 936, 'string', None, "Text"), ) type14 = ( spare(2), bitfield("text", 968, 'string', None, "Text"), ) type15 = ( spare(2), bitfield("mmsi1", 30, 'unsigned', 0, "First interrogated MMSI"), bitfield("type1_1", 6, 'unsigned', 0, "First message type"), bitfield("offset1_1", 12, 'unsigned', 0, "First slot offset"), spare(2), bitfield("type1_2", 6, 'unsigned', 0, "Second message type"), bitfield("offset1_2", 12, 'unsigned', 0, "Second slot offset"), spare(2), bitfield("mmsi2", 30, 'unsigned', 0, "Second interrogated MMSI"), bitfield("type2_1", 6, 'unsigned', 0, "Message type"), bitfield("offset2_1", 12, 'unsifned', 0, "Slot offset"), spare(2), ) type16 = ( spare(2), bitfield("mmsi1", 30, 'unsigned', 0, "Interrogated MMSI 1"), bitfield("offset1", 12, 'unsigned', 0, "First slot offset"), bitfield("increment1", 10, 'unsigned', 0, "First slot increment"), bitfield("mmsi2", 30, 'unsigned', 0, "Interrogated MMSI 2"), bitfield("offset2", 12, 'unsigned', 0, "Second slot offset"), bitfield("increment2", 10, 'unsigned', 0, "Second slot increment"), spare(2), ) def short_latlon_format(n): return str(n / 600.0) type17 = ( spare(2), bitfield("lon", 18, 'signed', 0x1a838, "Longitude", formatter=short_latlon_format), bitfield("lat", 17, 'signed', 0xd548, "Latitude", formatter=short_latlon_format), spare(5), bitfield("data", 736, 'raw', None, "DGNSS data"), ) type18 = ( bitfield("reserved", 8, 'unsigned', None, "Regional reserved"), bitfield("speed", 10, 'unsigned', 1023, "Speed Over Ground", formatter=cnb_speed_format), bitfield("accuracy", 1, 'unsigned', None, "Position Accuracy"), bitfield("lon", 28, 'signed', 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, 'signed', 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("course", 12, 'unsigned', 0xE10, "Course Over Ground", formatter=cnb_course_format), bitfield("heading", 9, 'unsigned', 511, "True Heading"), bitfield("second", 6, 'unsigned', None, "Time Stamp", formatter=cnb_second_format), bitfield("regional", 2, 'unsigned', None, "Regional reserved"), bitfield("cs", 1, 'unsigned', None, "CS Unit"), bitfield("display", 1, 'unsigned', None, "Display flag"), bitfield("dsc", 1, 'unsigned', None, "DSC flag"), bitfield("band", 1, 'unsigned', None, "Band flag"), bitfield("msg22", 1, 'unsigned', None, "Message 22 flag"), bitfield("assigned", 1, 'unsigned', None, "Assigned"), bitfield("raim", 1, 'unsigned', None, "RAIM flag"), bitfield("radio", 20, 'unsigned', None, "Radio status"), ) type19 = ( bitfield("reserved", 8, 'unsigned', None, "Regional reserved"), bitfield("speed", 10, 'unsigned', 1023, "Speed Over Ground", formatter=cnb_speed_format), bitfield("accuracy", 1, 'unsigned', None, "Position Accuracy"), bitfield("lon", 28, 'signed', 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, 'signed', 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("course", 12, 'unsigned', 0xE10, "Course Over Ground", formatter=cnb_course_format), bitfield("heading", 9, 'unsigned', 511, "True Heading"), bitfield("second", 6, 'unsigned', None, "Time Stamp", formatter=cnb_second_format), bitfield("regional", 4, 'unsigned', None, "Regional reserved"), bitfield("shipname", 120, 'string', None, "Vessel Name"), bitfield("shiptype", 8, 'unsigned', None, "Ship Type", # validator=lambda n: n >= 0 and n <= 99, formatter=ship_type_legends), bitfield("to_bow", 9, 'unsigned', 0, "Dimension to Bow"), bitfield("to_stern", 9, 'unsigned', 0, "Dimension to Stern"), bitfield("to_port", 6, 'unsigned', 0, "Dimension to Port"), bitfield("to_starbord", 6, 'unsigned', 0, "Dimension to Starboard"), bitfield("epfd", 4, 'unsigned', 0, "Position Fix Type", validator=lambda n: n >= 0 and n <= 8 or n == 15, formatter=epfd_type_legends), bitfield("assigned", 1, 'unsigned', None, "Assigned"), bitfield("raim", 1, 'unsigned', None, "RAIM flag"), bitfield("radio", 20, 'unsigned', None, "Radio status"), ) type20 = ( spare(2), bitfield("offset1", 12, 'unsigned', 0, "Offset number"), bitfield("number1", 4, 'unsigned', 0, "Reserved slots"), bitfield("timeout1", 3, 'unsigned', 0, "Time-out"), bitfield("increment1", 11, 'unsigned', 0, "Increment"), bitfield("offset2", 12, 'unsigned', 0, "Offset number 2"), bitfield("number2", 4, 'unsigned', 0, "Reserved slots"), bitfield("timeout2", 3, 'unsigned', 0, "Time-out"), bitfield("increment2", 11, 'unsigned', 0, "Increment"), bitfield("offset3", 12, 'unsigned', 0, "Offset number 3"), bitfield("number3", 4, 'unsigned', 0, "Reserved slots"), bitfield("timeout3", 3, 'unsigned', 0, "Time-out"), bitfield("increment3", 11, 'unsigned', 0, "Increment"), bitfield("offset4", 12, 'unsigned', 0, "Offset number 4"), bitfield("number4", 4, 'unsigned', 0, "Reserved slots"), bitfield("timeout4", 3, 'unsigned', 0, "Time-out"), bitfield("increment4", 11, 'unsigned', 0, "Increment"), ) aide_type_legends = ( "Unspecified", "Reference point", "RACON", "Fixed offshore structure", "Spare, Reserved for future use.", "Light, without sectors", "Light, with sectors", "Leading Light Front", "Leading Light Rear", "Beacon, Cardinal N", "Beacon, Cardinal E", "Beacon, Cardinal S", "Beacon, Cardinal W", "Beacon, Port hand", "Beacon, Starboard hand", "Beacon, Preferred Channel port hand", "Beacon, Preferred Channel starboard hand", "Beacon, Isolated danger", "Beacon, Safe water", "Beacon, Special mark", "Cardinal Mark N", "Cardinal Mark E", "Cardinal Mark S", "Cardinal Mark W", "Port hand Mark", "Starboard hand Mark", "Preferred Channel Port hand", "Preferred Channel Starboard hand", "Isolated danger", "Safe Water", "Special Mark", "Light Vessel / LANBY / Rigs", ) type21 = ( bitfield("aid_type", 5, 'unsigned', 0, "Aid type", formatter=aide_type_legends), bitfield("name", 120, 'string', None, "Name"), bitfield("accuracy", 1, 'unsigned', 0, "Position Accuracy"), bitfield("lon", 28, 'signed', 0x6791AC0, "Longitude", formatter=cnb_latlon_format), bitfield("lat", 27, 'signed', 0x3412140, "Latitude", formatter=cnb_latlon_format), bitfield("to_bow", 9, 'unsigned', 0, "Dimension to Bow"), bitfield("to_stern", 9, 'unsigned', 0, "Dimension to Stern"), bitfield("to_port", 6, 'unsigned', 0, "Dimension to Port"), bitfield("to_starboard", 6, 'unsigned', 0, "Dimension to Starboard"), bitfield("epfd", 4, 'unsigned', 0, "Position Fix Type", validator=lambda n: n >= 0 and n <= 8 or n == 15, formatter=epfd_type_legends), bitfield("second", 6, 'unsigned', 0, "UTC Second"), bitfield("off_position", 1, 'unsigned', 0, "Off-Position Indicator"), bitfield("regional", 8, 'unsigned', 0, "Regional reserved"), bitfield("raim", 1, 'unsigned', 0, "RAIM flag"), bitfield("virtual_aid", 1, 'unsigned', 0, "Virtual-aid flag"), bitfield("assigned", 1, 'unsigned', 0, "Assigned-mode flag"), spare(1), bitfield("name", 88, 'string', 0, "Name Extension"), ) type22 = ( spare(2), bitfield("channel_a", 12, 'unsigned', 0, "Channel A"), bitfield("channel_b", 12, 'unsigned', 0, "Channel B"), bitfield("txrx", 4, 'unsigned', 0, "Tx/Rx mode"), bitfield("power", 1, 'unsigned', 0, "Power"), bitfield("ne_lon", 18, 'signed', 0x1a838, "NE Longitude", formatter=short_latlon_format), bitfield("ne_lat", 17, 'signed', 0xd548, "NE Latitude", formatter=short_latlon_format), bitfield("sw_lon", 18, 'signed', 0x1a838, "SW Longitude", formatter=short_latlon_format), bitfield("sw_lat", 17, 'signed', 0xd548, "SW Latitude", formatter=short_latlon_format), bitfield("addressed", 1, 'unsigned', 0, "Addressed"), bitfield("band_a", 1, 'unsigned', 0, "Channel A Band"), bitfield("band_a", 1, 'unsigned', 0, "Channel A Band"), bitfield("zonesize", 3, 'unsigned', 0, "Zone size"), spare(23), ) station_type_legends = ( "All types of mobiles", "Reserved for future use", "All types of Class B mobile stations", "SAR airborne mobile station", "Aid to Navigation station", "Class B shipborne mobile station", "Regional use and inland waterways", "Regional use and inland waterways", "Regional use and inland waterways", "Regional use and inland waterways", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", "Reserved for future use", ) type23 = ( spare(2), bitfield("ne_lon", 18, 'signed', 0x1a838, "NE Longitude", formatter=short_latlon_format), bitfield("ne_lat", 17, 'signed', 0xd548, "NE Latitude", formatter=short_latlon_format), bitfield("sw_lon", 18, 'signed', 0x1a838, "SW Longitude", formatter=short_latlon_format), bitfield("sw_lat", 17, 'signed', 0xd548, "SW Latitude", formatter=short_latlon_format), bitfield("stationtype", 4, 'unsigned', 0, "Station Type", validator=lambda n: n >= 0 and n <= 31, formatter=station_type_legends), bitfield("shiptype", 8, 'unsigned', 0, "Ship Type", # validator=lambda n: n >= 0 and n <= 99, formatter=ship_type_legends), spare(22), bitfield("txrx", 2, 'unsigned', 0, "Tx/Rx mode"), bitfield("interval", 4, 'unsigned', 0, "Reporting interval"), bitfield("txrx", 4, 'unsigned', 0, "Quiet time"), ) type24a = ( bitfield("shipname", 120, 'string', None, "Vessel Name"), spare(8), ) type24b1 = ( bitfield("callsign", 42, 'string', None, "Call Sign"), bitfield("to_bow", 9, 'unsigned', 0, "Dimension to Bow"), bitfield("to_stern", 9, 'unsigned', 0, "Dimension to Stern"), bitfield("to_port", 6, 'unsigned', 0, "Dimension to Port"), bitfield("to_starbord", 6, 'unsigned', 0, "Dimension to Starboard"), spare(8), ) type24b2 = ( bitfield('mothership_mmsi', 30, 'unsigned', 0, "Mothership MMSI"), spare(8), ) type24b = ( bitfield("shiptype", 8, 'unsigned', None, "Ship Type", validator=lambda n: n >= 0 and n <= 99, formatter=ship_type_legends), bitfield("vendorid", 42, 'string', None, "Vendor ID"), dispatch("mmsi", {0: type24b1, 1: type24b2}, lambda m: 1 if repr(m)[:2] == '98' else 0), ) type24 = ( bitfield('partno', 2, 'unsigned', None, "Part Number"), dispatch('partno', {0: type24a, 1: type24b}), ) type25 = ( bitfield("addressed", 1, 'unsigned', None, "Addressing flag"), bitfield("structured", 1, 'unsigned', None, "Dimension to Bow"), bitfield("dest_mmsi", 30, 'unsigned', 0, "Destinstion MMSI", conditional=lambda i, v: v["addressed"]), bitfield("app_id", 16, 'unsigned', 0, "Application ID", conditional=lambda i, v: v["structured"]), bitfield("data", 0, 'raw', None, "Data"), ) # No type 26 handling yet, type27 = ( bitfield("accuracy", 1, 'unsigned', None, "Position Accuracy"), bitfield("raim", 1, 'unsigned', None, "RAIM flag"), bitfield("status", 4, 'unsigned', 0, "Navigation Status", formatter=cnb_status_legends), bitfield("lon", 18, 'signed', 0x1a838, "Longitude", formatter=short_latlon_format), bitfield("lat", 17, 'signed', 0xd548, "Latitude", formatter=short_latlon_format), bitfield("speed", 6, 'unsigned', 63, "Speed Over Ground", formatter=cnb_speed_format), bitfield("course", 9, 'unsigned', 511, "Course Over Ground"), bitfield("GNSS", 1, 'unsigned', None, "GNSS flag"), spare(1), ) aivdm_decode = ( bitfield('msgtype', 6, 'unsigned', 0, "Message Type", validator=lambda n: n > 0 and n <= 27), bitfield('repeat', 2, 'unsigned', None, "Repeat Indicator"), bitfield('mmsi', 30, 'unsigned', 0, "MMSI"), # This is the master dispatch on AIS message type dispatch('msgtype', {0: None, 1: cnb, 2: cnb, 3: cnb, 4: type4, 5: type5, 6: type6, 7: type7, 8: type8, 9: type9, 10: type10, 11: type4, 12: type12, 13: type7, 14: type14, 15: type15, 16: type16, 17: type17, 18: type18, 19: type19, 20: type20, 21: type21, 22: type22, 23: type23, 24: type24, 25: type25, 26: None, 27: type27}), ) # Length ranges. We use this for integrity checking. # When a range is a tuple, it's (minimum, maximum). lengths = { 1: 168, 2: 168, 3: 168, 4: 168, 5: 424, 6: (88, 1008), 7: (72, 168), 8: (56, 1008), 9: 168, 10: 72, 11: 168, 12: (72, 1008), 13: (72, 168), 14: (40, 1008), 15: (88, 168), 16: (96, 144), 17: (80, 816), 18: 168, 19: 312, 20: (72, 160), 21: (272, 368), 22: 168, 23: 160, 24: (160, 168), 25: 168, 26: (60, 1004), 27: 96, } field_groups = ( # This one occurs in message type 4 (3, ["year", "month", "day", "hour", "minute", "second"], "time", "Timestamp", lambda y, m, d, h, n, s: "%02d-%02d-%02dT%02d:%02d:%02dZ" % (y, m, d, h, n, s)), # This one is in message 5 (13, ["month", "day", "hour", "minute", "second"], "eta", "Estimated Time of Arrival", lambda m, d, h, n, s: "%02d-%02dT%02d:%02d:%02dZ" % (m, d, h, n, s)), ) # Message-type-specific information ends here. # # Next, the execution machinery for the pseudolanguage. There isn't much of # this: the whole point of the design is to embody most of the information # about the AIS format in the pseudoinstruction tables. BITS_PER_BYTE = 8 class BitVector(object): "Fast bit-vector class based on Python built-in array type." def __init__(self, data=None, length=None): self.bits = array('B') self.bitlen = 0 if data is not None: self.bits.extend(data) if length is None: self.bitlen = len(data) * 8 else: self.bitlen = length def extend_to(self, length): "Extend vector to given bitlength." if length > self.bitlen: self.bits.extend([0] * ((length - self.bitlen + 7) // 8)) self.bitlen = length def from_sixbit(self, data, pad=0): "Initialize bit vector from AIVDM-style six-bit armoring." self.bits.extend([0] * len(data)) for ch in data: ch = ord(ch) - 48 if ch > 40: ch -= 8 for i in (5, 4, 3, 2, 1, 0): if (ch >> i) & 0x01: self.bits[self.bitlen // 8] |= (1 << (7 - self.bitlen % 8)) self.bitlen += 1 self.bitlen -= pad def ubits(self, start, width): "Extract a (zero-origin) bitfield from the buffer as an unsigned int." fld = 0 for i in range(start // BITS_PER_BYTE, (start + width + BITS_PER_BYTE - 1) // BITS_PER_BYTE): fld <<= BITS_PER_BYTE fld |= self.bits[i] end = (start + width) % BITS_PER_BYTE if end != 0: fld >>= (BITS_PER_BYTE - end) fld &= ~(-1 << width) return fld def sbits(self, start, width): "Extract a (zero-origin) bitfield from the buffer as a signed int." fld = self.ubits(start, width) if fld & (1 << (width-1)): fld = -(2 ** width - fld) return fld def __len__(self): return self.bitlen def __repr__(self): "Used for dumping binary data." return (str(self.bitlen) + ":" + "".join(["%02x" % d for d in self.bits[:(self.bitlen + 7) // 8]])) class AISUnpackingException(exceptions.Exception): def __init__(self, lc, fieldname, value): self.lc = lc self.fieldname = fieldname self.value = value def __repr__(self): return ("%d: validation on fieldname %s failed (value %s)" % (self.lc, self.fieldname, self.value)) def aivdm_unpack(lc, data, offset, values, instructions): "Unpack fields from data according to instructions." cooked = [] for inst in instructions: if offset >= len(data): break elif (inst.conditional is not None and not inst.conditional(inst, values)): continue elif isinstance(inst, spare): offset += inst.width elif isinstance(inst, dispatch): i = inst.compute(values[inst.fieldname]) # This is the recursion that lets us handle variant types cooked += aivdm_unpack(lc, data, offset, values, inst.subtypes[i]) elif isinstance(inst, bitfield): if inst.type == 'unsigned': value = data.ubits(offset, inst.width) elif inst.type == 'signed': value = data.sbits(offset, inst.width) elif inst.type == 'string': value = '' # The try/catch error here is in case we run off the end # of a variable-length string field, as in messages 12 and 14 try: for i in range(inst.width // 6): newchar = ("@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\\]^- !\"" "#$%&'()*+,-./0123456789:;<=>?" [data.ubits(offset + 6 * i, 6)]) if newchar == '@': break else: value += newchar except IndexError: pass value = value.replace("@", " ").rstrip() elif inst.type == 'raw': # Note: Doesn't rely on the length. value = BitVector(data.bits[offset // 8:], len(data)-offset) values[inst.name] = value if inst.validator and not inst.validator(value): raise AISUnpackingException(lc, inst.name, value) offset += inst.width # An important thing about the unpacked representation this # generates is that it carries forward the meta-information from # the field type definition. This stuff is then available for # use by report-generating code. cooked.append([inst, value]) return cooked def packet_scanner(source): "Get a span of AIVDM packets with contiguous fragment numbers." payloads = {'A': '', 'B': ''} raw = '' well_formed = False lc = 0 while True: lc += 1 line = source.readline() if not line: return raw += line line = line.strip() # Strip off USCG metadata line = re.sub(r"(?<=\*[0-9A-F][0-9A-F]),.*", "", line) # Compute CRC-16 checksum packet = line[1:-3] # Strip leading !, trailing * and CRC csum = 0 for c in packet: csum ^= ord(c) csum = "%02X" % csum # Ignore comments if not line.startswith("!"): continue # Assemble fragments from single- and multi-line payloads fields = line.split(",") try: expect = fields[1] fragment = fields[2] channel = fields[4] if fragment == '1': payloads[channel] = '' well_formed = True payloads[channel] += fields[5] try: # This works because a mangled pad literal means # a malformed packet that will be caught by the CRC check. pad = int(fields[6].split('*')[0]) except ValueError: pad = 0 crc = fields[6].split('*')[1].strip() except IndexError: if skiperr: sys.stderr.write("%d: malformed line: %s\n" % (lc, line.strip())) well_formed = False else: raise AISUnpackingException(lc, "checksum", crc) if csum != crc: if skiperr: sys.stderr.write("%d: bad checksum %s, expecting %s: %s\n" % (lc, repr(crc), csum, line.strip())) well_formed = False else: raise AISUnpackingException(lc, "checksum", crc) if fragment < expect or not well_formed: continue # Render assembled payload to packed bytes bits = BitVector() bits.from_sixbit(payloads[channel], pad) yield (lc, raw, bits) raw = '' def postprocess(cooked): "Postprocess cooked fields from a message." # Handle type 21 name extension if cooked[0][1] == 21 and len(cooked) > 19: cooked[4][1] += cooked[19][1] cooked.pop(-1) return cooked def parse_ais_messages(source, scaled=False, skiperr=False, verbose=0): "Generator code - read forever from source stream, parsing AIS messages." values = {} for (lc, raw, bits) in packet_scanner(source): values['length'] = bits.bitlen # Without the following magic, we'd have a subtle problem near # certain variable-length messages: DSV reports would # sometimes have fewer fields than expected, because the # unpacker would never generate cooked tuples for the omitted # part of the message. Presently a known issue for types 15 # and 16 only. (Would never affect variable-length messages in # which the last field type is 'string' or 'raw'). bits.extend_to(168) # Magic recursive unpacking operation try: cooked = aivdm_unpack(lc, bits, 0, values, aivdm_decode) # We now have a list of tuples containing unpacked fields # Collect some field groups into ISO8601 format for (offset, template, label, legend, formatter) in field_groups: segment = cooked[offset:offset+len(template)] if [x[0] for x in segment] == template: group = formatter(*[x[1] for x in segment]) group = (label, group, 'string', legend, None) cooked = (cooked[:offset] + [group] + cooked[offset+len(template):]) # Apply the postprocessor stage cooked = postprocess(cooked) # Now apply custom formatting hooks. if scaled: for (i, (inst, value)) in enumerate(cooked): if value == inst.oob: cooked[i][1] = "n/a" elif inst.formatter: if isinstance(inst.formatter, tuple): # Assumes 0 is the legend for the "undefined" value if value >= len(inst.formatter): value = 0 cooked[i][1] = inst.formatter[value] elif isinstance(formatter, function): cooked[i][1] = inst.formatter(value) expected = lengths.get(values['msgtype'], None) # Check length; has to be done after so we have the type field bogon = False if expected is not None: if isinstance(expected, int): expected_range = (expected, expected) else: expected_range = expected actual = values['length'] if not (actual >= expected_range[0] and actual <= expected_range[1]): bogon = True if skiperr: sys.stderr.write( "%d: type %d expected %s bits but saw %s: %s\n" % (lc, values['msgtype'], expected, actual, raw.strip().split())) else: raise AISUnpackingException(lc, "length", actual) # We're done, hand back a decoding values = {} yield (raw, cooked, bogon) raw = '' except KeyboardInterrupt: raise KeyboardInterrupt except GeneratorExit: raise GeneratorExit except AISUnpackingException as e: if skiperr: sys.stderr.write("%s: %s\n" % (repr(e), raw.strip().split())) continue else: raise except Exception: (exc_type, exc_value, exc_traceback) = sys.exc_info() sys.stderr.write("%d: Unknown exception: %s\n" % (lc, raw.strip().split())) if skiperr: continue else: reraise_with_traceback(exc_type, exc_value, exc_traceback) # The rest is just sequencing and report generation. if __name__ == "__main__": import getopt import sys try: (options, arguments) = getopt.getopt(sys.argv[1:], "cdhjmqst:vx") except getopt.GetoptError as msg: print("ais.py: " + str(msg)) raise SystemExit(1) dsv = False dump = False histogram = False json = False malformed = False quiet = False scaled = False types = [] frequencies = {} verbose = 0 skiperr = True for (switch, val) in options: if switch == '-c': # Report in DSV format rather than JSON dsv = True elif switch == '-d': # Dump in a more human-readable format dump = True elif switch == '-h': # Make a histogram of type frequencies histogram = True elif switch == '-j': # Dump JSON json = True elif switch == '-m': # Dump malformed AIVDM/AIVDO packets raw malformed = True elif switch == '-q': # Suppress output quiet = True elif switch == '-s': # Report AIS in scaled form scaled = True elif switch == '-t': # Filter for a comma-separated list of types types = list(map(int, val.split(","))) elif switch == '-v': # Dump raw packet before JSON or DSV. verbose += 1 elif switch == '-x': # Skip decoding errors skiperr = False if not dsv and not histogram and not json and not malformed and not quiet: dump = True try: for ((raw, parsed, bogon) in parse_ais_messages(sys.stdin, scaled, skiperr, verbose)): msgtype = parsed[0][1] if types and msgtype not in types: continue if verbose >= 1 or (bogon and malformed): sys.stdout.write(raw) if not bogon: if json: def quotify(x): if isinstance(x, str): return '"' + str(x) + '"' else: return str(x) print("{" + ",".join(['"' + x[0].name + '":' + quotify(x[1]) for x in parsed]) + "}") elif dsv: print("|".join([str(x[1]) for x in parsed])) elif histogram: key = "%02d" % msgtype frequencies[key] = frequencies.get(key, 0) + 1 if msgtype == 6 or msgtype == 8: dac = 0 fid = 0 if msgtype == 8: dac = parsed[3][1] fid = parsed[4][1] elif msgtype == 6: dac = parsed[6][1] fid = parsed[7][1] key = "%02d_%04d_%02d" % (msgtype, dac, fid) frequencies[key] = frequencies.get(key, 0) + 1 elif dump: for (inst, value) in parsed: print("%-25s: %s" % (inst.legend, value)) print("%%") sys.stdout.flush() if histogram: keys = list(frequencies.keys()) keys.sort() for msgtype in keys: print("%-33s\t%d" % (msgtype, frequencies[msgtype])) except KeyboardInterrupt: pass # End