1
|
Liu SZ, Herbst M, Schaefer J, Weber T, Vogt S, Ritschl L, Kappler S, Kawcak CE, Stewart HL, Siewerdsen JH, Zbijewski W. Feasibility of bone marrow edema detection using dual-energy cone-beam computed tomography. Med Phys 2024; 51:1653-1673. [PMID: 38323878 DOI: 10.1002/mp.16962] [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: 05/04/2023] [Revised: 12/17/2023] [Accepted: 01/16/2024] [Indexed: 02/08/2024] Open
Abstract
BACKGROUND Dual-energy (DE) detection of bone marrow edema (BME) would be a valuable new diagnostic capability for the emerging orthopedic cone-beam computed tomography (CBCT) systems. However, this imaging task is inherently challenging because of the narrow energy separation between water (edematous fluid) and fat (health yellow marrow), requiring precise artifact correction and dedicated material decomposition approaches. PURPOSE We investigate the feasibility of BME assessment using kV-switching DE CBCT with a comprehensive CBCT artifact correction framework and a two-stage projection- and image-domain three-material decomposition algorithm. METHODS DE CBCT projections of quantitative BME phantoms (water containers 100-165 mm in size with inserts presenting various degrees of edema) and an animal cadaver model of BME were acquired on a CBCT test bench emulating the standard wrist imaging configuration of a Multitom Rax twin robotic x-ray system. The slow kV-switching scan protocol involved a 60 kV low energy (LE) beam and a 120 kV high energy (HE) beam switched every 0.5° over a 200° angular span. The DE CBCT data preprocessing and artifact correction framework consisted of (i) projection interpolation onto matched LE and HE projections views, (ii) lag and glare deconvolutions, and (iii) efficient Monte Carlo (MC)-based scatter correction. Virtual non-calcium (VNCa) images for BME detection were then generated by projection-domain decomposition into an Aluminium (Al) and polyethylene basis set (to remove beam hardening) followed by three-material image-domain decomposition into water, Ca, and fat. Feasibility of BME detection was quantified in terms of VNCa image contrast and receiver operating characteristic (ROC) curves. Robustness to object size, position in the field of view (FOV) and beam collimation (varied 20-160 mm) was investigated. RESULTS The MC-based scatter correction delivered > 69% reduction of cupping artifacts for moderate to wide collimations (> 80 mm beam width), which was essential to achieve accurate DE material decomposition. In a forearm-sized object, a 20% increase in water concentration (edema) of a trabecular bone-mimicking mixture presented as ∼15 HU VNCa contrast using 80-160 mm beam collimations. The variability with respect to object position in the FOV was modest (< 15% coefficient of variation). The areas under the ROC curve were > 0.9. A femur-sized object presented a somewhat more challenging task, resulting in increased sensitivity to object positioning at 160 mm collimation. In animal cadaver specimens, areas of VNCa enhancement consistent with BME were observed in DE CBCT images in regions of MRI-confirmed edema. CONCLUSION Our results indicate that the proposed artifact correction and material decomposition pipeline can overcome the challenges of scatter and limited spectral separation to achieve relatively accurate and sensitive BME detection in DE CBCT. This study provides an important baseline for clinical translation of musculoskeletal DE CBCT to quantitative, point-of-care bone health assessment.
Collapse
Affiliation(s)
- Stephen Z Liu
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
| | | | | | | | | | | | | | - Christopher E Kawcak
- Department of Clinical Sciences, Colorado State University College of Veterinary Medicine and Biomedical Sciences, Fort Collins, Colorado, USA
| | - Holly L Stewart
- Department of Clinical Sciences, Colorado State University College of Veterinary Medicine and Biomedical Sciences, Fort Collins, Colorado, USA
| | - Jeffrey H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
- Department of Imaging Physics, MD Anderson Cancer Center, Houston, Texas, USA
| | - Wojciech Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
| |
Collapse
|
2
|
Martins JC, Maier J, Gianoli C, Neppl S, Dedes G, Alhazmi A, Veloza S, Reiner M, Belka C, Kachelrieß M, Parodi K. Towards real-time EPID-based 3D in vivo dosimetry for IMRT with Deep Neural Networks: A feasibility study. Phys Med 2023; 114:103148. [PMID: 37801811 DOI: 10.1016/j.ejmp.2023.103148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/12/2023] [Revised: 08/17/2023] [Accepted: 09/22/2023] [Indexed: 10/08/2023] Open
Abstract
We investigate the potential of the Deep Dose Estimate (DDE) neural network to predict 3D dose distributions inside patients with Monte Carlo (MC) accuracy, based on transmitted EPID signals and patient CTs. The network was trained using as input patient CTs and first-order dose approximations (FOD). Accurate dose distributions (ADD) simulated with MC were given as training targets. 83 pelvic CTs were used to simulate ADDs and respective EPID signals for subfields of prostate IMRT plans (gantry at 0∘). FODs were produced as backprojections from the EPID signals. 581 ADD-FOD sets were produced and divided into training and test sets. An additional dataset simulated with gantry at 90∘ (lateral set) was used for evaluating the performance of the DDE at different beam directions. The quality of the FODs and DDE-predicted dose distributions (DDEP) with respect to ADDs, from the test and lateral sets, was evaluated with gamma analysis (3%,2 mm). The passing rates between FODs and ADDs were as low as 46%, while for DDEPs the passing rates were above 97% for the test set. Meaningful improvements were also observed for the lateral set. The high passing rates for DDEPs indicate that the DDE is able to convert FODs into ADDs. Moreover, the trained DDE predicts the dose inside a patient CT within 0.6 s/subfield (GPU), in contrast to 14 h needed for MC (CPU-cluster). 3D in vivo dose distributions due to clinical patient irradiation can be obtained within seconds, with MC-like accuracy, potentially paving the way towards real-time EPID-based in vivo dosimetry.
Collapse
Affiliation(s)
- Juliana Cristina Martins
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| | - Joscha Maier
- German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg, 69120, Germany.
| | - Chiara Gianoli
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| | - Sebastian Neppl
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, Munich, 81377, Germany.
| | - George Dedes
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| | - Abdulaziz Alhazmi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| | - Stella Veloza
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| | - Michael Reiner
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, Munich, 81377, Germany.
| | - Claus Belka
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, Munich, 81377, Germany.
| | - Marc Kachelrieß
- German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg, 69120, Germany; Heidelberg University, Grabengasse 1, Heidelberg, 69117, Germany.
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, Garching b. München, 85748, Germany.
| |
Collapse
|
3
|
Zhang X, Jiang Y, Luo C, Li D, Niu T, Yu G. Image-based scatter correction for cone-beam CT using flip swin transformer U-shape network. Med Phys 2023; 50:5002-5019. [PMID: 36734321 DOI: 10.1002/mp.16277] [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: 05/03/2022] [Revised: 12/23/2022] [Accepted: 01/23/2023] [Indexed: 02/04/2023] Open
Abstract
BACKGROUND Cone beam computed tomography (CBCT) plays an increasingly important role in image-guided radiation therapy. However, the image quality of CBCT is severely degraded by excessive scatter contamination, especially in the abdominal region, hindering its further applications in radiation therapy. PURPOSE To restore low-quality CBCT images contaminated by scatter signals, a scatter correction algorithm combining the advantages of convolutional neural networks (CNN) and Swin Transformer is proposed. METHODS In this paper a scatter correction model for CBCT image, the Flip Swin Transformer U-shape network (FSTUNet) model, is proposed. In this model, the advantages of CNN in texture detail and Swin Transformer in global correlation are used to accurately extract shallow and deep features, respectively. Instead of using the original Swin Transformer tandem structure, we build the Flip Swin Transformer Block to achieve a more powerful inter-window association extraction. The validity and clinical relevance of the method is demonstrated through extensive experiments on a Monte Carlo (MC) simulation dataset and frequency split dataset generated by a validated method, respectively. RESULT Experimental results on the MC simulated dataset show that the root mean square error of images corrected by the method is reduced from over 100 HU to about 7 HU. Both the structural similarity index measure (SSIM) and the universal quality index (UQI) are close to 1. Experimental results on the frequency split dataset demonstrate that the method not only corrects shading artifacts but also exhibits a high degree of structural consistency. In addition, comparison experiments show that FSTUNet outperforms UNet, Deep Residual Convolutional Neural Network (DRCNN), DSENet, Pix2pixGAN, and 3DUnet methods in both qualitative and quantitative metrics. CONCLUSIONS Accurately capturing the features at different levels is greatly beneficial for reconstructing high-quality scatter-free images. The proposed FSTUNet method is an effective solution to CBCT scatter correction and has the potential to improve the accuracy of CBCT image-guided radiation therapy.
Collapse
Affiliation(s)
- Xueren Zhang
- Shandong Key Laboratory of Medical Physics and Image Processing, Shandong Institute of Industrial Technology for Health Sciences and Precision Medicine, School of Physics and Electronics, Shandong Normal University, Jinan, Shandong, China
| | - Yangkang Jiang
- Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Institute of Translational Medicine, Zhejiang University, Hangzhou, Zhejiang, China
| | - Chen Luo
- Shenzhen Bay Laboratory, Shenzhen, China
- School of Automation, Zhejiang Institute of Mechanical & Electrical Engineering, Hangzhou, China
| | - Dengwang Li
- Shandong Key Laboratory of Medical Physics and Image Processing, Shandong Institute of Industrial Technology for Health Sciences and Precision Medicine, School of Physics and Electronics, Shandong Normal University, Jinan, Shandong, China
| | - Tianye Niu
- Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Institute of Translational Medicine, Zhejiang University, Hangzhou, Zhejiang, China
- Shenzhen Bay Laboratory, Shenzhen, China
| | - Gang Yu
- Shandong Key Laboratory of Medical Physics and Image Processing, Shandong Institute of Industrial Technology for Health Sciences and Precision Medicine, School of Physics and Electronics, Shandong Normal University, Jinan, Shandong, China
| |
Collapse
|
4
|
Itoh T, Noguchi K. Evaluation of the quantitative performance of non-enhanced dual-energy CT X-map in detecting acute ischemic brain stroke: A model observer study using computer simulation. Phys Med 2022; 104:85-92. [PMID: 36371946 DOI: 10.1016/j.ejmp.2022.10.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 09/02/2022] [Accepted: 10/30/2022] [Indexed: 11/11/2022] Open
Abstract
PURPOSE A simulation study was performed to evaluate the quantitative performance of X-map images-derived from non-enhanced (NE) dual-energy computed tomography (DECT)-in detecting acute ischemic stroke (AIS) compared with that of NE-DECT mixed images. METHODS A virtual phantom, 150 mm in diameter, filled with tissues comprising various gray- and white-matter proportions was used to generate pairs of NE-head images at 80 kV and Sn150 kV at three dose levels (20, 40, and 60 mGy). The phantom included an inserted low-contrast object, 15 mm in diameter, with four densities (0%, 5%, 10%, and 15%) mimicking ischemic edema. Mixed and X-map images were generated from these sets of images and compared in terms of detectability of ischemic edema using a channelized Hotelling observer (CHO). The area under the curve (AUC) of the receiver operating characteristic that generated CHO for each condition was used as a figure of merit. RESULTS The AUCs of X-map images were always significantly higher than those of mixed images (P < 0.001). The improvement in AUC for X-map images compared with that for mixed images at edema densities was 9.2%-12.6% at 20 mGy, 10.1%-17.7% at 40 mGy, and 14.0%-19.4% at 60 mGy. At any edema density, X-map images at 20 mGy resulted in higher AUCs than mixed images acquired at any other dose level (P < 0.001), which corresponded to a 66% dose reduction on X-map images. CONCLUSIONS The simulation study confirmed that NE-DECT X-map images have superior capability of detecting AIS than NE-DECT mixed images.
Collapse
Affiliation(s)
- Toshihide Itoh
- Department of CT Research and Collaboration, Siemens Healthineers, 1-11-1 Osaki, Shinagawa, Tokyo 141-8644, Japan.
| | - Kyo Noguchi
- Department of Radiology, Graduate School of Medicine and Pharmaceutical Science, University of Toyama, 2630 Sugitani, Toyama city, Toyama 930-0194, Japan
| |
Collapse
|
5
|
Liu SZ, Tivnan M, Osgood GM, Siewerdsen JH, Stayman JW, Zbijewski W. Model-based three-material decomposition in dual-energy CT using the volume conservation constraint. Phys Med Biol 2022; 67:10.1088/1361-6560/ac7a8b. [PMID: 35724658 PMCID: PMC9297826 DOI: 10.1088/1361-6560/ac7a8b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Accepted: 06/20/2022] [Indexed: 01/13/2023]
Abstract
Objective. We develop a model-based optimization algorithm for 'one-step' dual-energy (DE) CT decomposition of three materials directly from projection measurements.Approach.Since the three-material problem is inherently undetermined, we incorporate the volume conservation principle (VCP) as a pair of equality and nonnegativity constraints into the objective function of the recently reported model-based material decomposition (MBMD). An optimization algorithm (constrained MBMD, CMBMD) is derived that utilizes voxel-wise separability to partition the volume into a VCP-constrained region solved using interior-point iterations, and an unconstrained region (air surrounding the object, where VCP is violated) solved with conventional two-material MBMD. Constrained MBMD (CMBMD) is validated in simulations and experiments in application to bone composition measurements in the presence of metal hardware using DE cone-beam CT (CBCT). A kV-switching protocol with non-coinciding low- and high-energy (LE and HE) projections was assumed. CMBMD with decomposed base materials of cortical bone, fat, and metal (titanium, Ti) is compared to MBMD with (i) fat-bone and (ii) fat-Ti bases.Main results.Three-material CMBMD exhibits a substantial reduction in metal artifacts relative to the two-material MBMD implementations. The accuracies of cortical bone volume fraction estimates are markedly improved using CMBMD, with ∼5-10× lower normalized root mean squared error in simulations with anthropomorphic knee phantoms (depending on the complexity of the metal component) and ∼2-2.5× lower in an experimental test-bench study.Significance.In conclusion, we demonstrated one-step three-material decomposition of DE CT using volume conservation as an optimization constraint. The proposed method might be applicable to DE applications such as bone marrow edema imaging (fat-bone-water decomposition) or multi-contrast imaging, especially on CT/CBCT systems that do not provide coinciding LE and HE ray paths required for conventional projection-domain DE decomposition.
Collapse
Affiliation(s)
- Stephen Z. Liu
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Matthew Tivnan
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Greg M. Osgood
- Department of Orthopaedic Surgery, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Jeffrey H. Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - J. Webster Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Wojciech Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| |
Collapse
|
6
|
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
|