Next Article in Journal
Machine Learning Approaches for Detecting Tropical Cyclone Formation Using Satellite Data
Next Article in Special Issue
Large-Area, High Spatial Resolution Land Cover Mapping Using Random Forests, GEOBIA, and NAIP Orthophotography: Findings and Recommendations
Previous Article in Journal
Analysis of Vegetation Red Edge with Different Illuminated/Shaded Canopy Proportions and to Construct Normalized Difference Canopy Shadow Index
Previous Article in Special Issue
GEOBIA 2016: Advances in Object-Based Image Analysis—Linking with Computer Vision and Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Mapping and Monitoring of Whiting in the Semi-Enclosed Gulf Using Moderate Resolution Imaging Spectroradiometer (MODIS) Time Series Images and a Generic Ensemble Tree-Based Model

by
Abdallah Shanableh
1,2,
Rami Al-Ruzouq
1,2,*,
Mohamed Barakat A. Gibril
2,
Cristina Flesia
3 and
Saeed AL-Mansoori
4
1
Department of Civil and Environmental Engineering, University of Sharjah, Sharjah 27272, UAE
2
Research Institute of Sciences and Engineering, University of Sharjah, Sharjah 27272, UAE
3
Department of Earth and Environmental Sciences, University of Milano Bicocca, Piazza Della Scienza 4, 20126 Milano, Italy
4
Applications Development and Analysis Section (ADAS), Mohammed Bin Rashid Space Centre (MBRSC), Dubai 211833, UAE
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(10), 1193; https://doi.org/10.3390/rs11101193
Submission received: 31 March 2019 / Revised: 7 May 2019 / Accepted: 13 May 2019 / Published: 20 May 2019
(This article belongs to the Special Issue Geographic Object-Based Image Analysis (GEOBIA))

Abstract

:
Whiting events in seas and lakes are a natural phenomenon caused by suspended calcium carbonate (CaCO3) particles. The Arabian Gulf, which is a semi-enclosed sea, is prone to extensive whiting that covers tens of thousands of square kilometres. Despite the extent and frequency of whiting events in the Gulf, studies documenting the whiting phenomenon are lacking. Therefore, the primary objective of this study was to detect, map and document the spatial and temporal distributions of whiting events in the Gulf using daily images acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA’s Terra and Aqua satellites from 2002 to 2018. A method integrating a geographic object-based image analysis, the correlation-based feature selection technique (CFS), the adaptive boosting decision tree (AdaBoost DT) and the rule-based classification were used in the study to detect, quantify and assess whiting events in the Gulf from the MODIS data. Firstly, a multiresolution segmentation was optimised using unsupervised quality measures. Secondly, a set of spectral bands and indices were investigated using the CFS to select the most relevant feature(s). Thirdly, a generic AdaBoost DT model and a rule-based classification were adopted to classify the MODIS time series data. Finally, the developed classification model was compared with various tree-based classifiers such as random forest, a single DT and gradient boosted DT. Results showed that both the combination of the mean of the green spectral band and the normalised difference index between the green and blue bands (NDGB), or the combination of the NDGB and the colour index for estimating the concentrations of calcium carbonates (CI) of the image objects, were the most significant features for detecting whiting. Moreover, the generic AdaBoost DT classification model outperformed the other tested tree-based classifiers with an overall accuracy of 97.86% and a kappa coefficient of 0.97. The whiting events during the study period (2002–2018) occurred exclusively during the winter season (November to March) and mostly in February. Geographically, the whiting events covered areas ranging from 12,000 km2 to 60,000 km2 and were mainly located along the southwest coast of the Gulf. The duration of most whiting events was 2 to 6 days, with some events extending as long as 8 to 11 days. The study documented the spatiotemporal distribution of whiting events in the Gulf from 2002 to 2018 and presented an effective tool for detecting and motoring whiting events.

Graphical Abstract

1. Introduction

1.1. Background

Whiting is a short-lived phenomenon of milky parcels of water or bright in-water features, which has been reported globally in lacustrine, marine and freshwater environments and in semi-enclosed areas [1,2,3,4,5,6,7]. The ephemeral patches in whiting are turbid water with high levels of suspended fine-grained calcium carbonate mineral particles [8,9,10,11]. Whiting events last from days to weeks and can be visualised via satellite images as extended milky-white water varying from a few meters to square kilometres long [1,12].
Although whiting has been studied for decades to determine valid explanations for its occurrences and causes, the event remains controversial [10]. Various assumptions have been made to explain the causes of whiting. These assumptions include (1) the resuspension of fine-grained sediments caused by fish activities, microturbulent bursts and wind [8,12,13,14,15]; (2) bio-induced precipitation from the removal of CO2 by photosynthesis [11,16,17,18,19] and (3) abiotic precipitation initiated by fluctuations in water temperature and ion activities related to climate change [8,20]. Carbonates, which are produced by the physical and biological disintegration of animal and algal bioclasts, blooms of microscopic algae during photosynthesis and abiotic precipitation or calcification of suspended picoplankton and organic matter, may be the possible sources of suspended carbonate minerals, such as aragonite and high and low magnesium calcite [8].
The occurrence of whiting has prompted researchers to establish a scientific description of the phenomena and to determine its association with climate change and oil deposition [21,22,23]. Whiting events are ephemeral, thus sample collection using traditional field measurements can be challenging, time consuming and costly. Remote sensing with ocean colour satellite instruments provides a set of high temporal-resolution data with various scales and records of satellite images and derivatives. This method enables the spatiotemporal mapping of whiting events; however, limited studies have adopted remote sensing technology to study these occurrences in marine environments [10]. The majority of the studies used satellite data to map whiting in the Bahama Banks [12,24], in the Ten Thousand Islands in southwest Florida [10,25], in the Great Lakes in North America [26] and in the Feldberg Lake District, Klocksin Lake Chain and Rheinsberg Lake regions in Germany [27]. Whiting events based on general properties are varyingly recognised with remote sensing data from visual identification and manual delineation [10,24,28,29,30,31]. Dierssen et al. [12] studied the spectral behaviour of whiting to identify the spatial extent of whiting patches in the Bahama Banks. Considering the shallowness of the southwest Florida coast and the spectral similarity between whiting occurrences and bright shallow bottom sediments, Long et al. [25] delineated whiting patches manually by relying on visual inspection and spatial contrast. Recently, Long et al. [10] used contrast enhancement and floating algae index (FAI) images to differentiate in-water whiting features from clouds to map the spatiotemporal variability of the southwest Florida whiting events from 2003 to 2015. Other studies focused on detecting and estimating the concentration of particulate inorganic carbon (PIC), or calcium carbonate particles, in the surface layer of the water column from the water-leaving radiance and reflectance differences concept through the computation of chlorophyll-a concentration [32,33,34].
The whiting phenomenon in the Gulf was initially reported in the 1960s. In 1962, Wells and Illing [35] observed whiting in numerous places in the Gulf such as in the eastern part of the Qatar Peninsula towards the coast of Abu Dhabi and off the coast of Saudi Arabia between Ras Tanura and Ras Safaniya. Although more than 50 years have passed since whiting events were initially reported in the Gulf, this phenomenon remains unclear [36,37]. To the best of the authors’ knowledge, limited effort has been exerted to map whiting in the Gulf using remote sensing techniques. The objectives of the present study were: (1) to adopt the correlation-based feature selection (CFS) to identify the most significant features for whiting extraction, (2) to develop a semi-automated framework to detect the whiting coverage in the Gulf from the Moderate Resolution Imaging Spectroradiometer (MODIS) images, using adaptive boosting (AdaBoost) and rule-based classification approach, (3) to compare and assess the performance of various applied tree-based machine learning methods, namely, the single decision tree (DT), random forest and the gradient boosted decision tree (GBDT)and (4) to document the frequency, duration, seasonality, spatial coverage and distribution of whiting occurrences in the Gulf between 2002 and 2018 using satellite observations.

1.2. Study Region

The Arabian Gulf is located in the Middle East and surrounded by the coasts of eight countries, namely, United Arab Emirates (UAE), Saudi Arabia, Oman, Kuwait, Bahrain, Iraq, Qatar and Iran (Figure 1). The Gulf is a semi-enclosed marginal sea positioned in a subtropical hyperarid region (between the latitudes of 24° to 30°N and the latitudes of 48° to 57 °E), with an average annual rainfall of less than 5 cm in the coastal areas [33,34]. The Gulf is nearly 990 km long and 56–370 km wide. It has an average depth of 36 m and a maximum depth of nearly 100 m and occupies a surface area of 239 km2. The deepest region of the Gulf (more than 40 m deep) is near the Iranian coast and continues into the Strait of Hormuz. Meanwhile, the shallowest regions (less than 20 m deep) are located along the coasts of the UAE, Qatar, Bahrain and around the head of the Gulf. Owing to its shallowness, sea surface temperatures fluctuate significantly and the Gulf is considered as the hottest sea in the world during summer [35,38]. The seawater temperature ranges from less than 20 °C in winter to over 34 °C in summer. As a result of high evaporation rates during the hot and long summers in the region and the lack of precipitation, water in the Gulf is characterised by high salinity greater than 39 psu [33,39,40]. The Gulf is subjected to strong winds and often associated with dust storms, with the most extreme occurring in summer and late spring and moderate dust storms occurring in winter.

2. Methodology

2.1. Overview

Daily satellite data from MODIS for the period of 2002 to 2018 were obtained and visually inspected for whiting. The dates of whiting events were then classified in accordance with the month of the occurrence and the number of consecutive days during which whiting persisted. The data were then used to analyse frequency, seasonality and duration of whiting events in the Gulf.
The identification, mapping and classification of whiting in the Gulf required the extensive analysis of the acquired MODIS data. The analysis framework, which is illustrated in Figure 2, is summarised as follows: (1) preprocessing of the MODIS satellite data, (2) multiresolution image segmentation and parameter optimisation, (3) selection of the significant attributes using CFS, (4) classification of the daily MODIS time series images using AdaBoost DT and rule-based classification and (5) identification of the whiting spatiotemporal pattern (time series frequency, duration and seasonality, spatial coverage and distribution) in the Gulf. The aforementioned steps are further discussed in the following subsections.

2.2. MODIS Datasets for Whiting Exploration

The exploration of whiting events in the Gulf as a short-lived, repetitive phenomenon by field-based studies is challenging. Therefore, high temporal resolution satellite images offered excellent sources of data for mapping and monitoring whiting in the Gulf. The daily high temporal-resolution Terra/Aqua MODIS surface reflectance products (MOD09GA and MYD09GA) with a coarse spatial resolution of 500 m were used in this study. The products, which covered the entire Gulf in one scene for the period of 2002–2018, were downloaded from NASA’s Earthdata website (https://search.earthdata.nasa.gov/search). These products were atmospherically corrected by the MODIS Land Science Team for aerosols, thin cirrus clouds and gases.
Figure 3 illustrates the initiation and disappearance of a whiting event from February to March 2003. The consecutive MODIS images show the dramatic changes during the 6 consecutive days of observation. High-resolution satellite images, such as from the Landsat, with a 30 m spatial resolution and 16-day temporal resolutions, or from the Sentinel-2, with a 5-day revisit time and 10 m resolution, may be insufficient to capture such short-term whiting events given their limited temporal resolution.
The MODIS products (MOD09GA and MYD09GA) contained seven spectral bands with a spatial resolution of 0.5 km. Their reflectance bands were band 1 (red: 0.620–0.670 µm), band 2 (near infrared (NIR): 0.841–0.876 µm), band 3 (blue: 0.459–0.479 µm), band 4 (green: 0.545–0.565 µm), band 5 (shortwave infrared (SWIR1): 230–1.250 µm), band 6 (shortwave infrared (SWIR2): 1.628–1.652 µm) and band 7 (shortwave infrared (SWIR3): 2.105–2.155 µm). The data were reprojected from the sinusoidal coordinate system to the World Geodetic System 1984 and then converted from heirarchical data format to georeferenced tagged image file format. The images were then cropped on the basis of the boundaries of the Gulf.
The whiting events appeared in the MODIS satellite images as turbid and milky features. The spectral response of various surface water features at different locations in the Gulf appeared in a sample MOD09GA product during a whiting event and shown in Figure 4a–h. Figure 4f–h show whiting events with high intensities of bands 3 (blue region of the spectra) and 4 (green region of the spectra). Figure 4e shows a clear water sample with a low intensities of all bands. Shallow areas at the northwest section of the gulf show high reflections in bands 1 (red region of the spectra) and 2 (NIR region of the spectra), as shown in Figure 4a. Overall, the images show variations in terms of relative intensities of bands 3 and 4 and bands 1 and 2. Meanwhile, the intensities of bands 5, 6 and 7 did not exhibit significant variations.

2.3. Object-Based Analysis and Image Segmentation Optimisation

A traditional per-pixel classification approach is useful for feature extraction when the targets of interest are smaller than the spatial resolution of the remotely sensed data [41]. This approach only considers the spectral properties of each pixel and disregards any spatial or contextual information related to the classified pixel [42]. Geographic object-based image analysis (GEOBIA) has been extensively used in classifying very-high spatial resolution data as an alternative to a pixel-based approach. GEOBIA works by assessing spatially neighbouring groups of pixels rather than individual pixels [43]. GEOBIA is not only limited to high resolution images because the approach is not spatial-resolution dependent; therefore, it can be applied to different resolutions if the sizes of the intended objects are compatible with the spatial resolution of the images [44,45]. Thus, GEOBIA has been successfully adopted and implemented to classify MODIS time series data in different applications [44,45,46,47,48,49,50].
The generic GEOBIA framework can be divided into (1) image segmentation, which is the process of generating homogenous and nonoverlapping image objects/segments from image pixels and (2) image object/segments classification [51]. The multiresolution image segmentation algorithm (MRS) [52], which is one of the most used algorithms, is applied to the time series MODIS data. The MRS is a bottom-up region-growing algorithm that commences with pixels as individual segments. MRS is goverened by three main parameters, namely, (a) scale, (b) shape/color weight and (c) compactness/smoothness weight. This algorithm merges neighbouring pixels in each successive step on the basis of homogeneity (shape and compactness), which describes the similarity of contiguous objects. The degree of fitting, which is a value determined by the scale parameter defined by an analyst, is measured in each merging procedure. Moreover, the merge is performed if the degree of fitting is less than the minimum degree of fitting [53]. The scale parameter is one of the most critical parameters in the segmentation process. It profoundly influences resultant image objects and subsequent classification steps because it controls the size of image-generated objects [54,55]. Selecting high-scale values generate large image objects (undersegmentation), whereas selecting small-scale values yield small image objects (oversegmentation). Thus, the utilisation of an optimisation technique to find the optimum scale parameters for delineating whiting is vital to avoid the subjectivity of using a trial and error visual approach.
Several unsupervised segmentation quality measures have been utilised in the literature to identify the best scale parameters [56,57,58,59]. In the current study, the performance of two unsupervised segmentation quality measures, namely, the objective function (OF) [60] and the F-measure [57] were compared to find an optimum scale value that accurately delineates whiting. Both measures adopt oversegmentation and undersegmentation metrics by using the values of weighted variance and spatial autocorrelations (Moran’s I). They are expressed in Equations (1)–(5).
OF   =   WV norm +   MI norm   ,
F - measure   =   2   MI norm   ×   WV norm MI norm   +   WV norm ,  
WV   = i = 1 n a i . v i   i = 1 n a i   ,
MI   = n i = 1 n j = 1 n w i , j ( z i z ¯ ) ( z j z ¯ ) ( i = 1 n ( z i z ¯ ) 2 )   ( i j w i ,   j ) ,
where MInorm and WVnorm are the normalised Moran’s I and weighted variance, respectively. In Equation (3), WV denotes the weighted variance and a i and v i are the area and variance of image object/segment (generated by MRS) i, respectively. In Equation (4), n symbolises the total number of objects, z i and z j are the means of the spectral value of image objects i (Oi) and j (Oj), respectively, z ¯ is the mean spectral value of the total objects in a specific band and w i , j is a spatial proximity measure between image objects i (Oi) and j (Oj) in which nearby image objects are defined as 1 and other objects are considered as 0. The normalised function for WV and MI can be expressed using the following equation
F ( X )   =   X max ·   X   X max ·   X min ,
where X, Xmax and Xmin are the original, maximum and minimum values of the weighted variance or Moran’s I for a spectral band, respectively. The highest values of the OF or F-measure indicate high segmentation quality and the optimum scale is identified as the scale that achieves image objects with the highest OF or F-measure values. After the selection of an appropriate scale to delineate whiting, various spectral features and indices were computed and investigated to find the most relevant features to be used in the classification phase.

2.4. Feature Selection

One of the strongest characteristics of GEOBIA is it enables the generation of hundreds of features or attributes for analysis (e.g., spectral bands and indices, textural, geometric and contextual attributes). The extraction and utilisation of a considerable number of features (e.g., variables or attributes) in the analysis are computationally intensive; hence, they can negatively affect classification accuracy [61]. Therefore, FS, which involves the selection of an essential feature subset from an enormous amount of generated features, is a decisive step used in image analysis procedures. In addition, FS can achieve an equivalent or higher classification accuracy than the original feature space and improve the efficiency of GEOBIA [62,63,64,65,66]. FS methods may be generally grouped into three classes, namely, filter, wrapper and embedded methods [64,65]. The filter method is considered the simplest and fastest method among the three classes. It utilises certain statistical measures (e.g., correlation coefficients, variance, chi-square test measures and ANOVA F-values) to rank and select relevant features without using any learning algorithms [64,67,68,69]. By contrast, the wrapper method adopts a classification algorithm as part of the evaluation process to classify the training data and assess the results. The wrapper method selects the most significant feature subset that produces the highest classification accuracy [70,71]. Finally, the embedded method exhibits a trade-off between the filter and the wrapper methods. This method is considered feature ranking because features are selected during the construction of a classification model without further evaluating the selected feature subset [65]. Comparisons of various FS techniques are available in the literature [64,72,73,74].
A wrapper method that combines CFS and the naïve Bayes classifier was used in this study to find the most relevant features for extracting whiting features. CFS has recently been successfully applied to FS and has outperformed various FS techniques with GEOBIA [63,75].

2.4.1. CFS

CFS is a popular approach that uses a search algorithm with a heuristic evaluation function to assess the merit of feature subsets [76,77]. It measures the worth of each feature to predict the class label along with the intercorrelation level among features [78]. A heuristic evaluation function is designed on the basis of the hypothesis that superior feature subsets encompass correlated features with classes though they remain uncorrelated with one another [76,79]. Merits (heuristics) can be formalised using the following formula
Merits = k r ¯ cf k + k ( k   1 ) r ¯ ff ,
where k denotes the number of features, f indicates the feature, r ¯ cf symbolises the mean feature correlation with a class and r ¯ ff is the average intercorrelation among subset features.

2.4.2. Feature Acquisition and Computation

In general, various spectral indices and bands may be used to map water surface features and water bodies. In this study, numerous attributes, such as mean spectral reflectance, standard deviation and spectral indices, which have been previously reported in the literature, were computed for the image objects of several images for FS. Table 1 lists the 32 attributes that were examined using CFS.
The present study used CFS and the naïve Bayes classifier to determine the most relevant feature subset for extracting whiting events from the multitemporal MODIS data. Then, training and testing samples were prepared from the image objects of various whiting events. Table 2 shows the dates and number of samples that were selected for FS. These samples were normalised to a scale from 0 to 1 and then split 70% for training and 30% for testing.
The optimised image objects with the selected significant features via CFS were ultimately used to classify whiting using various tree-based classification algorithms.

2.5. Boosting Decision Tree Classification

Boosting is an ensemble machine learning algorithm that is used to improve the accuracy of a classifier by decreasing the classification algorithm’s sensitivity to noise and labelling errors in the training datasets [92,93]. The performance of numerous classifiers (weak learners) that are learned by resampled versions of the training samples is combined to improve classification accuracy. Boosting DTs are ensemble methods that use multiple iterations of DT classifiers. A boosted DT called the AdaBoost algorithm was adopted in this study.
AdaBoost [92], also known as adaptive boosting, is a generic iterative supervised learning algorithm that combines multiple classifiers (weak classifiers) to obtain high accuracy. The AdaBoost algorithm chooses training samples on the basis of adaptive resampling by selecting misclassified datasets produced by a previous classifier. The erroneously classified samples in a prior iteration are selected more often than correctly classified samples. Furthermore, new DT models are forced to focus on the misclassified samples and minimise the errors of the former trees [94,95]. Misclassified training samples are given increased weights in each iteration; thus, the classifier can improve its performance in new datasets. Ultimately, all trees (not only a final tree or trees) are incorporated because the additive model is designed such that combinations of all trees can give the optimum solution [96].
Given surface water features with diverse spectral values, training image objects were selected to represent five classes with different spectral responses, namely, whiting, low whiting, clear water, interlayer (a class with subtle differences between low whiting and clear water), apparent sediment (shallow areas next to the shore) and apparent green (potential true green-coloured algae). The prepared samples with the most significant features selected by CFS were split 70% for training and 30% for assessing the performance of the AdaBoost classification model. Using the verified prediction model developed by the AdaBoost algorithm, the multitemporal image objects were classified into the aforementioned five classes by rule-based classification. The performance of the AdaBoost classification technique was compared with that of various tree-based classification algorithms, such as random forest, the single DT and the gradient boosted decision tree. Classification results were assessed in this study by computing the overall accuracy (OA), the Kappa coefficient (KC) from the confusion matrix. While the OA of the classification indicates the percentage of correctly classified image objects, the KC statistically measures and analyses the degree of agreement between the classified and reference image objects [97].

3. Results & Discussion

The results of the spatiotemporal mapping of whiting in the Gulf using the MODIS time series images and a generic ensemble tree-based model are presented in this section. Furthermore, the frequency, seasonality, duration and geographic distribution and extent of whiting events in the Gulf are also documented and analysed.

3.1. Whiting Temporal Pattern in the Gulf

The primary objective of this study was to map the spatial and temporal extents of whiting in the Gulf for 16 years by using satellite images. Therefore, a total of 5800 image scenes of the Gulf were inspected prior to the analysis to detect the existence of whiting. Cloud-free image sets obtained for the period of 2002 to 2018 were used to generate generic statistics on the seasonality, frequency and duration of whiting events in the Gulf, as shown in Table 3.
The results showed that whiting events reoccurred in the region exclusively during the winter season (November to March). Figure 5 enumerates the number of events per month/year and the number of days per month/year where whiting events occurred during the study period (2002–2018). The frequency of whiting events during the study period ranged from one to two events per winter month to two to eight events per year. The highest frequency of whiting events was observed in February. The total number of whiting days ranged from 2 to 11 days per month to 8 to 34 days per year. Furthermore, the duration of individual events ranged from 2 to 8 days. The total number whiting days during the past 16 years was approximately 289 days or approximately 7% of the total number of days or 16% of the total days (November to March).

3.2. Results of the Integrated GEOBIA Approach

3.2.1. Results of Image Segmentation

Unsupervised image segmentation quality measures based on the OF and the F-measure were utilised to find the optimal scale for whiting event delineation by varying the scale parameter. The mean values of NIR image objects are frequently used in image segmentation assessments to compute undersegmentation and oversegmentation metrics. However, the NIR reflectance band was unsuitable for assessing the segmentation of images with whiting occurrences because the whiting phenomenon does not exist on this band. Therefore, the mean values and standard deviation of the mean blue spectral image objects were used to compute the OF and the F-score values (Equations (1)–(5)). Table 4 provides a sample of the OF and the F-measure computations for a single image acquired on 2 March 2003. The highest OF and F-measure values were on a scale value of 40. Thus, scale 40 was selected as the optimum scale value. The values of the shape and compactness parameters were set as 0.5 and 0.1, respectively.

3.2.2. Results of FS and Analysis

FS reduces the dimensionality of data, diminishes the complexity of classification models, minimises overfitting and accelerates the process. CFS was selected among various FS methods because of its successful implementation in several remote sensing applications. As listed in Table 1, various spectral bands and indices were examined to determine the most significant features for detecting whiting event occurrences from a large number of MODIS images. Given the considerable variation between the examined features in terms of range and numerical value, all the data were normalised to a scale ranging from 0 to 1. The CFS algorithm with the best-first search strategy [98] was implemented along with the naïve Bayes classifier to evaluate the worth of the selected features. Accuracy was assessed via tenfold stratified cross-validation on the selected training data. According to the results, the best-first eight significant features were the green (G), NDGB, FAI, CI, CI2, CI869, BRI and SD-R. Therefore, combinations of each pair of the best-first selected features were selected for further examination. Figure 6 shows the possible combinations of the selected features and their corresponding overall accuracy and kappa coefficient. These tests reveal that the utilisation of the NDGB and the mean green attribute values or the NDGB and CI of the image objects were similarly excellent features for accurately mapping whiting using the MODIS data. This finding was typically consistent with the spectral analysis of various surface features shown in Figure 4. In this study, the NDGB and the green band were selected for further analysis.

3.2.3. Classification Results

Given the lack of historical records of whiting events in the Gulf, the AdaBoost model was trained and tested with samples selected on the basis of satellite image inspection and the recommendation of Wells and Illing (1962). Figure 7 depicts the generated AdaBoost DT model from the most significant features, namely, the mean green and the NDGB, to extract whiting event features from satellite images with moderate spatial resolution. Optimised time series image objects were eventually classified into six classes (namely, clear water, interlayer, low whiting, whiting, apparent green and apparent sediments) by the developed AdaBoost model and GEOBIA rule-based classification. Whiting water with a high intensity can be simply identified with the generic developed model by rule-based classification when the slope (NDGB) is less than 0.1 and the intensity of the green band is greater than 1500 in the MODIS scene. Whiting with a low intensity can be recognised when the slope (NDGB) is less than 0.1 and the mean of the green band range is between 750–1500.
Figure 8a–h present four classification results together with the original images acquired on 2 March 2003, 4 March 2012, 10 February 2016 and 3 October 2018. Figure 8a–d show extensive whiting along the coasts of UAE and Qatar, whereas Figure 8e,f show low levels of whiting. As expected from the classification model, Figure 8g,h represent clear water without any significant whiting.
The performance of the developed AdaBoost classification model was compared with that of the three tree-based classification algorithms. The AdaBoost classification model gave superior classification outcomes compared to the other models, as shown in Table 5.

3.3. Spatial Distribution of Whiting in the Gulf

Generalisation of the occurrence pattern of whiting in the Gulf can be challenging because of the existence of clouds in the MODIS satellite images, especially during winter. Cloud-free images representing the highest concentration of whiting events from each year were selected to compute the peak/maximum area covered by whiting and to determine the spatial distribution of whiting in the region. A total of 17 MODIS images (one image per year) were classified with the developed model. The extent of the areas covered by whiting were then computed, as shown in Table 6. The maximum coverage of whiting occurred in March 2012, followed by March 2003, whereas the minimum coverage of assessed whiting events was recorded in December 2010. The classification results were used to identify the common areas subjected to whiting. To statistically identify significant hot spots (High spatial frequency of whiting events overs 17 years), statistical analysis of spatial clustering was carried out using Optimized Hot Spot Analysis (Getis-Ord Gi*) [99,100]. This tool automatically aggregates whiting events and identifies statistically significant spatial clusters where the focus is on presence or absence of each whiting event rather than a measured attribute associated with whiting events. Figure 9 shows significant whiting spatiotemporal clusters, in the semi-enclosed gulf, with various levels of confidence, areas that are statistically significant at the 99 percent confidence level showed that whiting events were the frequent in the southwest sections of the Gulf and along the coasts of UAE, Qatar, Bahrain and opposite the coast of Al Jubail in Saudi Arabia.

4. Conclusions

Previous studies mapping whiting events from satellite images, specifically in the Bahama Banks and in the coast of Southwest Florida, relied on the manual delineation of whiting on the basis of spatial contrast, spectral behaviour of whiting and derived spectral indices. Considering the limited studies on whiting event mapping in the Gulf, the present study aimed to document the spatial extent and the seasonal variability of whiting events in the Gulf between 2002 and 2018 using MODIS data. This study acquired and analysed extensive daily data for mapping and documenting, documented spatiotemporal distribution and presented an effective model (integrated CFS, Adaboost and rule-based classification) for detecting and classifying whiting in the Gulf using the MODIS data. The results of FS showed that the combination of the mean of the green band and the NDGB or the combination of the NDGB and CI were the most significant feature for detecting the brightness of in-water features compared with all the examined features in the classification. This study used various tree-based machine learning classifiers, namely, rule-based classification based on a single DT, GBDT, RF and AdaBoost, to classify the optimised multitemporal image objects. The results showed that the rule-based classification based on AdaBoost DT outperformed the supervised tree-based GEOBIA classifiers. Therefore, this study adopted the AdaBoost classification model to find a generic model for distinguishing objects of whiting water directly and classifying time series image objects by rule-based classification.
The adopted model showed an outstanding and expeditious approach to extracting and characterising whiting events quantitatively from time series images. Whiting events in the Gulf occurred during the winter season (November to March) and were extensively located in the southwestern section of the Gulf, mainly along the UAE coast. During the study period (2002–2018), the whiting events occurred exclusively for 5 to 34 days per year and covered areas ranging from 12,000 km2 to 60,000 km2. These events require further investigations for in-situ measurements and laboratory analysis, on the basis of the common spatial distribution of whiting. Therefore, whiting in the Arabian Gulf merits further attention from the scientific community to examine biophysical, biogeochemical and environmental factors that may reveal the causes of the whiting occurrences.

Author Contributions

Conceptualization, A.S., R.A.-R. and M.B.A.G.; Investigation, A.S., R.A.-R. and M.B.A.G.; Methodology, A.S., R.A.-R., M.B.A.G. and C.F.; Supervision, A.S., R.A.-R. and M.B.A.G.; Writing—original draft, A.S., R.A.-R. and M.B.A.G.; Writing—review & editing, A.S., R.A.-R., M.B.A.G., C.F. and S.A.-M.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Long, J.S.; Hu, C.; Robbins, L.L.; Byrne, R.H.; Paul, J.H.; Wolny, J.L. Optical and biochemical properties of a southwest Florida whiting event. Estuar. Coast. Shelf Sci. 2017, 196, 258–268. [Google Scholar] [CrossRef]
  2. Morse, J.W.; He, S. Influence of T, S and PCO2 on the homogeneous nucleation of calcium carbonate from seawater. Implications for whiting formation. Mar. Chem. 1993, 41, 291–298. [Google Scholar] [CrossRef]
  3. Watkins, J.M.; Rudstam, L.G.; Crabtree, D.L.; Walsh, M.G. Is reduced benthic flux related to the Diporeia decline? Analysis of spring blooms and whiting events in Lake Ontario. J. Great Lakes Res. 2013, 39, 395–403. [Google Scholar] [CrossRef]
  4. Thompson, J.B.; Schultze-Lam, S.; Beveridge, T.J.; Des Marais, D.J. Whiting events: Biogenic origin due to the photosynthetic activity of cyanobacterial picoplankton. Limnol. Oceanogr. 1997, 42, 133–141. [Google Scholar] [CrossRef] [PubMed]
  5. Wurgaft, E.; Steiner, Z.; Luz, B.; Lazar, B. Evidence for Inorganic Precipitation of CaCO3 on Suspended Solids in the Open Water of the Red Sea. Mar. Chem. 2016, 186, 145–155. [Google Scholar] [CrossRef]
  6. Friedman, G.M. Biochemical and ultrastructural evidence for the origin of whitings: A biologically induced calcium carbonate precipitation mechanism: Comment and reply. Geology 1993, 21, 287–288. [Google Scholar] [CrossRef]
  7. Bloch, R.; Littman, H.Z.; Elazari-Volcani, B. Occasional whiteness of the dead sea. Nature 1944, 154, 402–403. [Google Scholar] [CrossRef]
  8. Bustos-Serrano, H.; Morse, J.W.; Millero, F.J. The formation of whitings on the Little Bahama Bank. Mar. Chem. 2009, 113, 1–8. [Google Scholar] [CrossRef]
  9. Bathurst, R.G.C. Carbonate Sediments and Their Diagenesis, 2nd ed.; Developments in Sedimentology; Elsevier: Amsterdam, The Netherlands, 1975; Volume 12, p. 658. [Google Scholar]
  10. Long, J.S.; Hu, C.; Wang, M. Long-term spatiotemporal variability of southwest Florida whiting events from MODIS observations. Int. J. Remote Sens. 2018, 39, 906–923. [Google Scholar] [CrossRef]
  11. Robbins, L.L.; Blackwelder, P.L. Biochemical and ultrastructural evidence for the origin of whiting: A biologically induced calcium carbonate precipitation mechanism. Geology 1992, 20, 464–468. [Google Scholar] [CrossRef]
  12. Dierssen, H.M.; Zimmerman, R.C.; Burdige, D.J. Optics and remote sensing of Bahamian carbonate sediment whitings and potential relationship to wind-driven Langmuir circulation. Biogeosciences 2009, 6, 487–500. [Google Scholar] [CrossRef] [Green Version]
  13. Broecker, W.S.; Takahashi, T. Calcium carbonate precipitation on the Bahama Banks. J. Geophys. Res. 1966, 71, 1575–1602. [Google Scholar] [CrossRef]
  14. Boss, S.K.; Neumann, A.C. Physical versus chemical processes of “whiting” formation in the Bahamas. Carbonates Evaporites 1993, 8, 135–148. [Google Scholar] [CrossRef]
  15. Shinn, E.A.; Steinen, R.P.; Lidz, B.H.; Swart, P.K. Whitings, a Sedimentologic Dilemma. J. Sediment. Petrol. 1989, 59, 147–161. [Google Scholar] [CrossRef]
  16. Cloud, P.E. Environment of Calcium Carbonate Deposition West of Andros Island Bahamas. Geol. Surv. Prof. Pap. 1962, 350, 170. [Google Scholar]
  17. Riding, R.E.; Awramik, S.M. Microbial Sediments; Springer: Berlin, Germany, 2000; ISBN 978-3-642-08275-7. [Google Scholar]
  18. Ohlendorf, C.; Sturm, M. Precipitation and Dissolution of Calcite in a Swiss High Alpine Lake. Arctic Antarct. Alp. Res. 2001, 33, 410–417. [Google Scholar] [CrossRef]
  19. Otsuki, A.; Wetzel, R.G. Calcium and total alkalinity budgets and calcium carbonate precipitation of a small hard-water lake. Arch. Hydrobiol. 1974, 73, 14–30. [Google Scholar] [CrossRef]
  20. Morse, J.W.; Gledhill, D.K.; Millero, F.J. CaCO3 precipitation kinetics in waters from the Great Bahama Bank: Implications for the relationship between Bank hydrochemistry and whitings. Geochim. Cosmochim. Acta 2003, 67, 2819–2826. [Google Scholar] [CrossRef]
  21. Shinn, E.A.; Kendall, C.G.S.C. Back to the Future. Sediment. Rec. 2011, 9, 4–9. [Google Scholar] [CrossRef]
  22. Lidz, B.; Gibbons, H. Research on Whitings (Floating Patches of Calcium Carbonate Mud) Leads to Possible Explanation of Immense Middle East Oil Deposits. Available online: https://soundwaves.usgs.gov/2008/07/research.html (accessed on 30 September 2018).
  23. Whitton, B.A. Ecology of Cyanobacteria II: Their Diversity in Space and Time; Springer: Dordrecht, The Netherlands, 2012; ISBN 978-94-007-3854-6. [Google Scholar]
  24. Robbins, L.L.; Tao, Y.; Evans, C.A. Temporal and spatial distribution of whitings on Great Bahama Bank and a new lime mud budget. Geology 1997, 25, 947–950. [Google Scholar] [CrossRef]
  25. Long, J.; Hu, C.; Robbins, L. Whiting events in SW Florida coastal waters: A case study using MODIS medium-resolution data. Remote Sens. Lett. 2014, 5, 539–547. [Google Scholar] [CrossRef]
  26. Strong, A.E.; Eadie, B.J. Satellite observations of calcium carbonate precipitation in the Great Lakes. Limnol. Ocean. 1978, 23, 877–887. [Google Scholar] [CrossRef]
  27. Heine, I.; Brauer, A.; Heim, B.; Itzerott, S.; Kasprzak, P.; Kienel, U.; Kleinschmit, B. Monitoring of calcite precipitation in hardwater lakes with multi-spectral remote sensing archives. Water 2017, 9, 15. [Google Scholar] [CrossRef]
  28. Millero, J. The carbonate chemistry of grand bahama bank waters: After 18 years another look. J. Geophys. Res. 1984, 89, 3604–3614. [Google Scholar]
  29. Long, J.S. Whiting Events Off Southwest Florida: Remote Sensing and Field Observations. Ph.D. Dissertation, University of South Florida, Tampa, FL, USA, 2016. [Google Scholar]
  30. Lloyd, R.A. Remote Sensing of Whitings in the Bahamas. Master’s Thesis, University of South Florida, Tampa, FL, USA, 2012. [Google Scholar]
  31. Tao, Y. Whitings on the Great Bahama Bank: Distribution in Space and Time Using Space Shuttle Photographs. Ph.D. Dissertation, University of South Florida, Tampa, FL, USA, 1994. [Google Scholar]
  32. Balch, W.M.; Gordon, H.R.; Bowler, B.C.; Drapeau, D.T.; Booth, E.S. Calcium carbonate measurements in the surface global ocean based on Moderate-Resolution Imaging Spectroradiometer data. J. Geophys. Res. C Ocean. 2005, 110, 1–21. [Google Scholar] [CrossRef]
  33. Gordon, R.; Boynton, G.C.; Balch, W.M.; Harbour, D.S.; Smyth, T.J.; Bay, W.B. Retrieval of Coccolithophore from SeaWiFS Imagery Calcite Concentration radiance. Geophys. Res. Lett. 2001, 28, 1587–1590. [Google Scholar] [CrossRef]
  34. Mitchell, C.; Hu, C.; Bowler, B.; Drapeau, D.; Balch, W.M. Estimating Particulate Inorganic Carbon Concentrations of the Global Ocean From Ocean Color Measurements Using a Reflectance Difference Approach. J. Geophys. Res. Ocean. 2017, 122, 8707–8720. [Google Scholar] [CrossRef]
  35. Wells, A.J.; Illing, L.V. Present-day precipitation of calcium carbonate in the Persian Gulf. Dev. Sedimentol. 1964, 1, 429–435. [Google Scholar] [CrossRef]
  36. Shanableh, A.; Imteaz, M.; Hamad, K.; Omar, M.; Merabtene, T.; Siddique, M. Potential impact of global warming on whiting in a semi-enclosed gulf. Int. J. Glob. Warm. 2017, 13, 411–425. [Google Scholar]
  37. Shanableh, A.; Al-Ruzouq, R.; Al-Khayyat, G. Assessing the Spatial and Temporal Capacity of a Semi-Enclosed Gulf to Absorb and Release CO2 Using GIS and Remote Sensing. In Proceedings of the Global Civil Engineering Conference, Kuala Lumpur, Malaysia, 25–28 July 2017; pp. 1161–1173. [Google Scholar]
  38. Sheppard, C.; Al-Husiani, M.; Al-Jamali, F.; Al-Yamani, F.; Baldwin, R.; Bishop, J.; Benzoni, F.; Dutrieux, E.; Dulvy, N.K.; Durvasula, S.R.V.; et al. Environmental Concerns for the Future of Gulf Coral Reefs; Springer: Dordrecht, The Netherlands, 2012; p. E1. ISBN 978-94-007-3007-6. [Google Scholar]
  39. Kaempf, J.; Sadrinasab, M. The circualtion of the Persian Gulf: A numerical study. Ocean Sci. 2006, 2, 27–41. [Google Scholar] [CrossRef]
  40. Ben-Hasan, A.; Walters, C.; Christensen, V.; Al-Husaini, M.; Al-Foudari, H. Is reduced freshwater flow in Tigris-Euphrates rivers driving fish recruitment changes in the Northwestern Arabian Gulf? Mar. Pollut. Bull. 2018, 129, 1–7. [Google Scholar] [CrossRef] [PubMed]
  41. Blaschke, T.; Hay, G.J.; Kelly, M.; Lang, S.; Hofmann, P.; Addink, E.; Queiroz Feitosa, R.; van der Meer, F.; van der Werff, H.; van Coillie, F.; et al. Geographic Object-Based Image Analysis—Towards a new paradigm. ISPRS J. Photogramm. Remote Sens. 2014, 87, 180–191. [Google Scholar] [CrossRef]
  42. Makinde, E.O.; Salami, A.T.; Olaleye, J.B.; Okewusi, O.C. Object Based and Pixel Based Classification Using Rapideye Satellite Imager of ETI-OSA, Lagos, Nigeria. Geoinform. FCE CTU 2016, 15, 59. [Google Scholar] [CrossRef]
  43. Petropoulos, G.P.; Vadrevu, K.P.; Kalaitzidis, C. Spectral angle mapper and object-based classification combined with hyperspectral remote sensing imagery for obtaining land use/cover mapping in a Mediterranean region. Geocarto Int. 2013, 28, 114–129. [Google Scholar] [CrossRef]
  44. da Silva Junior, C.A.; Nanni, M.R.; de Oliveira-Júnior, J.F.; Cezar, E.; Teodoro, P.E.; Delgado, R.C.; Shiratsuchi, L.S.; Shakir, M.; Chicati, M.L. Object-based image analysis supported by data mining to discriminate large areas of soybean. Int. J. Digit. Earth 2019, 12, 270–292. [Google Scholar] [CrossRef]
  45. Bisquert, M.; Bégué, A.; Deshayes, M. Object-based delineation of homogeneous landscape units at regional scale based on modis time series. Int. J. Appl. Earth Obs. Geoinf. 2015, 37, 72–82. [Google Scholar] [CrossRef]
  46. Cano, E.; Denux, J.P.; Bisquert, M.; Hubert-Moy, L.; Chéret, V. Improved forest-cover mapping based on MODIS time series and landscape stratification. Int. J. Remote Sens. 2017, 38, 1865–1888. [Google Scholar] [CrossRef]
  47. Vintrou, E.; Desbrosse, A.; Bégué, A.; Traoré, S.; Baron, C.; Lo Seen, D. Crop area mapping in West Africa using landscape stratification of MODIS time series and comparison with existing global land products. Int. J. Appl. Earth Obs. Geoinf. 2012, 14, 83–93. [Google Scholar] [CrossRef]
  48. Bontemps, S.; Bogaert, P.; Titeux, N.; Defourny, P. An object-based change detection method accounting for temporal dependences in time series with medium to coarse spatial resolution. Remote Sens. Environ. 2008, 112, 3181–3191. [Google Scholar] [CrossRef]
  49. Bisquert, M.; Bégué, A.; Deshayes, M.; Ducrot, D. Environmental evaluation of MODIS-derived land units. GISci. Remote Sens. 2017, 54, 64–77. [Google Scholar] [CrossRef]
  50. Bellón, B.; Bégué, A.; Lo Seen, D.; de Almeida, C.A.; Simões, M. A remote sensing approach for regional-scale mapping of agricultural land-use systems based on NDVI time series. Remote Sens. 2017, 9, 600. [Google Scholar] [CrossRef]
  51. Kim, M.; Madden, M.; Warner, T.T. Forest Type Mapping using Object-specific Texture Measures from Multispectral Ikonos Imagery: Segmentation Quality and Image Classification Issues. Photogramm. Eng. Remote Sens. 2009, 75, 819–829. [Google Scholar] [CrossRef]
  52. Baatz, M.; Sch, A. Multiresolution Segmentation: An optimization approach for high quality multi-scale image segmentation. J. Photogramm. Remote Sens. 2004, 58, 239–258. [Google Scholar]
  53. Grybas, H.; Melendy, L.; Congalton, R.G. A comparison of unsupervised segmentation parameter optimization approaches using moderate- and high-resolution imagery. GISci. Remote Sens. 2017, 54, 515–533. [Google Scholar] [CrossRef]
  54. Kim, M.; Warner, T.A.; Madden, M.; Atkinson, D.S. Multi-scale GEOBIA with very high spatial resolution digital aerial imagery: Scale, texture and image objects. Int. J. Remote Sens. 2011, 32, 2825–2850. [Google Scholar] [CrossRef]
  55. Hussain, M.; Chen, D.; Cheng, A.; Wei, H.; Stanley, D. Change detection from remotely sensed images: From pixel-based to object-based approaches. ISPRS J. Photogramm. Remote Sens. 2013, 80, 91–106. [Google Scholar] [CrossRef]
  56. Gao, Y.; Mas, J.F.; Kerle, N.; Pacheco, J.A.N. Optimal region growing segmentation and its effect on classification accuracy. Int. J. Remote Sens. 2011, 32, 3747–3763. [Google Scholar] [CrossRef]
  57. Johnson, B.; Bragais, M.; Endo, I.; Magcale-Macandog, D.; Macandog, P. Image Segmentation Parameter Optimization Considering Within- and Between-Segment Heterogeneity at Multiple Scale Levels: Test Case for Mapping Residential Areas Using Landsat Imagery. ISPRS Int. J. Geo-Inf. 2015, 4, 2292–2305. [Google Scholar] [CrossRef] [Green Version]
  58. Drǎguţ, L.; Csillik, O.; Eisank, C.; Tiede, D. Automated parameterisation for multi-scale image segmentation on multiple layers. ISPRS J. Photogramm. Remote Sens. 2014, 88, 119–127. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, X.; Xiao, P.; Feng, X. An unsupervised evaluation method for remotely sensed imagery segmentation. IEEE Geosci. Remote Sens. Lett. 2012, 9, 156–160. [Google Scholar] [CrossRef]
  60. Espindola, G.M.; Camara, G.; Reis, I.A.; Bins, L.S.; Monteiro, A.M. Parameter selection for region-growing image segmentation algorithms using spatial autocorrelation. Int. J. Remote Sens. 2006, 27, 3035–3040. [Google Scholar] [CrossRef]
  61. Laliberte, A.S.; Browning, D.M.; Rango, A. A comparison of three feature selection methods for object-based classification of sub-decimeter resolution UltraCam-L imagery. Int. J. Appl. Earth Obs. Geoinf. 2012, 15, 70–78. [Google Scholar] [CrossRef]
  62. Ma, L.; Cheng, L.; Li, M.; Liu, Y.; Ma, X. Training set size, scale, and features in Geographic Object-Based Image Analysis of very high resolution unmanned aerial vehicle imagery. ISPRS J. Photogramm. Remote Sens. 2015. [Google Scholar] [CrossRef]
  63. Al-Ruzouq, R.; Shanableh, A.; Barakat, A.; Gibril, M.; AL-Mansoori, S. Image Segmentation Parameter Selection and Ant Colony Optimization for Date Palm Tree Detection and Mapping from Very-High-Spatial-Resolution Aerial Imagery. Remote Sens. 2018, 10, 1413. [Google Scholar] [CrossRef]
  64. Georganos, S.; Grippa, T.; Vanhuysse, S.; Lennert, M.; Shimoni, M.; Kalogirou, S.; Wolff, E. Less is more: Optimizing classification performance through feature selection in a very-high-resolution remote sensing object-based urban application. GISci. Remote Sens. 2018, 55, 221–242. [Google Scholar] [CrossRef]
  65. Zhou, Y.; Chen, Y.; Feng, L.; Zhang, X.; Shen, Z.; Zhou, X. Supervised and Adaptive Feature Weighting for Object-Based Classification on Satellite Images. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3224–3234. [Google Scholar] [CrossRef]
  66. Benediktsson, J.A.; Pesaresi, M.; Arnason, K. Classification and feature extraction for remote sensing images from urban areas based on morphological transformations. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1940–1949. [Google Scholar] [CrossRef] [Green Version]
  67. Dash, M.; Liu, H. Feature selection for classification. Intell. Data Anal. 1997, 1, 131–156. [Google Scholar] [CrossRef]
  68. Löw, F.; Michel, U.; Dech, S.; Conrad, C. Impact of feature selection on the accuracy and spatial uncertainty of per-field crop classification using Support Vector Machines. ISPRS J. Photogramm. Remote Sens. 2013, 85, 102–119. [Google Scholar] [CrossRef]
  69. Guyon, I.; Elisseeff, A. An introduction to variable and feature selection. J. Mach. Learn. Res. 2003, 3, 1157–1182. [Google Scholar]
  70. Duro, D.C.; Franklin, S.E.; Dubé, M.G. Multi-scale object-based image analysis and feature selection of multi-sensor earth observation imagery using random forests. Int. J. Remote Sens. 2012, 33, 4502–4526. [Google Scholar] [CrossRef]
  71. Hu, M.; Wu, F. Filter-wrapper hybrid method on feature selection. In Proceedings of the 2010 Second WRI Global Congress on Intelligent Systems, Wuhan, China, 16–17 December 2010; pp. 98–101. [Google Scholar] [CrossRef]
  72. Ma, L.; Fu, T.; Blaschke, T.; Li, M.; Tiede, D.; Zhou, Z.; Ma, X.; Chen, D. Evaluation of Feature Selection Methods for Object-Based Land Cover Mapping of Unmanned Aerial Vehicle Imagery Using Random Forest and Support Vector Machine Classifiers. ISPRS Int. J. Geo-Inf. 2017, 6, 51. [Google Scholar] [CrossRef]
  73. Gumus, E.; Kirci, P. Selection of spectral features for land cover type classification. Expert Syst. Appl. 2018, 102, 27–35. [Google Scholar] [CrossRef]
  74. Hamedianfar, A.; Barakat, A.M. Gibril Large-scale urban mapping using integrated geographic object-based image analysis and artificial bee colony optimization from worldview-3 data. Int. J. Remote Sens. 2019. [Google Scholar] [CrossRef]
  75. Ridha, M.; Pradhan, B. Catena An improved algorithm for identifying shallow and deep-seated landslides in dense tropical forest from airborne laser scanning data. Catena 2018, 167, 147–159. [Google Scholar] [CrossRef]
  76. Hall, M.A.; Smith, L.A. Feature subset selection: A correlation based filter approach. In International Conference on Neural Information Processing and Intelligent Information Systems; Springer: Berlin, Germany, 1997; pp. 855–858. [Google Scholar]
  77. Hall, M.A. Correlation-based feature selection of discrete and numeric class machine learning. In Proceedings of the 17th International Conference on Machine Learning, San Francisco, CA, USA, 29 June–2 July 2000; pp. 359–366. [Google Scholar]
  78. Hall, M.A.; Smith, L.A. Practical feature subset selection for machine learning. In Proceedings of the 21st Australasian Computer Science Conference ACSC’98, Perth, Australia, 4–6 February 1998; Volume 98, pp. 181–191. [Google Scholar]
  79. Hall, M.A.; Holmes, G. Benchmarking Attribute Selection Techniques for Discrete Class Data Mining. IEEE Trans. Knowl. Data Eng. 2003, 15, 1437–1447. [Google Scholar] [CrossRef]
  80. Trimble, T. ECognition Developer 8.7 Reference Book; Trimble Germany GmbH: Munich, Germany, 2011; pp. 319–328. [Google Scholar]
  81. Curran, P. Multispectral remote sensing of vegetation amount. Prog. Phys. Geogr. 1980, 4, 315–341. [Google Scholar] [CrossRef]
  82. Mcfeeters, S.K. The use of the Normalized Difference Water Index ( NDWI ) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  83. Xu, H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 2006, 27, 3025–3033. [Google Scholar] [CrossRef]
  84. Hu, C. A novel ocean color index to detect floating algae in the global oceans. Remote Sens. Environ. 2009, 113, 2118–2129. [Google Scholar] [CrossRef]
  85. Hu, C.; Lee, Z.; Franz, B. Chlorophyll a algorithms for oligotrophic oceans: A novel approach based on three-band reflectance difference. J. Geophys. Res. Ocean. 2012, 117, 1–25. [Google Scholar] [CrossRef]
  86. Xue, K.; Zhang, Y.; Duan, H.; Ma, R.; Loiselle, S.; Zhang, M. A remote sensing approach to estimate vertical profile classes of phytoplankton in a Eutrophic lake. Remote Sens. 2015, 7, 14403–14427. [Google Scholar] [CrossRef]
  87. Fensholt, R.; Sandholt, I. Derivation of a shortwave infrared water stress index from MODIS near- and shortwave infrared data in a semiarid environment. Remote Sens. Environ. 2003, 87, 111–121. [Google Scholar] [CrossRef]
  88. Jordan, C.F. Derivation of Leaf-Area Index from Quality of Light on the Forest Floor. Ecology 1969, 50, 663–666. [Google Scholar] [CrossRef]
  89. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  90. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  91. Zarco-Tejada, P.J.; Berjón, A.; López-Lozano, R.; Miller, J.R.; Martín, P.; Cachorro, V.; González, M.R.; De Frutos, A. Assessing vineyard condition with hyperspectral indices: Leaf and canopy reflectance simulation in a row-structured discontinuous canopy. Remote Sens. Environ. 2005, 99, 271–287. [Google Scholar] [CrossRef]
  92. Freund, Y.; Schapire, R.E. A Decision-Theoretic Generalization of On-Line Learning and an Application to Boosting. J. Comput. Syst. Sci. 1997, 55, 119–139. [Google Scholar] [CrossRef] [Green Version]
  93. Friedl, M.A.; McIver, D.K.; Hodges, J.C.F.; Zhang, X.Y.; Muchoney, D.; Strahler, A.H.; Woodcock, C.E.; Gopal, S.; Schneider, A.; Cooper, A.; et al. Global land cover mapping from MODIS: Algorithms and early results. Remote Sens. Environ. 2002, 83, 287–302. [Google Scholar] [CrossRef]
  94. DeFries, R. Multiple Criteria for Evaluating Machine Learning Algorithms for Land Cover Classification from Satellite Data. Remote Sens. Environ. 2000, 74, 503–515. [Google Scholar] [CrossRef]
  95. Zhou, Z.; Huang, J.; Wang, J.; Zhang, K.; Kuang, Z.; Zhong, S.; Song, X. Object-oriented classification of sugarcane using time-series middle-resolution remote sensing data based on AdaBoost. PLoS ONE 2015, 10, e0142069. [Google Scholar] [CrossRef] [PubMed]
  96. Maxwell, A.E.; Warner, T.A.; Fang, F. Implementation of machine-learning classification in remote sensing: An applied review. Int. J. Remote Sens. 2018, 39, 2784–2817. [Google Scholar] [CrossRef]
  97. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices, 2nd ed.; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2009; p. 183. ISBN 978-1-4200-5512-2. [Google Scholar]
  98. Xu, L.; Yan, P.; Chang, T. Best first strategy for feature selection. In Proceedings of the 9th International Conference on Pattern Recognition, Rome, Italy, 14 May–17 November 1988; pp. 706–708. [Google Scholar]
  99. Cheng, Z.; Zu, Z.; Lu, J. Traffic crash evolution characteristic analysis and spatiotemporal hotspot identification of urban road intersections. Sustainability 2018, 11, 161. [Google Scholar] [CrossRef]
  100. Al-Ruzouq, R.; Hamad, K.; Abu Dabous, S.; Zeiada, W.; Khalil, M.A.; Voigt, T. Weighted Multi-attribute Framework to Identify Freeway Incident Hot Spots in a Spatiotemporal Context. Arab. J. Sci. Eng. 2019. [Google Scholar] [CrossRef]
Figure 1. Map of the Gulf; (a) location map and (b) Moderate resolution imaging spectroradiometer (MODIS) satellite image with true colour captured on 1 March 2003 showing milky patches.
Figure 1. Map of the Gulf; (a) location map and (b) Moderate resolution imaging spectroradiometer (MODIS) satellite image with true colour captured on 1 March 2003 showing milky patches.
Remotesensing 11 01193 g001
Figure 2. Flowchart of the overall procedure.
Figure 2. Flowchart of the overall procedure.
Remotesensing 11 01193 g002
Figure 3. Initiation and disappearance of a whiting event in 2003 based on acquired MODIS images.
Figure 3. Initiation and disappearance of a whiting event in 2003 based on acquired MODIS images.
Remotesensing 11 01193 g003
Figure 4. Spectral response of various surface features in a sample MODIS product (MOD09GA) for the Gulf.
Figure 4. Spectral response of various surface features in a sample MODIS product (MOD09GA) for the Gulf.
Remotesensing 11 01193 g004
Figure 5. Generic statistics of whiting occurrence: (a) total frequency per year, (b) total frequency per month, (c) total number of days per year and (d) total number of days per month.
Figure 5. Generic statistics of whiting occurrence: (a) total frequency per year, (b) total frequency per month, (c) total number of days per year and (d) total number of days per month.
Remotesensing 11 01193 g005
Figure 6. Performance of the combinations of multiple selected features.
Figure 6. Performance of the combinations of multiple selected features.
Remotesensing 11 01193 g006
Figure 7. AdaBoost DT-generated model for developing rule sets for rule-based classification.
Figure 7. AdaBoost DT-generated model for developing rule sets for rule-based classification.
Remotesensing 11 01193 g007
Figure 8. MODIS images and their classification results: (a,b) 2 March 2003, (c,d) 4 March 2012, (e,f) 10 February 2016 and (g,h) 3 October 2018.
Figure 8. MODIS images and their classification results: (a,b) 2 March 2003, (c,d) 4 March 2012, (e,f) 10 February 2016 and (g,h) 3 October 2018.
Remotesensing 11 01193 g008
Figure 9. High spatial frequency of whiting events overs 17 years.
Figure 9. High spatial frequency of whiting events overs 17 years.
Remotesensing 11 01193 g009
Table 1. Description of the examined features derived from MODIS images.
Table 1. Description of the examined features derived from MODIS images.
No.Examined Feature NameAbbreviationsDescriptionMODIS BandsRef.
1–7Mean values of an image object of MODIS reflectance (ref.) bandsRef. 1–7Mean of bands 1–7 (Red, NIR, Blue, Green, SWIR1, SWIR2 and SWIR3)B1–B7[80]
8–14Standard deviation of an image object of ref. bandsSD 1–7Standard deviations of individual bands 1–7B1–B7[80]
15Normalised difference vegetation index NDVI ( Ref . N I R Ref . R e d ) ( Ref . N I R + Ref . R e d ) B2, B1[81]
16Normalised difference water index NDWI ( Ref . G r e e n Ref . N I R ) ( Ref . G r e e n + Ref . N I R ) B4, B2[82]
17Modified normalised difference water indexMNDWI ( Ref . G r e e n Ref . S W I R 1 ) ( Ref . G r e e n + Ref . S W I R 1 ) B4, B7[83]
18Floating algae indexFAI Ref . NIR - [ Ref . Red +   ( Ref . SWIR 1 - Ref . Red )   ( 859 - 645 ) ( 1240 - 645 ) ] B1, B2, B5[84]
19Color index for estimating PICCI Ref . Green - [ Ref . blue +   ( Ref . Red - Ref . Blue )   ( 555 - 443 ) ( 670 - 443 ) ] B1, B3, B4[85]
20Color index using 547, 667 and 869 nm for estimating PICCI869 Ref . Red - [ Ref . Green +   ( Ref . NIR - Ref . Green )   ( 667 - 547 ) ( 869 - 547 ) ] B1, B2, B4[34]
21Color index using 547 and 667 nm for estimating PICCI2 Ref . Green   -   Ref . Red B1, B4[34]
22Normalised difference algal bloom indexNDBI ( Ref . G r e e n Ref . R e d ) ( Ref . G r e e n + Ref . R e d ) B4, B1[86]
23Shortwave infrared water stress indexSIWS ( Ref . S W I R 2 Ref . N I R ) ( Ref . S W I R 2 + Ref . N I R ) B6, B2[87]
24Ratio vegetation index 1RVI 1 Ref . N I R   Ref . R e d B2, B1[88]
25Ratio vegetation index 2RVI 2 Ref . R e d   Ref . N I R B1, B2[88]
26Enhanced vegetation indexEVI2.5 ( Ref . N I R Ref . R e d ) 1 + Ref . N I R +   6 Ref . R e d   7.5 Ref . B l u e ) B2, B1, B3[89]
27Ratio of the reflectance values of red and green bandsRatio RG ( Ref . R e d Ref . G r e e n ) ( Ref . R e d + Ref . G r e e n ) B1, B4[90]
28Blue/red index BRI Ref . B l u e   Ref . R e d B3, B1[91]
29Blue/green indexBGI Ref . B l u e   Ref . G r e e n B3, B4[91]
30Normalised difference between green and red bands NDGR ( Ref . G r e e n Ref . R e d ) ( Ref . G r e e n + Ref . R e d ) B4, B1-
31Normalised difference between green and blue bands NDGB ( Ref . G r e e n Ref . B l u e ) ( Ref . G r e e n + Ref . B l u e ) B4, B3-
32Normalised difference between blue and green bands NDBG ( Ref . B l u e Ref . G r e e n ) ( Ref . B l u e + Ref . G r e e n ) B3, B4-
Table 2. Dates and number of image object samples selected from various images.
Table 2. Dates and number of image object samples selected from various images.
DatesWhiting SamplesClear WaterOther Segments
28 February 20037414974
2 March 2003100125126
3 March 20031213674
26 February 20046912466
2 February 20183616072
Sum400594412
Table 3. Statistics of occurrence of whiting in the Gulf for 16 years (July 2002 to July 2018).
Table 3. Statistics of occurrence of whiting in the Gulf for 16 years (July 2002 to July 2018).
YearMonthDatesPeriod (d)Frequency EventYearMonthDatesPeriod (d)Frequency Event
2002February1–551 November10–1122
December23–3081November24–252
2003February28–12-December1–662
March1–331December27–304
December9–11312012January15–1622
2004January29–3131January21–266
February7–1043February3–752
February16–172February21–244
February25–273March4–13102
March21–2331March18–192
November24–2961November2–322
2005February10–1121November12–143
November1–111122013January11–1881
November24–318February4–522
2006January15–2061February14–152
December9–1022March12–1321
December25–262December11–1882
2007January1–441December20–278
March4–5222014February11–1332
March11–133February19–213
April19–2131November7–1042
November25–3061November26–305
December11–1332December2–762
December24–252December25–262
2008February3–7522015January19–2461
February21–277February27–2821
March6–721November13–1531
December17–19312016Jan4–631
2009January4–742Jan–Feb29–241
January14–163February10–1341
February3–6412017February4–851
November7–821November10–1452
2010January27–2821November29–302
November23–2861December4–521
December16–17212018January2–541
2011January12–1431January29–3131
February4–522February1–441
Table 4. Optimum scale parameter values selected for a single scale using the objective function (OF) and the F-measure methods (Equations (1)–(2).
Table 4. Optimum scale parameter values selected for a single scale using the objective function (OF) and the F-measure methods (Equations (1)–(2).
SPNo. of ObjectsWeighted VarianceMoran’s IWV NormMI NormOFF-Measure
1032,227326.85100.41181010
2012,4461012.34730.25970.88830.37471.26300.5271
3066951784.30100.20130.76250.51871.28120.6174
4042082606.25370.14250.62850.66361.29210.6456
5028913394.25220.10220.50010.76301.26310.6042
6021974173.26230.07890.37320.82041.19360.5130
7017554780.72610.04850.27420.89531.16950.4198
8014595468.11060.03160.16210.93701.09910.2764
9012026039.13180.01780.06910.97081.03990.1290
10010396463.04620.00600110
Table 5. Accuracy assessment of various tree-based classifiers.
Table 5. Accuracy assessment of various tree-based classifiers.
ClassifierOAKC
AdaBoost DT97.86%0.97
Gradient boosted DT97.12%0.96
Single DT96.19%0.95
Random forest95.00%0.93
Table 6. Cloud-free MODIS images that show peak whiting patches between 2002 and 2018 and the computed whiting area coverage.
Table 6. Cloud-free MODIS images that show peak whiting patches between 2002 and 2018 and the computed whiting area coverage.
YearDateWhiting Area (km2)Percentage of Gulf Area
20025 February15,6556.55
20032 March53,68722.46
20049 February29,87412.50
20056 April44,89418.78
200630 January30,54912.78
200711 December20,4818.57
200822 February47,88720.04
20095 February22,5849.45
201016 December12,1005.06
20115 December17,3407.26
20124 March60,84725.46
20134 February39,54416.55
201410 November15,1376.33
201515 November19,2018.03
201610 February30,48012.75
20176 January45,75319.14
20183 February22,1599.27

Share and Cite

MDPI and ACS Style

Shanableh, A.; Al-Ruzouq, R.; Gibril, M.B.A.; Flesia, C.; AL-Mansoori, S. Spatiotemporal Mapping and Monitoring of Whiting in the Semi-Enclosed Gulf Using Moderate Resolution Imaging Spectroradiometer (MODIS) Time Series Images and a Generic Ensemble Tree-Based Model. Remote Sens. 2019, 11, 1193. https://doi.org/10.3390/rs11101193

AMA Style

Shanableh A, Al-Ruzouq R, Gibril MBA, Flesia C, AL-Mansoori S. Spatiotemporal Mapping and Monitoring of Whiting in the Semi-Enclosed Gulf Using Moderate Resolution Imaging Spectroradiometer (MODIS) Time Series Images and a Generic Ensemble Tree-Based Model. Remote Sensing. 2019; 11(10):1193. https://doi.org/10.3390/rs11101193

Chicago/Turabian Style

Shanableh, Abdallah, Rami Al-Ruzouq, Mohamed Barakat A. Gibril, Cristina Flesia, and Saeed AL-Mansoori. 2019. "Spatiotemporal Mapping and Monitoring of Whiting in the Semi-Enclosed Gulf Using Moderate Resolution Imaging Spectroradiometer (MODIS) Time Series Images and a Generic Ensemble Tree-Based Model" Remote Sensing 11, no. 10: 1193. https://doi.org/10.3390/rs11101193

APA Style

Shanableh, A., Al-Ruzouq, R., Gibril, M. B. A., Flesia, C., & AL-Mansoori, S. (2019). Spatiotemporal Mapping and Monitoring of Whiting in the Semi-Enclosed Gulf Using Moderate Resolution Imaging Spectroradiometer (MODIS) Time Series Images and a Generic Ensemble Tree-Based Model. Remote Sensing, 11(10), 1193. https://doi.org/10.3390/rs11101193

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