Track interpolation with AISdb involves generating estimated positions of vessels at specific intervals when actual AIS data points are unavailable. This process is important for filling in gaps in the vessel's trajectory, which can occur due to signal loss, data filtering, or other disruptions.
In this tutorial, we introduce different types of track interpolation implemented in AISdb with usage examples.
Example data preparation
First, we defined functions to transform and visualize the track data (a generator object), with options to view the data points or the tracks:
import plotly.graph_objects as goimport plotly.express as pximport pandas as pddeftrack2dataframe(tracks): data =[]for track in tracks: times = track['time'] mmsi = track['mmsi'] lons = track['lon'] lats = track['lat']# Iterate over the dynamic arrays and create a row for each time pointfor i inrange(len(times)): data.append({'mmsi': mmsi,'time': times[i],'longitude': lons[i],'latitude': lats[i],})return pd.DataFrame(data)defplotly_visualize(data,visual_type='lines'):if(visual_type=='scatter'):# Create a scatter plot for the vessel data points using scatter_geo fig = px.scatter_geo( data,lat="latitude",lon="longitude",color="mmsi",# Color by vessel identifierhover_name="mmsi",hover_data={"time":True},title="Vessel Data Points")else:# Create a line plot for the vessel trajectory using line_geo fig = px.line_geo( data,lat="latitude",lon="longitude",color="mmsi",# Color by vessel identifierhover_name="mmsi",hover_data={"time":True},title="Vessel Trajectory")# Set the map style and projection fig.update_geos(projection_type="azimuthal equal area",# Change this to 'natural earth', 'azimuthal equal area', etc.showland=True,landcolor="rgb(243, 243, 243)",countrycolor="rgb(204, 204, 204)",)# Set the layout to focus on a specific area or zoom level fig.update_layout(geo=dict(projection_type="azimuthal equal area",),width=1200,# Increase the width of the plotheight=800,# Increase the height of the plot) fig.show()
We will use an actual track retrieved from the database for the examples in this tutorial and interpolate additional data points based on this track. The visualization will show the original track data points:
Original data points of the vessel track queried from database
Linear interpolation
Linear interpolation estimates the vessel's position by drawing a straight line between two known points and calculating the positions at intermediate times. It is simple, fast, and straightforward but may not accurately represent complex movements.
With equal time window intervals
This method estimates the position of a vessel at regular time intervals (e.g., every 10 minutes). To perform linear interpolation with an equal time window on the track defined above:
Linear interpolation on the vessel track with equal time intervals
With equal distance intervals
This method estimates the position of a vessel at regular spatial intervals (e.g., every 1 km along its path). To perform linear interpolation with equal distance intervals on the pseudo track defined above:
Linear interpolation with equal distance intervals
Geodesic Track Interpolation
This method estimates the positions of a vessel along a curved path using the principles of geometry, particularly involving great-circle routes.
Linear interpolation of the vessel track along the geodesic curve
Cubic Spline Interpolation
Given a set of data points, cubic spline interpolation fits a smooth curve through these points. The curve is represented as a series of cubic polynomials between each pair of data points. Each polynomial ensures a smooth curve at the data points (i.e., the first and second derivatives are continuous).
Cubic spline interpolation with equal time intervals
Custom Track Interpolation
In addition to the standard interpolation methods provided by AISdb, users can implement other interpolation techniques tailored to their specific analytical needs. For instance, B-spline (Basis Spline) interpolation is a mathematical technique that creates a smooth curve through data points. This smoothness is important in trajectory analysis as it avoids sharp, unrealistic turns and maintains a natural flow.
Here is an implementation and example of using B-splines interpolation:
Then, we can apply the function just implemented on the vessel tracks generator:
The visualization of the interpolation shows as:
B-spline interpolation with equal distance intervals of 1 km
import aisdb
import numpy as np
import nest_asyncio
from aisdb import DBConn, DBQuery
from datetime import timedelta, datetime
nest_asyncio.apply()
dbpath='YOUR_DATABASE.db' # Define the path to your database
MMSI = 636017611 # MMSI of the vessel
# Set the start and end times for the query
start_time = datetime.strptime("2018-03-10 00:00:00", '%Y-%m-%d %H:%M:%S')
end_time = datetime.strptime("2018-03-31 00:00:00", '%Y-%m-%d %H:%M:%S')
with aisdb.SQLiteDBConn(dbpath=dbpath) as dbconn:
qry = aisdb.DBQuery(
dbconn=dbconn, start=start_time, end=end_time, mmsi = MMSI,
callback=aisdb.database.sqlfcn_callbacks.in_timerange_hasmmsi,
)
rowgen = qry.gen_qry()
tracks = aisdb.track_gen.TrackGen(rowgen, decimate=False)
# Visualize the original track data points
df = track2dataframe(tracks)
plotly_visualize(df, 'scatter')
import numpy as np
from scipy.interpolate import splrep, splev
def bspline_interpolation(track, key, intervals):
"""
Perform B-Spline interpolation for a specific key on the track data.
Parameters:
- track: Dictionary containing vessel track data (time, lat, lon, etc.).
- key: The dynamic key (e.g., 'lat', 'lon') for which interpolation is performed.
- intervals: The equal time or distance intervals at which interpolation is required.
Returns:
- Interpolated values for the specified key.
"""
# Get time and the key values (e.g., lat/lon) for interpolation
times = track['time']
values = track[key]
# Create the B-Spline representation of the curve
tck = splrep(times, values, s=0) # s=0 means no smoothing, exact fit to data
# Interpolate the values at the given intervals
interpolated_values = splev(intervals, tck)
return interpolated_values
def interp_bspline(tracks, step=1000):
"""
Perform B-Spline interpolation on vessel trajectory data at equal time intervals.
Parameters:
- tracks: List of vessel track dictionaries.
- step: Step for interpolation (can be time or distance-based).
Yields:
- Dictionary containing interpolated lat and lon values for each track.
"""
for track in tracks:
if len(track['time']) <= 1:
warnings.warn('Cannot interpolate track of length 1, skipping...')
continue
# Generate equal time intervals based on the first and last time points
intervals = np.arange(track['time'][0], track['time'][-1], step)
# Perform B-Spline interpolation for lat and lon
interpolated_lat = bspline_interpolation(track, 'lat', intervals)
interpolated_lon = bspline_interpolation(track, 'lon', intervals)
# Yield interpolated track
itr = dict(
mmsi=track['mmsi'],
lat=interpolated_lat,
lon=interpolated_lon,
time=intervals # Including interpolated time intervals for reference
)
yield itr
