Comparison of accuracies for different water extraction methods, with the best results given in bold text.
An improved water extraction method using a morphological linear enhancement technique is proposed to improve the delineation of narrow water features for the modified normalized difference water index (MNDWI) derived from remote sensing images. This method introduces a morphological white top-hat (WTH) transforming operation on the MNDWI to extract multi-scale and multidirectional differential morphological profiles and constructs a morphological narrow water index (MNWI). The MNWI can effectively enhance the local contrast of linear objects, allowing narrow water bodies to be easily separated from mountain shadows and other features. Furthermore, to accurately delineate surface water bodies, a dual-threshold segmentation method was also developed by combining an empirical threshold segmentation with the MNDWI for wide water bodies and an automatic threshold segmentation with the MNWI for narrow water bodies. This method was validated using three experimental datasets, which were taken from two different Landsat images. Our results demonstrate that narrow water bodies can be sufficiently identified, with an overall accuracy of over 90%. Most narrow streams or rivers keep a continuous shape in space, and the boundaries of the water bodies are accurately delineated as compared with the MNDWI method. Finally, the proposed method was used to extract the entire inland surface water of Fujian province, China.
- narrow water extraction
- white top-hat transform
- dual-threshold segmentation
Surface water is one the most vital earth resources undergoing changes in time and space as a consequence of land use/cover (LULC) changes, climate change, and other forms of environmental changes in many parts of the world. Timely and accurate monitoring and delivery of data of the dynamics of surface water are, therefore, critically important in various scientific disciplines, such as the assessment of present and future water resources, climate models, agricultural suitability, river dynamics, wetland inventory, watershed analysis, surface water surveys, and environmental monitoring .
Remote sensing at different spatial, spectral, and temporal resolutions provides an enormous amount of data for mapping water resources and its dynamics at local to global scales. At a result, it has become a routine approach for the monitoring of land surface water bodies, since the acquired data can provide macroscopic, real-time, dynamic, and cost-effective information, which is substantially different from conventional in situ measurements. Various approaches for water body extraction from multispectral images have been developed in the past decades [2, 3, 4], which can be broadly grouped into three categories: spectral band segmentation, image supervised classification, and water indices. Among all these methods, of particular interest is the spectral water index-based method, as it is a reliable and cost-effective method. This type of method takes advantage of reflectivity differences of each involved band for water body extraction based on the analysis of signature differences between water and other surfaces.
One of the most widely used indexes is the normalized difference water index (NDWI) , which utilizes the green (band 2) and near-infrared (band 4) of Landsat TM to delineate open water features. However, Xu found that the NDWI cannot efficiently suppress the signal from built-up surfaces and therefore proposed an improved one, called modified normalized difference water index (MNDWI) , where the NDWI was modified by replacing band 4 with band 5 of Landsat TM/ETM. The MNDWI has been validated as one of the most widely used water indices for various applications, though it is still difficult to obtain a high accuracy of water extraction in complex circumstances. Carleer and Wolff , among others, have found that the land cover classifications of water and shadow can often be confused. This issue often arises in environments where a large amount of shadow and water regions exist, such as urban and mountainous landscapes. The identification of narrow water bodies (such as narrow streams, canals, ponds, small reservoirs, etc.) can be a difficult task when using NDWI or MNDWI images, because the shallow and narrow water pixels may generate unstable spectral profiles or characteristics, due to the mixed reflectance caused by sediment and/or adjacent land covers. Narrow water is typically defined as a water body with an apparent width less than or equal to three pixels in an image. Therefore, narrow water features often contain mixed pixels, and the extraction of them from NDWI or MNDWI images generally exhibits a discontinuous shape in space.
To remedy this problem, past studies have attempted to identify narrow water features by combining different procedures. Yang et al. proposed a method of extracting initial water information via a user-defined specified water index, and then they performed a series of operations. These operations include morphological dilation, image filtering, and thinning techniques applied to the water index image to recover the continuity of narrow rivers or streams . This method can be effective in the extraction of narrow water bodies if the water disruption is short; however, it may increase false water identification when water disruption is large, since it is dependent on the morphological dilation operation to reconnect the narrow rivers. Such simple threshold techniques are not often a sufficient solution to identify narrow water bodies; therefore, Li et al. suggested an object-oriented method of small water body extraction . They first extracted textural and shape-related features from images as supplementary information to spectral bands and then performed a segmentation operation on the images using an optimal scale to identify the potential water bodies. Yet, their method is not an automatic process, since it involves multiple user-defined parameters in image segmentation, which prohibits its use in large areas. An alternative approach was performed by Jiang et al. who extracted narrow water features via the enhancement of linear features in NDWI images . However, their procedure involves multiple empirical thresholds, so it is not a cost-effective method for water feature extraction on large scale.
Attempting to improve on these previous approaches, here we propose an automatic water extraction method that constructs a novel narrow water index, denoted as morphological narrow water index (MNWI). The MNWI is constructed using multi-scale and multidirectional differential morphological profiles (DMPs) on a MNDWI image, and then water bodies are automatically extracted using a dual-threshold segmentation. The successful use of DMPs to extract various thematic information from images has been sufficiently validated, such as buildings in urban areas [10, 11], rare earth mining areas , and mapping of mangrove forests . In this paper, we introduce a DMP technique to highlight the contrast of bright features as a way of narrow water recognition in the MNDWI images. This can be accomplished because water pixels have higher value than surrounding pixels in MNDWI images. Our approach is expected to improve the ability of narrow water feature identification by enhancing its spatially implicit characteristics using multi-scale morphological features, e.g., other land cover uses that have similar values in a MNDWI image, such as bare patches and shadows, can be easily identified. Our approach also involves a dual-threshold strategy which is adopted for wide and narrow water body extraction, since a simple threshold is not often an adequate solution . An empirical threshold is used first to obtain possible water areas from a MNDWI image, followed by an automatic threshold which is determined by the maximum interclass variance criterion  used for extracting narrow water features from a MNWI image. Finally, a logical operation is performed by combining the two potential water features to identify the true water body boundary.
The remainder of this chapter is organized as follows. In Section 2, we describe the MNWI method of narrow water extraction. Experimental results are shown in Section 3 using multiple experiments on TM images, and we use the MNWI method in a practical application for extracting inland water features in Fujian province, China, in Section 4. Finally, Section 5 presents conclusions.
2. The proposed MNWI method
It is well accepted that open and wide water features can be easily separated from other land cover features by using the NDWI or the MNDWI methods, but extracting narrow water boundaries is generally a more difficult task due to its being confused with built-up areas, roads, hill shadows, etc. Therefore, the goal of our proposed MNWI method is the separation of narrow water and other land cover features by depicting the implicit spectral and structural characteristics of MNDWI. Narrow water usually exhibits strong linear shapes and continuous spatial curves; therefore, we propose a linear enhancement operation on a MNDWI image to form a MNWI image, according to the following steps:
where Green and MIR are the image reflectance of the green band and medium-wave infrared band (which correspond to the TM/ETM+ band 5), respectively.
where is the output image yielded by performing a closing operator to a MNDWI image. A closing operator can suppress smaller, darker objects and join adjacent objects together; thus we expect that small water bodies smaller than the structural elements will be highlighted after this reconstruction, and open water bodies and background objects larger than the structural elements should be suppressed. We previously defined a narrow water body as no larger than three pixels; thus, the linear structural elements in our algorithm are extended in four directions (0°, 45°, 90°, and 135°), and each direction was applied with three scales (smin = 1, smax = 3, Δs = 1) to generate the WTH image.
where WTHMax and WTHMin are the maximum and the minimum values of the reconstruction of the white top-hat procedure in all directions, respectively.
3. Method validations
To evaluate our method, three study regions are taken from two Landsat ETM+/OLI images with different water body types and different terrains, as shown in Figure 1. Study area 1 is a sub-scene image with 1000 by 1000 pixel size chosen from a Landsat 7 ETM+ image acquired on September 1, 2001, centered on the Panjiakou Reservoir near Tangshan and Chengde cities in Hebei province in northern China, which covers multiple branches of Luanhe River.
Study area 2 is also a sub-image of 1000 by 1000 pixels, which is located in Luoyuan County, Fujian province, which is acquired from the Landsat 8 operational Landsat image (OLI) on December 13, 2014. This region contains very narrow streams (Figure 2b), and it is used to test the extraction ability from mixed pixels. Study area 3 (Figure 3c) is a 1000 by 1000 pixel scene selected from the same OLI image as #2, which is located in southern Youxi County, Fujian province, and covers one of main branches of Minjiang River as well as other narrow streams (Qingyin, Qing, Wenjiang, etc.). Youxi County and several villages are in this study area, which contain multiple sources of possible background noise.
Actual water feature information for these study areas are not available, so the water bodies in the three images are manually digitized using high resolution spatial images to provide a basis map for comparison. High-resolution Google Earth™ images were also used as a complementary reference to assist in distinguishing water pixels that might be confused with background noise, such as mountain shadows, trails and built-up areas.
3.1 Validation of MNWI features
A comparison between MNDWI and MNWI methods was first described. The two types of water index images for study area 1 are shown in Figure 3a and b, with focused areas highlighted in red boxes shown in Figure 3c and d. The pixels’ values of water typically have higher values (white areas) than those of the other land cover features in these images, since both the MNDWI and the MNWI water indexes strongly enhance the water body signals. It is also observed in Figure 3a and c that the brightness difference between wide water and narrow water in the MNDWI image is larger, suggesting that a threshold segmentation on the MNDWI image cannot resolve entire water bodies. In contrast, the difference in the MNWI image is significantly reduced, though they maintain relatively higher values than other land cover features, as shown in Figure 3b and d. By looking at Figure 3d, it is seen that the local contrast between narrow water and other land cover features is significantly enhanced. Additionally, narrow water such as small rivers and branches maintains a continuous spatial shape, suggesting that narrow water features can be accurately extracted with a threshold implementation.
Furthermore, 1200 samples of typical land cover types from study area 3 were randomly selected from each category and were analyzed by calculating the maximal value, the minimal value, the mean value, and the deviation. The criteria for the sample selection are the following: (1) Each land cover has ~200 samples to keep a sample balance; and (2) each land cover contains several small patches from different locations to maintain a spectral variety. Figure 4 reports the spatial distribution of the samples for study area 3 and their statistical information for six typical land cover types, i.e., wide open water, narrow water, vegetation, built-up area, roads, and shadow. As can be seen in Figure 4, the values of the wide water are very high for MNDWI. However, it is difficult to discriminate narrow water from other land cover types, especially for shadows, roads, and built-up areas. In contrast, the values of the narrow water in the MNWI method are relatively high compared with other land cover types, indicating that it is relatively easy to separate the narrow water bodies from other land cover types.
It can also be seen in Figure 5 that the values of objects with polygon-like shapes, such as wide water, built-up areas, and shadows, are heavily suppressed in the MNWI image since they do not exhibit a linear structure. However, roads also show linear structural characteristics; thus they have high values in a MNWI image. This demonstrates that neither MNDWI nor MNWI can effectively identify entire water bodies using a threshold segmentation. To remedy this, we adopted a dual-threshold strategy. The first threshold is used for wide water extraction from MNDWI, and the second is employed to extract narrow water features from MNWI.
3.2 Validation of dual-threshold segmentation
An empirical threshold (0.2) was first used to perform a rough extraction of potential water features from a MNDWI image; then a second threshold determined by Otsu was determined from the MNWI image for possible narrow water features. Next, a combing procedure was carried out to extract entire water bodies using an “if-then” logic calculation according to the following rule: If two potential water regions are spatially connected, then they are determined to be water bodies; otherwise they are determined to be other land cover types.
As a demonstration of the dual-threshold segmentation method, comparisons between single threshold segmentation of a MNDWI image and dual-threshold segmentation of a MNWI image are shown in Figure 6(a)–(c). It is seen that when a smaller threshold (T = 20) is adopted, a narrow stream keeps a relatively complete spatial shape, yet it also contains a trail (road) at the bottom of the image (Figure 6a). However, if we increase the threshold T to 40, this trail is no longer extracted, but the stream exhibits discontinuities along the stream. Thus, the use of a dual-threshold segmentation can better identify this stream as continuous and avoid the identification of the trail.
3.3 Visual assessment
Figure 7 presents the extracted water features using our proposed method for each of the three study areas. As a comparison, water information derived from the MNDWI image using an optimal threshold segmentation is also listed. For clarity, the corrected, misclassified, and omitted water information is labeled with different color schemes. Visual inspection shows that our method significantly outperforms the MNDWI method when using an optimal threshold segmentation. It can be seen that the majority of narrow rivers in each of the study areas are successfully extracted by our method. Six branches of the Luanhe River are clearly delineated in study area 1, and most narrow streams in study areas 2 and 3 are clearly extracted. However, a few narrow tributaries were misclassified by the MNDWI method, highlighted in red in Figure 7. A closer inspection of study area 2 shows that there are still two omissions which are highlighted with white rectangles, due to the width of two streams being too narrow (less than 10 m) to occupy a footprint, and the reflectance of these pixels are strongly mixed with other land cover types. Conversely, the results derived from the MNDWI image are less effective, as only small portions of the narrow rivers were extracted correctly. This is especially true for the narrow streams in the top region of the image, as most of them are ignored.
Another misclassification issue is the delineation of the sides of rivers, due to mixed pixel effects. Another experiment in the study area 3 was conducted to demonstrate this. These results are reported in Figure 8, where the water information that was corrected, omitted, and misclassified is shown in cyan, magenta, and red, respectively. It can be found that the boundary of the Youxi River can be accurately extracted with the proposed method as shown in Figure 8b. However, using the MNDWI method with an optimal threshold, the omitted water pixels were along both sides of the river boundary (Figure 8c).
3.4 Quantitative evaluation
We now quantitatively evaluate our extracted results. Four measurements are used for comparison; the user and producer accuracy, the kappa coefficient, and the overall accuracy. Additionally, two recently developed methods for narrow water extraction, i.e., the method developed by Yang et al.  and the linear feature enhancement (LFE) developed by Jiang et al. , are also included for comparison. Table 1 gives a pixel-by-pixel analysis of the classification accuracies for all datasets, with the best results highlighted in bold.
|Study area 1||Producer accuracy||53.8%||75.1%||92.1%|
|Study area 2||Producer accuracy||35.2%||43.6%||82.8%|
|Study area 3||User accuracy||87.9%||73.8%||96.6%|
It can be seen in Table 1 that the optimal threshold segmentation method was the least accurate, as it failed nearly completely in study area 2 with a 35% product accuracy and a kappa coefficient of 0.477. The method of Yang has some similar effects in narrow water extraction to our datasets, especially for mixed water features. The method of Jiang significantly improves the accuracy of narrow water extraction as compared to the optimal threshold segmentation method, as most of the narrow streams are well extracted in each of the study areas and had a comparable accuracy to our method. However, it should also be noted that this method is not an automatic method, since many parameters need to be tuned, which prevents it from being effectively used in larger areas. Overall, our method outperforms all others in terms of measurements except for producer accuracy in study area 3, where Yang’s method achieves a relatively higher accuracy.
4. Extraction of inland water of Fujian province
The aforementioned experiments demonstrate that the MNWI is the most efficient algorithm for narrow water extraction, but it is still interest to address whether it is applicable to large-volume data in an actual scenario. Therefore, we apply our method to extract inland water features in Fujian province, China, which is a relatively large region that covers an area of about 121,000 km2. Fujian is a mountainous province, located on the southeast coast of China and facing Taiwan across the Taiwan Strait. It has significant vegetation cover because of high mean precipitation and warm annual temperatures. Topographically, Fujian is a very mountainous region, having abundant water resources, rivers, lakes, and reservoirs. Many rivers run through these mountains, of which the most important is the Min River, as its drainage area covers over 50% of the province. The upstream Jin River, Futun River, and Shaowu River all converge into the Min River. The Jiulong River flows south of the Min River, reaching the sea at Xiamen city, and the Ting River runs across Fujian’s southwestern border. It is thus an appropriate region to test our methods in a large area. To cover the entirety of Fujian province, 13 Landsat 8 OLI images were collected. Fujian is usually cloudy and rainy in spring and summer, so we collected all the images in winter to avoid cloud interference. The acquisition information is summarized in Table 2. Note that the quality of all acquired data is relatively good, with cloud cover less than 10%.
|Acquisition time||November 17||December 3||October 23||October 23||October 23||December 1||December 1|
|Acquisition time||December 1||December 1||December 1||October 5||October 5||October 5|
All the images were matched and stitched without altering their spectral color (Figure 9, left), and the final inland water information for the Fujian province (Figure 9, right) shows that 2,494,988 pixels were classified into inland surface water. It can be calculated that the total inland water area of Fujian province is about 2245.49 km2 in the winter of 2013. Visually, our method can extract the most of perceptible water bodies with a high accuracy, where the main rivers, such as Min River, Jiurong River, Ting River, etc., are all correctly delineated with clear river boundaries. Moreover, most of the spatial shapes of small tributaries were continuous and so were the lakes and reservoirs. To evaluate the quantitative classification accuracy, 6800 samples of water and non-water features were randomly selected. The water samples included 3340 pixels, of which wide water to narrow water ratio was about 2:1, as the non-water samples were 3460 pixels, chosen from possibly confused land cover types such as forest, built-up areas, and rare soil. The producer accuracy, user accuracy, overall accuracy, and kappa coefficient calculated from these samples are 94.33%, 98.7%, 96.61%, and 0.932, respectively, indicating that the proposed method can achieve a high accuracy for water extraction in a large and complex area and it is an effective optional tool for practical water extraction.
Most water indices can perform well for the extraction of wide water features from remote sensing images, but they are normally ineffective in the extraction of narrow water features. This chapter has described a new method using a morphological top-hat transform to form a novel narrow water index, denoted as MNWI, which improves the extraction accuracy of narrow water from Landsat images. Experimental results demonstrated impressive performances of our method of the narrow water extraction. A case study in Fujian province suggests that it is an effective and practical tool for large area inland water.
Funding was provided by the Natural Science Foundation of China (Grant No. 41571330, 41830108).
Sun F, Sun W, Chen J, Gong P. Comparison and improvement of methods for identifying water bodies in remotely sensed imagery. International Journal of Remote Sensing. 2012; 33:6854-6875
Frazier PS, Page KJ. Water body detection and delineation with Landsat TM data. Photogrammetric Engineering and Remote Sensing. 2000; 66:1461-1468
Rokni K, Ahmad A, Solaimani K, Hazini S. A new approach for surface water change detection: Integration of pixel level image fusion and image classification techniques. International Journal of Applied Earth Observation and Geoinformation. 2015; 34:226-234
McFeeters S. The use of the normalized difference water index (NDWI) in the delineation of open water features. International Journal of Remote Sensing. 1996; 17:1425-1432
Xu H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. International Journal of Remote Sensing. 2006; 27:3025-3033
Carleer AP, Wolff E. Urban land cover multi-level region-based classification of VHR data by selecting relevant features. International Journal of Remote Sensing. 2006; 27(6):1035-1051
Yang S, Xue S, Liu T. Method for automatically extracting fine water body by using TM image. Journal of Surveying and Mapping (in Chinese with English Abstract). 2010; 39(6):611-617
Li Y, Ding J, Yan R. Study on extraction method of small water bodies in mountainous area based on GF-1 remote sensing images. Resource Science (in Chinese with English Abstract). 2015; 37(2):408-416
Jiang H, Feng M, Zhu Y. An automated method for extracting rivers and lakes from Landsat imagery. Remote Sensing. 2014; 6(6):5067-5089
Beneditsson JA, Pesaresi M, Arnason K. Classification and feature extraction for remote sensing images from urban area based on morphological transformations. IEEE Transaction on Geosciences and Remote Sensing. 2003; 41(9):1940-1949
Huang X, Zhang L. A multidirectional and multi-scale morphological index for automatic building extraction from multispectral GeoEye-1 imagery. Photogrammetric Engineering & Remote Sensing. 2011; 77(7):721-732
Wu B, Fang C, Yu L, Huang X, Zhang Q . A fully automatic method to extract rare earth mining area from Landsat image. Photogrammetric Engineering and Remote Sensing. 2016; 82(9):55-64
Huang X, Zhang L, Wang L. Evaluation of morphological texture features for mangrove forest mapping and species discrimination using multispectral IKONOS imagery. IEEE Geoscience and Remote Sensing Letters. 2009; 6(3):393-397
Ji L, Zhang L, Wylie B. Analysis of dynamic thresholds for the normalized difference water index. Photogrammetric Engineering and Remote Sensing. 2009; 75(11):1307-1317
Otsu N. A threshold selection method from gray-level histograms. IEEE Transactions on Systems, Man and Cybernetics. 1979; 9(1):62-66
Zha Y, Gao Y, Ni S. Use of normalized difference built-up index in automatically mapping urban areas from TM imagery. International Journal of Remote Sensing. 2003; 24(3):583-594