About this Journal Submit a Manuscript Table of Contents
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 463760, 7 pages
http://dx.doi.org/10.1155/2013/463760
Research Article

An Image Filter Based on Multiobjective Genetic Algorithm and Shearlet Transformation

1Nanjing University of Science and Technology, Nanjing 210000, China
2Nanjing University of Information Science and Technology, Nanjing 210044, China

Received 4 June 2013; Revised 22 July 2013; Accepted 26 August 2013

Academic Editor: Marco Perez-Cisneros

Copyright © 2013 Zhi-yong Fan 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

Rician noise pollutes magnetic resonance imaging (MRI) data, making data’s postprocessing difficult. In order to remove this noise and avoid loss of details as much as possible, we proposed a filter algorithm using both multiobjective genetic algorithm (MOGA) and Shearlet transformation. Firstly, the multiscale wavelet decomposition is applied to the target image. Secondly, the MOGA target function is constructed by evaluation methods, such as signal-to-noise ratio (SNR) and mean square error (MSE). Thirdly, MOGA is used with optimal coefficients of Shearlet wavelet threshold value in a different scale and a different orientation. Finally, the noise-free image could be obtained through inverse wavelet transform. At the end of the paper, experimental results show that this proposed algorithm eliminates Rician noise more effectively and yields better peak signal-to-noise ratio (PSNR) gains compared with other traditional filters.

1. Introduction

Magnetic resonance imaging (MRI) devices are important imaging equipment, and the image quality directly impacts the diagnosis accuracy. However, MRI images are frequently contaminated by Rician noise during image gaining or transmission [1]. This phenomenon makes noise reduction to be one of the most important problems in image processing. Preservation of image details and attenuation of noise are both critical, but they are contradictory in nature. Therefore, this research is focused on Rician noise elimination and data details preservation at the same time.

Because of its good performance in both time domain and frequency domain, wavelet transform has become one of the most active research fields in image processing. It provides better results and preserves more details compared with traditional algorithms. However, wavelet transform cannot achieve optimal sparse for images containing higher-dimension singularity. To overcome the limitation, multiscale geometric analysis theory is proposed, and, based on it, a series of methods sprang out, for example, ridgelet [2], curvelet [3], contourlet [4], and bandlet [4]. One of the most successful ideas is the curvelets of Candes and Donoho, which achieve an (almost) optimal approximation for 2D piecewise smooth functions with discontinuities along with curves.

Recently, Labate et al. described a new class of multidimensional representation systems, which is called Shearlet. One advantage of this approach is that these systems can be constructed using generalized multiresolution analysis and implemented efficiently using a classical cascade algorithm [511].

Simple threshold denoising method of classical Shearlet transform could yield good performance because of the method’s multiscale and multidirection characteristics and image sparse representation. However, there is still room for improvement because classical Shearlet algorithm does not take energy distribution of different scales and different directions into consideration; as a result, it kills the coefficient excessively; therefore, image details could be lost. In order to solve the problem, Sun and Zhao [12] proposed a particle swarm optimization; it uses adaptive algorithm to search for optimal threshold of the highest PSNR values.

Based on these previous achievements, this paper proposed a new image-filtering algorithm. It has three characteristics: it uses soft threshold in Shearlet, it builds target function in MOGA by several evaluation methods, and it uses the MOGA to optimize coefficients of Shearlet wavelet threshold value in different scale and a different orientation.

The rest of this paper is organized as follows. Section 2 introduces related theories. Section 3 explains our algorithm, including workflow, Section 4 presents the experiment results of proposed algorithm, and Section 5 concludes this paper.

2. Related Theories

2.1. Rician Noise

Noised MRI image can be defined as ; here, represent original image pixels, and represent is noised pixels. When MR images are computed by using the magnitude of single-complex raw data, its distribution can be modeled as a Rician model [1315]. Consider the following: Here, is the standard deviation (STD) of Gaussian noise, is the amplitude of the signal without noise, is the value of the magnitude image, and is the 0th-order modified Bessel function. This model is used by the majority of the noise estimation methods.

When SNR is small enough (i.e., ), the Rician distribution is considered as a Rayleigh distribution. Consider the following:

When SNR is high (i.e.,), the Rician distribution is approximated as a Gaussian distribution.

2.2. Shearlet Transform

Labate et al. [511] proposed Shearlet transform based on wavelet. With dimension , consider the following affine system: Here, , and , are 2*2 invertible matrices with .

If satisfied Parseval , then, those elements of are called composite wavelets.

Shearlet is a special example of only when Here, is the anisotropic dilation matrix, and is the shear matrix.

For ,, when ,, and satisfy

Then, we get

Then, form a tiling of the set

From the condition on the support of and, it is easily deduced that have frequency support contained in the set as follows:

Thus, every element in is supported on a pair of trapezoids of approximate size , oriented along lines of slope.

For , here, is the vertical cone when the following formula was satisfied:

Then, collection is a Parseval frame for.

2.3. Multiobjective Genetic Algorithm (MOGA)

Multiobjective genetic algorithm seeks feasible solutions to problems comprising multiple objectives which are often in conflict with each other. A general minimization problem ofobjectives can be mathematically stated as , whereis the dimension of the decision variable space. Consider the following.

Minimize which satisfies where is theth objective function, is theth inequality constraint, and is theth equality constraint. The multiobjective optimization problem then reduces to finding an, such thatis optimized.

3. Proposed Algorithm

3.1. Threshold Rule

Threshold rule is the most important problem in image denoising of transform domain, and the hard-threshold and the soft-threshold approaches are two options. Donoho and Johnstone [16] proposed the following threshold rule: Here,is the pixels number of image, and is the noise level.

Research shows that Donoho threshold is the optimal threshold limit not the optimal threshold. With this considered, Donoho and Johnstone [16] proposed an improved threshold rule as follows:

As many researchers point out [12, 13], (12) did not consider energies of subwavelets in a different direction while being in the same scale, and this imperfection will make coefficients too much stifled.

Considering the variability of image content and Shearlet transformation of multiscale and multidirection characteristics, a novel threshold selection rule is proposed based on Shearlet transform multiscale and multidirection; this rule is the following.

Comprehensively considering complexity of image and the multiscale and multidirection characteristics of Shearlet transform, this paper proposed the following adaptive threshold rule:

Here, Sigmoid is adopted to build our rules. The Sigmoid curve is a mathematical concept which has been widely used to model the natural life cycle of many things, for its derivative is continuous and with higher accuracy.is the scale level, andis theth direction under theth scale level.

3.2. Target Function

We build MOGA target function by the signal-to-noise ratio (SNR) and the mean square error (MSE).

Signal-to-noise ratio (SNR) can be defined as

Here,is original image with size of pixels,is filtered image of noised image, andare coordinates of pixels.

Mean square error (MSE) expressed the correlation between images, and it is defined as follows:

Here,is original image with size of pixels,is filtered image of noised image, andare coordinates of pixels.

Our target function is defined as follows:

Here, , are weight coefficients of SNR and MSE.

3.3. Proposed Model

The most critical problem which lies in our optimal filtering performance study is, under optimization criterion, how to decide coefficients, considering energy of subwavelets not only in different scale but also in different direction.

Here, we proposed our algorithm which adopts MOGA algorithm to decide coefficients, of each subwavelet in different scale and direction of Shearlet transform, intending to get optimal filtering performance.

Our algorithm works as follows [17, 18]; see Figure 1.

463760.fig.001
Figure 1: Workflow of proposed algorithm.

Step 1 (initialization). Generate an initial population containing strings, where is the number of strings in each population. These strings contain weight coefficients of SNR, MSE, weight coefficients of Shearlet subwavelets,offunction, and other parameters in MOGA; thus, we need the following.

Step 2 (evaluation). (1)Use Shearlet transform to decompose target image.(2)Multiply subwavelets by weight coefficients .(3)Filter subwavelets by threshold rule.(4)Reconstruct image by filtered subwavelets.(5)Calculate the values of the objective functions (16) for the generated strings.(6)Update a tentative setoff Pareto optimal solution.

Step 3 (selection). Calculate the fitness value of each string using the random weights in (3). Select a pair of strings from the current population according to the following selection probability.

Step 4 (crossover). For each selected pair, apply a crossover operation to generate two new strings. new strings are generated by the crossover.

Step 5 (mutation). For each bit value of the strings generated by the crossover, apply a mutation with a prespecified mutation probability.

Step 6 (elitist strategy). Randomly remove strings from the set of strings generated by previous operations, and replace them with strings randomly selected from tentative set of Pareto optimal solutions.

Step 7 (termination test). If one stopping condition in the following is satisfied, go to Step 8; if not, return to Step 2.(i)Maximum iterations are exceeded.(ii)The optimal target value is achieved.

Step 8 (algorithm termination). Exit optimal algorithm.

4. Experimental Results and Analysis

4.1. Evaluation Index

Peak signal-to-noise ratio (PSNR) is defined as

Here,is original image with size of pixels,is filtered image of noised image, andare coordinates of pixels.

4.2. Experimental Results

To verify the validity of the algorithm, this paper designed two kinds of experimental methods to verify its effectiveness. One is use of objective data such as PNSR and MSE to objectively analyze its performance; and the other is making us able to obverse filtering performance directly by naked eyes [1921].

Experiment 1. We did filtering experiments on standard images Lena and Barbara in different noise level and listed results in Table 1. As we have seen from Table 1, PSNR of proposed algorithm (Shearlet-MOGA) is higher than PSNR of classical Shearlet algorithm, and its performance will be better with noise level increased.
Figure 2 is the original MR image we adopted to do experiments. Adding different noise level to Figure 2, we did filtering work by classical Shearlet and proposed algorithm in this paper and showed the statics data of MSE and PSNR as Tables 2 and 3.
In Table 2, the excellent effect of our algorithm is dramatic, shown in and our proposed MSE is smaller than classical Shearlet algorithm. Similar good results were found when the same experiment was repeated on PSNR. In Table 3, the PSNR of proposed algorithm is greater than that of classical Shearlet algorithm.

tab1
Table 1: Filtering results for Lena and Barbara.
tab2
Table 2: MSE in different (%) and different algorithm to that in Figure 2.
tab3
Table 3: PSNR in different (%) and different algorithm to that in Figure 2.
463760.fig.002
Figure 2: Original MR image.

Experiment 2. To evaluate the performance of proposed algorithm by naked eyes directly, several classical images such as Lena, Baboon, Barbara, and MRI are adopted to do filtering work, and all relative images are shown in Figure 3.
Figure 3(a) is the original Lena. Adding 5% Rician noise level to Lena, we get Figure 3(b).
Filtering Figure 3(b) by classical Shearlet algorithm, we got Figure 3(c). Figure 3(d) is the output of the filtering work we did to Figure 3(b) by proposed algorithm.
We did similar experiment to the image of Baboon. Add 10% Rician noise level to Baboon, we get Figure 3(e). Filtering Figure 3(f) by classical Shearlet algorithm, we got Figure 3(g). Figure 3(h) is the output of the filtering work we did to Figure 3(f) by proposed algorithm.
The image of Barbara is also adopted by us to test our algorithm. Figure 3(i) is the original Barbara. Figure 3(j) is Barbara noised by 15% Rician noise level. Filtering Figure 3(j) by classical Shearlet algorithm, we got Figure 3(k). Figure 3(l) is the output of the filtering work we did to Figure 3(j) by proposed algorithm.
At last, we measured our algorithm performance on MRI image. Figure 3(m) is the original MRI. Figure 3(n) is the MRI noised by 20% Rician noise level. Filtering Figure 3(n) by classical Shearlet algorithm, we got Figure 3(o). Figure 3(p) is the output of the filtering work we did to Figure 3(n) by proposed algorithm.
Through simple comparison, we can see that our proposed algorithm could effectively remove the noise from the degraded image of Rician noise with unknown intensity level and protect the image details better at the same time. To MRI image, experiments Paying particular attention data show that our algorithm has excellent performance in background. After strict analysis, we concluded that our algorithm retained the consistent component of low frequency in frequency domain by low-pass filtering, and background of MRI has this nature.

fig3
Figure 3: Experiment images in different noise levels and different algorithms.

5. Conclusions

In order to eliminate Rician noise and preserve image details as much as possible, this paper proposed a new image-filtering algorithm based on MOGA and classical Shearlet transform. It builds target functions in MOGA by several evaluation methods such as SNR and MSE. It also uses MOGA to find optimal Shearlet wavelet threshold coefficients in a different scale and different orientation. Computer simulations results are given to verify the effectiveness of this algorithm. At last, experiments data show that our algorithm has excellent performance in MRI imaging.

Acknowledgments

Whatever implementation the authors come up with needs to be based on classical Shearlet algorithm, and they cannot execute their work without Professor Labate’s earlier work and his open-source code online. The authors would like to thank Professor Labate. The authors would also like to express their thanks to the anonymous reviewers whose comments greatly improved the paper. This paper was funded under a Grant from the National Natural Science Foundation of China (no. 60773172) and a Grant from the Jiangsu Province Natural Science Foundation (BK2008411).

References

  1. H. Gudbjartsson and S. Patz, “The Rician distribution of noisy MRI data,” Magnetic Resonance in Medicine, vol. 34, no. 6, pp. 910–914, 1995. View at Publisher · View at Google Scholar · View at Scopus
  2. E. J. Candes, Monoscaleridgelets for the Representation of Images with Edges, Stanford University, 1999.
  3. E. J. Candes and D. L. Donoho, “Curvelets: a surprisingly effective non-adaptive representation for objects with edges,” in Curves and Surfaces Fitting, Saint-Malo 1999, pp. 105–120, Nashville, Tenn, USA, 2000.
  4. M. N. Do and M. Vetterli, “Contourlets,” in Beyond Wavelets, G. V. Welland, Ed., Academic Press, 2003.
  5. D. Labate, W. Q. Lim, G. Kutyniok, and G. Weiss, “Sparse multidimensional representation using shearlets,” in Wavelets XI, vol. 5914 of Proceedings of SPIE, pp. 59140U–59140U-9, International Society for Optics and Photonics, August 2005. View at Publisher · View at Google Scholar · View at Scopus
  6. S. Yi, D. Labate, G. R. Easley, and H. Krim, “A shearlet approach to edge analysis and detection,” IEEE Transactions on Image Processing, vol. 18, no. 5, pp. 929–941, 2009. View at Publisher · View at Google Scholar · View at Scopus
  7. S. Yi, D. Labate, G. R. Easley, and H. Krim, “Edge detection and processing using shearlets,” in Proceedings of 15th IEEE International Conference on Image Processing (ICIP '08), pp. 1148–1151, IEEE, October 2008. View at Publisher · View at Google Scholar · View at Scopus
  8. G. R. Easley, D. Labate, and F. Colonna, “Shearlet-based total variation diffusion for denoising,” IEEE Transactions on Image Processing, vol. 18, no. 2, pp. 260–268, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  9. S. Yi, D. Labate, G. R. Easley, and H. Krim, “A shearlet approach to edge analysis and detection,” IEEE Transactions on Image Processing, vol. 18, no. 5, pp. 929–941, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  10. G. R. Easley, D. Labate, and W.-Q. Lim, “Optimally sparse image representations using shearlets,” in Proceedings of the 40th Asilomar Conference on Signals, Systems, and Computers (ACSSC '06), pp. 974–978, October-November 2006. View at Publisher · View at Google Scholar · View at Scopus
  11. X. Chen, C. Deng, and S. Wang, “Shearlet-based adaptive shrinkage threshold for image denoising,” in Proceedings of the 1st International Conference on E-Business and E-Government (ICEE '10), pp. 1616–1619, May 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. H. Sun and J. Zhao, “Shearlet threshold denoising method based on two sub-swarm exchange particle swarm optimization,” in Proceedings of IEEE International Conference on Granular Computing (GrC '10), pp. 449–452, IEEE, August 2010. View at Publisher · View at Google Scholar · View at Scopus
  13. A. Firouzmanesh and P. Boulanger, “Image De-blurring using shearlets,” in Proceedings of the 9th Conference on Computer and Robot Vision (CRV '12), pp. 167–173, May 2012.
  14. Q. Guo, S. Yu, X. Chen, C. Liu, and W. Wei, “Shearlet-based image denoising using bivariate shrinkage with intra-band and opposite orientation dependencies,” in Proceedings of the International Joint Conference on Computational Sciences and Optimization (CSO '09), vol. 1, pp. 863–866, IEEE, April 2009. View at Publisher · View at Google Scholar · View at Scopus
  15. T. Sun and Y. Neuvo, “Detail-preserving median based filters in image processing,” Pattern Recognition Letters, vol. 15, no. 4, pp. 341–347, 1994. View at Publisher · View at Google Scholar · View at Scopus
  16. D. L. Donoho and I. M. Johnstone, “Adapting to unknownsmoothness via wavelet shrinkage,” Journal of the American Statistical Association, vol. 12, no. 90, pp. 1200–1224., 1995.
  17. T. Murata and H. Ishibuchi, “MOGA: multi-objective genetic algorithms,” in Proceedings of IEEE International Conference on Evolutionary Computation, vol. 1, pp. 289–294, IEEE, November-December 1995. View at Scopus
  18. K. Deb, “Multi-objective genetic algorithms: problem difficulties and construction of test problems,” Evolutionary Computation, vol. 7, no. 3, pp. 205–230, 1999. View at Scopus
  19. T. Li and X. Ye, “Improved stability criteria of neural networks with time-varying delays: an augmented LKF approach,” Neurocomputing, vol. 73, no. 4–6, pp. 1038–1047, 2010. View at Publisher · View at Google Scholar · View at Scopus
  20. J. Xu, Y.-Y. Cao, D. Pi, and Y. Sun, “An estimation of the domain of attraction for recurrent neural networks with time-varying delays,” Neurocomputing, vol. 71, no. 7–9, pp. 1566–1577, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. M. Xia, J. Fang, Y. Tang, and Z. Wang, “Dynamic depression control of chaotic neural networks for associative memory,” Neurocomputing, vol. 73, no. 4–6, pp. 776–783, 2010. View at Publisher · View at Google Scholar · View at Scopus