Application of machine learning classifiers to X-ray diffraction imaging with medically relevant phantoms.

PURPOSE: Recent studies have demonstrated the ability to rapidly produce large field of view X-ray diffraction (XRD) images, which provide rich new data relevant to the understanding and analysis of disease. However, work has only just begun on developing algorithms that maximize the performance toward decision-making and diagnostic tasks. In this study, we present the implementation of and comparison between rules-based and machine learning (ML) classifiers on XRD images of medically relevant phantoms to explore the potential for increased classification performance. METHODS: Medically relevant phantoms were utilized to provide well-characterized ground-truths for comparing classifier performance. Water and polylactic acid (PLA) plastic were used as surrogates for cancerous and healthy tissue, respectively, and phantoms were created with varying levels of spatial complexity and biologically relevant features for quantitative testing of classifier performance. Our previously developed X-ray scanner was used to acquire co-registered X-ray transmission and diffraction images of the phantoms. For classification algorithms, we explored and compared two rules-based classifiers (cross-correlation, or matched-filter, and linear least-squares unmixing) and two ML classifiers (support vector machines and shallow neural networks). Reference XRD spectra (measured by a commercial diffractometer) were provided to the rules-based algorithms, while 60% of the measured XRD pixels were used for training of the ML algorithms. The area under the receiver operating characteristic curve (AUC) was used as a comparative metric between the classification algorithms, along with the accuracy performance at the midpoint threshold for each classifier. RESULTS: The AUC values for material classification were 0.994 (cross-correlation [CC]), 0.994 (least-squares [LS]), 0.995 (support vector machine [SVM]), and 0.999 (shallow neural network [SNN]). Setting the classification threshold to the midpoint for each classifier resulted in accuracy values of CC = 96.48%, LS = 96.48%, SVM = 97.36%, and SNN = 98.94%. If only considering pixels ±3 mm from water-PLA boundaries (where partial volume effects could occur due to imaging resolution limits), the classification accuracies were CC = 89.32%, LS = 89.32%, SVM = 92.03%, and SNN = 96.79%, demonstrating an even larger improvement produced by the machine-learned algorithms in spatial regions critical for imaging tasks. Classification by transmission data alone produced an AUC of 0.773 and accuracy of 85.45%, well below the performance levels of any of the classifiers applied to XRD image data. CONCLUSIONS: We demonstrated that ML-based classifiers outperformed rules-based approaches in terms of overall classification accuracy and improved the spatially resolved classification performance on XRD images of medical phantoms. In particular, the ML algorithms demonstrated considerably improved performance whenever multiple materials existed in a single voxel. The quantitative performance gains demonstrate an avenue to extract and harness XRD imaging data to improve material analysis for research, industrial, and clinical applications.

X-ray fan beam coded aperture transmission and diffraction imaging for fast material analysis.

X-ray transmission imaging has been used in a variety of applications for high-resolution measurements based on shape and density. Similarly, X-ray diffraction (XRD) imaging has been used widely for molecular structure-based identification of materials. Combining these X-ray methods has the potential to provide high-resolution material identification, exceeding the capabilities of either modality alone. However, XRD imaging methods have been limited in application by their long measurement times and poor spatial resolution, which has generally precluded combined, rapid measurements of X-ray transmission and diffraction. In this work, we present a novel X-ray fan beam coded aperture transmission and diffraction imaging system, developed using commercially available components, for rapid and accurate non-destructive imaging of industrial and biomedical specimens. The imaging system uses a 160 kV Bremsstrahlung X-ray source while achieving a spatial resolution of ≈ 1 × 1 mm2 and a spectral accuracy of > 95% with only 15 s exposures per 150 mm fan beam slice. Applications of this technology are reported in geological imaging, pharmaceutical inspection, and medical diagnosis. The performance of the imaging system indicates improved material differentiation relative to transmission imaging alone at scan times suitable for a variety of industrial and biomedical applications.

Simulation based evaluation of a fan beam coded aperture x-ray diffraction imaging system for biospecimen analysis.

X-ray diffraction (XRD) imaging yields spatially resolved, material-specific information, which can aid medical diagnosis and inform treatment. In this work we used simulations to analyze the utility of fan beam coded aperture XRD imaging for fast, high-resolution scatter imaging of biospecimens for tissue assessment. To evaluate the proposed system's utility in a specific task, we employed a deterministic model to produce simulated data from biologically realistic breast tissue phantoms and model-based reconstruction to recover a spatial map of the XRD signatures throughout the phantoms. We found an XRD spatial resolution of ≈1 mm with a mean reconstructed spectral accuracy of 0.98 ± 0.01 for a simulated 1 × 150 mm2 fan beam operating at 160 kVp, 10 mA, and 4.5 s exposures. A classifier for cancer detection was developed utilizing cross-correlation of XRD spectra against a spectral library, with a receiver operating characteristic curve with an area under the curve value of 0.972. Our results indicated a potential diagnostic modality that could aid in tasks ranging from analysis of ex-vivo pathology biospecimens to intraoperative cancer margin assessment, motivating future work to develop an experimental system while enabling the development of improved algorithms for imaging and tissue analysis-based classification performance.

Monte Carlo dosimetry modeling of focused kV x-ray radiotherapy of eye diseases with potential nanoparticle dose enhancement.

PURPOSE: Eye plaque brachytherapy is the most common approach for intraocular cancer treatment. It is, however, invasive and subject to large setup uncertainty due to the surgical operation. We propose a novel-focused kV x-ray technique with potential nanoparticle (NP) enhancement and evaluate its application in treating choroidal melanoma and iris melanoma by Monte Carlo (MC) dosimetry modeling. METHODS: A polycapillary x-ray lens was used to focus 45 kVp x rays to achieve pinpoint accuracy of dose delivery to small tumors near critical structures. In addition to allowing for beam focusing, the use of kV x rays takes advantage of the strong photoelectric absorption of metallic NPs in that energy regime and hence strong radiosensitization. We constructed an MC simulation program that takes into account the x-ray optic modeling and used GEANT4 for dosimetric calculation. Extensive phantom measurements using a prototype-focused x-ray system were carried out. The MC simulation of simple geometry phantom irradiation was first compared to measurements to verify the x-ray optic lens modeling in conjunction with the Geant4 dosimetric calculation. To simulate tumor treatment, a geometric eye model and two tumor models were constructed. Dose distributions of the simulated treatments were then calculated. NP radiosensitization was also simulated for two concentrations of 2 nm gold NP (AuNP) uniformly distributed in the tumor. RESULTS: The MC-simulated full width at half maximum (FWHM) and central-axis depth dose of the focused kV x-ray beam match those measured on EBT3 films within ~10% around the depth of focus of the beam. Dose distributions of the simulated ocular tumor treatments show that focused x-ray beams can concentrate the high-dose region in or close to the tumor plus margin. For the simulated posterior choroidal tumor treatment, with sufficient tumor coverage, the doses to the optic disc and fovea are substantially reduced with focused x-ray therapy compared to eye plaque treatment (3.8 vs 39.8 Gy and 11.1 vs 53.8 Gy, respectively). The eye plaque treatment was calculated using an Eye Physics plaque with I-125 seeds under TG43 assumption. For the energy spectrum used in this study, the average simulated dose enhancement ratios (DERs) are roughly 2.1 and 1.1 for 1.0% and 0.1% AuNP mass concentration in the tumor, respectively. CONCLUSION: Compared to eye plaque brachytherapy, the proposed focused kV x-ray technique is noninvasive and shows great advantage in sparing healthy critical organs without sacrificing the tumor control. The NP radiation dose enhancement is considerable at our proposed kV range even with a low NP concentration in the tumor, providing better critical structure protection and more flexibility for treatment planning.

On the use of bolus for pacemaker dose measurement and reduction in radiation therapy.

Special attention is required in planning and administering radiation therapy to patients with cardiac implantable electronic devices (CIEDs), such as pacemaker and defibrillator. The range of dose to CIEDs that can induce malfunction is large among CIEDs. Clinically significant defects have been reported at dose as low as 0.15 Gy. Therefore, accurate estimation of dose to CIED and dose reduction are both important even if the dose is expected to be less than the often-used 2-Gy limit. We investigated the use of bolus in in vivo dosimetry for CIEDs. Solid water phantom measurements of out-of-field dose for a 6-MV beam were performed using parallel plate chamber with and without 1- to 2-cm bolus covering the chamber. In vivo dosimetry at skin surface above the CIED was performed with and without bolus covering the CIED for three patients with the CIED <5 cm from the field edge. Chamber measured dose at depth ~0.5-1.5 cm below the skin surface, where the CIED is normally located, was reduced by ~7-48% with bolus. The dose reduction became smaller at deeper depths and with smaller field size. In vivo dosimetry at skin surface also indicated ~20%-60% lower dose when using bolus for the three patients. The dose measured with bolus more accurately reflects the dose to CIED and is less affected by contaminant electrons and linac head scatter. In general, the treatment planning system (TPS) calculation underestimated the dose to CIED, but it predicts the CIED dose more accurately when bolus is used. We recommend the use of 1- to 2-cm bolus to cover the CIED during in vivo CIED dose measurements for more accurate CIED dose estimation. If the CIED is placed <2 cm in depth and its dose is mainly from anterior beams, we recommend using the bolus during the entire course of radiation delivery to reduce the dose to CIED.