Statistical cone-beam CT noise reduction with multiscale decomposition and penalized weighted least squares in the projection domain

Abstract

Purposes:

Suppressing noise can effectively promote image quality and save radiation dose in clinical imaging with x-ray computed tomography (CT). To date, numerous statistical noise reduction approaches have ever been proposed in image domain, projection domain or both domains. Especially, a multiscale decomposition strategy can be exploited to enhance the performance of noise suppression while preserving image sharpness. Recognizing the inherent advantage of noise suppression in the projection domain, we have previously proposed a projection domain multiscale penalized weighted least squares (PWLS) method for fan-beam CT imaging, wherein the sampling intervals are explicitly taken into account for the possible variation of sampling rates. In this work, we extend our previous method into cone-beam (CB) CT imaging, which is more relevant to practical imaging applications.

Methods:

The projection domain multiscale PWLS method is derived for CBCT imaging by converting an isotropic diffusion partial differential equation (PDE) in the three-dimensional (3D) image domain into its counterpart in the CB projection domain. With adoption of the Markov random field (MRF) objective function, the CB projection domain multiscale PWLS method suppresses noise at each scale. The performance of the proposed method for statistical noise reduction in CBCT imaging is experimentally evaluated and verified using the projection data acquired by an actual micro-CT scanner.

Results:

The preliminary result shows that the proposed CB projection domain multiscale PWLS method outperforms the CB projection domain single-scale PWLS, the 3D image domain discriminative feature representation (DFR), and the 3D image domain multiscale nonlinear diffusion methods in noise reduction. Moreover, the proposed method can preserve image sharpness effectively while avoiding generation of novel artifacts.

Conclusions:

Since the sampling intervals are explicitly taken into account in the projection domain multiscale decomposition, the proposed method would be beneficial to advanced applications where the CBCT imaging is employed and the sampling rates vary.

Keywords

CT cone-beam noise reduction multiscale decomposition penalized weighted least squares

Get full access to this article

View all access options for this article.

References

Brenner

Hall

. Computed tomography — an increasing source of radiation exposure. N Engl J Med 2007; 357: 2277–2289.

Zhao

Hoffman

, et al. Ultra-low dose lung CT perfusion regularized by a previous scan. Acad Radiol 2009; 16: 363–373.

Tang

Mou

, et al. Noise reduction by projection direction dependent diffusion for low dose x-ray cone beam computed tomography. In: Proceedings of the 10th international meeting on fully three-dimensional image reconstruction in radiology and nuclear medicine, Beijing, China, Sept. 5–10, 2009, pp.381–84.

Maier

Wigstrom

Hofmann

, et al. Three-dimensional anisotropic adaptive filtering of projection data for noise reduction in cone beam CT. Med Phys 2011; 38: 5896–5909.

Elbakri

Fessler

. Statistical image reconstruction for polyenergetic X-ray computed tomography. IEEE Trans Med Imag 2002; 21: 89–99.

Wang

Eremina

, et al. An experimental study on the noise properties of x-ray CT sinogram data in the Radon space. Phys Med Biol 2008; 53: 3327–3341.

Wang

, et al. Nonlinear sinogram smoothing for low-dose x-ray CT. IEEE Trans Nucl Sci 2004; 51: 2505–2513.

Wang

, et al. Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose x-ray computed tomography. IEEE Trans Med Imaging 2006; 25: 1272–1283.

Wang

, et al. Noise reduction for low-dose single-slice helical CT sinograms. IEEE Trans Nucl Sci 2006; 53: 1230–1237.

10.

Wang

Liang

, et al. Dose reduction for kilovoltage cone-beam computed tomography in radiation therapy. Phys Med Biol 2008; 53: 2897–2909.

11.

Wang

Wen

, et al. Multiscale penalized weighted least-squares sinogram restoration for low-dose x-ray computed tomography. IEEE Trans Biomed Eng 2008; 55: 1022–1031.

12.

Zhu

Wang

Xing

. Noise suppression in scatter correction for cone-beam CT. Med Phys 2009; 36: 741–752.

13.

Forthmann

Koehler

Defrise

, et al. Comparing implementations of penalized weighted least-squares sinogram restoration. Med Phys 2010; 37: 5929–5938.

14.

Yang

Silver

Noshi

. Adaptive weighted nonlinear diffusion for Computed Tomography denoising. In: Proceedings of the 11th interntional meeting on fully three-dimensional image reconstruction in radiology and nuclear medicine, Potsdam, Germany, July 11–15, 2011, pp. 210–213.

15.

Prakash

Kalra

Ackman

, et al. Diffuse lung diseases: cT of the chest with adaptive statistical iterative reconstruction technique. Radiology 2010; 256: 261–269.

16.

Tang

Yang

Gong

, et al. Image domain noise reduction with multiscale decomposition and nonlinear diffusion. In: IEEE nuclear science symposium and medical imaging conference, 2015.

17.

Chen

Liu

, et al. Discriminative feature representation: an effective postprocessing solution to low dose CT imaging. Phys Med Biol 2017; 62: 2103–2131.

18.

Chen

Shi

Feng

, et al. Artifact suppressed dictionary learning for low-dose CT image processing. IEEE Trans Med Imaging 2014; 33: 2271–2292.

19.

Chen

Zhang

Yang

, et al. Curve-like structure extraction using minimal path propagation with backtracking. IEEE Trans Image Processing 2016; 25: 988–1003.

20.

Chen

Feng

, et al. Nonlocal prior Bayesian tomographic reconstruction. J Math Imaging Vis 2008; 30: 133–146.

21.

Chen

Yin

, et al. Improving low-dose abdominal CT images by weighted intensity averaging over large-scale neighborhoods. Eur J Radiol 2011; 80: e42–e49.

22.

Chen

Yin

Shi

, et al. Improving abdomen tumor low-dose CT images using a fast dictionary learning based processing. Phys Med Biol 2013; 58: 5803–5820.

23.

Chen

Yang

, et al. Thoracic low-dose CT image processing using an artifact suppressed large-scale nonlocal means. Phys Med Biol 2012; 57: 2667–2688.

24.

Chen

Gao

Nie

, et al. Bayesian Statistical reconstruction for low-dose X-ray computed tomography using an adaptive-weighting nonlocal prior. Comput Med Imaging Graph 2009; 33: 495–500.

25.

Tang

. Statistical CT noise reduction with multiscale decomposition and penalized weighted least squares in the projection domain. Med Phys 2012; 39: 5498–5512.

26.

Liu

Zhang

, et al. Discriminative feature representation to improve projection data inconsistency for low dose CT imaging. IEEE Trans Med Imaging 2017; 36: 2499–2509.

27.

Bai

Yan

Jia

, et al. Z-Index parameterization (ZIP) for volumetric CT image reconstruction via 3D dictionary learning. IEEE Trans Med Imaging 2017; 36: 2466–2478.

28.

Liu

Yang

, et al. 3D Feature constrained reconstruction for low dose CT imaging. IEEE Trans Circuits Syst Video Technol 2018; 28: 1232–1247.

29.

Mallat

. A theory for multiresolution signal decomposition: the wavelet representation. IEEE Trans Pattern Anal Mach Intell 1989; 11: 674–693.

30.

Burt

Adelson

. The Laplacian pyramid as a compact image code. IEEE Trans Commun 1983; 31: 532–540.

31.

Lindeberg

. Scale-Space theory in computer vision. Springer Scale-Space Theory Comp Vis 1997; 1252: 1–165.

32.

Perona

Malik

. Scale-space and edge detection using anisotropic diffusion. IEEE Trans Pattern Anal Mach Intell 1990; 12: 629–639.

33.

Lampert

Wirjadi

. An optimal nonorthogonal separation of the anisotropic Gaussian convolution filter. IEEE Trans Image Processing 2006; 15: 3501–3513.

34.

Abd-Elmoniem

Youssef

Kadah

. Real-time speckle reduction and coherence enhancement in ultrasound imaging via nonlinear anisotropic diffusion. IEEE Trans Biomed Imaging 2002; 49: 997–1014.

35.

Gerig

Kubler

Kikinis

, et al. Nonlinear anisotropic filtering of MRI data. IEEE Trans Med Imaging 1992; 11: 221–231.

36.

Black

Sapiro

Marimont

, et al. Robust anisotropic diffusion. IEEE Trans Image Processing 1998; 7: 421–432.

37.

Yang

Zamyatin

Akino

. Effective data-domain noise and streak reduction for x-ray CT. In: Proceedings of the 11th interntional meeting on fully three-dimensional image reconstruction in radiology and nuclear medicine, Potsdam, Germany, July 11–15, 2011, pp. 290–293.

38.

Hsieh

. Adaptive streak artifact reduction in computed tomography resulting from excessive x-ray photon noise. Med Phys 1998; 25: 2139–2147.

39.

Grangeat

. Mathematical framework of cone-beam 3D reconstruction via the ﬁrst derivative of the Radon transform mathematical methods in tomography. Springer Lect Notes Math 1991; 1497: 66–97.

40.

Defrise

Noo

Kudo

. Improved two-dimensional rebinning of helical cone-beam computerized tomography data using John’s equation. Inverse Prob 2003; 19: S41–S54.

41.

Tang

. Axial cone beam reconstruction by weighted BPF/DBPF and orthogonal butterfly filtering. IEEE Trans Biomed Eng 2016; 63: 1895–1903.

42.

Tang

Hsieh

Dong

, et al. Minimization of over-ranging in helical CT scan via hybrid cone beam image reconstruction – benefits in dose efficiency. Med Phys 2008; 35: 3232–3238.

43.

Kozek

Crounse

Roska

, et al. Multiscale image analysis on the CNN universal machine. In: Fourth IEEE international workshop on cellular neural networks and their applications, Seville, Spain, 1996, pp. 24–26.

44.

Pauwels

Van Gool

Fiddelaers

, et al. An extended class of scale-invariant and recursive scale space filters. IEEE Trans Pattern Anal Mach Intell 1995; 17: 691–701.

45.

Natterer

. The mathematics of computerized tomography. New York: John Wiley & Sons, 1986.

46.

Thirion

. Direct extraction of boundaries from computed tomography scans. IEEE Trans Med Imaging 1994; 13: 322–328.

47.

Patch

. Computation of unmeasured third-generation VCT views from measured views. IEEE Trans Med Imaging 2002; 21: 801–813.

48.

Weickert

ter Haar Romeny

Viergever

. Efﬁcient and reliable schemes for nonlinear diffusion filtering. IEEE Trans Image Process 1998; 7: 398–410.

49.

Yan

Mou

Tang

, et al. Projection correlation based view interpolation for cone beam CT: primary ﬂuence restoration in scatter measurement with a moving beam stop array. Phys Med Biol 2010; 55: 6353–6375.