Next Article in Journal
Drivers’ Psychomotor Reaction Times Tested with a Test Station Method
Next Article in Special Issue
Special Issue on Interdisciplinary Artificial Intelligence: Methods and Applications of Nature-Inspired Computing
Previous Article in Journal
A Complete and Fast Analysis Procedure for Three-Phase Induction Motors Using Finite Element, Considering Skewing and Iron Losses
Previous Article in Special Issue
Integration of Ordinal Optimization with Ant Lion Optimization for Solving the Computationally Expensive Simulation Optimization Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation

by
Rogelio Bautista-Sánchez
1,†,
Liliana Ibeth Barbosa-Santillan
2,*,† and
Juan Jaime Sánchez-Escobar
3
1
Data Mining Engineering Group, Guadalajara 44630, Mexico
2
Department of Computer Science, University of Guadalajara, Guadalajara 44100, Mexico
3
Centro de Enseñanza Técnica Industrial, Guadalajara 44638, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2021, 11(5), 2429; https://doi.org/10.3390/app11052429
Submission received: 1 February 2021 / Revised: 2 March 2021 / Accepted: 2 March 2021 / Published: 9 March 2021

Abstract

:
The prediction of vessel maritime navigation has become an exciting topic in the last years, especially considering economics, commercial exchange, and security. In addition, vessel monitoring requires better systems and techniques that help enterprises and governments to protect their interests. Specifically, the prediction of vessel movements is essential for safety and tracking. However, the applications of prediction techniques have a high cost related to computational efficiency and low resource saving. This article presents a sample method to select historical data on vessel-specific routes to optimize the computational performance of the prediction of vessel positions and route estimation in real-time. These historical navigation data can help to estimate a complete path and perform vessel position predictions through time. This Select Best AIS Data in Prediction Vessel Movements and Route Estimation (PreMovEst) method works in a Vessel Traffic Service database to save computational resources when predictions or route estimations are executed. This article discusses AIS data and the artificial neural network. This work aims to present a prediction model that correctly predicts the physical movement in the route. It supports path planning for the Vessel Traffic Service. After testing the method, the results obtained for route estimation have a precision of 76.15%, and those for vessel position predictions through time have an accuracy of 81.043%.

1. Introduction

Maritime navigation is an essential part of international trade. Global maritime exchange increased exponentially in 2017 [1], mainly driven by Asia and Europe’s economies since many countries export products around the world in large volumes. At the beginning of 2018, the merchant fleet was estimated to be 58,329 vessels [1], which led to more commercial exchange. It has also led to the amount of maritime traffic increasing, thus it is necessary to analyze the information to improve the monitoring processes for different purposes. Naval navigation has experienced a boom. One main task is monitoring the seas to describe and predict moving through them, and, since e-commerce is growing exponentially, vessels are the primary mode of transport for products. Today, governments and enterprises are more interested in it. Data mining techniques (e.g., [2]) aim to discover patterns of movement of the ships from the delimitation of an area.
Statistical analysis and clustering techniques are used. The principles of vessels with similar properties and association rules are used to analyze the movement of the ships and have been implemented to improve the maritime description such as classification, clustering, and predictions [3]. Moreover, where transit routes have been established for navigation, the knowledge of vessel positioning’s historical data allows for the predictive analysis. There are different ways to interpret the information generated, for maritime traffic, density, route estimation, and arrival times.
Likewise, to find vessel movements and prevent incidents, techniques such as maritime tracking and vessel monitoring systems have been proposed. For the latter, in 2002, the International Maritime Organization (IMO) approved the Automatic Identification System (AIS), which has now become an indispensable tool in terms of safety and vessel monitoring around the world.
The vessel industry is volatile, and so commerce and the carriers collaborate and cooperate in alliances. The industry is characterized by having very different freight rates depending on the direction of the trip. Several carriers should deliver the products in several journeys based on the origin country and destination. In 2020, with the COVID-19 pandemic, the number of vessels carrying goods to satisfy client demand has been growing a lot.
The origin and destination of vessels imply costs related to the physical routes; however, the costs should be reduced by selecting the optimal routes.
Besides, to ensure driver adherence to routes and other requirements for the vessel, there should always be a degree of maneuverability in cases that the real time ocean conditions cause changes to origins, destinations, and physical routes.
A vessel can use specific routes many times based on its origin and destination concerning historical vessel positions. Sometimes ships have to change course for different reasons such as adverse climatic agents, political disputes, or prohibitions; however, their trajectories are similar.
Vessel traffic services save historical records to have evidence of the vessel’s location. Industries such as marine traffic, whose system manages large amounts of navigational information for nearly 800 million position records per month, are an example of this.
We define route as the vessel trajectory from the departure port to the arrival port. The position is one set of the coordinates in time of the vessel. The series of positions is all trajectory points of one vessel. Speed is the vessel’s velocity. Course means the direction of the vessel’s grades.
One of the principal goals is the prediction and route estimation of vessel movements. It is a difficult task, however, because systems require high-performance computing. In addition, it requires historical information that can describe patterns in the flow of vessels.
The inferential statistic has generated a significant number of tools to contribute to making scientific judgments on the uncertainty data. Besides, localization measures are designed to provide quantitative values of the central location of the sample. These values will be helpful to identify differences between several of the vessel’s routes. This paper proposes an approach to data selection to reduce computational cost. This work does not focus on precision, but it shows a precision of 76.15% and an accuracy of 81.043% for vessel position predictions through time.
The paper is organized as follows. Section 2 shows a brief description of the related works. Section 3 describes the PreMovEst method for selecting the best AIS data in predicting vessel movements. Section 4 presents preliminary experiments to show the data behavior and adjust the data before selection. Section 5 describes the operations of selecting routes by clustering using the best statistical values. Section 6 presents the results of the experiments. Section 7 summarizes the conclusions drawn from this work.

2. Related Works

The prediction field of vessel trajectory covers different concepts: direction, speed, registered locations, and statistical analysis to identify routes. One of the most critical topics is the way to analyze data to identify patterns. Pallotta et al. [4] proposed a framework with a progressive learning approach and without supervision to extract maritime movement patterns. Using DBSCAN removes outliers that may be location points for other routes. It is a basis for automatically detecting anomalies and projecting current trajectories and patterns into the future.
Perera et al. [5] provided an Extended Kalman Filter (EKF) for the estimation of vessel states and used it for the prediction of vessel trajectories to integrate intelligent features into Vessel Traffic Monitoring and Information Systems (VTMISs); however, the method has high computational complexity for the data selection procedures.
Zissis et al. [6] proposed a simple prediction system using a multilayer function for a web-based system capable of real-time learning and accurately predicting any vessel’s future behavior in low computational time. The selection of the historical AIS data is only related to obtaining the current data. To predict the next location in the future, it uses variables such as direction, course, latitude, longitude, and speed.
Mazzarella et al. [7] proposed a system for the detection and monitoring of multiple vessels. They also showed a proposal for estimating and predicting the navigation trajectory where the historical AIS selection is based on the similarity of the routes around the positions through the clustering technique. The authors used a Bayesian vessel prediction algorithm based on a Particle Filter (PF) to enhance the vessel position prediction quality.
Other works focus on predicting the next position of the vessel based on their neighborhood points and using recursively historical AIS data [8]. Vanneschi et al. [9] predicted the future position of a vessel through a genetic algorithm, and the construction of the dataset is done through the historical information of the same route to be processed. In all cases, the first step is data selection. Zhang et al. [10] selected the dataset based on the destination of the historical path. The most significant similarity to the travel path is predicted as the ship’s destination in the forest model. The essential part of measuring two trajectories’ similarity is the spatial measurement of the historical route (latitude and longitude) for practical applications. DBSCAN is used to find the spatial similarity.
The principal goal to obtain better results in the prediction vessel movements, independent of prediction techniques, is the selection of the correct historical information. For example, Liu considered AIS record’s similarity from the point of origin to the point of arrival and used vector support machines to predict ships’ movements.
Liu [11] used different techniques from Wang et al. [12], who considered a recurring neural network based on two-way recurring gate units. In this work, we use a recurrent neural network based on long short-term memory.
Dobrkovic et al. [13] considered short- and long-term forecasting approaches. The problem is to estimate the arrival times of ships, and they concluded that improvements are need in the quality of the data, the volume of the data, and the mining of these distributed data, as the precision of arrival time estimates, movement predictions, and route estimates depend on these issues. It is necessary to mention that the AIS historical data quality is sometimes not adequate, as mentioned by Dobrkovic et al. A process is needed to improve the quality of the datasets for use in route predictions. It highlights that the non-generalized segments in the routes are considered as trajectory outliers. These are irregular movements of the vessels when they evade barriers, collisions, and erratic traffic. Once the anomaly is determined, the authors proposed to remove it and reconstruct the segment to standardize the data [14].
Alizadeh et al. [15] proposed a vessel trajectory prediction to avoid a collision that reduced the error by 40.85% using neural networks with LSTM. The selection of the historical routes’ data is made considering only with those that show historical movements and are filtered based on the MMSI. Besides, a filter made on the historical data’s similarity based on the current route is predicted. Data selection can be improved to estimate vessel movements, even maritime traffic.
Ramin et al. [16] proposed a procedure based on predicting maritime traffic density using different time series models. They considered selecting the historical AIS data and the temporary labeling of the data for four distinct seasons in the year.
Young [17] estimated the future vessel location, which was tested with the validation of experts. The prediction is done independently for latitude and longitude with neural networks. The data selection is through the area of interest, defining the minimum and maximum latitudes and longitudes. Its MMSI differentiates this block of data for each route. The port makes the selection of the routes of departure and arrival. Similar routes used for prediction are selected with a clustering process to locate the similar routes using the Partition Around the Medoids (PAM) algorithm.
Filipiak et al. [18] developed a system that detects predefined maritime anomalies to support the maritime surveillance system for tanker vessels and specifies the difficulty of analyzing large volumes of data. The filter targets vessels on a large volume of data containing historical and erroneous information. The trajectory is considered an anomaly in the study of circular movements over an area detected by changes in angle and speed.
Xin et al. [19] and Daranda et al. [20] proposed the generation of navigation simulations on a route using statistical analysis. The historical AIS selection is based on only the vessels that travel the route from start to finish and discards those ships that have no similarity in the trajectory. The method highlights using statistical analysis to identify the similarity of the routes that can simulate the trajectory.
Alessandrini et al. [21] built a model that estimates a vessel’s time to arrive in a port. It highlights the AIS historical data selection based on the routes with similarities to the starting location and the destination location related to the arrival time.
Finally, Gao and Shi [22] proposed a model to predict the ships’ movements by identifying patterns in the AIS historical data. It highlights the selection of routes by similar patterns using clustering analysis and statistical classification by samples or indicators in three different ways: using the complete routes, using segments of the routes, and using each of the routes’ positions. They are highlighting as the best options for the segmentation of the routes for later selection by clustering.

3. PreMovEst Method for Select Best AIS Data in Prediction Vessel Movements

The PreMovEst method consists of four components, which are divided into two stages: training and discovering. For the first stage, there are three components: (i) GetAIS data; (ii) historical data collection; and (iii) the selection process of the best-routes collection. For the prediction and route estimation, there is one component: (iv) the process of finding and predicting the position of the vessel movements. We describe each of these steps in the following.
Supervised learning is used to build the knowledge base. It has information linked to the sample containing almost 158,274 records for training using Artificial Neural Networks (ANN). For the application of Multivariate Imputation by Chained Equations (MICE), all records were used, including records of the actual vessel target. In addition, both techniques were used to discover the longitude and latitude in the maritime area [23].
The implementation used Flask and the prediction model used Keras framework, both in Python.
The process of establishing a knowledge base is essential for training. It consists of three stages: (a) obtain the GetAIS data for vessels; (b) obtain the historical data collection; and (c) obtain the best routes collection through a Chi-squared selection process.
i
GetAIS data. Our first assumption is as follows. For the development of this work, the datasets were obtained from MarineCadastre [24] in Zones 15 and 16 in the maritime area with 30 GB volume.
ii
Obtain historical data collection. The method involves accessing MarineCadastre [24] and selecting the year and data segment to download the metadata and their content. The resulting instances are in simplified Dublin Core format. For each obtained instance, a transformation of the Dublin Core format to text file is performed. Historical AIS data can be obtained in different ways. Our process obtains the data through a polygon calculated by departure and destination points, on which our method can be built. After the maritime area is filtered, different AIS data are included. This considerable amount of data can be reduced using a diverse time range in one-month intervals or more. In this way, it is essential to preserve the behavior of our target. If the method considered all AIS data in the maritime area, it could break the general behavior seen in other data.
iii
The best-routes collection through statistical behavior selection process. The method finds values with a higher weight of the mean absolute deviation to be used for classification. They match the statistical target route to detect the best routes. It is not strict because it allows having information about routes with different performance levels and precision. Indices such as precision and recall also give a good measure of performance. Precision refers to the dispersion of the set of values obtained from repeated measurements of a magnitude. The smaller is the dispersion, the greater is the accuracy. Precision is calculated with (1)
Precision = t p t p + f p
Recall is the number of true positives (tp) divided by adding the number of true positives and the number of false positives (fp). True positives are data points classified as positive by the model (meaning they are correct) and false negatives are data points that the model identifies as negative that are positive (incorrect) (2).
Recall = t p t p + f n
This technique allows selecting those characteristics of mean absolute deviation, mean, median, and standard deviation dependent on each other, based on an expected value in the target path that has the absolute values of the prediction of the complete path. This can limit the current use case method since it does not know which is the expected case for each route. The process has to find suitable routes without a sustainable basis of classification or a predefined example.
iv
The process of finding and predicting the position of the vessel movements. Two techniques for prediction are proposed in this paper to show the accuracy of the position of vessel movements prediction for route estimation after selecting all routes that resemble the current one. The first is the use of an Artificial Recurrent Neural Network (ARNN) [25] with Long Short-Term Memory (LSTM) [26,27] using historical data as continual input streams. The second technique employs Multivariate Imputation by Chained Equations (MICE), a statistical method for handling missing data [28].
The network is made up of three layers, with 64 neurons in the input layer, 32 in the intermediate layer, and two neurons for the output layer. The ANN application is based on time-series forecasting, and it is auto-regressive to forecast multiple steps.
MICE offers a significant advantage over other missing data techniques in terms of its flexibility. However, a primary disadvantage is that MICE does not have the same theoretical justification as different imputation approaches. In particular, fitting a series of conditional distributions is done using a series of regression models that are not consistent with proper joint distribution. The purpose of using MICE is to generate the route and an approximation concerning the full path’s prediction based on it.

4. Data Analysis and Preliminary Experiments

The vessel prediction accuracy depends extensively on the existent information’s behavior, what the process obtains, and the appropriate selection of data. Historical AIS data can be obtained in different ways; our process obtains the data through a polygon calculated by the departure and destination points. The data are used to design a polygon, but this can include different AIS data. After filtering the maritime area, a considerable amount of data can be reduced by using a diverse time range, such as one month. In this way, it is crucial to preserve the behavior of our target. If all AIS data are considered in the maritime area, then it could break the general behavior from a different route.

4.1. Statistic Analysis

To understand how the method works, we use AIS data of a specific maritime area. Figure 1 shows some routes, one incomplete from a specific target, which helps to deploy our method.
The GetAIS data process consists of two stages: (a) statistic analysis; and (b) data adjustment. Figure 2 shows several data that exist for our specific target route.
An analysis of each route’s behavior was undertaken using this method to select the correct data. One drawback is the number of samples that a vessel can record on a whole journey. To understand what is happening in this process, we can describe the statistical analysis identifying central tendency and dispersion measures. Our method uses the Freedman–Diaconis rule to obtain a minimal difference between the area under the empirical probability distribution and the area under the theoretical probability distribution. It will help the method to know the behavior of data, so we can see a gap that can significantly differentiate the routes. The second step is to find visual differences among each route using a histogram and Kernel Density Estimation (KDE). If the KDE or the histogram follows a similar pattern for each route, the data behave similarly. However, that does not occur in real environments.

4.2. Data Adjustment

The graph of the histogram with its probability density curve is in two perspectives. The first uses the calculation of bins through the square root of the amount of data. The second is by using the Freedman–Diaconis rule, as shown in Figure 3 and Figure 4. The last graph is a quantile–quantile one to verify the non-normality of the distribution. Figure 3 and Figure 4 show the sample’s behavior around its histogram and its KDE.
Some route data change in latitude or longitude, depending on the destination, starting point, or route. The rest of the data are very similar in latitude and longitude. The data are different in the number of observations, or they have small deviations in the routes. At the same time, the historical data contain complete routes. The previous cases (the feature) affect the search and selection of statistical values that serve as discriminant properties to select the appropriate routes that are the basis for predicting and estimating routes.
The first step is to cut the data between the origin and the last actual position values for each route to adjust the sample route to the objective route.
In the target route, the method approximates all AIS data’s statistical behavior to our target and discards those routes with different behavior because some statistical values can perform as outliers, as shown in Figure 5. The data in the red section are discarded. This step helps approximate the historical data to the current route in statistical terms. Some routes have a different direction, and similar cases help to perform the selection of routes in the next step.
Once the data have been filtered, the method computes the statistical values of each route. All sets of the route with a behavior similar to the target route will maintain a similar statistical behavior despite the lack of information and regardless of the number of observations. This stability is what allows the method to select the most appropriate datasets. It is possible to find statistical values that have a higher weight for classification to detect the best routes that statistically match our route with the Chi-squared selection method. In this case, the method wants to find values with minor and significant statistical differences. The statistics with less dispersion and those with much more dispersion are kept, as shown in Table 1.
The statistical metrics calculated for the vessel position’s target routes (longitude and latitude) are as follows:
The mean is obtained from the tracking of each route. A central descriptive value is obtained from the distribution of the routes. The main disadvantage is the mean is susceptible to outliers’ influence, a property that the method can use as a discriminant feature.
Standard deviation measures the dispersion around the mean. It is highly influenced by the mean and allows the identification of outliers. This property is used to select the closest routes to the route that will generate the prediction.
The median is selected because outliers do not influence it. In addition, it provides an approximate location of the center point to the distribution when it is not normal to the target path distribution.
The trimmed mean is less sensitive to outliers. It allows discarding them and uses the distribution data to get closer to the target route.
MAD is less sensitive to outliers. It is an important variable to relate the routes that are statistically similar to the target route.
The range is susceptible to outliers. This property allows discriminating those routes that have large deviations outside the target route.
Skew identifies the degree of distortion of the distribution and is a characteristic most used as a discriminant.
Kurtosis is a statistical measure to discriminate by focusing on the tails of the distribution. It identifies outliers and routes that have deviations.
Standard Error of Mean (SE) rates the uncertainty in estimating the mean of the target route data compared to the other historical routes.
The selection operation is not strict because the routes can have a statistical similarity from the point of origin to the set point. The consequence is selecting routes with different weights for training and influencing the precision of the method.
Of the previous variables, MAD is one of the most stable variables depending on its nature, regardless of the number of observations, as shown in Figure 6.
Therefore, it is the variable that determines the selection of significant routes. The range is quite sensitive to the characteristics of the dataset because it is altered significantly by the existence of outliers.
As mentioned above, Table 1 shows the statistical values of each route as features. The data have a marked difference in most of the statistical values, specifically those related to latitude. The method applies the cluster technique to those statistical values to help us choose the best prediction routes. The cluster is on densities. The dense region of objects contains a target path and similar paths. The low-density section is those routes with a statistical difference.

5. Experiments

This section uses input data analysis and the selection of those statistical values that represent the best way for each route. A clustering technique was used and the data were selected to predict the route.
The process of finding and predicting the position of the vessel movements consists of three stages: (a) route selection by clustering; (b) using an artificial neural network; and (c) using multivariate imputation by chained equations.

5.1. Route Selection by Clustering

Density-Based Spatial Clustering of Applications with Noise (DBSCAN) [29] was applied to discover and select the best data. DBSCAN was proposed by Albert Ester and can identify clusters and outliers. A clear difference before the application of DBSCAN is found in the sample.
The first step is to cut the 146 data (for each route) between the origin and the last actual position values (147 longitudes and latitudes approximately) of the target route to adjust the sample routes to the objective route.
The route selection helps approximate all 148 AIS data’s statistical behavior to the target and discards those routes with different behavior. Some of the 149 statistical values are outliers, as shown in Figure 7.
The configuration parameters of DBSCAN are related to the minimum distance found at each point, and the minimum elements of clusters are defined by 3. The method applies DBSCAN to see several groups within its routes, as shown in Figure 8.
After applying DBSCAN, the method discards all routes that do not resemble the route to predict. For better results, it is necessary to have enough data to get an accurate approximation of its best statistical values to achieve the classification’s performance according to the theorem of limit central [30]. The algorithm chooses the information that justifies the use of specific routes that share similarities with the target.
The method’s comparison is with the methodological approach for extraction of the characteristics of biological signals [31], which uses the histogram class marks to make the selection of characteristics using the behavior of the distribution of biological signals.
The classmark allows a representative point of each histogram interval to be obtained and then a frequency polygon to be created. A frequency polygon is built with a dynamic prototype and biological signals to get a graphic form to represent the original data.
Table 2 shows the comparative results of route selection for each method, which are almost all similar. The differences are in the features. PreMovEst can use a small number of features to select routes. In contrast, the methodological approach for the extraction of biological signals’ characteristics needs all features, multiplying the number of latitude and longitude features.
For prediction and route estimation, the selection of all routes that resemble the current route was made. First, an Artificial Recurrent Neural Network (ARNN) with Long Short-Term Memory (LSTM) was used with historical data as continual input streams. As the second technique, multivariate imputation by chained equations used a statistical method for handling missing data.
The sample contains almost 158,274 records for training in the case of an Artificial Neural Network (ANN). For the application of MICE, all records were used, including records of the actual vessel target.
In addition, for both techniques, only longitude and latitude were used to save computational power.

5.2. Using an Artificial Neural Network

The type of ANN used before was an ARNN with a Long Short-Term Memory layer using Keras high-level framework [32].
First, some variables were defined: the network is made up of three layers, with 64 neurons in the input layer, 32 in the intermediate layer, and two neurons for the output layer, as shown in Figure 9.
The method learns through an adaptive prediction model and automatically determines when the next position on a vessel’s route occurs. A knowledge base was generated with all routes of the target vessel to achieve this goal. The first part of the data was used to create a model to predict the routes. Besides, the filtered historical data were used to train, validate, and test the neural network model to predict the current route.

5.3. Using Multivariate Imputation by Chained Equations

MICE is widely used to search for missing dataset values to get the best cases about data behavior [33].
For the full application of this technique, it was necessary to carry out a stratified sampling. The algorithm obtains a representative matrix of the data to get the most important information that contributes to generating the estimation route.

6. Results

The purpose was to show that, after obtaining the journey knowledge base, our method helps in the selection of the routes based on the statistical behavior of each one. ANN achieved the prediction of an interval of time and the estimation of the route using the MICE Algorithm. The PreMovEst method obtained the following results: Table 3 shows a contrast between the actual data of an interval of time and the ANN’s vessel movement predictions.
Table 3 presents ten positions of Route 25 of Segments 15 and 16. The results show an accuracy of 80.5% in the prediction. However, if the method does not have enough information about the route to be estimated or it is too short, it loses precision. That is why we combine the two processes: prediction using neural networks and route estimation with MICE. Besides, where there is not enough information, it is necessary to obtain more historical data, which increases the processing time that depends on the amount of AIS.
The amount of information used in the second example was more than in the first one. That allowed the accuracy to be better, at 84%. The historical AIS data samples are similar from the block to be estimated, out of 10 predictions. The absolute mean difference for latitude is 0.0015655 degrees, while the absolute mean difference for longitude is 0.00211949 degrees, as shown in Table 4.
Concerning the MICE technique’s application, as shown in Figure 10, the prediction generates the route approximation with 76.15% accuracy, using a small sample of routes. Finally, the prediction obtained an accuracy of 81.043%. However, it allows for approximating the results with the ten routes as the sample.
The PreMovEst method shows that the prediction is almost the same as the real one after selecting the best routes.
The best routes (series of positions) have a similarity in statistical properties after the clustering process.
The routes’ selection through their statistical data in the clustering allows discriminating those routes whose navigation course differs from the current route predicted. That is why they work as outliers, which is evident when clustering carried out, as they remain outside the cluster to which the route to be estimated belongs.
Figure 10 compares the PreMovEst method to the Zissis method that considers a cloud infrastructure to support marine traffic. Furthermore, their system stores the trained vessel and recalls it 24 h as input for the following predictions. Our method uses data in real-time and is capable of running on a computer with GPU. Computing power is demanded, depending on the application. Our method contributes to the search for processes that allow computational savings. It is understood that a dataset can be better known if it has the largest number of characteristics.
In the second example (see Figure 11), the method used more historical AIS data samples, i.e., 28% more than the first example. There is more information to predict the next movements of the route. In addition, the complete route was estimated. In this example, the historical samples are quite similar to the target route. It is also favorable because the route block to be estimated is smaller than in the previous example.
A segment of the original route with no information was observed. It occurred for technical reasons and was compensated by the historical information extracted in the selection process; however, our objective is not to recover missing data, which can lead to missing historical information. Thus, the neural network’s prediction process would be more uncertain when there are missing data in the segment to be estimated.

7. Conclusions

The present proposal, called Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation (PreMovEst), is related to obtaining reports on vessels’ status that allows for accurate and real-time monitoring.
It applies to all vessels navigating within a designated maritime area and predicting future positions of a vessel based on positions already traveled within the designated maritime area. The information is shown to the user through general and specific vessel status reports.
The PreMovEst method was tested using a small dataset. However, the correct selection of data for analysis and prediction in maritime routes is essential before applying the method.
The PreMovEst method uses a filter for the routes whose behavior are similar to obtain the best routes collection through a Chi-squared selection process. The execution time was reduced from 16 to 11 min on average. The experiments were run on a Core i7 processor with 16 GB of RAM.
The results show a prediction accuracy of 80.5–84%. The amount of information used in the second example was more than the first example, which allowed the accuracy to be better, at 84%. The historical AIS data samples are similar to the block to be estimated, out of 10 predictions. The absolute mean difference for latitude is 0.0015655 degrees, while the absolute mean difference for longitude is 0.00211949 degrees.
The results concern the assets of the vessels that move in the seas, especially in economic terms. The paper shows the PreMovEst method’s application and testing of its accuracy using ARNN and MICE techniques. Everyday, vessels move billions of items across the oceans from one country to another, a drawback of which is handling all of this information. The process of selecting data to train the neural network is a function of the number of trips made by the vessel from the port of departure to the arrival port.
One of the challenges to be undertaken is to obtain better accuracy while reducing the execution time.
Future work will be to accelerate the prediction model with Graphical Process Units (GPU). This requires image processing, which implies supercomputers to process such data, limiting its effectiveness for real-time monitoring of vessels navigating within a designated maritime area and increasing the amount of processed information.

Author Contributions

Conceptualization, R.B.-S., and J.J.S.-E.; software, R.B.-S., L.I.B.-S. and J.J.S.-E. investigation, R.B.-S., L.I.B.-S. and J.J.S.-E. writing—original draft preparation, R.B.-S., L.I.B.-S. and J.J.S.-E.; writing R.B.-S., L.I.B.-S. and J.J.S.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was founding by Research Council (CONACyT), Mexico.

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author.

Acknowledgments

This research was supported by the Sciences Research Council (CONACyT) through the research project number 262756 called “The use of GNSS data for tracking maritime flow for sea security” http://navigationgnssproject.net/index.html accessed on 28 February 2021.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ANAVE. Merchant Marine and Maritime Transport 2017/2018. 2018. Available online: https://unctad.org/system/files/official-document/rmt2018_en.pdf (accessed on 9 December 2020).
  2. Deng, F.; Guo, S.; Deng, Y.; Chu, H.; Zhu, Q.; Sun, F. Vessel track information mining using AIS data. In Proceedings of the 2014 International Conference on Multisensor Fusion and Information Integration for Intelligent Systems (MFI), Beijing, China, 28–29 September 2014; pp. 1–6. [Google Scholar] [CrossRef]
  3. Alessandrini, A.; Alvarez, M.; Greidanus, H.; Gammieri, V.; Arguedas, V.F.; Mazzarella, F.; Santamaria, C.; Stasolla, M.; Tarchi, D.; Vespe, M. Mining Vessel Tracking Data for Maritime Domain Applications. In Proceedings of the IEEE International Conference on Data Mining Workshops (ICDM 2016), Barcelona, Spain, 12–15 December 2016; pp. 361–367. [Google Scholar] [CrossRef]
  4. Pallotta, G.; Vespe, M.; Bryan, K. Vessel pattern knowledge discovery from AIS data: A framework for anomaly detection and route prediction. Entropy 2013, 15, 2218–2245. [Google Scholar] [CrossRef] [Green Version]
  5. Perera, L.P.; Oliveira, P.; Guedes Soares, C. Maritime Traffic Monitoring Based on Vessel Detection, Tracking, State Estimation, and Trajectory Prediction. IEEE Trans. Intell. Transp. Syst. 2012. [Google Scholar] [CrossRef]
  6. Zissis, D.; Xidias, E.K.; Lekkas, D. Real-time vessel behavior prediction. Evol. Syst. 2016, 7, 29–40. [Google Scholar] [CrossRef]
  7. Mazzarella, F.; Arguedas, V.F.; Vespe, M. Knowledge-based vessel position prediction using historical AIS data. In Proceedings of the 2015 Workshop on Sensor Data Fusion: Trends, Solutions, Applications (SDF), Bonn, Germany, 6–8 October 2015. [Google Scholar] [CrossRef]
  8. Hexeberg, S.; Flåten, A.L.; Eriksen, B.H.; Brekke, E.F. AIS-based vessel trajectory prediction. In Proceedings of the 2017 20th International Conference on Information Fusion (Fusion), Xi’an, China, 10–13 July 2017; pp. 1–8. [Google Scholar] [CrossRef]
  9. Vanneschi, L.; Castelli, M.; Re, A. Prediction of ships’ position by analysing AIS data: An artificial intelligence approach. Int. J. Web Eng. Technol. 2017, 12, 253. [Google Scholar] [CrossRef]
  10. Zhang, C.; Bin, J.; Wang, W.; Peng, X.; Wang, R.; Halldearn, R.; Liu, Z. AIS data driven general vessel destination prediction: A random forest based approach. Transp. Res. Part C Emerg. Technol. 2020, 118, 102729. [Google Scholar] [CrossRef]
  11. Liu, J.; Shi, G.Z.K. Vessel Trajectory Prediction Model Based on AIS Sensor Data and Adaptive Chaos Differential Evolution Support Vector Regression (ACDE-SVR). Appl. Sci. 2020, 9, 15. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, C.; Ren, H.; Li, H. Vessel trajectory prediction based on AIS data and bidirectional GRU. In Proceedings of the 2020 International Conference on Computer Vision, Image and Deep Learning (CVIDL), Nanchang, China, 15–17 May 2020; pp. 260–264. [Google Scholar] [CrossRef]
  13. Dobrkovic, A.; Iacob, M.; van Hillegersberg, J.; Mes, M.; Glandrup, M. Towards an Approach for Long Term AIS-Based Prediction of Vessel Arrival Times. In Logistics and Supply Chain Innovation: Bridging the Gap between Theory and Practice; Lecture Notes in Logistics; Zijm, H., Klumpp, M., Clausen, U., ten Hompel, M., Eds.; Springer: Berlin, Germany, 2015; pp. 281–294. [Google Scholar] [CrossRef]
  14. Tu, E.; Zhang, G.; Mao, S.; Rachmawati, L.; Huang, G.B. Modeling Historical AIS Data For Vessel Path Prediction: A Comprehensive Treatment. arXiv 2020, arXiv:2001.01592. [Google Scholar]
  15. Alizadeh, D.; Alesheikh, A.A.; Sharif, M. Vessel Trajectory Prediction Using Historical Automatic Identification System Data. J. Navig. 2021, 74, 156–174. [Google Scholar] [CrossRef]
  16. Ramin, A.; Mustaffa, M.; Ahmad, S. Prediction of Marine Traffic Density Using Different Time Series Model From AIS data of Port Klang and Straits of Malacca. Trans. Marit. Sci. 2020, 9, 217–223. [Google Scholar] [CrossRef]
  17. Young, B.L. Predicting Vessel Trajectories from AIS Data Using R. Master’s Thesis, Naval Postgraduate School, Monterey, CA, USA, 2017. [Google Scholar]
  18. Filipiak, D.; Stróżyna, M.; Węcel, K.; Abramowicz, W. Big Data for Anomaly Detection in Maritime Surveillance: Spatial AIS Data Analysis for Tankers. Zesz. Nauk. Akad. Mar. Wojennej 2018, 215, 5–28. [Google Scholar] [CrossRef] [Green Version]
  19. Xin, X.; Liu, K.; Yang, X.; Yuan, Z.; Zhang, J. A simulation model for ship navigation in the “Xiazhimen” waterway based on statistical analysis of AIS data. Ocean. Eng. 2019, 180, 279–289. [Google Scholar] [CrossRef]
  20. Daranda, A.; Dzemyda, G. Navigation Decision Support: Discover of Vessel Traffic Anomaly According to the Historic Marine Data. Int. J. Comput. Commun. Control. 2020, 15. [Google Scholar] [CrossRef] [Green Version]
  21. Alessandrini, A.; Mazzarella, F.; Vespe, M. Estimated Time of Arrival Using Historical Vessel Tracking Data. IEEE Trans. Intell. Transp. Syst. 2019, 20, 7–15. [Google Scholar] [CrossRef]
  22. Gao, M.; Shi, G.Y. Ship-handling behavior pattern recognition using AIS sub-trajectory clustering analysis based on the T-SNE and spectral clustering algorithms. Ocean Eng. 2020, 205, 106919. [Google Scholar] [CrossRef]
  23. Bautista, R.; Barbosa-Santillán, L.; Sánchez-Escobar, J. Statistical Approach in Data Filtering for Prediction Vessel Movements Through Time and Estimation Route Using Historical AIS Data. In Proceedings of the Mexican International Conference on Artificial Intelligence, Xalapa, Mexico, 27 October–2 November 2019; pp. 28–38. [Google Scholar] [CrossRef]
  24. U.S. Department of Commerce’s National Oceanic and Atmospheric Administration (NOAA) Office for Coastal Management and the U.S. Department of the Interior’s Bureau of Ocean Energy Management (BOEM). MarineCadastre. 2019. Available online: https://marinecadastre.gov/ais/ (accessed on 9 December 2020).
  25. Williams, R.J.; Zipser, D. A Learning Algorithm for Continually Running Fully Recurrent Neural Networks. Neural Comput. 1989, 1, 270–280. [Google Scholar] [CrossRef]
  26. Hochreiter, S.; Urgen Schmidhuber, J. LONG SHORT-TERM MEMORY. Neural Comput. 1997, 9, 1735–1780. [Google Scholar] [CrossRef] [PubMed]
  27. Gers, F.A.; Schmidhuber, J.; Cummins, F. Learning to forget: Continual prediction with LSTM. Neural Comput. 2000. [Google Scholar] [CrossRef] [PubMed]
  28. White, I.R.; Royston, P.; Wood, A.M. Multiple imputation using chained equations: Issues and guidance for practice. Stat. Med. 2011, 30, 377–399. [Google Scholar] [CrossRef] [PubMed]
  29. Ester, M.; Kriegel, H.P.; Sander, J.; Xu, X. A Density-based Algorithm for Discovering Clusters a Density-based Algorithm for Discovering Clusters in Large Spatial Databases with Noise. In Proceedings of the Second International Conference on Knowledge Discovery and Data Mining (KDD’96), Portland, OR, USA, 2–4 August 1996; pp. 226–231. [Google Scholar]
  30. Peck, R.; Devore, J.L. Statistics: The Exploration & Analysis of Data; Cengage Learning: Boston, MA, USA, 2011. [Google Scholar]
  31. Granados-Ruiz, J.; Hernandez, D.A.G.; Ramírez, C.L.; Zamudio, V.; Pérez-Mata, P. Methodological approach for extraction of characteristics of biological signals. Compusoft 2019, 8. [Google Scholar]
  32. Chollet, F. Keras Documentation. 2015. Available online: https://github.com/fchollet/keras (accessed on 9 December 2020).
  33. Azur, M.J.; Stuart, E.A.; Frangakis, C.; Leaf, P.J. Multiple imputation by chained equations: What is it and how does it work? Int. J. Methods Psychiatr. Res. 2011. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Routes in the selected area with different historical ways. The actual vessel location is at the center.
Figure 1. Routes in the selected area with different historical ways. The actual vessel location is at the center.
Applsci 11 02429 g001
Figure 2. The actual route for a specific target-destination is marked with X.
Figure 2. The actual route for a specific target-destination is marked with X.
Applsci 11 02429 g002
Figure 3. The Kernel Density Estimation (KDE) and histograms for each sample route (one and two) to visualize the differences. Freedman–Diaconis rule was used to calculate bins.
Figure 3. The Kernel Density Estimation (KDE) and histograms for each sample route (one and two) to visualize the differences. Freedman–Diaconis rule was used to calculate bins.
Applsci 11 02429 g003
Figure 4. Route 3 and 4 Kernel Density Estimation (KDE) and histograms for each sample route to visualize the differences. Freedman–Diaconis rule is used as a complement of Figure 3.
Figure 4. Route 3 and 4 Kernel Density Estimation (KDE) and histograms for each sample route to visualize the differences. Freedman–Diaconis rule is used as a complement of Figure 3.
Applsci 11 02429 g004
Figure 5. The data in the red section are discarded.
Figure 5. The data in the red section are discarded.
Applsci 11 02429 g005
Figure 6. A three-dimensional perspective of three features of the routes, two of which are MAD.
Figure 6. A three-dimensional perspective of three features of the routes, two of which are MAD.
Applsci 11 02429 g006
Figure 7. The actual and predicted routes.
Figure 7. The actual and predicted routes.
Applsci 11 02429 g007
Figure 8. Cluster discovered using DBSCAN. Those routes in the group (tagged as Similar) have a similarity in statistical properties and help predict and estimate the route-target.
Figure 8. Cluster discovered using DBSCAN. Those routes in the group (tagged as Similar) have a similarity in statistical properties and help predict and estimate the route-target.
Applsci 11 02429 g008
Figure 9. The neural network architecture.
Figure 9. The neural network architecture.
Applsci 11 02429 g009
Figure 10. The route prediction results with MICE.
Figure 10. The route prediction results with MICE.
Applsci 11 02429 g010
Figure 11. The route prediction results with MICE.
Figure 11. The route prediction results with MICE.
Applsci 11 02429 g011
Table 1. All statistic values are computed for latitude and longitude, and Route 25 is our route-target.
Table 1. All statistic values are computed for latitude and longitude, and Route 25 is our route-target.
RoutesLat MeanLat SDLAT MedianLat TrimmedLat MadLat_Range RangeLat SkewLat KurtosisLat SELon MeanLon SDLon MedianLon TrimmedLon MADLon RangeLon SkewLon KurtosisLon SE
Route 146.74520.369646.720646.7430.42751.3270.1046−1.23370.0129−85.18770.9905−84.8746−85.10191.01523.2322−0.5702−1.06050.0344
Route 246.7490.381.04346.705646.74550.43581.35470.1278−1.23770.0124−85.17370.9745−84.8458−85.0830.93563.2322−0.6222−0.99580.0318
Route 346.50890.006546.512446.50940.00450.0261−0.4288−1.37781 × 10 4 −84.37570.0169−84.3752−84.37610.00460.09230.16080.5053 × 10 4
Route 446.44350.113346.503246.47200.4378−1.86312.24550.0034−84.281.040.0884−84.327−84.300800.32881.81422.04840.0026
Route 546.81.040.348146.852246.80650.47481.3192−0.181.04−1.12270.01−85.30.9396−85.1586−85.24751.20063.2307−0.347−1.14760.0271
Route 646.42540.125646.503246.451500.4374−1.39480.5380.0026−84.26430.0986−84.327−84.283700.33081.32050.30280.0021
Route 746.48370.068746.501946.50191 × 10 4 0.4375−4.1181.016.64747 × 10 4 −84.38010.0729−84.4039−84.40156 × 10 4 0.40353.454911.55517 × 10 4
Route 846.89390.48146.878646.88440.59431.69870.1644−1.24340.0158−85.26110.9821−85.019−85.19141.11283.2296−0.4605−1.13270.0322
Route 946.75120.423246.594546.72960.46331.83480.461−0.7810.011−84.9810.8673−84.73−84.860.72283.2326−1.0074−0.05120.0225
Route 1046.7780.40646.81446.77390.5061.43110.0993−1.30.0138−85.23061.0034−84.9672−85.15041.11483.2241−0.4968−1.181.04330.0342
Route 1146.92460.522146.846646.92440.67561.63680.0926−1.61750.0159−85.14240.9692−84.881.04−85.04520.93623.2351−0.6855−0.87720.0296
Route 1246.76350.411146.571246.7470.4321.50040.3806−1.11570.0123−85.11270.9484−84.6995−85.00810.71913.2364−0.7724−0.781.040.0285
Route 1346.49070.056946.504146.504100.4487−5.0681.0426.41594 × 10 4 −84.3170.0461−84.3269−84.326900.39944.503122.23073 × 10 4
Route 1546.62360.249546.618546.63580.28270.9644−0.3162−0.93310.0083−85.11460.9591−84.7612−85.01260.81023.2341−0.7423−0.81970.0318
Route 1647.25690.073147.257647.25710.09370.2516−0.0208−1.21840.0055−86.75520.2824−86.7573−86.7560.36080.97160.0199−1.21920.0212
Route 1746.56890.184746.580146.58640.18820.7522−0.6377−0.2190.0058−85.2571.0155−84.9607−85.18481.12443.2318−0.4296−1.2650.0318
Route 1846.5039046.503946.503906 × 10 4 0.34262.10180−84.3270−84.327−84.32704 × 10 4 −0.62511.64650
Route 2146.71510.352246.529546.70470.38311.32610.3276−1.07410.0116−85.07430.9629−84.6554−84.96440.70263.2318−0.81.043−0.77220.0316
Route 2246.84930.522646.748146.80320.45952.14140.80190.10720.0155−84.93080.7366−84.8472−84.83970.73673.2328−1.20271.49060.0218
Route 2346.74840.368646.723246.74620.42311.32520.1032−1.24290.0129−85.19560.9924−84.8829−85.11051.01583.232−0.5567−1.08340.0346
Route 2446.59170.175146.603446.60970.16270.7872−0.81530.46080.0048−85.09130.8171−84.9781−84.9850.75543.2322−1.04730.30760.0225
Route 2546.78960.344146.845646.81.040.47131.2992−0.1852−1.12930.01−85.27610.9208−85.1329−85.22591.16763.1519−0.3408−1.1630.0267
Table 2. The results of route selection for all methods are similar.
Table 2. The results of route selection for all methods are similar.
Clustering after Using PreMovTestClustering after Using Granados-Ruiz 2019Route
11Route 1
11Route 2
00Route 3
00Route 4
11Route 5
00Route 6
00Route 7
10Route 8
10Route 9
11Route 10
10Route 11
10Route 12
00Route 13
10Route 15
00Route 16
00Route 17
00Route 18
11Route 21
00Route 22
11Route 23
00Route 24
11Route 25
Table 3. Real data versus the predicted path segment. The predicted latitude has better accuracy than the predicted longitude.
Table 3. Real data versus the predicted path segment. The predicted latitude has better accuracy than the predicted longitude.
Original LatitudePredicted LatitudeOriginal LongitudePredicted Longitude
47.3676747.88345−87.16083−86.99632
47.3688747.99324−87.16616−86.54012
47.3701848.01275−87.17136−86.62076
47.371548.02375−87.17667−86.56543
47.3728248.11034−87.18205−87.18984
47.3739648.11498−87.1872−87.22098
47.3753248.12876−87.19258−87.35465
47.3766748.12076−87.19767−87.65703
47.37848.13776−87.203−87.98583
47.381.0433348.22021−87.20817−87.99232
Table 4. Real data contrasted with the previously predicted path segment on the second sample.
Table 4. Real data contrasted with the previously predicted path segment on the second sample.
Original LatitudePredicted LatitudeOriginal LongitudePredicted Longitude
47.3676747.88345−87.16083−86.99632
47.3688747.99324−87.16616−86.54012
47.3701848.01275−87.17136−86.62076
47.371548.02375−87.17667−86.56543
47.3728248.11034−87.18205−87.18984
47.3739648.11498−87.1872−87.22098
47.3753248.12876−87.19258−87.35465
47.3766748.12076−87.19767−87.65703
47.37848.13776−87.203−87.98583
47.381.0433348.22021−87.20817−87.99232
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bautista-Sánchez, R.; Barbosa-Santillan, L.I.; Sánchez-Escobar, J.J. Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation. Appl. Sci. 2021, 11, 2429. https://doi.org/10.3390/app11052429

AMA Style

Bautista-Sánchez R, Barbosa-Santillan LI, Sánchez-Escobar JJ. Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation. Applied Sciences. 2021; 11(5):2429. https://doi.org/10.3390/app11052429

Chicago/Turabian Style

Bautista-Sánchez, Rogelio, Liliana Ibeth Barbosa-Santillan, and Juan Jaime Sánchez-Escobar. 2021. "Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation" Applied Sciences 11, no. 5: 2429. https://doi.org/10.3390/app11052429

APA Style

Bautista-Sánchez, R., Barbosa-Santillan, L. I., & Sánchez-Escobar, J. J. (2021). Method for Select Best AIS Data in Prediction Vessel Movements and Route Estimation. Applied Sciences, 11(5), 2429. https://doi.org/10.3390/app11052429

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop