Multienergy element-resolved cone beam CT (MEER-CBCT) realized on a conventional CBCT platform

Machine Learning Approach and Model for Predicting Proton Stopping Power Ratio and Other Parameters Using Computed Tomography Images

Purpose

The purpose of this study was to accurately estimate proton stopping power ratio (SPR), relative electron density ρ e, effective atomic number (Z eff), and mean excitation energy (I) using one simple robust model and design a machine learning algorithm that will lead to automation.

Methods

Empirical relationships between computed tomography (CT) number and SPR, ρ e (Z eff) and I were used to formulate a model that predicts all the four parameters using linear attenuation coefficients which can be converted to CT numbers. The results of these models were compared with the results of other existing models. Thirty-three ICRU human tissues were used as modeling data and 12 Gammex inserts as testing data for the machine learning algorithm designed. More ways of tissue classification were introduced to improve accuracy. In the examples, the dual energy methods were implemented using 80 kVp and 150 kVP/Sn.

Results

The proposed method gave modeling root mean square error (RMSE) near 1% at maximum for the case of SPR and ρ e for both single and dual-energy CT approaches considered with modeling RMSE of 0.32% for ρ e and 0.38% for SPR as modeling RMSE with room for improvement (this can be done by adjusting the model number of terms as well as the parameters). The method was able to achieve modeling RMSE of 1.11% for I and 1.66% for Z ef f. The mean error for all the estimated quantities was near 0.00%. In most cases, the proposed method has lower testing RMSE and mean error compare to the other methods presented in the study.

Conclusion

The proposed method proves to be more flexible and robust among all presented methods since it has lower testing error in most cases and can be improved based on data using the machine learning algorithm. The algorithm can also improve estimation by adjusting the model as well as aid in automation and it's easy to implement.

Advancing the Collaboration Between Imaging and Radiation Oncology

The fusion of cutting-edge imaging technologies with radiation therapy (RT) has catalyzed transformative breakthroughs in cancer treatment in recent decades. It is critical for us to review our achievements and preview into the next phase for future synergy between imaging and RT. This paper serves as a review and preview for fostering collaboration between these two domains in the forthcoming decade. Firstly, it delineates ten prospective directions ranging from technological innovations to leveraging imaging data in RT planning, execution, and preclinical research. Secondly, it presents major directions for infrastructure and team development in facilitating interdisciplinary synergy and clinical translation. We envision a future where seamless integration of imaging technologies into RT will not only meet the demands of RT but also unlock novel functionalities, enhancing accuracy, efficiency, safety, and ultimately, the standard of care for patients worldwide.

Estimation of Proton Stopping Power Ratio and Mean Excitation Energy Using Electron Density and Its Applications via Machine Learning Approach

Purpose

The purpose of this study was to develop a simple flexible method for accurate estimation of stopping power ratio (SPR) and mean excitation energy (I) using relative electron density (ρ e).

Materials and Methods

The model was formulated using empirical relationships between SPR, mean excitation energy I, and relative electron density. Some examples were implemented, and a comparison was carried out using other existing methods. The needed coefficients in the model were estimated using optimization tools. Basis vector method (BVM) and Hunemohr and Saito (H-S) method were applied to estimate the ρ e used in the application section. 80 kVp and 150 kVpSn were used as low and high energy, respectively, for the implementation of dual-energy methods.

Results

All the examples of the proposed method considered have modeling error that is ≤0.32% and testing root mean square error (RMSE) ≤0.92% for SPR with a mean error close to 0.00%. The method was able to achieve modeling RMSE of 2.12% for mean excitation energy with room for improvement. Similar or better results were achieved in application to BVM.

Conclusion

The method showed robustness in application by achieving lower testing error than other presented methods in most cases. It achieved accurate estimation which can be improved using the machine learning algorithm since it is flexible to implement in terms of the function (model) degree and tissue classification.

Super-resolution dual-layer CBCT imaging with model-guided deep learning

Objective.This study aims at investigating a novel super resolution CBCT imaging approach with a dual-layer flat panel detector (DL-FPD).Approach.With DL-FPD, the low-energy and high-energy projections acquired from the top and bottom detector layers contain over-sampled spatial information, from which super-resolution CT images can be reconstructed. A simple mathematical model is proposed to explain the signal formation procedure in DL-FPD, and a dedicated recurrent neural network, named suRi-Net, is developed based upon the above imaging model to nonlinearly retrieve the high-resolution dual-energy information. Physical benchtop experiments are conducted to validate the performance of this newly developed super-resolution CBCT imaging method.Main Results.The results demonstrate that the proposed suRi-Net can accurately retrieve high spatial resolution information from the low-energy and high-energy projections of low spatial resolution. Quantitatively, the spatial resolution of the reconstructed CBCT images from the top and bottom detector layers is increased by about 45% and 54%, respectively.Significance.In the future, suRi-Net will provide a new approach to perform high spatial resolution dual-energy imaging in DL-FPD-based CBCT systems.

Virtual cone-beam computed tomography simulator with human phantom library and its application to the elemental material decomposition

PURPOSE

The purpose of this study is to develop a virtual CBCT simulator with a head and neck (HN) human phantom library and to demonstrate the feasibility of elemental material decomposition (EMD) for quantitative CBCT imaging using this virtual simulator.

METHODS

The library of 36 HN human phantoms were developed by extending the ICRP 110 adult phantoms based on human age, height, and weight statistics. To create the CBCT database for the library, a virtual CBCT simulator that simulated the direct and scattered X-ray on a flat panel detector using ray-tracing and deep-learning (DL) models was used. Gaussian distributed noise was also included on the flat panel detector, which was evaluated using a real CBCT system. The usefulness of the virtual CBCT system was demonstrated through the application of the developed DL-based EMD model for case involving virtual phantom and real patient.

RESULTS

The virtual simulator could generate various virtual CBCT images based on the human phantom library, and the prediction of the EMD could be successfully performed by preparing the CBCT database from the proposed virtual system, even for a real patient. The CBCT image degradation owing to the scattered X-ray and the statistical noise affected the prediction accuracy, although these effects were minimal. Furthermore, the elemental distribution using the real CBCT image was also predictable.

CONCLUSIONS

This study demonstrated the potential of using computer vision for medical data preparation and analysis, which could have important implications for improving patient outcomes, especially in adaptive radiation therapy.

Dual- and multi-energy CT for particle stopping-power estimation: current state, challenges and potential

Range uncertainty has been a key factor preventing particle radiotherapy from reaching its full physical potential. One of the main contributing sources is the uncertainty in estimating particle stopping power (ρ_s) within patients. Currently, theρ_sdistribution in a patient is derived from a single-energy CT (SECT) scan acquired for treatment planning by converting CT number expressed in Hounsfield units (HU) of each voxel toρ_susing a Hounsfield look-up table (HLUT), also known as the CT calibration curve. HU andρ_sshare a linear relationship with electron density but differ in their additional dependence on elemental composition through different physical properties, i.e. effective atomic number and mean excitation energy, respectively. Because of that, the HLUT approach is particularly sensitive to differences in elemental composition between real human tissues and tissue surrogates as well as tissue variations within and among individual patients. The use of dual-energy CT (DECT) forρ_sprediction has been shown to be effective in reducing the uncertainty inρ_sestimation compared to SECT. The acquisition of CT data over different x-ray spectra yields additional information on the material elemental composition. Recently, multi-energy CT (MECT) has been explored to deduct material-specific information with higher dimensionality, which has the potential to further improve the accuracy ofρ_sestimation. Even though various DECT and MECT methods have been proposed and evaluated over the years, these approaches are still only scarcely implemented in routine clinical practice. In this topical review, we aim at accelerating this translation process by providing: (1) a comprehensive review of the existing DECT/MECT methods forρ_sestimation with their respective strengths and weaknesses; (2) a general review of uncertainties associated with DECT/MECT methods; (3) a general review of different aspects related to clinical implementation of DECT/MECT methods; (4) other potential advanced DECT/MECT applications beyondρ_sestimation.

Virtual computed-tomography system for deep-learning-based material decomposition

Objective. Material decomposition (MD) evaluates the elemental composition of human tissues and organs via computed tomography (CT) and is indispensable in correlating anatomical images with functional ones. A major issue in MD is inaccurate elemental information about the real human body. To overcome this problem, we developed a virtual CT system model, by which various reconstructed images can be generated based on ICRP110 human phantoms with information about six major elements (H, C, N, O, P, and Ca). Approach. We generated CT datasets labelled with accurate elemental information using the proposed generative CT model and trained a deep learning (DL)-based model to estimate the material distribution with the ICRP110 based human phantom as well as the digital Shepp–Logan phantom. The accuracy in quad-, dual-, and single-energy CT cases was investigated. The influence of beam-hardening artefacts, noise, and spectrum variations were analysed with testing datasets including elemental density and anatomical shape variations. Main results. The results indicated that this DL approach can realise precise MD, even with single-energy CT images. Moreover, noise, beam-hardening artefacts, and spectrum variations were shown to have minimal impact on the MD. Significance. Present results suggest that the difficulty to prepare a large CT database can be solved by introducing the virtual CT system and the proposed technique can be applied to clinical radiodiagnosis and radiotherapy.

Simultaneous Image Reconstruction and Element Decomposition for Iodine Contrast Agent Visualization in Multienergy Element-Resolved Cone Beam CT

Iodine contrast agent is widely used in liver cancer radiotherapy at CT simulation stage to enhance detectability of tumor. However, its application in cone beam CT (CBCT) for image guidance before treatment delivery is still limited because of poor image quality and excessive dose of contrast agent during multiple treatment fractions. We previously developed a multienergy element-resolved (MEER) CBCT framework that included x-ray projection data acquisition on a conventional CBCT platform in a kVp-switching model and a dictionary-based image reconstruction algorithm that simultaneously reconstructed x-ray attenuation images at each kilovoltage peak (kVp), an electron density image, and elemental composition images. In this study, we investigated feasibility using MEER-CBCT for low-concentration iodine contrast agent visualization. We performed simulation and experimental studies using a phantom with inserts containing water and different concentrations of iodine solution and the MEER-CBCT scan with 600 projections in a full gantry rotation, in which the kVp level sequentially changed among 80, 100, and 120 kVps. We included iodine material in the dictionary of the reconstruction algorithm. We analyzed iodine detectability as quantified by contrast-to-noise ratio (CNR) and compared results with those of CBCT images reconstructed by the standard filter back projection (FBP) method with 600 projections. MEER-CBCT achieved similar contrast enhancement as FBP method but significantly higher CNR. At 2.5% iodine solution concentration, FBP method achieved 170 HU enhancement and CNR of 2.0, considered the standard CNR for successful tumor visualization. MEER-CBCT achieved the same CNR but at ~6.3 times lower iodine concentration of 0.4%.

Improving dose calculation accuracy in preclinical radiation experiments using multi-energy element resolved cone beam CT

Cone-beam CT (CBCT) in modern pre-clinical small-animal radiation research platforms provides volumetric images for image guidance and experiment planning purposes. In this work, we implemented multi-energy element-resolved (MEER) CBCT using three scans with different kVps on a SmART platform (Precision X-ray Inc.) We performed comprehensive calibration tasks achieve sufficient accuracy for this quantitative imaging purpose. For geometry calibration, we scanned a ball bearing phantom and used an analytical method together with an optimization approach to derive gantry-angle specific geometry parameters. Intensity calibration and correction included the corrections for detector lag, glare, and beam hardening. The corrected CBCT projection images acquired at 30, 40 and 60 kVp in multiple scans were used to reconstruct CBCT images using the Feldkamp-Davis-Kress reconstruction algorithm. After that, an optimization problem was solved to determine images of relative electron density (rED) and elemental composition (EC) that are needed for Monte Carlo-based radiation dose calculation. We demonstrated effectiveness of our CBCT calibration steps by showing improvements in image quality and successful material decomposition in cases with a small animal CT calibration phantom and a plastinated mouse phantom. It was found that artifacts induced by geometry inaccuracy, detector lag, glare and beam hardening were visually reduced. CT number mean errors were reduced from 19\% to 5\%. In the CT calibration phantom case, median errors in H, O, and Ca fractions for all the inserts were below 1\%, 2\%, and 4\% respectively, and median error in rED was less than 5\%. Compared to standard approach deriving material type and rED via CT number conversion, our approach improved Monte Carlo simulation-based dose calculation accuracy in bone regions. Mean dose error was reduced from 47.5\% to 10.9\%.

A prior image constraint robust principal component analysis reconstruction method for sparse segmental multi-energy computed tomography

Background

Multi-energy computed tomography (MECT) is a promising technique in medical imaging, especially for quantitative imaging. However, high technical requirements and system costs barrier its step into clinical practice.

Methods

We propose a novel sparse segmental MECT (SSMECT) scheme and corresponding reconstruction method, which is a cost-efficient way to realize MECT on a conventional single-source CT. For the data acquisition, the X-ray source is controlled to maintain an energy within a segmental arc, and then switch alternately to another energy level. This scan only needs to switch tube voltage a few times to acquire multi-energy data, but leads to sparse-view and limited-angle issues in image reconstruction. To solve this problem, we propose a prior image constraint robust principal component analysis (PIC-RPCA) reconstruction method, which introduces structural similarity and spectral correlation into the reconstruction.

Results

A numerical simulation and a real phantom experiment were conducted to demonstrate the efficacy and robustness of the scan scheme and reconstruction method. The results showed that our proposed reconstruction method could have achieved better multi-energy images than other competing methods both quantitatively and qualitatively.

Conclusions

Our proposed SSMECT scan with PIC-RPCA reconstruction method could lower kVp switching frequency while achieving satisfactory reconstruction accuracy and image quality.

Multi-Material Decomposition for Single Energy CT Using Material Sparsity Constraint

Multi-material decomposition (MMD) decomposes CT images into basis material images, and is a promising technique in clinical diagnostic CT to identify material compositions within the human body. MMD could be implemented on measurements obtained from spectral CT protocol, although spectral CT data acquisition is not readily available in most clinical environments. MMD methods using single energy CT (SECT), broadly applied in radiological departments of most hospitals, have been proposed in the literature while challenged by the inferior decomposition accuracy and the limited number of material bases due to the constrained material information in the SECT measurement. In this paper, we propose an image-domain SECT MMD method using material sparsity as an assistance under the condition that each voxel of the CT image contains at most two different elemental materials. L₀ norm represents the material sparsity constraint (MSC) and is integrated into the decomposition objective function with a least-square data fidelity term, total variation term, and a sum-to-one constraint of material volume fractions. An accelerated primal-dual (APD) algorithm with line-search scheme is applied to solve the problem. The pixelwise direct inversion method with the two-material assumption (TMA) is applied to estimate the initials. We validate the proposed method on phantom and patient data. Compared with the TMA method, the proposed MSC method increases the volume fraction accuracy (VFA) from 92.0% to 98.5% in the phantom study. In the patient study, the calcification area can be clearly visualized in the virtual non-contrast image generated by the proposed method, and has a similar shape to that in the ground-truth contrast-free CT image. The high decomposition image quality from the proposed method substantially facilitates the SECT-based MMD clinical applications.

DECT-MULTRA: Dual-Energy CT Image Decomposition With Learned Mixed Material Models and Efficient Clustering

Dual-energy computed tomography (DECT) imaging plays an important role in advanced imaging applications due to its material decomposition capability. Image-domain decomposition operates directly on CT images using linear matrix inversion, but the decomposed material images can be severely degraded by noise and artifacts. This paper proposes a new method dubbed DECT-MULTRA for image-domain DECT material decomposition that combines conventional penalized weighted-least squares (PWLS) estimation with regularization based on a mixed union of learned transforms (MULTRA) model. Our proposed approach pre-learns a union of common-material sparsifying transforms from patches extracted from all the basis materials, and a union of cross-material sparsifying transforms from multi-material patches. The common-material transforms capture the common properties among different material images, while the cross-material transforms capture the cross-dependencies. The proposed PWLS formulation is optimized efficiently by alternating between an image update step and a sparse coding and clustering step, with both of these steps having closed-form solutions. The effectiveness of our method is validated with both XCAT phantom and clinical head data. The results demonstrate that our proposed method provides superior material image quality and decomposition accuracy compared to other competing methods.

Model-based material decomposition with a penalized nonlinear least-squares CT reconstruction algorithm

Spectral information in CT may be used for material decomposition to produce accurate reconstructions of material density and to separate materials with similar overall attenuation. Traditional methods separate the reconstruction and decomposition steps, often resulting in undesirable trade-offs (e.g. sampling constraints, a simplified spectral model). In this work, we present a model-based material decomposition algorithm which performs the reconstruction and decomposition simultaneously using a multienergy forward model. In a kV-switching simulation study, the presented method is capable of reconstructing iodine at 0.5 mg ml^-1 with a contrast-to-noise ratio greater than two, as compared to 3.0 mg ml^-1 for image domain decomposition. The presented method also enables novel acquisition methods, which was demonstrated in this work with a combined kV-switching/split-filter acquisition explored in simulation and physical test bench studies. This novel design used four spectral channels to decompose three materials: water, iodine, and gadolinium. In simulation, the presented method accurately reconstructed concentration value estimates with RMSE values of 4.86 mg ml^-1 for water, 0.108 mg ml^-1 for iodine and 0.170 mg ml^-1 for gadolinium. In test-bench data, the RMSE values were 134 mg ml^-1, 5.26 mg ml^-1 and 1.85 mg ml^-1, respectively. These studies demonstrate the ability of model-based material decomposition to produce accurate concentration estimates in challenging spatial/spectral sampling acquisitions.