Oil Spill Detection with Dual-Polarimetric Sentinel-1 SAR Using Superpixel-Level Image Stretching and Deep Convolutional Neural Network

: Synthetic aperture radar (SAR) has been widely applied in oil spill detection on the sea surface due to the advantages of wide area coverage, all-weather operation, and multi-polarization characteristics. Sentinel-1 satellites can provide dual-polarized SAR data, and they have high potential for successful application to oil spill detection. However, the characteristics of the sea surface and oil film on different images are not the same when imaging at different locations and in different conditions, which leads to the inconsistent accuracy of these images with the application of the current oil spill detection methods. In order to avoid the above limitation, we propose an oil spill detection method using image stretching based on superpixels and a convolutional neural network. Experiments were carried out on eight Sentinel-1 dual-pol data, and the optimal superpixel number and image stretching parameters are discussed. Mean intersection over union (MIoU) was used to evaluate classification accuracy. The proposed method could effectively improve the classification accuracy; when the expansion and inhibition coefficients of image stretching were set to 1.6 and 1.2 respectively, the experiments achieved a maximum MIoU of 85.4%, 7.3% higher than that without image stretching.


Introduction
Crude oil leakage is one of the most influential marine pollution forms due to its difficult elimination.On 20 April 2010, an explosion accident at the drilling platform of British Petroleum occurred in the Gulf of Mexico.The sunken platform leaked about 5000 barrels of crude oil every day, representing one of the most serious ecological disasters in the history of the United States [1].The rapid and accurate detection of oil spill is of great significance to protect marine environments and reduce economic loss.
Synthetic aperture radar (SAR) can work in all weather conditions and all day; thus, it has been widely applied in oil spill detection.The SAR sensor emits electromagnetic waves and receives echoes.The oil slick floating on the sea surface will suppress capillary gravity waves and short gravity waves, which affect the Bragg scattering of microwaves and reduce the normalized radar backscatter cross-section (NRCS), representing the intensity of radar backscattering [2].Therefore, oil films are usually observed as dark spots on SAR images.However, many other factors also appear as dark spots in SAR images, and they could be confused with oil spill areas.These include low wind speed, biogenic oil film, plankton, waves, and rain cells, which can result in false detection.Thus, the challenge of oil spill detection is to distinguish oil spill areas from these look-alike areas [3].
Research on oil spill detection using SAR sensors began in the 1990s.Early oil spill detection is usually based on texture analysis of single polarimetric SAR images.Solberg et al. summarized oil spill detection in three steps: (1) dark spots recognition, (2) feature extraction, and (3) oil spill detection [4].They extracted 12 texture features from the dark spots obtained from ENVISAT, Radarsat, and ERS SAR data, and the oil spill areas were successfully separated from other similar areas using a Bayesian classifier [4,5].In 2013, Pineda et al. proposed a texture classification neural network (TCCNA) and achieved high accuracy in detecting the Gulf of Mexico oil spill [6].Cheng et al. [7] used VV channel data acquired from COSMO-SkyMed to perform oil spill monitoring and tracking through model simulations.In 2014, Diego et al. [8] developed a tool that could offer an integrated framework for the detection and localization of marine spills with fuzzy clustering and wavelets.The PRIMI pilot project of Italy [9] developed a set of oil spill detection algorithms using texture analysis and deep learning, which achieved good monitoring results.Ajadi et al. [10,11] introduced Lipschitz regularity and multiscale change detection techniques into oil spill detection.These methods mainly focused on the extraction of geometric and texture features of dark spots in SAR images.
The main shortcoming of oil spill detection based on single-pol SAR images is that there is little chance to distinguish look-alike areas from oil spills.By contrast, the polarized information of SAR images can avoid this problem.Several polarization decomposition methods have been applied in oil spill detection with quad-or dual-polarimetric SAR data, such as Cloude decomposition and single-bounce eigenvalue relative difference (SERD) [12][13][14][15].Skrunes et al. [16] proved that the geometric intensity and the real part of the co-polarization cross-product of C-band SAR data could distinguish well between biogenic slicks and mineral oil types.Li et al. analyzed the polarimetric scattering properties of oil-covered water using ellipticity parameters and the degree of polarization.They revealed that the main scattering mechanism of oil seep is even bounce scattering, whereas oil spill is dominated by Bragg scattering [17].Salberg et al. [18,19] evaluated several polarized parameters and classifiers with the use of hybrid polarimetric SAR for oil spill detection.Ivonin et al. [20] generalized C-band Radarsat-2 and TerraSAR-X data using resonant to non-resonant signal damping (RND), and the method could separate mineral oil slicks from other oil slicks.Espeseth et al. [21] examined system additive and multiplicative noise on the X-, C-, and L-band SAR data covering oil slicks.The results proved that the SNR from TerraSAR-X and Radarsat-2 below 0 dB was not suitable for analysis of the scattering properties.In 2021, Li et al. [22] designed improved polarimetric feature parameters on the basis of polarimetric scattering entropy and anisotropy, and the proposed features could indicate the oil slick information and effectively suppress sea clutter and look-alike information.
With the rise of machine learning in recent years, deep learning and computer vision algorithms have gradually been applied to polarimetric SAR oil spill detection.Oil spill detection is essentially a classification problem.Accordingly, current studies on oil spill detection based on deep learning methods typically divide SAR images into fixed sizes, and then apply a convolutional neural network (CNN) model or derived semantic segmentation for classification.Many classical model structures, including fully convolutional network (FCN), U-Net, PSPNet, DeepLab series, and object-contextual representations (ORC) [23][24][25][26][27][28], are widely applied in image classification.On this basis, in oil spill detection, Zeng et al. [29] proposed a classification network named OSCNet based on vgg-16, improving the accuracy to 92.5% on specific datasets.In 2018, Yu et al. [30] introduced the f-divergence loss function into a generative adversarial network model, and it accurately segmented different shapes of oil spills in SAR images.Bianchi et al. [31] designed a deep learning framework using dense blocks to distinguish oil and non-oil areas.Zhang et al. [32] proposed a method using an advanced convolutional neural network based on the superpixel model applied to quad-polarimetric SAR images, and it led to the highest MIoU of 90.5%.Since the size and area of oil spill dark spots are different, multiscale learning methods have attracted more attention in recent years [33,34], which use multiple convolution kernels with different receptive fields to process input images for feature extraction.Zhu et al. [35] proposed the oil spill contextual and boundary-supervised detection network (CBD-Net) to extract refined oil spill regions by fusing multiscale features, and the spatial and channel squeeze excitation (scSE) block was also introduced to further improve the accuracy.
The Sentinel-1 satellite carries C-band SAR sensors.It provides single-or dual-polarimetric SAR data with large swath width.In recent years, it was widely used, and many achievements were made in geohazard detection [36,37], infrastructure deformation monitoring [38][39][40], etc.In the field of oil spill detection, Setiani et al. [41] analyzed the Balikpapan oil spill accidents that happened in 2018.Using several types of index, oil spill areas were separated from clean sea.Rajendran et al. used the parallelepiped supervised algorithm and confusion matrix to detect oil spill patches from the Wakashio oil spill accidents [42].Naz et al. [43] analyzed four oil spill events over the Indian Ocean using Sentinel-1 data, and the spread of oil spills could be observed.Moreover, Diego et al. [44] proved that the convolutional neural network achieved higher accuracy with lower computing time, compared with clustering and logistic regression algorithms.Chaturvedi et al. [45] detected oil slicks near the Al Khafji zone using Sentinel-1 C-band SAR data.Corcione et al. [46] found that the largest deviation was observed over oil slicks using the azimuth auto correlation function (AACF).Krestenitis et al. [47] proposed a benchmark based on deep neural networks for future oil spill detection, and they found that DeepLabv3+ presented the best performance in experiments.Ma et al. [48] combined the polarized information of Sentinel-1 dual-pol SAR data with the deep learning algorithm, took the pseudo color image calculated from the polarized information as the input of the neural network, and improved the DeepLabV3+ model.The final oil spill detection accuracy rate was increased to 98.92%.In 2022, Dong et al. [49] analyzed 563,705 Sentinel-1 images from 2014 to 2019.They found that the proportion of anthropogenic discharges was an order-ofmagnitude greater than natural seepages by quantifying oil slick area, and they provided a detailed global inventory of static and persistent sources.
However, the final accuracy of oil spill detection algorithms based on deep learning is highly dependent on the consistency of the training and test sets.In real applications, there are imaging differences when using different SAR data due to different imaging time and locations, leading to the inconsistent accuracy of oil spill detection among different data, which is usually much lower than normal accuracy [31,47].In order to further improve accuracy and address the above problem, we propose an image stretching method based on SLIC superpixels and a convolutional neural network, which can enlarge the features of oil spill areas, making it possible to achieve better detection accuracy in datasets with imaging differences.Moreover, we designed a semantic segmentation model for classification.Experiments were carried out on eight Sentinel-1 dual-pol data.Two polarization decomposition methods were applied, namely, H/Alpha decomposition and Stokes parameters.The experimental results showed that the maximum MIoU achieved was 85.4% when using polarized parameters with the image stretching method, representing a 7.3% improvement in accuracy compared with the method without image stretching.The proposed method could significantly improve the accuracy of oil spill detection and alleviate the inconsistent accuracy caused by imaging differences.

Materials and Methods
The flowchart of the proposed method is illustrated in Figure 1.The SAR data were firstly processed using the refined Lee filter, and the polarimetric parameters were calculated.The Pauli pseudo-color image was applied to calculate the SLIC superpixels.Then, the image stretching method was applied to process the SLIC superpixel results.Subsequently, the polarimetric parameters and the image stretching results were used as inputs to the constructed neural network model for training and testing.Finally, the classification results were extracted as the outputs of the model.

Polarization Decomposition
The Sentinel-1 satellite provides two types of dual-pol SAR data: VV/VH and HH/HV.This paper used VV/VH data to extract two groups of polarized parameters: H/Alpha decomposition and Stokes parameters, as shown in Figure 2.For VV/VH data, the scattering matrix of polarimetric SAR represents the statistics of actual values of multiple scattering centers.The covariance matrix C and coherence matrix T can be used to measure the statistical characteristics and calculate polarized parameters for more accurate data analysis.The scattering matrix can be vectorized as follows: K2L and K2P are treated by the inner product with their conjugate transpose matrix, and the covariance matrix C and coherence matrix T can be expressed as By calculating eigenvalues λ1 and λ2 of matrix T2, the components of each eigenvector can be defined as Therefore, entropy (H) is defined as The mean scattering angle (α) is calculated by H, α, λ1, and λ2 are the four components in H/Alpha decomposition, used as H/Alpha/λ1/λ2 parameters.
The Stokes vector is a decomposition method for describing partially polarized waves [50,51].The polarization state of plane waves is usually represented by wave parameter sets (g0, g1, g2, g3) as follows: where Eh and Ev represent the electromagnetic vector components in the H and V directions, respectively; it can be seen from the formula that g0 2 = g1 2 + g2 2 + g3 2 .Physically, g0 is equal to the total power of plane wave and is proportional to the total amplitude, g1 represents the amplitude difference between H and V direction components, and g2 and g3 are the phase difference between them.The plane wave contains a polarization component if any parameter is not equal to 0. This paper defined a new parameter contrast = g1/g0, with contrast/g0/g1/g2/g3 parameters set as the group of Stokes parameters.

Neural Network Model
The construction of the convolutional neural network incorporated Resblock.The idea of Resblock is to add a part to the basic neural network structure, so that it can learn the residual information between the input and output of the convolution layer [52].Assuming that the convolution layer is represented by H and the characteristic layer of the input convolution layer is x, and the output of the convolution layer is denoted by H(x); hence, the operation process of convolution layer is to establish the mapping from X to H(x).The residual of the convolution layer is defined as ResNet adopts Resblock to realize residual learning.Figure 3 shows the Resblock structure used in this paper, where the input tensor was calculated using two convolutional layers and a depthwise separable convolutional layer, and the calculated tensor was finally connected with the input tensor.The basic structure of the semantic segmentation model is shown in Figure 4.The input training data are first processed by a feature extraction step.The extracted feature maps pass through four different scales of dilated convolution layers, before being further calculated by multiple Resblocks to obtain highly extracted feature maps.These feature maps are inversed to the same size as the inputs via several transposed convolution layers, which are connected to the above Resblocks with the same size.Figure 5 shows the detailed information of the data input and output of this model when multi-input training data are used.Different input images are processed by the feature extraction (FE) module, which includes two convolution layers and a ResBlock to generalize feature maps.At the end, feature maps extracted by the FE module are processed by 1 × 1 convolution and stacked with the output of the transposed convolution layers to generate classification results.Table 1 shows the detailed parameters in this model.

SLIC Superpixel Segmentation and Image Stretching
SAR images are affected by speckle noises, and the weakening of the edge of the oil spill area can also have a great impact on the accuracy.Although the convolutional kernel of the CNN can calculate input data in an area of several pixels, it is difficult to consider the context information of the image in a larger range.In order to solve this problem, the simple linear iteration cluster (SLIC) superpixel [53] was introduced for oil spill detection, which divides adjacent pixels with similar characteristics into a region to form an irregular superpixel block.
The calculation of SLIC superpixel segmentation depends on the superpixel number, and each superpixel block contains approximately the same number of pixels.Firstly, the clustering centers of all superpixel blocks are evenly distributed across the image with S pixels as intervals, and the image is divided into k regular grids.The formula of interval S is as follows, where N represents the total number of pixels: For each pixel i to be classified, the measurement distance D' is defined as the distance from the pixel to the cluster center, expressed as where dc is usually called the pixel distance, l, a, and b represent the values of pixels in three different layers of the CIE lab color space, ds is the spatial distance, (xj,yj) and (xi,yi) are the coordinates of the pixel and cluster center, and Nc and Ns are constants.The distance between points and all cluster centers in surrounding 2S × 2S pixels is calculated, and every point is classified to the center with the smallest distance.The position of the cluster center is moved to the midpoint of the superpixel block after all points are classified, and then the above steps are repeated for several iterations.
Superpixel segmentation was calculated on the Pauli pseudo-color image in this paper.The Pauli image is an RGB color image containing three channels based on matrix C2, where the R channel is C11, B channel is C22, and the G channel is defined as 2 Re( ) The RGB color image is converted into the CIE lab space before superpixel segmentation.
A threshold calculation method based on the SLIC superpixel is proposed to separate the oil spill area from other parts of the image.The experiments average the pixels in each superpixel block, use the average value to represent the superpixel block, and count the average, median, maximum, and minimum value of all superpixel blocks in the image.The threshold formula can be summarized as follows: where mean represents the average value, min represents the minimum value, and mid is the median value.The logmeanmid function can prevent the deviation of the threshold calculation caused by different regions in different images.After calculating the numerical threshold, the pixel threshold at each superpixel block is calculated as follows: If the pixel value of the block is higher than the threshold, it is multiplied by the expansion coefficient.On the contrary, if the pixel value is lower, it is divided by the inhibition coefficient.Through this method, the characteristics of the oil spill area in the SAR image can be highlighted, and the imaging differences between different images can be alleviated.The process pseudo-code of image stretching is presented in Algorithm 1.  end for 25: end for

Experiments and Results
This section describes the experiments and results for the classification of polarized parameters, SLIC superpixels, and image stretching, as well as the comparative experiments with and without image stretching, and the selection of polarized parameters.

Dataset
The SAR dataset was composed of Sentinel-1 dual-pol SAR data from the Persian Gulf region.The Persian Gulf is located between the Iranian Plateau and the Arabian Peninsula.Pollution is steadily increasing due to intensive oil production, serious tanker transportation, pipeline leakage, ship navigation, and maritime accidents in this area.According to Scanex [54], an average of 100,000-160,000 tons of oil and petroleum products flow into the Persian Gulf every year, and the degree of oil pollution is 47 times higher than the world average.In this section, eight VV/VH oil spill SAR data from Sentinel-1A in IW-SLC (interferometric wide swath, single-look complex) mode were selected for experiments.All images were numbered 1-8 in chronological order.Figure 6 presents the pseudo-color images of these eight images, where the dark spots within red dotted boxes are oil spill areas, while the blue dotted boxes represent look-alike areas.Oil spill areas were verified on the basis of in situ investigations and recorded by Scanex.Look-alike areas were present as dark spots on the images due to environmental or biological factors, as can be seen in Images 3-6 (see Figure 6c-f).The look-alike area marked in Image 6 (see Figure 6f) was recorded in previous studies [32,33].Experts were invited to interpret the other unconfirmed dark spots in Images 3-5 (see Figure 6c-e).The dark spots in Figure 6d were identified as low-wind-speed areas, whereas Figure 6c,e represent biogenic oil films.According to the above information, we manually interpreted and labeled different areas on the image as the ground truth.Table 2 shows the basic information of the images, with a central frequency of 5.405 × 10 9 Hz.Each data point was marked manually, and the targets were divided into five classes: clean sea, oil spill, look-alike, land, and ship areas.Images were divided into small areas of 320 × 320 pixels for training and testing.In order to increase the number of samples, data enhancement operations such as rotation and resampling were carried out in some areas.The ratio of training to test data was about 4:1.Tables 3 and 4 show the number of samples belonging to each class.In this paper, the mean intersection over union (MIoU) was taken as the main indicator to evaluate the classification effect.MIoU calculates the average value of the ratio of intersection and union of all classes, defined as where k is the number of classes, i represents the real value, j represents the predicted value, and Pij represents the prediction of i as j.In the semantic segmentation task, it is defined in the form of a confusion matrix as follows: where TP means true positive, FN means false negative, and FP means false negative.For each individual category, the intersection over union (IoU) is calculated as Other indices, such as the accuracy, precision, recall, and F1-score, are also used to indicate classification accuracy.They are defined as follows:

TP TN FN FP TP TN
2 Precision Recall F1_score Precision Recall

Polarization Decomposition and SLIC Superpixels
We applied polarization decomposition and SLIC superpixels to the SAR data according to the method introduced in Section 2. Two groups of polarimetric parameters were extracted: Group 1 (H/Alpha/λ1/λ2) parameters and Group 2 (Stokes) parameters.Table 5 calculates the average and variance of different polarized parameters and different areas, reflecting their differences in numerical characteristics.As an example, Figure 7 presents the ground truth and several samples of extraction results from the dataset.SLIC superpixels with different settings were calculated for all eight images, and then input into the convolutional neural network model together with the polarized parameters for training and testing.The MIoU results are listed in Table 6.The MIoU obtained using H/Alpha/λ1/λ2 parameters only was 70.6%, Stokes parameters achieved an MIoU of 72.0%.When the superpixel block contained 50 × 50, 75 × 75, 100 × 100, and 125 × 125 pixels, SLIC superpixels improved the classification accuracy.In the case of 150 × 150 and 200 × 200, the increase in MIoU was relatively low.Stokes parameters achieved the highest MIoU for the 125 × 125 setting (78.1%), whereas H/Alpha/λ1/λ2 parameters achieved an MIoU of 77.3% using the 100 × 100 superpixel setting.Table 7 lists the IoU of the test results in different categories with and without superpixels, when the superpixels were set as 125 × 125.The IoU values of clean sea and land areas were higher, while those of look-alike areas were relatively lower.The IoU of oil spill areas increased to 78.5% in both groups using superpixels.This indicates that look alike areas could not be distinguished precisely despite superpixel segmentation being applied.Figure 9 presents the test results of several samples.When only polarized parameters were used, there was extensive misclassification of look alike and oil spill areas.This phenomenon was improved after introducing superpixels.

Image Stretching Based on Superpixels
This section introduces the image stretching method based on SLIC superpixels for oil spill detection.Because the result of image stretching depends on the setting of the inhibition coefficient and expansion coefficient, the effects of these coefficients on oil spill detection was explored through experiments.For example, Figure 10

Oil Spill Detection as a Function of Inhibition Coefficient
We carried out further experiments while varying the inhibition coefficient.Firstly, we set the expansion coefficient to 1, and then we used inhibition coefficients of 1.05, 1.1, 1.15, and 1.2 to stretch the superpixel results.The image stretching results are presented in Figure 11.The results after image stretching and polarized parameters were used as input to the neural network model.Table 8 lists   The results in Table 8 prove that the inhibition coefficient could improve the MIoU.When the inhibition coefficient was set as 1.15 or 1.2, the accuracy of oil spill detection was much better than for other settings.When the inhibition coefficient was set as 1.2 and the Stokes parameter was applied, the MIoU was the highest at 80%.H/Alpha parameters achieved MIoU values of 78.8% and 78.7% with inhibition coefficients set as 1.15 and 1.2, respectively.Subsequent experiments were conducted with the inhibition coefficient set as 1.2.Table 9 presents the IoU value of each category when the inhibition coefficient was set as 1.2.Similar to when polarized parameters and superpixels were used, the accuracy was highest in the clean sea and land areas, slightly lower in the oil spill area, and lowest in the look alike and ship areas, as outlined in Table 9.Experiments evaluating the effect of different expansion coefficients were carried out.We fixed the inhibition coefficient as 1.2 on the basis of the results from the previous section, and we varied the expansion coefficients as 1.2, 1.4, 1.6, and 1.8. Figure 12 illustrates the image stretching results for different expansion coefficients.As the expansion coefficient was changed from 1.2 to 1.8, the characteristics of the oil spill area became more obvious.The stretching results and polarized parameters were combined as input into the neural network model for classification.The MIoU values with different expansion coefficients are listed in Table 10.The use of expansion coefficients greatly improved the accuracy of oil spill and lookalike areas, with the highest IOU of the oil spill area reaching 88.6% with Stokes parameters.When the expansion coefficient was set as 1.2, 1.4, and 1.6, the accuracy was improved sequentially, with the rate of improvement being greatest for look-alike areas.When the expansion coefficient was increased from 1.6 to 1.8, the MIoU value of each region decreased slightly, while the oil spill region still maintained a high accuracy, as outlined in Table 11.When the expansion coefficient was set too high, the pixel values of land, clean sea, look alike, and ship areas on the superpixel segmentation images were over-enlarged could not be easily distinguished.An accuracy of 0.933 was achieved using Stokes parameters and a 1.6 expansion coefficient.The precision, recall, and F1-score were also highest at 0.869, 0.926, and 0.897, respectively.When setting the expansion coefficient as 1.6, the accuracy was highest at 0.925 using H/Alpha parameters.Figure 13 presents the test results using polarized parameters and the image stretching method.Compared with only applying the inhibition coefficient, the classification accuracy of look alike and oil spill area was further improved when the expansion coefficient was applied.13 shows the IoU and MIoU results of oil spill detection with and without image stretching.On the basis of the previous experiment, the inhibition coefficient and expansion coefficient were set as 1.2 and 1.6.'-' indicates no samples of this category on the current image.
Table 13 shows that all the MIoU of the results of test images were higher with image stretching than without.The accuracy of each category was improved.The IoU of lookalike and oil spill areas improved by 6.6% and 4.9% on average, respectively.Oil spill achieved an MIoU of 81.6%, while clean sea and land areas could also be distinguished effectively.The average MIoU of the three images increased from 75.0% to 79.4%.Furthermore, with the exception of the land area, the mean square deviation of the IoU of different classes decreased, and the mean square deviation of MIoU values of different images also decreased from 3.49% to 2.72%.This proves that the image stretching method could significantly improve the accuracy of oil spill detection, as well as alleviate the inconsistency of accuracy caused by imaging differences.

Discussion
The experiments in this section focused on superpixel segmentation, the image stretching method, and polarized parameters.The number of superpixels determined the size of each superpixel block and affects the subsequent image stretching and classification results.The experimental results showed that, after 1:4 multi-look processing, the recognition accuracy was highest when each superpixel block contained 125 × 125 pixels, with the MIoU reaching 78.1% on the Stokes parameter.Image stretching involved setting the expansion and suppression coefficients and calculating the adaptive threshold to further process the segmentation results on the basis of superpixel segmentation.The purpose of this method was to divide the oil spill area and other areas in the image, so as to ensure the accurate classification of oil spill areas in the case of a large amount of data.When the expansion coefficient was set to 1.6 and the inhibition coefficient was set to 1.2, MIoU was increased to a maximum of 85.4%, achieving the highest accuracy.Subsequently, five images were used as the training set, and three images were used as the verification set.The experimental results showed that the average MIoU of different images increased from 75.0% to 79.4%, while the mean square deviation decreased from 3.49% to 2.72%, proving that this method could partially overcome the imaging differences and inconsistency of accuracy of different images while improving the accuracy.
The image stretching method used in this paper depends on the threshold calculation of Equation ( 16).Its purpose is to characterize the boundary between oil spill dark spots and other areas, as well as retain the characteristics of the oil spill area through the inhibition coefficient and expansion coefficient, so as to overcome the distribution differences between different SAR data.Therefore, oil spill areas on different SAR images were classified within the same pixel distribution, while other similar regions were classified in another pixel distribution.In fact, it was found that the threshold is not an exact numerical value but usually an interval, and the threshold calculated by the formula fell into this interval.Taking Image 2 as example, Figure 14a,b show the processing results when the threshold was manually set to 0.67 and 0.73, respectively, while Figure 14c shows the processing results when the threshold was 0.69 according to the formula, and these three results were incredibly similar.

Conclusions
Polarized SAR data are widely used in oil spill detection.With the development of deep learning, more and more neural network models have been applied to oil spill detection and achieved good results.In this paper, we proposed an adaptive threshold calculation and image stretching method based on SLIC superpixels to overcome the influence of imaging differences between different SAR data, and we designed a semantic segmentation model to identify oil spill areas.Experiments were carried out on eight Sentinel-1 dual-pol SAR images.The experimental results showed that polarized parameters combined with SLIC superpixel segmentation could not effectively distinguish oil spill areas and similar areas, with the highest MIoU value only reaching 78.1%.After image stretching using the inhibition coefficient and expansion coefficient, the recognition accuracy of the two methods was greatly improved.The highest MIoU of the test set improved to 85.4%, while the accuracy and other indicators were also higher than those in other groups.Experiments on single images showed that this method partly overcame the imaging difference and inconsistent accuracy of Sentinel-1 dual-pol SAR data.
There were still some limitations of our current research.The optimal superpixel number used in the SLIC superpixel model depends on actual situations of the study site.This paper only discussed the impact of superpixel number on the current dataset, which may need to be adjusted according to the specific situation.Thus, future research can focus on adaptive optimization of the superpixel number.The current threshold of image stretching is set by experiments, and more data experiments and theoretical analyses are required.Future research can focus on expanding the data experiments and the adaptive threshold calculation, including numerical fitting and machine learning.From the perspective of the neural network, further research can be conducted on the improvement of the model structure explored in this paper.Introducing an appropriate attention mechanism is one option for enhancing the extraction and utilization of texture information in oil spill areas.
Author Contributions: Funding acquisition, Q.L. and Y.L.; methodology, J.Z. and Q.L.; supervision, Q.L., H.F., Y.Z. and J.L.; resources, J.L. and Z.Z.; writing-original draft, J.Z.; writing-review and editing, Q.L. and Y.L.All authors have read and agreed to the published version of the manuscript.

Figure 1 .
Figure 1.Overall flowchart of the proposed oil spill detection algorithm.

Figure 9 .
Figure 9. Classification results of some typical samples using polarized parameters and SLIC superpixels.Red represents oil spill, blue represents look alike, while green represents land areas.(a) oil spill, (b) look alike, (c) oil spill, (d) look alike, (e) oil spill, (f) look alike, (g) land.

Figure 13 .
Figure 13.Classification results of several samples using polarized parameters and image stretching.Red represents oil spill, blue represents look alike, while green represents land areas.(a) oil spill, (b) look alike, (c) oil spill, (d) look alike, (e) oil spill, (f) look alike, (g) land.

Funding:
This research was partially funded by the Key Project of Tianjin Natural Science Foundation (21JCZDJC00670), by the National Engineering Laboratory for Digital Construction and Evaluation Technology of Urban Rail Transit (grant No. 2021ZH04), by the Tianjin Transportation Science and Technology Development Project (grant No. 2022-40, 2020-02), and by the National Natural Science Foundation of China Grant (grant No. 41601446).

Table 1 .
Detailed parameters of neural network model.
* Numbers in parentheses represent dilation rates.

Table 3 .
Number of pixel samples of different classes.

Table 4 .
Number of samples of different classes.

Table 5 .
Average and variance of different polarized parameters.

Table 6 .
MIoU of classification results based on polarimetric parameters with different SLIC superpixel numbers.

Table 7 .
IoU values of oil spill detection with and without superpixels.

Table 8 .
MIoU value of different inhibition coefficients.

Table 9 .
Characteristics of polarized parameters in experiments.

Table 10 .
MIoU value of different inhibition coefficients.When setting the expansion coefficient to 1.2, 1.4, 1.6, and 1.8, all MIoU values of the results were improved compared with using the inhibition coefficient only.The highest MIoU was achieved with the expansion coefficient set as 1.6, followed by 1.4, 1.8, and 1.2.The results for an expansion coefficient of 1.4 were very close to those for 1.8.The MIoU based on Stokes parameters reached 85.4%, while H/Alpha parameters achieved an MIoU of 84.5% with the expansion coefficient.Tables11 and 12list the IoU values of different classes and multiple indicators, which were used to evaluate the model with different expansion coefficients.

Table 11 .
IoU Values of polarized parameters with 125 × 125 superpixels and different expansion coefficients.

Table 12 .
Indicators of Polarized parameters with 125 × 125 superpixels and different expansion coefficients.

Table 13 .
IoU and MIoU values of validation experiments with and without image stretching.