Next Article in Journal
Using Imagery Collected by an Unmanned Aerial System to Monitor Cyanobacteria in New Hampshire, USA, Lakes
Previous Article in Journal
Methodology for Mapping the Ecological Security Pattern and Ecological Network in the Arid Region of Xinjiang, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ocean Clutter Characterization Based on PolSAR Data and Second-Order Statistics of Elementary Scatterers

by
Georgia Koukiou
* and
Vassilis Anastassopoulos
*
Electronics Laboratory, Physics Department, University of Patras, 26504 Patras, Greece
*
Authors to whom correspondence should be addressed.
Remote Sens. 2023, 15(11), 2837; https://doi.org/10.3390/rs15112837
Submission received: 2 April 2023 / Revised: 21 May 2023 / Accepted: 27 May 2023 / Published: 30 May 2023
(This article belongs to the Section Ocean Remote Sensing)

Abstract

:
New features are proposed for sea clutter characterization when PolSAR data are employed. Cameron coherent target decomposition is applied to characterize each SAR pixel by means of the 8 basic elementary scatterers obtained by this decomposition. Since the examined SAR pixel does not match exactly to a specific ideal elementary scatterer, the closest scatterer is assigned to this pixel along with its distance (closeness) from the ideal one. The features proposed are (a) the percentage of each elementary scatterer in the investigated region; (b) the distribution of the closeness of each declared elementary scattering mechanism to its ideal counterpart; (c) the co-occurrence of the same scattering mechanism, taking into consideration its distance from the ideal one; and (d) The co-occurrence of the different scattering mechanisms in pairs, taking again into consideration their distance from the ideal ones. Simultaneously, the skewness and the kurtosis and their 2D versions of the previously mentioned probabilities are employed to further improve ocean clutter characterization. The above features correspond to the sea state condition in each separate region, i.e., the wave height and the wind speed. A clear correspondence between the proposed features and the sea state conditions is established. Data were available from RADARSAT-2 and ALOS-1 PALSAR systems for four different ocean regions on Earth.

1. Introduction

Radar technology has been used to support remote sensing procedures since the Second World War. With the launching of the first satellites, radar and especially SAR contributed to monitoring and recording the surface of the Earth and, simultaneously, to surveillance and detection of targets of interest. Nowadays, Fully Polarimetric SAR imaging (PolSAR) supports detailed surveillance of the Earth’s surface, either land or sea, based mainly on the well-established pixel’s scattering matrix decomposition methods.
The statistical characterization of the surface of the sea started more than three decades ago, discriminating radar and SAR echoes according to sea state conditions. The K–pdf [1,2,3,4,5] is a compound statistical model describing the radar echoes from the sea surface when the mean value of the pixels is modulated by a different pdf (gamma pdf) than the one describing each separate pixel (Rayleigh pdf). However, for high-resolution radar, each cell contains a small number of scatterers and, consequently, the Rayleigh pdf is not an appropriate distribution to represent each separate pixel due to the weakness of the Central Limit Theorem for application in the case of small number of scatterers. Accordingly, the Generalized Gamma pdf is used instead of the Rayleigh distribution, and, additionally, the same pdf is employed to model the modulation component, leading in this way to the Generalized Compound Model (GC-pdf) [1,2,3]. In general, the statistics of each separate pixel deviate from Gaussian behavior (Rayleigh amplitude) in the case that, in the area of the pixel, we have a small number of scatterers or one of them is much stronger than the rest [5,6].
With the advancement of PolSAR imaging, Earth’s surface was extensively studied, while numerous mathematical approaches for sea state characterization and target detection were developed. In [7], a method is proposed for multidimensional speckle rejection. The method is based on two new matrix models which arise from applying the stochastic summation approach to PolSAR. Methods were investigated in [8] which use PolSAR image data to measure ocean slopes and wave spectra. Another method to measure ocean wave slope spectra using PolSAR data was developed in [9] without the need for a complex hydrodynamic modulation transform function. In [10], the authors considered hybrid PolSAR data of ocean surface slicks and hypothesized that we can design a system that is able to discriminate between mineral oil, plant oil, and clean sea. The authors in [11] reported the results for two polarimetric methods implemented to evaluate their ship detection capabilities. In [12], an efficient algorithm is presented for retrieving the ocean surface wind vector from C-band Radar Satellite RADARSAT-2 PolSAR measurements based upon the co-polarized geophysical model function. A novel method is proposed in [13] to estimate the seawater surface salinity and temperature from the data collected by PolSAR imagery. In [14], PolSAR data are used to estimate the energy dissipation rate under different wind and sea conditions. A shallow-sea topography detection method was proposed based on fully polarimetric Gaofen-3 SAR data [15].
PolSAR images are employed in [16] for discriminating between oil spills and seawater based on the polarization ratio. This polarization ratio method is based on the difference between the scattering mechanism and the dielectric constant for the oil spills compared to that of seawater. For oil spill monitoring, the satellite revisit time needs to be as short as possible to identify minor spills before they can cause widespread damage [17]. Simultaneously, it is necessary to capture a sufficient amount of information about the surface to clearly distinguish between oil-spilled and oil-free sea regions. The hybrid-polarimetry synthetic aperture radar (SAR) system is exploited in [17] for such capabilities. Hybrid-polarized-mode compact polarization (CP) SAR imagery [18] will soon be available with the launch of the RADARSAT Constellation Mission. A methodology is presented in [18] to retrieve the oil–water mixture ratio at the ocean surface using CP SAR imagery. The advantage of the proposed new SAR system is that CP images will have wider swath and shorter revisit time compared to quad-polarization images, which are presently available from space-borne and air-borne SAR.
Millimeter and centimeter electromagnetic waves employed in radar technology can provide information regarding the types of electromagnetic scatterers on the Earth’s surface. In particular, coherent and fully polarimetric acquisitions can give detailed information regarding the Earth’s surface composition. Accordingly, PolSAR technology is widely used in land cover classification, security tasks, and surveillance [19,20,21,22,23].
Second-order statistics such as the transition matrices corresponding to Markov chains can give better insight into the correlation properties among the PolSAR pixels. Accordingly, the use of discrete Markov processes [3,24,25] outperforms, in detection and classification tasks, conventional distributions for clutter characterization.
Lee and Pottier provide, in [26], a comprehensive interpretation of the polarimetry-analyzing methods for target decomposition that are mainly employed. In [27], PolSAR data are analyzed so that the polarimetry properties of scatterers can contribute to man-made objects detection. The scattering mechanisms characterizing the properties of PolSAR pixels are automatically classified in [28]. Three different matrices, namely, the polarimetric covariance, the coherence, and the Muller matrix, are investigated for each PolSAR pixel in [29], so that covariance symmetries in the PolSAR images are established. Nonnegative eigenvalues are confirmed by means of simple modification when the decomposition of the covariance matrices is performed [30]. Coherent Target Decomposition (CTD) techniques were developed [31,32,33] for ship detection and identification making use of Cameron’s CTD [34]. This was feasible since it was possible to characterize the dominant scattering mechanisms on the ship deck as well as on the surface of the sea for various sea states. Kouroupis and Anastassopoulos in [35] re-investigated the co-diagonalization of the Sinclair backscattering matrix in order to overcome issues imposed by Huynen decomposition. In this way, the predominant scattering mechanisms are correctly selected. The work presented in [36] employed elementary scattering mechanisms to characterize the pixels belonging to the surface of a ship. These mechanisms are derived from Cameron CTD when the scattering matrix of PolSAR pixels is analyzed. Additionally, the paper introduces a robust feature to discriminate pixels belonging to a ship which is formed by the transition matrix of the first-order Markov chains [37] derived by the alternation of the elementary scatterers on the surface of the ship. Thus, second-order statistics regarding the appearance of the elementary scatterers are used. The work in [38] presents a CFAR detection approach which is based on the features obtained in [36].
In order to classify each PolSAR pixel from the sea surface into a specific elementary scatterer, Cameron decomposition is employed. The pixels, as well as the region of the Earth’s surface, are categorized according to the percentage of the elementary scatterers’ appearance, their closeness to the corresponding ideal scatterers, and the co-occurrence of the specific elementary scatterers with similar scatterers or with scatterers of a different kind. The closeness of the declared scattering behavior of each pixel to the corresponding ideal elementary scatterer is presented for the whole sea region as a histogram. This gives the capability for evaluating the first four statistical moments of this histogram as well as the 2D skewness and kurtosis of the above-mentioned co-occurrence matrices. It is proved that the sea state conditions affect the statistics and the second-order statistics of the elementary scatterers. A correspondence between the statistical properties of these histograms and the dominating sea conditions is established.
The SNAP software, which is used in our classification procedure, is ideal for analyzing Earth Observation images. It is provided by ESA as open-source software to support its Toolboxes [39]. The SNAP software is appropriate for geocoding and rectification of the remote-sensed images given that ground control points are available. The elementary scatterers map is created when the PolSAR images are geocoded.
The paper is organized as follows. In Section 2, Cameron CTD is analytically presented. In Section 3, the sea PolSAR data are explained in detail. The proposed features as well as the statistical analysis of the PolSAR ocean data and the experimental results are discussed in Section 4. The conclusions given in Section 5 present the main advantages of the proposed method.

2. Cameron Decomposition

In recent decades, target decomposition methods were proposed [26,27,28,29,30,31,32,33,34,35] in order to analyze the scattering matrix of a PolSAR pixel. Among them, the Cameron CTD method analyzes the behavior of each PolSAR pixel and classifies it into one of the elementary scattering mechanisms [34] with physical meaning. Cameron’s coherent decomposition can resolve two essential properties of the PolSAR pixels, namely, reciprocity and symmetry. Reciprocity is an inherent property of monostatic SAR and implies equal non-diagonal elements of the pixel’s backscattering matrix. Furthermore, a reciprocal scatterer is symmetric if it has an axis of symmetry in the plane which is normal on the radar Line of Sight (LOS).
Cameron CTD classifies the elementary scatterers according to the properties of reciprocity and symmetry in six categories, namely: trihedral, dihedral, dipole, cylinder, narrow diplane, and ¼-wave devices. Two other types of elementary scatterers, left and right helix, represent non-symmetric scattering mechanisms, making the total number of Cameron CTD elementary scatterers eight.
As exposed analytically in [36], Cameron decomposition transforms the backscattering matrix S to the backscattering vector S , which is further decomposed into a reciprocal S ^ rec and a non-reciprocal component S ^ non   rec . Moreover, the reciprocal part is being further decomposed into a maximum symmetric component S ^ sym max and a minimum symmetric component S ^ sym min . Accordingly, the backscattering vector S is as follows:
S = g cos θ rec cos ( τ sym ) S ^ sym max + sin ( τ sym ) S ^ sym min + sin θ rec S ^ non   rec
where the tone S is replaced with S ^ for a normalized vector and g represents the total span of the matrix S. The degree that the scatterer deviated from the reciprocal space is determined by θ rec , while the degree of symmetry of the scatterers is described by τ sym . Due to the fact that the reciprocity theorem always applies, we consider θ rec =   0 , and Equation (1) becomes
S = g cos τ sym S ^ sym max + sin τ sym S ^ sym min
The deviation of S from S sym max is expressed by the symmetry degree of the scatterer. A fully symmetric scatterer is represented by S s y m m a x when τ sym = 0 , while the fully asymmetric scatterer corresponds to τ sym = 45 . The component S ^ sym max , which represents maximum symmetry, can be converted into the normalized complex vector Λ ^ ( z ) , which contains z, i.e., the complex parameter describing the scattering mechanism. The complex vector Λ ^ ( z ) is as follows:
Λ ^ z = 1 1 + z 2 1 0 0 z
Table 1 contains the elementary scattering mechanisms and the corresponding values of z. The distance given in the following was introduced by Cameron to determine the underlying scattering mechanisms [40].
d z , z ref = sin 1 min d z , z ref , d * z , z ref
where
d z , z ref = z z ref 2 1 + z 2 1 + z ref 2
and
d * z , z ref = z z ref * 2 + 1 z 2 1 z ref * 2 1 + z 2 1 + z ref 2
to estimate the closeness of the measured scatterer parameter z from the reference parameters z r e f , depicted in Table 1. The symbol |…| stands for the modulus of a complex number, while the superscript * stands for complex conjugation. The scatterer in a specific PolSAR pixel is labelled according to the smallest distance obtained from Equation (4). The distance d z , z ref is employed in Section 4 to formulate the statistical features for sea clutter characterization.
The knowledge to interpret PolSAR imagery based on Cameron’s elementary scatterers is essential for understanding land cover types and their electromagnetic echoes response. Accordingly, areas that are flat like sea or bare land correspond to low-entropy pixel groups, while echoes from high-entropy regions like urban regions are strong. One of the color palettes available in MATLAB (JET color map) was employed to map the elementary scatterers in a PolSAR scene. This mapping, which is perceptually quite distinctive, is shown in Table 2. In [41], other types of land cover decompositions are described.

3. Sea PolSAR Data

3.1. Data Properties

The high-resolution PolSAR imagery comes from two different sources. The first one is the RADARSAT-2 platform (C-band, 5.6 cm), with the main task of all-weather sea and land monitoring. The employed PolSAR imagery is from the Wide Fine Quad-Pol, Single Look Complex, in SLC products [42]. For every PolSAR pixel, with a size of 13.6 × 7.6 m, the four polarimetric echoes are provided in complex numbers, having their I and Q components represented by 16-bit accuracy. The Wide Fine Quad-Pol, Single Look Complex SLC RADARSAT-2 data are acquired at an incidence angle of 30–50°. The SNAP software converts the SLC data into geocoded data which possess all corresponding kinds of polarizations, including their I and Q components. Truth maps needed for training and testing purposes are built by employing the geo-referenced images and their corresponding Google Earth maps. An example of such is the fully polarimetric SAR data containing the Vancouver region, which is depicted in Figure 1. The image in Figure 1a corresponds to the amplitude of the I and Q components of the HV part of the fully polarimetric data, while Figure 1b presents the geocoded SAR image registered using the SNAP software on the Google Earth map.
The second PolSAR data source was the Phase Array L-band Synthetic Aperture Radar (PALSAR) aboard ALOS-1, using the L-band in the polarimetric mode that provides a nominal resolution of approximately 9.4 × 3.6 ( m 2 ) (LSC data). ALOS-1 PALSAR acquires data at an incidence angle of 41.5°. Data from three different ocean regions were acquired so that it was possible to characterize, in each case, the sea state based on the statistical properties of the underlying elementary scatterers. The three regions depicted in Figure 2, are from the central Atlantic, Atlantic equator, and north Atlantic. Sea state conditions can be found from the Copernicus Climate Data Store [43].
In Figure 3 that follows statistical measurements are given for the ocean region Atlantic Center. Specifically, the distribution of the closeness (distance) of those SAR pixels characterized as trihedral scatterers from the ideal trihedral is given in Figure 3a. The pixels close to ideal trihedral contribute to the distribution near zero. In Figure 3b is given the co-occurrence matrix for trihedral and cylinder scatterers, with the axes representing the closeness to the corresponding ideal scatterer. Finally, in Figure 3c the same information is provided in a 3D graph.

3.2. SLC Fully Polarimetric Data Preprocessing Stage

SLC fully polarimetric SAR images require special preprocessing [44] in order to convert the data into a form suitable for decomposition and feature extraction for Earth surface classification. The stages followed in this work were the following:
-
Radiometric Calibration
Calibration radiometrically corrects a SAR image so that the pixel values truly represent the radar backscatter of the reflecting surface [45]. The corrections that get applied during calibration are mission-specific; therefore, the software automatically determines what the kind of input product is and what corrections need to be applied based on the product’s metadata. Calibration is essential for quantitative use of SAR data. During calibration, it is important to save the output product in complex format since Cameron decomposition requires the complex inputs of the Sinclair matrix.
-
Polarimetric Matrix Generation
All the polarimetric tools work with either Coherency or Covariance matrices as input. Starting from a Quad Pol SLC product, the Matrix Generation operator converts the product into these types of matrices.
-
Polarimetric Speckle Filtering
SAR images are deteriorated by the speckle noise during image capture and transmission. It is a type of multiplicative noise which appears due to the variation in the backscatter from non-homogeneous cells which gives a grainy appearance to SAR images [46]. The backscatter is the result of constructive and destructive interference of the microwave radiation in resolution cells. Many different methods have been discussed in the literature to eliminate speckle noise. The error in the captured image is mainly related to geometry and intensity of pixels. These errors are reformed using definite or statistical mathematical models. The classification of images becomes difficult as the speckle noise manipulates the statistics of the backscatter. It has undesirable effects on the Automatic Recognition System (ATR); hence, speckle filtering has become the first essential step for SAR image processing, such as segmentation and edge detection. Multi-look processing and spatial filtering reduce the speckle noise from images. Multi-look processing is performed at the time of acquisition of low-resolution images. It averages out the speckle noise by taking several “looks” at a target in a single RADAR radius. The average is incoherent. The single-look processing is quite noisy. Spatial filtering is performed after the acquisition of data.
-
Polarimetric Decomposition
In this work, feature extraction utilizes Cameron decomposition. Therefore, each PolSAR cell now corresponds to one elemental scattering mechanism.
-
Geometric Correction utilizing Range Doppler Terrain Correction
This is correcting of geometric distortions that lead to geolocation errors. The distortions are induced by side-looking (rather than straight-down-looking or nadir) imaging and are compounded by rugged terrain.

3.3. Statistical Feature Extraction

The statistical features used in this work to characterize the distributions of the data used are the mean, the variance, the skewness, and the kurtosis. Second-order skewness and kurtosis are employed to characterize the shape of the co-occurrence matrices due to their 2D nature. For the statistical analysis of the 1D and 2D histograms, the relevant functions of the MATLAB are employed.
In probability theory and statistics [47,48], skewness is a measure of the asymmetry of the probability distribution of a real-valued random variable about its mean. The skewness value can be positive, zero, negative, or undefined. For a unimodal distribution, negative skewness commonly indicates that the tail is on the left side of the distribution, and positive skew indicates that the tail is on the right. In cases where one tail is long but the other tail is fat, skewness does not obey a simple rule. Accordingly:
Negative skew: The left tail is longer. The mass of the distribution is concentrated on the right. The distribution is said to be left-skewed, left-tailed, or skewed to the left, despite the fact that the curve itself appears to be skewed or leaning to the right. A left-skewed distribution usually appears as a right-leaning curve [47].
Positive skew: The right tail is longer. The mass of the distribution is concentrated on the left. The distribution is said to be right-skewed, right-tailed, or skewed to the right, despite the fact that the curve itself appears to be skewed or leaning to the left. A right-skewed distribution usually appears as a left-leaning curve [47].
Skewness in a data series may sometimes be observed not only graphically but by simple inspection of the values. For example, consider the numeric sequence (49, 50, 51), whose values are evenly distributed around a central value of 50. We can transform this sequence into a negatively skewed distribution by adding a value far below the mean, which is probably a negative outlier, e.g., (40, 49, 50, 51). Therefore, the mean of the sequence becomes 47.5, and the median is 49.5. Based on the formula of nonparametric skew, defined as ( μ ν ) / σ , the skew is negative. Similarly, we can make the sequence positively skewed by adding a value far above the mean, which is probably a positive outlier, e.g., (49, 50, 51, 60), where the mean is 52.5, and the median is 50.5.
As mentioned earlier, a unimodal distribution with zero value of skewness does not imply that this distribution is symmetric necessarily. However, a symmetric unimodal or multimodal distribution always has zero skewness.
The skewness [47,48] of a random variable X is the third standardized moment μ ~ 3 , defined as
μ ~ = Ε X μ σ 3 = μ 3 σ 3
where μ is the mean, σ is the standard deviation, E is the expectation operator, and μ 3 is the third central moment. It is sometimes referred to as Pearson’s moment coefficient of skewness, or simply the moment coefficient of skewness, but should not be confused with Pearson’s other skewness statistics.
If σ is finite, μ is finite too and skewness can be expressed in terms of the non-central moment E [ X 3 ] by expanding the previous formula:
μ ~ 3 = Ε X μ σ 3 = E X 3 3 μ σ 2 μ 3 σ 3
Skewness is a descriptive statistic that can be used in conjunction with the histogram and the normal quantile plot to characterize the data or distribution. Skewness indicates the direction and relative magnitude of a distribution’s deviation from the normal distribution. With pronounced skewness, standard statistical inference procedures such as a confidence interval for a mean will be not only incorrect, in the sense that the true coverage level will differ from the nominal (e.g., 95%) level, but will also result in unequal error probabilities on each side. Many models assume normal distribution, i.e., the data are symmetric about the mean. The normal distribution has a skewness of zero. But, in reality, data points may not be perfectly symmetric. So, an understanding of the skewness of the dataset indicates whether deviations from the mean are going to be positive or negative. D’Agostino’s K-squared test is a goodness-of-fit normality test based on sample skewness and sample kurtosis [48].
In probability theory and statistics, kurtosis is a measure of the “tailedness” of the probability distribution of a real-valued random variable. Like skewness, kurtosis describes a particular aspect of a probability distribution. There are different ways to quantify kurtosis for a theoretical distribution, and there are corresponding ways of estimating it using a sample from a population. Different measures of kurtosis may have different interpretations.
The standard measure of a distribution’s kurtosis, originating with Karl Pearson, is a scaled version of the fourth moment of the distribution. This number is related to the tails of the distribution, not its peak. For this measure, higher kurtosis corresponds to a greater extremity of deviations (or outliers) and not the configuration of data near the mean.
The kurtosis is the fourth standardized moment, defined as
K u r t X = E X μ σ 4 = μ 4 σ 4
where μ 4 is the fourth central moment and σ is the standard deviation.

4. Statistical Analysis of PolSAR Data and Experimental Results

The statistical analysis of the PolSAR data acquired from the sea surface is carried out based on four different features:
  • The percentage of the appearance of each dominating elementary scatterer in the specific sea region.
  • The statistical distribution of the closeness of the estimated dominating elementary scatterer on each separate pixel to the ideal (reference) elementary scatterers. The closeness is evaluated by means of the distance d z , z ref given in Equation (4) and analyzed in Section 2. Accordingly, in a specific region of the ocean, all pixels characterized by the same elementary scatterer contribute to the same distribution describing their distance from this ideal elementary scatterer.
  • The co-occurrence matrix for encountering the neighboring SAR pixels which correspond to the same elementary scatterer; the matrix axes reveal the closeness to this ideal common elementary scatterer.
  • The co-occurrence matrix for encountering the neighboring SAR pixels which correspond to different elementary scatterers; the matrix axes reveal the closeness to the corresponding ideal elementary scatterers.
An example is given in Figure 3a regarding the distribution of the closeness (distance) of those SAR pixels characterized as trihedral scatterers from the ideal trihedral. In Figure 3b is depicted the co-occurrence matrix for trihedral and cylinder scatterers where the axes represent the closeness to the corresponding ideal scatterer. In Figure 3c, the same information as in Figure 3b is provided in a 3D graph.
In the following, experimental measurements are carried out for the three regions in the Atlantic Ocean as well as for the region of Vancouver. Working on the region of Vancouver, measurements were carried out separately for sea and land. The sea near Vancouver is characterized by calm weather conditions compared to those in all three regions of Atlantic. In the following five Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, the statistical properties of the five regions are presented regarding the closeness of the more frequently met scatterers to the corresponding ideal ones. Additionally, in the same Figures is given information regarding the location and the time of the data acquisition. In Figure 8, where the statistics of the appearance of the elementary scatterers are shown for the land (mainly) of the Vancouver region, the appearance of the trihedral scatterer is quite small in comparison to that in flat areas like sea. This is obvious from the scatterers’ distribution given at the side of the PolSAR HV amplitude image of the region.
According to the material exposed so far and the data available from all sources, Table 3 was constructed. Commenting on the data in Table 3 as well as based on the previous Figures, we can reach the following remarks:
  • PolSAR pixels at sea are mainly represented by the trihedral elementary scatterer followed by the cylinder. The ¼-wave device is usually not negligible to around 5%.
  • For land cover, the percentage of trihedral scatterer is small (<20%) compared to that found on the sea surface (>50%). The percentages of the cylinder, the dipole, the narrow diplane, and the ¼-wave device are large compared to the percentages of the corresponding elementary scatterers on the sea surface.
  • In the case of calm sea state conditions, the percentage of the trihedral elementary scatterer becomes quite large (over 70%).
  • In the open sea, for moderate sea waves height, the percentage of cylinder elementary scatterers increases.
  • The closeness of each elementary scatterer to the corresponding ideal scatterer, evaluated for each separate PolSAR pixel, does not change significantly with the sea state conditions.
  • Co-occurrences for the same most frequent elementary scatterer, the trihedral, happens in the calm sea (Vancouver), most frequently between scatterers which are close to the ideal counterparts. Co-occurrences between the most frequent different elementary scatterers (trihedral and cylinder) show no significant difference between calm and rough sea (Figure 9).
In order to take advantage of other statistical features for characterizing the sea state of the four ocean regions tested, the four first statistical moments of the distributions presented in Figure 4, Figure 5, Figure 6 and Figure 7 were evaluated and are presented in Table 4. These distributions depict the closeness of the declared scatterer in a specific pixel to the corresponding ideal one and represent the most frequent scatterers i.e., the trihedral, the dipole, the cylinder and the ¼-wave devise. According to what is exposed in Table 4, we mention the following:
  • The mean for each of the four elementary scatterers does not change significantly from one sea region to another sea region, even if the sea state conditions change.
  • The variance for the trihedral is large mainly for rough sea state. The variance for the rest of the three elementary scatterers remains unchanged for various sea conditions.
  • The skewness and the kurtosis slightly increase for the trihedral in calm sea. For the dipole, both the third and fourth moments remain almost unchanged. The same happens for the cylinder and the ¼-wave device as well.
  • In general, skewness is of small values (near zero), so the distributions, especially for the dipole, are expected to be symmetric around the mean.
The 2D skewness and kurtosis of the co-occurrence matrices presented in Figure 9, for all tested ocean regions, are employed for characterizing the sea clutter in the four ocean regions. The co-occurrence matrices 1-1 and 4-4 correspond to transitions between the same scatterer (trihedral or cylinder) while 1-4 corresponds to transitions between trihedral and cylinder. The results given in Table 5 lead to the following conclusions:
  • For small wave heights, the skewness and kurtosis for 1-1 co-occurrence is relatively large.
  • For low wind speed, large values for skewness and kurtosis of the 2D co-occurrences are noted.
  • The 2D skewness and the kurtosis for the co-occurrence matrices in case of transitions from cylinder to trihedral and vice versa are relatively smaller for low wind speed and wave height.

5. Conclusions

The results demonstrated in this work prove that the sea state conditions affect the statistics as well as the second-order statistics of the PolSAR pixels obtained from satellites. Each pixel is represented by an elementary scattering mechanism obtained by Cameron’s CTD. A variety of features were employed for this purpose, starting from the simple percentage of the appearance of the elementary scatterers in a sea region, the closeness of each pixel-scatterer to its ideal counterpart, and the second-order statistics (co-occurrence) of each elementary scatterer with similar ones or with scatterers of a specific different kind. Specifically, the four first statistical moments were employed to characterize the distribution which reveals the closeness of each elementary scatterer from the ideal one. Moreover, since co-occurrence matrices of various kinds were evaluated, their 2D skewness and kurtosis were employed to characterize the transitions between various elementary scatterers. The results provided in Table 3, Table 4 and Table 5 prove that new features are now available for characterizing the sea state conditions. These features can be used as well for ship detection at sea since the statistical behavior of the PolSAR pixels of a ship are totally different from those of the sea surface.

Author Contributions

G.K. and V.A. have contributed to this work almost equally in all aspects of its conceptualization, realization, and presentation. All authors have read and agreed to the published version of the manuscript.

Funding

This research is partially funded by the Research Committee of the University of Patras.

Data Availability Statement

The data are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anastassopoulos, V.; Lampropoulos, G.A. A generalized compound model for radar clutter. In Proceedings of the IEEE National Radar Conference, Atlanta, GA, USA, 29–31 March 1994. [Google Scholar]
  2. Anastassopoulos, V.; Lampropoulos, G.A. High resolution radar clutter classification. In Proceedings of the IEEE International Radar Conference, Washington, DC, USA, 8–11 May 1995. [Google Scholar]
  3. Anastassopoulos, V.; Lampropoulos, G.A.; Drosopoulos, A.; Rey, M. High Resolution Radar Clutter Statistics. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 43–60. [Google Scholar] [CrossRef]
  4. Nohara, T.J.; Haykin, S. Canadian East Coast radar trials and the K-distribution. IEEE Process. Part F 1991, 138, 80–81. [Google Scholar] [CrossRef]
  5. Watts, S. Radar detection prediction in K-distributed sea clutter and thermal noise. IEEE Trans. Aerosp. Electron. Syst. 1987, 1, 40–45. [Google Scholar] [CrossRef]
  6. Trunk, G.V. Radar properties of Non-Rayleigh Sea clutter. IEEE Trans. Aerosp. Electron. Syst. 1972, 8, 196–204. [Google Scholar] [CrossRef]
  7. Ferreira, J.A.; Nascimento, A.D.C.; Frery, A.C. PolSAR Models with Multimodal Intensities. Remote Sens. 2022, 14, 5083. [Google Scholar] [CrossRef]
  8. Schuler, D.L.; Leea, J.S.; Kasilingamb, D.; Pottier, E. Measurement of ocean surface slopes and wave spectra using polarimetric SAR image data. Remote Sens. Environ. 2004, 91, 198–211. [Google Scholar] [CrossRef]
  9. He, Y.; Shen, H.; Perrie, W. Remote Sensing of Ocean Waves by Polarimetric SAR. J. Atmos. Ocean. Technol. 2006, 23, 1768–1773. [Google Scholar] [CrossRef]
  10. Salberg, A.B.; Larsen, S.O. Classification of Ocean Surface Slicks in Simulated Hybrid-Polarimetric SAR Data. IEEE Trans. Geosci. Remote Sens. 2018, 56, 7062–7073. [Google Scholar] [CrossRef]
  11. Yeremy, M.; Campbell, J.W.M.; Mattar, K.; Potter, T. Ocean Surveillance with Polarimetric SAR. Can. J. Remote Sens. 2001, 27, 328–344. [Google Scholar] [CrossRef]
  12. Zhang, B.; Perrie, W.; Vachon, P.W.; Li, X.; Pichel, W.G.; Guo, J.; He, Y. Ocean Vector Winds Retrieval From C-Band Fully Polarimetric SAR Measurements. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4252–4261. [Google Scholar] [CrossRef]
  13. Ibrahim, Z.; Bahnacy, A.I.; Hussein, K.F.A.; Farahat, A.E. PolSAR Target Decomposition Technique to Estimate Ocean Surface temperature and Salinity. In Proceedings of the 39th National Radio Science Conference (NRSC 2022), Cairo, Egypt, 29 November–1 December 2022. [Google Scholar]
  14. Viana, R.D.; Lorenzzetti, J.A.; Carvalho, J.T.; Nunziata, F. Estimating Energy Dissipation Rate from Breaking Waves Using Polarimetric SAR Images. Sensors 2020, 20, 6540. [Google Scholar] [CrossRef] [PubMed]
  15. Huang, L.; Fan, C.; Meng, J.; Yang, J.; Zhang, J. Shallow sea topography detection using fully polarimetric Gaofen-3 SAR data based on swell patterns. Acta Oceanol. Sin. 2022, 41, 150–162. [Google Scholar] [CrossRef]
  16. Xie, T.; Ouyang, R.; Perrie, W.; Zhao, L.; Zhang, X. Proof and Application of Discriminating Ocean Oil Spills and Seawater Based on Polarization Ratio Using Quad-Polarization Synthetic Aperture Radar. Remote Sens. 2023, 15, 1855. [Google Scholar] [CrossRef]
  17. Kumar, A.; Mishra, V.; Panigrahi, R.K.; Martorella, M. Application of Hybrid-Pol SAR in Oil-Spill Detection. IEEE Geosci. Remote Sens. Lett. 2023, 20, 1–5. [Google Scholar] [CrossRef]
  18. Li, H.; Perrie, W.; Wu, J. Retrieval of Oil–Water Mixture Ratio at Ocean Surface Using Compact Polarimetry Synthetic Aperture Radar. Remote Sens. 2019, 11, 816. [Google Scholar] [CrossRef]
  19. Kreithen, D.E.; Halversen, S.S.; Owirka, G.J. Discriminating targets from clutter. Linc. Lab. J. 1993, 6, 25–52. [Google Scholar]
  20. Chaney, R.D.; Burl, M.C.; Novak, L.M. On the performance of polarimetric target detection algorithms. In Proceedings of the IEEE 1990 International Radar Conference, Arlington, VA, USA, 7–10 May 1990; pp. 520–525. [Google Scholar]
  21. Pottier, E.; Saillard, J. Optimal polarimetric detection of radar target in a slowly fluctuating environment of clutter. In Proceedings of the IEEE 1990 International Radar Conference, Arlington, VA, USA, 7–10 May 1990; pp. 211–216. [Google Scholar]
  22. Novak, L.M.; Burl, M.C.; Chaney, R.D.; Owirka, G.J. Optimal processing of polarimetric synthetic-aperture radar imagery. Linc. Lab. J. 1990, 3, 273–290. [Google Scholar]
  23. Gu, D.; Xu, X. Multi-Feature Extraction of Ships from SAR Images. Image and Signal Processing (CISP). In Proceedings of the 6th International Congress, Hangzhou, China, 16–18 December 2013; Volume 1, pp. 454–458. [Google Scholar]
  24. Anastassopoulos, V.; Lampropoulos, G.A. Optimal CFAR Detection in Weibull Clutter. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 52–64. [Google Scholar] [CrossRef]
  25. Fei, C.; Liu, T.; Lampropoulos, G.A.; Anastassopoulos, V. Markov Chain CFAR Detection for Polarimetric Data Using Data Fusion. IEEE Trans. Geosci. Remote Sens. 2012, 50, 397–408. [Google Scholar] [CrossRef]
  26. Lee, J.S.; Pottier, E. Polarimetric Radar Radar Imaging, 1st ed.; CRC Press: New York, NY, USA, 2009. [Google Scholar]
  27. Nunziata, F.; Migliaccio, M.; Brown, C.E. Reflection Symmetry for Polarimetric Observation of Man-Made Metallic Targets at Sea. IEEE J. Ocean. Eng. 2012, 37, 384–394. [Google Scholar] [CrossRef]
  28. Addabbo, P.; Biondi, F.; Clemente, C.; Orlando, D.; Pallotta, L. Classification of covariance matrix eigenvalues in polarimetric SAR for environmental monitoring applications. IEEE Trans. Aerosp. Electron. Syst. Mag. 2019, 34, 28–43. [Google Scholar] [CrossRef]
  29. Pallotta, L.; Clemente, C.; De Maio, A.; Soraghan, J.J. Detecting Covariance Symmetries in Polarimetric SAR Images. IEEE Trans. Geosci. Remote Sens. 2017, 55, 80–95. [Google Scholar] [CrossRef]
  30. Van Zyl, J.J.; Arii, M.; Kim, Y. Model-based decomposition of polarimetric SAR covariance matrices constrained for nonnegative eigenvalues. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3452–3459. [Google Scholar] [CrossRef]
  31. Touzi, R.; Charbonneau, F.; Hawkins, R.; Murnaghan, K.; Kavoun, X. Ship–sea contrast optimization when using polarimetric SARs. In Proceedings of the Geoscience and Remote Sensing Symposium, IGARSS ‘01, Sydney, NSW, Australia, 9–13 July 2001; Volume 1, pp. 426–428. [Google Scholar]
  32. Touzi, R.; Charbonneau, F.; Hawkins, R.; Vachon, P. Ship detection and characterization using polarimetric SAR. Can. J. Remote Sens. 2004, 30, 552–559. [Google Scholar] [CrossRef]
  33. Ringrose, R.; Harris, N. Ship detection using polarimetric SAR data. In Proceedings of the SAR Workshop: CEOS Committee on Earth Observation Satellites, Toulouse, France, 26–29 October 1999; pp. 687–691. [Google Scholar]
  34. Cameron, W.L.; Youssef, N.N.; Leung, L.K. Simulated polarimetric signatures of primitive geometrical shapes. IEEE Trans. Geosci. Remote Sens. 1996, 34, 793–803. [Google Scholar] [CrossRef]
  35. Kouroupis, G.; Anastassopoulos, V. Scatterer Characterization Based on the Condiagonalization of the Sinclair Backscattering Matrix. Prog. Electromagn. Res. M 2019, 85, 59–69. [Google Scholar] [CrossRef]
  36. Kouroupis, G.; Anastassopoulos, V. A polarimetric CFAR ship detector based on the joint probability function of simulated first-order Markov chains. Int. J. Remote Sens. 2019, 40, 5121–5140. [Google Scholar] [CrossRef]
  37. Therrien, C.W. Random Processes; Prentice Hall Signal Processing Series; Prentice-Hall: Hoboken, NJ, USA, 1992. [Google Scholar]
  38. Kouroupis, G.; Anastassopoulos, V. A Markov Chain Model Based on Cameron’s CTD Ship Detection Scheme. In Proceedings of the IEEE Imaging Systems and Techniques, Chania, Greece, 4–6 October 2016; pp. 100–105. [Google Scholar]
  39. SNAP Toolboxes. Science Toolbox Exploitation Platform. Available online: https://step.esa.int/main/doc/ (accessed on 1 December 2022).
  40. Cameron, W.L.; Rais, H. Conservative polarimetric scatterers and their role in incorrect extensions of the Cameron decomposition. IEEE Trans. Geosci. Remote Sens. 2006, 44, 3506–3516. [Google Scholar] [CrossRef]
  41. Sarker, S.; Veremyev, A.; Boginski, V.; Arvind Singh, A. Critical nodes in river networks. Sci. Rep. Nat. Res. 2019, 9, 11178. [Google Scholar] [CrossRef]
  42. MDA, RADARSAT-2 Product Description, RN-SP-52-1238 Issue 1/14. Available online: https://earth.esa.int/eogateway/documents/20142/0/Radarsat-2-Product-description.pdf/f2783c7b-6a22-cbe4-f4c1-6992f9926dca (accessed on 1 December 2022).
  43. Copernicus Climate Data Store. Available online: https://cds.climate.copernicus.eu/cdsapp#!/home (accessed on 1 January 2023).
  44. Kaushik, P.; Jabin, S. A Comparative study of Pre-processing Techniques of SAR Images. In Proceedings of the 4th International Conference on Computing Communication and Automation (ICCCA), Greater Noida, India, 14–15 December 2018; pp. 1–4. [Google Scholar]
  45. Veci, L.; Lu, J.; Prats-Iraola, P.; Scheiber, R.; Collard, F.; Norman, F.; Engdahl, M. The Sentinel-1 Toolbox. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Quebec, QC, Canada, 13–18 July 2014; pp. 1–3. [Google Scholar]
  46. Lee, J.S.; Jurkevich, L.; Dewaele, P.; Wambacq, P.; Oosterlinck, A. Speckle filtering of synthetic aperture radar images: A review. Remote Sens. Lett. 2009, 8, 313–340. [Google Scholar] [CrossRef]
  47. Joanes, D.N.; Gill, C.A. Comparing measures of sample skewness and kurtosis. J. R. Stat. Soc. 1998, 47, 183–189. [Google Scholar] [CrossRef]
  48. Westfall, P.H. Kurtosis as Peakedness, 1905–2014. R.I.P. Am. Stat. 2014, 68, 191–195. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) The broader area of Vancouver, British Columbia, Canada. The geocoded data were obtained by means of the SNAP software. The image corresponds to the amplitude of the I and Q components of the HV part of the fully polarimetric data. (b) The geocoded SAR image registered using SNAP software on the Google Earth map of the broader area of Vancouver.
Figure 1. (a) The broader area of Vancouver, British Columbia, Canada. The geocoded data were obtained by means of the SNAP software. The image corresponds to the amplitude of the I and Q components of the HV part of the fully polarimetric data. (b) The geocoded SAR image registered using SNAP software on the Google Earth map of the broader area of Vancouver.
Remotesensing 15 02837 g001
Figure 2. PALSAR ALOS-1 data from the Atlantic Ocean: (a) center, (b) equator, and (c) north.
Figure 2. PALSAR ALOS-1 data from the Atlantic Ocean: (a) center, (b) equator, and (c) north.
Remotesensing 15 02837 g002
Figure 3. Ocean Region: Atlantic Center. (a) The distribution of the closeness (distance) of those SAR pixels characterized as trihedral scatterers from the ideal trihedral. The pixels close to ideal trihedral contribute to the distribution near zero. (b) Co-occurrence matrix for trihedral and cylinder scatterers, with the axes representing the closeness to the corresponding ideal scatterer. (c) The same information is provided in a 3D graph.
Figure 3. Ocean Region: Atlantic Center. (a) The distribution of the closeness (distance) of those SAR pixels characterized as trihedral scatterers from the ideal trihedral. The pixels close to ideal trihedral contribute to the distribution near zero. (b) Co-occurrence matrix for trihedral and cylinder scatterers, with the axes representing the closeness to the corresponding ideal scatterer. (c) The same information is provided in a 3D graph.
Remotesensing 15 02837 g003
Figure 4. Region: Atlantic Center, Latitude −5.75, Longitude −14.75, Date-Time 15 January 2011 23:00. (a) The location on the Earth’s surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Figure 4. Region: Atlantic Center, Latitude −5.75, Longitude −14.75, Date-Time 15 January 2011 23:00. (a) The location on the Earth’s surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Remotesensing 15 02837 g004
Figure 5. Region: Atlantic Equator, Latitude −13.25, Longitude −55, Date-Time 2 December 2010 02:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Figure 5. Region: Atlantic Equator, Latitude −13.25, Longitude −55, Date-Time 2 December 2010 02:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Remotesensing 15 02837 g005
Figure 6. Region: Atlantic North, Latitude 68.25, Longitude −6.25, Date-Time 17 November 2009 11:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Figure 6. Region: Atlantic North, Latitude 68.25, Longitude −6.25, Date-Time 17 November 2009 11:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Remotesensing 15 02837 g006aRemotesensing 15 02837 g006b
Figure 7. Region: Vancouver Sea, Latitude 49, Longitude –123.25, Date-Time 15 April 2008 15:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude (the four squared region in the sea were used), (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Figure 7. Region: Vancouver Sea, Latitude 49, Longitude –123.25, Date-Time 15 April 2008 15:00. (a) The location on the Earth surface, (b) image of the sea surface in PolSAR HV amplitude (the four squared region in the sea were used), (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder, (e) the distribution of the pixels that respond as dipole elementary scatterers based on the closeness to ideal dipole, and (f) the distribution of the pixels that respond as ¼-wave device elementary scatterers based on the closeness to ideal ¼-wave device.
Remotesensing 15 02837 g007
Figure 8. Region: Vancouver Land, Latitude 49, Longitude –123.25, Date-Time 15 April 2008 15:00. (a) Image of the land surface in PolSAR HV amplitude, (b) the distribution of the elementary scatterers in the area of Vancouver, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder.
Figure 8. Region: Vancouver Land, Latitude 49, Longitude –123.25, Date-Time 15 April 2008 15:00. (a) Image of the land surface in PolSAR HV amplitude, (b) the distribution of the elementary scatterers in the area of Vancouver, (c) the distribution of the pixels that respond as trihedral elementary scatterers based on the closeness to ideal trihedral, (d) the distribution of the pixels that respond as cylinder elementary scatterers based on the closeness to ideal cylinder.
Remotesensing 15 02837 g008
Figure 9. Co-occurrence matrices for two different sea regions. The upper row comes from Atlantic North while the second row is from the calm Vancouver sea. Co-occurrences for the same most frequent elementary scatterer, the trihedral, happens in the calm sea (Vancouver) most frequently between scatterers which are close to the ideal counterparts (c). Co-occurrences between the most frequent different elementary scatterers (trihedral and cylinder) show no significant difference between calm (d) and rough sea (b).
Figure 9. Co-occurrence matrices for two different sea regions. The upper row comes from Atlantic North while the second row is from the calm Vancouver sea. Co-occurrences for the same most frequent elementary scatterer, the trihedral, happens in the calm sea (Vancouver) most frequently between scatterers which are close to the ideal counterparts (c). Co-occurrences between the most frequent different elementary scatterers (trihedral and cylinder) show no significant difference between calm (d) and rough sea (b).
Remotesensing 15 02837 g009aRemotesensing 15 02837 g009b
Table 1. Complex parameter z corresponding to elementary scattering mechanisms.
Table 1. Complex parameter z corresponding to elementary scattering mechanisms.
Complex Parameter z Normalized   Complex   Vector   Λ ^ ( z ) Scattering Mechanism
1 Λ ^ ( 1 ) Trihedral
−1 Λ ^ ( 1 ) Diplane
0 Λ ^ ( 0 ) Dipole
+1/2 Λ ^ ( + 1 / 2 ) Cylinder
−1/2 Λ ^ ( 1 / 2 ) Narrow diplane
± j Λ ^ ( ± j ) ¼-wave device
Table 2. Cameron 8 Elemental Scatterers color-coding. JET color map is from MATLAB.
Table 2. Cameron 8 Elemental Scatterers color-coding. JET color map is from MATLAB.
Symmetric Elementary ScattererClassCameron Color Representation
Trihedral1Remotesensing 15 02837 i001
Diplane2
Dipole3
Cylinder4
Narrow Diplane5
1/4 Wave Device6
Left Helix7
Right Helix8
Table 3. The percentage of the appearance of various elementary scatterers in five different regions along with the Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws).
Table 3. The percentage of the appearance of various elementary scatterers in five different regions along with the Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws).
RegionTrihedralDiplaneDipoleCylinderNarrow Diplane1/4 Wave DeviceLeft HelixRight Helixswhws
Atlantic center51.20.22.9390.65.30.30.31.735
Atlantic Equator66.90.22.524.60.750.010.012.346.6
Atlantic North59.50.22.532.20.64.60.20.22.936.51
Vancouver Sea73.60.11.720.50.43.50.10.10.193.7
Vancouver Sea + Land19.2315.225.682332.5 0.193.7
Table 4. The four first moments of the distributions presented in Figure 4, Figure 5, Figure 6 and Figure 7. These distributions depict the closeness of the declared scatterer to the corresponding ideal one and correspond to the most frequent scatterers. The moments are presented for the four different sea regions along with the Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws).
Table 4. The four first moments of the distributions presented in Figure 4, Figure 5, Figure 6 and Figure 7. These distributions depict the closeness of the declared scatterer to the corresponding ideal one and correspond to the most frequent scatterers. The moments are presented for the four different sea regions along with the Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws).
RegionMomentTrihedralDipoleCylinder1/4 Wave Deviceswhws
Atlantic centerMean8.8913.269.4925.931.735
Variance15.2226.9020.12169.7
Skewness1.02−0.080.720.43
Kurtosis4.102.333.302.0
Atlantic EquatorMean7.4713.3210.2325.762.346.6
Variance17.1127.3322.47169.02
Skewness1.05−0.080.570.44
Kurtosis4.072.32.862.03
Atlantic NorthMean8.2613.39.6625.852.936.51
Variance14.826.7120.18168.22
Skewness1.08−0.090.70.44
Kurtosis4.362.343.252.03
Vancouver SeaMean7.0513.3510.125.840.193.7
Variance14.8327.9921.35170.09
Skewness1.16−0.10.620.44
Kurtosis4.582.263.022.01
Table 5. Two-dimensional skewness and kurtosis of the co-occurrence matrices presented in Figure 9, for all tested ocean regions. Co-occurrence matrices 1-1 and 4-4 correspond to transitions between the same scatterer (trihedral or cylinder) while 1-4 correspond to transitions between trihedral and cylinder. The Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws) are also given for the corresponding ocean locations.
Table 5. Two-dimensional skewness and kurtosis of the co-occurrence matrices presented in Figure 9, for all tested ocean regions. Co-occurrence matrices 1-1 and 4-4 correspond to transitions between the same scatterer (trihedral or cylinder) while 1-4 correspond to transitions between trihedral and cylinder. The Significant Height of combined wind Waves and swell in meters (swh) and the horizontal Wind Speed in m/s (ws) are also given for the corresponding ocean locations.
RegionCo-Occurrence 1-1
Skewness
Co-Occurrence 1-1
Kurtosis
Co-Occurrence 4-4
Skewness
Co-Occurrence 4-4
Kurtosis
Co-Occurrence 1-4
Skewness
Co-Occurrence 1-4
Kurtosis
swhws
Atlantic center3.02317.7562.73311.892.86212.2271.735
Atlantic Equator2.2457.2382.2839.4282.2398.0712.346.6
Atlantic North3.01911.9272.87913.0563.15214.6112.936.51
Vancouver Sea3.05512.3172.7112.3562.82912.2730.193.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Koukiou, G.; Anastassopoulos, V. Ocean Clutter Characterization Based on PolSAR Data and Second-Order Statistics of Elementary Scatterers. Remote Sens. 2023, 15, 2837. https://doi.org/10.3390/rs15112837

AMA Style

Koukiou G, Anastassopoulos V. Ocean Clutter Characterization Based on PolSAR Data and Second-Order Statistics of Elementary Scatterers. Remote Sensing. 2023; 15(11):2837. https://doi.org/10.3390/rs15112837

Chicago/Turabian Style

Koukiou, Georgia, and Vassilis Anastassopoulos. 2023. "Ocean Clutter Characterization Based on PolSAR Data and Second-Order Statistics of Elementary Scatterers" Remote Sensing 15, no. 11: 2837. https://doi.org/10.3390/rs15112837

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

Article Metrics

Back to TopTop