Open AccessArticle

A Wavelet Neural Network for SAR Image Segmentation

Xian-Bin Wen

^1,2,*,

Hua Zhang

^1,2 and

Fa-Yu Wang

^1,2

Key Laboratory of Computer Vision and System of Ministry of Education, Tianjin University of Technology, Tianjin 300191, China

Tianjin Key Laboratory of Intelligence Computing and Novel Software Technology, Tianjin, 300191, China

Author to whom correspondence should be addressed.

Sensors 2009, 9(9), 7509-7515; https://doi.org/10.3390/s90907509

Submission received: 29 July 2009 / Revised: 9 September 2009 / Accepted: 12 September 2009 / Published: 22 September 2009

(This article belongs to the Special Issue Neural Networks and Sensors)

Download

Browse Figures

Versions Notes

Abstract

This paper proposes a wavelet neural network (WNN) for SAR image segmentation by combining the wavelet transform and an artificial neural network. The WNN combines the multiscale analysis ability of the wavelet transform and the classification capability of the artificial neural network by setting the wavelet function as the transfer function of the neural network. Several SAR images are segmented by the network whose transfer functions are the Morlet and Mexihat functions, respectively. The experimental results show the proposed method is very effective and accurate.

Keywords:

synthetic aperture radar; image segmentation; Wavelet Neural Network

1. Introduction

The synthetic aperture radar (SAR) system is a powerful tool for observing the Earth under all weather conditions. In recent years, SAR imaging has been rapidly gaining prominence in applications such as remote sensing, surface surveillance and automatic target recognition. Segmentation of SAR images is a critical preliminary operation in various SAR images processing applications, such as target detection, recognition, and image compression.

SAR images characteristically have a particular kind of noise, called speckle, which occurs by random interferences, either constructive or destructive, between electromagnetic waves from different reflections in the imaged area. This makes SAR segmentation a difficult task, though several different segmentation methods designed specifically for SAR images have been proposed. Three common methods are optical image segmentation after speckle filter, the multiscale method [1–3], and the neural networks method [4,5].

Artificial neural networks (ANNs) are a class of computational architectures that are composed of interconnected, simple processing nodes with weighted interconnections. Neural networks have proven to be a popular tool for knowledge extraction, pattern matching, and classification due to their capability of learning from examples with both linear and nonlinear relationships between the input and output signals. However, ANNs have limited ability to characterize local features, such as discontinuities in curvature, jumps in value or other edges, so these algorithma are not well suited for speckled SAR images. The wavelet transform, on the other hand, is efficient in representing and detecting local features in images due to the spatial and frequency localization properties of wavelet bases [6]. With the detection of local features, an object can be easily recognized. Many new algorithms based on wavelet transform have been developed to solve SAR image segmentation problems [7,8]. However, the feature-matching of these algorithms have some shortcomings. In order to ensure the reliability of the matching results, they all require an enormous number of scales to construct the time-frequency features at various scales during the classification process. Each scale corresponds to convolving the signal with a wavelet function; hence a large number of convolutions are needed for these algorithms, which make them computationally inefficient.

In this paper, a wavelet neural network (WNN) method is proposed for SAR image segmentation, which takes full advantages of the partial-resolution characteristic of the wavelet transform and the nonlinear mapping behavior of artificial neural networks.

This paper is organized as follows. In the next section, we will describe quadtree interpretation of SAR imagery and its mixture multiscale autoregressive (MMAR) modeling. In Section 3, we will propose a hybrid method based on the GA algorithm and EM algorithm for the MMAR model. In Section 4, we will present the experimental results. In Section 5, we will present a short conclusion concerning our algorithm.

2. Wavelet Transform

In signal analysis applications, it is necessary to extract signal features with Fourier transformation, but it is only a time domain transform, which has no time-frequency localization features. The theory of wavelet transformation was first proposed in the field of multi-resolution analysis; among others, it has been applied to image and signal processing. A continuous wavelet transform (CWT) can decompose a signal into a set of finite basis functions, which can uncover transient characteristics in the signal. Wavelet analysis is the breaking up of a signal into dilations and translation versions of the original wavelet, referred to as the mother wavelet. The wavelet must be oscillatory, have amplitudes that quickly decay to zero, and have at least one vanishing moment. Wavelet coefficients W_x(a, τ) are produced through the convolution of a mother wavelet function ψ(t) with the analyzed signal f(t), it is:

W_{x} (a, τ) = \frac{1}{\sqrt{| a |}} \int f (t) ψ (\frac{t - τ}{a}) dt

(1)

where a and τ denote the scale and local center of the analyzing wavelet. By adjusting the scale, a, a series of different frequency components in the signal can be obtained.

Several wavelet mother functions have been proposed in the wavelet theory. Each mother function has its suitable application. In this work, the wavelet employed is the Morlet Wavelet, due to its directional selectiveness capability of detecting oriented features, fine tuning to specific frequencies and its good localization in time and frequency [9]. This is a sinusoidal signal modulated by a Gaussian wave. It is characterized for its narrow frequency response, which offers a higher spectral resolution than the Mexican Hat wavelet. This wavelet is particularly useful for filtering out the background noise of the images. In this paper, the Morlet wavelet is applied as:

ψ_{a, τ} (t) = e^{- t_{i}^{2} / 2} cos (5 t_{i})

(2)

where

t_{i} = \frac{t - τ}{a}

3. Wavelet Neural Network

3.1. Structure of Wavelet Neural Network

The reason for the application of WNN in case of such a problem as classification is that the feature extraction and representation properties of the wavelet transform are merged into the structure of the ANN to further extend the ability to approximate complicated patterns.

The WNN can be considered an expanded perceptron in which the neurons of the first layer are replaced by wavelet nodes [10,11]. The wavelet nodes allow the detection of the transient, as well as the extraction and selection of a small number of meaningful features; the obtained features are then regarded as inputs to the subsequent neurons used as a classifier.

The WNN employed in this paper is designed as a three-layer structure with an input layer, a wavelet layer, and an output layer. The topological structure of the WNN is illustrated in Figure 1. In this WNN model, the hidden neurons have wavelet activation functions of different resolutions and ω_i is the weight connecting the hidden layer and output layer. For an input vector x = [x₁, x₂, …., x_n], the output of the i th wavelet layer neuron is described as follows:

ψ_{k} (x) = \sum_{i = 1}^{n} exp (- {(\frac{x_{i} - d_{k}}{t_{k}})}^{2} / 2) cos (5 \cdot \frac{x_{i} - d_{k}}{t_{k}})

(3)

where x_i is the i th input vector and k is the number of wavelet node. d_k and t_k are translation parameter and the dilation parameter, respectively.

The output of the third layer is the weighted sum of ψ_k(x)

y (x) = \sum_{m = 1}^{k} ω_{m} ψ_{m} (x)

(4)

3.2. Training of WNN

Wavelet network training consists of minimizing the usual least-squares cost function:

E = \frac{1}{2} \sum_{j = 1}^{s} {(y_{j} - o_{j})}^{2}

(5)

where s is the number of training samples for each class and o_j is the optimal output of the j th input vector.

Due to the fact that wavelets are rapidly vanishing functions, a wavelet may be too local if its dilation parameter is too small and it may sit out of the domain of interest if the translation parameter is not chosen appropriately.

Therefore, it is inadvisable to initialize the dilations and translations randomly, as is usually the case for the weights of a standard neural network with sigmoid activation function. We use the following initialization procedure, setting.

The same value to dilation parameter d_k is given randomly at the beginning, and initializing the translation parameter t_k is as follows:

t_{k} = (k \times s) / K, k = 0, 1, 2 \dots K - 1

(6)

where s is the number of training samples for each class and K is the number of nodes in the wavelet layer.

The partial derivative of parameters d, t, ω are as follows:

\begin{array}{l} \frac{\partial E}{\partial d_{m}} & = \sum_{j = 1}^{s} 2 (y_{j} - o_{j}) \cdot (\sum_{m = 1}^{k} ω_{m} exp (- {(\frac{x - d_{m}}{t_{m}})}^{2} / 2) ((\frac{x - d_{m}}{t_{m}^{2}}) cos (5 \cdot \frac{x - d_{m}}{t_{m}}) + \frac{5}{t_{m}} sin (5 \cdot \frac{x - d_{m}}{t_{m}}))) \\ = \sum_{j = 1}^{s} 2 (y_{j} - o_{j}) \cdot (\sum_{m = 1}^{k} ω_{m} exp (- \frac{s_{m}^{2}}{2}) \frac{(s_{m} cos (5 s_{m}) + 5 sin (5 s_{m}))}{t_{m}}) \end{array}

(7)

\begin{array}{l} \frac{\partial E}{\partial t_{m}} & = \sum_{j = 1}^{s} 2 (y_{j} - o_{j}) \cdot (\sum_{m = 1}^{k} ω_{m} exp (- {(\frac{x - d_{m}}{t_{m}})}^{2} / 2) \frac{x - d_{m}}{t_{m}^{2}} (\frac{x - d_{m}}{t_{m}} cos (5 \cdot \frac{x - d_{m}}{t_{m}}) + 5 sin (5 \cdot \frac{x - d_{m}}{t_{m}}))) \\ = \sum_{j = 1}^{s} 2 (y_{j} - o_{j}) \cdot (\sum_{m = 1}^{k} ω_{m} exp (- \frac{s_{m}^{2}}{2}) \frac{s_{m}}{t_{m}} (s_{m} cos (5 s_{m}) + 5 sin (5 s_{m}))) \end{array}

(8)

\frac{\partial E}{\partial ω_{m}} = \sum_{j = 1}^{s} \sum_{m = 1}^{k} ψ_{m} 2 (y_{j} - o_{j})

(9)

where

s_{m} = \frac{x - d_{m}}{t_{m}}

We adjust the parameters by the following equation:

Θ^{n} = Θ^{n - 1} - α Δ Θ

(10)

where Θ = (d,t,ω)^T is vector of the parameters d, t and ω, a is learning rate between 0.1 and 0.9.

4. Experiments

We applied the proposed method to two SAR images sized 256 × 256 pixels [Figure 2(a)] to demonstrate the differences between the Morlet and Mexihat procedures; these images include two regions.

First, about 100 samples were selected as the training data. In our experiment, the initial value of t_k is decided by Equation (6)d_k = 2, and ω is a random value between −0.5∼0.5. When the neuron number of the wavelet layer is 25, the segmentation results are best. The compared segmentation results are shown in Figure 2, with (b) showing the Mexihat Wavelet mother function as the transfer function of the second layer and (c) showing the Morlet Wavelet as the transfer function of the second layer. Table 1 is a comparison of the mean square of the above images. It shows that the accuracy ratio of the WNN using Morlet as the transfer function is higher than that of the WNN using Mexihat as the transfer function. Figure 3 is the convergence curve of the WNN training algorithm, which shows the error is almost 10⁻² by the 25^th iteration.

5. Conclusions

In this paper, an effective wavelet neural network for SAR image segmentation is proposed. The method not only has the feature of multiscale analysis, but also has a good performance in classification. Experimental results show that using Morlet as the transfer function is better than using Mexihat. WNN is an effective and accurate method for SAR image segmentation.

Acknowledgments

The authors would like to thank the anonymous reviewers for their detailed comments and questions which improved the quality of the presentation of this paper. This work is supported in part by the National Natural Science Foundation of China (No. 60872064 No. 60375003), the Tianjin Natural Science Foundation (08JCYBJC12300) (08JCYBJC12200), the Science Foundation of Tianjin University of Technology (2006BA15).

References

Chang, C.M.; Liu, T.S. A wavelet network control method for disk drives. IEEE Trans. Contr. Syst. Technol 2006, 14, 63–68. [Google Scholar]
Lin, C.H.; Wang, C.H. Adaptive wavelet networks for power-quality detection and discrimination in a power system. IEEE Trans. Pow. Deliv 2006, 21, 1106–1113. [Google Scholar]
Grimaldi, D. Time-of-flight measurement of ultrasonic pulse echoes using wavelet networks. IEEE Trans. Instrum. Meas 2006, 55, 5–13. [Google Scholar]
Angrisani, L.; Daponte, P. wavelet network-based detection and classification of transients. IEEE Trans. Instrum. Meas 2001, 50, 1425–1435. [Google Scholar]
Subasi, A.; Alkan, A.; Koklukaya, E.; Kiymik, M.K. Wavelet neural network classification of EEG signals by using AR model with MLE preprocessing. Neural Network. IEEE Trans. Pattern Anal. Mach. Intel 2005, 18, 985–997. [Google Scholar]
Mallat, S.; Hwang, W.L. Singularity detection and processing with wavelets. IEEE Trans. Inf. Theory 1992, 38, 617–643. [Google Scholar]
Grandi, G.D.; Lee, J.S.; Schuler, D.L. Target detection and texture segmentation in polarimetric SAR images using a wavelet frame: theoretical aspects. IEEE Trans. Geosci. Remote Sens 2007, 45, 3437–3453. [Google Scholar]
Bustamante, I.T.; Mata, F.S.; Gonzáleza, N.H.; Gazapoa, R.G.; Palaua, J.; Ferreirab, M.M.C. Application of chemometric tools for automatic classification and profile extraction of DNA samples in forensic tasks. Analyt. Chim. Acta 2007, 595, 43–50. [Google Scholar]
Armande, N.; Montesinos, P. Thin nets and crest lines: Application to satellite data and medical images. Compu. Vis. Image Underst 1997, 67, 285–295. [Google Scholar]
Zhang, J.; Walter, G.; Miao, G.Y.; Lee, W.N.W. Wavelet neural networks for function learning. IEEE Trans. Signal Processing 1995, 43, 1485–1497. [Google Scholar]
Delyon, B.; Juditsky, A.; Benveniste, A. Accuracy analysis for wavelet approximations. IEEE Trans. Neural Networks 1995, 6, 332–348. [Google Scholar]

Figure 1. Wavelet Neural Network Structure.

Figure 2. (a) Original SAR image. (b) Segmented image obtained using WNN(Mexihat). (c) Segmented image obtained using WNN(Morlet).

Figure 3. Convergence curve for WNN.

Table 1. Comparison of mean square of the WNN(Mexihat) and WNN(Morlet).

**Table 1.** Comparison of mean square of the WNN(Mexihat) and WNN(Morlet).
	WNN(Mexihat)	WNN(Motlet)
Figure 2 (top)	26.256	21.0044
Figure 2 (bottom)	83.69	66.406

© 2009 by the authors; licensee MDPI, Basel, Switzerland This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Share and Cite

MDPI and ACS Style

Wen, X.-B.; Zhang, H.; Wang, F.-Y. A Wavelet Neural Network for SAR Image Segmentation. Sensors 2009, 9, 7509-7515. https://doi.org/10.3390/s90907509

AMA Style

Wen X-B, Zhang H, Wang F-Y. A Wavelet Neural Network for SAR Image Segmentation. Sensors. 2009; 9(9):7509-7515. https://doi.org/10.3390/s90907509

Chicago/Turabian Style

Wen, Xian-Bin, Hua Zhang, and Fa-Yu Wang. 2009. "A Wavelet Neural Network for SAR Image Segmentation" Sensors 9, no. 9: 7509-7515. https://doi.org/10.3390/s90907509

APA Style

Wen, X. -B., Zhang, H., & Wang, F. -Y. (2009). A Wavelet Neural Network for SAR Image Segmentation. Sensors, 9(9), 7509-7515. https://doi.org/10.3390/s90907509

Article Menu

A Wavelet Neural Network for SAR Image Segmentation

Abstract

1. Introduction

2. Wavelet Transform

3. Wavelet Neural Network

3.1. Structure of Wavelet Neural Network

3.2. Training of WNN

4. Experiments

5. Conclusions

Acknowledgments

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI