Machine learning based oxygen and carbon concentration derivation using dual-energy CT for PET-based dose verification in proton therapy

PURPOSE

Online dose verification based on proton-induced positron emitters requires high accuracy in the assignment of elemental composition (e.g. C and O). We developed a machine learning framework for deriving oxygen and carbon concentration based on dual-energy CT (DECT).

METHODS

Digital phantoms at the head site were constructed based on single-energy CT (SECT) and stoichiometric calibration. DECT images (80 kVp and 140 kVp) were synthesized using two methods: 1) theoretical CT numbers with Gaussian noise (method 1) and 2) forward/backward image reconstruction with poly-energetic energy spectrum and Poisson noise modeled (method 2). Two architectures of convolutional neural networks, UNet and ResNet, were investigated to map from DECT images to C/O weights. Four cases (UNet-1: Method 1+UNet, ResNet-1: Method 1+ResNet, UNet-2: Method 2+UNet, and ResNet-2: Method 2 +ResNet) were tested for different tissue types and different noise levels. Monte-Carlo simulation was employed to identify the impact of fluctuation in oxygen and carbon concentration on activity/dose distribution.

RESULTS

When no noise present, all four cases are able to obtain <2% mean absolute errors (MAE) and <4% root mean square error (RMSE). For the worst image quality (e.g. lowest image SNR), the RMSE for O among all tissue types is 3.02% (UNet-1), 4.46% (ResNet-1), 4.38% (UNet-2) and 6.31% (ResNet-2), respectively. For UNet-1 and ResNet-1, the model performed slightly better in terms of RMSE for skeletal tissue than soft tissues. Such trend is not observed for UNet-2 and ResNet-2. With regard to the comparison between UNet and ResNet, different accuracy and noise immunity are observed. The activity profiles exhibit 3-5% difference in terms of mean relative error (MRE) between the ground truth and machine learning outcome.

CONCLUSION

We explored the feasibility of a machine learning framework to derive elemental concentration of oxygen and carbon based on DECT images. Two machine learning models, UNet and ResNet, are able to utilize spatial correlation and obtain accurate carbon and oxygen concentration. This study lays a foundation for us to apply the proposed approach to clinical DECT images. This article is protected by copyright. All rights reserved.

Tungsten anode spectral model using interpolating cubic splines: unfiltered x-ray spectra from 20 kV to 640 kV

PURPOSE

Monte Carlo methods were used to generate lightly filtered high resolution x-ray spectra spanning from 20 kV to 640 kV.

METHODS

X-ray spectra were simulated for a conventional tungsten anode. The Monte Carlo N-Particle eXtended radiation transport code (MCNPX 2.6.0) was used to produce 35 spectra over the tube potential range from 20 kV to 640 kV, and cubic spline interpolation procedures were used to create piecewise polynomials characterizing the photon fluence per energy bin as a function of x-ray tube potential. Using these basis spectra and the cubic spline interpolation, 621 spectra were generated at 1 kV intervals from 20 to 640 kV. The tungsten anode spectral model using interpolating cubic splines (TASMICS) produces minimally filtered (0.8 mm Be) x-ray spectra with 1 keV energy resolution. The TASMICS spectra were compared mathematically with other, previously reported spectra.

RESULTS

Using pairedt-test analyses, no statistically significant difference (i.e., p > 0.05) was observed between compared spectra over energy bins above 1% of peak bremsstrahlung fluence. For all energy bins, the correlation of determination (R(2)) demonstrated good correlation for all spectral comparisons. The mean overall difference (MOD) and mean absolute difference (MAD) were computed over energy bins (above 1% of peak bremsstrahlung fluence) and over all the kV permutations compared. MOD and MAD comparisons with previously reported spectra were 2.7% and 9.7%, respectively (TASMIP), 0.1% and 12.0%, respectively [R. Birch and M. Marshall, "Computation of bremsstrahlung x-ray spectra and comparison with spectra measured with a Ge(Li) detector," Phys. Med. Biol. 24, 505-517 (1979)], 0.4% and 8.1%, respectively (Poludniowski), and 0.4% and 8.1%, respectively (AAPM TG 195). The effective energy of TASMICS spectra with 2.5 mm of added Al filtration ranged from 17 keV (at 20 kV) to 138 keV (at 640 kV); with 0.2 mm of added Cu filtration the effective energy was 9 keV at 20 kV and 169 keV at 640 kV.

CONCLUSIONS

Ranging from 20 kV to 640 kV, 621 x-ray spectra were produced and are available at 1 kV tube potential intervals. The spectra are tabulated at 1 keV intervals. TASMICS spectra were shown to be largely equivalent to published spectral models and are available in spreadsheet format for interested users by emailing the corresponding author (JMB).

Comparative study of a low-Z cone-beam computed tomography system

Computed tomography images have been acquired using an experimental (low atomic number (Z) insert) megavoltage cone-beam imaging system. These images have been compared with standard megavoltage and kilovoltage imaging systems. The experimental system requires a simple modification to the 4 MeV electron beam from an Elekta Precise linac. Low-energy photons are produced in the standard medium-Z electron window and a low-Z carbon electron absorber located after the window. The carbon electron absorber produces photons as well as ensuring that all remaining electrons from the source are removed. A detector sensitive to diagnostic x-ray energies is also employed. Quantitative assessment of cone-beam computed tomography (CBCT) contrast shows that the low-Z imaging system is an order of magnitude or more superior to a standard 6 MV imaging system. CBCT data with the same contrast-to-noise ratio as a kilovoltage imaging system (0.15 cGy) can be obtained in doses of 11 and 244 cGy for the experimental and standard 6 MV systems, respectively. Whilst these doses are high for everyday imaging, qualitative images indicate that kilovoltage like images suitable for patient positioning can be acquired in radiation doses of 1-8 cGy with the experimental low-Z system.