Next Article in Journal
Recent Deceleration of the Ice Elevation Change of Ecology Glacier (King George Island, Antarctica)
Next Article in Special Issue
Quantifying Streamflow Variations in Ungauged Lake Basins by Integrating Remote Sensing and Water Balance Modelling: A Case Study of the Erdos Larus relictus National Nature Reserve, China
Previous Article in Journal
Satellite-Based Sea Ice Navigation for Prydz Bay, East Antarctica
Previous Article in Special Issue
Estimating Snow Mass and Peak River Flows for the Mackenzie River Basin Using GRACE Satellite Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Monitoring of the Largest Freshwater Lake in China Using a New Water Index Derived from High Spatiotemporal Resolution Sentinel-1A Data

1
The State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, P.O. Box 9718, Datun Road, Chaoyang, Beijing 100101, China
2
College of Resource and Environment, University of Chinese Academy of Sciences, Yuquan Road 19, Shijingshan, Beijing 100049, China
3
College of Environment and Planning, Henan University, Jinmingdadao Road, Kaifeng 475004, China
*
Authors to whom correspondence should be addressed.
Remote Sens. 2017, 9(6), 521; https://doi.org/10.3390/rs9060521
Submission received: 20 February 2017 / Revised: 18 May 2017 / Accepted: 19 May 2017 / Published: 24 May 2017
(This article belongs to the Special Issue Remote Sensing of Climate Change and Water Resources)

Abstract

:
Poyang Lake is the largest freshwater lake in China and is well known for its ecological function and economic importance. However, due to the influence of clouds, it is difficult to dynamically monitor the changes in water surface areas using optical remote sensing. To address this problem, we propose a novel method to monitor these changes using Sentinel-1A data. First, the Sentinel-1A water index (SWI) was built using a linear model and a stepwise multiple regression analysis method with Sentinel-1A and Landsat-8 imagery acquired on the same day. Second, water surface areas of Poyang Lake from 24 May 2015 to 14 November 2016 were extracted by the threshold method utilizing time-series SWI data with an interval of 12 days. The results showed that the SWI threshold classification method could be applied to different regions during different periods with high quantity accuracy (approximately 99%). The water surface areas ranged between 1726.73 km2 and 3729.19 km2 during the study periods, indicating an extreme variability in the short term. The maximum and average values of the changed areas were 875.57 km2 (with a change rate of 35%) and 197.58 km2 (with a change rate of 8.2%), respectively, after 12 days. The changes in the mid-western region of Poyang Lake were more dramatic. These results provide baseline data for high-frequency monitoring of the ecological environment and wetland management in Poyang Lake.

Graphical Abstract

1. Introduction

Lakes are important water and wetland resources. Lake water surface areas are sensitive to both climate change and human activities and, therefore, serve as an excellent indicator to assess local ecosystem changes [1,2,3,4,5]. Tremendous changes in lake water surface areas have occurred in various lakes worldwide due to temperature increases, precipitation variations, and so on, over the last 50 years [1,6,7,8,9]. In addition, monitoring of these changes in lake water surface areas is the first phase in further research related to ecosystem services [10], as well as wetland [11] and biodiversity assessments [12] around lakes.
Satellite remote sensing is one of the most efficient and widely used methods for the dynamic monitoring of water surface areas due to its fast response, wide field of view, periodic observation, and low cost [6,13,14,15,16,17,18]. Optical remote sensing is a frequently used dataset for dynamic monitoring of changes in water surface areas due to its availability. Changes in Poyang Lake have been monitored via the moderate-resolution imaging method of spectroradiometer (MODIS) imagery at monthly and yearly intervals from 2000 to 2010 [6]. However, MODIS spatial resolution (250 m/500 m) cannot delineate the detailed outline of lakes because of the mixed pixel problem [6]. Landsat datasets have been used to monitor changes in lake distributions on the Mongolian Plateau [2]. However, due to the influence of clouds, only inter-annual changes were observed. It is difficult for optical remote sensing satellites to provide enough data for dynamic monitoring of changes in water surface areas of Poyang Lake. This is especially true for medium or high spatial resolution remote sensing, which has a longer return cycle because the Poyang Lake water surface areas experience rapid and significant short-term (one or two weeks) fluctuations. Although several methods have been introduced to fill the gaps in medium or high spatial resolution remote sensing data by combination with high temporal remote sensing data, these methods cannot fill all the gaps because of the issues related to cloud cover [19,20,21,22].
Conversely, synthetic aperture radar (SAR) remote sensing has the advantages of being an all-weather and all-time tool [23,24]. In particular, the open-access Sentinel-1A SAR imagery (S1A) has high spatiotemporal resolution (return cycle of 12 days and imagery pixel spacing of 10 m) [25] and great potential for the monitoring of these dynamic changes in water surface areas. Cazals et al. [26] detected the hydrological dynamics of a coastal marsh located in the Regional Natural Park (close to the French Atlantic coastline) using the threshold method with S1A data, producing an overall accuracy of 82%. Radarsat-2 and ENVISAT imagery were also utilized to determine the extent of flooding in the upper Zambezi basin in southern Africa using the threshold method with an estimated error of 8–23% compared to the classification results of Landsat-8 [24]. However, it is a difficult and complex task to determine appropriate thresholds in different periods and regions since the applicability of the threshold is unknown. Moreover, the accuracy of the above study does not meet the needs of investigating dynamic changes in Poyang Lake. Therefore, a novel threshold value classification method with high accuracy is required.
As a spectral water index from optical imagery, the modified normalized difference water index (MNDWI) is an effective way to extract a water body [2,27,28,29,30]. Generally, a pixel is classified into water if its MNDWI value is more than a certain threshold value [2,30]. Therefore, the technique process of this threshold-categorized method is simple and efficient. However, the MNDWI value is based on optical remote sensing data that are easily influenced by clouds. Therefore, there is a need to monitor water surface areas dynamically using an easy-to-operate index, such as MNDWI, based on SAR data. To address the problems mentioned above, the aims of the present study were: (1) to propose a MNDWI-like index based on S1A data for dynamic monitoring of water surface areas; and (2) to monitor the water areas of Poyang Lake using the MNDWI-like index based on S1A data.

2. Study Area

Poyang Lake, the largest freshwater lake in China, was selected as the study area. Poyang Lake is located in the northern Jiangxi Province (29°05′–29°17′N and 115°54′–116°08′E, southeastern China) and outflows into the Yangtze River at Hukou. The average water surface area of Poyang Lake is more than 2500 km2 during the wet season [6]. The lake is the primary component of the Poyang Lake wetland and has been included in the first set of the Ramsar Convention List of Wetlands of International Importance [7]. The Poyang Lake region experiences a subtropical monsoon climate that is characterized by a rainy season (April to September) and a dry season (October to March) [6]. The annual average temperature is 25 °C, and the annual average accumulative precipitation is approximately 1600 mm [31]. The main water sources into the Poyang Lake are local precipitation and rivers (including the Xiushui, Ganjiang, Fuhe, Xinjiang, and Raohe rivers, Figure 1) [6].
Numerous small water bodies in the periphery have been disconnected from the large lake system because more than 500 levees have been built in Poyang Lake to protect approximately 5082 km2 of farmland and settlements [32,33]. In addition, the boundary of the water surface of Poyang Lake is unstable and needs to be defined. The outline of the maximum water surface areas (on 17 July 2016) was defined as being approximately 4105 km2 based on a previous study (Figure 1) [6].
Another minor study area is Taihu Lake, which is located in Jiangsu Province, China (30°56′–31°33′N and 119°52′–120°36′E) and is one of China’s five major freshwater lakes. This area was used to examine the results of the SWI-fitted equation and to examine the classification accuracy of S1A.

3. Materials and Methodology

3.1. Datasets and Preprocessing

Datasets utilized for the present study included S1A images (https://scihub.copernicus.eu/), and Landsat-8 Operational Land Imager images (OLI, http://earthexplorer.usgs.gov/).
Sentinel-1 is a polar orbiting two-satellite constellation (S1A and -B) that was launched by the European Space Agency (ESA) [34]. S1A and -B, the next generation of C-band (center frequency 5.405 GHz) radar missions, were launched on 3 April 2014 and 25 April 2016, respectively, from Europe’s Spaceport in French Guiana. The revisit time of a single satellite is 12 days [25], whereas the two-satellite constellation provides a revisit time of 6 days [35]. The S1A utilized in the present study is a SAR standard L1 product, and belongs to the Interferometric Wide-swath mode (IW) with dual-polarization (VV/VH) [26]. This data product has a wide coverage of 250 km with medium resolutions of 5 m and 20 m in range and azimuth directions, respectively, with an equivalent number of looks of five, and imagery pixel spacing of 10 m in ground geometry [25]. For detailed descriptions regarding S1A, refer to reference [25].
In the present study, 71 S1A scenes from 41 observations by Sentinel-1A satellite were selected from 24 May 2015 to 14 November 2016. Preprocessing of S1A imagery was performed with SNAP software provided by ESA (https://sentinel.esa.int/web/sentinel/toolboxes/sentinel-1). The workflow included five steps. First, S1A images were allocated orbit files. Second, S1A images were radiometrically calibrated to output radar backscatter bands (σ°). Third, radar backscatter bands were orthorectified using the Range Doppler Terrain Correction algorithm with SRTM DEM and spatial resolution of 1 s. Fourth, the backscattering coefficient (in dB) was acquired from the orthorectified radar backscatter band by the equation 10 × log10 (σ°) [36]. Fifth, a median filter with window size of 5 × 5 pixel was utilized to remove speckle noises [23,37,38].
A Landsat-8 satellite equipped with two sensor payloads (OLI and Thermal Infrared Sensor) was launched on 11 February 2013. In the present study, one OLI scene (acquired on 27 September 2016) was selected as the training data and applied for building the regression equation, and two OLI scenes (acquired on 9 September 2015 and 2 February 2016) were selected as testing data and applied to assess the classification accuracy of S1A. Eight OLI scenes (acquired on: 9 September 2015; 11 October 2015; 30 December 2015; 16 February 2016; 3 March 2016; 23 June 2016; 25 July 2016; and 27 September 2016) were used to extract the water surface areas of Poyang Lake and were compared with the water surface areas of Poyang Lake extracted from S1A from 24 May 2015 to 14 November 2016. Preprocessing work included atmospheric correction using the Fast Line-of-sight Atmospheric Analysis of Spectral Hypercube module, geometric correction, and data clipping based on the operational boundary of Poyang Lake. At the end, all images (S1A and OLI) were re-projected to an Asia North Albers Equivalent Conical Projection with the first standard parallel 25°N, the second standard parallel 47°N, and the central meridian 105°E.

3.2. Sentinel-1A Water Index

Water bodies create smooth and homogenous surfaces and, thereby, the return of radar pulses will be weak and the radar imagery over the area will be black or very dark in contrast with the surrounding non-water areas that produce a stronger return and brighter color [39]. Optical water indexes, such as MNDWI, have the opposite characteristics. In MNDWI imagery, water bodies appear as brighter colors and non-water areas appear as darker colors. Thus, there is a negative correlation between S1A imagery and MNDWI imagery. The S1A water index (SWI) was built based on the linear regression equation between S1A (including VH and VV polarization) and MNDWI. Simultaneously, SWI has a characteristic similar to MNDWI computed from OLI based on the flowing equation [30]:
M N D W I = ( ρ g r e e n ρ s w i r ) / ( ρ g r e e n + ρ s w i r )
where ρ g r e e n represents the reflectance of the green band (B3 of OLI) and ρ s w i r represents the reflectance of the shortwave infrared band (B6 of OLI).
To match S1A pixels with MNDWI imagery pixels, the spatial resolution of S1A was resampled to 30 m (because the spatial resolution of MNDWI imagery is 30 m) using the nearest-neighbor method. Then, the MNDWI imagery and S1A (with spatial resolution of 30 m) were stacked using the stacking layer function in ENVI (version 4.8) software. Thus, a pixel on the stacked imagery had three attributes (e.g., VH (the abbreviation for the backscattering coefficient in VH polarization and VV has the same meaning), VV, and MNDWI values) and then a regression analysis was performed.
A total of 44,284 pairs of data were selected evenly and randomly from the stacked imagery to build the SWI. OLI and S1A were acquired on the same day (27 September 2016), which helped reduce the error from changes in water surface areas for different imaging dates when the SWI was built. These data cover the water body, grassland, forestland, farmland, bare land, and built-up land, with 42.31% being the water body. MNDWI and S1A were defined as dependent variables and independent variables, respectively. The SWI was built using the stepwise regression analysis method based on MNDWI, and six independent variables (VH, VV, VH × VV, VH2, VV2, and VH − VV) were derived from S1A.

3.3. Index Assessment

Several parameters including the root mean square error (RMSE), mean absolute error (MAE), root mean absolute percentage error (RMAPE), and R2 were used to assess the SWI. These indexes are regularly employed in model evaluation studies and their expressions are as follows [40,41,42]:
δ R M S E = 1 n i = 1 n ( y i y ^ i ) 2
δ M A E = 1 n i = 1 n | y i y ^ i |
δ R M A P E = 1 n i = 1 n [ ( y i y ^ i ) / y i ] 2
R 2 = ( i = 1 n ( y ^ i y ^ ¯ ) × ( y i y ¯ ) / ( i = 1 n ( y ^ i y ^ ¯ ) 2 × i = 1 n ( y i y ¯ ) 2 ) ) 2
where n is the total number of samples; y i represents the observation value; y ^ i represents the estimated value; y ¯ represents the average observation value; and y ^ ¯ represents the average estimated value. When RMSE, MAE, and RMAPE are low and R2 is high, then the best-fit equation is achieved.
To examine the applicability of the SWI-fitted equation in different periods and regions, 190 and 225 pairs of testing data (SWI and MNDWI, respectively) were selected randomly from S1A and OLI images acquired on 2 February 2016, in the Taihu Lake region, China, and 9 September 2015, in the Poyang Lake region, respectively, to examine the fitted SWI using the indexes mentioned above.

3.4. Extraction of the Water Area

The water area extracted from OLI was regarded as the ground-truth water area and was extracted using the threshold method. If the MNDWI value of a pixel is more than 0 and the normalized difference vegetation index (NDVI) value is less than 0.15, it is classified as water. The NDVI was computed from OLI based on the flowing equation:
N D V I = ( ρ n i r ρ r e d ) / ( ρ n i r + ρ r e d )
where ρ n i r represents the reflectance of the near-infrared band (B5 of OLI) and ρ r e d represents the reflectance of the red band (B4 of OLI).
The histogram frequency of the SWI imagery (on 27 September 2016) based on the ground-truth water and non-water is plotted in Figure 2 to facilitate visualization. Figure 2 demonstrates that the majority of the SWI water values are more than 0.40 and the majority of the SWI non-water values are less than 0.18 through visual judgment. However, an appropriate and unique threshold ranging from 0.18 to 0.40 needed to be determined by calculating the positive and negative errors. The threshold is regarded as appropriate when the classification errors of water and non-water are approximate and the overall accuracy derived from the confusion matrix is optimal (the sum of positive and negative errors comes to the minimal value). The result showed that the sum is minimal (0.02%) when the threshold is 0.20. In the present study, the water information of Poyang Lake was extracted using the threshold method from 24 May 2015 to 14 November 2016.

3.5. Evaluation of Quantity Accuracy

The water area extracted from OLI was regarded as the ground-truth water area [24] to examine the classification accuracy of S1A in three steps. First, the spatial resolution of S1A was resized to 30 m using the nearest-neighbor method to match the spatial resolution of OLI. Then, three periods of ground-truth water areas from OLI (two OLI scenes acquired on 9 September 2015 and 27 September 2016 in Poyang Lake region, and one OLI scene acquired on 2 February 2016 in the Taihu Lake region) were selected to assess the accuracy of the SWI category during different periods and in different regions using the confusion matrix method. Finally, quantity accuracy was computed with the following equation:
ζ = ( 1 - | A e x t r a c t e d A t r u t h | / A t r u t h ) × 100
where ζ represents the quantity accuracy; Aextracted represents the water surface area extracted from SWI imagery; and Atruth represents the ground-truth water area extracted from the OLI bands.

3.6. Dynamic Change Analysis Method

To understand the characteristics and influences of the water changes in Poyang Lake, the submergence ratio was computed as the following equation for each pixel:
δ r a t i o = ( n / N ) × 100 %
where δratio represents the submergence ratio; n represents the number of inundations observed by S1A; and N represents the total observations times of S1A (i.e., 41 in this study). The submergence ratio approximates the probability of being submerged for each pixel.
Poyang Lake is composed of a number of lakes with different areas. Lakes were grouped to facilitate the statistical change detection analyses into the following classes: 0.1–1.0 km2, >1.0–10 km2, and >10–100 km2. The change in the number of lakes with different areas was analyzed to explain the dynamic change in the Poyang Lake water surface areas in this study.

3.7. Extraction of Water Surface Areas from SWI, VH, and VV Images Using Different Methods

In order to prove that the SWI image is better than the original image of S1A (including the VH image and VV image) to extract water, four classification methods were used, i.e., the threshold classification method (TCM), IsoData unsupervised classification (IUC), Minimum Distance Classification (MDC), and Support Vector Machine (SVM). The test zone was a rectangle with an area of 110 km × 70 km location in the Poyang Lake region.
Similar to Figure 2 and its theory, the threshold values were confirmed as −23 dB and −17 dB based on VH and VV imagery, respectively. For a single band, a pixel was classified as water if its value was less than −23 dB in VH imagery or −17 dB in VV imagery, respectively. For double bands, a pixel was classified as water if its value was less than −23 dB in VH imagery and −17 dB in VV imagery.
IUC is barely affected by subjectivity [43,44] and the classification results utilizing IUC demonstrate the ability of the imagery to distinguish between water and non-water areas. Water was extracted from the original image of S1A and SWI image using the same default parameters of IUC model in ENVI software.
MDC and SVM are supervised classification methods that are used in many studies [23,45,46]. The resamples (1937 pixels) used in MDC and SVM were randomly selected based on the ground-truth water extracted from OLI. The resamples are the same when extracting water from the original image and the SWI image using MDC and SVM.

3.8. Extraction of Water Surface Areas from MODIS Data

To understand the ability to identify the details of the water, the classification accuracy of water extracted from the MODIS image was evaluated based on the ground-truth water extracted from OLI. The MODIS image was acquired on 27 September 2016. The green band (B4 of MODIS, center wavelength 555 nm) and the shortwave infrared band (B6 of MODIS, center wavelength 1640 nm) were used to compute the MNDWI according to Equation (1). The red band (B1 of MODIS, center wavelength 645 nm) and the near-infrared band (B2 of MODIS, center wavelength 860 nm) were used to compute the NDVI according to Equation (6). If the MNDWI value of a pixel is more than 0 and the NDVI value of a pixel is less than 0.15, it is classified as water. Thus, the water image form MODIS was obtained. Then, the water image form MODIS was resized to 30 m using the nearest-neighbor method when building the accuracy assessment confusion matrix.

4. Results

4.1. The SWI-Fitted Equation

VV, VH × VV, VV2, and VH2 were selected with the stepwise regression analysis method. The SWI fitted equation was as follows:
S W I = 0.1747 × β v v + 0.0082 × β v h × β v v + 0.0023 × β v v 2 0.0015 × β v h 2 + 0.1904
where β v h and β v v represent the backscattering coefficient in VH polarization and VV polarization, respectively.
The R2 of the SWI fitted equation was 0.83, the RMSE was 0.23, the F-value was 31,874, and the p-value was less than 0.01. Therefore, the SWI fitted equation was statistically significant.
Table 1 shows the examination results of the SWI fitted equation in different periods and regions. The data in Table 1 show that the SWI-fitted equation was statistically significant.

4.2. Water Area Extraction Using SWI

The best classification results were obtained using the threshold of 0.2 to distinguish water and non-water from the SWI images. Figure 3 and Figure 4 show the extracted water surface areas of Poyang Lake from 24 May 2015 to 14 November 2016, with an interval of 12 days. The average, maximum, and minimum areas of Poyang Lake were 2638.36 km2, 3729.19 km2 (on 17 July 2016), and 1726.73 km2 (on 7 March 2016), respectively.
To better understand the dynamic changes of the water surface areas of Poyang Lake, the changed area and the variation rate of the water surface areas in 12-day intervals were analyzed. During the period of analysis, there were 22 changes of more than 100 km2. The maximum changed area was 875.57 km2 with a variation rate of 35% from 12 April to 24 April 2016. The average changed area was 197.58 km2. The variation rate in each 12-day interval showed that there were 11 changes of >10%. The average variation rate was 8.2%. The water surface areas of Poyang Lake had several significant changes within a 12-day window over the course of the study period.

4.3. Accuracy Evolution

The classification accuracy of SWI is presented in Table 2 and plotted in Figure 5 to facilitate visualization. In the Poyang Lake region, the overall accuracy was more than 96%. On 27 September 2016, the producer’s accuracy and user’s accuracy of water were 91.51% and 92.58%, respectively. On 9 September 2015, the producer’s accuracy and user’s accuracy of water were 94.18% and 93.41%, respectively. The approximation between producer’s accuracy and user’s accuracy leads to higher quantity accuracy than it does overall accuracy because the rates of omission and commission counteract. The quantity accuracy of water was 99.18% and 98.84% on 9 September 2015 and 27 September 2016, respectively. The accuracy was also superior in the Taihu Lake region. Therefore, the SWI threshold classification method could be applied to different regions during different periods with high and reliable accuracy.

4.4. Submergence Ratio Analysis

The submergence ratio can document the spatiotemporal dynamic change of the water surface areas of Poyang Lake and is plotted in Figure 6. The area that was submerged 100% of the time was 1307.41 km2, or 31.85% of the entirety of Poyang Lake; 876.07 km2, or 29.04% was submerged between 80% and 100%; and 1605.61 km2, or 39.11%, was submerged for less than 60% of the study period.
The change characteristics of Poyang Lake had obvious regional differences. To expediently describe the change characteristics in different regions, Poyang Lake was roughly divided into four regions based on submergence ratios (Figure 6, Regions I–IV). As a whole, the submergence ratio was more than 90% for the majority of Regions I and IV, less than 40% for the majority of Region III, and Region II had a complex submergence ratio. In Region II, there were many isolated lakes whose submergence ratios were more than 90% and had an area of 430.36 km2 (27% of the total area in Region II), and there were many districts whose submergence ratios were less than 40% with an area of 717.88 km2 (44% of the total area in Region II). During this study period, the variation rate in the changed area of Region II was the greatest among the four regions. The maximum variation rate was 80% between 12 April and 24 April 2016, with an average variation rate of 13% in Region II.

4.5. The Number of Lake with Different Area

Poyang Lake is composed of a number of lakes (Figure 7). The number of lakes with an area of 0.1–1 km2 ranged from 166 to 465 (average 336) in this study. The number of lakes with an area of >1–10 km2 ranged from 30 to 121 (average 78). The number of lakes with an area of >10–100 km2 ranged from 8 to 27 (average 15).
During the rainy season, from April to September, the water surface areas can expand to >3000 km2 and form several big lakes, and during the dry season, from October to March, the water surface areas can shrink to <2000 km2 and a number of small and isolated lakes are then segregated from the big lakes. The number of lakes of 0.1–1 km2 drastically fluctuated. During the dry season, the average, maximum, and minimum number of lakes was approximately 396, 465, and 320, respectively. During the rainy season, the average, maximum, and minimum number of lakes was approximately 232, 338, and 166, respectively. The change in the number of lakes of 0.1–1 km2 in 12-day intervals was analyzed. The maximum change in number was 157 in the 12-day interval from April 12 to April 24, 2016. The average change in number was 31, which means that 31 lakes of 0.1–1 km2 disappeared or appeared within 12-day intervals.

4.6. Comparison of Water Area Extracted Using Different Method and Different Data

Table 3 shows that the quantity accuracy and overall accuracy of water extracted from SWI imagery was higher than that of the original imagery of S1A using the four methods based on the ground-truth water area that was extracted from OLI acquired on 27 September 2016.
The water surface areas of Poyang Lake extracted from OLI and SWI are plotted in Figure 8. There were only eight OLI scenes that were cloud-free from 24 May 2015 to 14 November 2016 (Figure 8), and many variation characteristics of the water were not detected by OLI compared with the results from the SWI images.
Water extracted from MODIS on 27 September 2016 is shown in Figure 9b. The classification accuracy was assessed based on the ground-truth water area from OLI. The quantity accuracy was 67.56% (omission 537 km2). The overall accuracy was 88.96% (kappa coefficient 0.64) and the producer’s accuracy was only 58.71%, which indicates that 41.29% of the water was omitted (see Figure 9c). The quantity accuracy and overall accuracy from the SWI image were 31.28% and 7.56% higher, respectively, than that from the MODIS image.

5. Discussion

The edges of the water body are the main zones where incorrect classification occurred since these environments are very complex (e.g., there were widely distributed shallows, sand, and plankton) in the present study (Figure 5). The zones of the commission of water mainly appeared as a light-green and bright white color in OLI (Figure 10a). There were sparse emergent aquatic plants or plankton in the light-green zones, whose spectral characteristics were similar to the vegetation (Figure 11). Thus, these zones might be categorized as vegetation in the OLI images. However, in these zones, the main return of the radar pulse came from the water due to the low density of vegetation. Thus, these zones were categorized as water in the SWI images, which should be regarded as correct. The bright white zones were sand with high water content and low surface roughness; therefore, these bright white zones were wrongly categorized as water in the SWI image.
We cannot determine whether the omission water (Figure 10b) was part of the water body because its spectral characteristics were different from those of normal water (Figure 11). However, they were categorized as water in the MNDWI image. Therefore, further study is required to investigate this.
Compared with the original images of S1A, SWI images have the ability to enhance the difference between water and non-water. Furthermore, high speed, efficiency, and classification accuracy are regarded as being advantages of the SWI threshold classification method for extracting water from SWI imagery, as shown in Table 3.
Table 3 also shows that the area of water can be extracted with the highest QA of 99.57% using the SVM method with SWI data. However, when compared to the similarly performing TCM method (QA: 98.84%, see Table 3), it is more time consuming to analyze inundation extent using SVM.
Compared to the OLI and MODIS images, SWI images have high effective spatiotemporal resolution. The revisit time of OLI is 16 days, and the frequency of cloudy or rainy weather is very high in the south of China. Therefore, there is a lack of useful OLI data when detecting the dynamic changes in Poyang Lake (Figure 8).
Although MODIS has a revisit time of one day, its spatial resolution is too coarse to identify the details of the water and small water bodies. Poyang Lake is composed of a number of lakes with areas of less than 1 km2 that are easily omitted from MODIS. The annual water surface area of Poyang Lake extracted from MODIS ranged from 700 km2 to 2700 km2 during 2000 to 2010 [6]; however, the water surface area ranged from 1700 km2 to 3700 km2 in the present study. Additionally, the quantity accuracy from MODIS was 67.56%, and meaning that water surface area of 537 km2 wasn’t detected by MODIS in sub-region in the Figure 9. It may be possible to improve the quantity accuracy by other advanced methods, although we still believe that MODIS cannot precisely identify details of the water and small water bodies, especially in areas where water changes are severe.

6. Conclusions

In the present study, a new water index (SWI) was constructed based on S1A and OLI-MNDWI using the multiple stepwise regression analysis method. The best classification results were obtained using the threshold of 0.2 to distinguish water and non-water from the SWI images. The advantages of extracting water using the SWI are that it is quick and more efficient than the other Machine Learning methods. In addition, the SWI threshold classification method can be applied to different regions during different periods with high quantity accuracy (approximately 99%) and overall accuracy (approximately 97%).
A total of 41 days of water surface areas for Poyang Lake were extracted from 24 May 2015 to 14 November 2016, with an interval of 12 days. The average area of Poyang Lake was 2638.36 km2 during the study period. The maximum and minimum areas were 3729.19 km2 (on 17 July 2016) and 1726.73 km2 (on 7 March 2016), respectively. The results showed that changes in water surface areas were dramatic within periods of 12 days. The maximum and average changed areas were 875.57 km2 (with a change rate of 35% from 12 April to 24 April 2016) and 197.58 km2 (with a change rate of 8.2%) in 12 days, respectively. Moreover, the change in the mid-western region of Poyang Lake was more marked, with a maximum change rate of 80% between 12 April and 24 April 2016, and an average change rate of 13%.
Poyang Lake is composed of a number of lakes, with the number drastically fluctuating with the changes in water surface area. The number of lakes with areas of 0.1–1 km2, >1–10 km2, and >10–100 km2 ranged from 166 to 465, from 30 to 121, and from 8 to 27, respectively, in the present study.
These results provide scientific guidance and baseline data for high-frequency monitoring of the ecological environment and wetland management in Poyang Lake.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (41301390 and U1404401) and the Youth Innovation Promotion Association CAS (2017089).

Author Contributions

Haifeng Tian, Mingquan Wu, and Zheng Niu conceived and designed the methodology; Haifeng Tian and Wang Li performed the methodology; Ni Huang, Guodong Li, and Xiang Li analyzed the data; and Haifeng Tian wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ma, R.H.; Duan, H.T.; Hu, C.M.; Feng, X.Z.; Li, A.N.; Ju, W.M.; Jiang, J.H.; Yang, G.S. A half-century of changes in China’s lakes: Global warming or human influence? Geophys. Res. Lett. 2010, 37, 6. [Google Scholar] [CrossRef]
  2. Tao, S.L.; Fang, J.Y.; Zhao, X.; Zhao, S.Q.; Shen, H.H.; Hu, H.F.; Tang, Z.Y.; Wang, Z.H.; Guo, Q.H. Rapid loss of lakes on the mongolian plateau. Proc. Natl. Acad. Sci. USA 2015, 112, 2281–2286. [Google Scholar] [CrossRef] [PubMed]
  3. Hui, J.; Yao, L. Analysis and inversion of the nutritional status of China’s Poyang Lake using MODIS data. J. Indian Soc. Remote Sens. 2016, 44, 837–842. [Google Scholar] [CrossRef]
  4. Yang, P.; Liu, X.H.; Xu, B. Spatiotemporal pattern of bird habitats in the Poyang Lake based on landsat images. Environ. Earth Sci. 2016, 75, 8. [Google Scholar] [CrossRef]
  5. Tang, X.G.; Li, H.P.; Xu, X.B.; Yang, G.S.; Liu, G.H.; Li, X.Y.; Chen, D.Q. Changing land use and its impact on the habitat suitability for wintering anseriformes in China’s Poyang Lake region. Sci. Total Environ. 2016, 557, 296–306. [Google Scholar] [CrossRef] [PubMed]
  6. Feng, L.; Hu, C.M.; Chen, X.L.; Cai, X.B.; Tian, L.Q.; Gan, W.X. Assessment of inundation changes of Poyang Lake using MODIS observations between 2000 and 2010. Remote Sens. Environ. 2012, 121, 80–92. [Google Scholar] [CrossRef]
  7. Wu, G.P.; Liu, Y.B. Capturing variations in inundation with satellite remote sensing in a morphologically complex, large lake. J. Hydrol. 2015, 523, 14–23. [Google Scholar] [CrossRef]
  8. Yu, G.; Shen, H.D. Lake water changes in response to climate change in northern China: Simulations and uncertainty analysis. Quat. Int. 2010, 212, 44–56. [Google Scholar] [CrossRef]
  9. Bryant, M.D. Global climate change and potential effects on pacific salmonids in freshwater ecosystems of southeast alaska. Clim. Chang. 2009, 95, 169–193. [Google Scholar] [CrossRef]
  10. Dronova, I.; Gong, P.; Wang, L. Object-based analysis and change detection of major wetland cover types and their classification uncertainty during the low water period at Poyang Lake, China. Remote Sens. Environ. 2011, 115, 3220–3236. [Google Scholar] [CrossRef]
  11. Sica, Y.V.; Quintana, R.D.; Radeloff, V.C.; Gavier-Pizarro, G.I. Wetland loss due to land use change in the lower parana river delta, argentina. Sci. Total Environ. 2016, 568, 967–978. [Google Scholar] [CrossRef] [PubMed]
  12. Feng, L.; Han, X.; Hu, C.; Chen, X. Four decades of wetland changes of the largest freshwater lake in China: Possible linkage to the three gorges dam? Remote Sens. Environ. 2016, 176, 43–55. [Google Scholar] [CrossRef]
  13. Kutser, T.; Pierson, D.C.; Kallio, K.Y.; Reinart, A.; Sobek, S. Mapping lake cdom by satellite remote sensing. Remote Sens. Environ. 2005, 94, 535–540. [Google Scholar] [CrossRef]
  14. Kallio, K. Remote sensing monitors lakes. Trac Trends Anal. Chem. 2012, 34, 3–4. [Google Scholar]
  15. Giardino, C.; Bresciani, M.; Stroppiana, D.; Oggioni, A.; Morabito, G. Optical remote sensing of lakes: An overview on lake maggiore. J. Limnol. 2014, 73, 201–214. [Google Scholar] [CrossRef]
  16. Wright, R.; Pilger, E. Satellite observations reveal little inter-annual variability in the radiant flux from the mount erebus lava lake. J. Volcanol. Geotherm. Res. 2008, 177, 687–694. [Google Scholar] [CrossRef]
  17. Veettil, B.K.; Bianchini, N.; de Andrade, A.M.; Bremer, U.F.; Simoes, J.C.; de Souza, E. Glacier changes and related glacial lake expansion in the bhutan himalaya, 1990–2010. Reg. Environ. Chang. 2016, 16, 1267–1278. [Google Scholar] [CrossRef]
  18. Hereher, M.E. Environmental monitoring and change assessment of toshka lakes in southern egypt using remote sensing. Environ. Earth Sci. 2015, 73, 3623–3632. [Google Scholar] [CrossRef]
  19. Zhu, X.L.; Chen, J.; Gao, F.; Chen, X.H.; Masek, J.G. An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions. Remote Sens. Environ. 2010, 114, 2610–2623. [Google Scholar] [CrossRef]
  20. Wu, M.Q.; Wu, C.Y.; Huang, W.J.; Niu, Z.; Wang, C.Y.; Li, W.; Hao, P.Y. An improved high spatial and temporal data fusion approach for combining Landsat and MODIS data to generate daily synthetic Landsat imagery. Inf. Fusion 2016, 31, 14–25. [Google Scholar] [CrossRef]
  21. Gao, F.; Masek, J.; Schwaller, M.; Hall, F. On the blending of the Landsat and MODIS surface reflectance: Predicting daily Landsat surface reflectance. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2207–2218. [Google Scholar]
  22. Wu, M.Q.; Niu, Z.; Wang, C.Y.; Wu, C.Y.; Wang, L. Use of MODIS and Landsat time series data to generate high-resolution temporal synthetic Landsat data using a spatial and temporal reflectance fusion model. J. Appl. Remote Sens. 2012, 6, 063507. [Google Scholar]
  23. Tian, H.; Wu, M.; Niu, Z.; Wang, C.; Zhao, X. Dryland crops recognition under complex planting structure based on radarsat-2 images. Trans. Chin. Soc. Agric. Eng. 2015, 31, 154–159. [Google Scholar]
  24. Long, S.; Fatoyinbo, T.E.; Policelli, F. Flood extent mapping for namibia using change detection and thresholding with SAR. Environ. Res. Lett. 2014, 9. [Google Scholar] [CrossRef]
  25. Torres, R.; Snoeij, P.; Geudtner, D.; Bibby, D.; Davidson, M.; Attema, E.; Potin, P.; Rommen, B.; Floury, N.; Brown, M.; et al. Gmes Sentinel-1 mission. Remote Sens. Environ. 2012, 120, 9–24. [Google Scholar] [CrossRef]
  26. Cazals, C.; Rapinel, S.; Frison, P.L.; Bonis, A.; Mercier, G.; Mallet, C.; Corgne, S.; Rudant, J.P. Mapping and characterization of hydrological dynamics in a coastal marsh using high temporal resolution Sentinel-1a images. Remote Sens. 2016, 8, 17. [Google Scholar] [CrossRef]
  27. Qiao, C.; Luo, J.C.; Sheng, Y.W.; Shen, Z.F.; Zhu, Z.W.; Ming, D.P. An adaptive water extraction method from remote sensing image based on NDWI. J. Indian Soc. Remote Sens. 2012, 40, 421–433. [Google Scholar] [CrossRef]
  28. Gu, Y.X.; Hunt, E.; Wardlow, B.; Basara, J.B.; Brown, J.F.; Verdin, J.P. Evaluation of MODIS NDVI and NDWI for vegetation drought monitoring using oklahoma mesonet soil moisture data. Geophys. Res. Lett. 2008, 35, 5. [Google Scholar] [CrossRef]
  29. Gao, B.C. NDWI-a normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  30. 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]
  31. Li, Y.; Zhang, Q.; Li, X.; Yao, J. Hydrolgical effects of Poyang Lake catchment in response to climate changes Resour. Environ. Yangtze Basin 2013, 22, 1339–1347. [Google Scholar]
  32. Qi, S.; Brown, D.G.; Tian, Q.; Jiang, L.; Zhao, T.; Bergen, K.A. Inundation extent and flood frequency mapping using Landsat imagery and digital elevation models. Gisci. Remote Sens. 2009, 46, 101–127. [Google Scholar] [CrossRef]
  33. Cai, X.; Gan, W.; Ji, W.; Zhao, X.; Wang, X.; Chen, X. Optimizing remote sensing-based level-area modeling of large lake wetlands: Case study of Poyang Lake. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 471–479. [Google Scholar] [CrossRef]
  34. Plank, S. Rapid damage assessment by means of multi-temporal SAR-a comprehensive review and outlook to Sentinel-1. Remote Sens. 2014, 6, 4870–4906. [Google Scholar] [CrossRef]
  35. Jung, H.S.; Lu, Z.; Zhang, L. Feasibility of along-track displacement measurement from Sentinel-1 interferometric wide-swath mode. IEEE Trans. Geosci. Remote Sens. 2013, 51, 573–578. [Google Scholar] [CrossRef]
  36. Li, K.; Shao, Y.; Zhang, F. Rice information extraction using multi-polarization airborne synthetic aperture radar data. J. Zhejiang Univ. Agric. Life Sci. 2011, 37, 181–186. [Google Scholar]
  37. Gupta, K.K.; Gupta, R. Despeckle and geographical feature extraction in SAR images by wavelet transform. Isprs J. Photogramm. Remote Sens. 2007, 62, 473–484. [Google Scholar] [CrossRef]
  38. Schmidt, K.S.; Skidmore, A.K. Smoothing vegetation spectra with wavelets. Int. J. Remote Sens. 2004, 25, 1167–1184. [Google Scholar] [CrossRef]
  39. Kyriou, A.; Nikolakopoulos, K. Flood mapping from Sentinel-1 and Landsat-8 data. A case study from river Evros, Greece. In Earth Resources and Environmental Remote Sensing/GIS Applications VI; Michel, U., Schulz, K., Ehlers, M., Nikolakopoulos, K.G., Civco, D., Eds.; SPIE—International Society for Optics and Photonics: Bellingham, WA, USA, 2015. [Google Scholar]
  40. Chai, T.; Draxler, R.R. Root mean square error (RMSE) or mean absolute error (MAE)?-arguments against avoiding RMSE in the literature. Geosci. Model Dev. 2014, 7, 1247–1250. [Google Scholar] [CrossRef]
  41. Mentaschi, L.; Besio, G.; Cassola, F.; Mazzino, A. Problems in RMSE-based wave model validations. Ocean Model. 2013, 72, 53–58. [Google Scholar] [CrossRef]
  42. De Myttenaere, A.; Golden, B.; Le Grand, B.; Rossi, F. Mean absolute percentage error for regression models. Neurocomputing 2016, 192, 38–48. [Google Scholar] [CrossRef]
  43. Zhong, Y.F.; Zhang, L.P.; Huang, B.; Li, P.X. An unsupervised artificial immune classifier for multi/hyperspectral remote sensing imagery. IEEE Trans. Geosci. Remote Sens. 2006, 44, 420–431. [Google Scholar] [CrossRef]
  44. Zhong, Y.F.; Zhang, L.P.; Gong, W. Unsupervised remote sensing image classification using an artificial immune network. Int. J. Remote Sens. 2011, 32, 5461–5483. [Google Scholar] [CrossRef]
  45. Yu, H.Y.; Gao, L.R.; Li, J.; Li, S.S.; Zhang, B.; Benediktsson, J.A. Spectral-spatial hyperspectral image classification using subspace-based support vector machines and adaptive markov random fields. Remote Sens. 2016, 8, 21. [Google Scholar] [CrossRef]
  46. Lin, Z.L.; Yan, L.M. A support vector machine classifier based on a new kernel function model for hyperspectral data. Gisci. Remote Sens. 2016, 53, 85–101. [Google Scholar] [CrossRef]
Figure 1. Location of Poyang Lake: (a) location of Jiangxi Province, China; (b) location of Poyang Lake in Jiangxi Province; and (c) Poyang Lake and its surrounding environment.
Figure 1. Location of Poyang Lake: (a) location of Jiangxi Province, China; (b) location of Poyang Lake in Jiangxi Province; and (c) Poyang Lake and its surrounding environment.
Remotesensing 09 00521 g001
Figure 2. Frequency histogram of SWI imagery on 27 September 2016. The dashed line represents water and the unbroken line represents non-water.
Figure 2. Frequency histogram of SWI imagery on 27 September 2016. The dashed line represents water and the unbroken line represents non-water.
Remotesensing 09 00521 g002
Figure 3. Spatiotemporal distribution of surface water in Poyang Lake based on the SWI from 24 May 2015 to 14 November 2016. The interval of the time-series was 12 days and there were 38 periods in total; however, five periods (4 August 2015; 8 November 2015; 23 June 2016; 29 July 2016; and 3 September 2016) were missing due to a lack of Sentinel-1A data. The yellow color represents non-water and the blue color represents water.
Figure 3. Spatiotemporal distribution of surface water in Poyang Lake based on the SWI from 24 May 2015 to 14 November 2016. The interval of the time-series was 12 days and there were 38 periods in total; however, five periods (4 August 2015; 8 November 2015; 23 June 2016; 29 July 2016; and 3 September 2016) were missing due to a lack of Sentinel-1A data. The yellow color represents non-water and the blue color represents water.
Remotesensing 09 00521 g003
Figure 4. The area of Poyang Lake from 24 May 2015 to 14 November 2016, with intervals of 12 days.
Figure 4. The area of Poyang Lake from 24 May 2015 to 14 November 2016, with intervals of 12 days.
Remotesensing 09 00521 g004
Figure 5. The examined SWI classification accuracy based on OLI: (a) OLI true color (R band 4, G band 3, B band 2) imagery acquired on 27 September 2016 located in Poyang Lake; (b) classification results based on Sentinel-1A imagery acquired on 27 September 2016; (c) SWI classification accuracy based on OLI on 27 September 2016; (d) OLI true color imagery acquired on 9 September 2015 located in Poyang Lake; (e) classification results based on Sentinel-1A imagery acquired on 9 September 2015; (f) SWI classification accuracy based on OLI on 9 September 2015; (g) OLI true color imagery acquired on 2 February 2016 located in Taihu Lake, China; (h) classification results based on Sentinel-1A imagery acquired on 2 February 2016; and (i) SWI classification accuracy based on OLI on 2 February 2016.
Figure 5. The examined SWI classification accuracy based on OLI: (a) OLI true color (R band 4, G band 3, B band 2) imagery acquired on 27 September 2016 located in Poyang Lake; (b) classification results based on Sentinel-1A imagery acquired on 27 September 2016; (c) SWI classification accuracy based on OLI on 27 September 2016; (d) OLI true color imagery acquired on 9 September 2015 located in Poyang Lake; (e) classification results based on Sentinel-1A imagery acquired on 9 September 2015; (f) SWI classification accuracy based on OLI on 9 September 2015; (g) OLI true color imagery acquired on 2 February 2016 located in Taihu Lake, China; (h) classification results based on Sentinel-1A imagery acquired on 2 February 2016; and (i) SWI classification accuracy based on OLI on 2 February 2016.
Remotesensing 09 00521 g005
Figure 6. Submergence ratios in Poyang Lake from 24 May 2015 to 14 November 2016. The red line shows the boundary of the sub-region. Poyang Lake was divided into Regions I–IV based on the submergence ratio to describe the changed characteristics in different regions.
Figure 6. Submergence ratios in Poyang Lake from 24 May 2015 to 14 November 2016. The red line shows the boundary of the sub-region. Poyang Lake was divided into Regions I–IV based on the submergence ratio to describe the changed characteristics in different regions.
Remotesensing 09 00521 g006
Figure 7. The number of lakes with different areas that constitute the Poyang Lake. The black line represents the lakes with an area of 0.1–1 km2, the red line represents the lakes with an area of >1–10 km2, and the blue line represents the lake with area of >10–100 km2. The background color is used to distinguish the rainy season and dry season: the light green color represents the rainy season, and the light yellow color represents the dry season.
Figure 7. The number of lakes with different areas that constitute the Poyang Lake. The black line represents the lakes with an area of 0.1–1 km2, the red line represents the lakes with an area of >1–10 km2, and the blue line represents the lake with area of >10–100 km2. The background color is used to distinguish the rainy season and dry season: the light green color represents the rainy season, and the light yellow color represents the dry season.
Remotesensing 09 00521 g007
Figure 8. Water surface area extracted from SWI and OLI from 24 May 2015 to 14 November 2016.
Figure 8. Water surface area extracted from SWI and OLI from 24 May 2015 to 14 November 2016.
Remotesensing 09 00521 g008
Figure 9. MODIS classification accuracy based on OLI.
Figure 9. MODIS classification accuracy based on OLI.
Remotesensing 09 00521 g009
Figure 10. The surface features in OLI of the commission and omission of water extracted from SWI images: (a) and (b) OLI true color (R band 4, G band 3, B band 2) imagery acquired on 9 September 2015.
Figure 10. The surface features in OLI of the commission and omission of water extracted from SWI images: (a) and (b) OLI true color (R band 4, G band 3, B band 2) imagery acquired on 9 September 2015.
Remotesensing 09 00521 g010
Figure 11. The spectral characteristic of misclassified zones and normal vegetation and water regarded as references in OLI on 9 September 2015.
Figure 11. The spectral characteristic of misclassified zones and normal vegetation and water regarded as references in OLI on 9 September 2015.
Remotesensing 09 00521 g011
Table 1. Examination of the Sentinel-1A water index equation using indexes.
Table 1. Examination of the Sentinel-1A water index equation using indexes.
RegionDateNumber of PixelsRMSEMAERMAPE (%)R2
Poyang Lake9 September 20152250.210.173.990.92
Taihu Lake2 February 20161900.320.215.320.92
Table 2. Confusion matrix for the SWI classification accuracy.
Table 2. Confusion matrix for the SWI classification accuracy.
Region and DateCategoryQA (%)PA (%)UA (%)OA (%)Kappa
Poyang Lake, 27 September 2016Water98.8491.5192.5896.520.8981
Non-water99.6797.9397.61
Poyang Lake, 9 September 2015Water99.1894.1893.4196.160.9102
Non-water99.6497.0597.40
Taihu Lake, 2 February 2016Water98.5197.1698.7398.620.9690
Non-water99.1999.3798.57
QA = quantity accuracy; PA = producer’s accuracy; UA = user’s accuracy; and OA = overall accuracy.
Table 3. The accuracy of water extraction using different methods and datasets acquired on 27 September 2016.
Table 3. The accuracy of water extraction using different methods and datasets acquired on 27 September 2016.
Category MethodDataAccuracy of Water (%)Running Time (s)
QAOAUAPAKappa
Threshold Classification Method (TCM)SWI98.8496.5292.5891.510.902
VH93.0294.5085.0791.010.842
VV94.3895.9993.3588.100.882
VH and VV87.0594.5283.2594.030.853
IsoData Unsupervised Classification (IUC)SWI98.0196.5593.0391.180.904
VH46.0187.3963.8998.380.694
VV79.3792.5777.4893.460.804
VH and VV72.2792.2775.4096.310.805
Minimum Distance Classification (MDC)SWI96.3896.5993.8590.450.903
VH82.0793.5880.0494.400.823
VV94.7895.1587.0691.610.863
VH and VV96.6696.0689.7492.730.893
Support Vector Machine (SVM)SWI99.5796.4291.6992.080.9013
VH97.1194.5886.6589.160.84193
VV94.6395.1386.9691.630.8652
VH and VV97.8496.1990.4992.440.8918
QA = quantity accuracy; OA = overall accuracy; UA = user’s accuracy; PA = producer’s accuracy. The test zone location in the Poyang Lake region with an area of 110 km × 70 km. The ground-truth water area was extracted from OLI acquired on 27 September 2016. The last column shows the time needed to extract the water extent using different methods based on the same hardware. The unit is seconds.

Share and Cite

MDPI and ACS Style

Tian, H.; Li, W.; Wu, M.; Huang, N.; Li, G.; Li, X.; Niu, Z. Dynamic Monitoring of the Largest Freshwater Lake in China Using a New Water Index Derived from High Spatiotemporal Resolution Sentinel-1A Data. Remote Sens. 2017, 9, 521. https://doi.org/10.3390/rs9060521

AMA Style

Tian H, Li W, Wu M, Huang N, Li G, Li X, Niu Z. Dynamic Monitoring of the Largest Freshwater Lake in China Using a New Water Index Derived from High Spatiotemporal Resolution Sentinel-1A Data. Remote Sensing. 2017; 9(6):521. https://doi.org/10.3390/rs9060521

Chicago/Turabian Style

Tian, Haifeng, Wang Li, Mingquan Wu, Ni Huang, Guodong Li, Xiang Li, and Zheng Niu. 2017. "Dynamic Monitoring of the Largest Freshwater Lake in China Using a New Water Index Derived from High Spatiotemporal Resolution Sentinel-1A Data" Remote Sensing 9, no. 6: 521. https://doi.org/10.3390/rs9060521

APA Style

Tian, H., Li, W., Wu, M., Huang, N., Li, G., Li, X., & Niu, Z. (2017). Dynamic Monitoring of the Largest Freshwater Lake in China Using a New Water Index Derived from High Spatiotemporal Resolution Sentinel-1A Data. Remote Sensing, 9(6), 521. https://doi.org/10.3390/rs9060521

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