Next Article in Journal
An Automatic Road Network Construction Method Using Massive GPS Trajectory Data
Previous Article in Journal
A Polygon and Point-Based Approach to Matching Geospatial Features
Previous Article in Special Issue
A New Endmember Preprocessing Method for the Hyperspectral Unmixing of Imagery Containing Marine Oil Spills
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploratory Data Analysis of Synthetic Aperture Radar (SAR) Measurements to Distinguish the Sea Surface Expressions of Naturally-Occurring Oil Seeps from Human-Related Oil Spills in Campeche Bay (Gulf of Mexico)

by
Gustavo De Araújo Carvalho
1,*,
Peter J. Minnett
2,
Fernando Pellon De Miranda
1,
Luiz Landau
1 and
Eduardo Tavares Paes
3
1
Programa de Engenharia Civil (PEC), Laboratório de Métodos Computacionais em Engenharia (LAMCE), Laboratório de Sensoriamento Remoto por Radar Aplicado à Indústria do Petróleo (LabSAR), Instituto Alberto Luiz Coimbra de Pós-Graduação e Pesquisa de Engenharia (COPPE), Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, RJ 21941-909, Brazil
2
Department of Ocean Sciences (OCE), Rosenstiel School of Marine and Atmospheric Science (RSMAS), University of Miami (UM), Miami, FL 33149, USA
3
Laboratório de Ecologia Marinha e Oceanografia Pesqueira da Amazônia (LEMOPA), Instituto Socioambiental e dos Recursos Hídricos (ISARH), Universidade Federal Rural da Amazônia (UFRA), Belém, PA 66077-830, Brazil
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2017, 6(12), 379; https://doi.org/10.3390/ijgi6120379
Submission received: 30 September 2017 / Revised: 27 October 2017 / Accepted: 13 November 2017 / Published: 6 December 2017
(This article belongs to the Special Issue Oil and Gas Applications of Remote Sensing and UAV Systems)

Abstract

:
An Exploratory Data Analysis (EDA) aims to use Synthetic Aperture Radar (SAR) measurements for discriminating between two oil slick types observed on the sea surface: naturally-occurring oil seeps versus human-related oil spills—the use of satellite sensors for this task is poorly documented in scientific literature. A long-term RADARSAT dataset (2008–2012) is exploited to investigate oil slicks in Campeche Bay (Gulf of Mexico). Simple Classification Algorithms to distinguish the oil slick type are designed based on standard multivariate data analysis techniques. Various attributes of geometry, shape, and dimension that describe the oil slick Size Information are combined with SAR-derived backscatter coefficients—sigma-(σo), beta-(βo), and gamma-(γo) naught. The combination of several of these characteristics is capable of distinguishing the oil slick type with ~70% of overall accuracy, however, the sole and simple use of two specific oil slick’s Size Information (i.e., area and perimeter) is equally capable of distinguishing seeps from spills. The data mining exercise of our EDA promotes a novel idea bridging petroleum pollution and remote sensing research, thus paving the way to further investigate the satellite synoptic view to express geophysical differences between seeped and spilled oil observed on the sea surface for systematic use.

1. Introduction

Oil floating on the surface of the ocean is commonly detected using instruments flying onboard airplanes or satellites [1]. Currently, the most useful systems to detect oil-contaminated areas at sea are active radars, i.e., Synthetic Aperture Radars (SAR) operating in the microwave region of the electromagnetic spectrum [2]. As some environmental phenomena can generate signatures similar to oil in SAR imagery, this technology may yield false positives, in which the non-unique signature of oil caused by false targets can induce to ambiguous interpretations [3]. The so-called “look-alike features (or radar look-alikes)” range from atmospheric phenomena (e.g., rain cells) and oceanographic features (e.g., upwelling zones), amongst others [4]. Many scholars have described the use of SAR measurements to distinguish oil from look-alike features, and indeed, a broad scientific basis is devoted to the detection and identification of oil in SAR imagery, e.g., [5].
A supplementary application for satellite measurements is the recognition of the oil slick type. Herein, for the sake of brevity and simplicity, the term oil slick indicates the sea surface expression of oil naturally seeped out of the seafloor (i.e., oil seep) or spilled after man-made activities (i.e., oil spill). The latter accounts for operational, accidental, or illegal spills from platforms or ships—i.e., discharge of water containing petroleum products. Information about the oil slick type is an improvement requirement for several purposes, for example, environmental monitoring, emergency response, etc. [6]. Reliable surveillance based on satellite remote sensors capable of distinguishing the oil slick type can add accuracy to oil spill forecast systems, as well as enhance tridimensional mathematical models that backtrack (i.e., hindcast) oil seeps observed on the sea surface to its seafloor origin [7].
Environmental and economical issues associated with oil slicks are a constant concern to the oil and gas exploration and production industry. As such, the satellite synoptic view is an attractive option to distinguish the oil slick type that can lead to considerable scientific advances. From the standpoint of environmental monitoring programs, it can reduce ambiguities about the source of the observed oil (i.e., seeped or spilled). This can develop the relationship between governmental agencies and oil and gas exploration and production industry, thus reducing political uproar. From an economic standpoint, it can lead to offshore discoveries bringing invaluable information to explore active petroleum systems. Therefore, if the distinguishing of seeps and spills becomes possible, satellite information can directly assist in the search to find new oil fields in offshore exploration frontiers, a constant goal of a crucial sector for the world’s economic development.
Oil seeps and oil spills are both products of mineral oil floating on the sea surface, therefore, it is expected that their surface signatures are fairly similar in satellite imagery. Peer-reviewed and grey literature present limited information available about oil slick type differentiation using satellite sensors; to this matter, the reader is encouraged to see [8] and the vast list of references therein. As a result, the objective of the present research focuses on performing an Exploratory Data Analysis (EDA) to evaluate the use of SAR measurements to distinguish the sea surface expressions of oil seeps from oil spills to a useful level of confidence for systematic use.
In the course of achieving this objective, three scientific questions are sought about the sea surface expression of oil slicks observed in the Campeche Bay region (Gulf of Mexico—Figure 1), as determined by digital classification of satellite imagery:
  • Does seeped oil floating on the ocean surface have a SAR backscatter signature distinctive enough to distinguish it from anthropogenically-spilled oil?
  • Can the oil slick Size Information be used to distinguish seeps from spills?
  • Which characteristics lead to the generation of a system capable of distinguishing between seeped and spilled oil?
To address these matters, a long-term RADARSAT-2 dataset (2008–2012) is exploited to perform a data mining of selected oil slicks’ characteristics. Therefore, an additional goal is devised to design innovative qualitative–quantitative Classification Algorithms to distinguish man-made from natural oil slicks.

2. Dataset

An environmental monitoring observed the occurrence of oil slicks for more than a decade (2000–2012) in Campeche Bay—this has been carried out by Pemex (Petróleos Mexicanos; Mexico City, Mexico) to support its decision-making processes, thus demonstrating the operational usefulness SAR measurements to execute effective oil slick mapping. Such monitoring was based on the interpretation of SAR scenes of both RADARSAT satellites analyzed by domain specialists with support of ancillary meteo-oceanographic data from Earth Observing System (EOS) sensors. All logged oil slicks underwent a post-processing validation, in which internal Pemex field reports were considered—i.e., the sea surface source recognized by the satellite image analyses were corroborated by additional field observations.
The long-term RADARSAT dataset exploited by Pemex was made available for this investigation and it is herein referred to as Campeche Bay Oil Slick Satellite Database (i.e., CBOS-Data). It is good to point out that the CBOS-Data does not account for look-alike features and that oil seepage sites on the sea surface were identified with at least three different observations about the same location, thus indicating oil seep cases. For a thorough description of the Pemex’s operational environmental monitoring strategy and its produced dataset (i.e., CBOS-Data), the reader is encouraged to see a D.Sc. Thesis [8], a peer-reviewed article [9], and a conference paper [10]—these three references also provide a description of the occurrence and spatio-temporal distribution of the observed oil seeps and oil spills in the region.
To guarantee the quality of circumstances to cope with the objectives of the present research, our EDA only considers oil slicks observed between 2008 and 2012 with RADARSAT-2 (16-bit and VV-polarized—i.e., both ScanSAR Narrow beam modes: SCNA and SCNB). This mostly aims to circumvent technical specification differences (e.g., radiometric depth and polarization) between the two single-polarized RADARSAT measurements, as all analyzed RADARSAT-1 scenes are 8-bit and HH-polarized. Even though RADARSAT-2 provides a range of polarimetric information (i.e., linear cross-, dual parallel cross-, and full quad-polarized), our EDA focuses on using single linear parallel polarimetric SAR measurements from Pemex’s provided information—i.e., CBOS-Data. As such, the dataset exploited in our EDA consists of 277 RADARSAT-2 scenes and 4916 oil slicks, of which 2021 (41%) have been identified as oil seeps and 2895 (59%) as oil spills. Herein, this dataset is referred to as Campeche Bay Oil Slick Modified Database: CBOS-MOD. The entire content of the CBOS-MOD has been used to perform the proposed data mining exercises, as well as to design the innovative qualitative-quantitative Classification Algorithms, thus ensuring best statistical reliability of our EDA.

3. Methods

The Methods’ structure of our EDA is divided in two main Sections: Workable-Database Preparation (Section 3.1) and Multivariate Data Analysis (Section 3.2). Key concepts are emphasized providing enough detailed information to enable any knowledgeable scientist to replicate the results of our EDA with any SAR dataset containing identified oil slicks [8].

3.1. Workable-Database Preparation

3.1.1. RADARSAT Re-Processing

The RADARSAT-2 images, provided in SGF-format (SAR Georeferenced Fine Resolution), are given in uncalibrated grey-level count, i.e., Digital Numbers (DNs). These are enough for qualitative usage, but they are not recommended to cross-compare time series of SAR images [11,12]. Conversely, SAR backscatter coefficients are essential for quantitative analyses. The use of sigma-naught (σo), beta-naught (βo), or gamma-naught (γo) permit the comparison of data acquired with the same sensor from different dates or beam modes, as well as the evaluation of data acquired with various sensors [13,14,15]. The conversion from DN to σo, βo, and γo has been completed using the PCI Geomatica (version 2014, PCI Geomatics; Markham, ON, Canada). Various radiometric-calibrated image products are used in our EDA to describe the radar signal backscatter strength of individual oil slicks. Table 1 lists the twelve forms of σo, βo, and γo that are given with the amplitude of the received radar beam (C1) and in dB units (C2), both with and without the application of the Frost filter with 3 × 3 window [16]. A thirteenth image product gives the incidence angle values per pixel: INC.ang.

3.1.2. New Slick-Feature Attributes

Most new slick-feature attributes are bridged from elemental characteristics found in the literature; however, these have been suggested to distinguish oil slicks from look-alike features. Herein, they are explored to distinguish seeps from spills. A particular set of basic statistical measures is experimentally used to describe the oil slicks’ characteristics. Four Attribute Types (Table 2 and Table 3) are investigated:
  • 1st) Contextual Information: The Category attribute informs the oil slick type of each polygon: seeps or spills. Two other attributes provide the geographical location of the oil slick’s centroid: latitude (cLAT) and longitude (cLONG). A pair of satellite overpass attributes provides imaging time (SARtime) and observation date (SARdate).
  • 2nd) SAR Scene Descriptors: The satellite scene is described by the imaging configuration (i.e., beam mode: Bmode) and the incidence angle (INC.ang). Regarding the latter, different basic statistical measures are also calculated, as shown when the 4th Attribute Type (i.e., SAR backscatter signature) are introduced.
  • 3rd) Size Information: From the two basic morphological features originally present in the CBOS-Data, i.e., area (Area) and perimeter (Per), seven ratios are derived to describe the oil slicks’ geometry, shape, and dimension (Table 2). The first is the area to perimeter ratio (AtoP). Fiscella et al. [17] suggest using the perimeter to area ratio (PtoA). Fiscella et al. [17] and Singha et al. [18] recommend using a dimensionless normalized perimeter to area ratio (PtoA.nor). While small PtoA.nor values are related to simple geometry, larger values come from oil slicks with more complex geometries [19]. A dimensionless complexity measure is given by [20]: COMPLEX.ind. In addition, Bentz et al. [21] uses another dimensionless descriptor to illustrate how compact (i.e., close to a circle) is a sea-surface feature: COMPACT.ind. Two other indices have been utilized by [22]: SHAPE.ind and FRACTAL.ind. While these indices yield values close to the unit for regular forms (i.e., circular or square, respectively), larger numbers represent form irregularity [23]. The total number of pixels inside each oil slick polygon (LEN) is also provided.
  • 4th) SAR Backscatter Signature: All pixels within individual polygons are utilized to calculate the different basic statistical measures experimentally used to characterize the oil slick SAR backscatter signature (Table 3). These are separately calculated for the twelve radiometric-calibrated image products (Table 1). As suggested by [18,21], an arithmetic mean (AVG) of all pixels of each oil slick is computed. Three other central tendency representations are also considered: median (MED), mode (MOD), and mid-mean (MDM). To analyze the SAR backscatter signature spread, six dispersion measures are considered: standard deviation (STD), coefficient of dispersion (COD), variance (VAR), total range (RNG), average absolute deviation (AAD), and median absolute deviation (MAD). The coefficient of variation (COV) evaluates the relative relationship between dispersion and central tendency, thus comparing the degree of variation of data with different units and different meanings. COV originally involves the ratio between STD and AVG, but herein other pair-values are also given (Table 3). Different authors recommend exploring COV, e.g., [18,24] refer to it as power-to-mean ratio, whereas [20] describe it as the oil slick’s homogeneity, and for [21] it depicts the oil slick’s heterogeneity. However, Miranda [25] emphasizes that oil slicks with different spatial structures (i.e., with completely different pixel configuration) can have identical average or median values and the same standard deviation. Minimum (MIN) and maximum (MAX) values of the pixels inside each polygon are also used as supplementary quantities.
Information from the area surrounding oil slicks may play an essential role in the classification between oil slicks and look-alike features [17,20,26]. Nonetheless, herein, no information from the background oil-free surface around oil slicks has been taken into consideration. This decision aims to evaluate as simple as possible range of variables solely accounting for information within the oil slicks polygons’ limits.

3.1.3. Data Treatment Processes

The new slick-feature attributes (Table 2 and Table 3) have both non-numeric values (i.e., qualitative variables) and numeric values (i.e., quantitative variables). The multivariate data analysis techniques applied herein (see Section 3.2.) requires variables to be quantitative descriptors. To guarantee that all slick-feature attributes have discrete or continuous numerical metric values with no categorical variables, four Data Treatment Processes are applied:
  • Log10 Transformation: Various non-linear transformations (e.g., square, cubic, or fourth root, as well as natural and base 10 logarithms) were tested to bring non-symmetric distributions to, or at least close to, a normal distribution pattern. Through the visual analyses of their histograms, Log10 had best results to stabilize the data, thus decreasing the effect of outliers and reducing skewness. The exception is FRACTAL.ind with its negative to positive range that was transformed as the cubic root.
  • Negative Values Scaling: Because log functions cannot be applied to negative values, all pixels inside the oil slicks were inspected before such transformation. This mostly relates to the SAR Backscatter Signature (4th Attribute Type: Table 3), as the other Attribute Types all have positive values (except FRACTAL.ind). This specific Data Treatment Process consists of subtracting the minimum negative pixel value within the oil slick (PIXmin) from each pixel inside the oil slick (PIX) and adding one to it: PIXpos = [PIX − PIXmin + 1], where PIXpos is the new positive value. This simply adds two constants to each pixel, i.e., PIXmin and 1. The pixel with the minimum negative value becomes equal to 1, as PIX = PIXmin gives PIXpos = 1. Although this changes the values of the pixels inside some oil slicks, the relationship among pixel values within individual oil slicks is preserved [27,28].
  • Ranging Standardization: While some quantitative attributes are dimensionless (e.g., dB), others have incompatible units. To compare the oil slick characteristics in the subsequent data mining exercise of our EDA, a common scale is recommended [29,30,31]. Of the various available methods, there is no dearth of standardization; however, Milligan and Cooper [32] advocated that Ranging could be more effective than other methods (e.g., z-score). The Ranging standardization, besides bounding the magnitude of the attributes between zero (0) and one (1), also equalizes the variability of the attributes to a common scale: Xranged = [X − Xmin]/[Xmax − Xmin], where Xranged is the new Ranged value, X is the actual attribute value for a particular oil slick, Xmin and Xmax are, respectively, the minimum and maximum values of this attribute among all oil slicks [27]. For instance, if taking the SIG.amp.AVG value of one oil slick: subtract the minimum AVG for this radiometric-calibrated image product among all 4916 oil slicks, then the result is divided by the maximum AVG of all 4916 oil slicks minus the minimum AVG of all 4916 oil slicks. The Xranged variables assume non-negative values, and only one oil slick has Xranged = 0 and another one Xranged = 1, respectively when X = Xmin and when X = Xmax.
  • Dummy Variables: Certain qualitative attributes, such as string variables (Category and Bmode) and spatio-temporal variables (SARtime, SARdate, cLAT, and cLONG), have been converted to a specific indicator type: “dummy variable” [30]. For convenience, these have values of one (1) or zero (0), thus referring to its presence (Yes) or absence (No), respectively. These types of binary-coded variables attempt to decompose complex variabilities (or hidden knowledge) into more useful information.

3.2. Multivariate Data Analysis

3.2.1. Attribute Selection Methods

The attributes’ relevance is explored for three reasons: to eliminate redundant (or irrelevant) variables, to select more representative variables, and to reduce the variable-hyperspace dimension. A lower degree of dependence from one attribute to another reduces the messiness of the subsequent data mining exercise of our EDA (see Section 3.2.). Two Attribute Selection Methods (i.e., R-mode), having variable selection strategies different in essence, are explored mostly based on their simplicity:
  • Unweighted Pair Group Method with Arithmetic Mean (UPGMA): The resemblance between variables is revealed on a semi-automatic manner, in which an equally weighted pair wise group relationship among variables is used to form groups [27]. A variable is attributed to a group that has a larger similarity measure, i.e., Pearson’s r correlation coefficient, with all other variables of that same group [30,31]. A free scientific data analysis software package (PAST: PAleontological STatistics version 2.17c; Oslo, Norway) was utilized [33]. Rooted-tree dendrograms (i.e., diagrams of relationships) with a bootstrapping of 100 replicates help to find hierarchical relationships (i.e., affinity) among variables. Due to the large subjectivity while using dendrograms to form groups, two user-defined thresholds are explored [34]. Their choice aims to guarantee repeatability, as different thresholds give different groups of variables. The first threshold subjectively uses the Cophenetic Correlation Coefficient (CCC) that is automatically calculated in PAST to specify the clusters; its value varies depending upon dataset [35,36]. The second threshold is simply the fixed similarity value of 0.5. These two thresholds are used to draw horizontal lines (i.e., phenon lines) across the dendrograms, and when such lines cross a branch, a group is formed; from each group of similar variables only one variable is selected [37,38]. The Attribute Selection occurs twice for each UPGMA implementation, one for each phenon line, i.e., CCC and 0.5. This variable selection gives preference to an arbitrary user-defined strategy based on the simplicity of their meanings and on their calculation form, in which simple variables are preferred over complicated ones—e.g., basic variables in lieu of ratios (e.g., Area versus AtoP), C1 (amplitude) compared with C2 (dB), SAR backscatter signature without rather than with Frost filter, and central tendency is preferable to dispersion metrics. The variables choice follows the order presented on Table 1, Table 2 and Table 3.
  • Correlation-Based Feature Selection (CFS): The CFS Attribute Selection is a fully automated process based on a heuristic variable selection strategy [39,40]. The “Merit” of various groups of selected variables is evaluated to select attributes with low inter-correlation but highly correlated to the categories being distinguished [41,42]. In essence, a Merit is calculated as a measure of the usefulness of the selected best possible group of attributes. A scientific machine learning open source package (WEKA: Waikato Environment for Knowledge Analysis version 3.6.12; Waikato, New Zealand) was utilized. The user specifies an evaluation function and a search algorithm [43,44,45]. The former is the method by which the groups of attributes are evaluated (Backward Sequential Selection: CfsSubsetEval), and the latter improves the evaluation function (best-first searching strategy: BestFirst).
A summary of the analyzed 11 original sets, as well as the 33 optimal subsets of selected variables, i.e., UPGMA (CCC and 0.5) and CFS, is presented in Figure 2. These are collectively referred to as 44 data sub-divisions. It is expected that these particularly fruitful attribute-wise evaluations reveal hidden complexities into information capable of distinguishing the oil slick type.
Although literature suggests not using images given in DNs to compare SAR images [11,12,46], herein, an assessment explores DNs to quantify the consequences of using DNs. Only Frost filter and Antenna Pattern Compensation (APC) were applied. DN signatures are expressed by the same set of basic statistical measures calculated experimentally used to describe SAR backscatter signature (4th Attribute Type: Table 3), with the same Data Treatment Processes.

3.2.2. Principal Components Analysis (PCA)

PCA reduces the variable-hyperspace dimension, helping on the interpretation of large multivariate datasets. It also assures the fulfillment of the Discriminant Analysis pre-requisite by avoiding multi-collinearity among the explored hypothetical variables, i.e., Principal Components (PCs). PAST (version 3.06; Oslo, Norway) was utilized to run the PCAs, which were applied to the 44 data sub-divisions.
Even though some guidelines are available to decide the number of PCs to account for, there is no “rule of thumb” to determine the best “stopping rule” [47,48]. Two cut-offs were explored herein [49]. The first one is simply referred to as “Scree Plot”, as it analyzes the scree plot, but the PC-selection occurs based on the random model curve of expected eigenvalues, i.e., “broken stick”. The broken stick is also plotted on the scree plot and PCs falling below this curve, i.e., to its right side, are not considered [50,51]. Furthermore, if the broken stick curve crosses the bootstrapping eigenvalue error bars, i.e., is inside the 95% confidence interval, this PC is also not considered [30,52]. The second is a simpler and direct method based on the so-called Kaiser-Guttman criterion, thus referred to as “Kaiser Criterion”, which discards any PC after a specific lower bound: eigenvalue of 1 [49,53]. Only one of these two cut-offs was used to select meaningful PCs. This occurred based on the analysis of the global picture of the percentages of variance concentrated on the PCs (see [8]). The score values of the selected meaningful PC were subjected to the subsequent data mining exercise of our EDA.
Loadings expressing the relationship between variables (rows) and PCs (columns) verify the importance (i.e., meaning) of the original variables on each selected PC, thus determining if there were variables with higher weight (i.e., loadings) influencing each PC. Bootstrapping error bars were plotted per variable, and those having bars not crossing the abscissa are deemed to influence the PC [52]. Peres-Neto et al. [47] suggest 1000 row-wise bootstrapping replicates for the eigenvalues 95% confidence interval; however, as some PCAs took more than 10 h to run, only 500 replicates were used.

3.2.3. Correlation Matrix

Correlation matrices have been used to verify residual inter-variable correlation that might still be present after the Attribute Selection and PCA. These matrices provide the parametric Pearson’s r correlation coefficient, p (uncorrelated), to determine the relationship among all original values of the CBOS-MOD variables (no PCA applied), as well as among the scores of the selected PCs (with PCA applied).

3.2.4. Discriminant Function

Discriminant Analyses are linear combinations among the selected characteristics that seek the highest discriminating power to minimize the probability of miss-discrimination [31,49]. Discriminant Functions are mathematical features that arrange objects among reported groups [33], in which the dependent variable is the oil slick type (i.e., seep or spill). Each analyzed sample is a member of only one oil slick type, i.e., simple binary discrimination. In PAST (version 2.17c; Oslo, Norway), Discriminant Functions were tested using the 44 data sub-divisions as input.
The accuracy of the Discriminant Functions is assessed via Confusion Matrices, i.e., simple 2-by-2-tables, in which one can quantify misclassification chances by using a cross validation estimate [27,54]. Table 4 illustrates adapted Confusion Matrices following a wide-open character design proposed by [55,56]. Even though the Confusion Matrix is a well-known tool for analyzing classification performance, a two-fold reason is given to introduce basic details about its usage. Firstly, usually the nomenclature of these metrics may vary, e.g., [57,58,59], etc. Secondly, usually only the Overall Accuracy is reported, but [60] demonstrated that this may lead to a quality gap on the information provided that can mislead the user [61,62]. To this matter, various metrics should come into play to report the effectiveness of the Discriminant Functions. As a starting point, one should answer a first set of four questions focusing on the lines of Table 4, frequently referred to as Producer’s Accuracy (PA) and Commission Error (CE). A second set of four questions is given to interpret the columns of Table 4, converging to what is usually called as User’s Accuracy (UA) and Omission Error (OE).

3.2.5. Oil Slick Classification Algorithm

The Classification Algorithm design implies not only the understanding of the Overall Accuracy (i.e., Confusion Matrix diagonal), but it also requires a full line- and column-wise inspection of Table 4, i.e., Sensitivity-Specificity balance and a trade-off between the Positive and Negative Predictive Values. The former is obtained by looking into lines (PA/CE), whereas the latter is achieved exploring columns (UA/OE). Steadiness is sought among line and column information. This is because lines give the reference frame of what was previously known (i.e., spills and seeps): how many of the known oil slick samples are correctly (or incorrectly) identified? Conversely, the column reference frame indicates how many oil slick samples identified by the algorithm are correctly (or incorrectly) identified? While the first question (about lines: PA/CE) measures the success of identifying known samples, the second question (about columns: UA/OE) provides a metric of the success of the algorithm to identify oil slick samples.

4. Results

The Results presented herein mirror the organization structure of the Methods.

4.1. Workable-Database Preparation

4.1.1. RADARSAT Re-Processing

Each processed RADARSAT-2 scene is about 27 GB, as it accounts for several image products (Table 1), and together, the 277 scenes occupy almost 7.5 TB of disk space.

4.1.2. New Slick-Feature Attributes

Typical values (i.e., represented by a set of basic qualitative-quantitative statistics: minimum, maximum, average, and standard deviation) of the slick-feature attributes are presented in [8], along with histograms with their frequency distributions.

4.1.3. Data Treatment Processes

The skewness of most data transformed data attributes is reduced after the application of the Negative Values Scaling and the Log10 Transformation. The exception is the FRACTAL.ind, which has a cubic root transformation applied to it, and cLAT and cLONG, which did not change much in comparison to the frequency distributions prior to the Data Treatment Processes.
The application of these processes has five main consequences: (1) The Negative Values Scaling gives all pixels inside individual oil slicks that had negative pixel values a positive value; (2) The MIN attribute of 9 radiometric-calibrated image products lost its meaning and were removed, as PIXpos of several oil slicks became 1, as PIX = PIXmin gives PIXpos = 1; (3) The frequency of distribution of all variables has been brought to a Gaussian distribution with the Log10 Transformation; (4) All qualitative variables are coded with values lying between 0 and 1 after the Ranging Standardization; and (5) The 7 qualitative attributes have been replaced by 33 new dummy variables that are binary-coded to 1 or 0.
Consequently, the CBOS-MOD that had 485 different oil slick descriptors (Table 2: 53 and Table 3: 432) has been further refined after these processes, thus accounting for 502 variables: 485 − 9 − 7 + 33.

4.2. Multivariate Data Analysis

4.2.1. Attribute Selection Methods

Table 5-1 presents the CCC (UPGMA) and Merit (CFS) values of the 11 original sets, whereas Table 5-2 reviews the number of selected variables of the 33 optimal subsets. Even though CCC values are quite similar, there is a small variation on the separate investigations of σo, βo, and γo (with and without Size Information) that gives a slight variation in their number of selected variables. A higher CCC value stands out when only Size Information is considered: 0.9143.
The CFS-Merit decreasing trend reveals important information about the oil slick type predictability (Table 5-1). The CFS implementation using the all CBOS-MOD variables (502) mostly selected categorical contextual dummy variables (13 out of 15) and its Merit is the highest as possible: 1.00. The considerable CFS-Merit drop (0.131) when such variables are removed (i.e., all SAR backscatter signature with Size Information: 433) depicts the complexity of the exploratory nature of this research: the use of radiometric and size variables to distinguish seeps from spills is indeed not an easy task.
If Size Information variables are involved, there is a higher CFS-Merit as compared to when they are not considered, e.g., all SAR backscatter signatures with (0.131) and without (0.103) size. The same pattern is observed on the separate analysis of each SAR backscatter signature, e.g., σo with size (0.127) versus only σo (0.099). When using only Size Information the CFS-Merit (0.112) is higher than when separately exploring the SAR backscatter signature (σo, βo, and γo): 0.099, 0.099, and 0.097, respectively. These are strong indications that the sole use of radiometric variables to differentiate the oil slick type is indeed intricate, as well as suggests that by using the Size Information one has a somewhat better chance in differentiating seeps from spills. The CFS-Merit resulted from the analysis of DNs showed the worse value (0.059), thus agreeing with literature that does not recommend the use of DN values to cross-compare time series of SAR images [11,12,46].
The CFS selected attributes are presented on Table 6 (separate analyses of Size Information and DN variables) and on Table 7 (each SAR backscatter coefficients with and without Size Information). A series of other tables with the results of the several CFS implementations is presented in [8].
The results of the dendrogram analysis are depicted in Figure 3. Figure 3 (1st row left) illustrates the UPGMA dendrogram exploring only Size Information: of the initial 10 attributes (Table 2), the CCC (0.9143) forms six groups and the fixed similarity gives five groups: 6 and 5 variables, respectively (Table 6). Similarly, Figure 3 (1st row right) presents assessment of the DN values: of the initial 35 attributes, the CCC (0.8020) and the fixed similarity thresholds found 5 and 3 groups (variables), respectively (Table 6). Figure 3 also depicts the dendrograms of the separate UPGMA analysis of σo (CCC = 0.8262), βo (CCC = 0.8031), and γo (CCC = 0.8240) together with Size Information, whereas Table 7 presents an inventory with the UPGMA selected attributes. Carvalho [8] presents other dendrograms and lists the selected attributes of the 44 data sub-divisions.

4.2.2. Principal Components Analysis (PCA)

A comparison of Table 5-2,3 confirms the PCA dimensionality reduction. After analyzing the global picture of the percentages of variance concentrated on particular PCs among the 44 data sub-divisions a user-defined criterion was used to decide which of the two cut-offs should be used to select the meaningful PCs: if the “Scree Plot” PC concentrated less than 80% of the variance, the “Kaiser Criterion” PC was selected; otherwise, the Scree Plot PC was the chosen cut-off (see [8]). The number of selected PCs varies per data sub-divisions (Table 5-3).
Figure 4 (left panels) demonstrates the Scree Plot PC selection strategy that explores the broken stick analyzing only the Size Information. The 4th row of Figure 4 depicts the bootstrapping eigenvalue error bars touching the broken stick curve for the PCA of the CFS; as such, no PC has been selected with the Scree Plot. Although the Kaiser Criterion selected one PC concentrating 90% of the variance, the Discriminant Function requires at least 2 variables (i.e., in this case, two PCs); therefore, the first two PCs have been selected (Table 5-3). Carvalho [8] illustrates the full proficiency of the broken stick PCs selection strategy using the Scree Plot.
Figure 4 (middle panels) exemplifies scatterplots illustrating the relationship between the first two PCs resulted from the PCA exploring only Size Information: 10, 6, 5, and 3 variables (Table 2 and Table 6). It is evident the considerable overlap between the scores of the two populations: seeps (green circles) and spills (red triangles). The same pattern is observed on all other scatterplots, not only for the relationship between the first two PCs, but also for all other combinations of PCs [8].
The top three rows of Figure 4 (middle panels) reveal a peculiar gap in the cloud of data points affecting both populations, i.e., spills and seeps. This is mostly due to the FRACTAL.ind and its bi-modal (multi-modal) frequency distribution [8]. Such gap does not exist in the CFS analysis (Figure 4: 4th row middle panel), in which FRACTAL.ind is not accounted for (Table 6). However, this gap is not as prominent on the scatterplots when the FRACTAL.ind is accounted for on the analysis of Size Information together with the original set of SAR backscatter signature variables, i.e., σo, βo, or γo. This is probably because, in this case, each PC is influenced by a much larger number of variables: 433 or 151 (Table 1 and Table 3).
Figure 4 (right panels) depicts the influence (i.e., meaning) of each variable to each PC quantified by the inspection of Loadings Plots considering only Size Information: 10, 6, 5, and 3 variables (Table 2 and Table 6). Almost all variables largely influence the first PC, the exception is on the CFS (Figure 4: 4th row right) that has the perimeter (Per) being the sole variable to account for the variance on the first PC; the second PC has the opposite occurring. In fact, the only PC of all PCAs that had a preferable influence was indeed this case, as the loadings expressions of all PCs of the 44 data sub-divisions have shown that each PC receives the influence of every single variable [8].

4.2.3. Correlation Matrix

The correlation matrices were computed in PAST (version 2.17c; Oslo, Norway). Just about all information from the original values of the CBOS-MOD variables (no PCA applied: Table 5-2) showed residual inter-variable correlation among the variables of the 44 data sub-divisions. On the other hand, and as expected, the information from the scores of the selected PCs (with PCA applied: Table 5-3) presented no inter-variable correlation. This confirms that only the PCA information is carried to design the Classification Algorithms.

4.2.4. Discriminant Function

Table 8 summarizes the Discriminant Analyses’ findings, where it is possible to notice the correspondence between Hotelling’s t2 and Overall Accuracy. The former represents a measure of equality of the means of the groups been discriminated (i.e., spills and seeps), in which larger values correspond to better discrimination of the a priori informed groups. Likewise, the latter is the total classification efficiency that reports the success rate in correctly telling apart a priori known spills and seeps samples. Initially, only the Overall Accuracy is explored; see Section 4.2.5 for the full assessment of the other metrics. The Overall Accuracy somewhat matches the Hotelling’s t2 effectiveness. Because the results about the complete exploration of the CBOS-MOD accounted for many categorical contextual dummy variables, its discrimination power almost reached 100%. This corroborates the quality of excellence of the explored dataset, i.e., CBOS-MOD. DN is again the worst, matching expectations [11,12,46]. A slightly better discrimination is observed while using only the Size Information.
Higher Hotelling’s t2 is observed for the complete exploration of the CBOS-MOD that had undergone a major data reduction (Table 5-2: 502, 59, and 46 variables, and Table 5-3: 10, 21, and 18 PCs). This was anticipated as many categorical contextual dummy variables are accounted for. On the other hand, DNs (e.g., UPGMA-CCC: 163.0) has the lower Hotelling’s t2 by far. The same result observed for its CFS-Merit (Table 5-1: 0.059); important to note that Hotelling’s t2 and CFS-Merit are not comparable.
Column-wise, the Hotelling’s t2 behavior is fairly similar to the one from the CFS-Merit (Table 5-1). The pattern shown by the sole analysis of Size Information presents higher Hotelling’s t2 as compared to when these attributes are not included. Line-wise, higher Hotelling’s t2 values are usually observed for the UPGMA-CCC optimal subset. Again, the sole use of Size Information has about the same effectiveness to discriminate seeps from spills, as when these attributes are combined with the SAR backscatter signature. This is an indication that the SAR backscatter signature is not adding much to the whole differentiation process. The UPGMA-Fixed optimal subsets tend to present lower Hotelling’s t2, line-wise speaking. With the UPGMA-CCC optimal subset, noted above, it is possible to confirm that the number of variables has some influence on the discrimination process. This might be related not only to the amount of variables, but also to their quality, i.e., whether there is redundant correlated information or not.
The UPGMA-CCC Hotelling’s t2 values are somewhat equivalent to those of the CFS. Even though their variables are similar, they are not the same (e.g., Table 6 and Table 7). This also corroborates the importance of reducing the dimensionality on the attribute-domain, for instance, by using Attribute Selection and PCA. However, these processes should be applied with caution so as not to lose information (i.e., underestimation) or to include noise (i.e., overestimation) [47,48,49].

4.2.5. Oil Slick Classification Algorithm

The full set of metrics shown in Table 4 is assessed to design the qualitative-quantitative Classification Algorithms, which are simple-to-use and simple in their formulation. Because the UPGMA-Fixed showed lower discriminating power than that of UPGMA-CCC and CFS on the initial Overall Accuracy analysis (Table 8), only the latter two are presented.
The Discriminant Function results using selected PCs from all SAR backscatter signature variables (Table 3: σo, βo, and γo) analyzed with and without Size Information (Table 2) are presented in Table 9. A two-fold encouraging outcome is revealed: (1) The influence of Size Information on the discrimination between spills and seeps: if comparing the top with the lower tables, better results are observed when these variables are accounted for (top) than when they are not (lower); and (2) It confirms that the Overall Average (64%) may not inform the real discrimination performance; for instance, Specificity (59%) and Positive Predictive Value (54%) are not very effective, thus agreeing with [60] that the user should be aware of the real accuracy of the algorithm by examining different metrics.
The results of the separate analyses of σo variables with and without Size Information are presented in Table 10. The improvement observed when Size Information is considered is quite similar to the one observed in Table 9. This reinforces the fact that there is a benefit to using these variables, as all metrics evaluated are improved when they are accounted for. Confusion Matrices for βo and γo with similar results are presented in [8].
The Discriminant Functions results for separate analysis of Size Information (best) and of DNs variables (worst) are presented in Table 11. A two-fold positive outcome is evident: (1) The relevancy of Size Information to distinguish oil seeps from oil spills on the sea surface of the Campeche Bay region; and (2) The corroboration with literature (e.g., [11,12,46] that DN images should not be used to compare time-series of SAR imagery.
Because the Multivariate Data Analysis revealed the use of Size Information (i.e., 3rd Attribute Type: Table 2) as best to distinguish the oil slick type, a separate analysis explores the only two basic morphological features originally present in the CBOS-Data: area and perimeter. The same Data Treatment Process has been applied preceding the PCA, and two PCs have been selected prior the Discriminant Function. Table 12 presents quite similar results to those shown in Table 9, Table 10 and Table 11, thus demonstrating that the sole and simple use of area and perimeter is equally capable of distinguishing human-related oil spills from naturally-occurring oil seeps to a useful accuracy. Therefore, there is no need to evoke more complicated variables for distinguishing the oil slick type.

5. Discussion

Even though the three SAR backscatter coefficients can easily be transferred from one to another [8], the pool of attributes selected by the UPGMA and CFS implementations indicate the utility of σo, βo, and γo is not exactly the same (Table 7). It is also true that because σo is directly related to the reflectivity of the ground horizontal range plane, it is a common SAR backscatter coefficient considered in oil slick detection; see [8] and references therein. Nonetheless, the data mining exercise of our EDA has demonstrated that σo, βo, and γo present differences among each other (Table 7) and that it is worth analyzing them separately.
Two aspects agreed by the remote-sensing scientific community regarding the use of geometry, shape, and dimension attributes in the use of SAR measurements to identify oil slicks are: (1) Ship spills have distinct Size Information from other oil spills [63]; and (2) These attributes are crucial to distinguish oil slicks from look-alike features [64,65]. Thus, considering the exploratory nature of our Classification Algorithms, it should be emphasized that only a small fraction of the CBOS-Data (3.9%) and of the CBOS-MOD (3.2%) represents ship spills. In addition, the sole use of only two basic oil slicks’ morphological features (i.e., area and perimeter) is an innovative way to distinguish oil from oil, i.e., natural from man-made oil slicks—a task poorly documented in the scientific literature—peer-reviewed, as well as grey literature; see [8] and references therein.
The data mining exercise of our EDA promotes a novel idea bridging petroleum pollution and remote sensing research, thus paving the way for further investigations using the satellite synoptic view to express geophysical differences between spilled and seeped oil observed on the sea surface. Despite the substantial amount of work performed in the data mining exercise of our EDA that uses standard multivariate data analysis techniques applied to SAR measurements for the successful distinguishment of sea surface expressions of naturally-occurring oil seeps from human-related oil spills in the Campeche Bay region [8], a list of nine main suggestions that may be further explored to achieve improved accuracies for systematic use is given:
  • Information from the area surrounding the oil slicks, i.e., background oil-free surface, could come into play, e.g., damping ratio [66];
  • Information about radar beam incidence angles and environmental configurations (e.g., wind conditions) could add knowledge to the differentiation process;
  • Fractal dimension (e.g., box-counting or dynamic fractal approaches) could be used to measure, analyze, and classify the oil slicks textures and shapes [67,68];
  • The number of individual non-contiguous parts forming oil slick polygons could be explored along with other contextual properties, such as the distance between centroid and the shoreline or the distance from a possible point source, i.e., oil rig complex or oil seep site location on the sea surface [21];
  • Other methods instead of Ranging Standardization could be applied during the Data Treatment Process, e.g., z-score [30];
  • As an alternative to the proposed automatic Attribute Selection Methods to select groups and variables, visual analyses of the UPGMA dendrograms could be instituted, or even other custom methods, e.g., Ward’s [30];
  • Different “stopping rule” methods could select meaningful PCs [47,48,69];
  • Artificial Neural Networks (ANN), largely used to distinguish between oil slicks and look-alike features (e.g., [70,71,72]), indeed could be tested to differentiate oil seeps from oil spills using variables from the Attribute Selection Methods used herein, e.g., area and perimeter; and;
  • Clustering Analyses (e.g., K-means) could also reveal natural groupings among oil slicks, thus supporting the oil slick type differentiation process.
In addition to the abovementioned, the data mining exercise of our EDA surely serves as an archetype for developing other ways to differentiate the oil slick type [73]. It also paves the way for full polarimetric investigations that could indeed search other ways to differentiate oil seeps from oil spills using SAR measurements (e.g., [74,75]).

6. Conclusions

The data mining exercise of our EDA has successfully demonstrated that it is feasible to use multivariate data analysis techniques, applied to SAR measurements, to distinguish the sea surface expressions of naturally-occurring oil seeps from human-related oil spills observed in Campeche Bay. Indeed, the use of an assorted pool of standard multivariate data analysis techniques (e.g., R-mode Correlation, Principal Components Analysis, and Discriminant Function) has been effectively applied to analyze the remote sensing library of a specific long-term dataset (2008–2012) that consists of 277 RADARSAT-2 scenes and 4916 oil slicks.
The answers to the three scientific questions about the sea surface expression of oil slicks are:
  • Yes. The seeped oil has SAR backscatter signatures (i.e., σo, βo, and γo) distinctive enough to distinguish it from anthropogenically-spilled oil;
  • Yes. The Size Information (i.e., attributes describing the geometry, shape, and dimension of the oil slicks) can be used to distinguish seeps from spills, and in fact, the sole and simple use of area and perimeter can also distinguish natural from man-made oil slicks with an overall accuracy of about 70%; and;
  • The synergistic combination of various oil slick characteristics leads to systems capable of distinguishing between seeped and spilled oil, in which different combinations of variables promote similar differentiation; however, the one leading to the most effective classification is represented by the sole and simple use of area and perimeter (70%), whereas the one with the worst capabilities have variables expressed in Digital Numbers (DNs).
The proposed oil slick type Classification Algorithms of our EDA have useful accuracies, and range from uncomplicated algorithms combining the synergy of different variables to very simple ones using only SAR backscatter signature or only Size Information. In fact, the effective differentiation between oil seeps and oil spills achieved using only basic morphological features, such as the area and the perimeter of oil slicks, as determined by digital classification of satellite imagery, is a relevant outcome directly contributing to the remote sensing research, as well as to the activities of the oil and gas exploration and production industry. The field of environment monitoring also benefits from the outcomes of the data mining exercise conducted during our EDA. Conclusively, if the techniques presented herein are capable of distinguishing two types of oil observed on the sea surface (i.e., naturally-occurring oil seeps from human-related oil spills) based on RADARSAT-2 measurements, it is likely they can also produce good results if applied to the problem of discriminating oil from look-alike features (for instance, algal blooms or low wind zone signatures), indeed, solely using SAR measurements.

Acknowledgments

We are thankful to Pemex and to MDA Geospatial Services for making available the high-quality RADARSAT dataset, as well as to the three anonymous reviewers that greatly contributed to improve the manuscript with thoughtful criticism. The Human Resource Program (PRH) of the National Agency of Petroleum, Natural Gas, and Biofuels (ANP) provided the financial support for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jackson, C.R.; Apel, J.R. Synthetic Aperture Radar Marine User’s Manual; NOAA/NESDIS; Office of Research and Applications: Washington, DC, USA, 2004. Available online: http://www.sarusersmanual.com (accessed on 13 November 2017).
  2. Chan, Y.K.; Koo, V.C. An introduction to synthetic aperture radar (SAR). Prog. Electromagn. Res. B 2008, 2, 27–60. [Google Scholar] [CrossRef]
  3. Espedal, H.A.; Johannessen, O.M.; Knulst, J. Satellite detection of natural films on the ocean surface. Geophys. Res. Lett. 1996, 23, 3151–3154. [Google Scholar] [CrossRef]
  4. Johannessen, O.M.; Espedal, H.A.; Jenkins, A.J.; Knulst, J. SAR surveillance of ocean surface slicks. In Proceedings of the Second ERS Application Workshop, London, UK, 6–8 December 1995; pp. 187–192. [Google Scholar]
  5. Topouzelis, K. Oil spill detection by SAR Images: Dark formation detection, feature extraction and classification algorithms. Sensors 2008, 8, 6642–6659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Mchugh, S.L. Satellite Synthetic Aperture radar in the Prosecution of Illegal Oil Discharges. Master’s Thesis, Memorial University of Newfoundland, St. John’s, NL, Canada, 2009; p. 122. [Google Scholar]
  7. Mano, M.F.; Beisl, C.H.; Siqueira, C.Y.S.; Pereira, J.S. Evaluation of remote technologies applied to natural seep mapping and their impact in oil exploration. In Proceedings of the Rio Oil & Gas Expo and Conference, Rio de Janeiro, Brazil, 15–18 September 2014; p. 7. [Google Scholar]
  8. Carvalho, G.A. Multivariate data Analysis of Satellite-Derived Measurements to Distinguish Natural from Man-Made Oil Slicks on the Sea Surface of Campeche Bay (Mexico). Ph.D. Dissertation, Universidade Federal do Rio de Janeiro (UFRJ), COPPE, Rio de Janeiro, Brazil, 2015; p. 285. Available online: http://www.coc.ufrj.br/index.php?option=com_content&view=article&id=5354:gustavo-de-araujo-carvalho (accessed on 13 November 2017).
  9. Carvalho, G.A.; Minnett, P.J.; Miranda, F.P.; Landau, L.; Moreira, F. The use of a RADARSAT-derived long-term dataset to investigate the sea surface expressions of human-related oil spills and naturally-occurring oil seeps in Campeche Bay, Gulf of Mexico. Can. J. Remote Sens. 2016, 42, 307–321. [Google Scholar] [CrossRef]
  10. Carvalho, G.A.; Landau, L.; Miranda, F.P.; Minnett, P.; Moreira, F.; Beisl, C. The use of RADARSAT-derived information to investigate oil slick occurrence in Campeche Bay, Gulf of Mexico. In Proceedings of the XVII Brazilian Remote Sensing Symposium (SBSR), INPE, João Pessoa, Brazil, 25–29 April 2015; pp. 1184–1191. Available online: http://www.dsr.inpe.br/sbsr2015/files/p0217.pdf (accessed on 13 November 2017).
  11. Freeman, A. Radiometric calibration of SAR image data. In Proceedings of the XVII Congress for Photogrammetry and Remote Sensing (ISPRS), Washington, DC, USA, 2–14 August 1992; pp. 212–222. [Google Scholar]
  12. El-darymli, K.; Mcguire, P.; Gill, E.; Power, D.; Moloney, C. Understanding the significance of radiometric calibration for synthetic aperture radar imagery. In Proceedings of the IEEE 27th Canadian Conference on Electrical and Computer Engineering (CCECE), Toronto, ON, Canada, 4–7 May 2014; p. 6. [Google Scholar] [CrossRef]
  13. ESA (European Space Agency). 2014. Available online: http://earth.eo.esa.int/polsarpro/ (accessed on 1 August 2016).
  14. Thakur, P.K. SAR data processing to extract backscatter response from various features. In Proceedings of the Symposium Tutorials on Polarimetric SAR Data Processing and Applications, International Society for Photogrametry and Remote Sensing (ISPRS), Hyderabad, India, 9–12 December 2014. [Google Scholar]
  15. ASF (Alaska Satellite Facility). MapReady User Manual Remote Sensing Tool Kit, Engineering Group. 2015. Available online: https://media.asf.alaska.edu/uploads/pdf/mapready_manual_3.1.22.pdf (accessed on 13 November 2017).
  16. Frost, V.S.; Stiles, J.A.; Shanmugan, K.S.; Holtzman, J.C. A model for radar images and its application to adaptive digital filtering of multiplicative noise. IEEE Trans. Pattern Anal. Mach. Intell. 1982, 4, 157–166. [Google Scholar] [CrossRef] [PubMed]
  17. Fiscella, B.; Giancaspro, A.; Nirchio, F.; Pavese, P.; Trivero, P. Oil spill monitoring in the Mediterranean Sea using ERS SAR data. In Proceedings of the Envisat Symposium, ESA, Göteborg, Sweden, 16–20 October 2010; p. 9. [Google Scholar]
  18. Singha, S.; Bellerby, T.J.; Trieschmann, O. Satellite Oil Spill Detection Using Artificial Neural Networks. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 9. [Google Scholar] [CrossRef]
  19. Calabresi, G.; Del Frate, F.; Lichtenegger, I.; Petrocchi, A.; Trivero, P. Neural networks for the oil spill detection using ERS–SAR data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS ’99), Hamburg, Germany, 28 June–2 July 1999; pp. 215–217. [Google Scholar] [CrossRef]
  20. Solberg, A.H.S.; Storvik, G.; Solberg, R.; Volden, E. Automatic detection of oil spills in ERS SAR images. IEEE Trans. Geosci. Remote Sens. 1999, 37, 1916–1924. [Google Scholar] [CrossRef]
  21. Bentz, C.M. Reconhecimento Automático de Eventos Ambientais Costeiros e Oceânicos em Imagens de Radares Orbitais. Ph.D. Dissertation, Universidade Federal do Rio de Janeiro (UFRJ), COPPE, Rio de Janeiro, Brazil, 2006; p. 115. [Google Scholar]
  22. Pisano, A. Development of Oil Spill Detection Techniques for Satellite Optical Sensors and Their Application to Monitor Oil Spill Discharge in the Mediterranean Sea. Ph.D. Dissertation, Università di Bologna, Bologna, Italy, 2011; p. 146. [Google Scholar]
  23. McgarigaL, K.; Marks, B.J. FRAGSTATS: Spatial Pattern Analysis Program for Quantifying Landscape Structure. USDA For. Serv. Gen. Tech. Rep. PNW-351. 1994, p. 134. Available online: http://www.umass.edu/landeco/pubs/mcgarigal.marks.1995.pdf (accessed on 13 November 2017).
  24. Solberg, A.H.S.; Volden, E. Incorporation of prior knowledge in automatic classification of oil spills in ERS SAR images. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS ’97), Singapore, 3–8 August 1997; pp. 157–159. [Google Scholar]
  25. Miranda, F.P. Reconnaissance Geologic Mapping of a Heavily-Forested Shield Area (Guiana Shield, Northwestern Brazil). Ph.D. Dissertation, University of Nevada, Reno, NV, USA, 1990; p. 176. [Google Scholar]
  26. Fiscella, B.; Giancaspro, A.; Nirchio, F.; Pavese, P.; Trivero, P. Oil spill detection using marine SAR images. Int. J. Remote Sens. 2000, 21, 3561–3566. [Google Scholar] [CrossRef]
  27. Sneath, P.H.A.; SokaL, R.R. Numerical Taxonomy—The principles and Practice of Numerical Classification; W.H. Freeman and Company: San Francisco, CA, USA, 2017; p. 573. ISBN 0-7167-0697-0. [Google Scholar]
  28. Lane, D.M.; Scott, D.; Hebl, M.; Guerra, R.; Osherson, D.; Ziemer, H. Introduction to Statistics. 2015. Available online: http://onlinestatbook.com/Online_Statistics_Education.pdf (accessed on 13 November 2017).
  29. Moita Neto, J.M.; Moita, G.C. Uma introdução à análise exploratória de dados multivariados. Quím. Nova 1998, 21, 467–469. [Google Scholar] [CrossRef]
  30. Legendre, P.; Legendre, L. Developments in Environmental Modelling. In Numerical Ecology, Third ed.; Elsevier Science B.V.: Amsterdam, The Netherlands, 2012; p. 990. ISBN 978-0444538680. [Google Scholar]
  31. Valentin, J.L. Ecologia Numérica–Uma Introdução à Análise Multivariada de Dados Ecológicos, Second ed.; Editora Interciência: Rio de Janeiro, Brazil, 2012; 153p, ISBN 978-85-7193-230-2. [Google Scholar]
  32. Milligan, G.W.; Cooper, M.C. A study of standardization of variables in cluster analysis. J. Classif. 1988, 5, 181–204. [Google Scholar] [CrossRef]
  33. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: PAleontological STatistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 9. Available online: http://palaeo-electronica.org/2001_1/past/issue1_01.htm (accessed on 13 November 2017).
  34. Kelley, L.A.; Gardner, S.P.; Sutcliffe, M.J. An automated approach for clustering an ensemble of NMR-derived protein structures into conformationally related subfamilies. Protein Eng. Des. Sel. 1996, 9, 1063–1065. [Google Scholar] [CrossRef]
  35. Farris, J.S. On the cophenetic correlation coefficient. Syst. Biol. 1969, 18, 279–285. [Google Scholar] [CrossRef]
  36. Saraçli, S.; Dogan, N.; Dogan, I. Comparison of hierarchical cluster analysis methods by cophenetic correlation. J. Inequal. Appl. 2013, 2013, 203. [Google Scholar] [CrossRef]
  37. Sokal, R.R.; Rohlf, F.J. The Comparison of dendrograms by objective methods. Taxon 1962, 11, 33–40. [Google Scholar] [CrossRef]
  38. NCSS (Number Cruncher Statistical System). Hierarchical Clustering and Dendrograms. NCSS Statistical Software. 2015. Chapter 445. p. 15. Available online: http://ncss.wpengine.netdna-cdn.com/wp-content/themes/ncss/pdf/Procedures/NCSS/Hierarchical_Clustering-Dendrograms.pdf (accessed on 13 November 2017).
  39. Hall, M.A.; Smith, L.A. Feature subset selection: a correlation based filter approach. In Proceedings of the Fourth International Conference on Neural Information and Intelligent Information Systems, Dunedin, New Zealand, 24–28 November 1997; pp. 855–858. [Google Scholar]
  40. Hall, M.A.; Smith, L.A. Feature selection for machine learning: comparing a correlation-based filter approach to the wrapper. In Proceedings of the 12th International FLAIRS Conference, Orlando, FL, USA, 1–5 May 1999. [Google Scholar]
  41. Hall, M.A. Correlation-Based Feature Selection for Machine Learning. Ph.D. Dissertation, The University of Waikato, Department of Computer Science, Hamilton, New Zealand, 1999; p. 178. [Google Scholar]
  42. Bouckaert, R.R.; Frank, E.; Hall, M.; Kirkby, R.; Reutemann, P.; Seewald, A.; Scuse, D. WEKA Manual for Version 3–6-0; The University of Waikato: Hamilton, New Zealand, 2008; p. 212. Available online: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.153.9743 (accessed on 13 November 2017).
  43. Aha, D.W.; Bankert, R.L. Feature selection for case-based classification of cloud types: An empirical comparison. In Case-Based Reasoning: Papers From the 1994 Workshop (Technical Report WS-94–01); Advancement of Artificial Intelligence (AAAI): Menlo Park, CA, USA, 1994; pp. 106–112. Available online: https://www.aaai.org/Papers/Workshops/1994/WS-94-01/WS94-01-019.pdf (accessed on 13 November 2017).
  44. Tetko, I.; Baskin, I.; Varnek, A. Tutorial on Machine Learning, Part 2: Descriptor Selection Bias. 2008. Available online: http://infochim.u-strasbg.fr/CS3/program/Tutorials/Tutorial2b.pdf (accessed on 13 November 2017).
  45. Wilt, C.; Thayer, J.; Ruml, W.A. Comparison of Greedy Search Algorithms. In Proceedings of the Third Annual Symposium on Combinatorial Search (SoCS-10), Atlanta, GA, USA, 8–10 July 2010; pp. 129–136. [Google Scholar]
  46. Thompson, A.A.; McLeod, I.H. The RADARSAT-2 SAR processor. Can. J. Remote Sens. 2004, 30, 336–344. [Google Scholar] [CrossRef]
  47. Peres-Neto, P.R.; Jackson, D.A.; Somers, K.M. Giving meaningful interpretation to ordination axes: Assessing loading significance in principal component analysis. Ecology 2003, 84, 2347–2363. [Google Scholar] [CrossRef]
  48. Peres-Neto, P.R.; Jackson, D.A.; Somers, K.M. How many principal components? Stopping rules for determining the number of non-trivial axes revisited. Comput. Stat. Data Anal. 2005, 49, 974–997. [Google Scholar] [CrossRef]
  49. Hair, J.F.; Anserson, R.E.; Tatham, R.L.; Black, W.C. Multivariate Data Analysis. In Análise Multivariada de Dados, Fifth ed.; Sant’Anna, A.S., Chaves Neto, A., Eds.; Bookman: Porto Alegre, RS, Brazil, 2005; ISBN 0-13-014406-7. [Google Scholar]
  50. Cattell, R.B. The Scree Test For The Number of Factors. Multivar. Behav. Res. 1966, 1, 245–276. [Google Scholar] [CrossRef] [PubMed]
  51. Jackson, D.A. Stopping Rules in Principal Components Analysis: A Comparison of Heuristical and Statistical Approaches. Ecology 1993, 74, 2204–2214. [Google Scholar] [CrossRef]
  52. Hammer, Ø. PAST: Multivariate Statistics. 2015. Available online: http://folk.uio.no/ohammer/past/multivar.html (accessed on 13 November 2017).
  53. Kaiser, H.F. A note on Guttman’s lower bound for the number of common factors. Br. J. Stat. Psychol. 1961, 14, 2. [Google Scholar] [CrossRef]
  54. Maimon, O.; Rokach, L. Data Mining and Knowledge Discovery Handbook, Second ed.; Springer Science: New York, NY, USA, 2010; p. 1285. [Google Scholar] [CrossRef]
  55. Carvalho, G.A.; Minnett, P.J.; Fleming, L.E.; Banzon, V.F.; Baringer, W. Satellite remote sensing of harmful algal blooms: A new multi-algorithm method for detecting the Florida Red Tide (Karenia brevis). Harmful Algae 2010, 9, 440–448. Available online: https://www.ncbi.nlm.nih.gov/pubmed/21037979 (accessed on 13 November 2017). [CrossRef] [PubMed]
  56. Carvalho, G.A.; Minnett, P.J.; Banzon, V.F.; Baringer, W.; Heil, C.A. Long-term evaluation of three satellite ocean color algorithms for identifying harmful algal blooms (Karenia brevis) along the west coast of Florida: A matchup assessment. Remote Sens. Environ. 2011, 115, 18. Available online: http://www.ncbi.nlm.nih.gov/pubmed/22180667 (accessed on 13 November 2017). [CrossRef] [PubMed]
  57. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 1991, 37, 35–46. [Google Scholar] [CrossRef]
  58. Nussmeier, N.A.; Miao, Y.; Roach, G.W.; Wolman, R.L.; Mora-Mangano, C.; Fox, M.; Szekely, A.; Tommasino, C.; Schwann, N.M.; Mangano, D.T. Predictive value of the National Institutes of Health Stroke Scale and the Mini-Mental State Examination for neurologic outcome after coronary artery bypass graft surgery. J. Thorac. Cardiovasc. Surg. 2010, 139, 901–912. [Google Scholar] [CrossRef] [PubMed]
  59. Lee, M.-J.; Yun, M.J.; Park, M.-S.; Cha, S.H.; Kim, M.-J.; Lee, J.D.; Kim, K.W. Paraaortic lymph node metastasis in patients with intraabdominal malignancies: CT vs PET. World J. Gastroenterol. 2009, 15, 4434–4438. [Google Scholar] [CrossRef] [PubMed]
  60. Carvalho, G.A. The Use of Satellite-Based Ocean Color Measurements for Detecting the Florida Red Tide (Karenia brevis). Ph.D. Thesis, University of Miami (UM/RSMAS/MPO), Miami, FL, USA, 2008; p. 156. Available online: http://scholarlyrepository.miami.edu/oa_theses/116/ (accessed on 13 November 2017).
  61. Alberg, A.J.; Park, J.W.; Hager, B.W.; Brock, M.V.; Diener-West, M. The use of “overall accuracy” to evaluate the validity of screening or diagnostic tests. J. Gen. Intern. Med. 2004, 19, 460–465. [Google Scholar] [CrossRef] [PubMed]
  62. Theriault, C.; Scheibling, R.; Hatcher, B.; Jones, W. Mapping the distribution of an invasive marine alga (Codium fragile spp. tomentosoides) in optically shallow coastal waters using the compact airborne spectrographic imager (CASI). Can. J. Remote Sens. 2006, 32, 315–329. [Google Scholar] [CrossRef]
  63. Espedal, H.A.; Johannessen, O.M. Cover: Detection of oil spills near offshore installations using synthetic aperture radar (SAR). Int. J. Remote Sens. 2000, 21, 2141–2144. [Google Scholar] [CrossRef]
  64. Espedal, H.A. Detection of oil spill and natural film in the marine environment by spaceborne synthetic aperture radar. Ph.D. Dissertation, Department of Physics, University of Bergen and Nansen Environmental and Remote Sensing Center (NERSC), Bergen, Norway, 1998; p. 200. [Google Scholar]
  65. Espedal, H. Detection of oil spill and natural film in the marine environment by spaceborne SAR. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS ’99), Hamburg, Germany, 28 June–2 July 1999; pp. 1478–1480. [Google Scholar] [CrossRef]
  66. Holt, B. SAR imaging of the ocean surface. In Synthetic Aperture Radar Marine User’s Manual; NOAA/NESDIS; Office of Research and Applications: Washington, DC, USA, 2004; Chapter 2; pp. 25–79. Available online: http://www.sarusersmanual.com (accessed on 13 November 2017).
  67. Klinkenberg, B. A review of methods used to determine the fractal dimension of linear features. Math. Geol. 1994, 26, 23–46. [Google Scholar] [CrossRef]
  68. Bevilacqua, L.; Barros, M.M.; Galeão, A.C.R.N. Geometry, dynamics and fractals. J. Br. Soc. Mech. Sci. Eng. 2008, 30, 11–21. [Google Scholar] [CrossRef]
  69. Ledesma, R.D.; Valero-Mora, P. Determining the number of factors to retain in EFA: an easy-to-use computer program for carrying out Parallel Analysis. Pract. Assess. Res Eval. 2007, 12, 11. [Google Scholar]
  70. Garcia-Pineda, O.; MacDonald, I.; Zimmer, B. Synthetic aperture radar image processing using the Supervised Textural-Neural Network Classification Algorithm. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS ’08), Boston, MA, USA, 8–11 July 2008; pp. IV:1265–IV:1268. [Google Scholar] [CrossRef]
  71. Garcia-Pineda, O.; MacDonald, I.; Zimmer, B.; Shedd, B.; Roberts, H. Remote-sensing evaluation of geophysical anomaly sites in the outer continental slope, northern Gulf of Mexico. Deep Sea Res. Part II: Top. Stud. Oceanogr. 2010, 57, 1859–1869. [Google Scholar] [CrossRef]
  72. Garcia-Pineda, O.; Zimmer, B.; Howard, M.; Pichel, W.; Li, X.; MacDonald, I.R. Using SAR images to delineate ocean oil slicks with a texture-classifying neural network algorithm (TCNNA). Can. J. Remote Sens. 2009, 35, 411–421. [Google Scholar] [CrossRef]
  73. Mityagina, M.; Lavrova, O. Satellite Survey of Inner Seas: Oil Pollution in the Black and Caspian Seas. Remote Sens. 2016, 8, 875. [Google Scholar] [CrossRef]
  74. Song, D.; Ding, Y.; Li, X.; Zhang, B.; Xu, M. Ocean Oil Spill Classification with RADARSAT-2 SAR Based on an Optimized Wavelet Neural Network. Remote Sens. 2017, 9, 799. [Google Scholar] [CrossRef]
  75. Migliaccio, M.; Nunziata, F.; Buono, A. SAR polarimetry for sea oil slick observation. Int. J. Remote Sen. 2015, 36, 3243–3273. [Google Scholar] [CrossRef]
Figure 1. Gulf of Mexico highlighting the Campeche Bay region. The rectangle illustrates the region from which most oil slicks (98%) have been observed [9,10]. Isobaths of 200 m, 1000 m, and 3000 m are also shown. Courtesy of Adriano Vasconcelos.
Figure 1. Gulf of Mexico highlighting the Campeche Bay region. The rectangle illustrates the region from which most oil slicks (98%) have been observed [9,10]. Isobaths of 200 m, 1000 m, and 3000 m are also shown. Courtesy of Adriano Vasconcelos.
Ijgi 06 00379 g001
Figure 2. Summary of 44 data sub-divisions: 11 original sets (shown in the middle along with the number of variables; see also Table 5-1) and 33 optimal subsets of selected variables (shown in both sides). These have undergone two Attribute Selection Methods: UPGMA (Unweighted Pair Group Method with Arithmetic Mean) and CFS (Correlation-Based Feature Selection). The UPGMA implementation employs two user-defined thresholds to form groups on the analysis of dendrograms: CCC (Cophenetic Correlation Coefficient) and 0.5 (fixed similarity value). The CFS automatically provides the selected variables.
Figure 2. Summary of 44 data sub-divisions: 11 original sets (shown in the middle along with the number of variables; see also Table 5-1) and 33 optimal subsets of selected variables (shown in both sides). These have undergone two Attribute Selection Methods: UPGMA (Unweighted Pair Group Method with Arithmetic Mean) and CFS (Correlation-Based Feature Selection). The UPGMA implementation employs two user-defined thresholds to form groups on the analysis of dendrograms: CCC (Cophenetic Correlation Coefficient) and 0.5 (fixed similarity value). The CFS automatically provides the selected variables.
Ijgi 06 00379 g002
Figure 3. Rooted-tree dendrograms (i.e., diagrams of relationships) of the UPGMA implementation. 1st row left: only Size Information (10). 1st row right: only Digital Number (DN) variables (35). Size Information (10) together with: sigma-naught (σo; 141—2nd row), beta-naught (βo; 141—3rd row), and gamma-naught (γo; 141—4th row). Bottom and top horizontal red lines respectively correspond to fixed similarity of 0.5 and Cophenetic Correlation Coefficient (CCC: 0.9143 (Size), 0.8020 (DN), 0.8262 (σo), 0.8031 (βo), and 0.8240 (γo)).
Figure 3. Rooted-tree dendrograms (i.e., diagrams of relationships) of the UPGMA implementation. 1st row left: only Size Information (10). 1st row right: only Digital Number (DN) variables (35). Size Information (10) together with: sigma-naught (σo; 141—2nd row), beta-naught (βo; 141—3rd row), and gamma-naught (γo; 141—4th row). Bottom and top horizontal red lines respectively correspond to fixed similarity of 0.5 and Cophenetic Correlation Coefficient (CCC: 0.9143 (Size), 0.8020 (DN), 0.8262 (σo), 0.8031 (βo), and 0.8240 (γo)).
Ijgi 06 00379 g003
Figure 4. Results of the PCA exploring only Size Information. From top to bottom: original set (Table 2: 10 variables), UPGMA-CCC (Table 6: 6 variables), UPGMA-Fixed (Table 6: 5 variables), and CFS (Table 6: 3 variables). Left panels: Scree Plots (Eigenvalue x Principal Components (PC)), where the dashed red line corresponds to the broken stick curve. Bootstrapping eigenvalue error bars are shown as black lines: 2.5% and 97.5%. Middle panels: Scatterplots of the first two PCs, in which red triangles correspond to oil spills and green circles to oil seeps. Right panels: Loadings Plots.
Figure 4. Results of the PCA exploring only Size Information. From top to bottom: original set (Table 2: 10 variables), UPGMA-CCC (Table 6: 6 variables), UPGMA-Fixed (Table 6: 5 variables), and CFS (Table 6: 3 variables). Left panels: Scree Plots (Eigenvalue x Principal Components (PC)), where the dashed red line corresponds to the broken stick curve. Bootstrapping eigenvalue error bars are shown as black lines: 2.5% and 97.5%. Middle panels: Scatterplots of the first two PCs, in which red triangles correspond to oil spills and green circles to oil seeps. Right panels: Loadings Plots.
Ijgi 06 00379 g004
Table 1. Radiometric-calibrated image products: twelve forms of SAR backscatter signature using sigma-naught (σo), beta-naught (βo), and gamma-naught (γo), and the incidence angle (INC.ang).
Table 1. Radiometric-calibrated image products: twelve forms of SAR backscatter signature using sigma-naught (σo), beta-naught (βo), and gamma-naught (γo), and the incidence angle (INC.ang).
Image ProductsFrost Filterσoβoγo
C1 in AmplitudeWithoutSIG.ampBET.ampGAM.amp
WithSIG.amp.FFBET.amp.FFGAM.amp.FF
C2 in dBWithoutSIG.dBBET.dBGAM.dB
WithSIG.dB.FFBET.dB.FFGAM.dB.FF
Incience angle (INC.ang)
C1 = ([{DN2} + B]/A). C2 = [(10 × cc) × Log10(C1)]. DN: Digital Number. A: Range-dependent gain that varies per pixel in range direction. B: Constant offset nominally set to zero for SGF products. cc: 2 for the amplitude (or 1 for the intensity) of the received radar beam.
Table 2. Oil slick descriptors (53): 1st Attribute Type (Contextual Information: 6), 2nd Attribute Type (SAR Scene Descriptors: 37), and 3rd Attribute Type (Size Information: 10).
Table 2. Oil slick descriptors (53): 1st Attribute Type (Contextual Information: 6), 2nd Attribute Type (SAR Scene Descriptors: 37), and 3rd Attribute Type (Size Information: 10).
1st Attribute Type: Contextual Information §
1CategorySpill or SeepOil Slick Type *
2ClassSee [8]
3cLATLatitude (° N)Spatial Location
4cLONGLongitude (° W)
5SARtimeOverpass TimeTemporal Location
6SARdateOverpass Date
2nd Attribute Type: SAR Scene Descriptor
1Bmode § *Beam mode
2–37 **INC.angIncidence angle of the radar beam
3rd Attribute Type: Size Information (Geometry, Shape, and Dimension)
1LENNumber of pixels inside the oil slick polygon.
2Area §km2
3Per §kmPerimeter
4AtoPkmArea to Per ratioArea/Per
5PtoAkm−1Per to Area ratioPer/Area
6PtoA.nor*Normalized PtoAPer/[(2 × (Pi × Area))^(1/2)]
7COMPLEX.ind*Complexity Index(Per2)/Area
8COMPACT.ind*Compact Index(4 × Pi × Area)/(Per2)
9SHAPE.indkm−1Shape Index[0.25 × Per]/[Area^(1/2)]
10FRACTAL.ind*Fractal Index[2 × ln(0.25 × Per)]/[ln(Area)]
§ Present in the Campeche Bay Oil Slick Satellite Database (i.e., CBOS-Data). * Dimensionless quantity. ** Includes the 36 statistical measures calculated for the 4th Attribute Type (Table 3).
Table 3. Oil slick descriptors (432): 4th Attribute Type (SAR Backscatter Signature). Basic statistical measures (36) used to describe the oil slicks’ characteristics. Separately calculated for the different forms (12) of σo, βo, and γo: shown in Table 1: SIG.amp, SIG.amp.FF, SIG.dB, SIG.dB.FF, BET.amp, BET.amp.FF, BET.dB, BET.dB.FF, GAM.amp, GAM.amp.FF, GAM.dB, GAM.dB.FF.
Table 3. Oil slick descriptors (432): 4th Attribute Type (SAR Backscatter Signature). Basic statistical measures (36) used to describe the oil slicks’ characteristics. Separately calculated for the different forms (12) of σo, βo, and γo: shown in Table 1: SIG.amp, SIG.amp.FF, SIG.dB, SIG.dB.FF, BET.amp, BET.amp.FF, BET.dB, BET.dB.FF, GAM.amp, GAM.amp.FF, GAM.dB, GAM.dB.FF.
4th Attribute Type: SAR Backscatter Signature36 × 12 = 432
Basic Statistical Measures12 × 12 = 144
11–12AVGAverageCentral Tendency (4)
213–24MEDMedian
325–36MODMode
437–48MDM *Mid-Mean
549–60STDStandard DeviationDispersion (6)
661–72COD **Coefficient of Dispersion
773–84VARVariance
885–96RNGTotal Range
997–108AAD ***Average Absolute Deviation
10109–120MAD ****Median Absolute Deviation
11121–132MIN §Minimum
12133–144MAXMaximum
COV = Coefficient of Variation24 × 12 = 288
13145–156COV.STD/AVG ¥1st combined COV set
(STD divided by Central Tendency)
14157–168COV.STD/MED
15169–180COV.STDMOD
16181–192COV.STD/MDM
17192–204COV.COD/AVG2nd combined COV set
(COD divided by Central Tendency)
18205–216COV.COD/MED
19217–228COV.COD/MOD
20229–240COV.COD/MDM
21241–252COV.VAR/AVG3rd combined COV set
(VAR divided by Central Tendency)
22253–264COV.VAR/MED
23265–276COV.VAR/MOD
24277–288COV.VAR/MDM
25289–300COV.RNG/AVG4th combined COV set
(RNG divided by Central Tendency)
26301–312COV.RNG/MED
27313–324COV.RNG/MOD
28325–336COV.RNG/MDM
29337–348COV.AAD/AVG5th combined COV set
(AAD divided by Central Tendency)
30349–360COV.AAD/MED
31361–372COV.AAD/MOD
32373–384COV.AAD/MDM
33385–396COV.MAD/AVG6th combined COV set
(MAD divided by Central Tendency)
34397–408COV.MAD/MED
35409–420COV.MAD/MOD
36421–432COV.MAD/MDM
* MDM: Mean of the interquartile range. ** COD: Subtract the 1st from the 3rd interquartile, divided by their sum. *** AAD: Mean of the absolute difference of each value minus the mean. **** MAD: Median of the absolute difference of each value minus the median. § Negative Values Scaling eliminates 9 MIN variables. ¥ Original COV ratio.
Table 4. Confusion Matrices utilized to assess the Discriminant Function accuracy. Adapted from [55,56]. Producer’s Accuracy (PA), Commission Error (CE), User’s Accuracy (UA), and Omission Error (OE) are given as frequency of occurrence (§) and in percentages (¥).
Table 4. Confusion Matrices utilized to assess the Discriminant Function accuracy. Adapted from [55,56]. Producer’s Accuracy (PA), Commission Error (CE), User’s Accuracy (UA), and Omission Error (OE) are given as frequency of occurrence (§) and in percentages (¥).
§Algorithm Outcome SeepsAlgorithm Outcome Spills §Algorithm Outcome SeepsAlgorithm Outcome Spills
Actual SeepsGOOD SeepsBAD SeepsAll Actual SeepsActual SeepsABA + B
Actual SpillsBAD SpillsGOOD SpillsAll Actual SpillsActual SpillsCDC + D
All Algorithm SeepsAll Algorithm SpillsAll Actual Slicks A+CB+DA + B + C + D
Overall Accuracy Overall Accuracy
¥Algorithm Outcome SeepsAlgorithm Outcome Spills ¥Algorithm Outcome SeepsAlgorithm Outcome Spills
Actual SeepsSensitivityActual Seeps100% Actual Seeps A   ×   100 A   +   B B   ×   100 A   +   B 100%
Actual SpillsFalse PositiveSpecificity100% Actual Spills C   ×   100 C   +   D D   ×   100 C   +   D 100%
¥Algorithm Outcome SeepsAlgorithm Outcome Spills ¥Algorithm Outcome SeepsAlgorithm Outcome Spills
Actual SeepsPositive Predictive ValueInvert Neg Predictive Value Actual Seeps A   ×   100 A   +   C B   ×   100 B   +   D
Actual SpillsInvert Pos Predictive ValueNegative Predictive Value Actual Spills C   ×   100 A   +   C D   ×   100 B   +   D
100%100% 100%100%
PA/CE Q1)How many known oil seep samples are correctly identified?
Sensitivity: A, also given in percentage by (A × 100)/[A + B].
PA/CE Q2)How many known oil spill samples are correctly identified?
Specificity: D, also given in percentage by (D × 100)/[C + D].
PA/CE Q3)How many known oil seep samples are misidentified?
False Negative cases: B, also given in percentage by (B × 100)/[A + B]. Coupled with Sensitivity.
PA/CE Q4)How many known oil spill samples are misidentified?
False Positive cases: C, also given in percentage by (C × 100)/[C + D]. Linked to Specificity.
UA/OE Q1)How many oil seeps identified by the Function are indeed known oil seeps?
Positive Predictive Value: A, also given in percentage by (A × 100)/[A + C].
UA/OE Q2)How many oil spills identified by the Function are indeed known oil spills?
Negative Predictive Value: D, also given in percentage by (D × 100)/[B + D].
UA/OE Q3)Of samples identified by the Function as oil seeps, how many are oil spills?
Inverse of the Positive Predictive Value: C, also given in percentage by (C × 100)/[A + C].
UA/OE Q4)Of samples identified by the Function as oil spills, how many are oil seeps?
Inverse of the Negative Predictive Value: B, also given in percentage by (B × 100)/[B + D].
Table 5. Data mining summary of the 11 original sets and the 33 optimal subsets of selected variables proposed on Figure 3, collectively referred to as to as 44 data sub-divisions. 1) Values of UPGMA-CCC and CFS-Merit. 2) Number of variables. 3) Number of selected Principal Components (PCs). See also Table 9 for full PCA results.
Table 5. Data mining summary of the 11 original sets and the 33 optimal subsets of selected variables proposed on Figure 3, collectively referred to as to as 44 data sub-divisions. 1) Values of UPGMA-CCC and CFS-Merit. 2) Number of variables. 3) Number of selected Principal Components (PCs). See also Table 9 for full PCA results.
11 Original Sets and
33 Optimal Subsets
1) Values of2) Number of Variables:3) Number of Selected PC’s:
UPGMA (CCC)CFS (Merit)Original SetsUPGMACFSOriginal SetsUPGMACFS
(CCC)(0.5)(CCC)(0.5)
1. Complete Exploration of the CBOS-MOD0.81751.000502594615 #102118§
2. ❖ With Size Information0.82270.1314332086
3. ❖ Only0.82530.103423††††367††††6
4. σo with Size Information0.82620.1271512612266737
5. βo with Size Information0.80310.1291512512156735
6. γo with Size Information0.82400.1261512612196735
7. Size Information only0.91430.112106532222 ¥
8. Only σo0.83190.099141207245524
9. Only βo0.81300.099141197215425
10. Only γo0.82600.097141207226524
11. DNs only0.80200.0593553112222
❖ All SAR Backscatter Signature (σo, βo, and γo; See also Table 1 and Table 3). σo Sigma-naught: SIG.amp, SIG.amp.FF, SIG.dB, and SIG.dB.FF. βo Beta-naught: BET.amp, BET.amp.FF, BET.dB, and BET.dB.FF. γo Gamma-naught: GAM.amp, GAM.amp.FF, GAM.dB, and GAM.dB.FF. † Same as σo with Size Information. †† Same as σo only. § PCA not performed because most categorical are dummy variables. # 13 are categorical contextual dummy variables (see [8]). ¥ Because the Discriminant Function requires at least 2 variables. DNs: Digital Numbers.
Table 6. Attribute Selection results. Top: Only Size Information (Table 2: 10). UPGMA implementation: Variables in bold (5) were selected using the fixed similarity value of 0.5. See also Figure 3 (1st row left: CCC = 0.9143). CFS implementation: Variables in bold (1) was also selected with the UPGMA implementation. Bottom: Only Digital Numbers (DNs) variables (35). UPGMA implementation: Variables in bold (3) were selected using the fixed similarity value of 0.5. See also Figure 3 (1st row right: CCC = 0.8020). CFS implementation: Variables in bold (3) are the same as those selected with the UPGMA implementation.
Table 6. Attribute Selection results. Top: Only Size Information (Table 2: 10). UPGMA implementation: Variables in bold (5) were selected using the fixed similarity value of 0.5. See also Figure 3 (1st row left: CCC = 0.9143). CFS implementation: Variables in bold (1) was also selected with the UPGMA implementation. Bottom: Only Digital Numbers (DNs) variables (35). UPGMA implementation: Variables in bold (3) were selected using the fixed similarity value of 0.5. See also Figure 3 (1st row right: CCC = 0.8020). CFS implementation: Variables in bold (3) are the same as those selected with the UPGMA implementation.
Only Size Information
UPGMA Selected Attributes
1LEN
2AtoP
3PtoA
4COMPACT.ind
5COMPLEX.ind
6FRACTAL.ind
CFS Selected Attributes
1LEN
2PER
3SHAPE
Only DNs
UPGMA Selected Attributes
1DN.AVG
2DN.COD
3DN.STD/MOD
4DN.COD/AVG
5DN.RNG/AVG
CFS Selected Attributes
1DN.AVG
2DN.MDM
3DN.STD
4DN.COD
5DN.RNG
6DN.COD/MOD
7DN.RNG/AVG
8DN.RNG/MED
9DN.RNG/MOD
10DN.RNG/MDM
11DN.MAD/MED
Table 7. Pool of attributes selected by the UPGMA and CFS implementations considering each SAR backscatter coefficient, i.e., sigma-naught (σo), beta-naught (βo), and gamma-naught (γo), with and without the Size Information (*). Double-checked variables (**) have also been selected with the UPGMA fixed similarity value of 0.5. Figure 3 portrays the UPGMA dendrograms.
Table 7. Pool of attributes selected by the UPGMA and CFS implementations considering each SAR backscatter coefficient, i.e., sigma-naught (σo), beta-naught (βo), and gamma-naught (γo), with and without the Size Information (*). Double-checked variables (**) have also been selected with the UPGMA fixed similarity value of 0.5. Figure 3 portrays the UPGMA dendrograms.
Selection MethodUPGMACFSCFS
Size InformationWithout §WithWithout
Attribute Inventoryσoβoγoσoβoγoσoβoγo
1LEN ***
2AREA ***
3PER ***
4AtoP
5PtoA * *
6PtoAnor *
7COMPACT.ind ***
8COMPLEX.ind * *
9SHAPE.ind ***
10FRACTAL.ind ***
1amp.AVG******
2amp.MOD * * *
3amp.STD******
4amp.COD************
5amp.RNG ******
6amp.MAX * **
7amp.STD/MOD *
8amp.STD/MDM *
9amp.RNG/MED * *
10amp.RNG/MOD *
11amp.RNG/MDM * *
12amp.MAD/AVG *
13amp.FF.AVG***
14amp.FF.COD******
15amp.FF.RNG * **
16amp.FF.VAR/AVG***
17amp.FF.RNG/AVG * * ***
18amp.FF.RNG/MDM *
19amp.FF.RNG/MOD *
20dB.AVG******
21dB.MED * ***
22dB.MOD***
23dB.STD***** ***
24dB.COD******* ****
25dB.MAD**** ****
26dB.STD/AVG * ***
27dB.STD/MED ***
28dB.STD/MOD***
29dB.VAR/AVG*** *****
30dB.VAR/MED * ***
31dB.VAR/MOD * ***
32dB.RNG/AVG***
33dB.RNG/MED * ** *
34dB.RNG/MOD*********
35dB.FF.AVG******
36dB.FF.COD* *
37dB.FF.RNG * ****
38dB.FF.MIN* *
39dB.FF.MAX***
40dB.FF.MAD ******
41dB.FF.COD/AVG * *
42dB.FF.COD/MED *
43dB.FF.VAR/AVG *
44dB.FF.RNG/AVG* *
45dB.FF.RNG/MOD ***
46dB.FF.MAD/AVG * *
47dB.FF.MAD/MOD *
57Total (**)20 (7)19 (7)20 (7)261519242122
§ The UPGMA with Size Information is the same but with the attributes of Table 6.
Table 8. Discriminant Function results: Hotelling’s t2 values and Overall Accuracy of the 44 data sub-divisions proposed on Figure 2.
Table 8. Discriminant Function results: Hotelling’s t2 values and Overall Accuracy of the 44 data sub-divisions proposed on Figure 2.
Hotelling’s t2 ValuesOriginal SetsUPGMA (CCC)UPGMA (0.5)CFS
Complete Exploration17,056.089,403.094,086.0§
All SAR Backscatter Signature * with Size Information1065.51044.5
All SAR Backscatter Signature * only598.9††††632.4
σowith Size Information891.01123.5769.81065.8
βowith Size Information908.71099.5722.31106.2
γowith Size Information863.11115.4753.41028.4
Size Information only1034.21064.71086.41085.9
σoOnly594.1589.5447.7640.9
βoOnly584.1464.5435.7612.6
γoOnly579.2574.7432.2606.7
Digital Numbers (DNs) only138.0163.0138.4210.5
Overall AccuracyOriginal SetsUPGMA (CCC)UPGMA (0.5)CFS
Complete Exploration90.04%99.98%99.96%§
All SAR Backscatter Signature * with Size Information68.61%68.71%
All SAR Backscatter Signature * only63.63%††††63.93%
σowith Size Information67.64%69.51%65.46%69.18%
βowith Size Information68.02%69.57%65.58%69.55%
γowith Size Information67.47%69.32%65.79%68.59%
Size Information only69.59%70.00%70.02%70.22%
σoOnly63.28%63.16%63.26%64.12%
βoOnly63.53%63.20%63.02%64.46%
γoOnly63.16%63.02%63.24%64.18%
Digital Numbers (DNs) only57.65%57.89%56.75%59.48%
σo Sigma-naught: SIG.amp, SIG.amp.FF, SIG.dB, and SIG.dB.FF. βo Beta-naught: BET.amp, BET.amp.FF, BET.dB, and BET.dB.FF. γo Gamma-naught: GAM.amp, GAM.amp.FF, GAM.dB, and GAM.dB.FF. * σo, βo, and γo. See also Table 1 and Table 3. † Same as σo with Size Information. †† Same as σo only. § Discriminant Function not performed because most categorical are dummy variables.
Table 9. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): all SAR backscatter signature variables (Table 3: σo, βo, and γo) analyzed together with (top) and without (bottom) the Size Information (Table 2). Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). See Table 4 for support.
Table 9. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): all SAR backscatter signature variables (Table 3: σo, βo, and γo) analyzed together with (top) and without (bottom) the Size Information (Table 2). Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). See Table 4 for support.
All SAR Backsc. Sig. with Size Info. (7 PCs)All SAR Backsc. Sig. with Size Info. (6 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep13956262021 Seep13246972021
Spill91719782895 Spill84120542895
Total231226044916Over. Acc.Total216527514916Over. Acc.
69% 69%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep69%31%100% Seep66%34%100%
Spill32%68%100% Spill29%71%100%
CCCSeepSpill CFSSeepSpill
Seep60%24% Seep61%25%
Spill40%76% Spill39%75%
Total100%100% Total100%100%
All SAR Backscatter Signature only (5 PCs)All SAR Backscatter Signature only (6 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep14305912021 Seep13756462021
Spill119716982895 Spill112717682895
Total262722894916Over. Acc.Total250224144916Over. Acc.
64% 64%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep71%29%100% Seep68%32%100%
Spill41%59%100% Spill39%61%100%
CCCSeepSpill CFSSeepSpill
Seep54%26% Seep55%27%
Spill46%74% Spill45%73%
Total100%100% Total100%100%
Table 10. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): only sigma-naught (σo) variables analyzed together with (top) and without (bottom) the Size Information (Table 2). Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). Carvalho [8] presents similar Confusion Matrices for the separate analysis of beta-naught (βo) and gamma-naught (γo). See Table 4 for support.
Table 10. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): only sigma-naught (σo) variables analyzed together with (top) and without (bottom) the Size Information (Table 2). Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). Carvalho [8] presents similar Confusion Matrices for the separate analysis of beta-naught (βo) and gamma-naught (γo). See Table 4 for support.
Sigma with Size Information (7 PCs)Sigma with Size Information (7 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep13506712021 Seep13506712021
Spill82820672895 Spill84420512895
Total217827384916Over. Acc.Total219427224916Over. Acc.
70% 69%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep67%33%100% Seep67%33%100%
Spill29%71%100% Spill29%71%100%
CCCSeepSpill CFSSeepSpill
Seep62%25% Seep62%25%
Spill38%75% Spill38%75%
Total100%100% Total100%100%
Sigma only (5 PCs)Sigma only (4 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep14196022021 Seep13976242021
Spill120916862895 Spill114017552895
Total262822884916Over. Acc.Total253723794916Over. Acc.
63% 64%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep70%30%100% Seep69%31%100%
Spill42%58%100% Spill39%61%100%
CCCSeepSpill CFSSeepSpill
Seep54%26% Seep55%26%
Spill46%74% Spill45%74%
Total100%100% Total100%100%
Table 11. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs). Top: Size Information only (Table 2). Bottom: Digital Numbers (DNs) only. Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). See Table 4 for support.
Table 11. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs). Top: Size Information only (Table 2). Bottom: Digital Numbers (DNs) only. Left: CCC (Cophenetic Correlation Coefficient). Right: CFS (Correlation-Based Feature Selection). See Table 4 for support.
Size Information only (2 PCs)Size Information only (2 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep13147072021 Seep13376842021
Spill76821272895 Spill78021152895
Total208228344916Over. Acc.Total211727994916Over. Acc.
70% 70%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep65%35%100% Seep66%34%100%
Spill27%73%100% Spill27%74%100%
CCCSeepSpill CFSSeepSpill
Seep63%25% Seep63%24%
Spill37%75% Spill37%76%
Total100%100% Total100%100%
Digital Numbers only (2 PCs)Digital Numbers only (2 PCs)
CCCSeepSpillTotal CFSSeepSpillTotal
Seep10749472021 Seep10739482021
Spill112317722895 Spill104418512895
Total219727194916Over. Acc.Total211727994916Over. Acc.
58% 59%
CCCSeepSpillTotal CFSSeepSpillTotal
Seep53%47%100% Seep53%47%100%
Spill39%61%100% Spill36%64%100%
CCCSeepSpill CFSSeepSpill
Seep49%35% Seep51%34%
Spill51%65% Spill49%66%
Total100%100% Total100%100%
Table 12. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): sole and simple use of two specific oil slicks Size Information, i.e., area and perimeter. See also Table 4 for further information about Confusion Matrixes. See Table 4 for support.
Table 12. Confusion Matrices of the Discriminant Functions using selected Principal Components (PCs): sole and simple use of two specific oil slicks Size Information, i.e., area and perimeter. See also Table 4 for further information about Confusion Matrixes. See Table 4 for support.
Area and Perimeter (2 PCs)
SizeSeepSpillTotal
Seep13147072021
Spill79421022895
Total210828084916Over. Acc.
69%
SizeSeepSpillTotal
Seep65%35%100%
Spill27%73%100%
SizeSeepSpill
Seep62%25%
Spill38%75%
Total100%100%
PA/CE Q1)65%Sensitivity
PA/CE Q2)73%Specificity
PA/CE Q3)35%False Negative
PA/CE Q4)27%False Positive
UA/OE Q1)62%Positive Predictive Value
UA/OE Q2)75%Negative Predictive Value
UA/OE Q3)38%Inv. of the Pos. Pred. Val.
UA/OE Q4)25%Inv. of the Neg. Pred. Val.

Share and Cite

MDPI and ACS Style

Carvalho, G.D.A.; Minnett, P.J.; De Miranda, F.P.; Landau, L.; Paes, E.T. Exploratory Data Analysis of Synthetic Aperture Radar (SAR) Measurements to Distinguish the Sea Surface Expressions of Naturally-Occurring Oil Seeps from Human-Related Oil Spills in Campeche Bay (Gulf of Mexico). ISPRS Int. J. Geo-Inf. 2017, 6, 379. https://doi.org/10.3390/ijgi6120379

AMA Style

Carvalho GDA, Minnett PJ, De Miranda FP, Landau L, Paes ET. Exploratory Data Analysis of Synthetic Aperture Radar (SAR) Measurements to Distinguish the Sea Surface Expressions of Naturally-Occurring Oil Seeps from Human-Related Oil Spills in Campeche Bay (Gulf of Mexico). ISPRS International Journal of Geo-Information. 2017; 6(12):379. https://doi.org/10.3390/ijgi6120379

Chicago/Turabian Style

Carvalho, Gustavo De Araújo, Peter J. Minnett, Fernando Pellon De Miranda, Luiz Landau, and Eduardo Tavares Paes. 2017. "Exploratory Data Analysis of Synthetic Aperture Radar (SAR) Measurements to Distinguish the Sea Surface Expressions of Naturally-Occurring Oil Seeps from Human-Related Oil Spills in Campeche Bay (Gulf of Mexico)" ISPRS International Journal of Geo-Information 6, no. 12: 379. https://doi.org/10.3390/ijgi6120379

APA Style

Carvalho, G. D. A., Minnett, P. J., De Miranda, F. P., Landau, L., & Paes, E. T. (2017). Exploratory Data Analysis of Synthetic Aperture Radar (SAR) Measurements to Distinguish the Sea Surface Expressions of Naturally-Occurring Oil Seeps from Human-Related Oil Spills in Campeche Bay (Gulf of Mexico). ISPRS International Journal of Geo-Information, 6(12), 379. https://doi.org/10.3390/ijgi6120379

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