Multiresolution Fuzzy Clustering For SAR Image Segmentation
Yuttapong Rangsanseri and Punya Thitimajshima
Department of Telecommunications Engineering, Faculty of Engineering
King Mongkut's Institute of Technology Ladkraband, Bangkok 10520, Thailand
E-mail: Kryuttha@kmitl.ac.th, ktpunya@kmitl.ac.th
Abstract
This paper describes a pyramid-based method for clustering single polarization synthetic Aperture Radar (SAR) images. A pyramid of the image is first constructed by using the wavelet transform. The fuzzy c-means (FCM) algorithm is then applied to the pyramid with the coarse-to-fine approach. Finally, clustering is followed by a majority filtering to obtain more homogeneous regions. The algorithm has been tested on JERS-1/SAR image and the results demonstrate its potential usefulness.
Introduction
The Japanese Earth Resources Satellite-1 (JERS-1) has provided synthetic Aperture Radar (SAR) image operative at 1.275 GHZ since its launch on Feb. 11, 1992. This radar remote sensing provides terrain image regardless of weather and operations time. The image itself is the result of microwave backscattering and is expected to provide different and important information on terrain, which cannot be obtained by any optical sensors.
Remote sensing image segmentation is often accomplished by clustering when ground truth is not available to provide samples to train a supervised classifier[1]. In this technique an image is segmented into unknown classes. It is the task of the user to label those classes afterwards. There are a great number of clustering methods. Many are accomplished by an iterative technique. An approach based on clustering a pyramid data structure has been proposed by Trivedi and Bezdek [2]. In that research, each layer in the pyramid is created by averaging groups of 2x2 pixels to produce a composite pixel of reduced resolution.
The segmentation of SAR images usually uses texture measures, describing the relation between gray levels of neighboring pixels, as feature for the classification algorithm feature, based on the multiresolution decomposition provided by the wavelt transform [3] [4]. The method also produces a pyramid structure similar to [2]. The image segmentation is performed by the multiresolution fuzzy clustering technique.
Multiresolution analysis
Multiresolution analysis [5] permits the analysis of the signal in different frequency bands or at different scales. In practice, mutiresolution analysis is carried out using 2 channel filter banks composed of a low-pass and a high-pass filter and each filter bank is then sampled at a half rate (1/2 down sampling) of the previous frequency.
Wavelet is type of filter family. They are constructed from a mother wavelet [6], which can be any admissible function. A family is constructed by dilating and translating the mother wavelet by different amount. In the case of an image, the filtering is implemented in a separable way by filtering the lines and columns. The low-pass filter of level I provides an approximation of the initial image at the scale of 2'. The high-pass filters give the image of the wavelet coefficients or the high frequency images.
Fuzzy clustering
The fuzzy clustering has been proved to be very well suited to deal with the imprecise nature of geographical information, including remote sensing data. According to the fuzzy clustering framework, each cluster is a fuzzy set, and each pixel in the image has a membership value associated to each cluster, ranging between 0 and 1, measuring how much the pixel belongs to that particular cluster. There have been many different families of fuzzy clustering algorithms proposed in the last decade [7]. The one used in this paper is the Fuzzy C-Means algorithm (FCM), which is an iterative technique based on the minimization of a generalized within group sum of squared (WGSS) err objective function.