Next Article in Journal
Normalized Radar Scattering Section Simulation and Numerical Calculation of Freak Wave
Previous Article in Journal
Quantifying Efforts to Mitigate Interference between an Unmanned Surface Vessel and Starfish 990F for the Identification of Underwater Features in the Littoral Zone
Previous Article in Special Issue
Experimental Characterisation and Field Experience of a Reusable, Modified Polyurethane Foam for the Mechanical Clean-Up of Oil Spills on the Sea Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection of Massive Oil Spills in Sun Glint Optical Imagery through Super-Pixel Segmentation

1
Institute of Estuarine and Coastal Research, School of Ocean Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, China
2
School of Marine Sciences, Sun Yat-sen University, Zhuhai 519082, China
3
Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519000, China
4
Guangdong Provincial Key Laboratory of Marine Resources and Coastal Engineering, Guangzhou 510275, China
5
Pearl River Estuary Marine Ecosystem Research Station, Ministry of Education, Zhuhai 519000, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(11), 1630; https://doi.org/10.3390/jmse10111630
Submission received: 27 September 2022 / Revised: 19 October 2022 / Accepted: 23 October 2022 / Published: 2 November 2022

Abstract

:
Large volumes of crude oil accidentally released into the sea may cause irreversible adverse impacts on marine and coastal environments. Large swath optical imagery, acquired using platforms such as the moderate-resolution imaging spectroradiometer (MODIS), is frequently used for massive oil spill detection, attributing to its large coverage and short global revisit, providing rich data for oil spill monitoring. The aim of this study was to develop a suitable approach for massive oil spill detection in sun glint optical imagery. Specifically, preprocessing procedures were conducted to mitigate the inhomogeneous light field over the spilled area caused by sun glint, enhance the target boundary contrast, and maintain the internal homogeneity within the target. The image was then segmented into super-pixels based on a simple linear clustering method with similar characteristics of color, brightness, and texture. The neighborhood super-pixels were merged into target objects through the region adjacency graph method based on the Euclidean distance of their colors with an adaptive termination threshold. Oil slicks from the generated bright/dark objects were discriminated through a decision tree with parameters based on spectral and spatial characteristics. The proposed approach was applied to oil spill detection in MODIS images acquired during the Montara oil spill in 2009, with an overall extraction precision of 0.8, recall of 0.838, and F1-score of 0.818. Such an approach is expected to provide timely and accurate oil spill detection for disaster emergency response and ecological impact assessment.

1. Introduction

Oil spills from oil production (e.g., from drilling rigs) or during oil transportation (e.g., from oil tankers) represent a large portion of large oil spill accidents. Such large volumes of oil discharge may cause irreversible adverse impacts on marine and coastal environments. For example, research related to the 2010 Deepwater Horizon oil spill demonstrated that the oil spill caused significant declines in the abundance and diversity of coastal wetlands, deep-sea coral reefs, fish species, sea turtles, and cetaceans [1,2,3,4,5].
Recent advancements in technology and sensor availability have enabled remote sensing to play a vital role in oil spill detection and monitoring, taking advantage of its synoptic observation over a wide area and frequent revisit capabilities compared to traditional observation methods [6,7,8,9,10]. Of the various remote sensing techniques, synthetic-aperture radar (SAR) and optical remote sensing are most frequently used in oil spill detection on the ocean surface. Oil films on the ocean surface would dampen surface capillary and short gravity waves, creating a “smoother” surface compared to the uncontaminated nearby water, reducing the Bragg scattering signal under optimal wind conditions, and appearing darker than the surrounding water in SAR imagery [11]. SAR has day and night capabilities and can penetrate through clouds, enabling highly effective detection in oil spill responses. While SAR is primarily used to detect the presence of oil, recent development has demonstrated the capabilities of discriminating thick oil emulsions using quad-polarized SAR images [11,12,13].
Optical remote sensing is also frequently used in oil spill detection, especially for large oil spill accidents. Wide-swath sensors such as the moderate-resolution imaging spectroradiometer (MODIS) and visible infrared imaging radiometer (VIIRS) provide frequent revisit globally (1–2 days in equatorial orbit), which would largely compensate for the weaknesses of lack of coverage under clouds [14,15]. There are two main principles behind the optical detection of oil spills. First, the oil-dampening effect can modulate reflected sunlight, which enhances the oil–water spatial contrast under sun glint. Distinct from the dark spots observed in SAR imagery, oil-covered regions in optical sun glint imagery display either positive or negative contrast with nearby water depending upon solar/viewing geometry and surface roughness [15,16,17,18]. However, it should be noted that when sun glint is negligible or absent (e.g., LGN < 10−6 sr−1 for MODIS), oil films may be undetectable [18]. Second, oil has unique optical properties, including high absorption in the blue spectral band, and exponential decay toward longer wavelengths in crude oil. Mixtures of oil–water emulsions enable high scattering, creating higher reflectance in the red, near-infrared (NIR), and shortwave infrared bands. As the multi-spectral information in optical remote sensing can be used to rule out oil look-alikes [19,20], classify oil types [20,21,22], and quantify oil concentration or volume [23,24,25,26], it shows high promise for future research on oil detection through optical remote sensing.
While many semi-/automatic oil spill extraction algorithms have been developed based on SAR [27,28,29,30,31,32], there are few semi-/automatic oil spill detection algorithms working for optical imagery. Currently, oil spill extent detection in optical remote sensing images largely relies on visual interpretation and manual delineation [33,34,35], which require extensive experience and expertise; and the derived oil spill map is also subject to interpretation. Other than visual interpretation, methods including pixel-based indices and deep learning have been developed to identify and classify oil slicks based on optical imagery [36,37,38,39,40,41,42,43,44]. Such methods have achieved good results in small- to mid-sized oil spills (~5–20 km in size). While sun glint reflectance varies depending on the surface roughness and solar/viewing geometry, the above methods may not apply to massive/large oil spills (hundreds of km, or covering a large part of the image swath), because of differences in the ambient light field caused by sun glint over the covered area. A general robust satellite technique based on the historical long-term time series of satellite data and local anomaly detection in the NIR band has been applied to massive oil spill detection in sun glint MODIS imagery during the Deepwater Horizon oil spill [43,44]. However, such an approach has difficulties in discriminating the negative-contrast oil slicks because of the low range of variability in reflectance.
Therefore, the objective of this study is to develop an oil spill detection approach for large oil spill accidents to overcome the large reflectance variability over the spilled area and to accomplish the identification of both the negative and positive contrast oil slicks in sun glint imagery. In this paper, an object-based method consisting of super-pixel segmentation and neighborhood merging is proposed to identify oil slicks through object feature analysis. The approach was tested and applied to the extraction of oil slicks with different spatial sizes and under different sun glint intensities during the Montara oil spill, and the extraction performance was evaluated objectively by comparing it with manually delineated oil slicks.

2. Data

2.1. Study Area

The Montara oil field is located in the Australian sector of the Timor Sea, ~260 km off the northern shoreline of Western Australia. A blowout of oil and gas happened on 21 August 2009, at a water depth of 77 m (Figure 1), due to accidental loss of control of the H1 well of the Montara wellhead. Subsequent oil slicks were found in the Timor Sea. The uncontrolled well continued to discharge oil into the Timor Sea until it was finally capped on 3 November 2009. An estimated 29,600 barrels of crude oil were discharged into the Timor Sea, making the Montara oil spill one of the worst oil disasters in Australian history [45,46]. The study area of this paper includes the main polluted area of the oil spill (14–11° S, 123.5–126.5° E, Figure 1).

2.2. Data and Preprocessing

MODIS L1A data were acquired from the NASA ocean color data archive. The data were processed to Rayleigh-corrected reflectance (Rrc, dimensionless) through the SeaWiFS data analysis system (SeaDAS, version 7.5.3, NASA Goddard Space Flight Center, Greenbelt, MD, USA), with bands centered at 469, 555, 645, 859, 1240, and 1640 nm. These bands were used for spectral analysis, where bands centered at 859, 645, and 555 nm were used for segmentation, and the 859-nm band was used for additional object analysis in this study. Here, the land bands rather than the ocean bands were selected because the ocean bands tend to saturate under strong sun glint. Sea surface temperature (SST, in units of °C) data were also acquired through the SeaDAS processing. In this study, sun glint strength was defined as the normalized sun glint radiance (LGN, in units of sr−1). LGN was also acquired through the SeaDAS processing, and calculated via the Cox and Munk model based on the input from the wind speed (to calculate surface roughness) and the solar/viewing geometry [47,48,49]. The datasets were then mapped to a rectangular projection and all the bands were resampled to 250-m resolution. Appropriate image preprocessing procedures were used to mitigate the inhomogeneous light field over the image, enhance the contrast at the target boundary, and maintain homogeneity inside the target. The series of procedures conducted for preprocessing after the generation of Rrc data are summarized in Figure 2. Firstly, a selected cloud-free image was used to mask the land pixels, which was achieved with the normalized difference water index (NDWI) calculated using the MODIS green and NIR bands centered at 555 and 859 nm, respectively, with a threshold of <−0.05 for land and >0.5 for submerged islands [50]. Land mask, SST, and LGN datasets were then used as auxiliary data to discriminate oil slicks.
Secondly, the false color composite (FRGB, centered at 859, 645, and 555 nm) was generated after Gaussian stretching [51,52]. Here, instead of Rrc data, the FRGB was used in order to increase the computation efficiency. As can be seen from the transect in Figure 3a, the dark oil area generally shows lower brightness in the NIR band than in the surrounding water. However, the ambient light field is not homogeneous across the image: the NIR band brightness of the water decreases gradually from the left to the right side of the image because of the decreased sun glint intensities (LGN in Figure 3d). Moreover, such a gradual change in magnitude due to sun glint overwhelms the local contrasts between the oil and water, making it difficult to extract oil slicks based on the fixed threshold method.
Thirdly, the gradual color intensity variation across the image caused by sun glint needed to be minimized. This was achieved by subtracting a gradient background, which was modeled using the SFIT function in Interactive Data Language (IDL, version 8.5, Exelis Visual Information Solutions, Inc., a subsidiary of Harris Corporation, Boulder, CO, USA) software to determine the polynomial fit to the three channels of FRGB [53]. Taking the first channel (converted from Rrc_859) as an example, the brightness of the water became more homogeneous after subtraction of the gradient background caused by sun glint (Figure 3b,e vs. Figure 3a,d). Finally, the target boundary contrast was enhanced, and the target internal homogeneity was achieved through additional subtraction of median filtering with the kernel size determined through trial and error [54,55]. While the water brightness over the transect became more homogeneous, the brightness contrast between the oil and water was enhanced, as displayed in Figure 3c,f. The parameter settings in the preprocessing procedures are listed in Table 1.

3. Methods

3.1. Work Flow

The approach proposed in this study employs the commonly used chain structure flow in oil spill detection (Figure 4), which generally has three steps [28].
The first step, which consisted of two stages, was to generate the various objects: first, the image was segmented into super-pixels with similar color, brightness, and texture through the simple linear clustering (SLIC) method [56]; then, the super-pixels were merged into objects with different scales based on the region adjacency graph (RAG) method [57,58]. This is further explained in Section 3.2.
The next step was to separate the color anomaly objects from the background by analyzing the relationship between the sun glint reflectance and the lightness (in the CIELAB color space) of the background. Then, it needed to be decided whether the object belonged to the bright, or dark, object category. This is further explained in Section 3.3.
Finally, oil look-alikes needed to be excluded from the derived objects. These false positives included cloud shadow and sea surface roughness anomalies in dark objects, and seawater under sun glint and clouds in bright objects. These oil-look-alike objects were excluded through a decision tree whose nodes use various parameters including object-based spectral and spatial characteristics. This is further explained in Section 3.4.

3.2. Object Generation Based on Segmentation

SLIC segmentation was conducted after the image preprocessing to generate super-pixels from the FRGB image. A super-pixel is a small region consisting of multiple pixels that are located next to each other and have similar characteristics including color, brightness, and texture. These small regions (super-pixels) retain effective information and generally do not destroy the boundary information of objects in the image [59,60]. The SLIC method consists of three parts: color space transformation, distance metric construction, and iterative clustering. Firstly, the pixels of the FRGB image are transformed into five-dimensional feature vectors in the CIELAB color space and XY coordinates. The CIELAB color model can be considered a re-projection of the RGB color model, whose colors are composed of lightness (L), the green-to-red component (A), and the yellow-to-blue component (B). CIELAB is designed to approximate human vision, with the L component closely matching the human perception of lightness, which defines black at 0 and white at 100. Secondly, the similarity of two pixels can be measured by the distance between their corresponding five-dimensional vectors. The larger the distance, the smaller their similarity. The algorithm first generates K seed points, then searches the space surrounding each point, and finally classifies similar pixels into the same class with the seed until all pixels are classified. Thirdly, the average vector value in these K super-pixels is calculated to derive a recalculated center of the K clusters. The new K centers are then used to search for the most similar pixels around them. The above procedures are conducted iteratively until convergence, where the K centers do not shift [56,61,62]. Figure 5b shows the segmentation results of the SLIC method. Here, the color of the super-pixel represents the average brightness of the corresponding pixels in the RGB color space. The water areas were divided into various rectangular super-pixels with similar size and color. However, water super-pixels adjacent to other targets (such as oil slicks or clouds) were irregularly shaped, with their boundaries fitting the shape of the other targets. The parameter settings of the SLIC segmentation are listed in Table 1.
RAG is a typical topological data structure to characterize the adjacency of super-pixels [57,58]. Based on graph theory, RAG mergence was then conducted to merge super-pixels of a similar color into one object based on the Euclidean distance of three-dimensional feature vectors in the RGB color space. The over-segmented super-pixels with a similar color were gradually merged until no similar pairs of regions were found, which was operated by judging if the color difference was larger than a termination threshold. In this study, an adaptive termination threshold of neighborhood mergence was used, which was determined as follows:
T N M = 0.35 × N m 2 + R m 2 + G m 2
where TNM represents the adaptive termination threshold for neighborhood merging, and Nm, Rm, and Gm represent the median values of the three channels of the FRGB image.

3.3. Bright/Dark Object Recognition

As mentioned, oil slicks show either a positive contrast or a negative contrast with the water under the influence of sun glint, depending on the surface roughness and solar/viewing geometries. Therefore, the derived objects were first separated from the background into bright objects and dark objects. As most pixels within the study area were seawater, the median value of L in the CIELAB of the image was used to define the category of background seawater. When LGN is weak, the contrast between the dark object lightness (Lo) and the seawater background lightness (Lbg) is small. The contrast becomes larger with increased LGN values. The same applies to the contrast of bright objects with background seawater. The process of bright/dark object recognition can be divided into two steps. Firstly, objects that meet the conditions set by the different LGN are classified as bright targets if their lightness is greater than Lbg. Secondly, objects that satisfy the conditions set by the different LGN are classified as dark targets if their lightness value is smaller than Lbg. The rules for separating bright and dark objects with different sun glint conditions are shown in Figure 6.
As shown in Figure 5d, most of the oil slick objects were successfully segmented into objects other than the background. However, oil slicks were divided into several different objects due to the non-homogeneity within the slick object. After the bright/dark object recognition, these connected parts of the oil slick were merged into one through an eight-neighborhood connectivity analysis.

3.4. Oil Spill Identification

The characteristics of the derived objects were used to discriminate oil slicks from the derived dark and bright objects, which included four parameters and one decision condition: (1) the mean value of auxiliary data (LGN, SST) of the object, (2) the mean difference between the outer boundary pixels and the inner boundary pixels of the objects in the NIR reflectance (Rrc at 859 nm), (3) the mean reflectance difference between the object and a calculated shift object in the NIR band, (4) the size of the object, and (5) whether or not the object was located near land.
The discrimination of oil slicks was separated for bright and dark objects. For the discrimination of dark targets, false positive targets included surface roughness anomalies and cloud shadows. Sea surface roughness anomalies were ruled out through the object contrast with the background water. While sea surface roughness anomaly regions in low-wind conditions usually display an extensive dark area with scattered dark spots in the boundary area in SAR imagery, oil-slick-covered areas tend to show a sharp contrast at the slick boundary [63]. A similar phenomenon can be observed in optical imagery under sun glint. Therefore, the object boundary tends to display a larger contrast at the oil-slick-covered area than at the surface roughness anomaly area. In this study, the NIR mean reflectance difference between the inner boundary pixels (brown pixel with a red frame in Figure 7b) and the outer boundary pixels (blue pixel with a black frame in Figure 7b) of the bright/dark objects, i.e., MDN (mean difference to neighborhood pixels), was used to discriminate oil slicks from the surface roughness anomaly area. Under the same sun glint condition, the contrast between cloud and cloud shadow was larger than that between oil slicks and the background water. Due to the different local overpass times of MODIS Terra and Aqua, cloud shadow in the study area is usually observed on the west side of the cloud with Terra and the east side of the cloud with Aqua. Such cloud shadow was recognized through the reflectance difference in the NIR band between the bright and dark objects. A potential cloud region was searched by shifting the dark objects. Taking the Aqua image as an example, a potential cloud area (purple pixels with a black frame in Figure 7a) was searched by shifting N pixels to the west of the dark object (gray pixels in Figure 7a). N is the average value of the width of the dark object. The maximum value of N was set to 5, as a statistical analysis showed that the mean width of cloud shadows in the study area rarely exceeded 5 pixels. The mean difference in the NIR reflectance between dark objects and the potential cloud area, i.e., MDC (mean difference to cloud object), was calculated. Through a statistical analysis, we found that the larger the LGN, the larger the threshold values of MDN and MDC. The MDN and MDC parameter threshold settings for dark oil detection under different LGN values are shown in Figure 8. The dark objects from shear zones and the leeward side of the island were discriminated by checking whether the dark object was the neighborhood of the land mask. If any pixel within the dark object was less than a two-pixel distance from the land mask, the dark object was identified as a non-oil object.
False-positive targets in bright objects include clouds and water under strong sun glint. Large-size clouds were identified through the auxiliary SST data, as the difference in the SST between cloud and seawater was significantly larger than that between oil slicks and seawater [64]. If the mean SST of the bright object (SSTo) differs from the median value of the auxiliary SST of the image by a threshold SSTd, the object would be identified as a non-oil object. However, such cloud discrimination cannot be applied to clouds of small size because SST data with a resolution of 1 km are not sufficient to detect small-sized clouds. Instead, small-sized clouds were discriminated using procedures similar to the cloud shadow discrimination. Taking the Aqua image as an example, a potential cloud shadow area (gray pixels with a black frame in Figure 7c) was searched by shifting the bright object (purple pixels in Figure 7c) N pixels eastward without intersecting with the bright object. The mean difference between the bright object and the potential cloud shadow in the NIR band reflectance, i.e., MDCS (mean difference to cloud shadow object), was calculated. The MDCS in bright objects was lower for the oil slick area than for the cloud area. False positives from water under high sun glint were discriminated from the bright oil through the MDN, with the same principle as that of the discrimination of the roughness anomaly seawater. Through a statistical analysis, we found that the larger the LGN, the larger the threshold values for the MDN and MDCS parameters. The positive-contrast oil slick identification rules and the threshold settings of the SSTd, MDN, and MDCS parameters under different LGN values are shown in Figure 9.
Finally, small-sized clouds or cloud shadow objects that were not identified through the above processes were excluded through a sliding window method. The sliding window length/width and step were set to 1/6 of the image length/width, and the threshold was set to 30 pixels. The average size of the objects in the sliding window was calculated. If the average size was less than the threshold, the objects in the window were discriminated as non-oil objects. By this means, the dispersed cloud and associated cloud shadows were removed.

3.5. Performance Evaluation

When evaluating the performance of the approach, both the positive contrast and negative contrast oil objects were regarded as oil targets, while other targets were regarded as background targets. To objectively evaluate the performance of the approach, the oil slicks in the images were visually interpreted and manually delineated in ArcMap software (version 10.8, Environment System Research Institute, Redlands, CA, USA) [33]. These interpreted oil slicks were used to assess the performance using the approach proposed in this study. Three evaluation metrics widely used for target recognition were chosen in this study for evaluation [65,66,67]: precision, recall, and F1-score, which are defined as
P r e c i s i o n = T P T P + F P
R e c a l l = T P T P + F N
F 1 = 2 × T P 2 × T P + F N + F P
where TP (true positives) denotes the number of target objects (pixels) classified into the correct category, TN (true negatives) is the number of non-target objects classified into the correct category, FN (false negatives) is the number of target objects classified into the wrong category, and FP (false positives or false alarms) is the number of objects from a non-target category incorrectly classified as a target. High precision indicates high detection accuracy, and high recall suggests a high number of positive objects detected. F1 is an indicator used in statistics to measure the performance of a binary model, which considers both the precision and recall of the classification results, and it can be regarded as a weighted average of model precision and recall with a maximum value of 1 and a minimum of 0.

4. Results

4.1. Optical Characteristics of Oil Spill

With the same glint strength at the oil–water boundary, the oil and water contrast can be clearly observed due to the oil’s modulation of roughness (Figure 10a–c). The images in Figure 10a,b were taken at different times on the same day. The sun glint strength decreases from left to right in Figure 10a (0.052~0.002 sr−1), while the opposite is found in Figure 10b (0.0002~0.025 sr−1). Different oil–water contrasts can be found at the same oil slick, such as dark oil (blue circle point in Figure 10a) and bright oil (green circle point in Figure 10b). Here, the oil slick or area of the slick in the negative contrast with the water is referred to as dark oil, while the positive contrast slick or area of the slick is referred to as bright oil in this paper.
Oil look-alikes can also display similar positive/negative contrast characteristics, such as clouds, water under high sun glint, cloud shadows, and roughness anomaly seawater. Oil slicks and oil look-alikes under different sun glint conditions are presented in Figure 10a–c, and the corresponding reflectance spectra of the selected targets are displayed in Figure 10d. When sun glint strength differs, reflectance spectra of different targets may resemble each other (Figure 10d). For example, the sea surface roughness anomaly (Figure 10c, pink square point) and cloud shadow (Figure 10a, pink triangle point) resembled oil (Figure 10a, blue circle point) spectra in shape, and their reflectance magnitude was also close to that of oil slicks. This dark roughness anomaly region was caused by the local low wind (2.2 m/s in the pink square point vs. 4.5 m/s in the green square point; reflectance spectra displayed in Figure 10d). Seawater under high sun glint (Figure 10b, red square point) and clouds (Figure 10c, red triangle point) resembled oil (Figure 10c, red circle point) spectra in shape, with the reflectance magnitude also being close to that of oil slicks (Figure 10d). The target reflectance magnitude and shape also changed with the sun glint strength. As can be seen from the selected water locations in Figure 10, the water reflectance spectra decrease in all spectral bands with the decreased sun glint strength (LGN = 0.04 sr−1 at the red square point vs. 0.03 sr−1 at the green square point vs. 0.02 sr−1 at the blue square point).

4.2. Extracted Oil Slicks

Figure 11 shows the results of the extracted oil slicks using the approach proposed in this study. Most of the detected oil slicks were within 100 km of the platform location for the duration of the spill, with almost no slicks observed 210 km away from the platform. The detected area of oil slicks ranged from 1131 to 5310 km2 in the MODIS images, with an average area of 3314 km2. It can be seen that most of the oil slicks were sufficiently extracted, with false positive targets in the figure (clouds, water under strong sun glint, cloud shadows, and sea surface roughness anomaly area) being sufficiently removed. The oil slicks in Figure 11a,b,e,g are characterized by many dispersed slick patches in the image. They are composed of either one major slick with small scattered oil patches (Figure 11e,g) or several large slicks with small scattered patches that are distributed some distance from each other (Figure 11a,b,e). This approach achieved a good performance in extracting these types of oil slicks, with some minor omissions, demonstrating the approach’s capabilities in dealing with oil slicks of multiple scales. Figure 11c,d,f,i present the negative contrast of oil slicks with seawater. Oil slicks in these images are composed of one or several large patches close to each other. As can be seen, the extraction results show a clean slick delineation and are very close to the manual interpretation results. An oil slick with positive contrasting features can be observed in Figure 11h, and both positive and negative contrasting features can be observed in Figure 11e,g. The overall extraction results are strong, with all major slicks being identified, but with some false positive detections from clouds and cloud shadows (Figure 11e,g,h).

4.3. Performance Evaluation

Table 2 displays the quantitative assessment of the extraction results against the manually delineated slicks through visual interpretation. The statistical results demonstrate an overall precision value of 0.8, which means that 80% of the oil pixels were detected correctly. The overall recall value of 0.838 demonstrates that 83.8% of the oil spill pixels in the image were successfully detected.
For individual images, the recall of the approach ranged from 0.780 to 0.879. The main factors affecting the recall included the existence of clouds in the image and inhomogeneity within the oil slick. Cloud cover on top of an oil slick causes it to be divided into several parts. However, the dark oil objects would be mistakenly identified as cloud shadow if the relationship between the dark oil and cloud objects meets the third rule in Section 3.4, as shown by the omitted oil slicks in the regions marked by red ellipses in Figure 11a,b,i,j. The coexistence of bright and dark oil objects would also lead to similar effects described above (recognized as clouds and cloud shadow), leading to the omission of oil slicks in the extraction results (Figure 11e,f,g).
The precision of the extracted results ranged from 0.423 to 0.969 for individual images. Among them, three of the images had a precision of <0.8. The main factor affecting the precision is the sea surface roughness anomaly adjacent to the clouds, such as the regions marked by a red rectangle in Figure 11. Although there were only a few misidentified objects, the impact on the precision of the approach was great because of the large area of misidentified objects.

5. Discussion and Conclusions

The detection performed excellently for situations where the oil slick distribution in the image is concentrated, and where the oil spill area is not covered by clouds, as shown in Figure 11a–d,f,i. The extraction results of these situations display a good performance, with F1-scores of over 0.85. In circumstances where there is a multiple-target distribution in the image, e.g., the coexistence of surface roughness anomalies and clouds, the extraction result shows a low F1-score due to a large number of detected false positives and the small number of undetected false negatives. As shown in Figure 11h,j, the extraction results of these two situations are the worst among the 10 images, with F1-scores of 0.576 and 0.571.
As presented in this study, precise cloud/shadow detection and masking before the segmentation can significantly improve the performance of the approach proposed in this study. However, the MODIS cloud mask products would mistakenly recognize the bright oil slicks as clouds, prohibiting the application of such a cloud mask product to the approach, which calls for more precise cloud mask algorithms under strong sun glint conditions.
In conclusion, an object-based approach was proposed in this study to extract massive oil slicks in sun glint MODIS imagery of the Montara oil spill that occurred in the Timor Sea in 2009. After appropriate imaging preprocessing procedures including Gaussian stretching, median filtering, and sun glint gradient modeling, the SLIC and RAG algorithms were used to segment the image into super-pixels and to merge the neighborhood super-pixels into target objects with varying sizes. False positive detections from the objects were ruled out through a decision tree using object spectral and spatial characteristics. The approach achieved a good performance in extracting the oil slicks from the Montara oil spill, with a slick extraction precision of 0.8, recall of 0.838, and F1-score of 0.818. The images used in this study cover a broad range of situations including different sun glint strengths and oil slicks of multiple sizes with various distribution patterns, accompanied by various oil slick look-alikes. The performance assessment proved the robustness of the approach, which represents a useful approach for oil slick extraction using sun glint optical images during massive oil spills. Moreover, as the approach proposed in this study is based mainly on reflectance data in the 859, 645, and 555 nm bands, which can be commonly found in other optical sensors, this approach may apply to other optical images as well.

Author Contributions

Conceptualization, Z.S. and S.S.; Funding acquisition, S.S.; Methodology, Z.S.; Writing—original draft, Z.S. and S.S.; Writing—review & editing, Z.S., S.S., J.Z., B.A. and Q.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 42106173, and the Guangdong Basic and Applied Basic Research Foundation, grant number 2020A1515110957.

Data Availability Statement

MODIS data are openly accessible from https://oceancolor.gsfc.nasa.gov/ (accessed on 25 October 2022). The result dataset will be made available on request.

Acknowledgments

The authors would like to thank NASA for providing MODIS data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beyer, J.; Trannum, H.C.; Bakke, T.; Hodson, P.V.; Collier, T.K. Environmental effects of the Deepwater Horizon oil spill: A review. Mar. Pollut. Bull. 2016, 110, 28–51. [Google Scholar] [CrossRef] [Green Version]
  2. Shapiro, K.; Khanna, S.; Ustin, S.L. Vegetation Impact and Recovery from Oil-Induced Stress on Three Ecologically Distinct Wetland Sites in the Gulf of Mexico. J. Mar. Sci. Eng. 2016, 4, 33. [Google Scholar] [CrossRef] [Green Version]
  3. Girard, F.; Fisher, C.R. Long-term impact of the Deepwater Horizon oil spill on deep-sea corals detected after seven years of monitoring. Biol. Conserv. 2018, 225, 117–127. [Google Scholar] [CrossRef]
  4. McClain, C.R.; Nunnally, C.; Benfield, M.C. Persistent and substantial impacts of the Deepwater Horizon oil spill on deep-sea megafauna. R. Soc. Open Sci. 2019, 6, 191164. [Google Scholar] [CrossRef] [Green Version]
  5. Asif, Z.; Chen, Z.; An, C.J.; Dong, J.X. Environmental Impacts and Challenges Associated with Oil Spills on Shorelines. J. Mar. Sci. Eng. 2022, 10, 762. [Google Scholar] [CrossRef]
  6. Hu, C.M.; Feng, L.; Holmes, J.; Swayze, G.A.; Leifer, I.; Melton, C.; Garcia, O.; MacDonald, I.; Hess, M.; Muller-Karger, F.; et al. Remote sensing estimation of surface oil volume during the 2010 Deepwater Horizon oil blowout in the Gulf of Mexico: Scaling up AVIRIS observations with MODIS measurements. J. Appl. Remote Sens. 2018, 12, 026008. [Google Scholar] [CrossRef] [Green Version]
  7. Sun, S.J.; Lu, Y.C.; Liu, Y.X.; Wang, M.Q.; Hu, C.M. Tracking an Oil Tanker Collision and Spilled Oils in the East China Sea Using Multisensor Day and Night Satellite Imagery. Geophys. Res. Lett. 2018, 45, 3212–3220. [Google Scholar] [CrossRef]
  8. Hebbar, A.A.; Dharmasiri, I.G. Management of marine oil spills: A case study of the Wakashio oil spill in Mauritius using a lens-actor-focus conceptual framework. Ocean Coast. Manag. 2022, 221, 106103. [Google Scholar] [CrossRef]
  9. Rajendran, S.; Vethamony, P.; Sadooni, F.N.; Al-Kuwari, H.A.; Al-Khayat, J.A.; Seegobin, V.O.; Govil, H.; Nasir, S. Detection of Wakashio oil spill off Mauritius using Sentinel-1 and 2 data: Capability of sensors, image transformation methods and mapping. Environ. Pollut. 2021, 274, 116618. [Google Scholar] [CrossRef]
  10. Guo, G.; Liu, B.X.; Liu, C.Y. Thermal Infrared Spectral Characteristics of Bunker Fuel Oil to Determine Oil-Film Thickness and API. J. Mar. Sci. Eng. 2020, 8, 135. [Google Scholar] [CrossRef]
  11. Alpers, W.; Holt, B.; Zeng, K. Oil spill detection by imaging radars: Challenges and pitfalls. Remote Sens. Environ. 2017, 201, 133–147. [Google Scholar] [CrossRef]
  12. Garcia-Pineda, O.; Staples, G.; Jones, C.E.; Hu, C.M.; Holt, B.; Kourafalou, V.; Graettinger, G.; DiPinto, L.; Ramirez, E.; Streett, D.; et al. Classification of oil spill by thicknesses using multiple remote sensors. Remote Sens. Environ. 2020, 236, 111421. [Google Scholar] [CrossRef]
  13. Espeseth, M.M.; Jones, C.E.; Holt, B.; Brekke, C.; Skrunes, S. Oil-Spill-Response-Oriented Information Products Derived From a Rapid-Repeat Time Series of SAR Images. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 3448–3461. [Google Scholar] [CrossRef]
  14. Leifer, I.; Murray, J.; Streett, D.; Stough, T.; Ramirez, E.; Gallegos, S. The Federal Oil Spill Team for Emergency Response Remote Sensing, FOSTERRS: Enabling Remote Sensing Technology for Marine Disaster Response. In Time-Sensitive Remote Sensing; Lippitt, C.D., Stow, D.A., Coulter, L.L., Eds.; Springer: New York, NY, USA, 2015; pp. 91–111. [Google Scholar]
  15. Sun, S.J.; Hu, C.M. The Challenges of Interpreting Oil–Water Spatial and Spectral Contrasts for the Estimation of Oil Thickness: Examples From Satellite and Airborne Measurements of the Deepwater Horizon Oil Spill. IEEE Trans. Geosci. Remote Sens. 2019, 57, 2643–2658. [Google Scholar] [CrossRef]
  16. Hu, C.; Lu, Y.; Sun, S.; Liu, Y. Optical Remote Sensing of Oil Spills in the Ocean: What Is Really Possible? J. Remote Sens. 2021, 2021, 9141902. [Google Scholar] [CrossRef]
  17. Jackson, C.R.; Alpers, W. The role of the critical angle in brightness reversals on sunglint images of the sea surface. J Geophys. Res.-Ocean. 2010, 115. [Google Scholar] [CrossRef]
  18. Sun, S.J.; Hu, C.M. Sun glint requirement for the remote detection of surface oil films. Geophys. Res. Lett. 2016, 43, 309–316. [Google Scholar] [CrossRef]
  19. Hu, C.M.; Feng, L.; Hardy, R.F.; Hochberg, E.J. Spectral and spatial requirements of remote measurements of pelagic Sargassum macroalgae. Remote Sens. Environ. 2015, 167, 229–246. [Google Scholar] [CrossRef]
  20. Sun, S.J.; Hu, C.M.; Tunnel, J.W. Surface oil footprint and trajectory of the Ixtoc-I oil spill determined from Landsat/MSS and CZCS observations. Mar. Pollut. Bull. 2015, 101, 632–641. [Google Scholar] [CrossRef]
  21. Lu, Y.C.; Shi, J.; Hu, C.M.; Zhang, M.W.; Sun, S.J.; Liu, Y.X. Optical interpretation of oil emulsions in the ocean—Part II: Applications to multi-band coarse-resolution imagery. Remote Sens. Environ. 2020, 242, 111778. [Google Scholar] [CrossRef]
  22. Fingas, M. Visual Appearance of Oil on the Sea. J. Mar. Sci. Eng. 2021, 9, 97. [Google Scholar] [CrossRef]
  23. Clark, R.N.; Swayze, G.A.; Leifer, I.; Livo, K.E.; Kokaly, R.; Hoefen, T.; Lundeen, S.; Eastwood, M.; Green, R.O.; Pearson, N. A method for quantitative mapping of thick oil spills using imaging spectroscopy. US Geol. Surv. Open-File Rep. 2010, 1167, 1–51. [Google Scholar]
  24. Svejkovsky, J.; Hess, M.; Muskat, J.; Nedwed, T.J.; McCall, J.; Garcia, O. Characterization of surface oil thickness distribution patterns observed during the Deepwater Horizon (MC-252) oil spill with aerial and satellite remote sensing. Mar. Pollut. Bull. 2016, 110, 162–176. [Google Scholar] [CrossRef] [PubMed]
  25. Lu, Y.C.; Shi, J.; Wen, Y.S.; Hu, C.M.; Zhou, Y.; Sun, S.J.; Zhang, M.W.; Mao, Z.H.; Liu, Y.X. Optical interpretation of oil emulsions in the ocean—Part I: Laboratory measurements and proof-of-concept with AVIRIS observations. Remote Sens. Environ. 2019, 230, 111183. [Google Scholar] [CrossRef]
  26. Jiang, Z.C.; Ma, Y.; Yang, J.F. Inversion of the Thickness of Crude Oil Film Based on an OG-CNN Model. J. Mar. Sci. Eng. 2020, 8, 653. [Google Scholar] [CrossRef]
  27. Karathanassi, V.; Topouzelis, K.; Pavlakis, P.; Rokos, D. An object-oriented methodology to detect oil spills. Int. J. Remote Sens. 2007, 27, 5235–5251. [Google Scholar] [CrossRef]
  28. Konik, M.; Bradtke, K. Object-oriented approach to oil spill detection using ENVISAT ASAR images. ISPRS J. Photogramm. Remote Sens. 2016, 118, 37–52. [Google Scholar] [CrossRef]
  29. Chen, F.; Zhang, A.H.; Balzter, H.; Ren, P.; Zhou, H.Y. Oil Spill SAR Image Segmentation via Probability Distribution Modeling. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 15, 533–554. [Google Scholar] [CrossRef]
  30. Garcia-Pineda, O.; MacDonald, I.R.; Li, X.F.; Jackson, C.R.; Pichel, W.G. Oil Spill Mapping and Measurement in the Gulf of Mexico With Textural Classifier Neural Network Algorithm (TCNNA). IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 2517–2525. [Google Scholar] [CrossRef]
  31. Ren, P.; Xu, M.; Yu, Y.H.; Chen, F.; Jiang, X.Y.; Yang, E.F. Energy Minimization With One Dot Fuzzy Initialization for Marine Oil Spill Segmentation. IEEE J. Ocean. Eng. 2019, 44, 1102–1115. [Google Scholar] [CrossRef] [Green Version]
  32. Tong, S.W.; Liu, X.G.; Chen, Q.H.; Zhang, Z.J.; Xie, G.Q. Multi-Feature Based Ocean Oil Spill Detection for Polarimetric SAR Data Using Random Forest and the Self-Similarity Parameter. Remote Sens. 2019, 11, 451. [Google Scholar] [CrossRef]
  33. Hong, X.R.; Chen, L.S.; Sun, S.J.; Sun, Z.; Chen, Y.; Mei, Q.; Chen, Z.C. Detection of Oil Spills in the Northern South China Sea Using Landsat-8 OLI. Remote Sens. 2022, 14, 3966. [Google Scholar] [CrossRef]
  34. Bayramov, E.; Bayramov, R.; Aliyeva, S. Optical and Radar Remote Sensing and Contamination Probability Modelling for the Advanced Quantitative Risk Assessment of Marine Petroleum and Gas Industry. IFAC-Pap. 2018, 51, 31–33. [Google Scholar] [CrossRef]
  35. Sun, S.; Hu, C.; Garcia-Pineda, O.; Kourafalou, V.; Le Henaff, M.; Androulidakis, Y. Remote sensing assessment of oil spills near a damaged platform in the Gulf of Mexico. Mar. Pollut. Bull. 2018, 136, 141–151. [Google Scholar] [CrossRef] [PubMed]
  36. Zhao, J.; Temimi, M.; Al Azhar, M.; Ghedira, H.; Marpu, P. Multi-Sensor Based Approach for Detection of Oil Pollution in the Arabian Gulf and the Sea of Oman. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015. [Google Scholar]
  37. Zhao, D.; Cheng, X.W.; Zhang, H.P.; Zhang, H.T. An Oil Slick Detection Index Based on Landsat8 Remote Sensing Images. In Proceedings of the International Workshop on Big Geospatial Data and Data Science (BGDDS), Wuhan, China, 22–23 September 2018. [Google Scholar]
  38. Rajendran, S.; Vethamony, P.; Sadooni, F.N.; Al-Kuwari, H.A.; Al-Khayat, J.A.; Govil, H.; Nasir, S. Sentinel-2 image transformation methods for mapping oil spill—A case study with Wakashio oil spill in the Indian Ocean, off Mauritius. MethodsX 2021, 8, 101327. [Google Scholar] [CrossRef] [PubMed]
  39. Liu, B.X.; Zhang, Q.; Li, Y.; Chang, W.; Zhou, M.R. Spatial-Spectral Jointed Stacked Auto-Encoder-Based Deep Learning for Oil Slick Extraction from Hyperspectral Images. J. Indian Soc. Remote Sens. 2019, 47, 1989–1997. [Google Scholar] [CrossRef]
  40. Ozigis, M.S.; Kaduk, J.D.; Jarvis, C.H. Mapping terrestrial oil spill impact using machine learning random forest and Landsat 8 OLI imagery: A case site within the Niger Delta region of Nigeria. Environ. Sci. Pollut. Res. 2019, 26, 3621–3635. [Google Scholar] [CrossRef] [Green Version]
  41. Taravat, A.; Del Frate, F. Development of band ratioing algorithms and neural networks to detection of oil spills using Landsat ETM plus data. EURASIP J. Adv. Signal Process. 2012, 2012, 1. [Google Scholar] [CrossRef]
  42. Kolokoussis, P.; Karathanassi, V. Oil Spill Detection and Mapping Using Sentinel 2 Imagery. J. Mar. Sci. Eng. 2018, 6, 4. [Google Scholar] [CrossRef] [Green Version]
  43. Lacava, T.; Ciancia, E.; Coviello, I.; Di Polito, C.; Grimaldi, C.S.L.; Pergola, N.; Satriano, V.; Temimi, M.; Zhao, J.; Tramutoli, V. A MODIS-Based Robust Satellite Technique (RST) for Timely Detection of Oil Spilled Areas. Remote Sens. 2017, 9, 128. [Google Scholar] [CrossRef] [Green Version]
  44. Satriano, V.; Ciancia, E.; Lacava, T.; Pergola, N.; Tramutoli, V. Improving the RST-OIL Algorithm for Oil Spill Detection under Severe Sun Glint Conditions. Remote Sens. 2019, 11, 2762. [Google Scholar] [CrossRef]
  45. Burns, K.A.; Jones, R. Assessment of sediment hydrocarbon contamination from the 2009 Montara oil blow out in the Timor Sea. Environ. Pollut. 2016, 211, 214–225. [Google Scholar] [CrossRef] [PubMed]
  46. Spies, R.B.; Mukhtasor, M.; Burns, K.A. The Montara Oil Spill: A 2009 Well Blowout in the Timor Sea. Arch. Environ. Contam. Toxicol. 2017, 73, 55–62. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, M.; Bailey, S.W. Correction of Sun glint Contamination on the SeaWiFS Ocean and Atmosphere Products. Appl. Opt. 2001, 40, 4790–4798. [Google Scholar] [CrossRef]
  48. Zhang, H.; Wang, M.H. Evaluation of sun glint models using MODIS measurements. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 492–506. [Google Scholar] [CrossRef]
  49. Cox, C.S.; Munk, W.H. Measurement of the Roughness of the Sea Surface from Photographs of the Sun’s Glitter. J. Opt. Soc. Am. 1954, 44, 838–850. [Google Scholar] [CrossRef]
  50. Liu, Y.B. Why Ndwi threshold varies in delineating water body from multi-temporal images? In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Munich, Germany, 22–27 July 2012; pp. 4375–4378. [Google Scholar]
  51. Wang, M.Q.; Hu, C.M. Extracting Oil Slick Features From VIIRS Nighttime Imagery Using a Gaussian Filter and Morphological Constraints. IEEE Geosci. Remote Sens. Lett. 2015, 12, 2051–2055. [Google Scholar] [CrossRef]
  52. Wang, X.Z.; Liu, J.X.; Zhang, S.; Deng, Q.W.; Wang, Z.; Li, Y.H.; Fan, J.C. Detection of Oil Spill Using SAR Imagery Based on AlexNet Model. Comput. Intell. Neurosci. 2021, 2021, 4812979. [Google Scholar] [CrossRef]
  53. Yang, M.M.; Zhao, P.Y.; Feng, B.; Zhao, F. Water Surface Sun Glint Suppression Method Based on Polarization Filtering and Polynomial Fitting. Laser Optoelectron. Prog. 2021, 58. [Google Scholar] [CrossRef]
  54. Cui, X.; Xin, Y. An Effective Method in the Detection of Infrared Dim Target. Acta Photonica Sin. 2014, 43, 0210003. [Google Scholar]
  55. Zheng, Y.G.; Zhang, X.R.; Hou, B.; Liu, G.C. Using Combined Difference Image and k-Means Clustering for SAR Image Change Detection. IEEE Geosci. Remote Sens. Lett. 2014, 11, 691–695. [Google Scholar] [CrossRef]
  56. Achanta, R.; Shaji, A.; Smith, K.; Lucchi, A.; Fua, P.; Susstrunk, S. SLIC Superpixels Compared to State-of-the-Art Superpixel Methods. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 34, 2274–2281. [Google Scholar] [CrossRef] [PubMed]
  57. Cui, Q.N.; Pan, H.W.; Li, X.K.; Zhang, K.J.; Chen, W.P. CMSuG: Competitive mechanism-based superpixel generation method for image segmentation. J. Intell. Fuzzy Syst. 2022, 43, 4409–4430. [Google Scholar] [CrossRef]
  58. Fu, Z.L.; Sun, Y.J.; Fan, L.; Han, Y.T. Multiscale and Multifeature Segmentation of High-Spatial Resolution Remote Sensing Images Using Superpixels with Mutual Optimal Strategy. Remote Sens. 2018, 10, 1289. [Google Scholar] [CrossRef] [Green Version]
  59. Zhijie, C.; Baolong, G.; Cheng, L.; Hongyan, L. Review on Superpixel Generation Algorithms Based on Clustering. In Proceedings of the 2020 IEEE 3rd International Conference on Information Systems and Computer Aided Education (ICISCAE), Dalian, China, 27–29 September 2020; pp. 532–537. [Google Scholar] [CrossRef]
  60. Zhang, J.; Feng, H.; Luo, Q.L.; Li, Y.; Wei, J.J.; Li, J. Oil Spill Detection in Quad-Polarimetric SAR Images Using an Advanced Convolutional Neural Network Based on SuperPixel Model. Remote Sens. 2020, 12, 944. [Google Scholar] [CrossRef] [Green Version]
  61. Malik, R. Learning a classification model for segmentation. In Proceedings of the Ninth IEEE International Conference on Computer Vision, Nice, France, 13–16 October 2003; Volume 11, pp. 10–17. [Google Scholar]
  62. Csillik, O. Fast Segmentation and Classification of Very High Resolution Remote Sensing Data Using SLIC Superpixels. Remote Sens. 2017, 9, 243. [Google Scholar] [CrossRef] [Green Version]
  63. Dong, Y.Z.; Liu, Y.X.; Hu, C.M.; MacDonald, I.R.; Lu, Y.C. Chronic oiling in global oceans. Science 2022, 376, 1300–1304. [Google Scholar] [CrossRef]
  64. Li, Y.; Lan, G.X.; Li, J.J.; Ma, L. Potential Analysis of Maritime Oil Spill Monitoring Based on MODIS Thermal Infrared Data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009. [Google Scholar]
  65. Balogun, A.; Yekeen, S.; Pradhan, B.; Althuwaynee, O. Spatio-Temporal Analysis of Oil Spill Impact and Recovery Pattern of Coastal Vegetation and Wetland Using Multispectral Satellite Landsat 8-OLI Imagery and Machine Learning Models. Remote Sens. 2020, 12, 1225. [Google Scholar] [CrossRef] [Green Version]
  66. Temitope Yekeen, S.; Balogun, A.L.; Wan Yusof, K.B. A novel deep learning instance segmentation model for automated marine oil spill detection. ISPRS J. Photogramm. Remote Sens. 2020, 167, 190–200. [Google Scholar] [CrossRef]
  67. Huang, X.; Zhang, B.; Perrie, W.; Lu, Y.; Wang, C. A novel deep learning method for marine oil spill detection from satellite synthetic aperture radar imagery. Mar. Pollut. Bull. 2022, 179, 113666. [Google Scholar] [CrossRef]
Figure 1. The Montara wellhead platform location (latitude 12.7° S, longitude 124.5° E) and the study area of this research (in the red rectangular region).
Figure 1. The Montara wellhead platform location (latitude 12.7° S, longitude 124.5° E) and the study area of this research (in the red rectangular region).
Jmse 10 01630 g001
Figure 2. Flowchart for image preprocessing.
Figure 2. Flowchart for image preprocessing.
Jmse 10 01630 g002
Figure 3. First channel of FRGB (converted from Rrc_859) from MODIS Aqua on 30 August 2009, coordinated universal time (UTC) 5:20. Result image after a series of preprocessing procedures: (a) Gaussian stretching, (b) subtraction of the polynomial fit, and (c) subtraction of the additional median filtering. The horizontal transect over the spill regions is marked as a red line in panels a, b, and c. (df) represent brightness values over the transect marked in panels a, b, and c, respectively.
Figure 3. First channel of FRGB (converted from Rrc_859) from MODIS Aqua on 30 August 2009, coordinated universal time (UTC) 5:20. Result image after a series of preprocessing procedures: (a) Gaussian stretching, (b) subtraction of the polynomial fit, and (c) subtraction of the additional median filtering. The horizontal transect over the spill regions is marked as a red line in panels a, b, and c. (df) represent brightness values over the transect marked in panels a, b, and c, respectively.
Jmse 10 01630 g003
Figure 4. Overview flow chart to detect massive oil spills in optical imagery.
Figure 4. Overview flow chart to detect massive oil spills in optical imagery.
Jmse 10 01630 g004
Figure 5. (a) MODIS Aqua FRGB image on 30 August 2009, UTC 5:20 after image preprocessing. (b) Results after the super-pixel segmentation. (ce) are the RAG mergence results with the termination thresholds set to 15, 30, and 55, respectively. The threshold of 30 was employed in this study (panel d), which was calculated based on the adaptive threshold equation (Equation (1)). If the termination threshold is set at a lower value (panel c), the objects of oil slicks would be over-segmented, resulting in too many objects (the different colors represent different objects). On the contrary, a larger termination threshold (panel e) would lead to under-segmentation, where some of the small-sized oil slicks displaying low contrast with the background would be merged into the background water category.
Figure 5. (a) MODIS Aqua FRGB image on 30 August 2009, UTC 5:20 after image preprocessing. (b) Results after the super-pixel segmentation. (ce) are the RAG mergence results with the termination thresholds set to 15, 30, and 55, respectively. The threshold of 30 was employed in this study (panel d), which was calculated based on the adaptive threshold equation (Equation (1)). If the termination threshold is set at a lower value (panel c), the objects of oil slicks would be over-segmented, resulting in too many objects (the different colors represent different objects). On the contrary, a larger termination threshold (panel e) would lead to under-segmentation, where some of the small-sized oil slicks displaying low contrast with the background would be merged into the background water category.
Jmse 10 01630 g005
Figure 6. Flow chart for recognizing bright/dark objects.
Figure 6. Flow chart for recognizing bright/dark objects.
Jmse 10 01630 g006
Figure 7. Scheme to illustrate searching for object boundary pixels and the shifted object. (ac) are used to describe the pixels used in the calculation of mean difference to cloud shadow object (MDCS), mean difference to neighborhood pixels (MDN), and mean difference to cloud object (MDC), respectively.
Figure 7. Scheme to illustrate searching for object boundary pixels and the shifted object. (ac) are used to describe the pixels used in the calculation of mean difference to cloud shadow object (MDCS), mean difference to neighborhood pixels (MDN), and mean difference to cloud object (MDC), respectively.
Jmse 10 01630 g007
Figure 8. The negative contrast oil slick discrimination process.
Figure 8. The negative contrast oil slick discrimination process.
Jmse 10 01630 g008
Figure 9. The positive contrast oil slick discrimination process.
Figure 9. The positive contrast oil slick discrimination process.
Jmse 10 01630 g009
Figure 10. (a) MODIS Terra on 3 October 2009, UTC 02:15, (b) MODIS Aqua on 3 October 2009, UTC 05:10, and (c) MODIS Terra on 19 October 2009, UTC 02:15. Images in panels a and b present the same oil slicks under different solar/viewing geometries, but with different contrasts with the water. (d) Rrc spectra of the selected points marked in panels a–c.
Figure 10. (a) MODIS Terra on 3 October 2009, UTC 02:15, (b) MODIS Aqua on 3 October 2009, UTC 05:10, and (c) MODIS Terra on 19 October 2009, UTC 02:15. Images in panels a and b present the same oil slicks under different solar/viewing geometries, but with different contrasts with the water. (d) Rrc spectra of the selected points marked in panels a–c.
Jmse 10 01630 g010
Figure 11. (aj) MODIS imagery and extraction results on different dates: FRGB image (Row 1 and Row 4), delineated oil slick (Row 2 and Row 5) from visual interpretation, and extracted slick results (Row 3 and Row 6) using the approach proposed in this study. Areas with red ellipses and squares outline the omitted oil slicks regions and the misidentified oil slicks regions, respectively.
Figure 11. (aj) MODIS imagery and extraction results on different dates: FRGB image (Row 1 and Row 4), delineated oil slick (Row 2 and Row 5) from visual interpretation, and extracted slick results (Row 3 and Row 6) using the approach proposed in this study. Areas with red ellipses and squares outline the omitted oil slicks regions and the misidentified oil slicks regions, respectively.
Jmse 10 01630 g011
Table 1. Parameter settings for data preprocessing and segmentation.
Table 1. Parameter settings for data preprocessing and segmentation.
Data preprocessingSFIT functionMaximum degree of fit4
Median filterKernel size233
SegmentationSLICSegmentation number K200,000
Maximum color distance Nc10
Maximum iteration10
Table 2. Evaluation of the extracted oil spills. Data and time are in coordinated universal time (UTC).
Table 2. Evaluation of the extracted oil spills. Data and time are in coordinated universal time (UTC).
Serial NumberDateTime
(hh:mm)
SensorArea of the Oil Spill
by the Approach (km2)
PrecisionRecallF1-Score
a30 Aug 200905:20Aqua1468.440.9690.8610.912
b17 Sep 200902:15Terra3338.060.9810.8060.885
c24 Sep 200905:15Aqua5310.130.930.8630.895
d24 Sep 200902:15Terra4055.060.9390.8730.905
e03 Oct 200902:15Terra3559.560.8480.8370.843
f10 Oct 200905:15Aqua4293.940.8770.8260.851
g17 Oct 200905:20Aqua1131.440.6390.7830.704
h19 Oct 200902:15Terra2521.190.4480.8090.576
i21 Oct 200902:05Terra2208.380.9920.780.873
j30 Oct 200901:55Terra5255.060.4230.8790.571
Overall////0.80.8380.818
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Z.; Sun, S.; Zhao, J.; Ai, B.; Yang, Q. Detection of Massive Oil Spills in Sun Glint Optical Imagery through Super-Pixel Segmentation. J. Mar. Sci. Eng. 2022, 10, 1630. https://doi.org/10.3390/jmse10111630

AMA Style

Sun Z, Sun S, Zhao J, Ai B, Yang Q. Detection of Massive Oil Spills in Sun Glint Optical Imagery through Super-Pixel Segmentation. Journal of Marine Science and Engineering. 2022; 10(11):1630. https://doi.org/10.3390/jmse10111630

Chicago/Turabian Style

Sun, Zhen, Shaojie Sun, Jun Zhao, Bin Ai, and Qingshu Yang. 2022. "Detection of Massive Oil Spills in Sun Glint Optical Imagery through Super-Pixel Segmentation" Journal of Marine Science and Engineering 10, no. 11: 1630. https://doi.org/10.3390/jmse10111630

APA Style

Sun, Z., Sun, S., Zhao, J., Ai, B., & Yang, Q. (2022). Detection of Massive Oil Spills in Sun Glint Optical Imagery through Super-Pixel Segmentation. Journal of Marine Science and Engineering, 10(11), 1630. https://doi.org/10.3390/jmse10111630

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