Next Article in Journal
Effects of Microtopography on Runoff Generation in Plain Farmland: New Insights from an Event-Based Rainfall–Runoff Model
Next Article in Special Issue
Identifying Alpine Lakes in the Eastern Himalayas Using Deep Learning
Previous Article in Journal
Role of Reef-Flat Plate on the Hydrogeology of an Atoll Island: Example of Rangiroa
Previous Article in Special Issue
Cyclicities in the Regime of Groundwater and of Meteorological Factors in the Basin of the Southern Bug River
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Different Water Indices and Band Downscaling Methods for Water Bodies Mapping from Sentinel-2 Imagery at 10-M Resolution

1
Guangdong Technology Center of Hydraulic and Hydropower Project, Guangzhou 510635, China
2
Research Centre for Oceanic Remote Sensing Big Data Applications, Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China
3
Guangdong Open Laboratory of Geospatial Information Technology and Application, Key Laboratory of Guangdong for Utilization of Remote Sensing and Geographical Information System, Guangzhou Institute of Geography, Guangdong Academy of Sciences, Guangzhou 510070, China
*
Authors to whom correspondence should be addressed.
Water 2022, 14(17), 2696; https://doi.org/10.3390/w14172696
Submission received: 4 May 2022 / Revised: 13 August 2022 / Accepted: 24 August 2022 / Published: 30 August 2022
(This article belongs to the Special Issue Inland Surface Water and Deep Learning)

Abstract

:
Satellite-based remote sensing is important for monitoring the spatial distribution of water resources. The water index is currently one of the most widely used water body extraction methods. Based on Sentinel-2 remote sensing image, this study combines area-to-point regression kriging interpolation, bilinear interpolation, and the Gram–Schmidt (GS) pan-sharpening method with the water indices MNDWI, AWEIsh and WI2015 to compare different water body extraction methods. The experimental results showed that all water indices have satisfactory extraction ability, with the kappa coefficient as an accuracy threshold above 0.8. Moreover, the GS downscaling method combined with the WI2015 yielded the best performance. This research demonstrates the efficacy of the WI2015 method to extract water bodies in urban areas and its ability to comprehensively describe river water bodies. The findings indicate that high-resolution band information is particularly important for improving low-resolution band downscaling results and can significantly minimize erroneous water body extraction.

1. Introduction

Urban surface water bodies are important factors influencing the urban ecological environment and exert a certain impact on urban public health and people′s quality of life [1,2]. At the same time, urban water bodies play a key role in urban planning, regional climate change, the heat island effect, and water resource utilization [3,4]. In recent years, with adverse factors such as rapid urbanization, environmental degradation, and extreme climate, the water body area in urban cities has been decreasing significantly [5]. Therefore, the accurate and dynamic monitoring of urban surface water bodies has become essential for water resource management and decision-making. As a large-scale and real-time Earth observation technology, remote sensing has been widely used in surface feature recognition and extraction. Accordingly, remote sensing provides feasible technical means for the automatic and precise extraction of urban surface water bodies [6,7]. The growing role of satellite remote sensing technology in water extraction applications is becoming increasingly remarkable. To identify water body features, remote sensing images mainly use spectral differences between water bodies and other ground objects in different wavelength bands. The development of water body remote sensing methods has progressed through several stages, from the initial manual visual interpretation technique to semi-automatic extraction and classification techniques based on spectral features, and then further to extraction methods which couple spectral features and spatial information [8,9,10,11]. Currently, automatic high-precision water body extraction methods based on deep learning represent the gold standard [12,13,14,15].
Water index-based algorithms have become important for implementing rapid water body mapping in large-scale regions. Water index- and threshold-based approaches have been widely used to identify water bodies due to their unique spectral characteristics in the visible and infrared regions. Both approaches have undergone significant evolution. In 1996, McFeeters [16] proposed the normalized difference water index (NDWI) using the value of the green band minus the near-infrared (NIR) band, divided by the sum of the two bands. Under this, water bodies have positive values, while non-water bodies have negative values. Although NDWI can suppress and remove non-water features to a large degree, it fails to efficiently suppress built-up land signals. Consequently, certain features may actually comprise a mixture of water and built-up land noise. In 2006, based on the NDWI, in 2006, Xu proposed a modified normalized difference water index (MNDWI), replacing the NIR band with the shortwave-infrared (SWIR) band, which helped to remove disturbances caused by built-up lands [17]. However, the optimal thresholds varied based on location and time, and the method could not effectively remove shadow noise in some areas. In 2014, Feyisa et al. proposed an automated water extraction index (AWEI) and used different AWEI formats for scenes with shadows (AWEIsh) and without shadows (AWEInsh) [18]. This technique was sufficiently separated and systematic so as to improve the accuracy of water body mapping. Related scholars have used the natural logarithm of each band of the Landsat 7 ETM + image as a proxy of the reflection coefficient and interaction conditions, creating the water index WI2006. Subsequently, the water index WI2015 is based on WI2006, using linear discriminant analysis classification (LDAC) to determine the coefficient of the best segmentation training area category, further improving water extraction accuracy [19].
Many researchers have used moderate resolution imaging spectroradiometer (MODIS) [20], Landsat [21,22], and Sentinel [23,24] multispectral remote sensing images to achieve large-scale water body extraction research based on various water indices, such as MNDWI. In recent years, high-resolution remote sensing technology has achieved significant development. However, these fine spatial solution images, such as Gaofen-2 Satellite (GF-2) and Satellite Pour l′Observation de la Terre (SPOT), have no SWIR band, making it impossible to use the water index method [25]. In contrast, Sentinel-2 provides publicly available images. The Sentinel-2 mission has been organized by the Global Monitoring for Environment and Security. Using a bi-satellite system, it acquires multispectral, high-resolution optical observations over global terrestrial surfaces with a high revisitation frequency, approximately five days. Such a system is important for dynamic land cover mapping and updating. Sentinel-2 carries a multispectral instrument with 13 spectral bands spanning the visible spectrum (VIS) and NIR to SWIR. These spatial resolutions range from 10 to 60 m with a 290 km field of view on the ground [26]. With high-frequency and high-spectral resolution imaging, Sentinel-2 allows intensively and continuously monitoring of the Earth’s surface. Sentinel-2 multispectral instrument imagery includes 20 m resolution SWIR bands and 10 m resolution green and NIR bands, rendering water mapping based on water indices at 10 m resolution possible.
A useful way to improve the performance of water body mapping using Sentinel-2 imagery is to produce water indices results by downscaling the SWIR bands from 20 to 10 m. Trivially, the key challenge lies in accurately increasing the spatial resolution of the SWIR band. Spatial interpolation (such as bilinear interpolation) and image fusion (such as pan-sharpening) are the two most widely used methods to increase the spatial resolution of remote sensing imagery [27]. The spatial interpolation method is directly applied to coarse spatial resolution images without requiring additional datasets. In contrast, image fusion, such as pan-sharpening, is premised on the availability of the fine spatial resolution panchromatic (PAN) band of the same scene, aiming to downscale coarse multispectral imagery to the spatial resolution of the PAN band. Pan-sharpening is widely applied to remote sensing images with coarse multispectral bands and a fine spatial resolution PAN band [28].
Nevertheless, most previous studies primarily evaluated the effect of band downscaling methods or water indices on remote sensing water body extraction. However, combinations of typical band downscaling methods and frequently-used water indices have not been synthetically compared and analyzed. From this perspective, this study aims to compare the results of remote sensing water body extraction based on combinations of different band downscaling methods and water indices. Specifically, the effects of three types of factors on water body extraction were evaluated: (1) the extraction capability of the water indices MNDWI, AWEIsh, and WI2015; (2) the effect of the SWIR band downscaling results based on area-to-point regression kriging interpolation (ATPRK), bilinear interpolation (BIL), and the Gram–Schmidt (GS) pan-sharpening method; and (3) the segment precision between water/non-water bodies based on the marker-controlled watershed (MCW) algorithm.

2. Materials and Methods

2.1. Study Area

This study is a preliminary exploration to evaluate combinations of several kinds of downscaling methods and water indices in urban cities. Therefore, a typical urban area in Guangzhou city has been chosen for this comparative experiment. The study area covers the northeastern part of Haizhu District and the northernmost area of Panyu District, Guangzhou, covering a total area of more than 100 km2 (Figure 1). This urban area mainly includes large water bodies, such as the Haizhu Wetland National Park, Huangpuyong, and Guanzhou waterways, as well as many small water bodies. The densely distributed buildings and their shadows in urban areas tend to significantly interfere with the remote sensing extraction of water bodies, causing poor accuracy.

2.2. Methods

Figure 2 shows the workflow for the remote sensing extraction of water bodies. First, the spatial resolutions of the Sentinel-2 SWIR bands were improved to 10 m through BIL, ATPRK interpolation [29], and GS pan-sharpening. Then, the MNDWI, AWEIsh, and WI2015 water indices were calculated using the 10 m VIS, NIR, and SWIR bands. Next, masker threshold training was performed for each water index image, and the marker-controlled watershed algorithm was used to realize the segmentation and extraction of water bodies. Finally, waterbody reference data were used to evaluate and compare different combinations of band downscaling methods and water indices in terms of extraction accuracy.

2.2.1. Water Indices

NDWI is the first water index proposed for water body remote sensing extraction, based on the combination of the green band ρ green and the NIR band ρ NIR of remote sensing images, and is calculated as follows:
NDWI = ρ green ρ NIR ρ green + ρ NIR
NDWI mainly takes advantage of the strong absorption of water in the NIR band and the absence of strong reflectivity from vegetation [16]. It is used to mine water information from an image by suppressing vegetation and highlighting the water body, enhancing the clarity of the resulting images. However, NDWI considers only vegetation factors, ignoring the two key features of buildings and soil. When extracting water information through NDWI, the reflectivity of the green band is much higher than that of the NIR band. As such, the extraction results are often confused with soil and building information. When NDWI was used to extract urban water bodies, there were water bodies with more shadows of buildings, and the effect was poor.
Based on NDWI, the MNDWI modified the band combination of the water index and replaced the NIR band in NDWI with the SWIR band ρ SWIR . The calculation formula is as follows:
MNDWI = ρ Green ρ SWIR ρ Green + ρ SWIR
The spectral characteristics of the building shadows in the green and NIR bands were similar to those of water. Using the SWIR band to replace the NIR band, the contrast between water bodies and building shadows was significantly enhanced, which greatly improved the contrast between these two ground features, ultimately promoting accurate extraction of water body information in cities and towns. Xu [17] conducted experiments using remote sensing images containing different types of water bodies. The analysis revealed that, relative to NDWI, MNDWI could better extract fine features of water bodies, including the distribution of suspended sediments and water quality. Feyisa et al. [18] conducted experiments on Landsat TM images and proposed AWEI for factors such as low classification accuracy and relatively unfixed threshold selection in previous water body information extraction. AWEInsh is suitable for settings without shadows, whereas AWEIsh is designed to eliminate shadows and other ground objects easily confused with water information in the AWEInsh extraction results. Therefore, AWEIsh is suitable for scenes with more shadows; its formula is as follows:
AWEI sh = ρ blue + 2.5 ρ Green 1.5 ρ NIR + ρ SWIR 1 0.25 ρ SWIR 2
WI2006 is a water index created using standard variables to analyze the emissivity of the atmospheric surface. The natural logarithm of each band of the Landsat 7 ETM+ images was used to determine the reflection coefficient and interaction conditions. It has been applied to the extraction research of eastern Australian wetlands. In 2015, based on WI2006, Fisher et al. [19] created a new water index, WI2015, which uses LDAC as a coefficient to determine the best segmentation training area category. The calculation formula is as follows:
WI 2015 = 1.7204 + 171 ρ Green + 3 ρ Red 70 ρ NIR 45 ρ SWIR 1 71 ρ SWIR 2

2.2.2. Gram–Schmidt Pan-Sharpening

Spatial interpolation and pan-sharpening algorithms are typically used to achieve spatial downscaling of low-resolution bands. The most widely used pan-sharpening algorithms include the principal component analysis (PCA) [29], hue–saturation–value (HSV) [30], high pass filter (HPF) [31], and GS techniques [32]. Specifically, GS exhibits the highest spectral fidelity and can maintain the consistency of the band spectral characteristics before and after pan-sharpening; that is, the high-resolution band data obtained by downscaling retains the spectral characteristics of the original low-resolution band. Therefore, this study uses the GS pan-sharpening algorithm to preserve the original Sentinel-2 SWIR spectral information as much as possible.
There is a lack of panchromatic band information in Sentinel-2 images, whereas GS pan-sharpening requires the incorporation of high-resolution band information, similar to the PAN band. The resolutions of all four VIS/NIR multispectral bands of the Sentinel-2 image have been determined to be 10 m. However, the correlation between the bands tends to lead to data redundancy. The classical PCA method is widely used for the dimensionality reduction of multispectral remote sensing image band information since PCA can compress VIS/NIR multispectral bands into a panchromatic-like band. By linearly transforming the four VIS/NIR bands, mutually orthogonal spectral spaces are generated, in which the first principal component (FPC) contains the most abundant information. Therefore, the FPC can be regarded as a 10 m panchromatic-like band. Through fusion with the FPC using the GS method, the spatial resolution of the SWIR band can be increased to 10 m.

2.2.3. Area-to-Point Regression Kriging (ATPRK)

ATPRK was used to perform remote sensing image band fusion. ATPRK combines the traditional regression kriging interpolation and quantitative remote sensing scale conversion theory. First, one must assume that the band reflectivity Z l x i is a random variable of the grid points x i   i = 1 , , M in the low-resolution band l   l = 1 , , L , M is the number of grid points. Z k v j is the random variable of the grid points v j   j = 1 , , MF 2 in the high-resolution band k   k = 1 , , K , and F is the ratio of the high- and low-resolution values. According to the regression kriging theory, the spatial downscaling results of the band l achieved by ATPRK in the high-resolution grid v are realized by the estimation and addition of the trend term and residual term m ^ l v and r ^ l v , respectively), as shown in the following formula:
Z ^ l v = m ^ l v + r ^ l v
At a specific high-resolution grid point v 0 , the estimated value m ^ l v 0 of the trend item was obtained by linear regression of the value Z k v 0 at the high-resolution band k :
m ^ l v 0 = k = 0 K a k l Z k v 0 ,   Z 0 v 0 = 1
According to the assumption of scale invariance, the regression model above is consistent with the regression model established by the value Z l x at the low-resolution band l and the value Z k x at the upscaling band k :
Z l x = k = 0 K a k l Z k x + r l x ,   Z 0 x = 1 ,   x
where r l x is the regression residual term at the corresponding band l , and the regression coefficient a k l is estimated using least squares.
After performing regression analysis on the trend item, ATPRK interpolation can be used to implement spatial downscaling of the residual item. The ATPRK interpolation downscaling results maintain the original spectral band information. The estimated value r ^ l v 0 of the residual term at the corresponding high-resolution band l , is the linear weighted average of the residual term r l x i of the low-resolution adjacent grid point:
r ^ l v 0 = i = 1 N λ i r l x i ,   s . t .   i = 1 N λ i = 1
where N is the number of adjacent grid points, and λ i is the corresponding weight value calculated by the following Kriging equations:
γ cc l x 1 , x 1 γ cc l x 1 , x N 1 γ cc l x N , x 1 1 γ cc l x N , x N 1 1 0 λ 1 λ N μ = γ fc l v 0 , x 1 γ fc l v 0 , x N 1
where γ cc l x i , x j is the area-to-area variogram between the low-resolution grid points on the band, γ fc l v 0 , x j is the area-to-point variogram between the high-resolution grid points to be estimated and the low-resolution neighboring grid points, and μ is the Lagrangian operator. s is assumed to be the distance between the centers of any two grid points. The variograms γ cc l s and γ fc l s can be calculated by a convolution between the point-to-point variogram γ ff l s and the point spread function h l s (* is the convolution operator) [33]:
γ fc l s = γ ff l s     h l s γ cc l s = γ ff l s     h l s     h l s
It is essential to determine the point-to-point variogram, calculate the low-resolution residual variogram by fitting, and then perform deconvolution inference. Notably, the point spread function selected in this project is a simple arithmetic average operation. Thus, the area-to-area and area-to-point variograms are converted into the mean values of multiple point-to-point variograms for calculation.
The estimated values of the trend item and residual item of the band at the high-resolution grid point were calculated, and the sum of these two was the final spatial downscaling result. The above calculation process was carried out for the low-resolution bands individually. Finally, the resolution of all bands was unified through the fusion of high-resolution and low-resolution bands.

2.2.4. Marker-Controlled Watershed Segmentation

Previous studies have shown that the MCW algorithm is particularly suitable for waterbody segmentation [34]. Compared with algorithms that use a single threshold to segment water/non-water bodies, such as maximum between-class variance, it performs better at the edges of water bodies. The typical process of using the MCW algorithm to conduct water body extraction includes three steps: (1) marking the water body/non-water body area for each water index image (including MNDWI, AWEIsh, and WI2015), marking the water body/non-water body area with high reliability; (2) gradient image generation: by applying the Sobel operator to each water index image, calculate and generate the corresponding gradient images, which are then used to determine the boundary and markers between water and non-water bodies; (3) performing water body segmentation based on water body/non-water area markers and gradient images. Here, the watershed algorithm iteratively expands each marker until all unmarked pixels are marked as water or non-water.
The last two steps of the MCW algorithm are relatively fixed, with no parameter involved. However, the water body segmentation result is more sensitive to the water/non-water area marking in the first step. Therefore, it was necessary to calibrate the marking selection parameters. In this study, the threshold method was used to automatically generate water/non-water markings, by combining the real distribution data of the water bodies and a pair of mask thresholds determined for each water index. The left range of the smaller threshold and the right range of the larger oner corresponds to the non-water and water areas, respectively. Notably, the transition range between thresholds is divided into water and non-water bodies by the watershed algorithm.

2.2.5. Accuracy Evaluation Indicators

Taking the waterbody reference data of high-resolution remote sensing images of the study area as the standard, the accuracy of the water body results extracted from Sentinel-2 remote sensing images were evaluated. The following four accuracy evaluation indicators were selected: the producer accuracy (PA), user accuracy (UA), overall accuracy (OA), and kappa coefficient. The calculation methods were as follows:
PA = TP TP + FN UA = TP TP + FP OA = TP + TN T Kappa = T × TP + TN T × T
where TP is the number of water pixels that are correctly extracted, FN is the number of water pixels that have not been extracted, FP is the number of water pixels that are incorrectly extracted, TN is the number of non-water pixels that are correctly extracted, T is the total number of image pixels, and = TP + FP × TP + FN + FN + TN × FP + TN .

3. Results

3.1. Band Downscaling Quality

An important criterion for evaluating the quality of a fusion image is its ability to maintain spectral characteristics (i.e., quality preservation). The 10 m band downscaling results generated by the three methods were upscaled to 20 m and then compared with the original 20 m band. Figure 3 shows the scatter plots comparing the upscaled results and the original coarse SWIR bands for the ATPRK, BIL, and GS methods. Notably, the spectral characteristics of each method remain unchanged over different coarse bands. Comparing the quality assurance of different methods, the GS band obtained by upscaling showed the least correlation with the original band. Compared with the GS method, the results of the BIL method had a stronger correlation with the original data but evidently underestimated the high-value area. A significant advantage of the ATPRK method for band fusion was its quality preservation. The comparison results show that ATPRK achieves non-destructive preservation of the original band spectrum information.

3.2. Water Indices Results

Figure 4 shows the three types of water index results obtained using the BIL, ATPRK, and GS downscaling methods. It is observed that by using the three water indices, the contrast between water bodies and land areas are better highlighted, and the boundaries of the water bodies are sufficiently clear. In the MNDWI image, the contrast between water and land is particularly strong. Conversely, the range of water index values calculated using the SWIR band by downscaling using the ATPRK method was wider than that of the GS and BIL methods. A characteristic of kriging interpolation compared with other interpolation methods was that its interpolation result might exceed the original data range.
Figure 5 shows the histograms of MNDWI, AWEIsh, and WI2015 water indices using ATPRK, BIL, and GS downscaling methods. All histograms have bimodal shapes, with threshold ranges between water and non-water bodies located at the bottom. The histogram results reveal that the technique of using single threshold values to exactly segment water/non-water bodies belongs to the category of theoretical cases for water index-based extraction methods. Comparatively, the AWEIsh results derived from the ATPRK, BIL, and GS band downscaling methods all had narrow threshold ranges. Furthermore, for all three water indices, the results derived from the ATPRK and BIL interpolation methods yielded similar numerical ranges, which were wider than those from the GS method.

3.3. Analysis of Water Body Commissions and Omissions

The MCW segmentation algorithm was used to segment each water index image to accurately extract the water body area. By superimposing and comparing the water extraction results with the reference data, misclassified and unidentified water body areas were obtained, as shown in Figure 6 and Figure 7, respectively.
In urban areas, the features most easily misclassified as water bodies are building shadows. Here, we focus on two high building density regions within the study area, which showed obvious shadows on the original image, characterized by similarly low-reflectivity features as water bodies. It can be observed that the water bodies extracted from the water index images were affected by the misclassification of building shadows to a significant extent (e.g., the misclassification in the yellow box). The MNDWI, AWEIsh, and WI2015 indices corresponding to the BIL method easily misclassified building shadows as water bodies. This was due to the fact that although the BIL interpolation method maintains the original spectral information of the low-resolution NIR band to a certain extent, it does not introduce other high-resolution band information. Therefore, the effect of improving the spatial details is limited. In comparison, the water body extraction results of the MNDWI_ATPRK and WI2015_GS combinations performed better; only a few building shadows were mistakenly identified as water bodies. The results show that the integration of high-resolution band information was particularly important for improving downscaling results of low-resolution bands, which significantly reduced the probability of water body identification. On the other hand, it is difficult to achieve the same effect by relying only on simple spatial interpolation methods (such as the BIL interpolation).
Many small water bodies in the study area were missed or poorly estimated by the three water indices since the ATPRK interpolation, BIL interpolation, and GS pan-sharpening methods unify the spatial resolution of all bands involved in the calculation of the water index to 10 m. Therefore, it is naturally impossible to extract sub-pixel-level small water bodies with a width of less than 10 m. Specifically, the MNDWI_GS and AWEIsh_GS combinations yielded significant omission errors of a small stream shown in the red box of Figure 7, indicating that the GS pan-sharpening method may be unsuitable for the extraction of small water bodies.

3.4. Quantitative Evaluation of Water Body Extraction Accuracy

In addition to the qualitative and intuitive expression of the results of inaccurate extractions and omissions of water bodies, a variety of precision indicators such as UA were used to quantitatively evaluate the water body extraction results of different downscaling methods and water index combinations. Table 1 shows the accuracy verification indices of water body extraction using the three water indices generated by different downscaling methods. Overall, each combination achieved good water extraction results; several accuracy evaluation indicators, such as UA, were above 80%. The MNDWI against ATPRK had the best UA accuracy (95.02%). The WI2015 for GS provided the best PA accuracy (89.10%), OA accuracy (96.79%), and kappa coefficient (0.897). Overall, the WI2015 calculated using GS downscaling band information yielded the best water extraction effect.
The ATPRK method is designed to preserve the spectral information of observed images. However, the water mapping results of AWEIsh and WI2015 from ATPRK were poorer than those from GS. The water mapping results are related to the ability of downscaling methods to maintain spectral information. Still, they are also affected by other properties, such as the preservation of spatial details. Similar results were found in a previous study [5]. Although the HPF downscaling method can better preserve the spectral information of the original image, it cannot produce water body maps with higher accuracy than the other methods.
The extraction effects of the three water indices in the study area were reflected by the average accuracy indices of the water body extraction results corresponding to different downscaling methods. MNDWI displayed the highest average UA accuracy (94.23%), implying the least commission errors. WI2015 has the highest average PA value (88.54%), meaning the lowest omission errors. Likewise, the OA and kappa coefficients returned the best results for water bodies extracted by WI2015. Therefore, in this study, it can be concluded that WI2015 performed best on the urban water bodies mapping, while MNDWI and AWEIsh provided similar performances.

4. Conclusions

Based on Sentinel-2 remote sensing images, this study used BIL interpolation, ATPRK interpolation, and panchromatic sharpening GS spatial downscaling methods to increase the spatial resolution of the SWIR bands to 10 m and calculate the corresponding water indices MNDWI, AWEIsh, and WI2015. The MCW water segmentation algorithm was used to segment the water/non-water body area on each water index image and qualitatively analyze the water body extraction, inaccurate extractions, and omission results. Finally, combined with UA, PA, OA, and the kappa coefficient, a quantitative evaluation of the water body extraction results by combining the different spatial downscaling methods with water indices was conducted.
Our results indicate that the water index based on Sentinel-2 remote sensing images can effectively extract water body information in urban areas, especially river water bodies. The water body extraction accuracy of different water indices remained above 0.8. The combination of the GS spatial downscaling method and the WI2015 water index yielded the best water body extraction efficiency, with a kappa coefficient of 0.897. The effectiveness and feasibility of using satellite remote sensing technology to monitor the distribution of water bodies were verified.
Limited by the spatial resolution of Sentinel-2 images, the extraction results of each water index did not have enough resolution to provide information on sub-pixel-level, small water bodies with a width less than 10 m. Extracting small water bodies has long been a challenge in remote sensing water body extraction research. In follow-up research, we will implement high-resolution remote sensing images (such as GF-2 images) to comparatively assess the extraction capabilities of various water indices for small water bodies in urban areas.

Author Contributions

Conceptualization, H.H. and X.L.; methodology, H.L.; validation, H.J. and W.L.; model construction, H.L. and H.H.; data curation, W.L. and X.Y.; writing—original draft preparation, H.L.; writing—review and editing, H.H. and X.L.; visualization, H.J. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly supported by the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (No. GML2019ZD0301), the Natural Science Foundation of Guangdong Province (No. 2020A1515010643 and 2020A1515011068), the GDAS′ Special Project of Science and Technology Development (No. 2020GDASYL-20200104004), 2019 Science and Technology Innovation Project of Provincial Special Fund for Economic Development (No. 2019B1), Guangdong Province Agricultural Science and Technology Innovation and Promotion Project (No. 2022KJ102) and the Guangdong Innovative and Entrepreneurial Research Team Program (No. 2016ZT06D336).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the Geographical Science Data Center of the Greater Bay Area and the Guangdong Provincial Engineering Laboratory of Geographic Spatio-temporal Big Data for the data support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xie, H.; Luo, X.; Xu, X.; Pan, H.; Tong, X. Automated Subpixel Surface Water Mapping from Heterogeneous Urban Environments Using Landsat 8 OLI Imagery. Remote Sens. 2016, 8, 584. [Google Scholar] [CrossRef]
  2. Zhou, Y.A.; Luo, J.C.; Shen, Z.F.; Hu, X.D.; Yang, H.P. Multiscale Water Body Extraction in Urban Environments From Satellite Images. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4301–4312. [Google Scholar] [CrossRef]
  3. Huang, X.; Xie, C.; Fang, X.; Zhang, L.P. Combining Pixel- and Object-Based Machine Learning for Identification of Water-Body Types From Urban High-Resolution Remote-Sensing Imagery. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 2097–2110. [Google Scholar] [CrossRef]
  4. Sun, X.; Tan, X.Y.; Chen, K.L.; Song, S.; Zhu, X.D.; Hou, D.L. Quantifying landscape-metrics impacts on urban green-spaces and water-bodies cooling effect: The study of Nanjing, China. Urban For. Urban Green. 2020, 55, 126838. [Google Scholar] [CrossRef]
  5. Du, N.; Ottens, H.; Sliuzas, R. Spatial impact of urban expansion on surface water bodies-A case study of Wuhan, China. Landsc. Urban Plann. 2010, 94, 175–185. [Google Scholar] [CrossRef]
  6. Pekel, J.-F.; Cottam, A.; Gorelick, N.; Belward, A.S. High-resolution mapping of global surface water and its long-term changes. Nature 2016, 540, 418–422. [Google Scholar] [CrossRef] [PubMed]
  7. Lehner, B.; Doll, P. Development and validation of a global database of lakes, reservoirs and wetlands. J. Hydrol. 2004, 296, 1–22. [Google Scholar] [CrossRef]
  8. Qi, B.; Zhuang, Y.; Chen, H.; Dong, S.; Li, L. Fusion feature multi-scale pooling for water body extraction from optical panchromatic images. Remote Sens. 2019, 11, 245. [Google Scholar] [CrossRef]
  9. Jiang, W.; Ni, Y.; Pang, Z.; Li, X.; Ju, H.; He, G.; Lv, J.; Yang, K.; Fu, J.; Qin, X. An effective water body extraction method with new water index for sentinel-2 imagery. Water 2021, 13, 1647. [Google Scholar] [CrossRef]
  10. Yue, H.; Li, Y.; Qian, J.; Liu, Y. A new accuracy evaluation method for water body extraction. Int. J. Remote Sens. 2020, 41, 7311–7342. [Google Scholar] [CrossRef]
  11. Chen, Y.; Tang, L.; Kan, Z.; Bilal, M.; Li, Q. A novel water body extraction neural network (WBE-NN) for optical high-resolution multispectral imagery. J. Hydrol. 2020, 588, 125092. [Google Scholar] [CrossRef]
  12. Li, L.; Yan, Z.; Shen, Q.; Cheng, G.; Gao, L.; Zhang, B. Water body extraction from very high spatial resolution remote sensing data based on fully convolutional networks. Remote Sens. 2019, 11, 1162. [Google Scholar] [CrossRef]
  13. Li, M.; Wu, P.; Wang, B.; Park, H.; Yang, H.; Wu, Y. A deep learning method of water body extraction from high resolution remote sensing images with multisensors. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 3120–3132. [Google Scholar] [CrossRef]
  14. Yu, Y.; Yao, Y.; Guan, H.; Li, D.; Liu, Z.; Wang, L.; Yu, C.; Xiao, S.; Wang, W.; Chang, L. A self-attention capsule feature pyramid network for water body extraction from remote sensing imagery. Int. J. Remote Sens. 2021, 42, 1801–1822. [Google Scholar] [CrossRef]
  15. Shi, W.; Sui, H. An effective superpixel-based graph convolutional network for small waterbody extraction from remotely sensed imagery. Int. J. Appl. Earth Obs. Geoinf. 2022, 109, 102777. [Google Scholar] [CrossRef]
  16. 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]
  17. Xu, H.Q. 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]
  18. Feyisa, G.L.; Meilby, H.; Fensholt, R.; Proud, S.R. Automated Water Extraction Index: A new technique for surface water mapping using Landsat imagery. Remote Sens. Environ. 2013, 140, 23–35. [Google Scholar] [CrossRef]
  19. Fisher, A.; Flood, N.; Danaher, T. Comparing Landsat water index methods for automated water classification in eastern Australia. Remote Sens. Environ. 2016, 175, 167–182. [Google Scholar] [CrossRef]
  20. Sharma, R.C.; Tateishi, R.; Hara, K.; Luong Viet, N. Developing Superfine Water Index (SWI) for Global Water Cover Mapping Using MODIS Data. Remote Sens. 2015, 7, 13807–13841. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, X.B.; Xie, S.P.; Zhang, X.L.; Chen, C.; Guo, H.; Du, J.K.; Duan, Z. A robust Multi-Band Water Index (MBWI) for automated extraction of surface water from Landsat 8 OLI imagery. Int. J. Appl. Earth Obs. Geoinf. 2018, 68, 73–91. [Google Scholar] [CrossRef]
  22. Liu, X.P.; Hu, G.H.; Chen, Y.M.; Li, X.; Xu, X.C.; Li, S.Y.; Pei, F.S.; Wang, S.J. High-resolution multi-temporal mapping of global urban land using Landsat images based on the Google Earth Engine Platform. Remote Sens. Environ. 2018, 209, 227–239. [Google Scholar] [CrossRef]
  23. Kaplan, G.; Avdan, U. Object-based water body extraction model using Sentinel-2 satellite imagery. Eur. J. Remote Sens. 2017, 50, 137–143. [Google Scholar] [CrossRef]
  24. Jiang, H.; Wang, M.; Hu, H.; Xu, J. Evaluating the Performance of Sentinel-1A and Sentinel-2 in Small Waterbody Mapping over Urban and Mountainous Regions. Water 2021, 13, 945. [Google Scholar] [CrossRef]
  25. Yang, X.; Zhao, S.; Qin, X.; Zhao, N.; Liang, L. Mapping of Urban Surface Water Bodies from Sentinel-2 MSI Imagery at 10 m Resolution via NDWI-Based Image Sharpening. Remote Sens. 2017, 9, 596. [Google Scholar] [CrossRef]
  26. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  27. Du, Y.; Zhang, Y.; Ling, F.; Wang, Q.; Li, W.; Li, X. Water Bodies’ Mapping from Sentinel-2 Imagery with Modified Normalized Difference Water Index at 10-m Spatial Resolution Produced by Sharpening the SWIR Band. Remote Sens. 2016, 8, 354. [Google Scholar] [CrossRef]
  28. Javan, F.D.; Samadzadegan, F.; Mehravar, S.; Toosi, A.; Khatami, R.; Stein, A. A review of image fusion techniques for pan-sharpening of high-resolution satellite imagery. ISPRS J. Photogramm. Remote Sens. 2021, 171, 101–117. [Google Scholar] [CrossRef]
  29. Chavez, P.S.; Kwarteng, A.Y. Extracting spectral contrast in Landsat thematic mapper image data using selective principal component analysis. Photogramm. Eng. Remote Sens. 1989, 55, 339–348. [Google Scholar]
  30. Carper, W.J.; Lillesand, T.M.; Kiefer, R.W. The use of intensity-hue-saturation transformations for merging spot panchromatic and multispectral image data. Photogramm. Eng. Remote Sens. 1990, 56, 459–467. [Google Scholar]
  31. Chavez, P.S.; Sides, S.C.; Anderson, J.A. Comparison of 3 different methods to merge multiresolution and multispectral data—Landsat TM and SPOT panchromatic. Photogramm. Eng. Remote Sens. 1991, 57, 295–303. [Google Scholar]
  32. Laben, C.A.; Brower, B.V. Process for Enhancing the Spatial Resolution of Multispectral Imagery Using PanSharpening. Available online: http://www.google.com/patents/US6011875 (accessed on 12 June 2017).
  33. Wang, Q.; Shi, W.; Atkinson, P.M.; Zhao, Y. Downscaling MODIS images with area-to-point regression kriging. Remote Sens. Environ. 2015, 166, 191–204. [Google Scholar] [CrossRef]
  34. Jiang, H.; Feng, M.; Zhu, Y.Q.; Lu, N.; Huang, J.X.; Xiao, T. An Automated Method for Extracting Rivers and Lakes from Landsat Imagery. Remote Sens. 2014, 6, 5067–5089. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study area: (a) location in Guangzhou city; (b) true color image; (c) false color image; (d) real water body distribution.
Figure 1. Study area: (a) location in Guangzhou city; (b) true color image; (c) false color image; (d) real water body distribution.
Water 14 02696 g001
Figure 2. Workflow of the remote sensing water body extraction.
Figure 2. Workflow of the remote sensing water body extraction.
Water 14 02696 g002
Figure 3. Preservation of spectral properties of the coarse SWIR bands for the ATPRK, BIL, and GS methods.
Figure 3. Preservation of spectral properties of the coarse SWIR bands for the ATPRK, BIL, and GS methods.
Water 14 02696 g003
Figure 4. Results of the MNDWI, AWEIsh, and WI2015 water indices using the ATPRK, BIL, and GS downscaling methods.
Figure 4. Results of the MNDWI, AWEIsh, and WI2015 water indices using the ATPRK, BIL, and GS downscaling methods.
Water 14 02696 g004
Figure 5. Histograms of the MNDWI, AWEIsh, and WI2015 water indices using the ATPRK, BIL, and GS downscaling methods.
Figure 5. Histograms of the MNDWI, AWEIsh, and WI2015 water indices using the ATPRK, BIL, and GS downscaling methods.
Water 14 02696 g005
Figure 6. Misclassified results of water bodies for different combinations of water indices and downscaling methods, e.g., misclassification of building shadows in the yellow box.
Figure 6. Misclassified results of water bodies for different combinations of water indices and downscaling methods, e.g., misclassification of building shadows in the yellow box.
Water 14 02696 g006
Figure 7. Unidentified results of water bodies for different combinations of water indices and downscaling methods, e.g., omission errors of a small stream in the red box.
Figure 7. Unidentified results of water bodies for different combinations of water indices and downscaling methods, e.g., omission errors of a small stream in the red box.
Water 14 02696 g007
Table 1. Quantitative evaluation of water body extraction results in different downscaling methods and water indices combinations.
Table 1. Quantitative evaluation of water body extraction results in different downscaling methods and water indices combinations.
Water
Indices
Downscaling
Methods
UAPAOAKappa
MNDWIATPRK95.02%86.79%96.45%0.885
BIL94.12%84.84%95.92%0.867
GS93.55%87.56%96.31%0.882
AWEIshATPRK92.32%87.81%96.11%0.876
BIL91.51%87.40%95.87%0.868
GS93.80%88.70%96.57%0.891
WI2015ATPRK93.65%87.99%96.41%0.885
BIL92.45%88.54%96.27%0.881
GS94.48%89.10%96.79%0.897
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, H.; Hu, H.; Liu, X.; Jiang, H.; Liu, W.; Yin, X. A Comparison of Different Water Indices and Band Downscaling Methods for Water Bodies Mapping from Sentinel-2 Imagery at 10-M Resolution. Water 2022, 14, 2696. https://doi.org/10.3390/w14172696

AMA Style

Liu H, Hu H, Liu X, Jiang H, Liu W, Yin X. A Comparison of Different Water Indices and Band Downscaling Methods for Water Bodies Mapping from Sentinel-2 Imagery at 10-M Resolution. Water. 2022; 14(17):2696. https://doi.org/10.3390/w14172696

Chicago/Turabian Style

Liu, Haiyang, Hongda Hu, Xulong Liu, Hao Jiang, Wanxia Liu, and Xiaoling Yin. 2022. "A Comparison of Different Water Indices and Band Downscaling Methods for Water Bodies Mapping from Sentinel-2 Imagery at 10-M Resolution" Water 14, no. 17: 2696. https://doi.org/10.3390/w14172696

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