Filters

CHANGELOG.txt

@@ -1,6 +1,14 @@
 ====================
  pyais CHANGELOG
 ====================
+-------------------------------------------------------------------------------
+ Version 2.5.9 26 Nov 2023
+-------------------------------------------------------------------------------
+* Initial release of the AIS Filters system.
+* Basic filtering classes: `NoneFilter`, `MessageTypeFilter`, and `DistanceFilter`.
+* `FilterChain` class to combine multiple filters into a sequence.
+* Utility functions `haversine` and `is_in_grid` for distance and grid calculations.
+
 -------------------------------------------------------------------------------
  Version 2.5.9 26 Nov 2023
 -------------------------------------------------------------------------------

README.md

@@ -327,6 +327,93 @@ Such details include information used by the slot allocation algorithm (either S
 
 Refer to [readthedocs](https://pyais.readthedocs.io/en/latest/messages.html#communication-state) for more information.
 
+
+# AIS Filters
+
+The filtering system is built around a series of filter classes, each designed to filter messages based on specific criteria. Filters are chained together using the `FilterChain` class, which allows combining multiple filters into a single, sequential filtering process. The system is flexible, allowing for the easy addition or removal of filters from the chain.
+
+### How It Works
+
+1. **AIS Stream**: Messages are provided as a stream to the filters.
+2. **Filter Application**: Each filter in the chain applies its criteria to the stream, passing the messages that meet the criteria to the next filter.
+3. **Filter Chain**: The `FilterChain` class orchestrates the passing of messages through each filter, from the first to the last.
+
+## Filters
+
+### 1. AttributeFilter
+
+- **Description**: Filters messages based on a user-defined function.
+- **Usage**: Initialize with a function that takes an AIS message and returns `True` if the message should be kept.
+
+### 2. NoneFilter
+
+- **Description**: Filters out messages where specified attributes are `None`.
+- **Usage**: Initialize with the names of attributes that should not be `None` in the messages.
+
+### 3. MessageTypeFilter
+
+- **Description**: Filters messages based on their type.
+- **Usage**: Initialize with message types to include.
+
+### 4. DistanceFilter
+
+- **Description**: Filters messages based on distance from a reference point.
+- **Usage**: Initialize with a reference point (latitude and longitude) and a distance threshold in kilometers.
+
+### 5. GridFilter
+
+- **Description**: Filters messages based on whether they fall within a specified geographical grid.
+- **Usage**: Initialize with the boundaries of the grid (minimum and maximum latitude and longitude).
+
+## Utility Functions
+
+### 1. Haversine
+
+- **Description**: Calculates the great circle distance between two points on the Earth.
+- **Parameters**: Takes two tuples representing the latitude and longitude of each point.
+- **Returns**: Distance between the points in kilometers.
+
+### 2. Is In Grid
+
+- **Description**: Checks if a point is within a defined geographical grid.
+- **Parameters**: Latitude and longitude of the point and the boundaries of the grid.
+- **Returns**: `True` if the point is within the grid, `False` otherwise.
+
+## FilterChain
+
+- **Description**: Chains multiple filters together into a single filtering process.
+- **Usage**: Initialize with a list of filters to be applied in order. The chain can be used to filter a stream of AIS messages.
+
+## Example Usage
+
+```python
+from pyais import decode
+# ... (importing necessary classes)
+
+# Define and initialize filters
+attribute_filter = AttributeFilter(lambda x: not hasattr(x, 'turn') or x.turn == -128.0)
+none_filter = NoneFilter('lon', 'lat', 'mmsi2')
+message_type_filter = MessageTypeFilter(1, 2, 3)
+distance_filter = DistanceFilter((51.900, 5.320), distance_km=1000)
+grid_filter = GridFilter(lat_min=50, lon_min=0, lat_max=52, lon_max=5)
+
+# Create a filter chain
+chain = FilterChain([
+    attribute_filter,
+    none_filter,
+    message_type_filter,
+    distance_filter,
+    grid_filter,
+])
+
+# Decode AIS data and filter
+data = [decode(b"!AIVDM..."), ...]
+filtered_data = list(chain.filter(data))
+
+for msg in filtered_data:
+    print(msg.lat, msg.lon)
+```
+
 # AIS tracker
 
 **pyais** comes with the the ability to collect and maintain the state of individual vessels over time.

docs/filters.rst

@@ -0,0 +1,94 @@
+=========================
+AIS Filters Documentation
+=========================
+
+AIS filters provide a modular and flexible way to filter AIS (Automatic Identification System) messages based on various criteria such as attributes, message types, geographical boundaries, and more. Filters can be chained together to create complex filtering logic.
+
+Overview
+********
+
+The filtering system is built around a series of filter classes, each designed to filter messages based on specific criteria. Filters are chained together using the ``FilterChain`` class, which allows combining multiple filters into a single, sequential filtering process. The system is flexible, allowing for the easy addition or removal of filters from the chain.
+
+**How It Works**
+
+#. **AIS Stream**: Messages are provided as a stream to the filters.
+#. **Filter Application**: Each filter in the chain applies its criteria to the stream, passing the messages that meet the criteria to the next filter.
+#. **Filter Chain**: The ``FilterChain`` class orchestrates the passing of messages through each filter, from the first to the last.
+
+Filters
+*******
+
+1. AttributeFilter
+   - **Description**: Filters messages based on a user-defined function.
+   - **Usage**: Initialize with a function that takes an AIS message and returns ``True`` if the message should be kept.
+
+2. NoneFilter
+   - **Description**: Filters out messages where specified attributes are ``None``.
+   - **Usage**: Initialize with the names of attributes that should not be ``None`` in the messages.
+
+3. MessageTypeFilter
+   - **Description**: Filters messages based on their type.
+   - **Usage**: Initialize with message types to include.
+
+4. DistanceFilter
+   - **Description**: Filters messages based on distance from a reference point.
+   - **Usage**: Initialize with a reference point (latitude and longitude) and a distance threshold in kilometers.
+
+5. GridFilter
+   - **Description**: Filters messages based on whether they fall within a specified geographical grid.
+   - **Usage**: Initialize with the boundaries of the grid (minimum and maximum latitude and longitude).
+
+Utility Functions
+*****************
+
+1. Haversine
+   - **Description**: Calculates the great circle distance between two points on the Earth.
+   - **Parameters**: Takes two tuples representing the latitude and longitude of each point.
+   - **Returns**: Distance between the points in kilometers.
+
+2. Is In Grid
+   - **Description**: Checks if a point is within a defined geographical grid.
+   - **Parameters**: Latitude and longitude of the point and the boundaries of the grid.
+   - **Returns**: ``True`` if the point is within the grid, ``False`` otherwise.
+
+FilterChain
+***********
+
+- **Description**: Chains multiple filters together into a single filtering process.
+- **Usage**: Initialize with a list of filters to be applied in order. The chain can be used to filter a stream of AIS messages.
+
+Example Usage
+*************
+
+.. code-block:: python
+
+    from pyais import decode
+    # ... (importing necessary classes)
+
+    # Define and initialize filters
+    attribute_filter = AttributeFilter(lambda x: not hasattr(x, 'turn') or x.turn == -128.0)
+    none_filter = NoneFilter('lon', 'lat', 'mmsi2')
+    message_type_filter = MessageTypeFilter(1, 2, 3)
+    distance_filter = DistanceFilter((51.900, 5.320), distance_km=1000)
+    grid_filter = GridFilter(lat_min=50, lon_min=0, lat_max=52, lon_max=5)
+
+    # Create a filter chain
+    chain = FilterChain([
+        attribute_filter,
+        none_filter,
+        message_type_filter,
+        distance_filter,
+        grid_filter,
+    ])
+
+    # Decode AIS data and filter
+    data = [decode(b"!AIVDM..."), ...]
+    filtered_data = list(chain.filter(data))
+
+    for msg in filtered_data:
+        print(msg.lat, msg.lon)
+
+Conclusion
+**********
+
+The AIS filter system provides a powerful and flexible way to filter AIS messages based on a wide range of criteria. By chaining together filters, you can create complex filtering logic tailored to your specific needs.

docs/index.rst

@@ -13,3 +13,4 @@ Welcome to Pyais's documentation!
    install
    examples
    messages
+   filters

examples/filters.py

@@ -0,0 +1,46 @@
+from pyais import decode
+from pyais.filter import (
+    AttributeFilter,
+    DistanceFilter,
+    FilterChain,
+    GridFilter,
+    MessageTypeFilter,
+    NoneFilter
+)
+
+# Define the filter chain with various criteria
+chain = FilterChain([
+    # Filter out messages based on the 'turn' attribute or lack thereof
+    AttributeFilter(lambda x: not hasattr(x, 'turn') or x.turn == -128.0),
+
+    # Ensure 'lon', 'lat', and 'mmsi2' attributes are not None
+    NoneFilter('lon', 'lat', 'mmsi2'),
+
+    # Include only messages of type 1, 2, or 3
+    MessageTypeFilter(1, 2, 3),
+
+    # Limit messages to within 1000 km of a specific point
+    DistanceFilter((51.900, 5.320), distance_km=1000),
+
+    # Restrict messages to a specific geographic grid
+    GridFilter(lat_min=50, lon_min=0, lat_max=52, lon_max=5),
+])
+
+# Example AIS data to filter
+data = [
+    decode(b"!AIVDM,1,1,,B,15NG6V0P01G?cFhE`R2IU?wn28R>,0*05"),
+    decode(b"!AIVDM,1,1,,A,13HOI:0P0000VOHLCnHQKwvL05Ip,0*23"),
+    decode(b"!AIVDM,1,1,,B,100h00PP0@PHFV`Mg5gTH?vNPUIp,0*3B"),
+    decode(b"!AIVDM,1,1,,A,133sVfPP00PD>hRMDH@jNOvN20S8,0*7F"),
+    decode(b"!AIVDM,1,1,,B,13eaJF0P00Qd388Eew6aagvH85Ip,0*45"),
+    decode(b"!AIVDM,1,1,,A,14eGrSPP00ncMJTO5C6aBwvP2D0?,0*7A"),
+    decode(b"!AIVDM,1,1,,A,15MrVH0000KH<:V:NtBLoqFP2H9:,0*2F"),
+    decode(b"!AIVDM,1,1,,A,702R5`hwCjq8,0*6B"),
+]
+
+# Filter the data using the defined chain
+filtered_data = list(chain.filter(data))
+
+# Print the latitude and longitude of each message that passed the filters
+for msg in filtered_data:
+    print(msg.lat, msg.lon)

pyais/__init__.py

@@ -5,7 +5,7 @@
 from pyais.tracker import AISTracker, AISTrack
 
 __license__ = 'MIT'
-__version__ = '2.5.9'
+__version__ = '2.6.0'
 __author__ = 'Leon Morten Richter'
 
 __all__ = (