Accuracy data for the considered test images corrupted by i.i.d. speckle (single-look case)
Blind estimation of noise characteristics (BENC), such as noise type, its statistics and spectrum, has become an actual practical task for various image processing applications (Vozel et al., 2009). There are several reasons for this. First, noise is one of the main factors degrading and determining the quality of images of different types: grayscale and color optical (Liu et al., 2008; Foi et al., 2007; Plataniotis&Venetsanopoulos, 2000), component images in certain sub-bands of hyperspectral remote sensing data (Aiazzi et al., 2006), radar and ultrasound medical images (Lin et al., 2010; Oliver&Quegan, 2004), etc. Second, information on noise characteristics is valuable and widely exploited in most of stages of image processing. For example, it is used in edge detection for threshold setting (Davies, 2000), image filtering (Touzi, 2002; Lee et al., 2009) including denoising techniques based on orthogonal transforms (Mallat, 1998; Sendur&Selesnick, 2002; Egiazarian et al., 1999), image reconstruction (Katsaggelos, 1991), lossy compression of noisy images (Bekhtin, 2011), non-reference assessment of image visual quality (Choi et al., 2009), etc. Third, although there can be initial assumptions on noise type and a range of variations of its statistical parameters, these parameters can be quite different even for a given imaging system depending upon conditions of its operation. The requirements to information accuracy on noise parameters are rather strict, e.g., variance of pure additive or pure multiplicative noise has to be known or pre-estimated with a relative error not larger than ±20% (Abramov et al., 2004). Thus, it is often desirable to estimate noise characteristics for a given image.
Besides, amount of images offered by various imaging systems increases enormously. Therefore, it becomes difficult to evaluate noise characteristics in an interactive manner since this requires time, perfect skills, and availability of the corresponding software. Moreover, there are practical situations and applications for which it is impossible to find a highly qualified expert to perform the task of evaluation of the noise characteristics. The examples are estimation of noise characteristics in remote sensing images on-board satellites (Van Zyl et al., 2009). BENC can be also useful even if an expert is involved to analysis of the noise characteristics. This happens, e.g., if a newly designed and manufactured imaging system is verified to check do the main properties of the noise present in the formed images conform expected (forecasted) ones. Then, the output estimates of BENC can be compared to the outcomes of the expert analysis and support (control) each other.
There are quite many known methods of BENC designed so far. A few of them can operate on images corrupted by a general type of signal-dependent noise (Liu et al., 2008). Most of known BENC methods are able to deal only with a particular type of noise under assumption that the noise type is known a priori or pre-determined in an automatic manner (Vozel et al., 2009). The case of pure additive noise has been studied more thoroughly in literature (see Vozel et al. 2009; Zoran&Weiss, 2009; Abramov et al., 2008; Lukin et al. 2007, and references therein). Some of these methods can be, after certain modifications, applied to estimation of multiplicative noise variance. These modifications basically relate to either application of logarithmic type homomorphic transform or a special approach to form local estimates of multiplicative noise relative variance as a normalization of local variance estimates by squared local mean (Vozel et al. 2009). However, quite many BENC methods exploit local estimate scatter-plots and line (curve) fitting into them to evaluate multiplicative noise variance (Lee et al., 1992; Ramponi&d’Alvise, 1999). Note that a multiplicative noise is typical for radar imagery, in particular, images acquired by synthetic aperture radars (SARs) where coherent principles of image forming are employed (Solbo&Eltoft, 2004; Oliver&Quegan, 2004). Speckle is a specific noise-like phenomenon arising in formed images and it is known to be the dominant factor degrading their quality (Oliver&Quegan, 2004). For many operations of radar (and ultrasound) image processing, the characteristics of the speckle are to be known in advance or pre-estimated (Lee et al., 2009; Solbo&Eltoft, 2008).
One can argue that there are many practical situations when speckle characteristics such as the (relative) variance of the multiplicative noise (or the efficient number of looks) and the speckle distribution law are known in advance or can be predicted from theory (Oliver&Quegan, 2004). This holds if a given SAR operates in a known mode (e.g. forms one-look amplitude images) and the operation parameters are stable. Then, it is enough to carry out a preliminary analysis of several images acquired by this SAR manually (in interactive mode) to be sure that the aforementioned characteristics (parameters) conform theory and are stable enough.
However, in many practical situations, it is worth applying BENC, sometimes in addition to an interactive analysis. First, suppose that a new SAR is tested and it is desirable to know whether or not it provides the desired (forecasted, expected) characteristics. Second, one might deal with SAR images for which full description of the imaging mode used is absent (Lee et al., 1992; Ramponi&d’Alvise, 1999). Third, although it is assumed that the multi-look mode of image formation allows decreasing the speckle variance by the number of looks, this is not absolutely true and, in practice, noise reduction is not as efficient as ideally predicted (Anfinsen et al., 2009; Foucher et al., 2000).
Therefore, two important questions arise: what is the accuracy of the existing blind estimation methods and what BENC to apply? To our best knowledge, there are no studies dealing with intensive testing of BENC with application to speckle (our conference paper (Lukin et al., 2011) seems to be one of the first attempts in this direction). By intensive testing we mean the use of tens of different images having different content and/or many realizations of speckle for both single and multi-look modes. There are several reasons why such testing has not been carried out yet. The main reason is the absence of the test radar images commonly accepted by the radar data processing community. We have to stress here that it is quite difficult to create test SAR images since one has to find an answer to many particular questions as what terrain and objects to simulate, what model of the carrier trajectory and its instabilities to use, to consider moving objects or not, what is SNR in radar receiver input, what kind of received signal processing is used (Dogan&Kartal, 2010; Di Martino et al., 2012), etc. Another reason is that, maybe, designers of BENC for speckle have been satisfied by accuracy of the obtained estimates for a limited set of processed images and have not tried applying their methods to a wider variety of data.
Experience obtained recently in testing BENCs for additive and signal dependent noise cases (Vozel et al., 2009; Abramov et al., 2011; Lukin et al., 2009b) clearly demonstrates the following. First, whilst a given method can produce an acceptable accuracy for many tested images, there can be a few test images (usually highly textural ones and/or with clipping effects) for which abnormal (unacceptable) estimates are obtained. Just to these images one has to pay more attention in attempts to improve a methods’ performance. Second, a spatial correlation of noise present in most of real life images and often ignored in a design and testing of many BENC techniques can considerably influence an accuracy of estimation methods (Abramov et al., 2008). Recall that a spatial correlation of speckle is a feature typical for SAR images (Solbo&Eltoft, 2008, Lukin et al., 2008; Lukin et al., 2009; Ponomarenko et al., 2011) which is not often taken into account in SAR image simulations.
Thus, we come to a necessity to perform intensive testing of BENC methods without having a set of standard test images. Our idea then is to create a set of test SAR images with a priori known characteristics of the speckle similar to those ones observed in practice. In this sense, TerraSAR-X images can be a good choice (in Section 2, we explain this in detail). Note that quite many of them are now available in the convenient form and their amount is rapidly growing (see http://www.infoterra.de/free-sample-data). Then, it becomes possible to test BENCs for simulated data (Section 3) and to predict what could happen in practice. These predictions are then verified for the considered methods for high quality data provided by TerraSAR-X data (Section 4) to offer practical recommendations on the BENC method selection and setting its parameters. Finally, conclusions follow.
2. Basic properties of speckle and its modeling
Speckle is a typical example for which pure multiplicative model is usually exploited (Touzi, 2002; Oliver&Quegan, 2004). This means that a dependence of signal dependent noise variance on true value is monotonically increasing proportionally to squared (true value). Speckle is not Gaussian and its probability density function (PDF) depends upon a way of image forming (amplitude or intensity) and number of looks (Oliver&Quegan, 2004). PDF of the speckle considerably differs from Gaussian if a single-look imaging mode is used and it is either Rayleigh (for amplitude images) or negative exponential (for intensity images) for the case of fully developed speckle. If multi-look imaging mode is applied, the speckle PDF becomes closer to Gaussian and depends upon the number of looks.
To get an imagination on fully developed speckle PDF, consider real-life data produced by TerraSAR-X imager. Its attractive feature is that data (images) are freely available at the aforementioned site. These data have full description of parameters of the imaging system operation mode used for obtaining each presented image. Large size images (thousands to thousands pixels) for many different areas of the Earth are offered. Furthermore, a brief description of a territory, observed effects and cover types is given. This allows selecting and processing data with different numbers of looks, properties of a sensed terrain, a desired polarization, etc. Fragments of certain size as 512x512 pixels can be easily cut from large size data arrays and studied more thoroughly. Another positive feature is that single-look images are presented in the complex valued form. This allows obtaining single-look images in aforementioned forms (representations). It also makes possible to analyze distributions of real and imaginary part values for image fragments, etc. While considering single-look SAR images in this paper, we use amplitude images since it is the most common form and it provides convenient representation for visual analysis.
One more advantage is that TerraSAR-X is a high quality system designed by specialists from German Aerospace Agency DLR who have large experience in creation and management of spaceborne SAR systems (Herrmann et al., 2005). The TerraSAR-X imager provides a stability of noise characteristics and practical absence of an additive noise in formed images. Later, it will be explained why this is so important for further analysis.
To partly corroborate conformity of theory and practice for single-look SAR images, we have manually cropped several sub-arrays of complex valued data that correspond to homogeneous terrain regions. The histogram of real (in-phase) part values for one such a fragment is presented in Fig. 1(a) where sample data mean is close to zero. The histogram for imaginary (quadrature) part is very similar. Gaussianity tests hold for both data sub-arrays. The histogram of the amplitude single-look image for the same fragment is represented in Fig. 1(b). Since both components of complex valued data are Gaussian with approximately the same variance, the amplitude values obey Rayleigh distribution. This one more time shows that for single-look amplitude images speckle can be simulated as pure multiplicative noise having Rayleigh PDF. Using modern simulation tools as, e.g., Matlab, this can be done easily, at least, for the case of independent identically distributed (i.i.d.), i.e. spatially uncorrelated speckle.
The histogram in Fig. 1(b) also shows one more aspect important for simulations. Speckle image values can be by 3...4 times larger than mean (which is close to in homogeneous image regions). Then, if the values of are modelled as 8-bit data, the noisy values can be outside the limits 0…255 and, therefore, 16-bit representation of a simulated noisy image is to be used to preserve statistics of the speckle. In Section 3, we will show what might happen to the estimates provided by BENCs if clipping effects take place for simulated noisy SAR images, i.e. if they are represented as 8-bit data.
To additionally analyze statistics of the speckle, we have also tested several manually cropped homogeneous regions in different single-look images that correspond to either rather large (about 70x70 pixels) agricultural fields and water surface. The estimated speckle variance ( is an
To consider and simulate speckle more adequately, we have also analyzed spatial correlation of speckle using TerraSAR-X data. This has been done in three different ways. First, standard 2D autocorrelation function (ACF) estimates have been obtained for 32x32 pixels size homogeneous fragments. They have been inspected visually and have demonstrated the absence of far correlation and the presence of essential correlation for neighboring pixels in single-look amplitude images (see example in Fig. 2(a)). It is interesting that even higher correlation for neighboring pixels has been observed for multi-look images (see example in Fig. 2(b)). There are also ACF side lobes for azimuth direction that, most probably, arise due to peculiarities of the SAR response to a point target.
Second, we have analyzed a normalized 8x8 DCT spectrum estimates obtained in a blind manner (Ponomarenko et al., 2010) for several considered images. These estimates also clearly indicate that speckle is spatially correlated, i.e., not i.i.d. (Lukin et al., 2011).
Third, we have also calculated a parameter
and the second estimate is calculated in the DCT domain as
where are DCT coefficients of
To simulate single- and multi-look SAR images, we have used four aerial optical images as (all test images are of size 512x512pixels). These four images are presented in Fig. 3. Positive features of these images allowing to use them in simulation of SAR data are the following. First, these images have practically no self-noise that could later influence blind estimation of speckle statistics.
Second, these are the images of natural scenes and, thus, they contain large quasi-homogeneous regions, edges with different contrasts, various textures and small-sized targets.
Note that we have simulated speckle with the same statistics for all pixels ignoring the fact that for small-sized targets it might differ from speckle in homogeneous image regions. This simplification is explained by the following two reasons. First, more complicated models of speckle are required for small-sized targets. Second, the percentage of pixels occupied by small-sized targets is quite small in real-life images (Lee et al., 2009) and local estimates of noise statistics in the corresponding scanning windows are anyway abnormal. Thus, these local estimates are “ignored” by the BENCs considered below which are robust (see Section 3 for more details).
Fig. 4 gives two examples of noisy test image with fully developed speckle (single look). For one of them (Fig. 4(a)) speckle is i.i.d. whilst for the second (Fig. 4(b)) speckle is spatially correlated (see details of its simulation below). Even visual analysis of these two noisy images allows noticing the difference in speckle spatial correlation. As it will become clear from the visual analysis of real-life SAR images presented later in Section 4, the case shown in Fig. 4(b) is much closer to practice.
Thus, from a practical point of view, it is more reasonable to simulate spatially correlated speckle. This can be done in different ways. In our study, we have employed the following simulation algorithm:
Generate 2D array of a required size for Gaussian zero mean spatially correlated noise (GCN – Gaussian correlated noise) with a desired spatial spectrum (this is a standard task solution of which is omitted).
Transform the 2D GCN data array into 1D array C of size K= in a pre-selected way, e.g., by row-by-row scanning.
Generate 1D array B of i.i.d. Rayleigh distributed unity mean random variables of size K=.
For the array C, form an array of indices CI in such a manner that CI(1) is the index of the element in C which is the largest, СI(2) is the index of the element in C which is the second largest, and so on. Finally, CI(K) is the last element of the array CI which is the index of the smallest element of C.
Similarly, form an index array BI for the array B.
For k=1..K make valid the condition C(CI(i))=B(BI(i)). Then, noise with Gaussian distribution is replaced by noise with the required distribution (Rayleigh in our case).
The obtained array C is transformed to 2D array RES of size K= in the way inverse to it has been done in step 2.
The source code in Matlab that realizes the described algorithm is presented below:
Here M, N correspond to IIm and JIm (that is to the simulated image size), and all other notations are the same as described above. The obtained 2D array has a Rayleigh distribution and has practically the same spatial correlation properties as GCN. Then the values of RES(i,j) are pixel-wise multiplied by to obtain the corresponding speckle values .
The image presented in Fig. 4(b) has been obtained in the way described above. Moreover, this allows getting multi-look images if several realizations of the speckle with desired spectrum are generated and then averaged.
3. Considered blind estimation techniques and their accuracy analysis for simulated data
Describing the considered BENC methods, one should keep in mind that blind estimates of speckle characteristics obtained for a given method can differ from each other due to the following factors:
properties and parameters (if they can be varied or user defined) of a method applied;
method robustness with respect to outliers;
content of an analyzed image;
an observed speckle realization in the considered image;
clipping effects (if they take place).
Because of this, we first describe BENCs used in our studies and the main principles put into their basis. Then, simulation results are presented for simulated single- and multi-look SAR images, and the analysis of these results is performed.
3.1. Considered BENCs
As it has been mentioned in Introduction, there are two basic approaches to blind estimation of . The first approach presumes forming local estimates of speckle variance and robust processing of the obtained local estimates. The second approach is based on obtaining a scatter-plot and robust regression line fitting into it.
Let us start from considering the former approach. It consists of the following stages. At the first stage, an analyzed image is divided into non-overlapping or overlapping blocks and local estimates are obtained as
To understand the operation principles of the first group of methods, it can be useful to look at distributions of the local estimates (3). As examples, the two distributions of local estimates for the two real-life (TerraSAR-X) single-look elementary images (presented in Fig. 5(a) and
Fig. 5(b)) are shown in Figs. 6(a) and 6(b), respectively (
Besides, the distributions in Fig. 6 differ by heaviness of the right-hand tail. Recall that this tail stems from the presence of the so-called “abnormal” local estimates (3) that are obtained in heterogeneous image blocks (Vozel et al., 2009). For the elementary images that have a simpler structure (Figures 5(b) and 7), the tail heaviness is considerably less (see Fig. 7(b)).
The property that the distributions have maxima with the mode close to the true value of has been put into the basis of several BENC methods for estimation of noise variance (Vozel et al., 2009). The task is then to find the distribution mode automatically, robustly and with a high enough accuracy. For this purpose, it is possible to exploit robust mode finders such as a sample myriad, bootstrapping and minimal inter-quantile distance with properly (adaptively) set parameters. Since an improved minimal inter-quantile distance estimator provides the best accuracy (Lukin et al., 2007), we use it in our further studies. The technique based on obtaining the set of local estimates according to (3) and estimation of its mode by the improved minimal inter-quantile distance estimator (Lukin et al., 2007) is further referred as
The second group of BENC methods, as it has been mentioned above, is based on scatter-plots. A traditional way of scatter-plot representation for signal-dependent noise is the following. For each block, a point in Cartesian system is obtained where its Y coordinate corresponds to a local variance estimate and a local mean estimate is its argument (X axis coordinate ). An example of such a scatter-plot for the image in Fig. 5(b) is presented in Fig. 8(a). A curve is depicted in this scatter-plot. It is seen that it goes through the centers of the main clusters of this scatter-plot (the cluster centers are indicated by red squares) where the clusters are formed by normal local estimates (see details below). However, there are also quite many points that are located far away from this curve and cluster centers. These points correspond to abnormal local estimates obtained in heterogeneous blocks. This means that if one presumes to fit a polynomial type curve and then to obtain , where is the parameter of the fitted curve, the method of curve fitting should be robust with respect to outliers.
There are also other ways to obtain a scatter-plot. One variant is that a point Y coordinate corresponds to a local variance estimate and a squared local mean estimate is its argument (X axis coordinate ). An example of such a scatter-plot obtained for the same single-look image is represented in Fig. 8(b). Then one has to fit a curve
i.e. straight line where the estimate , is the parameter of the fitted line. Another option is to obtain a scatter-plot in such a way that a point coordinate Y relates to a local standard deviation estimate where its argument (X axis coordinate ) is the corresponding local mean estimate. Then one has where is the fitted straight line parameter that serves as the estimate . This kind of a scatter-plot is shown in Fig. 8(c) for the same single-look SAR image. Visual analysis of the scatter-plots in Figs. 8(b) and 8(c) shows that for them there are also some clusters of normal local estimates whereas abnormal estimates are present as well. An advantage of the latter two approaches is that it is, in general, simpler to fit a straight line than a higher order polynomial. In particular, there are standard means for this purpose as, e.g., the Matlab version of robustfit method (DuMouchel&O'Brien, 1989). For methods analyzed below, we have used the approach based on (6).
Finally, there are also methods that exploit scatter-plot data to find cluster centers and curve (line) fitting using these scatter-plot centers (Zabrodina et al., 2011; Abramov et al., 2011). Cluster centers are indicated by red color dots in scatter-plots in Fig. 8. The cluster center coordinates relate to where is the estimate of distribution mode of local variance estimates for a
There is the straight line fitted into the cluster centers in Figs. 8(b) and 8(c). In this case, robustness with respect to abnormal local estimates is provided indirectly due to robust methods used for finding cluster centers. However, there can be also abnormal cluster centers. To reduce their influence, special techniques as RANSAC or double weighting (DW) LMSE fit can be applied (Zabrodina et al., 2011). Taking into account the comparison results (Abramov et al., 2011; Zabrodina et al., 2011), below we consider only the DW curve fitting to scatter-plot since this method, on the average, provides the best results. It is possible to use different sizes of blocks for local variance and local mean estimation in blocks. Below we study 5x5 and 7x7 pixel blocks. The technique based on forming a scatter-plot, its division into fixed number of clusters, finding cluster centers using mode estimation and DW line fitting is referred below as
There are also other techniques based on curve fitting into cluster centers with improved robustness with respect to outliers. First, cluster centers can be determined without image pre-segmentation (as for the
3.2. Analysis of simulation results
Let us analyze the obtained simulation results. The main properties and accuracy characteristics of the aforementioned methods based on finding a distribution mode have been intensively studied for the case of additive noise (Lukin et al., 2007). Although the multiplicative noise case is considered here, the conclusions drawn for the additive case might be still valid for
Another conclusion is that the estimation bias (denoted as for the multiplicative noise case) usually contributes more to aggregate error , where denotes the variance of blind estimation of . Here and where notation means averaging by realizations.
Let us check are these conclusions valid for the multiplicative noise case. Usually variance is determined for a large number of realizations of the artificially added noise that corrupts a given test noise-free image. Thus, we have simulated 200 realizations of i.i.d. speckle with Rayleigh distribution. The obtained simulation results are presented in Table 1. Analysis shows that estimation bias is also positive for all four test images and for both studied sizes of blocks. The values of are of the order 10-6. Thus, they are two magnitude order less than squared bias and have negligible contribution to . This shows that, in fact, it is possible to analyze only the estimation bias or even the estimates obtained for only one realization of the speckle. At least, this is possible for the test images of the considered size of 512x512 pixels or larger ( decreases if a processed image size increases).
One more conclusion that follows from data analysis for
Consider now data for
Finally, let us analyse data for
We have also obtained simulation results for 4-look test images corrupted by i.i.d. speckle (theoretical is equal to 0.273/4=0.068). They are the following. For the first test image, estimation bias is 0.0101, 0.0100 and 0.0005 for
Consider now the case of spatially correlated noise. We have carried out preliminary simulations and established that estimation bias contributes considerably more than estimation variance to the . Thus, below we present only the errors determined as the difference between the obtained estimate and the true value for single (only one) realization. The simulation results for single-look images are collected in Table 2.
An interesting observation that follows from data analysis in Table 2 is that the differences are mostly negative, at least, for 5x5 block size, i.e. speckle variance is underestimated. This can be explained as follows. One factor that influences blind estimation is distribution mode position. Normal local estimates in blocks that form this mode are mostly smaller than (Lukin et al., 2011b). Because of this, speckle variance estimates tend to smaller values for
As it can be also seen from analysis of data in Table 2, the estimates for 7x7 blocks are larger than the corresponding estimates for 5x5 blocks. This is because mode position for normal local estimates shifts to right (to larger values) if the block size increases. This effects have been illustrated for spatially correlated speckle (Lukin et al., 2011b) and for spatially correlated additive noise (Abramov et al., 2008). Then, the final estimates for all BENCs also increase. Analysis shows that it is worth using the block size 7x7 pixels for
Finally, simulation results for four-look test images corrupted by spatially correlated speckle are represented in Table 3. The data are presented as - similarly to the previous case =0.068. Overestimation is observed for
4. Verification results for real-life SAR images
First, we will verify our BENCs for the single-look real-life TerraSAR-X images presented in Figures 5 and 7. The obtained data will be considered in subsection 4.1. Besides, in subsection 4.2, we will verify our BENCs for multi-look SAR images of urban area in Canada (Toronto) (these images are presented later). All of them are acquired for HH polarization. As it is stated in file description, approximate number of looks is about 6. Thus, the expected . Similarly, assuming for multi-look data, we can get the limits 0.036…0.054 for blind estimates that can be considered appropriate in practice. Let us keep these limits in mind in further analysis.
4.1. Verification results for single-look SAR images
Let us start from data obtained for
Concerning other blind estimates, they all are within the required limits. For three of four real-life images, 7x7 block size is the best choice from the viewpoint of estimation accuracy.
Let us consider the results for two other BENC methods, both based on scatter-plots. Here we consider only the case of 7x7 blocks according to recommendations given in the previous Section. The estimates obtained by the
Table 5 also contains blind estimates obtained by
One more advantage of
An example of such scatter-plot obtained as (6) for image with clipping effects is given in Fig. 9. Straight line shows the true position of the line to be fitted. As it is seen, there are three clusters (that correspond to large means) positions which are erroneous (vertical coordinates are considerably smaller than they should be). Although line fitting method is robust, the presence of a large percentage of such clusters can lead to essential errors in blind estimation.
4.2. Verification results for multi-look SAR images
The real-life six-look SAR images used in verification tests are given in Fig. 10. From visual inspection, the image in Fig. 10(d) seems to have more complex structures whilst other three images have quite large quasi-homogeneous regions. Let us see how this will influence blind estimates.
The obtained blind estimates for
The presented results clearly show that for estimation techniques based on scatter-plots and robust fitting it is often not enough to carry out robust fitting. Image pre-processing able to partly remove local estimates expected to be abnormal (due to block heterogeneity or to presence of clipping effects) is desirable. Such pre-processing might include image pre-segmentation which in our experiments has been performed by unsupervised variational classification through image multi-thresholding (Klaine et al., 2005). Its advantage is that pre-processing is quite fast. This allows obtaining blind estimates quite quickly since other operations (obtaining of local estimates and robust regression) are also very fast.
5. Conclusions and future work
Some aspects of SAR image simulation have been considered. In particular, it has been stressed that spatial correlation of speckle is to be taken into account. One algorithm to do this is described.
Three methods for blind estimation of noise statistical characteristics in SAR images have been first tested for simulated images. It has been shown that there are several factors influencing their performance. These factors are image content (complexity), the method used and its parameters. It is not always possible to provide blind estimates within desired limits especially for highly textural (complex structure) images. Then, these methods have been verified for real life TerraSAR-X images of limited size of 512x512 pixels. Preliminary tests have clearly demonstrated the presence of essential spatial correlation of speckle, especially for multi-look images. This is taken into account in setting parameters of BENC methods. The block size of 7x7 pixels is recommended for practical use.
The BENC methods based on scatter-plots without image pre-processing produce, on the average, worse accuracy than the method based on mode determination for local estimates’ distribution. If pre-processing is applied, BENC methods (as
This work has been partly supported by French-Ukrainian program Dnipro (PHC DNIPRO 2013, PROJET N° 28370QL).
Abramov S., Lukin V., Ponomarenko N., Egiazarian K., & Pogrebnyak O. (2004). Influence of multiplicative noise variance evaluation accuracy on MM-band SLAR image filtering efficiency. , Vol. 1, pp. 250-252, Kharkov, Ukraine, June 2004 Proceedings of MSMW 2004
Abramov S., Lukin V., Vozel B., Chehdi K., & Astola J. (2008). Segmentation-based method for blind evaluation of noise variance in images, Vol. 2(1), No. 023533, (August 13, 2008), DOI:10.1117/1.2977788 Journal of Applied Remote Sensing,
Abramov S., Zabrodina V., Lukin V., Vozel B., Chehdi K., & Astola J. (2011). Methods for Blind Estimation of the Variance of Mixed Noise and Their Performance Analysis, In: , Ed. J. Awrejcewicz, pp. 49-70, In-Tech, Austria, ISBN 978-953-307-389-7 Numerical Analysis – Theory and Applications
Aiazzi B., Alparone L., Barducci A., Baronti S., Marcoinni P., Pippi I., & Selva M. (2006). Noise modelling and estimation of hyperspectral data from airborne imaging spectrometers. , Vol. 49, No. 1, February 2006 Annals of Geophysics
Anfinsen S.N., Doulgeris A.P., & Eltoft T. (2009). Estimation of the Equivalent Number of Looks in Polarimetric Synthetic Aperture Radar Imagery, , Vol. 47, No. 11, pp. 3795-3809 IEEE Transactions on Geoscience and Remote Sensing
Bekhtin Yu. S. (2011). Adaptive Wavelet Codec for Noisy Image Compression, , Sevastopol, Ukraine, Sept., 2011, pp. 184-188 Proc. of the 9-th East-West Design and Test Symp.
Choi M.G., Jung J.H., & Jeon J.W. (2009). No-reference Image Quality Assessment Using Blur and Noise, , Vol. 50, pp. 163-167 World Academy of Science, Engineering and Technology
Davies E.R. (2000). , World Scientific, ISBN 9810240228 Image Processing for the Food Industry
Di Martino G., Poderico M., Poggi G., Riccio D., & Verdoliva L. (2012). SAR Image Simulation for the Assessment of Despeckling Techniques, , Munich, Germany, July 2012, pp. 1797-1800 Proceedings of IGARSS
Dogan O., & Kartal M. (2010). Time Domain SAR Raw Data Simulation of Distributed Targets, , Article ID 784815 EURASIP Journal on Advances in Signal Processing
DuMouchel W. & O'Brien F. (1989). Integrating a Robust Option into a Multiple Regression Computing Environment in Computing Science and Statistics. , pp. 297-301, American Statistical Association, Alexandria, VA Proc. of the 21st Symposium on the Interface
Egiazarian K., Astola J., Helsingius M., & Kuosmanen P. (1999). Adaptive denoising and lossy compression of images in transform domain. , Vol. 8(3), pp. 233-245, DOI:10.1117/1.482673 Journal of Electronic Imaging
Foi A., Trimeche M., Katkovnik V., & Egiazarian K. (2007). Practical Poissonian-Gaussian Noise Modeling and Fitting for Single Image Raw Data. , Vol. 17, No. 10, pp. 1737-1754 IEEE Transactions on Image Processing
Foi A. (2009). Clipped Noisy Images: Heteroskedastic Modeling and Practical Denoising. , Vol. 89, No. 12, pp. 2609-2629 Signal Processing
Foucher S., Boucher J.-M., & Benie G. B. (2000). Maximum likelihood estimation of the number of looks in SAR images, Proc. of Int. Conf. Microwave, Radar Wireless Communication, Wroclaw, Poland, May 2000, Vol. 2, pp. 657–660
Herrmann J., Faller N., Kern A., & Weber M. (2005). INFOTERRA GMBH Initiatives Commercial Exploitation of TerraSAR-X, , Hannover, Germany, May 2005 Proc. of ISPRS Hannover Workshop
Katsaggelos A.K. (Ed.). (1991). , Springer-Verlag, New York Digital Image Restoration
Klaine L., Vozel B., & Chehdi K. (2005). Unsupervised Variational Classification Through Image Multi-Thresholding. , Antalya, Turkey Proc. of the 13th EUSIPCO Conference
Lee J.-S., Hoppel K., & Mango S.A. (1992). Unsupervised Estimation of Speckle Noise in Radar Images, Vol. 4, pp. 298-305 Int. Journal of Imaging Systems and Technology.
Lee J.-S., Wen J.H., Ainsworth T.I., Chen K.S., & Chen A.J. (2009). Improved sigma filter for speckle filtering of SAR imagery, IEEE Transactions on Geoscience and Remote Sensing Vol. 47(1), pp. 202–213
Lin C.H., Sun Y.N., & Lin C.J. (2010). A Motion Compounding Technique for Speckle Reduction in Ultrasound Images, , Vol. 23(3), pp. 246–257 Journal of Digital Imaging
Liu C., Szeliski R., Kang S.B., Zitnick C.L., & Freeman W.T. (2008). Automatic estimation and removal of noise from a single image. , Vol. 30, No 2, pp. 299-314 IEEE Transactions on Pattern Analysis and Machine Intelligence
Lukin V., Abramov S., Zelensky A., Astola J., Vozel B., & Chehdi K. (2007). Improved minimal inter-quantile distance method for blind estimation of noise variance in images, Proc. SPIE 6748 of Image and Signal Processing for Remote Sensing XIII, 67481I October 24, 2007, DOI:10.1117/12.738006
Lukin V., Ponomarenko N., Egiazarian K., & Astola J. (2008). Adaptive DCT-based filtering of images corrupted by spatially correlated noise, Proc. SPIE 6812 of Image Processing: Algorithms and Systems VI, 68120W, San Jose, USA, January 2008, DOI:10.1117/12.764893
Lukin V., Abramov S., Ponomarenko N., Uss M., Vozel B., Chehdi K., & Astola J. (2009a). Processing of images based on blind evaluation of noise type and characteristics. Proceedings of SPIE Symposium on Remote Sensing, Vol. 7477, Berlin, Germany, September 2009
Lukin V.V., Abramov S.K., Uss M.L., Marusiy I.A., Ponomarenko N.N., Zelensky A.A., Vozel B., & Chehdi K. (2009b). Testing of methods for blind estimation of noise variance on large image database, In: Practical Aspects of Digital Signal Processing, Shahty, Russia, Retrieved from <http://k504.xai.edu.ua/html/prepods/lukin/BookCh1.pdf>
Lukin V., Abramov S., Popov A., Eltsov P., Vozel B., & Chehdi K. (2010). A method for automatic blind estimation of additive noise variance in digital images, Telecomunicaions and Radio Engineering, Vol. 69(19), pp. 1681-1702
Lukin V., Abramov S., Ponomarenko N., Uss M., Zriakhov M., Vozel B., Chehdi K., & Astola J. (2011). Methods and Automatic Procedures for Processing Images Based on Blind Evaluation of Noise Type and Characteristics. SPIE Journal on Advances in Remote Sensing, DOI: 10.1117/1.3539768
Lukin V.V., Abramov S.K., Fevralev D.V., Ponomarenko N.N., Egiazarian K.O., Astola J.T., Vozel B., & Chehdi K. (2011b). Performance evaluation for Blind Methods of Noise Characteristics Estimation for TerraSAR-X Images, Proc. SPIE 8180 of Image and Signal Processing for Remote Sensing XVII, 81800X, Prague, Czech Republic, September 2011, DOI:10.1117/12.897730
Mallat S. (1998). A Wavelet tour of signal processing,Academic Press, San Diego
Oliver C. & Quegan S. (2004). Understanding Synthetic Aperture Radar Images, SciTech Publishing
Plataniotis K.N. & Venetsanopoulos A.N. (2000). Color Image Processing and Applications, Springer-Verlag, NY
Ponomarenko N.N., Lukin V.V., Egiazarian K.O., & Astola J.T. (2010). A method for blind estimation of spatially correlated noise characteristics, Proc. SPIE 7532 of Image Processing: Algorithms and Systems VIII, 753208, San Jose, USA, January 2010, DOI:10.1117/12.847986
Ponomarenko N.N., Lukin V.V., & Egiazarian K.O. (2011). Visually Lossless Compression of Synthetic Aperture Radar Images, Proceedings of ICATT, Kiev, Ukraine, September 2011, pp. 263-265
Ramponi G. & D’Alvise R. (1999). Automatic Estimation of the Noise Variance in SAR Images for Use in Speckle Filtering, Proceedings of EEE-EURASIP Workshop on Nonlinear Signal and Image Processing, Vol. 2, pp. 835-838, Antalya, Turkey
Sendur L. & Selesnick I.W. (2002). Bivariate shrinkage with local variance estimation. IEEE Signal Processing Letters, Vol. 9, No. 12, pp. 438-441
Solbo S. & Eltoft T. (2004). Homomorphic Wavelet-based Statistical Despeckling of SAR Images. IEEE Trans. on Geosc. and Remote Sensing, Vol. GRS-42, No. 4, pp. 711-721
Solbo S. & Eltoft T. (2008). A Stationary Wavelet-Domain Wiener Filter for Correlated Speckle, IEEE Trans. on Geoscience and Remote Sensing, Vol. 46(4), pp. 1219-1230
Touzi R. (2002). A Review of Speckle Filtering in the Context of Estimation Theory. IEEE Transactions on Geoscience and Remote Sensing, Vol. 40, No. 11, pp. 2392-2404
Uss M., Vozel B., Lukin V., & Chehdi K. (2011). Local Signal-Dependent Noise Variance Estimation from Hyperspectral Textural Images. IEEE Journal of Selected Topics in Signal Processing, Vol. 5, No. 2, DOI: 10.1109/JSTSP.2010.2104312
Uss M., Vozel B., Lukin V., & Chehdi K. (2012). Maximum Likelihood Estimation of Spartially Correlated Signal-Dependent Noise in Hyperspectral Images, Optical Engineering, Vol. 51, No 11, DOI: 10.1117/1.OE.51.11.111712
Van Zyl Marais I., Steyn W.H., & du Preez J.A. (2009). On-board image quality assessment for a small low Earth orbit satellite, Proc. of the 7th IAA Symp. on Small Satellites for Earth Observation, Berlin, Germany, May 2009
Vozel B., Abramov S., Chehdi K., Lukin V., Ponomarenko N., Uss M., & Astola J. (2009). Blind methods for noise evaluation in multi-component images, In: Multivariate Image Processing,pp. 263-295, France
Zabrodina V., Abramov S., Lukin V., Astola J., Vozel B., & Chehdi K. (2011). Blind Estimation of mixed noise parameters in images using robust regression curve fitting, Proc. of 19th European Signal Processing Conference EUSIPCO2011, Barcelona, Spain, August 2009, pp. 1135 – 1139, ISSN 2076-1465
Zoran D. & Weiss Y. (2009). Scale Invariance and Noise in Natural Images, Proc. of IEEE 12th International Conference on Computer Vision ICCV, Kyoto, Japan, September 2009, pp. 2209-2216, DOI:10.1109/ICCV.2009.5459476