Research Article  Open Access
Eftychios Protopapadakis, Athanasios Voulodimos, Anastasios Doulamis, Nikolaos Doulamis, Dimitrios Dres, Matthaios Bimpas, "Stacked Autoencoders for Outlier Detection in OvertheHorizon Radar Signals", Computational Intelligence and Neuroscience, vol. 2017, Article ID 5891417, 11 pages, 2017. https://doi.org/10.1155/2017/5891417
Stacked Autoencoders for Outlier Detection in OvertheHorizon Radar Signals
Abstract
Detection of outliers in radar signals is a considerable challenge in maritime surveillance applications. HighFrequency SurfaceWave (HFSW) radars have attracted significant interest as potential tools for longrange target identification and outlier detection at overthehorizon (OTH) distances. However, a number of disadvantages, such as their low spatial resolution and presence of clutter, have a negative impact on their accuracy. In this paper, we explore the applicability of deep learning techniques for detecting deviations from the norm in behavioral patterns of vessels (outliers) as they are tracked from an OTH radar. The proposed methodology exploits the nonlinear mapping capabilities of deep stacked autoencoders in combination with densitybased clustering. A comparative experimental evaluation of the approach shows promising results in terms of the proposed methodology’s performance.
1. Introduction
Detection of targets and outliers in radar signals is a research issue that has gained significant attention in the academic and industrial research community, mainly because of the important associated impact of relevant applications in surveying of large areas. HighFrequency SurfaceWave (HFSW) radars are a category of radars that operate at the frequency band 3–30 MHz and, in contrast with other radars, use ground wave or sky wave propagation and ionospheric reflections of the electromagnetic waves for target detection, which allows for achieving longer ranges, where microwave radars cannot perform [1], but to the detriment of the attained accuracy. For many years, HFSW radars, or overthehorizon (OTH) radars, as they are commonly known, have been used to remotely measure oceanographic parameters, providing information about surface currents, wave spectra, wind direction and intensity, and so on [2]. Their extraordinary range (up to 200 nautical miles) combined with their continuous mode of operation make for an ideal candidate tool for longrange oceanic surveillance. However, many associated weaknesses, for example, low spatial resolution, high nonlinearity, and important presence of clutter, negatively impact their performance as earlywarning tools for detection, tracking, and identification of vessels.
The promising capabilities of OTH radars have attracted significant interest from the research community and have already resulted in various approaches (e.g., [3, 4]). Nevertheless, related research issues continue to present significant challenges, which can be attributed to few reasons, briefly described below:(i)Different targets may present similar dielectric and frequency properties thus making it hard to make a clear distinction among them.(ii)Given multipath propagation effects of rough surfaces, scattering from some objects tends to overwhelm the weak backscattering of targets.(iii)Due to the changes in atmosphere and ground conditions, noise is added which can confuse the analysis of a radar signal.(iv)Ocean and ionospheric clutter generate noise especially for HFSW radars.On a different note, the surge of deep learning and the great results it has produced in other signal analysis domains, such as computer vision, speech recognition, and natural language processing, create certain expectations regarding its potential efficacy in radar signal analysis applications. Deep learning allows computational models of multiple processing layers to learn and represent data with multiple levels of abstraction mimicking how the brain perceives and processes multimodal information, thereby implicitly capturing intricate structures of large‐scale data. Complex abstractions are learnt at a given level based on relatively simpler abstractions formulated in the preceding layer in the hierarchy.
The goal of this paper is to present a framework for detecting deviations from the norm in behavioral patterns of vessels (henceforth called outliers), as they are tracked from an OTH radar. The proposed methodology exploits the nonlinear mapping capabilities of deep stacked autoencoders (SAs) [5] in combination with densitybased clustering. Stacked autoencoders are used in an unsupervised way to map the track history of any vessel into a compact and informative feature vector. Then, at any moment all tracked ships are projected into a new feature space and clustered using densitybased algorithms, such as OPTICS [6]. The outcome of the clustering stage then indicates possible outliers.
The remainder of this paper is structured as follows: Section 2 presents an overview of the related work. In Section 3 we describe in detail the proposed methodology for outlier detection in OTH radar signals, which is followed by the experimental evaluation of the methodology in Section 4. Finally, Section 5 concludes the paper.
2. Related Work
In the literature, several signal processing and machine learning methods have been investigated and proposed to acquire more reliable data with lower noise and extract semantic information from radar signals. Kouemou and Opitz [7] introduced a waveletbased feature analysis combined with Hidden Markov Models (HMM) to classify real radar signals into predefined categories. Spectral analysis [8] is used by GarbanzoSalas and Hocking [9] for detecting small objects from harmonic pulse radar data. The use of online bootstrapping machine learning tools to improve target detection rate of radar signals is also one major research area [10]. Radar data can be analyzed using the concepts of transfer learning since often we have only a small number of labelled data available while the majority of signals captured are unlabelled (nonannotated) [11]. Other works focus on modeling of ionospheric disturbances on spaceborne interferometric synthetic aperture radar (SAR) via EchoState Networks [12, 13] or ensemble classifiers [14].
Denoising techniques for radar signals include low level processing such as the median filter or other nonlinear convolution schemes [15]. Other approaches spatially or temporally decompose radar signals by wavelet transforms [16, 17]. This way, we can find patterns distributed on space and time domain to improve targets detection efficiency. These methodologies can be extended to the analysis of synthetic aperture radar (SAR) images [11], or by incorporating sparsitybased signal analysis [18]. A neural network based scheme for detecting salient objects in SAR images is recently presented [19]. The goal is to identify changes in SAR content. A similar approach for detecting changes using nonlinear stacked restricted Boltzmann machines is given in the work of Liu et al. [20], while multilayered feature learning to improve detection accuracy of SAR images is described by Xie et al. [21]. Furthermore, lowpower HF surfacewave (HFSW) radars have demonstrated being a costeffective longrange earlywarning sensor for ship detection and tracking [22, 23]. A detailed description of various ways in which HFSW radar technology can be used for maritime surveillance is provided by Braca et al. [24].
Regarding deep versus “shallow” learning schemes, traditional machine learning techniques exploit shallow architectures; that is, they use a single layer for data/feature transformation, even in a highly nonlinear space. Shallowness refers here to the simplicity of these architectures that use only one (or few) layer(s) of processing, responsible for transforming the raw input signals or features into the problemspecific feature space. Instead, in a deep learning paradigm, the architectures are composed of many (deep) nonlinear processing stages [25]. Deep learning has been extensively applied in many fields, such as computer vision [26] (e.g., behavior recognition [27] and human tracking [28]) and speech recognition [29]. However, its applicability in radar signal processing had not being investigated until very recently [30]. Even so, most of the proposed works pertain to object detection in SAR image data [31], essentially resembling visual analysis approaches.
3. The Proposed Methodology
The proposed methodology exploits the nonlinear mapping abilities of stacked autoencoders (SAs) [5] in combination with densitybased clustering, to identify irregular occurrences, using overthehorizon radar data. Such an approach is based on two main assumptions:(1)The history of a naval vessel, in terms of speed, position, course, signal frequency, or other related data, provided by a ground radar, suffices to extract meaningful features.(2)Unexpected deviation from the norm is observed for a few ships, denoted henceforth as outliers.The approach is relatively straightforward: Given a set of OTH data entries, SAs are used in an unsupervised way to map the track history of any vessel into a compact and informative feature vector. Then, at any moment all tracked ships are projected into a new feature space and clustered using OPTICS [6], a widely used densitybased algorithm. The clustering outcome informs about possible outliers. In the following subsections, the different stages of the methodology are presented, after a brief description of the data involved. Figure 1 provides a highlevel view of the proposed approach.
3.1. OTH and AIS Data
Heterogeneous data, such as automatic identification system (AIS) data, highfrequency surface wave (HFSW) radar data, and synthetic aperture radar (SAR) data, have been exploited in research for maritime surveillance purposes [32]. In our case, two sources of information were fused to support the outlier detection process: OTH radar and AIS data.
The OTH radar data used for the setting and evaluation of the presented work was acquired by the HFSW STRADIVARIUS radar by Diginext [33]. OTH radar detection (plot) and tracking (track) data are the output of the OTH radar for a given period. The plot and track data provided include estimated position coordinates, velocity, course, Doppler frequency, global and local noise levels, azimuth, and other parameters, appropriately timestamped.
On a different note, AIS is an automatic tracking system used for collision avoidance on ships and by vessel traffic services. AIS information supplements marine radar, which continues to be the primary method of collision avoidance for water transport. Vessels equipped with AIS transceivers can be tracked by AIS base stations located along coast lines. The International Maritime Organization’s International Convention for the Safety of Life at Sea requires AIS to be present aboard international voyaging ships with gross tonnage of 300 or more and all passenger ships regardless of size [34]. AIS reports contain both dynamic information (e.g., latitude, longitude, course over ground, speed over ground, and time) and static information (e.g., vessel type and dimension information).
3.2. DensityBased Clustering as a Basis for Outlier Detection
Clustering refers to the task of identifying groups or clusters in a dataset. In densitybased clustering, a cluster is a set of data objects spread in the data space over a contiguous region of high density of objects. Densitybased clusters are separated from each other by contiguous regions of low density of objects. Data objects located in lowdensity regions are typically considered noise or outliers [35]. OPTICS algorithm [6], as one among various approaches for hierarchical densitybased clustering, includes ordering points to identify the clustering structure. OPTICS is based on DBSCAN [36] and the work of Stuetzle [37].
OPTICS computes a Minimum Spanning Tree (MST) of the data, where edge weights represent pairwise distances. These distances are smoothed by a density estimator, called core distance. The core distance of a point is the smallest threshold such that is still considered a core object by the DBSCAN algorithm; that is, has at least objects in its neighborhood within radius . The resulting distance, which is used to construct the MST, is called reachability distance (RD). Taking as input parameter for smoothing the density estimation, the reachability distance of point is defined relative to a reference object as the minimum of the core distance of and the actual distance between and . The outcome of the algorithm can provide us information about the clustering of the objects (see Section 3.4).
3.3. Using Stacked Autoencoders for Data Representation
Densitybased algorithms, traditionally, use the Euclidian distance metric [38]. Such distance metrics are prone to high dimensionality related problems. If we have a feature space of many dimensions, that is, the tracked course of a ship, clustering performance decreases.
Let and be points drawn from a ddimensional Gaussian distribution, so that and . Then their expected distance satisfies [39]Thus, the term , where is a scalar denoting the dimensions of the Gaussian distribution, overshadows the informative term . At this point, the need of robust lowdimension features becomes apparent. In such cases the use of autoencoders is advised [5].
An autoencoder is a neural network that is trained to attempt to copy its input to its output. Internally, it has a hidden layer that describes a code used to represent the input. The network may be viewed as consisting of two parts: an encoder function and a decoder that produces a reconstruction . Autoencoders are designed to be unable to learn to copy perfectly, since they are trained such that instead of . The model often learns useful properties of the data, because it is forced to prioritize which aspects of the input should be copied.
Usually, training the autoencoder to perform the input copying task will result in taking on useful properties, constraining to have smaller dimension than . An autoencoder whose code dimension is less than the input dimension is called undercomplete. Learning an undercomplete representation forces the autoencoder to capture the most salient features of the training data
The learning process is described simply as minimizing a loss function, for example, , where is a loss function penalizing or being dissimilar from , such as the mean squared error. When the decoder is linear and is the mean squared error, an undercomplete autoencoder learns to span the same subspace as PCA. In this case, an autoencoder trained to perform the copying task has learnt the principal subspace of the training data as a side effect
A sparse autoencoder is simply an autoencoder whose training criterion involves a sparsity penalty on the code layer , in addition to the reconstruction error, that is, . Sparse autoencoders are typically used to learn features for another task such as classification. An autoencoder that has been regularized to be sparse must respond to unique statistical features of the dataset it has been trained on, rather than simply acting as an identity function.
The core idea of our work lies in using stacked autoencoders to capture a representation of the main patterns present in the data. By doing so, any outlier in data samples will not be explained well using that representation. In other words, outliers will have significant variations from the rest of the data.
3.4. Identifying Outliers
The outlier detection is a combinatory thresholdbased approach built on the interquartile range rule, as in [40], OPTICS output (see Section 3.2), and AIS/OTH matched data (see Section 3.5).
OPTICS outputs (i.e., reachability distances of the ordered ships) are treated as a continuous signal, over which we identify the peaks. Peaks correspond to significant changes between the closest compared vehicles. As such, anything that varies from the norm has a peak, allowing the easy identification of a possible outlier. Then, we calculate a threshold value defined as , , where denotes the number of ships at a time and is the reachability distances vector, in a descending order.
In case that an outlier provides AIS data, the detection regarding that ship is ignored. At first, for a specific time instance, ships are ordered in a densityreachable way (Figure 2). Points close to each other should belong to the same cluster, unless there is a significant change in RD value. Then, the outlier RD value threshold is defined over 10% of highest RDs.
3.5. Matching OTH Data to AIS
As explained in Section 3.1, AIS data contain, among others, ships’ trajectory points. These coordinates are compared to the radar ones, to identify the similarity among the trajectories. Let us denote as the available discrete time instances, created from the ground radar for ship , . The equivalent case for AIS data is for any ship , , that provides AIS data.
Figure 3 illustrates the available trajectories over a specified area for both radar and AIS data. At this point, we should note that trajectories are calculated for various time intervals, which do not, usually, coincide among the two systems. Typically, for the same ship , in a ratio of four radar time instances to one AIS time instance. Also, note that , so that a 1to1 match among radar and AIS tracked ships is not feasible. Therefore, we should consider both the temporal and the spatial information, to find the matches. The algorithm (presented in pseudocode in Algorithm 1) performs the vessel matching (Figure 4), given OTH and AIS information recorded at the same time (for further details about used data see Section 4.2).

(a)
(b)
The matching process is based on a voting mechanism. For each of the radar tracked ships , at a time instance , , we calculate the closest ships , according to their AIS position at the specific time. In order to identify the corresponding (closest) AIS time instance of ship to radar entry , at a time , we calculate the time difference ; then the corresponding time instance is given as . In case that , AIS entry is not matched to at time . Then, once we have a set of matched ship instances , , we find closest entries to according to their position (i.e., longitude, latitude), so that , .
4. Experimental Results
In the following subsections, we describe the dataset utilized for the experiments, the performance evaluation metrics employed, and the system setup details, before presenting the experimental evaluation of the proposed framework.
4.1. Computational Complexity
Data preprocessing creates a set of OTH data related entries, for a predefined set of past moments, for each one of the tracked ships, at a moment . Since both and are constants defined by the user, the required runtime is . The mapping process of a trained SA is per datum, since SAs are neural networks with a defined number of neurons. OPTICS processes each point once and performs one neighborhood query during this processing. Given a spatial index that grants a neighborhood query in runtime, an overall runtime of is obtained. The matching process between AIS and OTH entries requires an overall runtime of , , since we compare each of the OTH tracked ships to each of the ships equipped with AIS. Table 1 displays the computational complexity of the different processing steps.

4.2. Utilized Dataset
The utilized dataset pertains to approximately 6 hours of data captured from the Mediterranean coast of France by Diginext in July 2016 in the context of the RANGER EU Horizon 2020 project. AIS data for the same period were also obtained for use as ground truth.
A total of 556 ship entries were in this 6hour dataset. The following data provided entries are used:(1)Longitude and latitude: position values provided in degrees. The typical range is and , respectively(2)Course and speed: course is calculated in degrees, typically in the range , and speed in m/s(3)Doppler frequency: it is calculated in Hz, typically in the range (4)Raw Rx azimuth: azimuth angle from the Rx site in the raw spatial grid (equivalent to the reception beam), typically in the range (5)Local noise: noise level in the surrounding of the plot. It is calculated in dBm, in the range (6)Global noise: background noise level of all rangeDoppler map. It is calculated in dBm, in the range .
4.3. Performance Metrics
Formally, a cluster analysis can be described as the partitioning a number of classification objects in groups or clusters , . Given objects , where denotes the th element of . The grouping of all objects , , in clusters can be defined as follows: The above formulation ensures that the association of each object to a cluster is unique. A unique association is a valid case for both hierarchical and partitioning cluster analysis. Given matrix , various internal quality indices have been calculated, to determine an optimal clustering.
4.3.1. Calinski–Harabasz Index
The Calinski–Harabasz index (CHI) [41] is defined according to the following equation:where is defined asand is defined as starts at a comparably large value. With increasing number of clusters , approaching the optimal clustering solution in groups, the value should significantly decrease due to increasing compactness of each cluster. As soon as the optimal solution is exceeded an increase in compactness and thereby a decrease in value might still occur. However, any decrease in value should be notably smaller.
Calculated for each possible cluster solution, the maximum CHI value indicates the best cluster partitioning of the data.
4.3.2. Davies–Bouldin Index
The Davies–Bouldin index (DBI) [42] is an internal evaluation scheme, where the validation of how well the clustering has been done is made using quantities and features inherent to the dataset. DBI is defined as follows:where is defined as is a distance function, defined as , and is defined asAll the above equations assume that .
For each cluster an utmost similar cluster—regarding their intracluster error sum of squares—is searched, leading to . The index then defines the average over these values. In this case, the minimum index value corresponds to the best cluster solution.
4.3.3. Silhouette
The silhouette value is a measure of how similar an object is to its own cluster (cohesion) compared to other clusters (separation). The silhouette ranges from −1 to 1, where a high value indicates that the object is well matched to its own cluster and poorly matched to neighboring clusters. If most objects have a high value, then the clustering configuration is appropriate. If many points have a low or negative value, then the clustering configuration may have too many or too few clusters.
For each datum , let be the average dissimilarity (distance) of with all other data within the same cluster . Let be the lowest average dissimilarity of to any other cluster , , of which is not a member. We now define a silhouette asthus, . Values close to one indicate that the datum is appropriately clustered at . The average silhouette value over all data, that is, , is another measurement for the quality of the generated clusters.
4.4. Experimental Setup
The first step should be the definition of the feature space on which radar data are mapped. As a starting point, we investigated the dimensional space provided by PCA, maintaining 99.1% of the original variation. The adopted stacked autoencoder approach consists of three layers or four layers, depending on the PCA outcome. The loss function was the wellknown mean square error [43] with L2 and sparsity regularizers [44].
Ships track history is composed of 9 consecutive frames, each containing all data as described in Section 4.2. Data are normalized using minmax approach, prior to mapping or clustering approach. The system ignores ships with a narrow appearance span. Any ship that has no enough sufficient entries, that is, 3/4 of past moments tracks, is not taken under consideration.
4.5. Evaluation of Results
OPTICS algorithm outcomes depend on the selection of minimum cluster size. We have investigated the clustering outputs assuming at least 2, 5, 8, 11, 14, 17, 20, 23, and 26 members in each cluster. Clustering over SA mapped data performed better than using raw or PCA mapped data, for most of investigated cases.
According to CHI (Figure 5), highest scores are achieved when using 26 ships per cluster. It is intriguing that cluster performance scores over raw data outperform PCA mapped data scores. There is an increasing trend on the CHI as the minimum cluster size increases. The trend is clearly illustrated for SAs, less for raw data, and slightly for PCA projected data.
The next step was the investigation of DBI scores for the same minimum cluster size setup (Figure 6). This time, the best scores are achieved using 14 or 20 as the cluster size. SA mapping provides better clustering scores in five out of seven investigated cases. Regardless of the mapping method, CHI scores, over SA mapped data, improve as the number of clusters rises, but not in a monotonic way.
The last cluster performance metric was the average silhouette distance (Figure 7). Results suggest that accepting two ships as minimum cluster size is the best possible setup, for PCA mapped data. On the other hand, if we use SA for data mapping, the minimum cluster size should be set as 20.
Another significant performance metric is the average reachability distance itself. The smaller the reachability distance of a point is, the higher the density is around it. The core idea of the proposed approach is that only outliers should vary significantly from the norm, on the projected feature space. Thus, all the ships, minus the outliers, should have similar feature values, which results in reduced reachability distances.
Providing more training data allows SA to adjust the mapping process to the norm. As illustrated in Figure 8 the average reachability distance tends to one, at a slow pace, while increasing the number of training samples. The variance of the RD is, also, reduced when using more time instances for training, as shown in Figure 9. Furthermore, SA mapping allows for the creation of more clusters compared to PCA or raw data clustering (Figure 10).
Regardless of the adopted feature mapping approach, OPTICS outputs are at least four times less in value, compared to calculated RDs using raw data (see Figures 11(a) and 11(b), top). Additionally, SAs result in more clusters, in most of the cases (see Figures 11(a) and 11(b), bottom). Increasing the number of minimum ships per cluster, close objects have almost identical reachability distances, resulting in almost linear subregions, within RD curve.
(a)
(b)
The last step of the performance analysis provides empirical findings. In most of the cases, SAs mapped data results in detection of more outliers compared to the other approaches (Figure 12). The maximum number of detected outliers was three. PCA resulted in no detection at any time.
There was the possibility of unwanted outlier identification. In particular, ships providing AIS data were considered, a few times, possible outliers. Figure 13 illustrates the case. Typically, using SAs resulted in few possible outliers, which however were not accepted as valid detection, as explained in Section 3.4.
5. Conclusions
In our article, a novel approach that identifies unexpected behavior in ship plot and track patterns, as captured by an OTH radar, has been presented. The core idea is the unsupervised development of a mapping process, which can project the raw data in a compact, lower feature space. Outliers projected to the same space should have significantly different values. Stacked autoencoders and PCA were used for the mapping process and compared against the exploitation of raw data, for the identification of unusual ship behavior. Densitybased clustering algorithms (OPTICS) were employed for clusteringbased outlier detection. Experimental results suggest that the approach based on SAs outperforms the other approaches in both generated cluster quality and outliers’ identification.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research leading to these results has received funding from the European Commission’s H2020 Research and Innovation Programme, under Grant Agreement no. 700478 (RANGER project). The authors would like to thank all project partners for their collaboration and especially their partners from Diginext for the provision of the OTH and AIS data and documentation.
References
 L. Sevgi, A. Ponsford, and H. C. Chan, “An integrated maritime surveillance system based on highfrequency surfacewave radars, Part 1. Theoretical background and numerical simulations,” IEEE Antennas and Propagation Magazine, vol. 43, no. 4, pp. 28–43, 2001. View at: Publisher Site  Google Scholar
 S. Grosdidier, A. Baussard, and A. Khenchaf, “HFSW radar model: Simulation and measurement,” IEEE Transactions on Geoscience and Remote Sensing, vol. 48, no. 9, pp. 3539–3549, 2010. View at: Publisher Site  Google Scholar
 D. Nikolio, Z. Popovic, M. Borenovio et al., “Multiradar multitarget tracking algorithm for maritime surveillance at OTH distances,” in Proceedings of the 17th International Radar Symposium, IRS 2016, pp. 1–6, May 2016. View at: Publisher Site  Google Scholar
 G. Vivone, P. Braca, and J. Horstmann, “Knowledgebased multitarget ship tracking for HF surface wave radar systems,” IEEE Transactions on Geoscience and Remote Sensing, vol. 53, no. 7, pp. 3931–3949, 2015. View at: Publisher Site  Google Scholar
 I. Goodfellow, Y. Bengio, and A. Courville, Deep Learning, MIT Press.
 M. Ankerst, M. M. Breunig, H.P. Kriegel, and J. Sander, “OPTICS: ordering points to identify the clustering structure,” in Proceedings of the ACM SIGMOD International Conference on Management of Data (SIGMOD '99), pp. 49–60, New York, NY, USA, 1999. View at: Publisher Site  Google Scholar
 G. Kouemou and F. Opitz, “Impact of Wavelet based signal processing methods in radar classification systems using Hidden Markov Models,” in Proceedings of the 2008 International Radar Symposium, IRS, pp. 1–4, May 2008. View at: Publisher Site  Google Scholar
 M. GarbanzoSalas and W. K. Hocking, “Spectral analysis comparisons of Fouriertheorybased methods and minimum variance (Capon) methods,” Journal of Atmospheric and SolarTerrestrial Physics, vol. 132, pp. 92–100, 2015. View at: Publisher Site  Google Scholar
 P.H. Jau, Z.M. Tsai, N.C. Kuo et al., “Signal processing for harmonic pulse radar based on spread spectrum technology,” IET Radar, Sonar & Navigation, vol. 8, no. 3, pp. 242–250, 2014. View at: Publisher Site  Google Scholar
 L. Liu, D. Fu, Y. Zhu, D. Su, and M. Diao, “Applications of Bootstrap in Radar Signal Processing,” in Intelligent Computation in Big Data Era, vol. 503, pp. 474–479, 2015. View at: Publisher Site  Google Scholar
 D. Xu, L. Du, H. Liu et al., “Compressive sensing of steppedfrequency radar based on transfer learning,” IEEE Transactions on Signal Processing, vol. 63, no. 12, pp. 3076–3087, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 B. A. Massinas, A. Doulamis, N. Doulamis, and D. Paradissis, “An echo state network for ionospheric disturbances behavior modeling on Spaceborne Interferometric Synthetic Aperture Radar,” in Proceedings of the AIAA SPACE 2013 Conference and Exposition, September 2013. View at: Google Scholar
 B. A. Massinas, A. Doulamis, N. Doulamis, and D. Paradissis, “Applied optimal estimation for ionospheric disturbances behavior on spaceborne interferometric synthetic aperture radar systems,” in Proceedings of the AIAA SPACE Conference and Exposition, 2015, American Institute of Aeronautics and Astronautics, September 2015. View at: Google Scholar
 B. A. Massinas, A. Doulamis, N. Doulamis, and D. Paradissis, “Ensemble classifiers in optimal estimation for ionospheric disturbances behavior on spaceborne interferometric SAR systems,” in Proceedings of the AIAA Space and Astronautics Forum and Exposition, SPACE 2016, September 2016. View at: Google Scholar
 M. Cao, S. Li, R. Wang, and N. Li, “Interferometric phase denoising by median patchbased locally optimal wiener filter,” IEEE Geoscience and Remote Sensing Letters, vol. 12, no. 8, pp. 1730–1734, 2015. View at: Publisher Site  Google Scholar
 X. Meng, Z. He, G. Feng, and B. Xiao, “An improved wavelet denoising algorithm for wideband radar targets detection,” Circuits, Systems and Signal Processing, vol. 32, no. 4, pp. 2003–2026, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 R. Ahmed, N. Maheshwari, and P. Lalla, “Wavelet based iterative thresholding for denoising of remotely sensed optical and synthetic aperture radar images,” in Proceedings of the IEEE International Conference on Advanced Communication, Control and Computing Technologies (ICACCCT '14), pp. 1331–1335, May 2014. View at: Publisher Site  Google Scholar
 M. Astry, R. Narayanan, and M. Rangaswamy, “Sparsitybased signal processing for noise radar imaging,” IEEE Transactions on Aerospace and Electronic Systems, vol. 51, no. 1, pp. 314–325, 2015. View at: Publisher Site  Google Scholar
 M. Gong, J. Zhao, J. Liu, Q. Miao, and L. Jiao, “Change detection in synthetic aperture radar images based on deep neural networks,” IEEE Transactions on Neural Networks and Learning Systems, vol. 27, no. 1, pp. 125–138, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 J. Liu, M. Gong, J. Zhao, H. Li, and L. Jiao, “Difference representation learning using stacked restricted Boltzmann machines for change detection in SAR images,” Soft Computing, vol. 20, no. 12, pp. 4645–4657, 2016. View at: Publisher Site  Google Scholar
 H. Xie, S. Wang, K. Liu, S. Lin, and B. Hou, “Multilayer feature learning for polarimetric synthetic radar data classification,” in Proceedings of the Joint 2014 IEEE International Geoscience and Remote Sensing Symposium, IGARSS 2014 and the 35th Canadian Symposium on Remote Sensing, CSRS 2014, pp. 2818–2821, July 2014. View at: Publisher Site  Google Scholar
 S. Maresca, P. Braca, J. Horstmann, and R. Grasso, “A network of HF surface wave radars for maritime surveillance: Preliminary results in the German Bight,” in Proceedings of the 2014 IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2014, pp. 6077–6081, May 2014. View at: Publisher Site  Google Scholar
 S. Maresca, P. Braca, J. Horstmann, and R. Grasso, “Maritime surveillance using multiple highfrequency surfacewave radars,” IEEE Transactions on Geoscience and Remote Sensing, vol. 52, no. 8, pp. 5056–5071, 2014. View at: Publisher Site  Google Scholar
 P. Braca, S. Maresca, R. Grasso, K. Bryan, and J. Horstmann, “Maritime surveillance with multiple overthehorizon HFSW radars: An overview of recent experimentation,” IEEE Aerospace and Electronic Systems Magazine, vol. 30, no. 12, pp. 4–18, 2015. View at: Publisher Site  Google Scholar
 G. E. Hinton and R. R. Salakhutdinov, “Reducing the dimensionality of data with neural networks,” Science, vol. 313, no. 5786, pp. 504–507, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Guo, Y. Liu, A. Oerlemans, S. Lao, S. Wu, and M. S. Lew, “Deep learning for visual understanding: a review,” Neurocomputing, vol. 187, pp. 27–48, 2016. View at: Publisher Site  Google Scholar
 K. Makantasis, A. Doulamis, N. Doulamis, and K. Psychas, “Deep learning based human behavior recognition in industrial workflows,” in Proceedings of the 23rd IEEE International Conference on Image Processing, ICIP 2016, pp. 1609–1613, September 2016. View at: Publisher Site  Google Scholar
 N. Doulamis and A. Voulodimos, “FASTMDL: Fast Adaptive Supervised Training of multilayered deep learning models for consistent object tracking and classification,” in Proceedings of the 2016 IEEE International Conference on Imaging Systems and Techniques, IST 2016, pp. 318–323, October 2016. View at: Publisher Site  Google Scholar
 G. Hinton, L. Deng, D. Yu et al., “Deep neural networks for acoustic modeling in speech recognition: the shared views of four research groups,” IEEE Signal Processing Magazine, vol. 29, no. 6, pp. 82–97, 2012. View at: Publisher Site  Google Scholar
 E. Mason, B. Yonel, and B. Yazici, “Deep learning for radar,” in Proceedings of the 2017 IEEE Radar Conference (RadarConf17), pp. 1703–1708, Seattle, WA, USA, May 2017. View at: Publisher Site  Google Scholar
 S. Deng, L. Du, C. Li, J. Ding, and H. Liu, “SAR automatic target recognition based on euclidean distance restricted autoencoder,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, no. 99, pp. 1–11, 2017. View at: Publisher Site  Google Scholar
 G. Vivone, L. M. Millefiori, P. Braca, and P. Willett, “Performance assessment of vessel dynamic models for longterm prediction using heterogeneous data,” IEEE Transactions on Geoscience and Remote Sensing, no. 99, pp. 1–14, 2017. View at: Publisher Site  Google Scholar
 “HF Surface Wave Radar,” [Accessed: 16Jun2017] Available: http://www.diginext.fr/en/offer/criticaloperationsupportsystems/hfsurfacewaveradar. View at: Google Scholar
 “Automatic Identification Systems (AIS),” [Accessed: 04Jun2017] Available: http://www.imo.org/en/OurWork/safety/navigation/pages/ais.aspx. View at: Google Scholar
 H.P. Kriegel, P. Kröger, J. Sander, and A. Zimek, “Densitybased clustering,” Wiley Interdisciplinary Reviews: Data Mining and Knowledge Discovery, vol. 1, no. 3, pp. 231–240, 2011. View at: Publisher Site  Google Scholar
 M. Ester, H.P. Kriegel, J. Sander, and X. Xu, “A densitybased algorithm for discovering clusters in large spatial databases with noise,” in Proceedings of the 2nd International Conference on Knowledge Discovery and Data Mining (KDD '96), pp. 226–231, 1996. View at: Google Scholar
 W. Stuetzle, “Estimating the cluster type of a density by analyzing the minimal spanning tree of a sample,” Journal of Classification, vol. 20, no. 1, pp. 25–47, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 E. Protopapadakis and A. Doulamis, “Semisupervised image metafiltering using relevance feedback in cultural heritage applications,” International Journal of Heritage in the Digital Era, vol. 3, no. 4, pp. 613–627, 2014. View at: Publisher Site  Google Scholar
 M. Hein, J.Y. Audibert, and U. von Luxburg, “From graphs to manifolds – weak and strong pointwise consistency of graph Laplacians,” in Learning Theory, P. Auer and R. Meir, Eds., vol. 3559, pp. 470–485, Springer, Berlin, Germany, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 F. Xiao and C. Fan, “Data mining in building automation system for improving building operational performance,” Energy and Buildings, vol. 75, pp. 109–118, 2014. View at: Publisher Site  Google Scholar
 T. Caliński and J. Harabasz, “A dendrite method for cluster analysis,” Communications in Statistics, vol. 3, no. 1, pp. 1–27, 1974. View at: Google Scholar  MathSciNet
 D. L. Davies and D. W. Bouldin, “A cluster separation measure,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. PAMI1, no. 2, pp. 224–227, 1978. View at: Publisher Site  Google Scholar
 E. Protopapadakis, M. Schauer, E. Pierri et al., “A genetically optimized neural classifier applied to numerical pile integrity tests considering concrete piles,” Computers & Structures, vol. 162, pp. 68–79, 2016. View at: Publisher Site  Google Scholar
 D. W. Winters, B. D. Van Veen, and S. . Hagness, “A sparsity regularization approach to the electromagnetic inverse scattering problem,” Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation, vol. 58, no. 1, pp. 145–154, 2010. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Eftychios Protopapadakis et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.