Denoising of PET images by combining wavelets and curvelets for improved preservation of resolution and quantitation

Denoising of Positron Emission Tomography (PET) images is a challenging task due to the inherent low signal-to-noise ratio (SNR) of the acquired data. A pre-processing denoising step may facilitate and improve the results of further steps such as segmentation, quantification or textural features characterization. Different recent denoising techniques have been introduced and most state-of-the-art methods are based on filtering in the wavelet domain. However, the wavelet transform suffers from some limitations due to its non-optimal processing of edge discontinuities. More recently, a new multi scale geometric approach has been proposed, namely the curvelet transform. It extends the wavelet transform to account for directional properties in the image. In order to address the issue of resolution loss associated with standard denoising, we considered a strategy combining the complementary wavelet and curvelet transforms. We compared different figures of merit (e.g. SNR increase, noise decrease in homogeneous regions, resolution loss, and intensity bias) on simulated and clinical datasets with the proposed combined approach and the wavelet-only and curvelet-only filtering techniques. The three methods led to an increase of the SNR. Regarding the quantitative accuracy however, the wavelet and curvelet only denoising approaches led to larger biases in the intensity and the contrast than the proposed combined algorithm. This approach could become an alternative solution to filters currently used after image reconstruction in clinical systems such as the Gaussian filter.

[1]  C. L. Le Rest,et al.  Validation of a Monte Carlo simulation of the Philips Allegro/GEMINI PET systems using GATE , 2006, Physics in medicine and biology.

[2]  I. Selesnick,et al.  Bivariate shrinkage with local variance estimation , 2002, IEEE Signal Processing Letters.

[3]  Michael B. Wakin Sparse Image and Signal Processing: Wavelets, Curvelets, Morphological Diversity (Starck, J.-L., et al; 2010) [Book Reviews] , 2011, IEEE Signal Processing Magazine.

[4]  I. Johnstone,et al.  Adapting to Unknown Smoothness via Wavelet Shrinkage , 1995 .

[5]  Kris Thielemans,et al.  Study of direct and indirect parametric estimation methods of linear models in dynamic positron emission tomography. , 2008, Medical physics.

[6]  Christian Roux,et al.  A Fuzzy Locally Adaptive Bayesian Segmentation Approach for Volume Determination in PET , 2009, IEEE Transactions on Medical Imaging.

[7]  Martin A Lodge,et al.  Simultaneous measurement of noise and spatial resolution in PET phantom images. , 2010, Physics in medicine and biology.

[8]  Andrew J. Reader,et al.  One-pass list-mode EM algorithm for high-resolution 3-D PET image reconstruction into large arrays , 2002 .

[9]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[10]  R Lu,et al.  Noisy Image Segmentation by Modified Snake Model , 2006 .

[11]  P. Dutilleux An Implementation of the “algorithme à trous” to Compute the Wavelet Transform , 1989 .

[12]  A. Webb,et al.  Introduction to biomedical imaging , 2002 .

[13]  G. A Theory for Multiresolution Signal Decomposition : The Wavelet Representation , 2004 .

[14]  Jyh-Cheng Chen,et al.  Development of wavelet de-noising technique for PET images. , 2005, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[15]  Jeff Kershaw,et al.  PET kinetic analysis: wavelet denoising of dynamic PET data with application to parametric imaging , 2007, Annals of nuclear medicine.

[16]  P B Hoffer,et al.  Computerized three-dimensional segmented human anatomy. , 1994, Medical physics.

[17]  Laurent Demanet,et al.  Fast Discrete Curvelet Transforms , 2006, Multiscale Model. Simul..

[18]  Xiaoming Huo,et al.  Sparse image representation via combined transforms , 1999 .

[19]  Robert D. Nowak,et al.  Wavelet-based statistical signal processing using hidden Markov models , 1998, IEEE Trans. Signal Process..

[20]  DAVID K. MASLEN,et al.  COMPUTATIONAL HARMONIC ANALYSIS FOR TENSOR FIELDS , .

[21]  Mohamed-Jalal Fadili,et al.  The Undecimated Wavelet Decomposition and its Reconstruction , 2007, IEEE Transactions on Image Processing.

[22]  Paul E. Kinahan,et al.  PET Image Reconstruction , 2005 .

[23]  F. Turkheimer,et al.  Evaluation of a 3D local multiresolution algorithm for the correction of partial volume effects in positron emission tomography. , 2011, Medical physics.

[24]  John Dickson,et al.  PET imaging for prediction of response to therapy and outcome in oesophageal carcinoma , 2011, European Journal of Nuclear Medicine and Molecular Imaging.

[25]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[26]  E. Candès New tight frames of curvelets and optimal representations of objects with C² singularities , 2002 .

[27]  Anthonin Reilhac,et al.  Optimization of dynamic measurement of receptor kinetics by wavelet denoising , 2006, NeuroImage.

[28]  Emmanuel J. Candès,et al.  The curvelet transform for image denoising , 2001, Proceedings 2001 International Conference on Image Processing (Cat. No.01CH37205).

[29]  H R Schelbert,et al.  Blood flow and metabolism by PET. , 1994, Cardiology clinics.

[30]  Martin Vetterli,et al.  Adaptive wavelet thresholding for image denoising and compression , 2000, IEEE Trans. Image Process..

[31]  Jean-Luc Starck,et al.  Very high quality image restoration by combining wavelets and curvelets , 2001, SPIE Optics + Photonics.

[32]  D. Donoho,et al.  Atomic Decomposition by Basis Pursuit , 2001 .

[33]  Carole Lartizien,et al.  Incorporating Patient-Specific Variability in the Simulation of Realistic Whole-Body $^{18}{\hbox{F-FDG}}$ Distributions for Oncology Applications , 2009, Proceedings of the IEEE.

[34]  S Ali Hyder,et al.  An efficient algorithm for denoising MR and CT images using digital curvelet transform. , 2011, Advances in experimental medicine and biology.

[35]  Dimitris Visvikis,et al.  PET Image Denoising Using a Synergistic Multiresolution Analysis of Structural (MRI/CT) and Functional Datasets , 2008, Journal of Nuclear Medicine.

[36]  Stéphane Mallat,et al.  A Wavelet Tour of Signal Processing - The Sparse Way, 3rd Edition , 2008 .

[37]  F E Turkheimer,et al.  Modeling Dynamic PET-SPECT Studies in the Wavelet Domain , 2000, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[38]  Eric Andres,et al.  3-D Discrete Analytical Ridgelet Transform , 2006, IEEE Transactions on Image Processing.

[39]  Frank Wübbeling PET Image Reconstruction , 2012 .

[40]  E. Candès,et al.  New tight frames of curvelets and optimal representations of objects with piecewise C2 singularities , 2004 .

[41]  F E Turkheimer,et al.  Multiresolution Analysis of Emission Tomography Images in the Wavelet Domain , 1999, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[42]  Hamid Soltanian-Zadeh,et al.  Assessment of the Wavelet Transform in Reduction of Noise from Simulated PET Images , 2009, Journal of Nuclear Medicine Technology.

[43]  G. Easley,et al.  Sparse directional image representations using the discrete shearlet transform , 2008 .

[44]  M. Hatt,et al.  Reproducibility of Tumor Uptake Heterogeneity Characterization Through Textural Feature Analysis in 18F-FDG PET , 2012, The Journal of Nuclear Medicine.

[45]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[46]  Eero P. Simoncelli Modeling the joint statistics of images in the wavelet domain , 1999, Optics & Photonics.

[47]  D. Visvikis,et al.  Incorporation of wavelet-based denoising in iterative deconvolution for partial volume correction in whole-body PET imaging , 2009, European Journal of Nuclear Medicine and Molecular Imaging.