Land Surface Water Mapping Using Multi-Scale Level Sets and a Visual Saliency Model from SAR Images

Land surface water mapping is one of the most basic classification tasks to distinguish water bodies from dry land surfaces. In this paper, a water mapping method was proposed based on multi-scale level sets and a visual saliency model (MLSVS), to overcome the lack of an operational solution for automatically, rapidly and reliably extracting water from large-area and fine spatial resolution Synthetic Aperture Radar (SAR) images. This paper has two main contributions, as follows: (1) The method integrated the advantages of both level sets and the visual saliency model. First, the visual saliency map was applied to detect the suspected water regions (SWR), and then the level set method only needed to be applied to the SWR regions to accurately extract the water bodies, thereby yielding a simultaneous reduction in time cost and increase in accuracy; (2) In order to make the classical Itti model more suitable for extracting water in SAR imagery, an improved texture weighted with the Itti model (TW-Itti) is employed to detect those suspected water regions, which take into account texture features generated by the Gray Level Co-occurrence Matrix (GLCM) algorithm, Furthermore, a novel calculation method for center-surround differences was merged into this model. The proposed method was tested on both Radarsat-2 and TerraSAR-X images, and experiments demonstrated the effectiveness of the proposed method, the overall accuracy of water mapping is 98.48% and the Kappa coefficient is 0.856.


Introduction
Water resources are an irreplaceable strategic resource for human survival.Land surface water (LSW) is an important part of the water cycle.LSW mapping, using remote sensing techniques, plays an important role in wetland monitoring, flood monitoring, flood disaster assessment, surface water area estimation, and water resources management.Although multi-spectral optical images have been found to be effective when the sky is clear, there are some limitations in overcast weather conditions.However, due to their capabilities of large-area coverage, cloud penetration and all-weather acquisition, the SAR image data could compensate for the weakness in optical images in such conditions for LSW mapping [1].This paper proposes a fast and reliable land surface water mapping method of large-area and high-resolution SAR images.Only calm water was considered including calm rivers and lakes, not including puddles and paddy fields.
Water classification is a particular case of SAR image classification, with only two classes to assign; it is widely used in applications such as wetland monitoring, flood monitoring, surface water area estimation and so on.However, due to the presence of coherent speckle, which can be modeled as proposed a method for detecting multiple objects in cluttered scenes based on a bottom-up visual attention model.Li et al. [22] proposed a wavelet transform saliency model for generating a saliency map to detect natural objects.Gao et al. [23] presented a novel detection model based on a visual attention mechanism to detect ships in remote sensing images.Li and Itti [24] integrated saliency features acquired by a visual salience model, and with gist features, to detect objects, such as ships, buildings and airplanes in remote sensing images.However, these models are aimed at analyzing images of nature or at extracting ships, airplanes, etc. from optical satellite images, and they are not designed for SAR images which make it difficult to apply them directly.
The main purpose of this paper is to propose a fast and accurate land surface water mapping method, and most important is that it can solve some complex situation, such as shadow, farmland confusion, etc.The remainder of this paper is structured as follows.In Section 2, our method is formulated and applied to SAR image water extraction.Firstly, a TW-Itti model is proposed.Secondly, we describe the adaptive multi-scale level set method based on Gamma model.Thirdly, a post-processing method is described.Some experimental results are shown in Section 3, before we conclude the paper in Section 4.

Proposed Method
Given the above problems, a novel water segmentation method was proposed based on the integration of multi-scale level sets and a visual salience model, as follows, (1) To improve the accuracy and efficiency of water extraction, the classic Itti model is first introduced to detect the SWR.After that, the level set method can be applied only to these regions, and at the same time, the visual salience map can be used for initializing the zero level set function.Thus, we can detect water regions rapidly using Itti model and we can accurately segment water from those regions using level sets; (2) The classical Itti model is not designed for SAR water images, thus, an improved visual saliency attention model was presented.First, we replace the color images in the Itti model with textured images, which are generated using the gray-level co-occurrence matrix (GLCM) algorithm.Based on the grey values of the water body, a weight coefficient is introduced to calculate the center-surround differences; (3) In order to better avoid "salt and pepper" effects, multi-scale technology is introduced in the level sets method.A multi-scale level set method was proposed by Sui et al. [25], the difference with our approach being that we use OTSU algorithms [26] to adaptively define the initial level set function; (4) To detect objects that may be confused with calm water, a post-processing method was proposed based on object-oriented geometrical features.
In this section, the framework of our rapid water extraction method for SAR imagery was presented, as shown in Figure 1.
The principal steps of our method are as follows.
Step 1: Generate a visual salience map using the improved Itti model.
Step 2: Decompose the original image and visual saliency map into L scales by the bilinear interpolation algorithm.Let K = L.
Step 3: SWR regions are obtained: segment the saliency map into binary images at each scale using the OTSU method, and then objects in these binary images can be marked using the connecting area algorithm.Next, SWR regions (black rectangles in the SAR image shown in Figure 1) can be obtained by calculating the maximum external rectangle of each marked object.
Step 4: Use the OTSU algorithms to initialize the level set function at scale L. Obtain the scale K segmentation results using the level set method with the Gamma model.
Step 5: K = K ´1.If K ě 0, return to Step 4, until the final segmentation result is obtained.
Step 6: Post-processing: the weighting of the smoothness within the shape is applied to the final binary image to distinguish paddy fields from rivers or lakes (see Section 2.3).

differences;
(3)In order to better avoid "salt and pepper" effects, multi-scale technology is introduced in the level sets method.A multi-scale level set method was proposed by Sui et al. [25], the difference with our approach being that we use OTSU algorithms [26] to adaptively define the initial level set function; (4) To detect objects that may be confused with calm water, a post-processing method was proposed based on object-oriented geometrical features.
In this section, the framework of our rapid water extraction method for SAR imagery was presented, as shown in Figure 1.

Improved TW-Itti Model
The model proposed in this paper draws its inspiration from the bottom-up attention model proposed by Itti et al. [20] (see the left part of Figure 1).Our contribution is focused mainly on the extraction of early visual features.Instead of using color information, the textures generated by the GLCM are applied to obtain the feature maps.Since we used single-polarized and single-band SAR images, no colors can be acquired; thus, we fused the texture feature to generate color images.The benefit of doing this is to enhance the visibility of the water body and the background.Moreover, in view of the dark grey value of water, a new calculation method for center-surround differences is presented.Center-surround is implemented in the model as the difference between fine and coarse scales, The center is a pixel at scale c P t2, 3, 4u, and the surround is the corresponding pixel at scale s " c `δ, with δ P t3, 4u.The across-scale difference between two maps, denoted "a" below, is obtained by interpolation to the finer scale and point-by-point subtraction.The input is provided as a static gray image, and nine spatial scales are generated using dyadic Gaussian pyramids, which progressively low-pass filter and subsample the input image.I is used to create a Gaussian pyramid I(σ), where σP(0 . . .8) is the scale.The center-surround differences (a defined previously) between a"center" fine scale c and a "surround" coarser scale s yield the feature maps.
The first set of feature maps is concerned with intensity contrast, which, in mammals, is detected by neurons sensitive either to dark centers on bright surrounds or to bright centers on dark surrounds.Here, both types of sensitivities are simultaneously computed (using a rectification) in a set of six maps I pc, sq.
I pc, sq " |I pcq a I psq|, c P t2, 3, 4u , s " c `σ, σ P t3, 4u Considering the dark grey value of water, a weight coefficient is introduced to calculate the center-surround differences, mainly to maintain the low reflective value of water.Thus, Equation (1) can be rewritten as below: where ω is the weight coefficient and ω P r0, 1s.
A second set of maps is similarly constructed for the texture channels.A texture image T is generated using GLCM and principal components analysis (PCA) algorithm.Here, T is comprised of five conventional textures, i.e., energy, correlation, inertia, entropy and inverse difference moment).The window size is 7 ˆ7, the distance value is 1 pixel and there are 16 grey levels.After acquiring these textures, the PCA method is used to obtain the first principal component that represents texture image T. The texture feature maps T(c,s) can be represented as follows: T pc, sq " |T pcq a T psq| (3) The GLCM-based features, proposed by Haralick et al. [27], have been used widely in texture segmentation.They are defined for the neighborhood of an input image and represented by using the second order joint condition probability density p ij , which is the relative frequency at which two image pixels are related.The construction of the GLCM depends on three parameters, i.e., the relative distance (d), the window size (w), and the orientation (θ) between pixels in a pair.In this paper, the distanced is set to one pixel, the window size is set to 7 ˆ7 pixels, while the orientation is equal to the average value of four directions (0 ˝, 45 ˝, 90 ˝and 135 ˝).Next, we extract five features from the GLCM, as follows: Energy: ASM " Correlation: Inertia: Entropy: ENT " ´ÿ i ÿ j ppi, jqlogppi, jq Inverse Difference Moment IDM " where u x , u y , σ x , σ y are the means and standard deviations of the marginal distribution associated with P pi, jq {R, and R is a normalizing constant.
After acquiring these five texture features, they are fused using PCA where the principal components are then used for segmentation purposes.Finally, the fuzzy-C means algorithm is performed with the first three principal components during image segmentation.
The third set of maps is the local orientation O, here O is obtained from I using oriented Gabor pyramids O pσ, θq, where σ P p0 . . .8q represents the scale and θ P t0 ˝, 45 ˝, 90 ˝, 135 ˝u is the preferred orientation.The orientation feature maps, O pc, s, θq, are encoded as a group, with a local orientation contrast between the center and surround scales, as follows: In total, 36 feature maps were computed: six for intensity, six for texture, and 24 for orientation.After processing of across-scale combinations, normalization and linear combinations with six for intensity, six for texture, and 24 for orientation.As mentioned in [20], the improved salience map is generated.
After obtaining the saliency map, suspected water regions can be found by using the Otsu method.The Otsu algorithm is a nonparametric and unsupervised method of automatic threshold selection for image segmentation, and it is a classical method for optical image segmentation.However, due to the strong speckle noise of SAR images, especially for fine spatial resolution SAR images, Otsu method could be invalidated.Thus, it was applied to the visual saliency map, which has no speckle noise, and it is used to segment the visual saliency map and obtain the suspected water regions.At the same time, we can also use this segmentation result of the visual saliency map to initialize the zero level set function, supposing the segmentation result of the visual saliency map is t pxq.The details for using the visual saliency map to initialize the zero level set function can be seen in the last part of Section 2.2.

Adaptive Multi-Scale Level Set Method Based on the Gamma Model
In order to remove the influence of speckle and preserve important structural information, a multi-scale level set algorithm is introduced.We acquire multi-scale images at several scales by decomposing the SAR image using the block averaging algorithm.Because the Gamma model can well represent the distribution of SAR images if the speckle is fully developed and spatially uncorrelated, and the radar backscatter of the region is constant, it is employed to define the energy functional of the level set function for the original SAR image.Even if there is correlation between neighboring pixels, or the radar backscatter features texture, the Gamma distribution can still be used as an approximation by adjusting the number of looks L. For scale images decomposed by the block averaging algorithm, the probability density function (PDF) of the pixel intensity can be also given by a Gamma distribution.For a given image u pxq, r u pxq is the decomposed image generated by block averaging algorithm.If u pxq is modeled by a Gamma distribution as u pxq "Γ pµ, Lq, then r u pxq "Γ pµ, 4Lq.After using the threshold level set initialization method at the coarsest scale, the initial segmentation result, which is marked by a label, is generated.This label is used not only for the level set initialization at finer scales but also as an additional constraint to guide contour evolution.Consequently, the finest scale image is extracted as the final result.
Chan and Vese [17] proposed a model that implemented the Mumford-Shah functional via the level set function for the purpose of bimodal segmentation.The segmentation is performed using an active contour model without boundaries.Let Ω be a bounded open subset of R 2 , with BΩ being its boundary.Let u 0 px, yq : Ω Ñ R be a given image, and C be a curve in the image domain Ω.Segmentation is achieved via the evolution of curve C, which is the basic idea of the active contour model.In the level set method, C Ă Ω is represented by the zero level set of a Lipschitz function φ : Ω Ñ R, and we replace the unknown variable C with the unknown variable φ.
In this paper, assume that the image is partitioned into two classes Ω 1 and Ω 2 , which are separated by a curve C, and that classes Ω 1 and Ω 2 are modeled using the probability density functions (pdfs) P 1 and P 2 , respectively.
The partition is obtained by minimizing the following energy function: A gamma model was used for high resolution SAR image segmentation.Suppose u SAR px, yq is a SAR image, we model the image in each region R i by a Gamma distribution of mean intensity u i and number of looks L: Using the Heaviside function H, and the one-dimensional Dirac function δ 0 , which are defined, respectively, by The segmentation is performed via the evolution of φ by minimizing the following energy function: where µ, λ 1 , λ 2 are non-negative weighted parameters.Function ∅ represents class Ω 1 for ∅ ą 0, and The evolution of φ is governed by the following partial differential motion equation: where δ ε pφq is a regularized version of the Dirac function.
Next, we aim to estimate the Gamma parameters θ = {µ i }.Using maximum likelihood estimation θ ˚" argmax θ logp py|θq.Assuming that the samples y j , j " 1, . . ., N in each region are independent and identically distributed, the log likelihood is logp py|θq " log ś N j"1 p `yj ˇˇθ ˘" ř N j"1 logp `yj ˇˇθ ˘.Taking the derivative of logp py|θq with respect to θ and setting them equal to zero, µ i " ř N i j"1 y j {N i was obtained, where N i is the pixel number in Ω i .
As we know that the segmentation result of the visual saliency map is t pxq, then we simply initialize the level set function φ as follows: The main step of the level sets method is as follows (also see the details for Step 4 of our method): (1) Initialize the level set function using Equation ( 15); (2) Evolve the level set function φ according to Equations ( 13) and ( 14); (3) Check whether the evolution is stationary.If not, return to Step 2.

Post-Processing: Using Object-Oriented Geometrical Feature
Water objects in SAR images are continuous regions with low light, and some objects with a low radar cross-section similar to calm water, such as urban roads, airstrips, building shadows or paddy fields, may influence the results of water extraction.
Compared with a large water area, urban roads, airstrips and building shadows are not significant, thus, these areas are ignored after we acquire the SWR regions using the visual salience model.However, some paddy fields remain as salient objects because of their occupation of large areas.In this study, geometrical features are used for post-processing to distinguish paddy fields from water bodies.First, a fast pixel labeling method is applied to search for connected areas in binary images (more details can be found in [28]).Then, the attributes of each object can be acquired, such as the area, perimeter and number of objects.In SAR images, the geometrical features of paddy fields usually appear as regular rectangles, whereas, water bodies such as rivers and lakes, mostly appear as irregular geometrical shape.Thus, inspired by the user guide of eCognition software (Definiens AG, Munich, Germany) [29], the weighting of smoothness (WS) within the shape was introduced to distinguish paddy fields from water bodies, as shown in Equation (16).

WS "
Area of Minimum Enclosing Rectangle Area of the object (16) ISPRS Int.J. Geo-Inf.2016, 5, 58 We obtained good segmentation results using the proposed method, and most paddy fields had geometrical shapes such are as shown in Figure 2a.However, after the evolution of the level sets, some paddy fields were connected each other, such as the shape shown in Figure 2b. Figure 2c shows the object shape of a river and Figure 2d shows the object shape of a lake.The red dashed rectangles are the minimum enclosing rectangles for each object.During image segmentation, the ratio of the minimum enclosing rectangle's area and the object's area is used to make the segmentation object's shape more regular, and the value of the ratio is very close to 1, so the object's shape is very close to a rectangle.The ratio values for Figure 2a-d are 1.344, 2.143, 6.524 and 2.5869, respectively.Thus, based on the experience value (3000 objects with 15 images were tested), we assume that an object is "water class", if the value of smoothness weighting in an object is greater than 2, and all others water classes can be removed during post-processing.

Experimental Data Set
SAR amplitude images (Radarsat2 and TerraSAR-X) of river basin regions in China were used in the experiment, detailed descriptions of the experiment data can be seen in Table 1.

Results and Discussion
The following three methods were compared to evaluate the performance of the proposed method.
ALG1: A MRF segmentation method proposed by Picco and Placio [30].ALG2: A method for the separation of land and water in SAR images which was proposed by Silveira and Heleno [14].
ALG3: A multi-scale level set method for SAR image segmentation which was proposed by Sui et al. [25].
Note that all of our images were tested without a filter when using the proposed method and the scale level was 2. The value of the smooth index WS = 2, i.e., if the value of the object smoothing index was less than 2, this object could be removed.The weight coefficient of the center-surround differences ω is equal to 0.1.

Experimental Data Set
SAR amplitude images (Radarsat2 and TerraSAR-X) of river basin regions in China were used in the experiment, detailed descriptions of the experiment data can be seen in Table 1.

Results and Discussion
The following three methods were compared to evaluate the performance of the proposed method.
ALG1: A MRF segmentation method proposed by Picco and Placio [30].ALG2: A method for the separation of land and water in SAR images which was proposed by Silveira and Heleno [14].
ALG3: A multi-scale level set method for SAR image segmentation which was proposed by Sui et al. [25].
Note that all of our images were tested without a filter when using the proposed method and the scale level was 2. The value of the smooth index WS = 2, i.e., if the value of the object smoothing index was less than 2, this object could be removed.The weight coefficient of the center-surround differences ω is equal to 0.1.

ISPRS Int. J. Geo-Inf. 2016, 5, 58 9 of 19
China, which was acquired by Radarsat-2 in 2009 in descending VV polarization mode, with a resolution of 3 m and a size of 3024 × 2263 pixels.Figure 3b,c shows the saliency map generated using Itti model and our TW-Itti model, respectively.The experiment validates our model performance much better than the Itti model.In Figure 3b, except for water body, other features are also incorrectly considered to be salient, which made the water body discontinuous.Figure 3d shows the Otsu segmentation result after the post-processing of Figure 3c.It can be seen from Figure 3d that segmentation using the saliency map can acquire main outlines, however, missing edge details.It can be seen that the edges are jagged and irregular.As such, we need to apply the level set method to refine the edge and make it smooth.Figure 4 shows the water extraction results for the Huai River. Figure 4a is the SWR generated by the saliency map (see Figure 3d) after post-processing.In general, the actual positions of boundaries within a SAR scene are unknown.Thus, the segmentation quality measures were modified to allow a comparison of the automatic segmentation approaches and manual segmentation, ISPRS Int.J. Geo-Inf.2016, 5, 58 9 of 19 China, which was acquired by Radarsat-2 in 2009 in descending VV polarization mode, with a resolution of 3 m and a size of 3024 × 2263 pixels.Figure 3b,c shows the saliency map generated using Itti model and our TW-Itti model, respectively.The experiment validates our model performance much better than the Itti model.In Figure 3b, except for water body, other features are also incorrectly considered to be salient, which made the water body discontinuous.Figure 3d shows the Otsu segmentation result after the post-processing of Figure 3c.It can be seen from Figure 3d that segmentation using the saliency map can acquire main outlines, however, missing edge details.It can be seen that the edges are jagged and irregular.As such, we need to apply the level set method to refine the edge and make it smooth.Figure 4 shows the water extraction results for the Huai River. Figure 4a is the SWR generated by the saliency map (see Figure 3d) after post-processing.In general, the actual positions of boundaries within a SAR scene are unknown.Thus, the segmentation quality measures were modified to allow a comparison of the automatic segmentation approaches and manual segmentation,

Experiment on Radarsat-2 Image
Figure 3 compares the salience maps generated using the Itti model and our model.This shows clearly that our model readily distinguished areas of water.Figure 3a shows the Huai River in eastern China, which was acquired by Radarsat-2 in 2009 in descending VV polarization mode, with a resolution of 3 m and a size of 3024 ˆ2263 pixels.Figure 3b,c shows the saliency map generated using Itti model and our TW-Itti model, respectively.The experiment validates our model performance much better than the Itti model.In Figure 3b, except for water body, other features are also incorrectly considered to be salient, which made the water body discontinuous.Figure 3d shows the Otsu segmentation result after the post-processing of Figure 3c.It can be seen from Figure 3d that segmentation using the saliency map can acquire main outlines, however, missing edge details.It can be seen that the edges are jagged and irregular.As such, we need to apply the level set method to refine the edge and make it smooth.as shown in Figure 4b.Figures 4c-f show segmentation results by ALG1, ALG2, ALG3 and our method, respectively.Figure 4g shows the post-processing results of the proposed method.The experiment validates that our method (see Figure 4g) performs best in distinguishing water and preventing speckle.Detailed comparisons of accuracy and efficiency are given in Table 1.

Experiment on TerraSAR-X Imagery
Figure 5 shows a comparative saliency map using a TerraSAR-X image.Figure 5a    Figure 4 shows the water extraction results for the Huai River. Figure 4a is the SWR generated by the saliency map (see Figure 3d) after post-processing.In general, the actual positions of boundaries within a SAR scene are unknown.Thus, the segmentation quality measures were modified to allow a comparison of the automatic segmentation approaches and manual segmentation, as shown in Figure 4b. Figure 4c-f shows segmentation results by ALG1, ALG2, ALG3 and our method, respectively.Figure 4g shows the post-processing results of the proposed method.The experiment validates that our method (see Figure 4g) performs best in distinguishing water and preventing speckle.Detailed comparisons of accuracy and efficiency are given in Table 1.

Experiment on TerraSAR-X Imagery
Figure 5 shows a comparative saliency map using a TerraSAR-X image.Figure 5a shows the Hanjiang River in central China, which was acquired by TSX-1 in 2008 in a descending VV polarization mode, with a resolution of 1 m and a size of 7153 ˆ1948 pixels.After comparing Figure 5b,c, the same conclusions about the model improvements in Radarsat-2 were applicable.
Figure 6 shows the water extraction results for the Hanjiang River. Figure 6a shows the SWR generated using the saliency map (see Figure 5c).Figure 6b is the manual segmentation result.Figure 6c-f show segmentation result by ALG1, ALG2, ALG3 and our method, respectively.Figure 6g shows the post-processing result of the proposed method.The experiment also validates that our method can achieve the best segment result.The speckle noise is very strong in this image, and we can see that there are many very small-scale objects (i.e., "salt and pepper" effects) in Figure 6c,d, and there are many other features misrepresented as water.Figure 6e is much better than the previous methods, however, compared with Figure 6f, there are still some noise effects.Detailed comparisons of accuracy and efficiency are given in Table 1.

Experiment on TerraSAR-X Imagery
Figure 5 shows a comparative saliency map using a TerraSAR-X image.Figure 5a   Figure 6 shows the water extraction results for the Hanjiang River. Figure 6a shows the SWR generated using the saliency map (see Figure 5c).Figure 6b is the manual segmentation result.Figure 6c-f show segmentation result by ALG1, ALG2, ALG3 and our method, respectively.Figure 6g shows the post-processing result of the proposed method.The experiment also validates that our method can achieve the best segment result.The speckle noise is very strong in this image, and we can see that there are many very small-scale objects (i.e., "salt and pepper" effects) in Figure 6c,d, and there are many other features misrepresented as water.Figure 6e is much better than the previous methods, however, compared with Figure 6f, there are still some noise effects.Detailed comparisons of accuracy and efficiency are given in Table 1. Figure 6 shows the water extraction results for the Hanjiang River. Figure 6a shows the SWR generated using the saliency map (see Figure 5c).Figure 6b is the manual segmentation result.Figure 6c-f show segmentation result by ALG1, ALG2, ALG3 and our method, respectively.Figure 6g shows the post-processing result of the proposed method.The experiment also validates that our method can achieve the best segment result.The speckle noise is very strong in this image, and we can see that there are many very small-scale objects (i.e., "salt and pepper" effects) in Figure 6c,d, and there are many other features misrepresented as water.Figure 6e is much better than the previous methods, however, compared with Figure 6f, there are still some noise effects.Detailed comparisons of accuracy and efficiency are given in Table 1.

Applications to Large-Area SAR Images
To demonstrate the speed and reliability of our method, we tested the proposed method using real large-area SAR images.The Changjiang River basin, which is one of the main river basins in China, was selected for this experiment.Figure 7 shows the water extraction results for the Changjiang River. Figure 7a shows the original SAR image with an image size of 15540 ×14550 pixels.Figure 7b shows the saliency map generated by Itti model.Figure 7c is the salience map generated using the TW-Itti model and Figure 7d shows the water segmentation result using our method.The time required to generate this image was 1209.463s.

Applications to Large-Area SAR Images
To demonstrate the speed and reliability of our method, we tested the proposed method using real large-area SAR images.The Changjiang River basin, which is one of the main river basins in China, was selected for this experiment.Figure 7 shows the water extraction results for the Changjiang River. Figure 7a shows the original SAR image with an image size of 15540 ˆ14550 pixels.Figure 7b shows the saliency map generated by Itti model.Figure 7c is the salience map generated using the TW-Itti model and Figure 7d shows the water segmentation result using our method.The time required to generate this image was 1209.463s.

Applications to Large-Area SAR Images
To demonstrate the speed and reliability of our method, we tested the proposed method using real large-area SAR images.The Changjiang River basin, which is one of the main river basins in China, was selected for this experiment.Figure 7 shows the water extraction results for the Changjiang River. Figure 7a shows the original SAR image with an image size of 15540 ×14550 pixels.Figure 7b shows the saliency map generated by Itti model.Figure 7c is the salience map generated using the TW-Itti model and Figure 7d shows the water segmentation result using our method.The time required to generate this image was 1209.463s.

Accuracy Analysis and Discussion
In order to perform quantitative analysis, overall accuracy (short for OA) and Kappa coefficient are adopted as the accuracy assessment tools in this paper.Assuming that the size of an original SAR image is  × , we denote the label image as ′, whose size is also  × .Correspondingly, R represents the label image from manual extraction (it also indicates an ideal water extraction) for the same original SAR image.The error image is therefore defined as  =  ′ − .
where l denotes the number of non-zero pixels in E. Kappa coefficient can be described as below: where is the row numbers of the confusion matrix,   is the value of principal diagonal,  + and  + are the sum value of row i and column i, respectively.Table 2 compares the overall accuracy and Kappa coefficient of the above methods.

Accuracy Analysis and Discussion
In order to perform quantitative analysis, overall accuracy (short for OA) and Kappa coefficient are adopted as the accuracy assessment tools in this paper.Assuming that the size of an original SAR image is M ˆN, we denote the label image as X 1 , whose size is also M ˆN.Correspondingly, R represents the label image from manual extraction (it also indicates an ideal water extraction) for the same original SAR image.The error image is therefore defined as E " X 1 ´R.
where l denotes the number of non-zero pixels in E. Kappa coefficient can be described as below: Kappa " pM ˆNq where r is the row numbers of the confusion matrix, x ii is the value of principal diagonal, x i`a nd x `i are the sum value of row i and column i, respectively.Table 2 compares the overall accuracy and Kappa coefficient of the above methods.
Based on the qualitative analysis (Experiments 1 and 2) and the quantitative analysis (Table 1), the following conclusions can be reached: (1) Compared with the classical Itti model, the improved visual saliency model provides better enhancement of water information and it suppresses other surface features (See Figure 2c,d for Experiment 1and Figure 4b,c for Experiment 2).( 2) Compared with the state-of-the-art methods (ALG1, ALG2 and ALG3), the method of integration using the multi-scale technique and visual saliency model produced better accuracy.(3) In all cases, our method frequently produced superior results and it produced the most accurate results in the shortest time (see Figure 4g for Experiment 1, Figure 6g for Experiment 2 and the bold line in Table 1).Moreover, in order to better evaluate the performance of the proposed method, the error matrix of the results for ALG1, ALG2, ALG3 and our method were calculated for validation; in this paper, we only take the dataset of Radarsat-2 (Experiment 1) as an example.An error matrix is a very effective way to assess accuracy in that the accuracies of each category are plainly described along with the errors of omission and commission present in the classification.The error matrixes of Experiment 1 are shown in Tables 3-6.Overall accuracies for ALG1, ALG2, ALG3 are 79.22% and 92.39%, 91.38%, and the Kappa coefficients are 0.262, 0.53 and 0.511.Individual class accuracies were highly variable, ranging from 20.27% to 99.51%.For our method, overall accuracy is 98.48% and the Kappa coefficient is 0.856.Moreover, individual class accuracies are all very high, above 80%.This proves that our classifier is accurate and robust.

Conclusions
The main purpose of this study was to devise a water mapping method that improves accuracy by using a visual attention-based method for multi-scale level sets, particularly in areas with roads, airstrips, radar shadows and paddy fields, which are often major causes of low classification accuracy.Experiments demonstrated that our method significantly improved accuracy in areas where shadow and other dark surfaces were the main sources of classification errors.Our method was aimed at calm water bodies and, therefore, did not consider dynamic water bodies, such as flood water with high wind.In future work, we will focus more on other types of water bodies, looking at, for instance, flood-induced backscatter changes in SAR data [31] and so on.
Although experiments have testified that our method gets satisfactory results, much work remains to be done.The multi-scale analysis framework is a new component in the level set method.In future work, multi-region level sets method can be merged with multi-scale analysis, and, as such, more land features can be considered.Moreover, we will develop our method to be used in full polarized SAR images for acquiring more accurate classification.

Figure 1 .
Figure 1.The basic principle of the proposed method.(Step1, saliency map generated from original image; Step 2, Decomposed the saliency map into L scales; Step3, Obtained SWR regions in different scales; Step4, Segmented the water in SWR regions; Step5, refine the segmented result; Step6, post-processing to acquire accurate water body.)

Figure 2 .
Figure 2. The demonstration of the weighting of smoothness within the objects.(a) Single paddy field; (b) connected paddy fields; (c) river; (d) lake.

Figure 2 .
Figure 2. The demonstration of the weighting of smoothness within the objects.(a) Single paddy field; (b) connected paddy fields; (c) river; (d) lake.
Figure5shows a comparative saliency map using a TerraSAR-X image.Figure5ashows the Hanjiang River in central China, which was acquired by TSX-1 in 2008 in a descending VV polarization mode, with a resolution of 1 m and a size of 7153 × 1948 pixels.After comparing Figure5b,c, the same conclusions about the model improvements in Radarsat-2 were applicable.
Figure5shows a comparative saliency map using a TerraSAR-X image.Figure5ashows the Hanjiang River in central China, which was acquired by TSX-1 in 2008 in a descending VV polarization mode, with a resolution of 1 m and a size of 7153 × 1948 pixels.After comparing Figure5b,c, the same conclusions about the model improvements in Radarsat-2 were applicable.

Figure 7 .
Figure 7. Salience map comparison and segmentation result (Changjiang River).(a) original SAR image; (b) the saliency map generated by Itti model; (c) the saliency map generated by TW-Itti model; (d) LSW map obtained using our method.

Table 1 .
Experimental data of different water bodies.

Study Area of Huai River: Huai
(2)er catchment around eastern China is taken as one of the study areas.The Radarsat 2 images (VV polarization) at a spatial resolution of 3 m were acquired on 7 December 2009.A Google map image of 0.5m resolution with the same region and same season is used as the true water class image, then manual water extraction in this google image is used for reference to check the accuracy of the water extraction process.(2)

Study Area of Hanjiang and Changjiang River:
Hanjiang River is the largest branch of Changjiang River, and they are located in Wuhan City, Hubei Province, China.The TerraSAR-X images (VV polarization) at a spatial resolution of 1 m were acquired on 9 October 2008.A vector map of the same region supported by Map Institute of Hubei Province is used for reference to check the accuracy of the water extraction process.

Table 1 .
Experimental data of different water bodies.
map of the same region supported by Map Institute of Hubei Province is used for reference to check the accuracy of the water extraction process.

Table 2 .
Overall Accuracy and Kappa Coefficient Comparison.

Table 2 .
Overall Accuracy and Kappa Coefficient Comparison.

Table 6 .
Error matrix of the proposed method classified image (Dataset 1).