A deep learning approach to radiation dose estimation.

Currently methods for predicting absorbed dose after administering a radiopharmaceutical are rather crude in daily clinical practice. Most importantly, individual tissue density distributions as well as local variations of the concentration of the radiopharmaceutical are commonly neglected. The current study proposes machine learning techniques like Green's function-based empirical mode decomposition and deep learning methods on U-net architectures in conjunction with soft tissue kernel Monte Carlo (MC) simulations to overcome current limitations in precision and reliability of dose estimations for clinical dosimetric applications. We present a hybrid method (DNN-EMD) based on deep neural networks (DNN) in combination with empirical mode decomposition (EMD) techniques. The algorithm receives x-ray computed tomography (CT) tissue density maps and dose maps, estimated according to the MIRD protocol, i.e. employing whole organ S-values and related time-integrated activities (TIAs), and from measured SPECT distributions of 177Lu radionuclei, and learns to predict individual absorbed dose distributions. In a second step, density maps are replaced by their intrinsic modes as deduced from an EMD analysis. The system is trained using individual full MC simulation results as reference. Data from a patient cohort of 26 subjects are reported in this study. The proposed methods were validated employing a leave-one-out cross-validation technique. Deviations of estimated dose from corresponding MC results corroborate a superior performance of the newly proposed hybrid DNN-EMD method compared to its related MIRD DVK dose calculation. Not only are the mean deviations much smaller with the new method, but also the related variances are much reduced. If intrinsic modes of the tissue density maps are input to the algorithm, variances become even further reduced though the mean deviations are less affected. The newly proposed hybrid DNN-EMD method for individualized radiation dose prediction outperforms the MIRD DVK dose calculation method. It is fast enough to be of use in daily clinical practice.

Factors affecting accuracy of S values and determination of time-integrated activity in clinical Lu-177 dosimetry.

INTRODUCTION: In any radiotherapy, the absorbed dose needs to be estimated based on two factors, the time-integrated activity of the administered radiopharmaceutical and the patient-specific dose kernel. In this study, we consider the uncertainty with which such absorbed dose estimation can be achieved in a clinical environment. METHODS: To calculate the total error of dose estimation we considered the following aspects: The error resulting from computing the time-integrated activity, the difference between the S-value and the patient specific full Monte Carlo simulation, the error from segmenting the volume-of-interest (kidney) and the intrinsic error of the activimeter. RESULTS: The total relative error in dose estimation can amount to 25.0% and is composed of the error of the time-integrated activity 17.1%, the error of the S-value 16.7%, the segmentation error 5.4% and the activimeter accuracy 5.0%. CONCLUSION: Errors from estimating the time-integrated activity and approximations applied to dose kernel computations contribute about equally and represent the dominant contributions far exceeding the contributions from VOI segmentation and activimeter accuracy.

A tool to automatically analyze electromagnetic tracking data from high dose rate brachytherapy of breast cancer patients.

During High Dose Rate Brachytherapy (HDR-BT) the spatial position of the radiation source inside catheters implanted into a female breast is determined via electromagnetic tracking (EMT). Dwell positions and dwell times of the radiation source are established, relative to the patient's anatomy, from an initial X-ray-CT-image. During the irradiation treatment, catheter displacements can occur due to patient movements. The current study develops an automatic analysis tool of EMT data sets recorded with a solenoid sensor to assure concordance of the source movement with the treatment plan. The tool combines machine learning techniques such as multi-dimensional scaling (MDS), ensemble empirical mode decomposition (EEMD), singular spectrum analysis (SSA) and particle filter (PF) to precisely detect and quantify any mismatch between the treatment plan and actual EMT measurements. We demonstrate that movement artifacts as well as technical signal distortions can be removed automatically and reliably, resulting in artifact-free reconstructed signals. This is a prerequisite for a highly accurate determination of any deviations of dwell positions from the treatment plan.

On the use of particle filters for electromagnetic tracking in high dose rate brachytherapy.

Modern radiotherapy of female breast cancers often employs high dose rate brachytherapy, where a radioactive source is moved inside catheters, implanted in the female breast, according to a prescribed treatment plan. Source localization relative to the patient's anatomy is determined with solenoid sensors whose spatial positions are measured with an electromagnetic tracking system. Precise sensor dwell position determination is of utmost importance to assure irradiation of the cancerous tissue according to the treatment plan. We present a hybrid data analysis system which combines multi-dimensional scaling with particle filters to precisely determine sensor dwell positions in the catheters during subsequent radiation treatment sessions. Both techniques are complemented with empirical mode decomposition for the removal of superimposed breathing artifacts. We show that the hybrid model robustly and reliably determines the spatial positions of all catheters used during the treatment and precisely determines any deviations of actual sensor dwell positions from the treatment plan. The hybrid system only relies on sensor positions measured with an EMT system and relates them to the spatial positions of the implanted catheters as initially determined with a computed x-ray tomography.

On the use of multi-dimensional scaling and electromagnetic tracking in high dose rate brachytherapy.

High dose rate brachytherapy affords a frequent reassurance of the precise dwell positions of the radiation source. The current investigation proposes a multi-dimensional scaling transformation of both data sets to estimate dwell positions without any external reference. Furthermore, the related distributions of dwell positions are characterized by uni-or bi-modal heavy-tailed distributions. The latter are well represented by α-stable distributions. The newly proposed data analysis provides dwell position deviations with high accuracy, and, furthermore, offers a convenient visualization of the actual shapes of the catheters which guide the radiation source during the treatment.