Next Article in Journal
Vulnerability Assessment as a Basis for Sanitary Zone Delineation of Karst Groundwater Sources—Blederija Spring Case Study
Previous Article in Journal
Implications of Hysteresis on the Horizontal Soil Water Redistribution after Infiltration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining Spectral Water Indices and Mathematical Morphology to Evaluate Surface Water Extraction in Taiwan

1
Department of Land Economics, National Chengchi University, No. 64, Sec. 2, ZhiNan Rd., Wenshan District, Taipei 11605, Taiwan
2
Fengyuan Land Office, Taichung City, Taichung 420211, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 4 August 2021 / Revised: 27 September 2021 / Accepted: 3 October 2021 / Published: 6 October 2021

Abstract

:
Rivers in Taiwan are characterised by steep slopes and high sediment concentrations. Moreover, with global climate change, the dynamics of channel meandering have become complicated and frequent. The primary task of river governance and disaster prevention is to analyse river changes. Spectral water indices are mostly used for surface water estimation, which separates the water from the background based on a threshold value, but it can be challenging in the case of environmental noise. Edge detection uses a canny edge detector and mathematical morphology for extracting geometrical features from the image and effective edge detection. This study combined spectral water indices and mathematical morphology to capture water bodies based on downloaded remote sensing images. From the findings, this study summarised the applicability of various spectral water body indices to the surface water extraction of different river channel patterns in Taiwan. The normalised difference water index and the modified normalised difference water index are suitable for braided rivers, whereas the automated water extraction index is ideal for meandering rivers.

1. Introduction

Taiwan is a narrow (north–south) island, with rivers following an east–west flow and having shorter flow path characteristics compared with other countries. Most of Taiwan’s rivers originate in high mountain areas; therefore, they have high flow rates and steep slopes [1]. Moreover, Taiwan is known for extreme weather events. A complex, seasonally shifting, monsoon-dominated weather system brings an abundance of rain, with a yearly average rainfall of 1000–3000 mm [2]. Because of the effects of extreme weather in recent years and the scale and frequency of rain events, the number of rainfall-related landslides and injuries has significantly increased [1]. Such weather causes the water flow to have a high sand content, which leads to intense riverbed silting. River sand mining and the construction of cross-river structures have also made river changes more frequent and complex. The prevention and control of river disasters have become critical tasks for the government. Therefore, analysing and investigating river changes must be prioritised [3,4].
In the past, river change patterns were the primary data source for river planning and management; traditional measurement techniques included control, topographic, and cross-sectional surveys [5]. River cross-sectional survey data have helped researchers to understand the actual changes in river troughs. Much effort, time, and money are required when short-term information on river changes cannot be obtained. Remote sensing data have the following characteristics: a synoptic view, the potential for remeasurement, and a wide range of spatial and temporal resolutions. The format of remote sensing data is suitable for computer processing; combined with their other advantages, this has made such data the primary source for change detection [6,7,8]. Remote sensing images include multispectral information and the band data of different ground objects can be collected for feature classification [9]. Hence, this study selected remote sensing images as a data source for river extraction. First, the inspection of river channel change is necessary for determining the characteristics of water bodies. Most water extraction methods are based on remote sensing images. Methods for determining water characteristics include supervised classification, nonsupervised classification [10,11], mixed cell decomposition [12,13,14], single-band threshold, and spectral water indices [15,16,17,18]. Spectral water indices have higher accuracy and convenience compared with other methods [19,20]. The fundamental water index normalized difference water index (NDWI) using green and near-infrared (NIR) bands of a Landsat thematic mapper (TM) image maximize water feature identification by McFeeters [16]. Xu replaced the near-infrared (NIR) band in McFeeters’s NDWI with the shortwave infrared (SWIR) band and derived the modified NDWI (MNDWI) [18] to resolve water bodies over built-up areas. However, there were still significant problems due to shadows in mountainous terrain. Feyisa et al. [15] proposed the automated water extraction index (AWEI) to identify water features, which has two conditions: AWEIsh is designed primarily to remove shadow pixels, whereas AWEInsh is designed for areas with an urban background. Extracting rivers and lakes by combining the water indices (NDWI, MNDWI, AWEIsh, and AWEInsh) was proposed as an automated method with digital image processing techniques in Landsat scenes in China [21]. Hence, in this study, we evaluated the performance of the most widely used water indices: NDWI, MNDWI, AWEIsh, and AWEInsh in a Landsat image of Taiwan. Mathematical morphology can be used to extract water body features, and it has the advantage of noise filtering and feature information retention. In particular, a continuous shape can be maintained for narrow water bodies in space [22].
The water body index and mathematical morphology have their advantages in water body extraction from images. This study combined the benefits of these two methods to extract water bodies. The spectral water body indices of the images in each period were calculated to enhance the differences between the water body and other objects. Then, a series of mathematical morphological operations were performed, and the characteristics of the water body were extracted. Subsequently, this study analysed the optimal water body index corresponding to different river channel patterns in Taiwan.

2. Materials and Methods

2.1. Data

The satellite optical images are free to download, and have regular fixed-point shooting and more spectral resolution options. This study used Landsat imagery from the United States as a source of information for detecting river changes. The data product, Level-1 TP, was downloaded from the U.S. Geological Survey (USGS). Landsat-7 was launched by NASA on 15 April 1999, with enhanced theme mapping (Enhanced Thematic Mapper Plus (ETM plus)) sensors in eight bands. Landsat-8 was launched by NASA on 11 February 2013 and is equipped with a land imager (OLI) and thermal infrared sensor (Thermal InfraRed Sensor (TIRS)) with 11 bands.
Taiwan has 26 main rivers. This study selected rivers of different types and locations: the Tamsui River in Northern Taiwan, the Beigang River in Central Taiwan, the Zengwen River and Qishan River in Southern Taiwan, and the Xiuguluan River in Eastern Taiwan. River change was the primary target of this study; for every testing area, images recorded at different times were used for detecting changes (Figure 1 and Table 1).

2.2. Method

This study’s primary goal was to identify water bodies from remote sensing images for use in flood prevention and river management. The research steps were as follows: First, the spectral water body indices of the remote sensing images of each period were calculated to enhance the difference between water bodies and other earth objects. Then, a series of mathematical morphological operations were conducted to extract the water body features. A flow chart of the research process is presented in Figure 2.

2.2.1. Spectral Water Indices

Spectral water indices are algebraic operations that combine two or more bands to enhance the difference between water bodies and nonwater bodies. Then, appropriate thresholds are established to distinguish between water bodies and other land covers; the process is simple, fast, and low-cost, especially for treating specific noise. Spectral water indices are widely used to identify water bodies [19,20]. We used four typical spectral water body indices for water body extraction in this study: the normalised difference water index (NDWI), the modified NDWI (MNDWI), the automated water extraction index shadow (AWEIsh), and the automated water extraction index no shadow (AWEInsh; see Table 2).
In the surface water extraction process, the threshold of 0 is used for all individual index maps. An optimum threshold is searched to give the highest overall accuracy and kappa coefficient based on the trial and error method. Furthermore, indices with better extraction abilities are combined to increase the separability between water and non-water surfaces based on the evaluation matrices [21].

2.2.2. Mathematical Morphology

Edges are important features in an image since they represent significant local intensity changes. They provide important clues to separate regions within an object or to identify changes. Most remote sensing applications, such as image registration, image segmentation, region separation, object description, and recognition, use edge detection as a preprocessing stage for feature extraction. Real images, such as remote sensing images, can be corrupted by point noise [23,24]. Edge detection uses a canny edge detector and mathematical morphology for extracting geometrical features from the image and for effective edge detection [25]. After data preprocessing and spectral water body index calculations, the images were processed using mathematical morphology. Mathematical morphology is an image processing technique that was proposed by Matheron and applied by Serra to image processing and analysis [26]. Compared with traditional classification, which only considers the grayscale value of pixels, mathematical morphology also considers the shape, neighbourhood, size, and other characteristics of an object to preserve its shape and edges. This is particularly suitable for water bodies with specific spatial features, resulting in feature extraction [27]. Mathematical morphology is highly effective at distinguishing between water and nonwater bodies. For example, KuPidura used mathematical morphology to distinguish between lakes and rivers and proved that mathematical morphology results are accurate [28]. Sghaier et al. combined image texture with mathematical morphology to form an algorithm for extracting rivers from high-spatial-resolution images [22]. Rishikeshan and Ramesh developed a mathematical-morphology-based semiautomated lake extraction method that preserves the actual size and shape of lakes and ponds [29]. Intervention studies involving animals or humans, and other studies that require ethical approval, must list the authority that provided approval and the corresponding ethical approval code.
This study used Matlab toolsto operate mathematical morphology [30]. Median filtering is the first step, which is a nonlinear operation often used in image processing to reduce salt and pepper noise. It makes the image smoother and intensifies pixels, mainly background pixels close to each other. The second step is to determine the size of the structural elements, which depends on satellite imagery at different spatial resolutions (15 m in this study). The third and fourth steps are dilation and erosion operating “lengthening”, “thickening”, “shrinks”, and “thins” in the binary image, which achieve control by a structural element. The last two steps are opening and closing operations. The opening operation makes the contour of objects smooth and disconnects narrow discontinuous, and removes thin protrusions. Similarly, with the opening operation, the closing operation also makes the outline smooth. The opposite is that it usually eliminates discontinuity and narrows long, thin gaps; clears up small holes; and fills the ruptures of the contour line. (For theoretical principles and foundations, please refer to [31]).

2.3. Accuracy Assessment

After all the research procedures, the water body results were compared with reference images (i.e., the ground truth data). To verify the accuracy of the results, the reference images were recorded at the same time as the result images. In this study, the reference data of water bodies and nonwater bodies in colour images were collected through random sampling to evaluate the accuracy of the classification results. Landsat 8 colour images were used as reference images, and random sampling was used to collect reference data on water and nonwater bodies. Sampling points were established using ArcGIS software, in which random sampling methods comprise the following options:
  • Simple random sampling: The samples are randomly distributed throughout the image.
  • Stratified random sampling: A random number of samples is distributed throughout the image; the number of samples for each category is proportional to the area of the image.
  • Equalised stratified random sampling: A randomly distributed sample count is created throughout the image, with each category having the same sample size.
The interest category of the topographic features was the water bodies. To ensure accurate and reliable classification, images were selected in which water bodies accounted for a lower proportion of the whole image compared with other objects. Hence, we used equalised stratified random sampling to produce sampling points. The water body and nonwater sample sizes were the same to avoid the sample number being concentrated in the nonwater body category. For each category, a sufficient number of samples must be collected to accurately assess the classification accuracy when using random sampling. However, a crucial consideration is the number of sampling points, which must be determined to balance accuracy with the cost of collecting reference data. In general, at least 50 random samples are required for each category when the study area is less than 1 million acres and there are fewer than 12 categories [32]. At least 75–100 random samples should be collected in each category when the study area is large or the surface distribution is more complex (100 samples in this study). These criteria were verified as statistically valid and helpful.
In this study, a confusion matrix was used as an evaluation indicator, the results of which can be divided into four categories, as presented in Table 3:
(1)
True positive (TP) represents correctly identified water pixels;
(2)
True negative (TN) indicates the correctly rejected nonwater pixels;
(3)
False positive (FP) indicates incorrectly identified water pixels;
(4)
False negative (FN) indicates undetected water pixels.
The values in the confusion matrix can be used to calculate the producer accuracy (PA), user accuracy (UA), overall accuracy (OA), and kappa coefficient, which are detailed in Table 4.
In Table 4, ∑ is the chance accuracy represented by (TP + FP)(TP + FN) + (FN + TN)(FP + TN), where T is the total number of pixels. The accuracy indices of PA, UA, and OA range from 0 to 1, with values closer to 1 indicating more accurate results. However, these indices do not account for random agreement between datasets. Hence, the kappa coefficient, a tool for controlling the random agreement factor, is often used with PA, UA, and OA. Usually, the kappa coefficient can range from −1 to +1, where 0 represents the agreement that is expected from random chance, and 1 represents perfect agreement between the raters [33].

3. Result and Discussion

3.1. Element Sizes in Water Body Extraction

An image’s spatial resolution affects the sizes of structural elements when the characteristics of water bodies are extracted using spectral water body indices combined with mathematical morphology; specifically, Landsat images with a spatial resolution of 15 m were used in this study. Because morphological operators are based on minimum/maximum computations, they require a total ordering relationship between the values of the pixels falling within the neighbourhood defined by the structuring element. This study used three consecutive years (2013 to 2015) of image data from the Beigang River area to find the most suitable structural element size, which also calculates OA of each water index; the results are presented in Table 5 and Figure 3.
With a structural element size of three, the Beigang River water body was entirely extracted using mathematical morphology. However, some smaller pixels were misjudged as water bodies, as indicated by the circles in Figure 3. The accuracy of each spectral water index was higher than those of other structural element sizes. The Beigang River could also be fully extracted when the structural element size was four. When the structural element size was five, however, water bodies were determined to be nonwater, as indicated by the circles in Figure 3. In summary, for extracting water bodies using water body indices combined with mathematical morphology, the most suitable structural element size was four.

3.2. Results of Water Body Extraction

The original images were processed using radiation correction, atmospheric correction, and image fusion. NDWI, MNDWI, AWEIsh, and AWEInsh were used to calculate and enhance the differences between water bodies and nonwater bodies in the images. Subsequently, mathematical morphological calculations were performed, and the water features were finally extracted. The accuracy of water body classification was supplemented by Google Earth images with an accuracy assessment and confusion matrix through random sampling. The total number of samples collected was 100. Then, the PA, UA, OA, and kappa coefficient were calculated based on the confusion matrix of the four spectral water body indices. The OA and kappa coefficient were then used as evaluation criteria to determine the quality of the classification results.
  • Tamsui River
As indicated in Table 6, the MNDWI performed the best in terms of OA for the Tamsui River watershed. MNDWI, using an MIR band instead of an NIR band, can considerably improve the enhancement in open water features. It can quickly and accurately discriminate water from nonwater features [18]. Based on the results, the MNDWI is particularly suitable for urban areas with both buildings and water bodies. It can enhance the differences between artificial objects and water bodies to avoid misjudgements, as shown in Figure 4.
2.
Beigang River
The AWEI indices were optimal for the Beigang River in terms of OA (Table 7). The studied area is located in the downstream plain, with less shadow-free problems. AWEIsh could not fully depict tributaries and buildings, which it misjudged as water bodies; therefore, the AWEInsh is more suitable for the Beigang River (Figure 5).
3.
Zengwen River
In terms of OA in the Zengwen River basin, the AWEI indices had the best performance (Table 8). The study area includes downstream plains and midstream areas of hillsides. When the water extraction results were examined, this study found that AWEIsh could show continuous water flows and contributed to the analysis of subsequent change detection (Figure 6).
4.
Qishan River
The Qishan River is located in an upstream area with mostly mountainous terrain and only a few artificial buildings. In terms of OA, the NDWI exhibited the best performance for the Qishan River (Table 9). Based on the water extraction results, the NDWI successfully distinguished between water bodies and other objects, and then fully extracted the river channel of the Qishan River, as detailed in Figure 7.
5.
Xiuguluan River
Similar to the Qishan River, the Xiuguluan River is mostly located in mountainous areas. In addition to artificial buildings, there are ponds and agricultural land. As detailed in Table 10, the NDWI had the best performance in terms of OA for the Xiuguluan River. Based on the water extraction results, the NDWI eliminated the mountain shadow problem, as detailed in Figure 8.

3.3. Result of River Changing

The findings of this study were used to detect river changes between two epochs in the testing areas. It is easy to see the difference between epochs, based upon which river management and government can be conducted. For the Tamsui River watershed, we used the MNDWI and mathematic morphology to extract the water body of two epochs. The difference between the two epochs is shown in Figure 9.
The AWEInsh and mathematic morphology were more suitable for the Beigang River to extract water bodies. The difference between the two epochs is shown in Figure 10.
The AWEIsh had the best performance in the Zengwen River, including downstream plains and midstream areas of hillsides. The difference between the two epochs is shown by the changes in the river channel (Figure 11).
The Qishan River is located in an upstream area with mostly mountainous terrain and only a few artificial buildings; the NDWI and I and mathematic morphology exhibited the best performance for the Qishan River. The difference between the two epochs is demonstrated by the changes in the river channel (Figure 12).
For the Xiuguluan River, using NDWI and mathematic morphology had the best performance in extracting the water bodies of the two epochs. The difference between the two epochs is shown in Figure 13.

4. Conclusions

Much effort, time, and money are required to determine river change patterns for river planning and management. Remote sensing data have the potential to be used for remeasurement with a wide range of spatial and temporal resolutions. This study integrated spectral water body indices and mathematical morphology to extract the water bodies of the watershed from remote sensing images of Northern, Central, Southern, and Eastern Taiwan. The results demonstrated the performance of each spectral water body index in every watershed as well as which index is most suitable for each type of river channel. Based on the results, this study drew the following main conclusions:
  • NDWI and MNDWI are ideal for braided rivers (e.g., Qishan River, Xiuguluan River, and Tamsui River).
  • AWEInsh should be applied to meandering rivers (e.g., Beigang River and Zengwen River).
  • The various spectral water body indices differ in their application. Therefore, the contribution of this study is mainly for Taiwan’s rivers, summarising the applicability of various spectral water bodies indexes in water body extraction.
  • In this study, pairs of remote sensing images were used to identify the water bodies of river channels to observe changes; the results may facilitate river management.

Author Contributions

Conceptualization, F.-S.N.; data curation, Y.-C.L.; formal analysis, F.-S.N. and Y.-C.L.; methodology, F.-S.N. and Y.-C.L.; project administration, F.-S.N.; validation, F.-S.N.; writing—original draft, F.-S.N.; writing—review and editing, F.-S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the TAIWAN STUDIES CENTER, NATIONAL CHENGCHI UNIVERSITY.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Water Resources Agency, MOEA. Water Environment in Taiwan. 2021. Available online: http://eng.wra.gov.tw/cl.aspx?n=5118 (accessed on 1 April 2021).
  2. Central Weather Berau. Weather FAQ. 2021. Available online: https://opendata.cwb.gov.tw/index (accessed on 1 April 2021).
  3. Louw, E.; Olanrewaju, C.C.; Olanrewaju, O.A.; Chitakira, M. Impacts of flood disasters in Nigeria: A critical evaluation of health implications and management. Jàmbá: J. Disaster Risk Stud. 2019, 11, 1–9. [Google Scholar]
  4. Tingsanchali, T. Urban flood disaster management. Procedia Eng. 2012, 32, 25–37. [Google Scholar] [CrossRef] [Green Version]
  5. Boon, P.; Paul, R. (Eds.) River Conservation and Management; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  6. Chen, G.; Hay, G.J.; Carvalho, L.M.; Wulder, M.A. Object-based change detection. Int. J. Remote Sens. 2012, 33, 4434–4457. [Google Scholar] [CrossRef]
  7. Coops, N.C.; Wulder, M.A.; White, J.C. Identifying and describing forest disturbance and spatial pattern: Data selection issues and methodological implications. In Understanding Forest Disturbance and Spatial Pattern; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  8. Lunetta, R.S.; Johnson, D.M.; Lyon, J.G.; Crotwell, J. Impacts of imagery temporal frequency on land-cover change detection monitoring. Remote Sens. Environ. 2004, 89, 444–454. [Google Scholar] [CrossRef]
  9. 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]
  10. Mciver, D.K.; Friedl, M.A. Using prior probabilities in decision-tree classification of remotely sensed data. Remote Sens. Environ. 2002, 81, 253–261. [Google Scholar] [CrossRef]
  11. Sivanpillai, R.; Miller, S.N. Improvements in mapping water bodies using ASTER data. Ecol. Inform. 2010, 5, 73–78. [Google Scholar] [CrossRef]
  12. Frazier, P.S.; Page, K.J. Water body detection and delineation with Landsat TM data. Photogramm. Eng. Remote Sens. 2000, 66, 1461–1468. [Google Scholar]
  13. Klein, I.; Dietz, A.J.; Gessner, U.; Galayeva, A.; Myrzakhmetov, A.; Kuenzer, C. Evaluation of seasonal water body extents in Central Asia over the past 27 years derived from medium-resolution remote sensing data. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 335–349. [Google Scholar] [CrossRef]
  14. Work, E.A.; Gilmer, D.S. Utilization of satellite data for inventorying prairie ponds and lakes. Photogramm. Eng. Remote Sens. 1976, 42, 685–694. [Google Scholar]
  15. Feyisa, G.L.; Meilby, H.; Fensholt, R.; Proud, S.R. Automated Water Extraction Index: A new technique for surface water mapping using Landsat imagery. Remot. Sens. Environ. 2014, 140, 23–35. [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. Rogers, A.S.; Kearney, M.S. Reducing signature variability in unmixing coastal marsh Thematic Mapper scenes using spectral indices. Int. J. Remote Sens. 2004, 25, 2317–2335. [Google Scholar] [CrossRef]
  18. 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]
  19. Acharya, T.D.; Subedi, A.; Lee, D.H. Evaluation of water indices for surface water extraction in a Landsat 8 scene of Nepal. Sensors 2018, 18, 2580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Jiang, H.; Feng, M.; Zhu, Y.; Lu, N.; Huang, J.; Xiao, T. An automated method for extracting rivers and lakes from Landsat imagery. Remote Sens. 2014, 6, 5067–5089. [Google Scholar] [CrossRef] [Green Version]
  21. Du, Z.; Bin, L.; Ling, F.; Li, W.; Tian, W.; Wang, H.; Zhang, X. Estimating surface water area changes using time-series Landsat data in the Qingjiang River Basin, China. J. Appl. Remote Sens. 2012, 6, 063609. [Google Scholar] [CrossRef]
  22. NagaRaju, C.; NagaMani, S.; Rakesh Prasad, G.; Sunitha, S. Morphological edge detection algorithm based on multi-structure elements of different directions. Int. J. Inf. Commun. Technol. Res. 2011, 1, 37–43. [Google Scholar]
  23. Sghaier, M.O.; Foucher, S.; Lepage, R. River extraction from high-resolution sar images combining a structural feature set and mathematical morphology. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 10, 1025–1038. [Google Scholar] [CrossRef]
  24. Ali, M.; David, C. Using the Canny edge detector for feature extraction and enhancement of remote sensing images. In Proceedings of the IGARSS 2001. Scanning the Present and Resolving the Future. Proceedings IEEE 2001 International Geoscience and Remote Sensing Symposium (Cat. No. 01CH37217), Sydney, NSW, Australia, 9–13 July 2001; pp. 2298–2300. [Google Scholar]
  25. Liu, Q.; Dong, J.; Liu, G.; Huang, C.; Xie, C. Using the canny edge detector and mathematical morphology operators to detect vegetation patches. In Proceedings of the Third International Conference on Digital Image Processing (ICDIP 2011), Chengdu, China, 15–17 April 2011. [Google Scholar]
  26. Dougherty, E. Mathematical Morphology in Image Processing; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  27. Rishikeshan, C.A.; Ramesh, H. An automated mathematical morphology driven algorithm for water body extraction from remotely sensed images. ISPRS J. Photogramm. Remote Sens. 2018, 146, 11–21. [Google Scholar] [CrossRef]
  28. Kupidura, P. Distinction of lakes and rivers on satellite images using mathematical morphology. Biul. Wojsk. Akad. Tech. 2013, 62, 57–69. [Google Scholar]
  29. Rishikeshan, C.A.; Ramesh, H. An ANN supported mathematical morphology based algorithm for lakes extraction from satellite images. ISH J. Hydraul. Eng. 2018, 24, 222–229. [Google Scholar] [CrossRef]
  30. Putera, S.H.I.; Ibrahim, Z. Printed circuit board defect detection using mathematical morphology and MATLAB image processing tools. In Proceedings of the 2010 2nd International Conference on Education Technology and Computer, Shangai, China, 22–24 June 2010. [Google Scholar]
  31. Suresha, D.; Ganesh, V. A Survey-Mathematical Morphology operations on Images in MATLAB. Int. J. Adv. Sci. Res. Technol. 2012, 2, 37–43. [Google Scholar]
  32. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 1991, 37, 35–46. [Google Scholar] [CrossRef]
  33. McHugh, M.L. Interrater reliability: The kappa statistic. Biochem. Med. 2012, 22, 276–282. [Google Scholar] [CrossRef]
Figure 1. Locations of testing rivers.
Figure 1. Locations of testing rivers.
Water 13 02774 g001
Figure 2. Research flow chart.
Figure 2. Research flow chart.
Water 13 02774 g002
Figure 3. Characteristics of the water body of Beigang River extracted using differently sized structural elements.
Figure 3. Characteristics of the water body of Beigang River extracted using differently sized structural elements.
Water 13 02774 g003
Figure 4. Tamsui River: Results of spectral water body indices for water body extraction at different times.
Figure 4. Tamsui River: Results of spectral water body indices for water body extraction at different times.
Water 13 02774 g004
Figure 5. Beigang River: Results of spectral water body indices in terms of water body extraction at different times.
Figure 5. Beigang River: Results of spectral water body indices in terms of water body extraction at different times.
Water 13 02774 g005
Figure 6. Zengwen River—Results of spectral water body indices in terms of water body extraction at different times.
Figure 6. Zengwen River—Results of spectral water body indices in terms of water body extraction at different times.
Water 13 02774 g006
Figure 7. Qishan River: Results of spectral water body indices in terms of water body extraction on different dates.
Figure 7. Qishan River: Results of spectral water body indices in terms of water body extraction on different dates.
Water 13 02774 g007
Figure 8. Xiuguluan River: Results of spectral water body indices in terms of water body extraction on different dates.
Figure 8. Xiuguluan River: Results of spectral water body indices in terms of water body extraction on different dates.
Water 13 02774 g008
Figure 9. The Tamsui River watershed changes between two epochs.
Figure 9. The Tamsui River watershed changes between two epochs.
Water 13 02774 g009
Figure 10. The Beigang River changes between two epochs.
Figure 10. The Beigang River changes between two epochs.
Water 13 02774 g010
Figure 11. The Zengwen River changes between two epochs.
Figure 11. The Zengwen River changes between two epochs.
Water 13 02774 g011
Figure 12. The Qishan River changes between two epochs.
Figure 12. The Qishan River changes between two epochs.
Water 13 02774 g012
Figure 13. The Xiuguluan River changes between two epochs.
Figure 13. The Xiuguluan River changes between two epochs.
Water 13 02774 g013
Table 1. Data of testing images.
Table 1. Data of testing images.
Testing AreaDateSatellite
Tamsui River3 December 2001Landsat-7 ETM+
23 December 2017Landsat-8 OLI
Beigang River17 November 2013Landsat-8 OLI
15 November 2018
Zengwen River3 December 2013Landsat-8 OLI
30 October 2018
Qishan River26 December 2001Landsat-7 ETM+
28 November 2017Landsat-8 OLI
Xiuguluan River14 November 2000Landsat-7 ETM+
13 December 2019Landsat-8 OLI
Table 2. Spectral water indices for water feature extraction.
Table 2. Spectral water indices for water feature extraction.
Spectral Water IndexEquation
Normalised difference water index (NDWI) NDWI = B g r e e n B N i r / B g r e e n + B N i r
Modified normal difference water index (MNDWI) M N D W I = ( B g r e e n B M I R ) / ( B g r e e n + B M I R )
Automated water extract index shadow (AWEIsh) A W E I s h = B b l u e + 2.5 × B g r e e n   1.5 × B N i r + B S W I R 1   0.25 × B S W I R 2
Automated water extract index no shadow (AWEInsh) A W E I n s h = 4 × B g r e e n B S W I R 1 0.25 × B N i r + 2.75 × B S W I R 2
Table 3. Confusion matrix.
Table 3. Confusion matrix.
Classified DataWaterNon-Water
WaterTPFP
Non-waterFNTN
Table 4. Indicators in the accuracy evaluation.
Table 4. Indicators in the accuracy evaluation.
Accuracy Evaluation MetricsEquationDescription
Producer accuracy P A = T P T P + F N The ratio of the correctly classified pixels to the actual total number of pixels for each category.
User accuracy U A = T P T P + F P The ratio of the correctly classified pixels to the total pixels categorised for each category.
Overall accuracy O A = T P + T N T The ratio of the correctly classified pixels to the total number of pixels.
Kappa coefficient K a p p a = T ( T P + T N ) T 2 The percentage of errors reduced by classification results compared with completely random classifications.
Table 5. Overall accuracy of the Beigang River area based on the different sizes of structural elements.
Table 5. Overall accuracy of the Beigang River area based on the different sizes of structural elements.
Spectral Water Index Accuracy(%)NDWIMNDWIAWEIshAWEInsh
Element SizeYear
3201373.977.888.996.7
201473.376.786.191.1
201578.378.385.091.1
4201373.376.192.896.1
201475.076.791.194.4
201582.277.887.293.9
5201371.175.094.495.0
201476.776.183.390.0
201577.477.892.292.2
Table 6. Tamsui River: Accuracy of spectral water body indices in water body extraction on different dates.
Table 6. Tamsui River: Accuracy of spectral water body indices in water body extraction on different dates.
Water
Index
3 December 200123 December 2017
PAUAOAKappaPAUAOAKappa
NDWI0.360.81820.640.280.510.740.4898
MNDWI11110.961510.980.96
AWEIsh0.5810.790.580.730810.860.7227
AWEInsh0.4410.720.440.461510.720.4514
Table 7. Beigang River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Table 7. Beigang River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Water
Index
17 November 201315 November 2018
PAUAOAKappaPAUAOAKappa
NDWI0.88890.72730.760.50740.7170.79170.750.5012
MNDWI0.9630.74290.80.58680.94340.80650.850.6956
AWEIsh0.90740.96080.930.85980.86790.95830.910.8204
AWEInsh0.90740.980.940.880.943410.970.94
Table 8. Zengwen River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Table 8. Zengwen River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Water
Index
3 December 201330 October 2018
PAUAOAKappaPAUAOAKappa
NDWI0.98250.75680.810.59260.96360.8030.850.6901
MNDWI0.98250.80.850.68220.96360.86890.90.7955
AWEIsh0.877210.930.860.909110.950.9
AWEInsh0.91230.98110.940.8790.92730.96230.940.8793
Table 9. Qishan River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Table 9. Qishan River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Water
Index
26 December 200128 November 2017
PAUAOAKappaPAUAOAKappa
NDWI0.892910.940.880.781310.860.72
MNDWI0.80360.95740.870.74190.79690.96230.850.6949
AWEIsh0.39290.91670.640.32230.57810.97370.720.4753
AWEInsh0.50.44440.37−0.30060.6250.63490.53−0.0138
Table 10. Xiuguluan River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Table 10. Xiuguluan River: Accuracy of spectral water body indices in terms of water body extraction on different dates.
Water
Index
14 November 200013 December 2019
PAUAOAKappaPAUAOAKappa
NDWI0.90790.97180.93570.87130.909110.950.9
MNDWI0.89470.97140.92860.85710.94550.9630.950.8992
AWEIsh0.80260.96830.88030.76240.69090.82610.750.504
AWEInsh0.855310.92140.84380.763610.870.7441
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ning, F.-S.; Lee, Y.-C. Combining Spectral Water Indices and Mathematical Morphology to Evaluate Surface Water Extraction in Taiwan. Water 2021, 13, 2774. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192774

AMA Style

Ning F-S, Lee Y-C. Combining Spectral Water Indices and Mathematical Morphology to Evaluate Surface Water Extraction in Taiwan. Water. 2021; 13(19):2774. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192774

Chicago/Turabian Style

Ning, Fang-Shii, and Yu-Chan Lee. 2021. "Combining Spectral Water Indices and Mathematical Morphology to Evaluate Surface Water Extraction in Taiwan" Water 13, no. 19: 2774. https://0-doi-org.brum.beds.ac.uk/10.3390/w13192774

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