Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 864019, 11 pages

http://dx.doi.org/10.1155/2015/864019

## Synthetic Aperture Radar Image Background Clutter Fitting Using SKS + MoM-Based Distribution

^{1}School of Information Engineering, Southwest University of Science and Technology, Mianyang 621010, China^{2}College of Information Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China^{3}Robot Technology Used for Special Environment Key Laboratory of Sichuan Province, Southwest University of Science and Technology, Mianyang, Sichuan 621010, China

Received 10 April 2015; Revised 18 August 2015; Accepted 6 September 2015

Academic Editor: Dan Simon

Copyright © 2015 Zhengwei Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

distribution can accurately model various background clutters in the single-look and multilook synthetic aperture radar (SAR) images and is one of the most important statistic models in the field of SAR image clutter modeling. However, the parameter estimation of distribution is difficult, which greatly limits the application of the distribution. In order to solve the problem, a fast and accurate distribution parameter estimation method, which combines second-kind statistics (SKS) technique with Freitas’ method of moment (MoM), is proposed. First we deduce the first and second second-kind characteristic functions of distribution based on Mellin transform, and then the logarithm moments and the logarithm cumulants corresponding to the above-mentioned characteristic functions are derived; finally combined with Freitas’ method of moment, a simple iterative equation which is used for estimating the distribution parameters is obtained. Experimental results show that the proposed method has fast estimation speed and high fitting precision for various measured SAR image clutters with different resolutions and different number of looks.

#### 1. Introduction

SAR image background clutter statistic modeling [1–6] is one of the basic subjects in the SAR image interpretation such as SAR target detection and discrimination. The research on the subject has important theoretical value and wide application prospects in the automatic target detection and discrimination in high-resolution SAR images [7, 8]. Among the SAR image clutter statistic models, the modeling ability of distribution is good for the heterogeneous clutters, but it cannot fit well for extremely heterogeneous clutters [8]. distribution has a good modeling ability for homogeneous, heterogeneous, and extremely heterogeneous clutters, but the distribution only can be applied to single-look SAR images [9]. distribution, which is proposed by Frery et al. [10], not only has very strong modeling capability for the clutters with wide intensity range and different intensity uniformity but also has very good compatibility for most of clutter models. The distribution can fit various background clutter regions in the single-look and multilook SAR images well and so far it is one of the optimal statistic models in the SAR image background clutter fitting [11], but the distribution also has the defect that the parameter estimation is difficult, which has become a main technical bottleneck in a practical application of the distribution. At present, the typical distribution parameter estimation methods mainly have Freitas’ moment method [1] and Frery’s moment estimation method [10]. The shape and scale parameters of distribution are estimated by using the first-order and second-order moments in Freitas’ method, and Frery’s method realizes the parameter estimation by using 1/2- and 1/4-order moments. The former of the above-mentioned two methods is simpler; the latter needs to be solved by using a numerical method and the calculation cost is larger. But both above-mentioned methods have several common disadvantages: a priori knowledge of the SAR image needs to be used for obtaining the parameter equivalent number of looks (ENL) [12], due to the constraint of convergence condition, the full range estimation of distribution parameters cannot be realized [13, 14], and for very homogeneous clutter, the parameter ENL usually is greater than one; the estimation results of the above-mentioned two methods might be wrong [13, 14]. These defects result in that the parameters of distribution obtained by the above-mentioned two methods are not very accurate and the fitting precision of distribution for the clutters is declined. Theoretically the maximum likelihood estimation (MLE) can also be used for estimating the parameters of distribution, but its equations that are used for estimating the distribution parameters are highly nonlinear, which might cause a sharp rise in the calculation and even a wrong estimated result, so the method is rarely applied in the parameter estimation of distribution. In order to solve the distribution parameter estimation problem, by studying and analyzing the modeling principle of distribution, a new, fast, and accurate distribution parameter estimation method, which combines second-kind statistics (SKS) technique with Freitas’ method of moment (MoM), is proposed in this paper. The proposed method starts from Mellin transform, by deducing the first and second second-kind characteristic functions of distribution, the logarithm moments (log-moments), and logarithm cumulants (log-cumulants) corresponding to the characteristic functions, and then, combining with Freitas’ moment estimation method, a simple iterative equation, which is used for estimating the distribution parameters, namely, the equivalent number of looks (ENL), the shape parameter, and the scale parameter, is obtained ultimately. The method not only solves the problems that “the MoM cannot realize the full range estimation of distribution parameters” and “the MLE needs a long iterative time and the accuracy of the iteration is poor,” but, more importantly, the method considers the parameter ENL as an unknown estimated parameter like the shape and scale parameters, which can ensure that distribution has a high fitting precision for the clutters. Finally, the performance of the proposed method is compared and evaluated by using the measured SAR images with different resolutions and different number of looks. Experimental results show that not only the ENL, shape, and scale parameters can be more quickly and accurately estimated by using the proposed method but also the estimated ENL, shape, and scale parameters are more reasonable than other methods. The proposed method has a high fitting precision for the measured SAR image clutter regions with wide changes in the uniformity.

#### 2. New Parameter Estimation Method for Distribution

For a SAR intensity image, the expression of distribution is where , , , and . is the intensity of pixels. is the equivalent number of looks. is the shape parameter, which essentially reflects the uniformity of a measured clutter region. The parameter indicates that distribution can model various clutter regions with wide uniformity changes. The bigger the parameter is, the more uniform the corresponding clutter region is. is the scale parameter, which is related to the average energy of a measured clutter region. The greater the parameter is, the greater the average energy of the corresponding clutter region is.

In order to fundamentally solve the parameter estimation problem of distribution, a new estimation method of the distribution parameters , , and will be proposed in the following. What needs to be pointed out is that all the parameters are estimated based on the intensity form of distribution in the following discussion.

##### 2.1. Second-Kind Statistics Method

By analyzing the modeling principle of distribution, we found that distribution derives from the speckle model and also integrates with some characteristics of the multiplicative model. We know that the speckle noise may be considered as Mellin convolution in the Mellin transform-based second-kind statistics [4–6], which coincides with the SAR image multiplicative model, so it may be expected that, by imitating the solving ideas that “the first and second second-kind characteristic functions may be generated by Fourier transform and then the cumulants can be obtained by taking the derivative of the characteristic functions at zero” [15], we can derive a new parameter estimation method of distribution by using Mellin transform, which can greatly reduce the difficulty of distribution parameter estimation.

For a function with the domain , the following integral is defined as Mellin transform [4–6]:where the two equals signs with a Delta over them denote a definition formula of Mellin transform given behind them. We know that the domain of the probability density function of a SAR intensity image just is , so Mellin transform may be applied to analyze the SAR intensity image.

Substituting formula (1) into formula (2), we have

According to formula (3) and [16], we may obtain

By carefully observing formula (4), we can find that if we take the logarithm on both sides of the formula, this will greatly simplify the computation:where is called the second second-kind characteristic function corresponding to distribution. By taking the derivative of the function at , we have where , , are called log-cumulants.

Substituting formula (5) into formula (6), we can obtain the log-cumulants corresponding to distribution: where is called digamma function, which is the first-order derivative of logarithm gamma function. is called -order polygamma function. Obviously, formulas (7) of log-cumulants are concise and easy to be computed. More fortunately, thereto the polygamma function is a monotonic function, so formulas (7) are very suitable to be used for estimating parameters , , and of distribution.

Instead of the log-cumulants, usually the sample log-cumulants are used in a practical application. Supposing that , and are grayscale observations of sample pixels of a certain clutter region, then the formulas of sample log-cumulants may be expressed as

Combining formula (7) with formula (8), we can obtain the ultimate estimation equations of distribution parameters , , and :where (9) are called the second-kind statistics (SKS) estimator and the corresponding parameter estimation method is called the SKS method.

##### 2.2. Advantages and Disadvantages of the SKS Method

By analyzing the above-mentioned estimation method of distribution parameters and compared with the moment estimator and the maximum likelihood estimator, we can draw the following conclusions.

The convergence condition of obtaining formula (4) from formula (3) by integral isThe purpose of obtaining formula (4) by integral is to get the log-cumulants corresponding to distribution. We can see from formula (6) that the log-cumulants are only concerned with the derivative of at the point . When , the convergence condition (10) is obviously satisfied. So the SKS estimator successfully solves the problem that “the moment estimator cannot realize the full range estimate of distribution parameters due to the constraint of convergence condition,” which can ensure that distribution has a high fitting precision for the clutter regions with wide uniformity changes.

The moment estimator and the maximum likelihood estimator usually use a priori knowledge of a SAR image to obtain the ENL parameter or take the ISNR (inherent signal to noise ratio) [13, 17] as the ENL parameter of a SAR image. In the above-mentioned two methods to obtain parameter , the partial correlation between the samples that participate in average processing will bring the bias to the estimated parameter , which will result in the decline of the fitting precision. But the ENL parameter is considered as an unknown estimated parameter like the shape and scale parameters in the SKS estimator, so it solves the problem of the estimation bias of the ENL parameter in the moment estimator and the maximum likelihood estimator well.

We can easily see that when , it has a strong regularity between the last two equations of (9) and moreover the polygamma functions that they contain are monotonic, so the solution of the equations, namely, the estimates of the shape, scale, and ENL parameters, can be quickly and accurately obtained by using a numerical iteration method, which overcomes the defects that the maximum likelihood estimator needs a long iteration time and has a poor iteration accuracy.

When a clutter region is very homogeneous, the estimates of distribution parameters are very difficult to be obtained by using the SKS method.

Usually, a good fitting effect can be achieved by the SKS method. But in practice we found that when a clutter region is very homogeneous, the solution of (9) is very difficult to be obtained. As mentioned above, it has a strong regularity between the last two equations of (9), and moreover the function is monotonically decreasing and is close to zero with increasing of in its definition domain and the function is monotonically increasing and is close to zero with increasing of in its definition domain, so the solution of the equations, namely, the estimates of parameters , , and , can be quickly and accurately determined by using a numerical method. For (9), letwhere represents the inverse function of . It may be proved that the function is monotonically decreasing with increasing of and the function is monotonically increasing with increasing of . So in theory, the intersection of the curves and just is the estimates of parameters and that we want to obtain, and then the estimate of parameter can be obtained by substituting the estimates of parameters and into (9). However, when a clutter region is very homogeneous, the curves and would never intersect; that is to say, the solution of (9) is very difficult to be obtained. By analyzing (9) and formula (11), we can find that when a clutter region is very homogeneous, the value of will be very large and the functions and will approach zero. When , both and are constants; if , then the curves and would never intersect; that is to say, the solution of (9) is very difficult to be obtained.

So we can conclude from the above-mentioned analysis that when a clutter region is very homogeneous, the estimates of distribution parameters will be difficult to be obtained by the SKS method. It is this reason that prompts us to study a new distribution parameter estimation method.

##### 2.3. New SKS + MoM Method

The ENL is considered as a parameter to be estimated in the new method. Combining the moment estimation method proposed by Freitas and the first-order logarithm cumulant that is derived from Mellin transform, a new distribution parameter estimation method may be obtained as follows:

Using the first-order origin moment of distribution, we can obtain

Substituting formula (13) and the third formula of (12) into the first equation of (12) and simultaneously considering the constraints , , and , the solution of can be obtained by a numerical method. Substituting the solution into the third formula of (12), we can get the numerical solution of . And finally substituting solution into the second formula of (12), we can get the numerical solution of . So parameters , , and of distribution are obtained quickly and accurately by using the new method. The new method only uses first-order and second-order origin moments and a first-order log-cumulant. In addition, we can see in the following experiment that the new method is more robust than the SKS method.

What needs to be pointed out is that just as Freitas’ moment method the shape parameter in (12) still exists in the constraint , which limits the application of the new method in the fitting of extremely heterogeneous clutters. But the main purpose of the new method is to solve the distribution parameter estimation problem when the clutters are very homogeneous (at this time, ) and the SKS method is not applicable. When the clutters are so extremely heterogeneous that their actual distribution parameter , we still use the SKS method to estimate the parameters of distribution.

The proposed SKS + MoM method not only has the advantage of the SKS method that can quickly and accurately estimate the parameters of distribution but also overcomes the disadvantage of the SKS and MoM methods that when a clutter region is very homogeneous, both the SKS method and the MoM method are not applicable.

#### 3. Experimental Results and Analysis

The processing flow of the proposed method is shown in Figure 1.