The performance comparison of the LWF-BUWPD and several state-of-the-art image denoising algorithms
Two-dimensional over-complete wavelet packet transform can better represent the texture and long oscillatory patterns in natural images.
- Adaptive wavelet packet transform
- Wiener cost function
- image denoising
Wavelet with vanishing moments are very effective for representing piecewise smooth images . However, two-dimensional separable wavelets are ill-suited to represent long oscillatory patterns in images with abundant textures, partly owing to their poor directional selectivity in frequency domain. These oscillatory variations of intensity can only be represented by small-scale wavelet coefficients. In some image-processing applications, such as image denoising or image compression, those small-scale coefficients are quantized to zero in the low bit rate image compression and are thresholded or shrunken to zero in image denoising, which degrades compression and denoising performance significantly. To overcome this circumstance, one way is to find the more suitable image representation techniques such as curvelet, contourlet , and ridgelet . But these methods need the researchers to design the new directional compact filter or multichannel filter banks, which is also challenging in filter design area. Another way is to improve the idea of wavelet design method to accommodate the new requirement where the over-complete wavelet packet decomposition is proposed.
Over-complete wavelet packet contains a mass of libraries of waveforms, from which the best wavelet packet base can be selected to efficiently represent long oscillatory patterns given the corresponding criterion. For example, in image compression application, the best wavelet packet base is usually found by pruning the full wavelet packet decomposition given a user predefined cost function. Recently, a large variety of cost function have been proposed, such as Shannon’s entropy cost function and vector entropy cost function, which are used in the rate distortion control strategy. On the other hand, different from finding the best tree structure of the wavelet packet decomposition, some researchers devoted themselves for finding the best wavelet packet decomposition base  which is equivalent to find the best filter banks of wavelet packet decomposition.
In image denoising application, due to unknown noiseless image as well as a great diversity of filtering methods in the transform domain, selecting the best wavelet packet base is always a difficult problem. Enlightened by the idea of doubly local Wiener filtering method  and spatially adaptive wavelet domain shrinkage image denoising algorithms, the best wavelet packet bases selection method using the local Wiener cost function and its application in image denoising is proposed. In this chapter, we will first give the detail of the best wavelet packet-based selection algorithm and then discuss its application in image denoising.
2. Best wavelet packet base selection
Let be the wavelet function and scaling functions, the basic idea of wavelet packet function can be defined as:
where are the orthonormal filters, respectively. Given a signal which is satisfied the Nyquist sampling criterion, subspace is usually assumed to be the first-level subspace of signal . The subspace with depth is defined as . To understand this subspace easily, the dynamic interval can be associated with the subspace . For example, the interval is equivalent to .
2.1. Two-dimensional wavelet packet bases
The two-dimensional separable wavelet packet functions are defined as the tensor products of two one-dimensional wavelet packet functions, that are,
Similar to the one-dimensional case, the subspace
is referred to as the subspace with depth and a wavelet packet function . Figure 1 gives the subspace of 2D wavelet packet base diagram. It can be seen that each node is divided into two branches. Let us associate the dyadic square with the subspace , for example, associates with the square , associates with the rectangle , and associates with the square . It has been proved that: if the squares , are a partition of the unit square , then the family of functions constitutes an orthogonal wavelet packet base of .
2.2. Best wavelet packet base selection under wiener cost function
Similar to the tree structure illustrated in Fig. 1, 2D wavelet packet decomposition can be easily expressed by the quad-tree structure with the root node . Each node (except the last level node) has four child nodes and . This quad-tree structure facilitates the best wavelet packet decomposition procedure. Given the image and the noise level (if unknown, the noise level can be estimated by the MAD estimator in the wavelet domain), the wavelet packet coefficient at the node under J-level wavelet packet decomposition is represented as . For each node, its Wiener cost function is computed by
Assuming all the node in the quad-tree decomposition forms a set , our task is to find a subset from set Ω which have the smallest total cost function (the summation of the Wiener cost function at all wavelet packet decomposition node).
Similar to most of the search algorithms of the best wavelet packet base, the best wavelet packet base under Wiener cost function is obtained by pruning a J-level full quad-tree from bottom to top. The searching algorithm can be described as follows:
Let represent 22J leaf nodes in the J-level full quad-tree.
For and each node in the j-th level, calculate the corresponding Wiener cost function
as well as the summation of the Wiener cost functions of its four child nodes by
If , then
Calculate the Wiener cost function of each set by
If , set and return the step (ii); otherwise, output .
The set is composed of all leaf nodes of the best quad-tree decomposition structure. In this way, given the input noisy image and the noise level, the optimal orthogonal wavelet packet decomposition structure and the minimal total Wiener cost function can be obtained. Figure 2 shows a demo of best wavelet packet decomposition where Fig. 2(a) is its square representation and Fig. 2 (b) is the corresponding tree structure of this best wavelet packet decomposition.
3. Image denoising algorithms based on best wavelet packet decomposition
This section gives two image denoising algorithms based on the best wavelet packet decomposition. In the above section, the optimal wavelet packet decomposition structure is obtained under the Wiener cost function. When computing the Wiener function, the noiseless image and the noise level are assumed as known. In practice, this assumption is unreasonable. It is well known that the noise level can be accurately estimated by the MAD estimator from the input noisy image. The key problem is how to estimate the noiseless image. Enlightened by the empirical Wiener filtering and the doubly local Wiener filtering image denoising algorithms, we first filter the noisy image to get the pilot image, and then the pilot image is used to get the best wavelet packet decomposition structure.
It is well known that undecimated wavelet packet decomposition can achieve better image denoising performance than the decimated wavelet packet decomposition. This is owing to the property that they are shift invariance and robustness. So, the undecimated wavelet packet decomposition can ameliorate some unpleasant phenomenon that appears in the maximally decimated wavelet packet decomposition, such as Gibbs-like ringing around edges and specks in smooth region.
In what follows, we will give the detail of the image denoising algorithm using the undecimated wavelet decomposition. Let be a input noisy image of size , the operator denotes circularly shifting the input image by indices in the vertical direction and indices in the horizontal direction, and the operator is a similar operation but in the opposite direction. In the proposed algorithm, we first use the decimated best wavelet packet decomposition algorithm to denoise all possible shift versions of the noisy image to get a set of denoised images and then unshifting these denoised images and average them to get the final denoised image. The new algorithm is referred to as the local Wiener filtering using the best undecimated wavelet packet decomposition (LWF-BUWPD), which can be summarized as follows:
Using the local Wiener filtering image denoising algorithm to each shift of input noisy image shift () to get a pilot image ;
Given the pilot image and noise level, the best wavelet packet decomposition structure can be obtained by using the searching algorithm in section 2.
Given the best wavelet packet decomposition structure, the empirical energy distribution of the pilot image can be estimated by
where and represent the directional window and the pilot image’s best wavelet packet decomposition coefficients, respectively.
Using the estimated energy distributions and noise level, the local Wiener filtering is operated on the base wavelet packet decomposition coefficients of the noisy image, that is,
where are the best wavelet packet decomposition coefficients of the noisy image in the subspace .
Unshift all the shifted denoised images and average them to obtain the final denoised image
4. Experimental results
We choose the 8-bit 512 × 512 grayscale images “Lena” and “Barbara,” and a 256 × 256 texture image, as the test images. In the proposed image denoising algorithm, the best wavelet packet decomposition structure is varied with different shift noisy images and different noise level. For better illustration, the best wavelet packet decomposition structure for different noise level is shown in Fig. 3 for “Barbara” image.
In Table 1 and Fig. 4, we give the denoising performance of the image denoising algorithms using the undecimated best wavelet packet decomposition. The experimental results show that for images of structural textures, for example, “Barbara” and texture images, the proposed algorithm greatly improves denoising performance as compared with the existing state-of-the-art algorithms.