Arthur Lalonde

The potential of dual‐energy CT to reduce proton beam range uncertainties

Purpose Dual‐energy CT (DECT) promises improvements in estimating stopping power ratios (SPRs) for proton therapy treatment planning. Although several comparable mathematical formalisms have been proposed in literature, the optimal techniques to characterize human tissue SPRs with DECT in a clinical environment are not fully established. The aim of this work is to compare the most robust DECT methods against conventional single‐energy CT (SECT) in conditions reproducing a clinical environment, where CT artifacts and noise play a major role on the accuracy of these techniques. Methods Available DECT tissue characterization methods are investigated and their ability to predict SPRs is compared in three contexts: (a) a theoretical environment using the XCOM cross section database; (b) experimental data using a dual‐source CT scanner on a calibration phantom; (c) simulations of a virtual humanoid phantom with the ImaSim software. The latter comparison accounts for uncertainties caused by CT artifacts and noise, but leaves aside other sources of uncertainties such as CT grid size and the I‐values. To evaluate the clinical impact, a beam range calculation model is used to predict errors from the probability distribution functions determined with ImaSim simulations. Range errors caused by SPR errors in soft tissues and bones are investigated. Results Range error estimations demonstrate that DECT has the potential of reducing proton beam range uncertainties by 0.4% in soft tissues using low noise levels of 12 and 8 HU in DECT, corresponding to 7 HU in SECT. For range uncertainties caused by the transport of protons through bones, the reduction in range uncertainties for the same levels of noise is found to be up to 0.6 to 1.1 mm for bone thicknesses ranging from 1 to 5 cm, respectively. We also show that for double the amount noise, i.e., 14 HU in SECT and 24 and 16 HU for DECT, the advantages of DECT in soft tissues are lost over SECT. In bones however, the reduction in range uncertainties is found to be between 0.5 and 0.9 mm for bone thicknesses ranging from 1 to 5 cm, respectively. Conclusion DECT has a clear potential to improve proton beam range predictions over SECT in proton therapy. However, in the current state high levels of noise remain problematic for DECT characterization methods and do not allow getting the full benefits of this technology. Future work should focus on adapting DECT methods to noise and investigate methods based on raw‐data to reduce CT artifacts.

A general method to derive tissue parameters for Monte Carlo dose calculation with multi-energy CT

To develop a general method for human tissue characterization with dual- and multi-energy CT and evaluate its performance in determining elemental compositions and quantities relevant to radiotherapy Monte Carlo dose calculation. Ideal materials to describe human tissue are obtained applying principal component analysis on elemental weight and density data available in literature. The theory is adapted to elemental composition for solving tissue information from CT data. A novel stoichiometric calibration method is integrated to the technique to make it suitable for a clinical environment. The performance of the method is compared with two techniques known in literature using theoretical CT data. In determining elemental weights with dual-energy CT, the method is shown to be systematically superior to the water-lipid-protein material decomposition and comparable to the parameterization technique. In determining proton stopping powers and energy absorption coefficients with dual-energy CT, the method generally shows better accuracy and unbiased results. The generality of the method is demonstrated simulating multi-energy CT data to show the potential to extract more information with multiple energies. The method proposed in this paper shows good performance to determine elemental compositions from dual-energy CT data and physical quantities relevant to radiotherapy dose calculation. The method is particularly suitable for Monte Carlo calculations and shows promise in using more than two energies to characterize human tissue with CT.

Robust quantitative contrast‐enhanced dual‐energy CT for radiotherapy applications

Purpose The purpose of this study was to develop and validate accurate methods for determining iodine content and virtual noncontrast maps of physical parameters, such as electron density, in the context of radiotherapy. Methods A simulation environment is developed to compare three methods allowing extracting iodine content and virtual noncontrast composition: (a) two‐material decomposition, (b) three‐material decomposition with the conservation of volume constraint, and (c) eigentissue decomposition. The simulation allows comparing the performance of the methods using iodine‐based contrast agent contents in tissues from a reference dataset with variable density and elemental composition. The comparison is performed in two ways: (a) with a priori knowledge on the composition of the targeted tissue, and (b) without a priori knowledge on the base tissue. The three methods are tested with patient images scanned with dual‐energy CT and iodine‐based contrast agent. An experimental calibration adapted to the presence of iodine is performed by imaging tissue equivalent materials and diluted contrast agent solutions with known atomic composition. Results Results show that in the case of known a priori on the composition of the targeted tissue, the two‐material decomposition is robust to variable densities and atomic compositions without biasing the results. In the absence of a priori knowledge on the target tissue composition, the eigentissue decomposition method yields minimal bias and higher robustness to variations. Results from the experimental calibration and the images of two patients show that the extracted quantities are accurate and the bias is negligible for both methods with respect to clinical applications in their respective scope of use. For the patient imaged with a contrast agent, virtual noncontrast electron densities are found in good agreement with values extracted from the scan without contrast agent. Conclusion This study identifies two accurate methods to quantify iodine‐based contrast agents and virtual noncontrast composition images with dual‐energy CT. One is the two‐material decomposition with a priori knowledge of the constituent components focused on organ‐specific applications, such as kidney or lung function assessment. The other method is the eigentissue decomposition and is useful for general radiotherapy applications, such as treatment planning where accurate dose calculations are needed in the absence of contrast agent.

A Bayesian approach to solve proton stopping powers from noisy multi‐energy CT data

Purpose: To propose a new formalism allowing the characterization of human tissues from multienergy computed tomography (MECT) data affected by noise and to evaluate its performance in estimating proton stopping powers (SPR). Methods: A recently published formalism based on principal component analysis called eigentissue decomposition (ETD) is adapted to the context of noise using a Bayesian estimator. The method, named Bayesian ETD, uses the maximum a posteriori fractions of eigentissues in each voxel to determine physical parameters relevant for proton beam dose calculation. Simulated dual‐energy computed tomography (DECT) data are used to evaluate the performance of the proposed method to estimate SPR and to compare it to the initially proposed maximum‐likelihood ETD and to a state‐of‐the‐art ρe − Z formalism. To test the robustness of each method towards clinical reality, three different levels of noise are implemented, as well as variations in elemental composition and density of reference tissues. The impact of using more than two energy bins to determine SPR is also investigated by simulating MECT data using two to five energy bins. Finally, the impact of using MECT over DECT for range prediction is evaluated using a probabilistic model. Results: For simulated DECT data of reference tissues, the Bayesian ETD approach systematically gives lower root‐mean‐square (RMS) errors with negligible bias. For a medium level of noise, the RMS errors on SPR are found to be 2.78%, 2.76% and 1.53% for ρe − Z, maximum‐likelihood ETD, and Bayesian ETD, respectively. When variations are introduced to the elemental composition and density, all implemented methods give similar performances at low noise. However, for a medium noise level, the proposed Bayesian method outperforms the two others with a RMS error of 1.94%, compared to 2.79% and 2.78% for ρe − Z and maximum‐likelihood ETD, respectively. When more than two energy spectra are used, the Bayesian ETD is able to reduce RMS error on SPR using up to five energy bins. In terms of range prediction, Bayesian ETD with four energy bins in realistic conditions reduces proton beam range uncertainties by a factor of up to 1.5 compared to ρe − Z. Conclusion: The Bayesian ETD is shown to be more robust against noise than similar methods and a promising approach to extract SPR from noisy DECT data. In the advent of commercially available multi‐energy CT or photon‐counting CT scanners, the Bayesian ETD is expected to allow extracting more information and improve the precision of proton therapy beyond DECT.

The impact of dual- and multi-energy CT on proton pencil beam range uncertainties: a Monte Carlo study

The purpose of this work is to evaluate the impact of single-, dual- and multi-energy CT (SECT, DECT and MECT) on proton range uncertainties in a patient like geometry and a full Monte Carlo environment. A virtual patient is generated from a real patient pelvis CT scan, where known mass densities and elemental compositions are overwritten in each voxel. Simulated CT images for SECT, DECT and MECT are generated for two limiting cases: (1) theoretical and idealistic CT numbers only affected by Gaussian noise (case A, the best scenario) and (2) reconstructed polyenergetic sinograms containing beam hardening, projection-based Poisson noise, and reconstruction artifacts (case B, the worst scenario). Conversion of the simulated SECT images into Monte Carlo inputs is done following the stoichiometric calibration method. For DECT and MECT, the Bayesian eigentissue decomposition method of Lalonde (2017 Med. Phys. 44 5293-302) is used. Pencil beams from seven different angles around the virtual patient are simulated using TOPAS to assess the performance of each method. Percentage depth doses curves (PDD) are compared to ground truth in order to determine the accuracy of range prediction of each imaging modality. For the idealistic images of case A, MECT and DECT slightly outperforms SECT. Root mean square (RMS) errors or 0.78 mm, 0.49 mm and 0.42 mm on R 80 mm, are observed for SECT, DECT and MECT respectively. In case B, PDD calculated in the MECT derived Monte Carlo inputs generally shows the best agreement with ground truth in both shape and position, with RMS errors of 2.03 mm, 1.38 mm and 0.86 mm for SECT, DECT and MECT respectively. Overall, the Bayesian eigentissue decomposition used with DECT systematically predicts proton ranges more accurately than the gold standard SECT-based approach. When CT numbers are severely affected by imaging artifacts, MECT with four energy bins becomes more reliable than both DECT and SECT.

Unsupervised classification of tissues composition for Monte Carlo dose calculation

The purpose of this study is to investigate the potential of k-means clustering to efficiently reduce the variety of materials needed in Monte Carlo (MC) dose calculation. A numerical phantom with 31 human tissues surrounded by water is created. K-means clustering is used to group the tissues in clusters of constant elemental composition. Four different distance measures are used to perform the clustering technique: Euclidean, Standardized Euclidean, Chi-Squared and Cityblock. Dose distributions are calculated with MC simulations for both low-kV photons and MeV protons using the clustered and reference elemental composition. Comparison between the dose distributions in the clustered and non-clustered phantom are made to assess the impact of clustering with each distance measure. The statistical significance of the differences observed between the four different metrics is determined by comparing the accuracy of energy absorption coefficients (EAC) of low-kV photons and proton stopping powers relative to water (SPR) for repeated clustering procedures. The performance of the proposed approach for a larger number of original materials is evaluated similarly by using a population of 62 000 statistically generated materials grouped into classes defined with supervised and unsupervised classification. In the phantom geometry, the Chi-Squared distance is the one introducing the smallest error on dose distribution and significant differences are observed between the EAC and SPR values predicted by each distance metric. The proposed approach is also shown to be equivalent to a state-of-the-art supervised classification method for proton therapy, but beneficial for low-kV photons applications. In conclusion, k-means clustering successfully reduces the variety of materials needed for accurate MC dose calculation. Based on the performance of four distance measures, we conclude that k-means clustering using the Chi-Squared distance introduces the smallest errors on dose distribution. The method is shown to yield similar or improved accuracy on key physical parameters compared to supervised classification.

Dosimetric impact of dual-energy CT tissue segmentation for low-energy prostate brachytherapy: a Monte Carlo study

The purpose of this study is to evaluate the impact of a novel tissue characterization method using dual-energy over single-energy computed tomography (DECT and SECT) on Monte Carlo (MC) dose calculations for low-dose rate (LDR) prostate brachytherapy performed in a patient like geometry. A virtual patient geometry is created using contours from a real patient pelvis CT scan, where known elemental compositions and varying densities are overwritten in each voxel. A second phantom is made with additional calcifications. Both phantoms are the ground truth with which all results are compared. Simulated CT images are generated from them using attenuation coefficients taken from the XCOM database with a 100 kVp spectrum for SECT and 80 and 140Sn kVp for DECT. Tissue segmentation for Monte Carlo dose calculation is made using a stoichiometric calibration method for the simulated SECT images. For the DECT images, Bayesian eigentissue decomposition is used. A LDR prostate brachytherapy plan is defined with 125I sources and then calculated using the EGSnrc user-code Brachydose for each case. Dose distributions and dose-volume histograms (DVH) are compared to ground truth to assess the accuracy of tissue segmentation. For noiseless images, DECT-based tissue segmentation outperforms the SECT procedure with a root mean square error (RMS) on relative errors on dose distributions respectively of 2.39% versus 7.77%, and provides DVHs closest to the reference DVHs for all tissues. For a medium level of CT noise, Bayesian eigentissue decomposition still performs better on the overall dose calculation as the RMS error is found to be of 7.83% compared to 9.15% for SECT. Both methods give a similar DVH for the prostate while the DECT segmentation remains more accurate for organs at risk and in presence of calcifications, with less than 5% of RMS errors within the calcifications versus up to 154% for SECT. In a patient-like geometry, DECT-based tissue segmentation provides dose distributions with the highest accuracy and the least bias compared to SECT. When imaging noise is considered, benefits of DECT are noticeable if important calcifications are found within the prostate.

Optimized I-values for use with the Bragg additivity rule and their impact on proton stopping power and range uncertainty

Novel imaging modalities can improve the estimation of patient elemental compositions for particle treatment planning. The mean excitation energy (I-value) is a main contributor to the proton range uncertainty. To minimize their impact on beam range errors and quantify their uncertainties, the currently used I-values proposed in 1982 are revisited. The study aims at proposing a new set of optimized elemental I-values for use with the Bragg additivity rule (BAR) and establishing uncertainties on the optimized I-values and the BAR. We optimize elemental I-values for the use in compounds based on measured material I-values. We gain a new set of elemental I-values and corresponding uncertainties, based on the experimental uncertainties and our uncertainty model. We evaluate uncertainties on I-values and relative stopping powers (RSP) of 70 human tissues, taking into account statistical correlations between tissues and water. The effect of new I-values on proton beam ranges is quantified using Monte Carlo simulations. Our elemental I-values describe measured material I-values with higher accuracy than ICRU-recommended I-values (RMSE: 6.17% (ICRU), 5.19% (this work)). Our uncertainty model estimates an uncertainty component from the BAR to 4.42%. Using our elemental I-values, we calculate the I-value of water as 78.73  ±  2.89 eV, being consistent with ICRU 90 (78  ±  2 eV). We observe uncertainties on tissue I-values between 1.82-3.38 eV, and RSP uncertainties between 0.002%-0.44%. With transport simulations of a proton beam in human tissues, we observe range uncertainties between 0.31% and 0.47%, as compared to current estimates of 1.5%. We propose a set of elemental I-values well suited for human tissues in combination with the BAR. Our model establishes uncertainties on elemental I-values and the BAR, enabling to quantify uncertainties on tissue I-values, RSP as well as particle range. This work is particularly relevant for Monte Carlo simulations where the interaction probabilities are reconstructed from elemental compositions.

TU-FG-BRB-02: The Impact of Using Dual-Energy CT for Determining Proton Stopping Powers: Comparison Between Theory and Experiments

PURPOSE To evaluate the clinical performance of dual-energy CT (DECT) in determining proton stopping power ratios (SPR) and demonstrate advantages over conventional single-energy CT (SECT). METHODS SECT and DECT scans of tissue-equivalent plastics as well as animal meat samples are performed with a Siemens SOMATOM Definition Flash. The methods of Schneider et al. (1996) and Bourque et al. (2014) are used to determine proton SPR on SECT and DECT images, respectively. Waterequivalent path length (WEPL) measurements of plastics and tissue samples are performed with a 195 MeV proton beam. WEPL values are determined experimentally using the depth-dose shift and dose extinction methods. RESULTS Comparison between CT-based and experimental WEPL is performed for 12 tissue-equivalent plastic as well as 6 meat boxes containing animal liver, kidney, heart, stomach, muscle and bones. For plastic materials, results show a systematic improvement in determining SPR with DECT, with a mean absolute error of 0.4% compared to 1.7% for SECT. For the meat samples, preliminary results show the ability for DECT to determine WEPL with a mean absolute value of 1.1% over all meat boxes. CONCLUSION This work demonstrates the potential in using DECT for determining proton SPR with plastic materials in a clinical context. Further work is required to show the benefits of DECT for tissue samples. While experimental uncertainties could be a limiting factor to show the benefits of DECT over SECT for the meat samples, further work is required to adapt the DECT formalism in the context of clinical use, where noise and artifacts play an important role.

TU-AB-BRC-03: Accurate Tissue Characterization for Monte Carlo Dose Calculation Using Dual-and Multi-Energy CT Data

PURPOSE To develop a general method for human tissue characterization with dual-and multi-energy CT and evaluate its performance in determining elemental compositions and the associated proton stopping power relative to water (SPR) and photon mass absorption coefficients (EAC). METHODS Principal component analysis is used to extract an optimal basis of virtual materials from a reference dataset of tissues. These principal components (PC) are used to perform two-material decomposition using simulated DECT data. The elemental mass fraction and the electron density in each tissue is retrieved by measuring the fraction of each PC. A stoichiometric calibration method is adapted to the technique to make it suitable for clinical use. The present approach is compared with two others: parametrization and three-material decomposition using the water-lipid-protein (WLP) triplet. RESULTS Monte Carlo simulations using TOPAS for four reference tissues shows that characterizing them with only two PC is enough to get a submillimetric precision on proton range prediction. Based on the simulated DECT data of 43 references tissues, the proposed method is in agreement with theoretical values of protons SPR and low-kV EAC with a RMS error of 0.11% and 0.35%, respectively. In comparison, parametrization and WLP respectively yield RMS errors of 0.13% and 0.29% on SPR, and 2.72% and 2.19% on EAC. Furthermore, the proposed approach shows potential applications for spectral CT. Using five PC and five energy bins reduces the SPR RMS error to 0.03%. CONCLUSION The proposed method shows good performance in determining elemental compositions from DECT data and physical quantities relevant to radiotherapy dose calculation and generally shows better accuracy and unbiased results compared to reference methods. The proposed method is particularly suitable for Monte Carlo calculations and shows promise in using more than two energies to characterize human tissue with CT.