1. Introduction
Efficient ship detection constitutes an important element of marine surveillance since it improves the supervision of activities such as illegal fishing, oil spillage, marine traffic management, smuggling and drug trafficking [
1,
2]. Standard ship detectors use Synthetic-Aperture Radar (SAR) imagery to find vessels by identifying the energy backscattered from the ship (Radar Cross Section, RCS) in the image [
3,
4]. In some cases however such a method fails due to the fact that the ships may not be clearly visible. The technique of direct vessel identification has several limitations that depend on the ship and radar characteristics, the sea state conditions and the image quality and resolution. Moreover, the ship’s signature may appear displaced with respect to the actual position of the ship, this phenomenon occurs when the sailing direction presents a component perpendicular to the radar’s direction of flight [
5].
Furthermore in some cases ships may be made of non-conducting materials such as fiberglass or wood which are transparent to microwaves leading to a poor or null signature in a radar image. Examples of fiberglass ships are the so-called
go-fast ships.
Go-fast ships are small, between 10 and 12 m, and very light. These characteristics allow them to reach very high speeds (up to 80 knots) which create characteristic wakes [
6]. The transparency of such kind of ships to radars has made them the preferred mean for illegal activities [
7].
An alternative methodology for ship detection is the identification of slicks on the ocean surface due to the release of engine or fish oils. However, specific conditions must occur for these ship-related slicks to be visible, in particular there must be a sufficient amount of oil spilled and moderate wind conditions. Finally, a ship may be identified through the detection of its wake. This paper describes and validates a new methodological approach to detecting ships by exclusively identifying their wakes.
Ship detection is a vastly researched topic due to the wide range of possible applications, from maritime security to traffic planning and illegal activity monitoring such as drug trafficking and illegal fishing. The studies of possible ship detection algorithms use a wide spectrum of radar imagery types, such as optical, infrared and SAR [
8]. References Elvidge et al. [
9] and Lebona et al. [
10] present methods for ship detection through optical imagery. Optical imagery’s advantages include the possibility of vessel identification and feature extraction regardless of their size, shape and material composition thanks to the high spatial resolution provided by optical sensors. The performance of these sensors is however limited due to the dependence on light to illuminate the scene under observation precluding the possibility of analysing the scene at night. Moreover, optical sensors cannot see through clouds therefore the weather conditions must be optimal, i.e., no clouds must cover the scene [
8].
Thermal infrared spectrum has also been studied for the application of vessel detection since such sensors do not depend on solar illumination. These sensors detect the emissions from the objects in the imaged area allowing the analysis of the area at nighttime. Thermal infrared sensors are however still affected by the weather not allowing visibility with the presence of clouds and of high moisture [
8]. References Mirghasemi et al. [
11] and Wu et al. [
12] amongst others, studied the use of thermal infrared imagery for vessel detection.
SAR imagery allows for the imaging of wide areas at relatively high spatial resolution, in all weather conditions and at all times. It presents however high levels of noise, moreover it does not allow for ship type identification and is more highly prone to false detections [
8]. The use of SAR imagery for vessel detection has a vast body of research, using single, double and fully polarimetric products. The approach of target detection using SAR imagery is varied. Papers such as Kouroupis et al. [
13], Allard et al. [
14] and Ringrose et al. [
15] use coherent decompositions of fully polarimetric SAR scattering matrices, such as Cameron’s coherent target decomposition (CTD) and the Odd-Even basis decomposition or use spectral analysis, such as the subaperture coherence, ref. [
14], to detect the targets.
When using single polarization SAR data, the most common method of ship target detection consists in determining a threshold that allows distinguishing between the visible target and the sea. The most common of such methods is the constant false alarm rate (CFAR) method. The CFAR method consists in creating an adaptive threshold within a certain area and analysing the
pdf of such area [
16]. In particular, three windows are selected, one with a size that allows containing a ship, depending on the image’s resolution, a buffer window and a background window that contains the first two. These windows move along the image analysing the
pdf within the windows and for each scan, a threshold is applied [
17]. The use of the CFAR method has also been used as a starting point to create algorithms that would take into account more factors. Leng et al. [
18] propose an improved CFAR algorithm that takes into account not only the pixels’ intensity but also their spatial distribution, allowing to achieve better overall results. Deng et al. [
19] have used the CFAR in combination with deep learning which increased the detection speed and its recognition rate. Other algorithms that use single polarization SAR imagery for target detection are also developed [
20]. Arnaud et al. [
21] uses SAR interferometry and coherency techniques, Liu et al. [
16] uses information from polarimetric SAR data in combination with the CFAR. More recently, the GLRT ship detector [
22] has gained attention for its better performance with respect to the state of the art, though the inclusion and consideration of a scattering model for the ship leads to an increase of the processing time; on the other hand Yang et al. [
23] have introduced a GPU based CFAR method that has significantly improved the computational time with respect to standard CFAR detectors.
Finally, deep learning models are used on SAR imagery for ship detection by adapting object detectors to the situation at hand [
16,
24]. The main hurdles of this method are, firstly, the presence of geometric distortions that generally present themselves in SAR imagery in the forms of shadowing, foreshortening and layover, which cause a high number of false positives along the coastal areas, consequently, sea-land segmentation must always be performed, and secondly, the lack of large, diversified and labelled datasets, necessary for these types of models to produce high quality performances [
2].
Overall the majority of the studies overviewed, whether using double or fully polarimetric SAR products, detect vessels solely through the identification of their signature. Hence the majority of the research focuses on the detection of ships with a clearly visible RCS.
The field of SAR imagery usage for vessel detection through wake identification is instead a less studied one. References [
4,
25,
26] are the few studies presenting an algorithm for wake detection through SAR imagery. The algorithm presented by M.D. Graziano et al. [
25] presents a novel approach for the wakes’ pattern recognition within an image section. The algorithm scans the section pixel by pixel and assigns a confidence index of wake presence in such pixel. The confidence index increases if the neighboring pixels’ index is also high as it increases the possibility of the image indeed containing a wake. The Radon transform is used to convert straight lines in the image domain into dots in a transformed domain. The detection algorithm was applied to SAR imagery from TerraSAR-X mission, which is characterised by geometric resolutions of 1.2 m in slant range and 6.6 m in azimuth. The data used was acquired in VV polarisation and at varying incidence angles. Only the turbulent wake and the narrow-V wake’s components were looked for. The presence of a wake is confirmed if both the narrow-V and the turbulent wake are imaged, therefore Ref. [
27] proposed an index that took into account the contribution of both wakes to the final identification of the wake. The results presented are at the preliminary stage, in fact results regarding only one image containing three wakes are presented. While within these results only one out of the three wakes is identified, the proposed method is very promising and easily applicable to methods such as ML algorithms and CNNs [
27].
The Radon transform however faces a problem when there is some speckle noise affecting the image since it uses the collinearity of points without analysing how these points are distributed. Therefore G. Yang et al. [
4] used the Parallel Coordinate Transform (PCT) for noisy images. Moreover in Ref. [
4] the original image was not binarised but was instead left in grayscale in order to lose as little information as possible. These elements have led to an improvement in performance. Reference [
4] however applies as a first step the detection of the ship’s signature and only then it proceeds to identify the wake, making the method not useful for
go-fast ship detections.
Patents regarding the subject are also present. Yang L. et al. [
28] propose an algorithm in which the identification of the vessel’s signature followed the application of a series of image processing functions and the subsequent application of the Radon Transform. A few older studies and patents on the topic can be found in the literature [
29,
30,
31]. Apart from Graziano et al. [
27], in general the algorithms proposed consist of a three steps approach in which the first step is the detection of the ship’s signature, the second step looks into the processing of the image through various techniques, and finally the third step consists of the detection of the lines composing the wakes. Finally, Karakus et al. [
26] propose an inverse approach to the wake detection problem by enhancing the information in the Radon transform and improving the linear features. The solution involves a Bayesian methodology which leads to a maximum a posteriori (MAP) estimator and generalised minimax concave (GMC) penalty of Selesnick for the cost function. This approach provides promising results while simultaneously avoiding ad hoc thresholding for the enhancement of the image.
The following paper seeks to contribute to the field of SAR ship detection by tackling the problem of the identification of ship wakes through a novel approach consisting of image processing and computer vision techniques applied to C-band imagery. The paper is organized as follows: in
Section 2 the theory behind the wake detection algorithm is explained and the dataset used for the development and testing of the algorithm is presented.
Section 3 describes the methodology of the algorithm proposed. In
Section 4 and
Section 5 the results are, respectively, presented and discussed. In
Section 6, the ship’s speed and heading results are validated and certain particular cases are presented. Finally in
Section 7 conclusions are drawn.
2. Theory & Dataset
Ships of small-to-medium sized dimensions, including
go-fast ships, create a wake that is a combination of two V-shaped wakes and a turbulent wake. The latter appears as a dark line behind the ship that follows its direction of movement [
5]. The two V-shaped wakes are characterised by different possible angle values; one of the V-shaped wakes has a varying half-angle, i.e., angle between one of its arms and the turbulent wake, that can reach a maximum value of ∼19.45°. The second V-shaped wake is characterized by a varying half-angle that can reach a maximum value of 3°. In the paper the wake with the wider angle will be referred to as a V-shaped wake, while that with the narrowest angle will be referred to as a narrow-V wake. The two V-shaped wake components, i.e., their arms, usually appear brighter compared to the background. In most cases only one of the arms of each of the two V-shaped wakes and/or the turbulent wake are visible. Moreover, V-shaped wakes are hardly detectable when the ship is travelling at speeds below 15 knots [
3,
5] when imaging the area with Sentinel-1 imagery (the products used in this study) due to the characteristics of the sensor, i.e., the frequency band used and its spatial resolution.
Figure 1 portrays the various components of the wake just described.
The V-shaped wake’s angle range depends on the ship’s characteristics, in particular its travelling speed and its length. Darmon et al. [
32] describe the link between the half-angle of the V-shaped wake and the
Froude number. The Froude number is defined as:
where
U is the relative flow velocity between the sea and the vessel,
g is the acceleration of gravity and
L is the length of the vessel at waterline level. References [
32,
33] demonstrate that for values of the Froude number bigger than 0.7 the V-shaped wake’s half-angle (
) can be approximated by:
Due to their small size, light weight and high speed capabilities
go-fast ships tend to fall within the higher spectrum of the Froude number and as such, their V-shaped wakes’ half angle can be expected to be approximated with Equation (
2). Relationships (
1) and (
2) are used in the following sections to retrieve the relative velocity of the identified ship when the V-shaped half angle leads to a Froude number greater than 0.7.
Dataset
Sentinel-1 is a mission carried out by the European Space Agency (ESA) within the Copernicus program. The mission images the Earth in 6 days. The sensors operate in C-band [
34] and offer data both in single and double polarization mode. The whole dataset used in this research consists of data acquired in Interferometric Wide Swath (IWS) mode and HH polarisation, processed at GRDH level [
35]. Part of the dataset is used for the development of the algorithm, part for its testing and validation.
The dataset used for the development comprises a total of 13 images, 10 representing non-go-fast ships and 3 representing go-fast ships. In the context of this paper, for go-fast ships the authors mean ships visible in the image only by means of their wakes.
The testing dataset comprises a total of 43 images, 28 representing non-
go-fast ships and 15 representing
go-fast ships. Out of the testing dataset images, five were used for the speed and heading validation later described. The data was acquired between July 2018 and May 2019 in the area of the Gulf of Naples (Italy), Singapore and Genoa (Italy).
Figure 2 shows an example of a non-
go-fast ship and a
go-fast ship signatures in SAR images taken from the dataset. For the validation of the results, two preparatory tasks were performed. Firstly, to determine the visibility of the ship, a standard ship detector available from the SNAP toolbox [
36] was applied to each image in the dataset. The parameters used were:
The ships not detected, hence visible only through the presence of their wake, were classified as
go-fast, while those detected were classified as non-
go-fast. Secondly, for the validation of the ships’ speed and direction of movement, data from the Marine Traffic Portal was used. The Marine Traffic Portal provides the AIS (Automatic Identification System) data that each ship passing within a specific geographical area of choice provides to the local authorities The data in the portal is updated every two seconds to three minutes depending on the vessel’s type and moving status [
38].
Table 1 and
Table 2 present the complete dataset, distinguishing between the elements used during the development and those used during the testing and validation phases.
3. Method Description
The scope of the algorithm is to efficiently detect small to medium sized ships and
go-fast ships exclusively through the identification of their wakes. Moreover, the algorithm must evaluate the RCS of the ship, if present, in order to distinguish between
go-fast and non-
go-fast ships. Finally, it must evaluate the backscattering of the wake in order to determine the position, speed and direction of movement of the vessel.
Figure 3 presents the flow chart of the wake detector algorithm developed. As described in
Section 2 the dataset used consisted in Sentinel-1 C-band imagery, however, the algorithm can be used with any type of SAR imagery. After a preliminary image processing phase, the algorithm is based on two main functions described in the following, a line detection and an error check function. Depending on whether a ship was found or not it finally determines the ship’s characteristics, namely its type, speed and direction of movement.The input image is divided into smaller sections each independently analysed. Each section goes through up to three steps depending on whether a wake is found or not. The dimension of the section was chosen based on the average dimension of a wake, and was such that, when present, it would contain the wake.
3.1. Image Processing
The main focus of the image processing step is to emphasize the particular features withing the image’s section under analysis whilst decreasing the effects of the sensor’s noise. This is done in order to highlight a possible wake within the section under analysis. The step includes a combination of global and adaptive histogram equalization and enhancement. Overall this first step of the algorithm is divided in four sub-steps:
Brightness Equalization
Despeckling using a Gaussian filter
Features enhancement since wakes appear as a particularity in an assumed uniform background, i.e., the sea. This step comprises the execution of three sub-steps:
- (a)
Adaptive Histogram Equalization: The function used is the Contrast Limited Adaptive Histogram Equalization. In this function contrast limiting is applied, i.e., if any histogram bin is above a certain contrast limit (defined by the user), those pixels are clipped and distributed uniformly to other bins before applying histogram equalization [
39]. The tile dimension has been set to
pixels and the contrast limit to 20. Such values were chosen after a period of testing on the testing dataset.
- (b)
Impose to pixel a non-linearly proportional value: Since an important feature will be characterized by a value of pixel different with respect to the background, this difference is wanted to be emphasized. The function used was customed for the specific problem at hand, each pixel value is used as an exponent to a base-two power. The pixel value is first divided by 31 in order to have a resulting value within the range (0–255).
- (c)
Global Histogram Equalization
Image Binarization
3.2. Line Detection
The Hough Transform [
39] is used for the detection of the lines. The wakes’ arms are expected to be brighter in the image compared to the background, while the turbulent wake is expected to be darker. Hence the algorithm looks for two bright lines and a dark line. The transform accepts a binary image in input. Two binary images are computed, one for the identification of bright lines (emphasizing the brightest pixels in the image) and another for the identification of dark lines (emphasizing the darkest pixels in the image).
For each line found, the Hough transform function returns the characteristics of the line and a score as output. If more than one line is found the line with the highest score is considered one of the arms of the wake. The following lines are then compared with such line in terms of angle aperture, i.e., the second arm must create an angle in the range 8°–40° (40°
°, where 19.45° is the maximum possible half-angle of a wake). The angle of 8° on the lower end of the spectrum leads to a half-angle of 4°. As explained in
Section 2 a narrow-V wake’s angle creates a maximum half-angle of 3° therefore an angle greater than 3° was chosen in order to avoid the possible identification of the two narrow-V wakes arms as that would lead to a wrong result in the ship’s speed determination. Moreover, two lines are considered components of a wake if they are within a maximum set distance from one another. On the other hand, if only one bright line is found the algorithm will immediately proceed with looking for a dark line. Depending on the number of bright lines detected (one or two), the line detection algorithm will then proceed differently in search of the turbulent wake.
3.2.1. One Bright Line Found
The bright line is assumed to be the arm of the V-shaped wake, hence the algorithm looks for dark lines in the angle range of 6° to 20°. If no turbulent wake is found, no wake is determined to be found.
3.2.2. Two Bright Lines Found
The possibility of one of them being the arm of the narrow-V wake must be investigated. The turbulent wake is looked for in the proximity of each of the bright lines found in the angle range 0°–3°. If a dark line is found, such a line is classified as the turbulent wake; the bright line detected close to it is classified as the narrow-V wake’s arm, while the other is classified as the arm of the V-shaped wake. In case no turbulent wake is found, the two bright lines are classified differently depending on the angle they create. If such angle is in the range of 8°–12°, then the two lines are considered an arm of the narrow-V wake and one of the V-shaped wake, otherwise, the two lines found are considered the two arms of the V-shaped wake. The angle range was chosen based on the typical structure of a wake witnessed during the initial training phase. In the first case two possible turbulent wake’s positions are determined hence two headings of the ship are presented in output.
3.3. Error Check Function
The sea generally represents a uniform background to the wake. However, due to environmental conditions occurring during the acquisition of the image, this may not always be the case. To avoid false identifications due to sea artifacts, the Error Check function is implemented. If bright lines are found in the image, before looking for dark lines the algorithm checks that the lines detected are indeed visible lines in the original non-binary image. Analogously this is done for the detection of dark lines.
For each detected line, its pixels are compared with the surrounding area. The comparison is performed by calculating the absolute value of the difference between the average value of the pixels composing the line and that of the pixels of a window of the surrounding area. If such value is smaller than a set threshold determined during the training phase then the line is deemed as being a sea artifact otherwise it is confirmed as a line. The surrounding area was chosen in order to avoid the pixels in the closest proximity of the detected line since these could also constitute the wake, but also in order to avoid considering an area too far from the detected line, since a change in the sea texture and luminosity could cause an unwanted line confirmation. The surrounding window to analyse was defined as the area far from the detected line in the range of 4 to 8 pixels.
If no wake is found and part of the input image must still be analysed, the window is scrolled of a certain amount of resolution cells defined by considering one order of magnitude smaller than the dimension of the window, and the process is repeated from the first step on the new region. If a wake is found the algorithm proceeds with determining the ship’s position and type.
3.4. Ship Type Determination
Based on the wake’s shape the ship’s expected position is determined as the intersection point between the lines composing the wake. An area around such point is then searched in order to find any ship signature, the presence or absence of which would allow distinguishing between go-fast and non-go-fast ships.
3.5. Ship Direction Determination
The direction of the ship is characterised by two main components: the direction of movement i.e., North, South, East or West and the orientation i.e., the angle with respect to the North.The ship’s direction is determined by comparing the position of the turbulent wake with respect to the estimated position of the ship (not the position of its signature), taking in consideration that the turbulent wake is found behind the travelling vessel. On the other hand, the ship’s orientation is determined by getting the inclination of the turbulent wake since, by definition of turbulent wake, this appears as a straight line behind the travelling vessel. If the turbulent wake was not detected its expected orientation is computed.
When the step is concluded, if part of the image still has to be analysed, the window is scrolled again and the process is repeated on the new section. Once the entire image’s analysis is completed, the algorithm provides the position of the ships found, together with their type, direction and speed based on the characteristics of their related wake. For every ship/wake found the ship speed and direction of movement are evaluated based on Equations (
1) and (
2).
7. Conclusions
In this paper, a novel algorithm able to detect ships exclusively from their wakes was presented. The application of a wake detector is of particular interest in situations where there is a good likelihood that wakes are visible even if the corresponding vessels are not. Contrary to similar previous state-of-the-art, the images analysed covered a widely extended area of ocean not allowing for a fair comparison to be made amongst different techniques. A total of 56 images portraying a wake and in thirty-eight of the cases the ship’s projection, taken from Sentinel-1 products, were used for the development and testing of the algorithm presented. The algorithm resulted in a 93% success rate, meaning that forty out of the forty-three ships composing the testing dataset were correctly identified. A wake is seen as composed of five elements, the arms of the three sub-wakes that compose a wake and that are imaged as either bright or dark lines in the sea. Out of these five elements, the algorithm looks for three, namely one of the arms of each sub-wake, and deems an object in the image to be a wake if at least two out of the three lines are detected at particular angle positions and below a set distance between one another. Out of the three wakes not found at least one of the arms was always found. In one case the wake corresponded to that of a medium to large sized ship, which caused the arms to be not well defined. Such ships and therefore wakes are not important for the scope of the study since they always present a very visible signature, however, this means that this algorithm cannot be considered a universal ship wake detector since certain indexes would have to be adjusted for this scenario.
Some further limitations in the algorithm have been discussed. The problem of the presence of ship signature artifacts and how they may cause problems in the detection of the wakes’ components has been presented. While this did not affect the detection of the ship itself, it did affect the identification of the wake’s components and hence of the ship’s heading and speed. It was then presented how the removal of the artifact from the image before the application of the wake detector eliminated the components’ false detection.
Based on the limitations just discussed, a future step could be to consider adopting the wake detector in combination with another ship detector which masks the ships’ signatures once detected. Considering the type of ships whose wakes were not found, this additional step may also help improve the algorithm’s overall ship detection performance and create a more universal ship detector.
A validation for the ship’s speed and heading determination was performed. Overall the ship’s velocity calculation resulted in a ∼3 knots error and the heading in a ∼4° error. While very limited in terms of the number of images used, this validation allowed to substantiate the methods used to assess such values.
The paper also presents the problem that sea artifacts may bring in causing False Positives. While the assumption of considering the sea surface uniform is valid, this is valid only in limited circumstances such as in open sea scenarios and in situations with low wind conditions as discussed also by Reference [
25]. Sea artifacts were found in a dataset composed of only thirteen sea images and it is therefore believed that they should not be considered as a rare event but instead be included in any algorithm dealing with wake detection through line identifications. Finally, the algorithm presented uses low computationally expensive techniques allowing for potential real-time analysis of the area under investigation.
The results presented pave the way for further studies relating the application of the Hough transform to wake detectors and have identified soft spots that should be further addressed with any type of wake detector involving line identifications. Since the wake detector has been tested on Sentinel-1 SAR images only, further validation on additional SAR datasets from different missions or with different settings of the radar is required. In particular, it would be of great advantage using radar imagery with higher spacial resolution since the elements of the wake would appear more clearly and therefore facilitate the identification.