Tags

Type your tag names separated by a space and hit enter

Ultra-low-dose CT image denoising using modified BM3D scheme tailored to data statistics.
Med Phys 2019; 46(1):190-198MP

Abstract

PURPOSE

It is important to enhance image quality for low-dose CT acquisitions to push the ALARA boundary. Current state-of-the-art block-matching three-dimensional (BM3D) denoising scheme assumes white Gaussian noise (WGN) model. This study proposes a novel filtering module to be incorporated into the BM3D framework for ultra-low-dose CT denoising, by accounting for its specific power spectral properties.

METHODS

In the current BM3D algorithm, the Wiener filtering is applied in the transform domain to a post-thresholding signal for enhanced denoising. However, unlike most natural/synthetic images, low-dose CTs do not obey the ideal Gaussian noise model. Based on the specific noise properties of ultra-low-dose CT, we derive the optimal transform-domain coefficients of Wiener filter based on the minimum mean-square-error (MMSE) criterion, taking the noise spectrum and the signal/noise cross spectrum into consideration. In the absence of ground-truth signal, the hard-thresholding denoising module in the previous stage is used as a plug-in estimator. We evaluate the denoising performance on thoracic CT image datasets containing paired full-dose and ultra-low-dose images simulated by a well-validated clinical engine (or pipeline). We also assess its clinical implication by applying the denoising methods to the emphysema quantification task. Our modified BM3D method is compared with the current one, using peak signal-to-noise ratio (PSNR) and emphysema scoring results as evaluation metrics.

RESULTS

The noise in ultra-low-dose CT presented distinct non-Gaussian characteristics and was correlated with image intensity. Performance evaluation showed that the current Wiener filter in basic BM3D algorithm yielded little denoising enhancement on ultra-low-dose CT images. In contrast, the proposed Wiener filter achieved (1.46, 1.91) dB performance gain in mean and median peak signal-to-noise ratio (PSNR) for 5%-dose image denoising and (0.93, 0.95) dB improvement for 10% dose. A paired t-test of the PSNRs between denoising using the current and the proposed Wiener filters demonstrated statistically significant improvement, yielding P-values of 1.45E-12 and 1.34E-7 on 5% and 10%-dose images, respectively. In addition, emphysema quantification on the denoised images using the modified BM3D method also had statistically significant advantage over that using the current BM3D scheme, resulting in a P-value of 6.30E-5 with the commonly used measure.

CONCLUSIONS

This work tailors the Wiener filter in BM3D algorithm to data statistics and demonstrates statistically significant performance improvement on ultra-low-dose CT image denoising and a subsequent emphysema quantification task. Such performance gain is more pronounced with a lower dose level. The development and rationale are generally enough for other image denoising tasks when the WGN assumption is violated.

Authors+Show Affiliations

Department of Radiation Oncology, University of California, Los Angeles, CA, 90095, USA.Departments of Biomedical Physics and Radiology, University of California, Los Angeles, CA, 90095, USA.Departments of Biomedical Physics and Radiology, University of California, Los Angeles, CA, 90095, USA.Department of Radiation Oncology, University of California, Los Angeles, CA, 90095, USA.

Pub Type(s)

Journal Article

Language

eng

PubMed ID

30351450

Citation

Zhao, Tingting, et al. "Ultra-low-dose CT Image Denoising Using Modified BM3D Scheme Tailored to Data Statistics." Medical Physics, vol. 46, no. 1, 2019, pp. 190-198.
Zhao T, Hoffman J, McNitt-Gray M, et al. Ultra-low-dose CT image denoising using modified BM3D scheme tailored to data statistics. Med Phys. 2019;46(1):190-198.
Zhao, T., Hoffman, J., McNitt-Gray, M., & Ruan, D. (2019). Ultra-low-dose CT image denoising using modified BM3D scheme tailored to data statistics. Medical Physics, 46(1), pp. 190-198. doi:10.1002/mp.13252.
Zhao T, et al. Ultra-low-dose CT Image Denoising Using Modified BM3D Scheme Tailored to Data Statistics. Med Phys. 2019;46(1):190-198. PubMed PMID: 30351450.
* Article titles in AMA citation format should be in sentence-case
TY - JOUR T1 - Ultra-low-dose CT image denoising using modified BM3D scheme tailored to data statistics. AU - Zhao,Tingting, AU - Hoffman,John, AU - McNitt-Gray,Michael, AU - Ruan,Dan, Y1 - 2018/11/19/ PY - 2018/04/09/received PY - 2018/10/03/revised PY - 2018/10/09/accepted PY - 2018/10/24/pubmed PY - 2019/1/22/medline PY - 2018/10/24/entrez KW - BM3D KW - data statistics KW - low-dose CT SP - 190 EP - 198 JF - Medical physics JO - Med Phys VL - 46 IS - 1 N2 - PURPOSE: It is important to enhance image quality for low-dose CT acquisitions to push the ALARA boundary. Current state-of-the-art block-matching three-dimensional (BM3D) denoising scheme assumes white Gaussian noise (WGN) model. This study proposes a novel filtering module to be incorporated into the BM3D framework for ultra-low-dose CT denoising, by accounting for its specific power spectral properties. METHODS: In the current BM3D algorithm, the Wiener filtering is applied in the transform domain to a post-thresholding signal for enhanced denoising. However, unlike most natural/synthetic images, low-dose CTs do not obey the ideal Gaussian noise model. Based on the specific noise properties of ultra-low-dose CT, we derive the optimal transform-domain coefficients of Wiener filter based on the minimum mean-square-error (MMSE) criterion, taking the noise spectrum and the signal/noise cross spectrum into consideration. In the absence of ground-truth signal, the hard-thresholding denoising module in the previous stage is used as a plug-in estimator. We evaluate the denoising performance on thoracic CT image datasets containing paired full-dose and ultra-low-dose images simulated by a well-validated clinical engine (or pipeline). We also assess its clinical implication by applying the denoising methods to the emphysema quantification task. Our modified BM3D method is compared with the current one, using peak signal-to-noise ratio (PSNR) and emphysema scoring results as evaluation metrics. RESULTS: The noise in ultra-low-dose CT presented distinct non-Gaussian characteristics and was correlated with image intensity. Performance evaluation showed that the current Wiener filter in basic BM3D algorithm yielded little denoising enhancement on ultra-low-dose CT images. In contrast, the proposed Wiener filter achieved (1.46, 1.91) dB performance gain in mean and median peak signal-to-noise ratio (PSNR) for 5%-dose image denoising and (0.93, 0.95) dB improvement for 10% dose. A paired t-test of the PSNRs between denoising using the current and the proposed Wiener filters demonstrated statistically significant improvement, yielding P-values of 1.45E-12 and 1.34E-7 on 5% and 10%-dose images, respectively. In addition, emphysema quantification on the denoised images using the modified BM3D method also had statistically significant advantage over that using the current BM3D scheme, resulting in a P-value of 6.30E-5 with the commonly used measure. CONCLUSIONS: This work tailors the Wiener filter in BM3D algorithm to data statistics and demonstrates statistically significant performance improvement on ultra-low-dose CT image denoising and a subsequent emphysema quantification task. Such performance gain is more pronounced with a lower dose level. The development and rationale are generally enough for other image denoising tasks when the WGN assumption is violated. SN - 2473-4209 UR - https://www.unboundmedicine.com/medline/citation/30351450/Ultra_low_dose_CT_image_denoising_using_modified_BM3D_scheme_tailored_to_data_statistics_ L2 - https://doi.org/10.1002/mp.13252 DB - PRIME DP - Unbound Medicine ER -