1
|
Ge T, Liao R, Medrano M, Politte DG, Williamson JF, O’Sullivan JA. MB-DECTNet: a model-based unrolling network for accurate 3D dual-energy CT reconstruction from clinically acquired helical scans. Phys Med Biol 2023; 68:245009. [PMID: 37802071 PMCID: PMC10714406 DOI: 10.1088/1361-6560/ad00fb] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 09/11/2023] [Accepted: 10/06/2023] [Indexed: 10/08/2023]
Abstract
Objective.Over the past several decades, dual-energy CT (DECT) imaging has seen significant advancements due to its ability to distinguish between materials. DECT statistical iterative reconstruction (SIR) has exhibited potential for noise reduction and enhanced accuracy. However, its slow convergence and substantial computational demands render the elapsed time for 3D DECT SIR often clinically unacceptable. The objective of this study is to accelerate 3D DECT SIR while maintaining subpercentage or near-subpercentage accuracy.Approach.We incorporate DECT SIR into a deep-learning model-based unrolling network for 3D DECT reconstruction (MB-DECTNet), which can be trained end-to-end. This deep learning-based approach is designed to learn shortcuts between initial conditions and the stationary points of iterative algorithms while preserving the unbiased estimation property of model-based algorithms. MB-DECTNet comprises multiple stacked update blocks, each containing a data consistency layer (DC) and a spatial mixer layer, with the DC layer functioning as a one-step update from any traditional iterative algorithm.Main results.The quantitative results indicate that our proposed MB-DECTNet surpasses both the traditional image-domain technique (MB-DECTNet reduces average bias by a factor of 10) and a pure deep learning method (MB-DECTNet reduces average bias by a factor of 8.8), offering the potential for accurate attenuation coefficient estimation, akin to traditional statistical algorithms, but with considerably reduced computational costs. This approach achieves 0.13% bias and 1.92% mean absolute error and reconstructs a full image of a head in less than 12 min. Additionally, we show that the MB-DECTNet output can serve as an initializer for DECT SIR, leading to further improvements in results.Significance.This study presents a model-based deep unrolling network for accurate 3D DECT reconstruction, achieving subpercentage error in estimating virtual monoenergetic images for a full head at 60 and 150 keV in 30 min, representing a 40-fold speedup compared to traditional approaches. These findings have significant implications for accelerating DECT SIR and making it more clinically feasible.
Collapse
Affiliation(s)
- Tao Ge
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| | - Rui Liao
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| | - Maria Medrano
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| | - David G Politte
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| | - Jeffrey F Williamson
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| | - Joseph A O’Sullivan
- Washington University in St. Louis, Saint Louis, MO 63130, United States of America
| |
Collapse
|
2
|
Zhang G, Wang Y, Chen W, Li T, Tian Y. Correction of Bowtie filter induced scatter signals based on air scan data and object scan data. Biomed Phys Eng Express 2022; 8. [PMID: 35276688 DOI: 10.1088/2057-1976/ac5d0c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 03/11/2022] [Indexed: 11/12/2022]
Abstract
In a cone beam CT system, a bowtie filter brings in additional scatter signals with respect to object induced scatter signals, which can degrade image quality and sometimes result in artifacts. This work aims to improve the image quality of CT scans by analyzing the contribution of bowtie filter induced scatter signals and removing them from projection data. Air calibration is a very useful preprocessing step to eliminate the response variations of detector pixels. Bowtie filter induced scattered x-ray signals of air scans are recorded in air calibration tables and therefore considered as a part of primary signals. However, scattered X-rays behave differently in scanned objects compared to primary x-rays. The difference should be corrected to eliminate the impact of bowtie filter induced scatter signals. A kernel based correction algorithm based on air scan data, named bowtie filter scatter correction algorithm, is applied to estimate and to eliminate the bowtie filter induced scatter signals in object scans. The scatter signals of air scans can be measured with air scans or retrieved from air calibration tables of a CT system, and can be used as input of the correction algorithm to estimate the change of scatter signals caused by the scanned objects in the scan field. Based on the assumption that the scatter signals in the projection data scanned with narrow collimation can be neglected, the difference signals between narrow and broad collimations can be used to estimate bowtie filter induced scatter signals for air scans with the correction of extra-focal radiations (EFRs). The calculated bowtie filter induced scatter signals have been compared with the results of Monte Carlo simulations, and the parameters of correction algorithm have been determined by fitting the measured scatter signal curves of phantom scans with calculated curves. Projection data have been reconstructed using Filtered BackProjection (FBP) method with and without bowtie filter correction to check whether the image quality is improved. Scatter signals can be well approximated with the bowtie filter scatter correction algorithm together with an existing object scatter correction algorithm. After removing the bowtie filter induced scatter signals, the dark bands in reconstructed images in the regions near the edges of scanned objects can be mostly eliminated. The difference signals of air scan data between narrow and broad collimations can be used to estimate the bowtie filter induced scatter for air scans. The proposed bowtie filter scatter correction algorithm using air scan data can be applied to estimate and to remove most of the bowtie filter induced scatter signals in object scans.
Collapse
Affiliation(s)
- Guoqing Zhang
- Siemens Shanghai Medical Equipment Ltd, 278 Zhouzhu Rd., Pudong District, Shanghai 201318, People's Republic of China
| | - Yang Wang
- Siemens Shanghai Medical Equipment Ltd, 278 Zhouzhu Rd., Pudong District, Shanghai 201318, People's Republic of China
| | - Wenhao Chen
- Siemens Shanghai Medical Equipment Ltd, 278 Zhouzhu Rd., Pudong District, Shanghai 201318, People's Republic of China
| | - Taotao Li
- Siemens Shanghai Medical Equipment Ltd, 278 Zhouzhu Rd., Pudong District, Shanghai 201318, People's Republic of China
| | - Yi Tian
- Siemens Shanghai Medical Equipment Ltd, 278 Zhouzhu Rd., Pudong District, Shanghai 201318, People's Republic of China
| |
Collapse
|
3
|
Haase V, Hahn K, Schöndube H, Stierstorfer K, Maier A, Noo F. Single material beam hardening correction via an analytical energy response model for diagnostic CT. Med Phys 2022; 49:5014-5037. [PMID: 35651302 PMCID: PMC9388575 DOI: 10.1002/mp.15787] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 05/17/2022] [Accepted: 05/18/2022] [Indexed: 11/07/2022] Open
Abstract
BACKGROUND Various clinical studies show the potential for a wider quantitative role of diagnostic X-ray computed tomography (CT) beyond size measurements. Currently, the clinical use of attenuation values is however limited due to their lack of robustness. This issue can be observed even on the same scanner across patient size and positioning. There are different causes for the lack of robustness in the attenuation values; one possible source of error is beam hardening of the X-ray source spectrum. The conventional and well-established approach to address this issue is a calibration-based single material beam hardening correction (BHC) using a water cylinder. PURPOSE We investigate an alternative approach for single material BHC with the aim of producing a more robust result for the attenuation values. The underlying hypothesis of this investigation is that calibration based BHC automatically corrects for scattered radiation in a manner that is sub-optimal in terms of bias as soon as the scanned object strongly deviates from the water cylinder used for calibration. METHODS The approach we propose performs BHC via an analytical energy response model that is embedded into a correction pipeline that efficiently estimates and subtracts scattered radiation in a patient-specific manner prior to BHC. The estimation of scattered radiation is based on minimizing, in average, the squared difference between our corrected data and the vendor-calibrated data. The used energy response model is considering the spectral effects of the detector response and of the pre-filtration of the source spectrum including a beam-shaping bowtie filter. The performance of the correction pipeline is first characterized with computer simulated data. Afterwards, it is tested using real 3-D CT data sets of two different phantoms, with various kV settings and phantom positions, assuming a circular data acquisition. The results are compared in the image domain to those from the scanner. RESULTS For experiments with a water cylinder, the proposed correction pipeline leads to similar results as the vendor. For reconstructions of a QRM liver phantom with extension ring, the proposed correction pipeline achieved a more uniform and stable outcome in the attenuation values of homogeneous materials within the phantom. For example, the root mean squared deviation between centered and off-centered phantom positioning was reduced from 6.6 HU to 1.8 HU in one profile. CONCLUSIONS We have introduced a patient-specific approach for single material BHC in diagnostic CT via the use of an analytical energy response model. This approach shows promising improvements in terms of robustness of attenuation values for large patient sizes. Our results contribute towards improving CT images so as to make CT attenuation values more reliable for use in clinical practice. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Viktor Haase
- Siemens Healthcare GmbH, Siemensstr. 3, Forchheim, 91301, Germany.,Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander-Universität Erlangen-Nürnberg, Martensstr. 3, Erlangen, 91058, Germany
| | - Katharina Hahn
- Siemens Healthcare GmbH, Siemensstr. 3, Forchheim, 91301, Germany
| | - Harald Schöndube
- Siemens Healthcare GmbH, Siemensstr. 3, Forchheim, 91301, Germany
| | | | - Andreas Maier
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander-Universität Erlangen-Nürnberg, Martensstr. 3, Erlangen, 91058, Germany
| | - Frédéric Noo
- Department of Radiology and Imaging Sciences, University of Utah, 729 Arapeen Drive, Salt Lake City, Utah, 84108, USA
| |
Collapse
|
4
|
Medrano M, Liu R, Zhao T, Webb T, Politte DG, Whiting BR, Liao R, Ge T, Porras-Chaverri MA, O'Sullivan JA, Williamson JF. Towards sub-percentage uncertainty proton stopping-power mapping via dual-energy CT: direct experimental validation and uncertainty analysis of a statistical iterative image reconstruction method. Med Phys 2022; 49:1599-1618. [PMID: 35029302 DOI: 10.1002/mp.15457] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 10/28/2021] [Accepted: 12/22/2021] [Indexed: 11/06/2022] Open
Abstract
PURPOSE To assess the potential of a joint dual-energy CT reconstruction process (Statistical Image Reconstruction method built on a Basis Vector Model (JSIR-BVM)) implemented on a 16-slice commercial CT scanner to measure high spatial-resolution stopping-power ratio (SPR) maps with uncertainties of less than 1%. METHODS JSIR-BVM was used to reconstruct images of effective electron density and mean excitation energy from dual-energy CT (DECT) sinograms for ten high-purity samples of known density and atomic composition inserted into head and body phantoms. The measured DECT data consisted of 90 kVp and 140 kVp axial sinograms serially acquired on a Philips Brilliance Big Bore CT scanner without beam-hardening corrections. The corresponding SPRs were subsequently measured directly via ion chamber measurements on a MEVION S250 superconducting synchrocyclotron and evaluated theoretically from the known sample compositions and densities. Deviations of JSIR-BVM SPR values from their theoretically calculated and directly measured ground-truth values were evaluated for our JSIR-BVM method and for our implementation of the Hünemohr-Saito (H-S) DECT image-domain decomposition technique for SPR imaging. A thorough uncertainty analysis was then performed for 5 different scenarios (comparison of JSIR-BVM SPR/SP to International Commission on Radiation Measurements and Units (ICRU) benchmarks; comparison of JSIR-BVM SPR to measured benchmarks; and uncertainties in JSIR-BVM SPR/SP maps for patients of unknown composition) per the Joint Committee for Guides in Metrology (JCGM) and the Guide to expression of Uncertainty in Measurement (GUM), including the impact of uncertainties in measured photon spectra, sample composition and density, photon cross-section and I-value models, and random measurement uncertainty. Estimated SPR uncertainty for three main tissue groups in patients of unknown composition and the weighted proportion of each tissue type for three proton treatment sites were then used to derive a composite range uncertainty for our method. RESULTS Mean JSIR-BVM SPR estimates deviated by less than 1% from their theoretical and directly measured ground-truth values for most inserts and phantom geometries except for high density Delrin and Teflon samples with SPR error relative to proton measurements of 1.1% and -1.0% (Head Phantom) and 1.1% and -1.1% (Body Phantom). The overall RMS deviations over all samples were 0.39% and 0.52% (head phantom) and 0.43% and 0.57% (body phantom) relative to theoretical and directly measured ground-truth SPRs, respectively. The corresponding RMS (maximum) errors for the image-domain decomposition method were 2.68% and 2.73% (4.68% and 4.99%) for the head phantom and 0.71% and 0.87% (1.37% and 1.66%) for the body phantom. Compared to H-S SPR maps, JSIR-BVM yielded 30% sharper and two-fold sharper images for soft tissues and bone-like surrogates, respectively, while reducing noise by factors of 6 and 3, respectively. The uncertainty (coverage factor k = 1) of the DECT-to-benchmark values comparison ranged from 0.5% to 1.5% and is dominated by scanning-beam photon-spectra uncertainties. An analysis of the SPR uncertainty for patients of unknown composition showed a JSIR-BVM uncertainty of 0.65%, 1.21%, and 0.77% for soft-, lung-, and bony-tissue groups which led to a composite range uncertainty of 0.6%-0.9%. CONCLUSIONS Observed JSIR-BVM SPR estimation errors were all less than 50% of the estimated k = 1 total uncertainty of our benchmarking experiment, demonstrating that JSIR-BVM high spatial-resolution, low-noise SPR mapping is feasible and is robust to variations in the geometry of the scanned object. In contrast, the much larger H-S SPR estimation errors are dominated by imaging noise and residual beam-hardening artifacts. While the uncertainties characteristic of our current JSIR-BVM implementation can be as large as 1.5%, achieving <1% total uncertainty is feasible by improving the accuracy of scanner-specific scatter-profile and photon-spectrum estimates. With its robustness to beam-hardening artifact, image noise and variations in phantom size and geometry, JSIR-BVM has the potential to achieve high spatial-resolution SPR mapping with sub-percentage accuracy and estimated uncertainty in the clinical setting. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Maria Medrano
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Ruirui Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Tianyu Zhao
- Department of Radiation Oncology, Washington University in St. Louis, St. Louis, MO, 63110, USA
| | - Tyler Webb
- Department of Physics, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - David G Politte
- Mallinckrodt Institute of Radiology, St. Louis, MO, 63110, USA
| | - Bruce R Whiting
- Department of Radiology, University of Pittsburgh, Pittsburgh, PA, 15213, USA
| | - Rui Liao
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Tao Ge
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Mariela A Porras-Chaverri
- Atomic, Nuclear and Molecular Sciences Research Center (CICANUM), University of Costa Rica, San Jose, Costa Rica
| | - Joseph A O'Sullivan
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Jeffrey F Williamson
- Department of Radiation Oncology, Washington University in St. Louis, St. Louis, MO, 63110, USA
| |
Collapse
|
5
|
Liu R, Lei Y, Wang T, Zhou J, Roper J, Lin L, McDonald MW, Bradley JD, Curran WJ, Liu T, Yang X. Synthetic dual-energy CT for MRI-only based proton therapy treatment planning using label-GAN. Phys Med Biol 2021; 66:065014. [PMID: 33596558 DOI: 10.1088/1361-6560/abe736] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
MRI-only treatment planning is highly desirable in the current proton radiation therapy workflow due to its appealing advantages such as bypassing MR-CT co-registration, avoiding x-ray CT exposure dose and reduced medical cost. However, MRI alone cannot provide stopping power ratio (SPR) information for dose calculations. Given that dual energy CT (DECT) can estimate SPR with higher accuracy than conventional single energy CT, we propose a deep learning-based method in this study to generate synthetic DECT (sDECT) from MRI to calculate SPR. Since the contrast difference between high-energy and low-energy CT (LECT) is important, and in order to accurately model this difference, we propose a novel label generative adversarial network-based model which can not only discriminate the realism of sDECT but also differentiate high-energy CT (HECT) and LECT from DECT. A cohort of 57 head-and-neck cancer patients with DECT and MRI pairs were used to validate the performance of the proposed framework. The results of sDECT and its derived SPR maps were compared with clinical DECT and the corresponding SPR, respectively. The mean absolute error for synthetic LECT and HECT were 79.98 ± 18.11 HU and 80.15 ± 16.27 HU, respectively. The corresponding SPR maps generated from sDECT showed a normalized mean absolute error as 5.22% ± 1.23%. By comparing with the traditional Cycle GANs, our proposed method significantly improves the accuracy of sDECT. The results indicate that on our dataset, the sDECT image form MRI is close to planning DECT, and thus shows promising potential for generating SPR maps for proton therapy.
Collapse
Affiliation(s)
- Ruirui Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Yang Lei
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Tonghe Wang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Jun Zhou
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Justin Roper
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Liyong Lin
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Mark W McDonald
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Jeffrey D Bradley
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Walter J Curran
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Tian Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| | - Xiaofeng Yang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322, United States of America
| |
Collapse
|