Wake Component Detection in X-Band SAR Images for Ship Heading and Velocity Estimation

A new algorithm for ship wake detection is developed with the aim of ship heading and velocity estimation. It exploits the Radon transform and utilizes merit indexes in the intensity domain to validate the detected linear features as real components of the ship wake. Finally, ship velocity is estimated by state-of-the-art techniques of azimuth shift and Kelvin arm wavelength. The algorithm is applied to 13 X-band SAR images from the TerraSAR-X and COSMO/SkyMed missions with different polarization and incidence angles. Results show that the vast majority of wake features are correctly detected and validated also in critical situations, i.e., when multiple wake appearances or dark areas not related to wake features are imaged. The ship route estimations are validated with truth-at-sea in seven cases. Finally, it is also verified that the algorithm does not detect wakes in the surroundings of 10 ships without wake appearances.


Introduction
Maritime traffic monitoring is gaining great advantages from spaceborne SAR imagery.Several international projects are investigating an operational ship monitoring system [1][2][3][4], and many experimental results are available on ship detection in SAR images [5][6][7][8][9][10].In this ambit, ship wake detection not only has strong interest as the validation of ship presence, but also generates data useful for ship route estimation.In fact, if wake direction provides information about ship heading, methodologies have also been developed to derive ship velocity from additional wake-based information [11][12][13][14][15][16].
The typical ship-wake appearance in high-resolution SAR images exhibits the assessed structures detailed in [17].A central dark line is the signature of the turbulent wake, which is aligned with the ship longitudinal axis.Two bright linear features (narrow-V wakes) can be seen within a half-angle ranging from 1.5 ˝to 4 ˝with respect to the turbulent dark line.The most external couple of wakes are representative of the Kelvin arms, which develop in an angle with an aperture of ˘19.5 ˝with respect to the turbulent dark line.Finally, it is rarely possible to observe the alternating dark and bright wave pattern roughly parallel to the ship track, called internal wave wakes [17].However, as it is more apparent in L-band rather than in X-band or in C-band for the longer L-band Bragg waves' persistence [17], it is not included in the proposed wake detection process.
Since wakes appear in SAR images as dark or bright straight lines, the existing algorithms for wake detection exploit the classical linear feature detection methods, such as Radon [18] or Hough [19] transforms, which convert bright (dark) lines in the input domain into peaks (troughs) in the transformed domain.It has been shown that the Hough transform is unsuited for wake detection due to a too heavy preprocessing phase [20], and the Radon transform utilization [21] is encouraged by the low computational cost of the inverse transform, which can be performed by FFT.Some authors [15,22,23] proposed techniques to enhance the peaks of the wake features and applied thresholding techniques to detect such peaks.Tunaley [22] identifies the portion of Radon transform interested by the wake feature under analysis as the area in the Radon domain comprised within two sinusoidal curves.Rey et al. [23] suggest to manually size the subsets of the SAR image accordingly to the wake length.Zilmann et al. [15] identify a region in the Radon domain where the transform is characterized by constant grey, i.e., intermediate level, and its mean value is used to detect troughs and peaks whose values are respectively below and above the mean value of the background domain.
With reference to ship route estimation, the most used and assessed technique is based on the evaluation of the azimuth displacement between ship and ship-wake vertex [11][12][13].A conventional SAR processor, in fact, erroneously focuses on moving targets, focused with an azimuth offset with respect to their true position, which depends on the target-radar radial velocity [12].Results on 250 wakes extracted from medium to low resolution SAR data (Radarsat-1, ERS, Seasat) [14] show that the accuracy of the computed velocity is within about ˘15 percent of the true speed for most of the images, provided that the wakes exhibit a clear appearance and are not closely aligned with the azimuth direction.A more recent technique is based on the analysis of the Kelvin arms, which exhibit a dark/bright pattern whose wavelength depends on the gravity acceleration and on the ship velocity [15,16].
In this paper, a wake detection algorithm based on the Radon transform is presented.Differently from the literature, the bright/dark dots in the transformed domain are detected considering their angular distance as foreseen by hydrodynamics [17].In addition, they are validated as representative of the ship wake components by means of novel merit indexes in the intensity domain.State-of-the-art techniques for heading and velocity estimation are also incorporated in the processing.Algorithm validation is performed on 13 X-band SAR images from the TerraSAR-X and COSMO/SkyMed missions, with different polarization and incidence angles.

Algorithm Overview and Objectives
The algorithm herein presented is aimed at estimating ship heading and velocity.In addition, solutions reducing the efforts of human operators are pursued, in view of fully-automatic procedures for future inclusion in a potentially operational environment.In this respect, the presented algorithm is designed to look for a wake structure, always formed by: one central dark line (turbulent wake), two pairs of bright linear features (narrow-V wake and Kelvin arms), with an azimuth offset between ship and wake vertex.Since it is well known that wake appearance can, in practice, have less features imaged than the above five ones, a successive analysis is devoted to the confirmation (rejection) of the identified linear structures as real (false) components of the ship wake.
The proposed algorithm introduces evolutions with respect the state-of-the-art in both the detection and the confirmation phases.In particular, detection is performed considering the angular relations among the components of the wake structures as foreseen by hydrodynamics theory [17].The confirmation phase is designed in order to mitigate the main limitations of state-of-the-art algorithms, such as the utilization of thresholding techniques [14,15,23].For low threshold levels, in fact, many peaks can come up, mostly caused by the integration of a noisy clutter (which some authors [21,22] try to mitigate with image preprocessing).If the threshold is set high, weak structures can be masked if: (1) shorter than image/subset dimension and/or (3) fainter than the noisy clutter.In addition, thresholding techniques typically require intense customization, which is difficult to automate.
The algorithm herein developed (Figure 1) still follows a wake-based approach because of its maturity and assessed capability.First, a Radon-based detection phase identifies a single trough relative to the turbulent wake and two pairs of peaks, representative of narrow-V and Kelvin wakes.Then, both dark and bright dots are validated as representative of real wake structures by merit indexes.
The presented main algorithm originalities are: • in the Radon domain, peaks and troughs corresponding to the typical five wake components are searched and always identified according to the hydrodynamic theory [17] (no threshold logic), exploiting the physical characteristics in terms of relative angular distance; • ship signature is masked and not summed in the Radon integral; • the turbulent and one narrow-V wakes are contemporaneously identified as a pair; • merit indexes are defined in the intensity domain and are utilized to get the final confirmation of the detected wake features.

Wake Component Detection
The ship locations are known as the center of the cluster of bright pixels affected by ship returns, which can be detected by the ship detection techniques [7].Here, the constant false alarm rate (CFAR) approach [24,25] has been applied.The algorithm input is a ship-centered tile, selected from the entire image whose sizes are set with reference to the physical expected length of wakes (3 to 5 km for the turbulent wake [19]) or less if the ship is closer to image borders.
The ship-centered tile shows the bright pixels not only relative to the ship return pixels, but also to the land signature, as well as the bright disturbances due to ship propellers, typically imaged within a distance of about one ship length aft [26].Such bright returns have to be removed.Differently from Rey et al. [23], who propose to detect ship bright returns by setting a threshold value and to replace them with the mean value of the amplitude distribution, hereinafter, the tile is masked with a rectangular area comprised between [−δrmax, δrmax] (range) and [−δamax, δamax] (azimuth), where δrmax It is worth noting that a Radon domain point only identifies a line.Therefore, the above five lines also extend before the bow, while the wake features are half-lines starting from the ship's true position; hence, the so-called 180 ˝ambiguity in ship heading results and has to be solved (see Sections 2.3 and 3 for details).After wake component identification and validation, ship velocity is estimated by both techniques of azimuth shift and Kelvin pattern wavelength.
The presented main algorithm originalities are: ‚ in the Radon domain, peaks and troughs corresponding to the typical five wake components are searched and always identified according to the hydrodynamic theory [17] (no threshold logic), exploiting the physical characteristics in terms of relative angular distance; ‚ ship signature is masked and not summed in the Radon integral; ‚ the turbulent and one narrow-V wakes are contemporaneously identified as a pair; ‚ merit indexes are defined in the intensity domain and are utilized to get the final confirmation of the detected wake features.

Wake Component Detection
The ship locations are known as the center of the cluster of bright pixels affected by ship returns, which can be detected by the ship detection techniques [7].Here, the constant false alarm rate (CFAR) approach [24,25] has been applied.The algorithm input is a ship-centered tile, selected from the entire image whose sizes are set with reference to the physical expected length of wakes (3 to 5 km for the turbulent wake [19]) or less if the ship is closer to image borders.
The ship-centered tile shows the bright pixels not only relative to the ship return pixels, but also to the land signature, as well as the bright disturbances due to ship propellers, typically imaged within a distance of about one ship length aft [26].Such bright returns have to be removed.Differently from Rey et al. [23], who propose to detect ship bright returns by setting a threshold value and to replace them with the mean value of the amplitude distribution, hereinafter, the tile is masked with a rectangular area comprised between [´δr max , δr max ] (range) and [´δa max , δa max ] (azimuth), where δr max is 1.5-times the length of the ship cluster projected onto the range direction and δa max is the maximum expected azimuth shift, computed from Equation (1) [12].
where R is the slant range, V the SAR antenna velocity, v g and v r are the ship velocity along the ground and slant range, respectively, β i the local incidence angle and α the angle between the ship route and ground range direction.To calculate δa max , the maximum registered ship speed (v gmax ) in the area is assumed, as provided by the information [27] on the typical local maritime traffic, whereas R, V and β i are retrieved from the ancillary data.The Radon transform [18] is then applied to the ship-centered, masked image tile (Step 1 in Figure 1).The Radon transform integrates image brightness along straight lines identified by their distance(s) from the origin of the image reference frame and by the angle (θ) between the normal to the straight line and the image x-axis.Therefore, one bright (dark) line-shaped feature in the input image leads to one peak (trough) in the Radon domain (θ,s), whose extent in (θ,s) depends on the line's width.Hereinafter, the Radon transform is computed for any (θ,s) by summing the intensity values of unmasked pixels only.Consequently, the integration path length usually considered in the standard Radon operator is not applicable for the lines passing through the mask.To take this effect into account, each summation is divided by the actual length of the integration path (portion of the line outside the mask) so that masked pixels contribute neither to summation, nor to path length.Finally, the Radon domain is restricted (Step 2.1 in Figure 1) between two sine curves [14] to only detect peaks relative to the wake of the ship under investigation (δa ď δa max ), imposing: |s| ď δa max sinθ.
The algorithm core consists of a sequential procedure to identify 4 peaks and one trough corresponding to the 5 potentially imaged wake structures of interest: the turbulent wake (the trough), the two narrow-V wakes and the two Kelvin arms.The peaks/trough are only searched where they are expected to be, opportunely limiting the search area in terms of (θ,s).First, a trough/peak pair is searched, corresponding to the turbulent wake and one narrow-V wake (herein defined as "first narrow-V wake"), which must be within a 4 ˝θ-interval (Figure 2 and Step 2.2 in Figure 1).To this end, the whole tile is scanned by a window sized with ∆θ v = 4 ˝+ δθ and ∆s v1 = δa max + δs, where δθ accounts for the angular discretization step in Radon transform derivation, and δs includes both the linear discretization step and the width of the wake.The trough/peak pair is selected as the one that maximizes their difference.When the turbulent wake direction (θ t , s t ) and the first narrow-V wake direction (θ v1 , s v1 ) are known, the second narrow-V wake is identified (Figure 2 and Step 2.3 in Figure 1) as the peak in (θ v2 , s v2 ) within ∆θ v from θ t on the opposite side with respect to the first narrow-V wake.Depending on θ t , s v2 can be smaller or larger than s t .In particular, s v2 is searched either from s t to s t + ∆s v2 or from s t ´∆s v2 to s t , from s t ´∆s v2 to s t + ∆s v2 , otherwise, with ∆s v2 = s t ∆θ v /tanθ t + δs, according to the value of θ t , as shown in Table 1.For the sake of clarity, being s v2 = δa¨sinθ v2 , if θ v2 is counterclockwise with respect to θ t , s v2 is smaller than s t if θ t < 90 ˝, s v2 is larger than s t if θ t > 90 ˝+ ∆θ and s v2 is comprised between s t ´∆s v2 and s t + ∆s v2 if 90 ˝< θ t < 90 ˝+ ∆θ.

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been introduced to: (1) assess whether the identified line-shaped features refer to the image of a wake structure; (2) define a merit index for the confirmed wake features.Validation is performed in the input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transform (Step 3 in Figure 1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in Figure 1).To this end, the lines of the wake features are intersected with the constant-range line passing through the focused ship location.Then, the turbulent wake is identified as the half-line with lower mean intensity.The narrow-V and Kelvin arm half-lines are identified as the ones within 90° from the unambiguous turbulent wake half-line.

Second Narrow-V search area First Kelvin arm search area
Turbulent wake/ First Narrow-V wake pair Δθv Δθk2

Δθk1
If θ v1 is clockwise with respect to θ t , θ v2 is counterclockwise with respect to θ t Table 1.Search logic for the second narrow-V wake.
θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been introduced to: (1) assess whether the identified line-shaped features refer to the image of a wake structure; (2) define a merit index for the confirmed wake features.Validation is performed in the input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transform (Step 3 in Figure 1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in Figure 1).To this end, the lines of the wake features are intersected with the constant-range line passing through the focused ship location.Then, the turbulent wake is identified as the half-line with lower mean intensity.The narrow-V and Kelvin arm half-lines are identified as the ones within 90° from the unambiguous turbulent wake half-line.

Second Narrow-V search area First Kelvin arm search area
Turbulent wake/ First Narrow-V wake pair Δθv Δθk2

Δθk1
If θ v1 is counterclockwise with respect to θ t , θ v2 is clockwise with respect to θ t Table 1.Search logic for the second narrow-V wake.
θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon dom each one identified by the coordinates (θ,s) and representing a candidate component of the typ wake appearance, which must be validated as a real wake arm.The validation phase has introduced to: (1) assess whether the identified line-shaped features refer to the image of a w structure; (2) define a merit index for the confirmed wake features.Validation is performed in input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transf (Step 3 in Figure 1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in Figure 1 this end, the lines of the wake features are intersected with the constant-range line passing thro the focused ship location.Then, the turbulent wake is identified as the half-line with lower m intensity.The narrow-V and Kelvin arm half-lines are identified as the ones within 90° from unambiguous turbulent wake half-line.Finally, in order to detect the Kelvin boundaries directions (θ k1 , s k1 ) and (θ k2 , s k2 ), it is assumed that the first (second) Kelvin arm is found clockwise (counterclockwise) with respect to the turbulent wake (Figure 2 and Step 2.4 in Figure 1).In terms of θ k1 (θ k2 ), the Kelvin arm is identified by the peak within ∆θ k = 19.5 ˝+ δθ from θ t clockwise (counterclockwise).In terms of s k1 (s k2 ), it has to be searched similarly to s v2 and according to Table 2 with ∆s k = s t ∆θ k /tanθ t + δs.

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been introduced to: (1) assess whether the identified line-shaped features refer to the image of a wake structure; (2) define a merit index for the confirmed wake features.Validation is performed in the input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transform (Step 3 in Figure 1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in Figure 1).To

Wake Component Validation
The The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been introduced to: (1) assess whether the identified line-shaped features refer to the image of a wake structure; (2) define a merit index for the confirmed wake features.Validation is performed in the input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transform (Step 3 in Figure 1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in Figure 1).To this end, the lines of the wake features are intersected with the constant-range line passing through the focused ship location.Then, the turbulent wake is identified as the half-line with lower mean intensity.The narrow-V and Kelvin arm half-lines are identified as the ones within 90° from the unambiguous turbulent wake half-line.

Wake Component Validation
The results of the previous peaks/trough detection phase are five points in the Radon domain, each one identified by the coordinates (θ,s) and representing a candidate component of the typical wake appearance, which must be validated as a real wake arm.The validation phase has been introduced to: (1) assess whether the identified line-shaped features refer to the image of a wake structure; (2) define a merit index for the confirmed wake features.Validation is performed in the input domain (intensity image).
The five points are converted into five lines in the input image by the inverse Radon transform (Step 3 in Figure 1), and the 180 ˝ambiguity in ship heading must be solved (Step 4.1 in Figure 1).To this end, the lines of the wake features are intersected with the constant-range line passing through the focused ship location.Then, the turbulent wake is identified as the half-line with lower mean intensity.The narrow-V and Kelvin arm half-lines are identified as the ones within 90 ˝from the unambiguous turbulent wake half-line.
The validation phase is carried out by means of merit indexes.Here, only the merit index in the intensity domain is included, since merit indexes in the Radon domain [28] have shown worse performance in terms of confirmation/rejection of wake components.For bright features, the merit index F I indicates to what extent the considered half-line is brighter than the mean intensity of the masked image subset.The merit index is calculated along each detected half-line representing narrow-V and Kelvin wakes, discarding the 5% brighter pixels of the half-line to filter out the effect of localized brighter spots along dimmer lines.It is estimated by Equation (2).
where CDF(x i ), the cumulative distribution function of the image intensity (x i ) along the identified half-line, is calculated as the summation between x 0 and x i of the histogram of the intensity along the half-line.CDF ´1(0.95) is the intensity value (x i ), which corresponds to CDF = 0.95.In addition, n halfline is the number of pixels of the half-line, and I is the intensity mean value over the whole, masked tile.As defined, the merit index limits the erroneous confirmation of the identified line-shaped features, which could show a mean value larger than the mean intensity of the masked image subset only because of the effect of passing through a bright cluster of pixels.It is worth noting that it has been chosen to discard the 5% brighter pixels as a compromise between the erroneous confirmations of features not related to wakes and the erroneous no-confirmation of bright, but short line-shaped features related to wakes.
A peak in the Radon domain is assumed to be representative of a real wake feature depending on its F I value.In Section 2.4, it will be clear that, for the narrow-V wakes (and for the turbulent wake, discussed later in this section), it is desirable to avoid missed confirmations, even if some false confirmations could arise.Therefore, a narrow-V wake arm is retained as part of the wake and defined as a "confirmed wake feature" if F I is positive and F I is used as a merit index.Otherwise, the line-shaped feature is discarded.On the contrary, it is better to limit the possibility to attain false confirmations of Kelvin wake arms.Therefore, a margin of 33% is assumed: a Kelvin arm is confirmed if F I > 0.33.
Differently, for the turbulent wake, the merit index F I t defines to what extent the half-line is darker than the mean intensity of the masked image subset.In particular, since the 95% cutting strategy does not make sense for the turbulent, dark wake, because it would lead to optimistic estimation, the following parameter is used: where I t is the intensity mean value over the turbulent half-line.A trough is assumed to be representative of a real wake feature and defined as a "confirmed wake feature" if F I t is negative.If the estimated turbulent and the first narrow-V wake are not confirmed, the whole wake structure is considered as not imaged, and the procedure is stopped.In details, for all of the ships whose wakes are not visible in the image, the algorithm should exclude wake appearances, providing positive F I t and/or negative F I relative to the first narrow-V wake.

Ship Heading and Velocity
The last algorithm module is performed to estimate ship heading and velocity.The former one is known when the turbulent wake direction is calculated, the azimuth direction w.r.t. the True North being included in ancillary data; whereas, ship velocity is estimated by state-of-the-art techniques of azimuth shift [11][12][13] and Kelvin pattern wavelength.
A procedure is herein introduced for wake vertex estimation with the aim of improving the accuracy of the azimuth displacement computation.In detail, the proposed merit index-based weighting strategy aims to mitigate the influence on the vertex position of wakes with low merit index.First of all, since the wake structures visible in SAR images develop after 1.5 ship lengths, the constant range line is selected at range δr max before the ship cluster center.Then, its intersections with the confirmed wake structures are determined.The wake vertex is then estimated by means of the weighted least squares approach, where merit indexes are used as weights.Finally, ship velocity slant-range component v r and/or ground component v g is obtained from Equation ( 1), where δa is now the imaged azimuth displacement, i.e., the distance between the estimated wake vertex and the ship (corresponding to the center of the tile) along the azimuth direction.
Moreover, another wake-based technique to estimate ship velocity exploits the measurements of the wavelength of the Kelvin waves.In particular, it can be derived that ship velocity is related to the wavelength of the Kelvin waves [14,15]: where g is gravity acceleration and q = 1 for Kelvin transverse waves [14] or q = ?3{2 for cusp waves [15].In [15], Fourier analysis is performed along the Kelvin arms to calculate λ as the wavelength at the spectrum peak, and finally, the velocity is derived through Equation ( 4), averaging three spectrograms along three parallel cuts (the middle one coincides with Kelvin arms detected by the Radon analysis) with 1-pixel separation.The application of the method showed that the velocity can be estimated with accuracy typically ranging from 4% to 9%.
It is worth underlining that false confirmations with low merit indexes do not pose serious drawbacks on velocity estimation based on the azimuth displacement method because they slightly impact wake vertex estimation thanks to the index-based weighting strategy.Nonetheless, false confirmations of Kelvin wake arms determine misleading velocity estimations based on the wavelength method because the merit index does not have any beneficial effect in this case.Consequently, it the above-mentioned margin of 33% on the merit index of Kelvin arms has been assumed (Section 2.3).

Results and Discussion
The algorithm has been tested over 13 SAR images of the Gulf of Naples, Italy, from COSMO/SkyMed (Table 3) and TerraSAR-X orbiting SARs (Table 4).The COSMO/SkyMed images are Stripmap products in HH or VV polarizations and with three-meter ground range, three-meter azimuth resolutions.The TerraSAR-X images are "TanDEM-X Co-registered Single Look Slant Range Complex Experimental bistatic acquisitions".The VV polarization, monostatic products only are used, whose geometric resolutions (1.2-meter slant range, 6.6-meter azimuth) coincide with those of the standard TerraSAR-X dual-polarization, Stripmap products.In these images, it is possible to identify 34 wakes by visual inspection, with a variable number of imaged features (turbulent, narrow-V, Kelvin), which are used to test the algorithm.
Considering typical turbulent wake length and image resolutions, processed image tiles should have a dimension of about 3000 pixels.Instead, due to the fact that detected ships associated with the visible wakes are closer to image borders, image tiles are selected as large as 2000 ˆ2000 pixels.With reference to masks, δa max is determined for any detected ship by using ancillary data for the satellite velocity and to estimate the local slant range by interpolation.In addition, 30 knots is selected as the maximum ship velocity to be expected on the grounds of the scheduled maritime traffic in the Gulf of Naples, Italy.As for δr max , it is estimated from the ship cluster dimension for each detected ship.For the processing validations, δa max and δr max vary in the ranges of 156 to 428 pixels (618 pixels for wake ID #29-30) and 18 to 67 pixels, respectively.The end-to-end wake detection process is presented for Wake #30 (Table 4), imaged by TerraSAR-X.The wake (Figure 3a) has been chosen because: (1) it is very close to another ship wake, to test the robustness of the proposed method with respect to critical situations where multiple wakes appear; (2) only one narrow-V wake and one Kelvin arm are distinguishable, even if very weak, in addition to the turbulent wake.This is useful to test the proposed method's capability to confirm only the imaged features.
Remote Sens. 2016, 8, 498 9 of 16 addition to the turbulent wake.This is useful to test the proposed method's capability to confirm only the imaged features.An intensity image subset is generated around the identified ship and properly masked (red rectangle in Figure 3a).The masked intensity subset is then Radon-transformed (Figure 3b).The input intensity image being square, the transform is available only in a butterfly-shaped region out of which (red area) it is not defined.The area enclosed by the red sinusoidal curves represents the sub-domain where wake structures shall be located, if imaged, considering the foreseen maximum azimuth shift.The peak/trough search phase is carried out, and the identified wake structures are shown in Figure 3b.The blue triangle is the candidate turbulent wake; the two green squares are the candidate narrow-V wake arms; and the two yellow circles are the candidate Kelvin wake arms.Finally, the inverse Radon transform is applied to identify the wake direction, plotted as lines superimposed on the intensity image (Figure 3c). Figure 3d shows the confirmed half-lines, after the 180° ambiguity solution and the estimation of the merit indexes.It is worth recalling that the turbulent wake is An intensity image subset is generated around the identified ship and properly masked (red rectangle in Figure 3a).The masked intensity subset is then Radon-transformed (Figure 3b).The input intensity image being square, the transform is available only in a butterfly-shaped region out of which (red area) it is not defined.The area enclosed by the red sinusoidal curves represents the sub-domain where wake structures shall be located, if imaged, considering the foreseen maximum azimuth shift.
The peak/trough search phase is carried out, and the identified wake structures are shown in Figure 3b.The blue triangle is the candidate turbulent wake; the two green squares are the candidate narrow-V wake arms; and the two yellow circles are the candidate Kelvin wake arms.Finally, the inverse Radon transform is applied to identify the wake direction, plotted as lines superimposed on the intensity image (Figure 3c). Figure 3d shows the confirmed half-lines, after the 180 ˝ambiguity solution and the estimation of the merit indexes.It is worth recalling that the turbulent wake is confirmed if F I t is negative; the narrow-V wakes are confirmed if F I is positive; and the Kelvin arms are confirmed if F I > 0.33.
The first result is that the algorithm correctly identifies the imaged ship wake, discarding the other, longer dark line, not related to the ship under investigation.In [15], the dark extremum of the Radon transform values is related to the turbulent wake.For the wake in Figure 3, the absolute minimum of the Radon transform values is not located in the blue triangle of Figure 3b.Consequently, any technique not using the pair-search logic for the turbulent and narrow-V wake components fails the detection in this case.
Then, the 180 ˝heading-ambiguity is solved, and merit indexes correctly validate all visible features.In detail, the merit index relative to turbulent wake is negative, whereas a positive index is shown for only one narrow-V wake.Since the first and second Kelvin arms show a negative and 0.02 merit index, respectively, they are not confirmed.Then, the ship heading is estimated at 7 ˝with respect to the azimuth direction.
Velocity cannot be evaluated because the azimuth displacement method is not applicable (the angle α between the ship moving direction and range is greater than 75 ˝), as well as the Kelvin pattern method, since no Kelvin arms are detected.
In order to have an overall view of wake detectability, ships related to the 34 visible wakes are detected in the images, and 34 image subsets are generated to verify whether the algorithm automatically detects existing wakes.The developed confirmation criterion in intensity domains is applied over 170 detected wake features of the 34 wakes in the Stripmap SAR images (Tables 3 and 4).The wake detection results are validated through image visual inspection.
First of all, no missed confirmations are obtained for turbulent and narrow-V wakes, i.e., all of the bright (dark) wake features really imaged show a positive (negative) merit index.Among the Kelvin wake arms, 22 cases show a positive merit index, even if in most cases, they should be discarded, because visual inspection revealed that the detected arms could not be confirmed.Therefore, no false confirmations occur thanks to the strategy of confirming only Kelvin arms when F I > 0.33.As a drawback, in one case, a real Kelvin wake (although with a weak appearance) is also discarded.Because it refers to the only case (Wake #1) for which the first Kelvin arm shows a merit index of 0.68, it does not impact on the velocity estimation capability.
The proposed method for wake detection is robust with respect to critical situations where multiple wakes appear as verified for Wake #30.The pair-search logic shows also an additional advantage when SAR images are characterized by dark areas not related to wake features.On some occasions, it is possible that these features restitute the absolute minimum in the restricted area of the Radon domain, which would lead to erroneous turbulent wake identification.It is worth noting that since the method assumes a wake structure composed of one straight turbulent wake, two straight narrow-V wakes and two straight Kelvin arms, it is not able to correctly reconstruct wakes with different structures, as the one shown in Figure 4, which shows a lack of straight features and more than one dark wake component.
With reference to wake appearance, it is worth noting that the case of the turbulent and one narrow-V wakes comprises the vast majority (24), whereas the turbulent and two narrow-V are identified in nine cases.In only one case (Wake #1), a Kelvin arm is detected, in addition to the turbulent and two narrow-V wakes.Validation of ship velocity with sea truth (Automated Identification System (AIS) data) is possible in seven cases.The AIS data have been provided by exactEarth Ltd. (Cambridge, ON, Canada) and interpolated to the SAR image acquisition time by using the Hermite cubit spline [29] and including the azimuth shift due to ship velocity components along the slant range.As an example, Wake #33 (Table 4) is shown in Figure 5a.The ship detected in a TerraSAR-X image is associated with the "BF Catania" (container ship, 102 m × 19 m) with available AIS information at the time of image acquisition (heading 258° w.r.t.True-North, velocity 9.6 knots).The confirmed wake components are shown in Figure 5b.SAR processing outputs are 258.4°as the heading and 10.2 knots as the velocity (azimuth displacement method), with a 6.25% overestimation for the latter.Additionally, for Wake #1 in Table 3 (Figure 5c), the wake features are detected and correctly confirmed (Figure 5d).Since Wake #1 is the only one included in the dataset, showing a high merit index (0.68) for a Kelvin arm, both azimuth displacement and Kelvin wake pattern methods can be applied.They lead to 27.5 knots and 24.1 knots, respectively, with a difference of 12%.It is worth noting that the absolute minimum of the Radon transform values is related to the dark stretched spot on the left part of the tile, and the turbulent wake direction is correctly identified thanks to the pairsearch logic.Furthermore, in such a case, the thresholding techniques [15] fail the wake reconstruction.
On the basis of the time of image acquisition and the geolocation of the ship cluster, it is possible to associate the detected ship with the "Don Francesco" (high speed craft, catamaran; 60 m × 16.5 m; maximum velocity 39 knots; en route from Castellammare to Capri), but unfortunately, no AIS data are available.
The last two columns of Tables 3 and 4 list the ship velocity along the heading direction estimated by the azimuth displacement method and provided by the AIS data, respectively.Generally, the ship velocity errors range from 2% to 12% with a single out-of-range value of 24% relative to Wake #11 (Table 3, Figure 6).The larger error could be due to the ship location, which is very close to the windy west coast of Capri Island and/or to a ship heading of 27° with respect to the azimuth direction.The last columns of Tables 3 and 4 list velocity estimation results, which have been performed 10 times by means of the azimuth displacement method, with the other 24 cases discarded because the angle between the ship heading and azimuth direction is below 15 Validation of ship velocity with sea truth (Automated Identification System (AIS) data) is possible in seven cases.The AIS data have been provided by exactEarth Ltd. (Cambridge, ON, Canada) and interpolated to the SAR image acquisition time by using the Hermite cubit spline [29] and including the azimuth shift due to ship velocity components along the slant range.As an example, Wake #33 (Table 4) is shown in Figure 5a.The ship detected in a TerraSAR-X image is associated with the "BF Catania" (container ship, 102 m ˆ19 m) with available AIS information at the time of image acquisition (heading 258 ˝w.r.t.True-North, velocity 9.6 knots).The confirmed wake components are shown in Figure 5b.SAR processing outputs are 258.4˝as the heading and 10.2 knots as the velocity (azimuth displacement method), with a 6.25% overestimation for the latter.Additionally, for Wake #1 in Table 3 (Figure 5c), the wake features are detected and correctly confirmed (Figure 5d).Since Wake #1 is the only one included in the dataset, showing a high merit index (0.68) for a Kelvin arm, both azimuth displacement and Kelvin wake pattern methods can be applied.They lead to 27.5 knots and 24.1 knots, respectively, with a difference of 12%.It is worth noting that the absolute minimum of the Radon transform values is related to the dark stretched spot on the left part of the tile, and the turbulent wake direction is correctly identified thanks to the pair-search logic.Furthermore, in such a case, the thresholding techniques [15] fail the wake reconstruction.
On the basis of the time of image acquisition and the geolocation of the ship cluster, it is possible to associate the detected ship with the "Don Francesco" (high speed craft, catamaran; 60 m ˆ16.5 m; maximum velocity 39 knots; en route from Castellammare to Capri), but unfortunately, no AIS data are available.
The last two columns of Tables 3 and 4 list the ship velocity along the heading direction estimated by the azimuth displacement method and provided by the AIS data, respectively.Generally, the ship velocity errors range from 2% to 12% with a single out-of-range value of 24% relative to Wake #11 (Table 3, Figure 6).The larger error could be due to the ship location, which is very close to the windy west coast of Capri Island and/or to a ship heading of 27 ˝with respect to the azimuth direction.4; (b) confirmed half lines for Wake #33 in Table 4; (c) masked image tile for Wake #1 in Table 3; (d) confirmed half lines for Wake #1 in Table 3.
Finally, 10 ships not related to wakes are detected, and the end-to-end processing is applied.The corresponding ten image subsets are generated to confirm that the algorithm excludes wake presence.As an example, Figure 7 portrays a ship signature that is not associated with a wake signature.The results show that for the 10 subsets, the figure of merit relative to the turbulent wake is positive and/or F I relative to the first narrow-V wake is negative.Consequently, the whole wake structure is correctly considered as not present.3.
situations, such as multiple wakes' appearance or dark areas not related to wake features.Indeed, the algorithm has been applied to an image subset including two dark turbulent wakes of two ships cruising in the bottom-up direction.The algorithm correctly identifies the wake of the ship under investigation, discarding the other, longer turbulent wakes.It is the bright/dark pair search logic that allows one to discard the additional dark features, whose bright narrow-V wake arms have been already dampened.Ship velocity is then estimated by state-of-the art techniques of azimuth shift and Kelvin pattern wavelength.Results' validation with truth-at-sea (Automated Identification System data) has been possible in seven cases, where it is verified that ship velocity is estimated by the azimuth shift technique with errors typically below 12%, with only one case of velocity error around 24% for a ship very close to the Capri Island coast.Finally, another 10 ships not related to wakes are detected, and 10 image subsets are generated.The processing is also applied to the 10 subsets, and the results show that all wake-like features are correctly rejected.

Figure 1 .
Figure 1.The flow chart of the proposed method.

Step 5 : 1 :Figure 1 .
Figure 1.The flow chart of the proposed method.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Table 1 .Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Table 1 .
Search logic for the second narrow-V wake.θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt If θv1 is counterclockwise with respect to θt, θv2 is clockwise with respect to θt sv2 Selection

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 2 .
Figure 2. Angular search logic for narrow-V wakes and Kelvin arms.

Figure 3 .
Figure 3. End-to-end wake detection process applied to Wake #30 of Table 4: (a) masked image tile; (b) Radon transform with detected wake components (blue triangle: turbulent wake, two green squares: narrow-V wake arms, two yellow circles: Kelvin wake arms); (c) detected lines superimposed to the intensity image; (d) confirmed half lines.

Figure 3 .
Figure 3. End-to-end wake detection process applied to Wake #30 of Table 4: (a) masked image tile; (b) Radon transform with detected wake components (blue triangle: turbulent wake, two green squares: narrow-V wake arms, two yellow circles: Kelvin wake arms); (c) detected lines superimposed to the intensity image; (d) confirmed half lines.

Table 1 .
Search logic for the second narrow-V wake.

Table 2 .
Search logic for the Kelvin arms.
θk1 Selection θk1 is clockwise with respect to θt sk1 Selection

Table 2 .
Search logic for the Kelvin arms.
θk1 Selection θk1 is clockwise with respect to θt sk1 Selection

Table 2 .
Search logic for the Kelvin arms.

Table 2 .
Search logic for the Kelvin arms.

Table 1 .
Search logic for the second narrow-V wake.θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt If θv1 is counterclockwise with respect to θt, θv2 is clockwise with respect to θt sv2 Selection

Table 2 .
Search logic for the Kelvin arms.

Table 1 .
Search logic for the second narrow-V wake.
θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt If θv1 is counterclockwi respect to θt, θv2 is clockw respect to θt sv2 Selection

Table 2 .
Search logic for the Kelvin arms.

Second Narrow-V search area First Kelvin arm search area Turbulent wake/ First Narrow-V wake pair Δθv Δθk2 Second Kelvin arm search area Δθk1 s k1 Selection if
results of the previous peaks/trough detection phase are five points in the Rad each one identified by the coordinates (θ,s) and representing a candidate component of wake appearance, which must be validated as a real wake arm.The validation phas introduced to: (1) assess whether the identified line-shaped features refer to the image structure; (2) define a merit index for the confirmed wake features.Validation is perfor input domain (intensity image).The five points are converted into five lines in the input image by the inverse Radon (Step 3 in Figure1), and the 180° ambiguity in ship heading must be solved (Step 4.1 in F θ t < 90 ˝´∆θ k s t < s k1 <s t + ∆s k if 90 ˝´∆θ k < θ t < 90 ˝st ´∆s k < s k1 < s t + ∆s k if θ t > 90 ˝st ´∆s k < s k1 < s t

Table 1 .
Search logic for the second narrow-V wake.

Table 2 .
Search logic for the Kelvin arms.

Table 1 .
Search logic for the second narrow-V wake.
θv2 Selection If θv1 is clockwise with respect to θt, θv2 is counterclockwise with respect to θt If θv1 is counterclockwise with respect to θt, θv2 is clockwise with respect to θt sv2 Selection

Table 2 .
Search logic for the Kelvin arms.

First Kelvin arm search area Turbulent wake/ First Narrow-V wake pair Δθv Δθk2 Second Kelvin arm search area Δθk1 s k2 Selection if
θ t < 90 ˝st ´∆s k < s k2 < s t if 90 ˝< θ t < 90 ˝+ ∆θ k s t ´∆s k < s k2 < s t + ∆s k if θ t > 90 ˝+ ∆θ k s t < s k2 < s t + ∆s k