Journal of Electrical and Computer Engineering

Journal of Electrical and Computer Engineering / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 382979 | https://doi.org/10.1155/2012/382979

Graham V. Weinberg, Daniel Finch, "Analysis of a Pareto Mixture Distribution for Maritime Surveillance Radar", Journal of Electrical and Computer Engineering, vol. 2012, Article ID 382979, 6 pages, 2012. https://doi.org/10.1155/2012/382979

Analysis of a Pareto Mixture Distribution for Maritime Surveillance Radar

Academic Editor: Igor Djuroviฤ‡
Received24 Apr 2012
Accepted06 Aug 2012
Published05 Sep 2012

Abstract

The Pareto distribution has been shown to be an excellent model for X-band high-resolution maritime surveillance radar clutter returns. Given the success of mixture distributions in radar, it is thus of interest to consider the effect of Pareto mixture models. This paper introduces a formulation of a Pareto intensity mixture distribution and investigates coherent multilook radar detector performance using this new clutter model. Clutter parameter estimates are derived from data sets produced by the Defence Science and Technology Organisation's Ingara maritime surveillance radar.

1. Introduction

In recent years at the Australian Defence Science and Technology Organisation (DSTO), clutter mixture models have been analysed, since they can provide tighter fits to spiky clutter in the distributionโ€™s upper tail region [1, 2]. The main application has been to improve the fit of the K-Distribution to real data, resulting in the KK-Distribution introduced in [1]. This mixture model involves two separate K-Distributions combined in the amplitude domain. The justification for this is that one component of the mixture models the Bragg/whitecap scatterers, while the other accounts for sea spikes. It was also shown in [1] that the KK-Distribution modelled real polarimetric high-resolution sea clutter returns better than a single K-Distribution, for the case of horizontally polarised clutter, which tends to be spikier than vertically polarised returns.

Recently, the Pareto distribution has been proposed as an alternative clutter model and has also been validated for high-resolution radar clutter returns [3, 4]. This intensity model has been found to fit both low and high grazing angle clutter returns, for X-band radar operating in a maritime surveillance environment. It has been found to match the performance of the K-Distribution and to be almost as accurate as the KK-model. There are examples where it matches the performance of the KK-Distribution exactly [4]. The advantage of using a Pareto clutter model is that it is simpler mathematically, resulting in more basic and efficient detection schemes [5]. Thus it is of interest to examine whether, in comparison to the success of the KK-Distribution, a Pareto mixture model yields any improvement in clutter modelling and detection performance.

The paper is structured by introducing the Pareto distribution in Section 2 and illustrating why it is becoming of interest to the radar signal processing community. Section 3 proposes a Pareto mixture model and describes the result of fitting it to DSTOโ€™s high-resolution radar clutter sets. Section 4 derives the Neyman-Pearson optimal detector, based upon this clutter model. In addition, a generalised likelihood ratio test (GLRT) detector is produced. Finally, Section 5 gives some examples of detector performance curves, with parameters estimated from real clutter returns.

2. The Pareto Distribution

As maritime surveillance radar resolution has improved over the years, the backscattering from the sea surface, known as sea clutter, has been found to deviate significantly from the traditional Gaussian model [6]. As a consequence of this, many models of clutter have been proposed and analysed, including the Lognormal [7], Weibull [8] and the K-Distribution [9, 10]. The KK-Distribution [11, 12] was proposed to see whether more accurate clutter fits could be achieved, especially for the very spiky horizontally polarised clutter returns. This model proved to be very successful but has a number of serious shortcomings. The first is that the model is a 5-parameter family, although [1] reduces this to 4, based upon DSTO X-band sea clutter data. Previous models studied in radar signal processing have been two parameter distributions. The complexity of the K- and KK-Distribution clutter models also results in difficult mathematical expressions for the coherent multilook detectors. A good example of this is the detector introduced in [13], which shows the KK-Distribution coherent multilook detector is analytically difficult to work with.

The Pareto distribution [14, 15] is named after the Italian economist Vilfredo Pareto (15 July 1848โ€“19, August 1923) [16, 17] and is a power law probability distribution that has been found to be an excellent model of long-tailed phenomena [18]. Its original application was in the modelling of income over a population [19]. It has been used in the modelling of actuarial data; an example is in excess of loss quotations in insurance [20]. Its usefulness as a model for long-tailed distributions has resulted in applications to a number of diverse areas including physics, hydrology, and seismology [18]. It has also found suitable applications in engineering, such as internet teletraffic modelling [21]. From a radar perspective, its simplicity has resulted in the construction of much simpler detection schemes [5, 22, 23].

In the intensity domain, the Pareto distribution has a power law density given by ๐‘“๐›ผp,๐›ฝp๐›ผ(๐‘ก)=p๐›ฝ๐›ผpp๐‘ก๐›ผp+1,(1) for ๐‘กโ‰ฅ๐›ฝp. The parameter ๐›ผp>0 is called the distributionโ€™s shape parameter, which governs the tail behaviour. Parameter ๐›ฝp>0 is the scale parameter, which determines where the distributionโ€™s support begins. Its mean is given by ๐”ผ(๐‘‹)=๐›ผp๐›ฝp/(๐›ผpโˆ’1), while its variance is var(๐‘‹)=๐›ฝp2๐›ผp/(๐›ผpโˆ’1)2(๐›ผpโˆ’2). Hence for finite first and second moments, it is required that ๐›ผp>2.

Figure 1 shows an example of correlated Pareto intensity clutter, generated using the compound Gaussian model defined in [5]. The figure shows the clutter in both the time and frequency domain. In this case, the Pareto shape parameter is ๐›ผp=25, while its scale parameter is ๐›ฝp=1. The compound Gaussian model has a clutter covariance matrix ฮฃ that is Toeplitz, with factor ๐œ…=0.01, resulting in weak clutter returns. This means the (๐‘–,๐‘—)th entry of ฮฃ is given by ๐œ…|๐‘–โˆ’๐‘—|, for ๐‘– and ๐‘— in the set {1,2,โ€ฆ,๐‘} with ๐‘ the number of rows/columns of ฮฃ.

Figure 2 shows the same Pareto distribution simulated with stronger clutter returns (Toeplitz factor ๐œ…=0.99). In this case we see the correlation having a significant effect on the clutter returns.

These figures show the difficulty of performing target detection in a spiky clutter environment. Both Figures 1 and 2 illustrate that many small targets would easily be missed because of the spikiness of the clutter. In addition, false detections are likely for the same reason. Hence it is important to try to find detectors which allow increased probability of detection, even if there are only minor gains in performance. Small gains are significant for detection of small targets in spiky clutter environments.

3. A Pareto Mixture Model

The extension of the Pareto model (1) to a mixture model can be based upon the definition of the KK-Distribution in [1]. The Pareto mixture distribution, in the intensity domain, is defined by ๐‘“(๐›ผpp1,๐›ผpp2,๐›ฝpp,๐‘˜)(๐›ผ๐‘ก)=(1โˆ’๐‘˜)pp1๐›ฝpp๐›ผpp1๐‘ก๐›ผpp1+1๐›ผ+๐‘˜pp2๐›ฝpp๐›ผpp2๐‘ก๐›ผpp2+1,(2) where ๐›ฝpp is the distributionโ€™s scale parameter (so that ๐‘กโ‰ฅ๐›ฝpp) and ๐›ผpp1 and ๐›ผpp2 are two shape parameters. These three parameters are all assumed to be positive real numbers. The mixing coefficient is ๐‘˜โˆˆ[0,1]. The standard Pareto distribution density (1) is recovered by setting ๐‘˜ to zero or one in (2). In contrast to the KK-Distribution, the Pareto mixture was chosen to have a single-scale parameter, and two shape parameters. The idea behind this was that since the scale parameter determines the support, it would be useful to have only one such parameter and allow for greater distributional fitting through two shape parameters. In MATLAB, an algorithm was created to optimise the fit of the Pareto mixture to real data, in a similar vein to MATLABโ€™s automatic Pareto fit function gpfit, using maximum likelihood estimation. This algorithm fitted the density (2) to individual data sets based upon the four distributional parameters (๐›ผpp1,๐›ผpp2,๐›ฝpp, and ๐‘˜).

Figure 3 shows an example of the fit of a Pareto mixture distribution to real data, using empirical distribution functions. The data used for this purpose is the DSTO Ingara maritime sea clutter trials data, which is fully polarised high-resolution X-band clutter returns. Details of the data and radar can be found in [1, 24]. This data was collected during a trial in 2004, in the Southern Ocean, and has been well documented. A description of the radar and clutter analysis can be found in [1], while the fitting of the Pareto distribution to this data is described in [4]. Key points about the trial are that the radar operated in a circular spotlight mode, scanning the same patch of sea surface at different azimuth angles. The Ingara radar is a fully polarimetric, operating in X-band [24]. No targets were present in the trial, so that the returns are pure clutter.

This figure shows the Ingara data corresponding to run 34690, at azimuth angle of 255โˆ˜, with horizontal polarisation, as considered in [4]. The fitted Pareto parameters are ๐›ผp=3.191 and ๐›ฝp=0.0047, while the Pareto mixture parameters are ๐›ผpp1=3.6637,๐›ผpp2=7.0837,๐›ฝpp=0.0083, and ๐‘˜=0.5402. For comparison, a K- and KK-Distribution were also fitted, using the parameter estimation guidelines in [1]. The fitted K-Distribution has scale parameter ๐‘=74.45 and shape parameter ๐œˆ=2.863, while the KK-Distribution has scale parameters ๐‘1=71.45 and ๐‘2=344.39, the same shape parameter and mixing coefficient ๐‘˜=0.01. In this case the K- and KK-Distributions are extremely close together, while the Pareto mixture (denoted PP in Figure 3) fits the data more closely. Under magnification, the KK-Distribution is closer to the data than the K-Distribution, especially in the upper tail region. The Pareto mixture clearly fits the data better than a single Pareto distribution. This was observed in all the Ingara data examined, implying that the Pareto mixture model gave a better fit to the data.

4. Detector Performance

Neyman-Pearson detectors [15] can be constructed using the Pareto mixture distribution as an intensity clutter model. The formulation of the basic detection problem is outlined briefly, following the set up in [5]. The radar return is denoted ๐ณ, which is a vector of length ๐‘, corresponding to the number of looks. The coherent multilook detection problem is stated as the statistical hypothesis test ๐ป0โˆถ๐ณ=๐œ against ๐ป1โˆถ๐ณ=๐‘…๐ฉ+๐œ, where all complex vectors are ๐‘ร—1. The vector ๐ฉ is the Doppler steering vector, whose components are given by ๐ฉ(๐‘—)=๐‘’๐‘—2๐œ‹๐‘–๐‘“๐ท, for ๐‘—โˆˆ{1,2,โ€ฆ,๐‘}, where ๐‘“๐ท is the normalised target Doppler frequency (โˆ’0.5โ‰ค๐‘“๐‘‘โ‰ค0.5) and is assumed to be completely known. The complex random variable ๐‘… represents the target, and |๐‘…| is its amplitude. The clutter vector ๐œ will be modelled as a spherically invariant random process (SIRP). This means we assume that the clutter takes the form ๐œ=๐‘†๐“–, where ๐“– is a zero mean complex Gaussian process with covariance matrix ฮฃ and ๐‘† is a univariate nonnegative random variable with density ๐‘“๐‘†. Useful references on the SIRP approach include [25, 26]. As in [5], it will be assumed that ฮฃ is semi-positive definite, implying that a Cholesky factorisation can be applied to its inverse. Consequently, there exists a matrix ๐ด such that ฮฃโˆ’1=๐ด๐ป๐ด, where ๐ด๐ป is the Hermitian transpose. This enables a whitening approach to be applied to the original hypothesis test, since SIRPs are unaffected by a linear transformation [25]. As a result of this, the original hypothesis test can be re-expressed as ๐ป0โˆถ๐ซ=๐ง versus ๐ป1โˆถ๐ซ=๐‘…๐ฎ+๐ง, where ๐ซ=๐ด๐ณ,๐ง=๐ด๐œ and ๐ฎ=๐ด๐ฉ. The advantage of this is that the clutter process ๐ง=๐‘†๐ด๐“– is, when conditioned on ๐‘†, a complex Gaussian process because ๐ด๐“– is multidimensional complex Gaussian process, with zero mean but covariance the ๐‘ร—๐‘ identity matrix. This enables the corresponding densities to be written down more simply. In particular, since ๐‘“๐งโˆฃ๐‘†(๐ฑ)=(1/๐œ‹๐‘๐‘ 2๐‘)๐‘’โˆ’๐‘ โˆ’2โ€–๐ฑโ€–2, it can be shown that ๐‘“๐ง1(๐ฑ)=๐œ‹๐‘โ„Ž๐‘๎€ทโ€–๐ฑโ€–2๎€ธ,(3) where โ„Ž๐‘๎€œ(๐‘ฅ)=โˆž0๐‘ โˆ’2๐‘๎‚€โˆ’๐‘ฅexp๐‘ 2๎‚๐‘“๐‘†(๐‘ )๐‘‘๐‘ (4) is called the characteristic function [25].

For the case of a fixed target model (so that ๐‘… is constant), the density under ๐ป1 can be shown to be ๐‘“๐ป11(๐ซ)=๐œ‹๐‘โ„Ž๐‘๎€ทโ€–๐ซโˆ’๐‘…๐ฎโ€–2๎€ธ.(5) To extend this to the case where ๐‘… is unknown the GLRT uses the approximation: โ€–๐ซโˆ’๐‘…๐ฎโ€–2โ‰ˆโ€–๐ซโ€–2โˆ’||๐ฎ๐ป๐ซ||2โ€–๐ฎโ€–2(6) in the density under ๐ป1.

In order to examine the performance of detectors based upon a Pareto mixture distribution assumption, it is necessary to specify the SIRP that generates the desired marginal intensity distributions. To generate a Pareto mixture SIRP, one uses the characteristic function: โ„Ž๐‘๎€ท(๐‘)=(1โˆ’๐‘˜)๐œ๐‘,๐›ผpp1,๐›ฝpp๎€ธ๎€ท,๐‘+๐‘˜๐œ๐‘,๐›ผpp2,๐›ฝpp๎€ธ,,๐‘(7) where ๐œ is defined by ๐›ฝ๐œ(๐‘,๐›ผ,๐›ฝ,๐‘)=๐›ผฮ“(๐‘+๐›ผ)[]ฮ“(๐›ผ)๐‘+๐›ฝ๐‘+๐›ผ.(8) It can be shown that the SIRP with characteristic function (7) has Pareto mixture marginal intensity distributions. The Neyman-Pearson optimal detector [15] is given by the ratio of densities under the respective hypotheses: ๐‘“๐ฟ(๐ซ)=๐ป1(๐ซ)๐‘“๐ป0(๐ซ)๐ป1โ‰ท๐ป0๐œ,(9) where ๐ฟ(๐ซ) is the likelihood function and ๐œ is the detection threshold. The notation (9) means we reject the null hypothesis ๐ป0 (return is pure clutter) if the likelihood exceeds the threshold (set by the probability of false alarm). When we reject ๐ป0 we accept the alternative hypothesis ๐ป1, implying that there is strong evidence of a target present in the radar return ๐ซ. The densities (3) and (5), together with (7), can be applied to (9) to construct the optimal detector explicitly. The GLRT is obtained explicitly by applying (6) to (5) prior to it being applied to (9).

5. Detector Performance Considerations

Using the DSTO Ingara data sets, the performance of the optimal detector, and GLRT, for Pareto mixture clutter models was assessed. Due to the fact that the Ingara data is pure clutter, a synthetic Gaussian target model has been used throughout, as in [5]. Pareto and Pareto mixture clutter models have had their parameters estimated from Ingara data sets. The detector performance curves have then been generated using Monte Carlo simulation, using 106 runs in each case.

The comparison of detectors based upon the Pareto and Pareto mixture clutter models produced some very interesting results. In most cases it was observed that the mixture distribution assumption did not improve detector performance significantly. This is illustrated in Figure 4, which is a detection performance, plotting the probability of detection (Pd) as a function of the target signal to clutter ratio (SCR). In this example, the mixture distribution has ๐›ผpp1=20,๐›ผpp2=100,๐›ฝpp=0.05, and ๐‘˜=0.8. The single Pareto distribution has ๐›ผp=2.347 and ๐›ฝp=0.001277. Additionally, ๐‘=10 looks and the normalised Doppler frequency is ๐‘“๐ท=0.5. The false alarm probability has been set to 10โˆ’6. The simulated clutter has been generated with a Toeplitz covariance matrix with ๐œ…=0.4. This example corresponds to a typical horizontally polarised Ingara data set. Figure 4 shows the optimum detector result as well as the GLRT detector. As can be observed, although the mixture model provided a better clutter fit, there is only a minor gain in detection performance.

Figure 5 is a second example of detector performance. Here, the same clutter parameters and normalised Doppler frequency have been used as for Figure 4. However, the number of looks is ๐‘=5 and the Toeplitz factor has been increased to 0.7. Additionally, the false alarm probability has been reduced to 10โˆ’7. Comparing the pairs of decision rules, we again observe that the mixture model does not provide a significant improvement in detection performance.

Varying the number of looks, normalised Doppler frequency and Toeplitz covariance factor did not significantly alter the results. Hence it appears that although the KK-mixture distribution provides significant detection performance relative to that provided by a detector based upon a single K-Distributed clutter assumption [11], the same is not repeated when the Pareto assumption is applied. This implies a single Pareto distribution is sufficient from a detection performance perspective.

6. Conclusions

Although a Pareto mixture model can provide a better fit to the upper tail region in real radar clutter returns, this has been observed to not significantly improve detection performance, which is the primary application of radar. Increasing the number of degrees of freedom in a model adds computational cost, slowing the radarโ€™s ability to perform its operation in real time. Hence it is computationally better to use a nonmixture model when designing radar detection schemes, under the basic Pareto clutter model assumption.

References

  1. Y. Dong, โ€œDistribution of X-band high resolution and high grazing angle sea clutter,โ€ DSTO Report RR-0316, 2006. View at: Google Scholar
  2. G. V. Weinberg, โ€œNeyman-pearson optimal and suboptimal detection for signals in general clutter mixture distributions,โ€ DSTO Report RR-0363, 2011. View at: Google Scholar
  3. M. Farshchian and F. L. Posner, โ€œThe pareto distribution for low grazing angle and high resolution X-band sea clutter,โ€ in Proceedings of the IEEE International Radar Conference 2010, RADAR 2010, pp. 789โ€“793, May 2010. View at: Publisher Site | Google Scholar
  4. G. V. Weinberg, โ€œAssessing Pareto fit to high-resolution high-grazing-angle sea clutter,โ€ Electronics Letters, vol. 47, no. 8, pp. 516โ€“517, 2011. View at: Publisher Site | Google Scholar
  5. G. V. Weinberg, โ€œCoherent multilook radar detection for targets in pareto distributed clutter,โ€ Electronics Letters, vol. 47, no. 14, pp. 822โ€“824, 2011. View at: Google Scholar
  6. N. Levanon, Radar Principles, Wiley, New York, NY, USA, 1988.
  7. G. V. Trunk and S. F. George, โ€œDetection of Targets in Non-Gaussian Sea Clutter,โ€ IEEE Transactions on Aerospace and Electronic Systems, vol. 6, no. 5, pp. 620โ€“628, 1970. View at: Google Scholar
  8. D. C. Schleher, โ€œRadar detection in weibull clutter,โ€ IEEE Transactions on Aerospace and Electronic Systems, vol. AES-12, no. 6, pp. 736โ€“743, 1976. View at: Google Scholar
  9. E. Jakeman and P. N. Pusey, โ€œA model for non-rayleigh sea echo,โ€ IEEE Transactions on Antennas and Propagation, vol. AP-24, no. 6, pp. 806โ€“818, 1976. View at: Google Scholar
  10. Y. Dong, โ€œOptimal coherent radar detection in a K-distributed clutter environment,โ€ IET Radar, Sonar and Navigation, vol. 6, no. 5, pp. 283โ€“292, 2012. View at: Google Scholar
  11. L. Rosenberg, D. J. Crisp, and N. J. Stacy, โ€œAnalysis of the KK-distribution with medium grazing angle sea-clutter,โ€ IET Radar, Sonar and Navigation, vol. 4, no. 2, pp. 209โ€“222, 2010. View at: Publisher Site | Google Scholar
  12. Y. Dong and B. Haywood, โ€œHigh grazing angle X-band sea clutter distributions,โ€ in Proceedings of the RADAR 2007—The Institution of Engineering and Technology International Conference on Radar Systems, October 2007. View at: Publisher Site | Google Scholar
  13. G. V. Weinberg, โ€œCoherent multilook radar detection for targets in KK-distributed clutter,โ€ in Digital Communication, C. Palanisamy, Ed., chapter 9, Intech, Rijeka, Croatia, 2012. View at: Google Scholar
  14. M. Evans, N. Hastings, and B. Peacock, Statistical Distributions, Wiley, New York, NY, USA, 3rd edition, 2000.
  15. G. P. Beaumont, Intermediate Mathematical Statistics, Chapman and Hall, London, UK, 1980.
  16. G. Bruno, โ€œPareto, Vilfredo,โ€ in The New Palgrave: A Dictionary of Economics, vol. 5, pp. 799โ€“804, 1987. View at: Google Scholar
  17. R. Cirillo, The Economics of Vilfredo Pareto, Frank Cass Publishers, London, UK, 1978.
  18. I. B. Aban, M. M. Meerschaert, and A. K. Panorska, โ€œParameter estimation for the truncated pareto distribution,โ€ Journal of the American Statistical Association, vol. 101, no. 473, pp. 270โ€“277, 2006. View at: Publisher Site | Google Scholar
  19. A. V. Asimit, E. Furman, and R. Vernic, โ€œOn a multivariate pareto distribution,โ€ Insurance, vol. 46, no. 2, pp. 308โ€“316, 2010. View at: Publisher Site | Google Scholar
  20. M. Rytgaard, โ€œEstimation in the pareto distribution,โ€ ASTIN Bulletin, vol. 20, no. 2, pp. 201โ€“216, 1990. View at: Google Scholar
  21. E. Chlebus and R. Ohri, โ€œEstimating parameters of the pareto distribution by means of zipf's law: application to internet research,โ€ in Proceedings of the IEEE Globecom, pp. 1039โ€“1043, December 2005. View at: Publisher Site | Google Scholar
  22. X. Shang and H. Song, โ€œRadar detection based on compound-gaussian model with inverse gamma texture,โ€ IET Radar, Sonar and Navigation, vol. 5, no. 3, pp. 315โ€“321, 2011. View at: Publisher Site | Google Scholar
  23. K. J. Sangston, F. Gini, and M. S. Greco, โ€œNew results on coherent radar target detection in heavy-tailed compound-gaussian clutter,โ€ in Proceedings of the IEEE International Radar Conference, pp. 779โ€“784, May 2010. View at: Publisher Site | Google Scholar
  24. N. J. S. Stacy, D. Crisp, A. Goh, D. Badger, and M. Preiss, โ€œPolarimetric analysis of fine resolution X-band SAR sea clutter data,โ€ in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '05), pp. 2787โ€“2790, July 2005. View at: Publisher Site | Google Scholar
  25. M. Rangaswamy, D. Weiner, and A. Oeztuerk, โ€œNon-Gaussian random vector identification using spherically invariant random processes,โ€ IEEE Transactions on Aerospace and Electronic Systems, vol. 29, no. 1, pp. 111โ€“124, 1993. View at: Publisher Site | Google Scholar
  26. K. Yao, โ€œRadar detection based on compound-Gaussian model with inverse gamma texture,โ€ IEEE Transactions on Information Theory, vol. IT-19, no. 5, pp. 600โ€“608, 1973. View at: Google Scholar

Copyright © 2012 Graham V. Weinberg and Daniel Finch. 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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views2778
Downloads810
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.