Three-dimensional Anisotropic Adaptive Filtering of Projection Data for Noise Reduction in Cone Beam CT

Overview

Journal Med Phys

Specialty Biophysics

Date 2011 Nov 4

PMID 22047354

Citations 9

Authors

Andreas Maier

Lars Wigstrom

Hannes G Hofmann

Joachim Hornegger

Lei Zhu

Norbert Strobel

Rebecca Fahrig

Affiliations

Soon will be listed here.

Abstract

Purpose: The combination of quickly rotating C-arm gantry with digital flat panel has enabled the acquisition of three-dimensional data (3D) in the interventional suite. However, image quality is still somewhat limited since the hardware has not been optimized for CT imaging. Adaptive anisotropic filtering has the ability to improve image quality by reducing the noise level and therewith the radiation dose without introducing noticeable blurring. By applying the filtering prior to 3D reconstruction, noise-induced streak artifacts are reduced as compared to processing in the image domain.

Methods: 3D anisotropic adaptive filtering was used to process an ensemble of 2D x-ray views acquired along a circular trajectory around an object. After arranging the input data into a 3D space (2D projections + angle), the orientation of structures was estimated using a set of differently oriented filters. The resulting tensor representation of local orientation was utilized to control the anisotropic filtering. Low-pass filtering is applied only along structures to maintain high spatial frequency components perpendicular to these. The evaluation of the proposed algorithm includes numerical simulations, phantom experiments, and in-vivo data which were acquired using an AXIOM Artis dTA C-arm system (Siemens AG, Healthcare Sector, Forchheim, Germany). Spatial resolution and noise levels were compared with and without adaptive filtering. A human observer study was carried out to evaluate low-contrast detectability.

Results: The adaptive anisotropic filtering algorithm was found to significantly improve low-contrast detectability by reducing the noise level by half (reduction of the standard deviation in certain areas from 74 to 30 HU). Virtually no degradation of high contrast spatial resolution was observed in the modulation transfer function (MTF) analysis. Although the algorithm is computationally intensive, hardware acceleration using Nvidia's CUDA Interface provided an 8.9-fold speed-up of the processing (from 1336 to 150 s).

Conclusions: Adaptive anisotropic filtering has the potential to substantially improve image quality and/or reduce the radiation dose required for obtaining 3D image data using cone beam CT.

Citing Articles

Trainable joint bilateral filters for enhanced prediction stability in low-dose CT.

Wagner F, Thies M, Denzinger F, Gu M, Patwari M, Ploner S Sci Rep. 2022; 12(1):17540.

PMID: 36266416 PMC: 9585057. DOI: 10.1038/s41598-022-22530-4.

Temporal and volumetric denoising via quantile sparse image prior.

Schirrmacher F, Kohler T, Endres J, Lindenberger T, Husvogt L, Fujimoto J Med Image Anal. 2018; 48:131-146.

PMID: 29913433 PMC: 6088387. DOI: 10.1016/j.media.2018.06.002.

Low-dose x-ray tomography through a deep convolutional neural network.

Yang X, DE Andrade V, Scullin W, Dyer E, Kasthuri N, De Carlo F Sci Rep. 2018; 8(1):2575.

PMID: 29416047 PMC: 5803233. DOI: 10.1038/s41598-018-19426-7.

Applications of nonlocal means algorithm in low-dose X-ray CT image processing and reconstruction: A review.

Zhang H, Zeng D, Zhang H, Wang J, Liang Z, Ma J Med Phys. 2017; 44(3):1168-1185.

PMID: 28303644 PMC: 5381744. DOI: 10.1002/mp.12097.

Dynamic detector offsets for field of view extension in C-arm computed tomography with application to weight-bearing imaging.

Herbst M, Schebesch F, Berger M, Choi J, Fahrig R, Hornegger J Med Phys. 2015; 42(5):2718-29.

PMID: 25979070 PMC: 4425729. DOI: 10.1118/1.4915542.

References

La Riviere P . Penalized-likelihood sinogram smoothing for low-dose CT. Med Phys. 2005; 32(6):1676-83. DOI: 10.1118/1.1915015. View

Kalra M, Rizzo S, Maher M, Halpern E, Toth T, Shepard J . Chest CT performed with z-axis modulation: scanning protocol and radiation dose. Radiology. 2005; 237(1):303-8. DOI: 10.1148/radiol.2371041227. View

Pflederer T, Rudofsky L, Ropers D, Bachmann S, Marwan M, Daniel W . Image quality in a low radiation exposure protocol for retrospectively ECG-gated coronary CT angiography. AJR Am J Roentgenol. 2009; 192(4):1045-50. DOI: 10.2214/AJR.08.1025. View

Hsieh J . Adaptive streak artifact reduction in computed tomography resulting from excessive x-ray photon noise. Med Phys. 1998; 25(11):2139-47. DOI: 10.1118/1.598410. View

Fahrig R, Dixon R, Payne T, Morin R, Ganguly A, Strobel N . Dose and image quality for a cone-beam C-arm CT system. Med Phys. 2007; 33(12):4541-50. DOI: 10.1118/1.2370508. View

Lauritsch G, Boese J, Wigstrom L, Kemeth H, Fahrig R . Towards cardiac C-arm computed tomography. IEEE Trans Med Imaging. 2006; 25(7):922-34. DOI: 10.1109/tmi.2006.876166. View

Borsdorf A, Kappler S, Raupach R, Hornegger J . Analytic noise-propagation in indirect fan-beam FBP reconstruction. Annu Int Conf IEEE Eng Med Biol Soc. 2009; 2008:2701-4. DOI: 10.1109/IEMBS.2008.4649759. View

Westin C, Richolt J, Moharir V, Kikinis R . Affine adaptive filtering of CT data. Med Image Anal. 2000; 4(2):161-77. DOI: 10.1016/s1361-8415(00)00011-6. View

Kachelriess M, Watzke O, Kalender W . Generalized multi-dimensional adaptive filtering for conventional and spiral single-slice, multi-slice, and cone-beam CT. Med Phys. 2001; 28(4):475-90. DOI: 10.1118/1.1358303. View

10.

Jakobs T, Becker C, Ohnesorge B, Flohr T, Suess C, Schoepf U . Multislice helical CT of the heart with retrospective ECG gating: reduction of radiation exposure by ECG-controlled tube current modulation. Eur Radiol. 2002; 12(5):1081-6. DOI: 10.1007/s00330-001-1278-x. View

11.

Westin C, Wigstrom L, Loock T, Sjoqvist L, Kikinis R, Knutsson H . Three-dimensional adaptive filtering in magnetic resonance angiography. J Magn Reson Imaging. 2001; 14(1):63-71. DOI: 10.1002/jmri.1152. View

12.

Wang J, Li T, Lu H, Liang Z . 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(10):1272-83. PMC: 1619874. DOI: 10.1109/42.896783. View

13.

Hausleiter J, Meyer T, Hadamitzky M, Huber E, Zankl M, Martinoff S . Radiation dose estimates from cardiac multislice computed tomography in daily practice: impact of different scanning protocols on effective dose estimates. Circulation. 2006; 113(10):1305-10. DOI: 10.1161/CIRCULATIONAHA.105.602490. View

14.

Hausleiter J, Meyer T, Hermann F, Hadamitzky M, Krebs M, Gerber T . Estimated radiation dose associated with cardiac CT angiography. JAMA. 2009; 301(5):500-7. DOI: 10.1001/jama.2009.54. View

15.

Gerig G, Kubler O, Kikinis R, Jolesz F . Nonlinear anisotropic filtering of MRI data. IEEE Trans Med Imaging. 1992; 11(2):221-32. DOI: 10.1109/42.141646. View

16.

Zhu L, Wang J, Xing L . Noise suppression in scatter correction for cone-beam CT. Med Phys. 2009; 36(3):741-52. PMC: 2736744. DOI: 10.1118/1.3063001. View

17.

Fahrig R, Fox A, Lownie S, Holdsworth D . Use of a C-arm system to generate true three-dimensional computed rotational angiograms: preliminary in vitro and in vivo results. AJNR Am J Neuroradiol. 1997; 18(8):1507-14. PMC: 8338156. View

18.

Wang J, Li T, Liang Z, Xing L . Dose reduction for kilovotage cone-beam computed tomography in radiation therapy. Phys Med Biol. 2008; 53(11):2897-909. DOI: 10.1088/0031-9155/53/11/009. View

19.

La Riviere P, Bian J, Vargas P . Penalized-likelihood sinogram restoration for computed tomography. IEEE Trans Med Imaging. 2006; 25(8):1022-36. DOI: 10.1109/tmi.2006.875429. View

20.

Borsdorf A, Raupach R, Flohr T, Hornegger J . Wavelet based noise reduction in CT-images using correlation analysis. IEEE Trans Med Imaging. 2008; 27(12):1685-703. DOI: 10.1109/TMI.2008.923983. View