The future of MRI in radiation therapy: Challenges and opportunities for the MR community.

Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field.

A qualitative investigation of resilience and well-being among medical physics residents.

PURPOSE: Medical physics residents (MPRs) will define and shape the future of physics in medicine. We sought to better understand the residency experience, as related to resilience and well-being, through the lens of current MPRs and medical physicists (MPs) working with residents. METHODS AND MATERIALS: From February-May 2019, we conducted 32, 1-h, confidential, semi-structured interviews with MPs either currently enrolled in an accredited residency (n = 16) or currently employed by a department with an accredited residency (n = 16). Interviews centered on the topics of mentorship, work/life integration, and discrimination. Qualitative analysis methods were used to derive key themes from the interview transcripts. RESULTS: With regard to the medical physics residency experience, four key themes emerged during qualitative analysis: the demanding nature of medical physics residencies, the negative impacts of residency on MPRs during training and beyond, strategies MPRs use to cope with residency stress, and the role of professional societies in addressing residency-related change. CONCLUSIONS: Residency training is a stress-inducing time in the path to becoming a board-certified MP. By uncovering several sources of this stress, we have identified opportunities to support the resiliency and well-being of MPs in training through recommendations by professional societies, programmatic changes, and interventions at the department and residency program director level for residency programs, as well as strategies that MPRs themselves can use to support well-being on their career journey.

Clinical implementation of MR-guided vaginal cylinder brachytherapy.

We present an institutional experience on the clinical implementation of magnetic resonance (MR)-guided vaginal brachytherapy using commercially available solid applicator models. To test the fidelity of solid applicator models to digitize vaginal cylinder applicators, three datasets were evaluated. The first included 15 patients who were simulated with CT alone. Next, a water phantom was used to evaluate vaginal cylinders ranging from 20 to 35 mm in diameter with CT and MR. Finally, three patients undergoing HDR brachytherapy with vaginal cylinders that were simulated with both CT and MR were evaluated. In these assessments, the solid applicator models were aligned based on the outline of the applicators on the corresponding volumetric image, and deviations between the central source positions defined based on X-ray markers (on CT) and solid applicator models (on CT and MR), and the percent dose difference between select reference points were calculated. The mean central source deviation defined based on X-ray markers (on CT) and solid applicator models (on CT and MR) for the 15-patient cohort, the phantom, and the 3-patient cohort is 0.6 mm, 0.6 mm, and 1.2 mm, respectively. The average absolute percent dose difference for the bladder, rectum, prescription, and inferior reference points were 2.2%, 2.3%, 2.2%, and 2.4%, respectively, for the 15 patient cohort. For the phantom study, the average, absolute percent dose difference for the prescription and inferior reference points are 2.0% and 2.1% for the CT, 2.3% and 2.2% for the T1W, and 2.8% and 3.0% for the T2W images. For the three patient cohort, the average absolute percent dose difference for the bladder, rectum, prescription, and inferior reference points are 2.9%, 2.6%, 3.0%, and 4.2% for the CT, 6.5%, 1.6%, 2.5%, and 4.7% for the T1W, and 6.0%, 7.4%, 2.6, and 2.0% for the T2W images. Based on the current study, aligning the applicator model to MR images provides a practical, efficient approach to perform MR-based brachytherapy planning.

Volumetric MRI with sparse sampling for MR-guided 3D motion tracking via sparse prior-augmented implicit neural representation learning.

BACKGROUND: Volumetric reconstruction of magnetic resonance imaging (MRI) from sparse samples is desirable for 3D motion tracking and promises to improve magnetic resonance (MR)-guided radiation treatment precision. Data-driven sparse MRI reconstruction, however, requires large-scale training datasets for prior learning, which is time-consuming and challenging to acquire in clinical settings. PURPOSE: To investigate volumetric reconstruction of MRI from sparse samples of two orthogonal slices aided by sparse priors of two static 3D MRI through implicit neural representation (NeRP) learning, in support of 3D motion tracking during MR-guided radiotherapy. METHODS: A multi-layer perceptron network was trained to parameterize the NeRP model of a patient-specific MRI dataset, where the network takes 4D data coordinates of voxel locations and motion states as inputs and outputs corresponding voxel intensities. By first training the network to learn the NeRP of two static 3D MRI with different breathing motion states, prior information of patient breathing motion was embedded into network weights through optimization. The prior information was then augmented from two motion states to 31 motion states by querying the optimized network at interpolated and extrapolated motion state coordinates. Starting from the prior-augmented NeRP model as an initialization point, we further trained the network to fit sparse samples of two orthogonal MRI slices and the final volumetric reconstruction was obtained by querying the trained network at 3D spatial locations. We evaluated the proposed method using 5-min volumetric MRI time series with 340 ms temporal resolution for seven abdominal patients with hepatocellular carcinoma, acquired using golden-angle radial MRI sequence and reconstructed through retrospective sorting. Two volumetric MRI with inhale and exhale states respectively were selected from the first 30 s of the time series for prior embedding and augmentation. The remaining 4.5-min time series was used for volumetric reconstruction evaluation, where we retrospectively subsampled each MRI to two orthogonal slices and compared model-reconstructed images to ground truth images in terms of image quality and the capability of supporting 3D target motion tracking. RESULTS: Across the seven patients evaluated, the peak signal-to-noise-ratio between model-reconstructed and ground truth MR images was 38.02 ± 2.60 dB and the structure similarity index measure was 0.98 ± 0.01. Throughout the 4.5-min time period, gross tumor volume (GTV) motion estimated by deforming a reference state MRI to model-reconstructed and ground truth MRI showed good consistency. The 95-percentile Hausdorff distance between GTV contours was 2.41 ± 0.77 mm, which is less than the voxel dimension. The mean GTV centroid position difference between ground truth and model estimation was less than 1 mm in all three orthogonal directions. CONCLUSION: A prior-augmented NeRP model has been developed to reconstruct volumetric MRI from sparse samples of orthogonal cine slices. Only one exhale and one inhale 3D MRI were needed to train the model to learn prior information of patient breathing motion for sparse image reconstruction. The proposed model has the potential of supporting 3D motion tracking during MR-guided radiotherapy for improved treatment precision and promises a major simplification of the workflow by eliminating the need for large-scale training datasets.

Does Vascular Collapse Occur After Treatment of Hepatocellular Cancer With Stereotactic Body Radiation Therapy?

There is debate about why stereotactic body radiation therapy (SBRT) produces superior control of hepatocellular cancer (HCC) compared to fractionated treatment. Both preclinical and clinical evidence has been presented to support a "classic" biological explanation: the greater BED of SBRT produces more DNA damage and tumor cell kill. More recently, preclinical evidence has supported the concept of a "new biology", particularly radiation-induced vascular collapse, which increases hypoxia and free radical activation. This is hypothesized to cause much greater tumor cell death than was produced by the initial radiation-induced DNA damage to the tumor. We decided to investigate if vascular collapse occurs after standard SBRT for patients with HCC. Eight patients with 10 lesions underwent dynamic contrast enhanced MRI at the time of simulation and either 48 or 96 hours after the first fraction. Only three of 10 tumors showed a decrease in blood flow. These findings suggest that vascular collapse does not typically occur after SBRT for HCC.

Gender-Based Discrimination and Sexual Harassment in Medical Physics.

PURPOSE: Gender-based discrimination and sexual harassment have been well-studied in the fields of science, technology, engineering, math, and medicine. However, less is known about these topics and their effect within the profession of medical physics. We aimed to better understand and clarify the views and experiences of practicing medical physicists and medical physics residents regarding gender-based discrimination and sexual harassment. METHODS AND MATERIALS: We conducted in-depth, semistructured, and confidential interviews with 32 practicing medical physicists and medical physics residents across the United States. The interviews were broad and covered the topics of discrimination, mentorship, and work/life integration. All participants were associated with a department with a residency program accredited by the Commission on Accreditation of Medical Physics Education Programs and had appointments with a clinical component. RESULTS: Participants shared views about gender-based discrimination and sexual harassment that were polarized. Some perceived that discrimination and harassment were a current concern within medical physics, while some either perceived that they were not a concern or that discrimination positively affected women and minoritized populations. Many participants shared personal experiences of discrimination and harassment, including those related to unequal compensation, discrimination against mothers, discrimination during the hiring process, gender-biased assumptions about behaviors or goals, communication biases, and overt and persistent sexual harassment. CONCLUSIONS: There is an urgent need to acknowledge, better understand, and address gender-based discrimination and sexual harassment in the field of medical physics.

Dosimetric measurements of Onyx embolization material for stereotactic radiosurgery.

PURPOSE: Arteriovenous malformations are often treated with a combination of embolization and stereotactic radiosurgery. Concern has been expressed in the past regarding the dosimetric properties of materials used in embolization and the effects that the introduction of these materials into the brain may have on the quality of the radiosurgery plan. To quantify these effects, the authors have taken large volumes of Onyx 34 and Onyx 18 (ethylene-vinyl alcohol copolymer doped with tantalum) and measured the attenuation and interface effects of these embolization materials. METHODS: The manufacturer provided large cured volumes (∼28 cc) of both Onyx materials. These samples were 8.5 cm in diameter with a nominal thickness of 5 mm. The samples were placed on a block tray above a stack of solid water with an Attix chamber at a depth of 5 cm within the stack. The Attix chamber was used to measure the attenuation. These measurements were made for both 6 and 16 MV beams. Placing the sample directly on the solid water stack and varying the thickness of solid water between the sample and the Attix chamber measured the interface effects. The computed tomography (CT) numbers for bulk material were measured in a phantom using a wide bore CT scanner. RESULTS: The transmission through the Onyx materials relative to solid water was approximately 98% and 97% for 16 and 6 MV beams, respectively. The interface effect shows an enhancement of approximately 2% and 1% downstream for 16 and 6 MV beams. CT numbers of approximately 2600-3000 were measured for both materials, which corresponded to an apparent relative electron density (RED) ρ(e) (w) to water of approximately 2.7-2.9 if calculated from the commissioning data of the CT scanner. CONCLUSIONS: We performed direct measurements of attenuation and interface effects of Onyx 34 and Onyx 18 embolization materials with large samples. The introduction of embolization materials affects the dose distribution of a MV therapeutic beam, but should be of negligible consequence for effective thicknesses of less than 8 mm. The measured interface effects are also small, particularly at 6 MV. Large areas of high-density artifacts and low-density artifacts can cause errors in dose calculations and need to be identified and resolved during planning.

Real time volumetric MRI for 3D motion tracking via geometry-informed deep learning.

PURPOSE: To develop a geometry-informed deep learning framework for volumetric MRI with sub-second acquisition time in support of 3D motion tracking, which is highly desirable for improved radiotherapy precision but hindered by the long image acquisition time. METHODS: A 2D-3D deep learning network with an explicitly defined geometry module that embeds geometric priors of the k-space encoding pattern was investigated, where a 2D generation network first augmented the sparsely sampled image dataset by generating new 2D representations of the underlying 3D subject. A geometry module then unfolded the 2D representations to the volumetric space. Finally, a 3D refinement network took the unfolded 3D data and outputted high-resolution volumetric images. Patient-specific models were trained for seven abdominal patients to reconstruct volumetric MRI from both orthogonal cine slices and sparse radial samples. To evaluate the robustness of the proposed method to longitudinal patient anatomy and position changes, we tested the trained model on separate datasets acquired more than one month later and evaluated 3D target motion tracking accuracy using the model-reconstructed images by deforming a reference MRI with gross tumor volume (GTV) contours to a 5-min time series of both ground truth and model-reconstructed volumetric images with a temporal resolution of 340 ms. RESULTS: Across the seven patients evaluated, the median distances between model-predicted and ground truth GTV centroids in the superior-inferior direction were 0.4 ± 0.3 mm and 0.5 ± 0.4 mm for cine and radial acquisitions, respectively. The 95-percentile Hausdorff distances between model-predicted and ground truth GTV contours were 4.7 ± 1.1 mm and 3.2 ± 1.5 mm for cine and radial acquisitions, which are of the same scale as cross-plane image resolution. CONCLUSION: Incorporating geometric priors into deep learning model enables volumetric imaging with high spatial and temporal resolution, which is particularly valuable for 3D motion tracking and has the potential of greatly improving MRI-guided radiotherapy precision.

Real-time prediction of stomach motions based upon gastric contraction and breathing models.

Objective.Precision radiation therapy requires managing motions of organs at risk that occur during treatment. While methods have been developed for real-time respiratory motion tracking, non-breathing intra-fractional variations (including gastric contractile motion) have seen little attention to date. The purpose of this study is to develop a cyclic gastric contractile motion prediction model to support real-time management during radiotherapy.Approach. The observed short-term reproducibility of gastric contractile motion permitted development of a prediction model that (1) extracts gastric contraction motion phases from few minutes of golden angle stack of stars scanning (at patient positioning), (2) estimate gastric phase of real-time sampled data acquired during treatment delivery to these reconstructed phases and (3) predicting future gastric phase by linear extrapolation using estimation results from step 2 to account for processing and system latency times. Model was evaluated on three parameters including training time window for step 1, number of spokes for real-time sampling data in step 2 and future prediction time. Mainresults. The model was tested on a population of 20 min data samples from 25 scans from 15 patients. The mean prediction error with 10 spokes and 2 min training was 0.3 ± 0.1 mm (0.1-0.7 mm) with 5.1 s future time, slowly rising to 0.6 ± 0.2 mm (0.2-1.1 mm) for 6.8 s future time and then increasing rapidly for longer forward predictions, for an average 3.6 ± 0.5 mm (2.8-4.7 mm) HD95 of gastric motion. Results showed that reducing of train time window (5-2 min) does not influence the prediction performance, while using 5 spokes increased prediction errors.Significance. The proposed gastric motion prediction model has sufficiently accurate prediction performance to allow for sub-millimeter accuracy while allowing sufficient time for data processing and machine interaction and shows the potential for clinical implementation to support stomach motion tracking during radiotherapy.

Contrastive self-supervised learning from 100 million medical images with optional supervision.

Purpose: Building accurate and robust artificial intelligence systems for medical image assessment requires the creation of large sets of annotated training examples. However, constructing such datasets is very costly due to the complex nature of annotation tasks, which often require expert knowledge (e.g., a radiologist). To counter this limitation, we propose a method to learn from medical images at scale in a self-supervised way. Approach: Our approach, based on contrastive learning and online feature clustering, leverages training datasets of over 100,000,000 medical images of various modalities, including radiography, computed tomography (CT), magnetic resonance (MR) imaging, and ultrasonography (US). We propose to use the learned features to guide model training in supervised and hybrid self-supervised/supervised regime on various downstream tasks. Results: We highlight a number of advantages of this strategy on challenging image assessment problems in radiography, CT, and MR: (1) significant increase in accuracy compared to the state-of-the-art (e.g., area under the curve boost of 3% to 7% for detection of abnormalities from chest radiography scans and hemorrhage detection on brain CT); (2) acceleration of model convergence during training by up to 85% compared with using no pretraining (e.g., 83% when training a model for detection of brain metastases in MR scans); and (3) increase in robustness to various image augmentations, such as intensity variations, rotations or scaling reflective of data variation seen in the field. Conclusions: The proposed approach enables large gains in accuracy and robustness on challenging image assessment problems. The improvement is significant compared with other state-of-the-art approaches trained on medical or vision images (e.g., ImageNet).

Gastrointestinal 4D MRI with respiratory motion correction.

PURPOSE: Gastrointestinal motion patterns such as peristalsis and segmental contractions can alter the shape and position of the stomach and intestines with respect to other irradiated organs during radiation therapy. Unfortunately, these deformations are concealed by conventional four-dimensional (4D)-MRI techniques, which were developed to visualize respiratory motion by binning acquired data into respiratory motion states without considering the phases of GI motion. We present a method to reconstruct breathing-compensated images showing the phases of periodic gastric motion and study the effect of this motion on regional anatomical structures. METHODS: Sixty-seven DCE-MRI examinations were performed on patients undergoing MRI simulation for hepatocellular carcinoma using a golden-angle stack-of-stars sequence that collected 2000 radial spokes over 5 min. The collected data were reconstructed using a method with integrated respiratory motion correction into a time series of 3D image volumes without visible breathing motion. From this series, a gastric motion signal was extracted by temporal filtering of time-intensity curves in the stomach. Using this motion signal, breathing-corrected back-projection images were sorted according to the gastric phase and reconstructed into 21 gastric motion state images showing the phases of gastric motion. RESULTS: Reconstructed image volumes showed gastric motion states clearly with no visible breathing motion or related artifacts. The mean frequency of the gastric motion signal was 3 cycles/min with a standard deviation of 0.27 cycles/min. CONCLUSIONS: Periodic gastrointestinal motion can be visualized without confounding respiratory motion using the presented GI 4D MRI technique. GI 4D MRIs may help define internal target volumes for treatment planning, aid in planning organ at risk volume definition, or support motion model development for gastrointestinal motion tracking algorithms for real-time MR-guided radiation therapy.

A hierarchical model of abdominal configuration changes extracted from golden angle radial magnetic resonance imaging.

Abdominal organs are subject to a variety of physiological forces that superimpose their effects to influence local motion and configuration. These forces not only include breathing, but can also arise from cyclic antral contractions and a range of slow configuration changes. To elucidate each individual motion pattern as well as their combined effects, a hierarchical motion model was built for characterization of these 3 motion modes (characterized as deformation maps between states) using golden angle radial MR signals. Breathing motions are characterized first. Antral contraction states are then reconstructed after breathing motion-induced deformation are corrected; slow configuration change states are further extracted from breathing motion-corrected image reconstructions. The hierarchical model is established based on these multimodal states, which can be either individually shown or combined to demonstrate any arbitrary composited motion patterns. The model was evaluated using 20 MR scans acquired from 9 subjects. Poor reproducibility of breathing motions both within as well as between scan sessions was observed, with an average intra-subject difference of 1.6 cycles min-1 for average breathing frequencies of 12.0 cycles min-1. Antral contraction frequency distributions were more stable than breathing, but also presented poor reproducibility between scans with an average difference of 0.3 cycles min-1 for average frequencies of 3.2 cycles min-1. The magnitudes of motions beyond breathing were found to be significant, with 14.4 and 33.8 mm maximal motions measured from antral contraction and slow configuration changes, respectively. Hierarchical motion models have potential in multiple applications in radiotherapy, including improving the accuracy of dose delivery estimation, providing guidance for margin creation, and supporting advanced decisions and strategies for immobilization, treatment monitoring and gating.

Volumetric prediction of breathing and slow drifting motion in the abdomen using radial MRI and multi-temporal resolution modeling.

Abdominal organ motions introduce geometric uncertainties to radiotherapy. This study investigates a multi-temporal resolution 3D motion prediction scheme that accounts for both breathing and slow drifting motion in the abdomen in support of MRI-guided radiotherapy. Ten-minute MRI scans were acquired for 8 patients using a volumetric golden-angle stack-of-stars sequence. The first five-minutes was used for patient-specific motion modeling. Fast breathing motion was modeled from high temporal resolution radial k-space samples, which served as a navigator signal to sort k-space data into different bins for high spatial resolution reconstruction of breathing motion states. Slow drifting motion was modeled from a lower temporal resolution image time series which was reconstructed by sequentially combining a large number of breathing-corrected k-space samples. Principal components analysis (PCA) was performed on deformation fields between different motion states. Gaussian kernel regression and linear extrapolation were used to predict PCA coefficients of future motion states for breathing motion (340 ms ahead of acquisition) and slow drifting motion (8.5 s ahead of acquisition) respectively. k-space data from the remaining five-minutes was used to compare ground truth motions states obtained from retrospective reconstruction/deformation with predictions. Median distances between predicted and ground truth centroid positions of gross tumor volume (GTV) and organs at risk (OARs) were less than 1 mm on average. 95- percentile Hausdorff distances between predicted and ground truth GTV contours of various breathing motions states were 2 mm on average, which was smaller than the imaging resolution and 95-percentile Hausdorff distances between predicted and ground truth OAR contours of different slow drifting motion states were less than 0.2 mm. These results suggest that multi-temporal resolution motion models are capable of volumetric predictions of breathing and slow drifting motion with sufficient accuracy and temporal resolution for MRI-based tracking, and thus have potential for supporting MRI-guided abdominal radiotherapy.

Mapping lung ventilation through stress maps derived from biomechanical models of the lung.

PURPOSE: Most existing computed tomography (CT)-ventilation imaging techniques are based on deformable image registration (DIR) of different respiratory phases of a four-dimensonal CT (4DCT) scan of the lung, followed by the quantification of local breathing-induced changes in Hounsfield Units (HU) or volume. To date, only moderate correlations have been reported between these CT-ventilation metrics and standard ventilation imaging modalities for adaptive lung radiation therapy. This study evaluates the use of stress maps derived from biomechanical model-based DIR as an alternative CT-ventilation metric. MATERIALS AND METHODS: Six patients treated for lung cancer with conventional radiation therapy were retrospectively analyzed. For each patient, a 4DCT scan and Tc-99m SPECT-V image acquired for treatment planning were collected. Biomechanical model-based DIR was applied between the inhale and exhale phase of the 4DCT scans and stress maps were calculated. The voxel-wise correlation between the reference SPECT-V image and map of the maximum principal stress was measured with a Spearman correlation coefficient. The overlap between high (above the 75th percentile) and low (below the 25th percentile) functioning volumes extracted from the reference SPECT-V and the stress maps was measured with Dice similarity coefficients (DSC). The results were compared to those obtained when using two classical CT-ventilation metrics: the change in HU and Jacobian determinant. RESULTS: The mean Spearman correlation coefficients were: 0.37 ± 18 and 0.39 ± 13 and 0.59 ± 0.13 considering the changes in HU, Jacobian and maximum principal stress, respectively. The corresponding mean DSC coefficients were 0.38 ± 0.09, 0.37 ± 0.07 and 0.52 ± 0.07 for the high ventilation function volumes and 0.48 ± 0.13, 0.44 ± 0.09 and 0.52 ± 0.07 for the low ventilation function volumes. CONCLUSION: For presenting a significantly stronger and more consistent correlation with standard SPECT-V images than previously proposed CT-ventilation metrics, stress maps derived with the proposed method appear to be a promising tool for incorporation into functional lung avoidance strategies.

Modeling intra-fractional abdominal configuration changes using breathing motion-corrected radial MRI.

Abdominal organ motions introduce geometric uncertainties to gastrointestinal radiotherapy. This study investigated slow drifting motion induced by changes of internal anatomic organ arrangements using a 3D radial MRI sequence with a scan length of 20 min. Breathing motion and cyclic GI motion were first removed through multi-temporal resolution image reconstruction. Slow drifting motion analysis was performed using an image time series consisting of 72 image volumes with a temporal sampling rate of 17 s. B-spline deformable registration was performed to align image volumes of the time series to a reference volume. The resulting deformation fields were used for motion velocity evaluation and patient-specific motion model construction through principal component analysis (PCA). Geometric uncertainties introduced by slow drifting motion were assessed by Hausdorff distances between unions of organs at risk (OARs) at different motion states and reference OAR contours as well as probabilistic distributions of OARs predicted using the PCA model. Thirteen examinations from 11 patients were included in this study. The averaged motion velocities ranged from 0.8 to 1.9 mm min-1, 0.7 to 1.6 mm min-1, 0.6 to 2.0 mm min-1and 0.7 to 1.4 mm min-1for the small bowel, colon, duodenum and stomach respectively; the averaged Hausdorff distances were 5.6 mm, 5.3 mm, 5.1 mm and 4.6 mm. On average, a margin larger than 4.5 mm was needed to cover a space with OAR occupancy probability higher than 55%. Temporal variations of geometric uncertainties were evaluated by comparing across four 5 min sub-scans extracted from the full scan. Standard deviations of Hausdorff distances across sub-scans were less than 1 mm for most examinations, indicating stability of relative margin estimates from separate time windows. These results suggested slow drifting motion of GI organs is significant and geometric uncertainties introduced by such motion should be accounted for during radiotherapy planning and delivery.

Evaluating deep learning methods in detecting and segmenting different sizes of brain metastases on 3D post-contrast T1-weighted images.

Purpose: We investigate the impact of various deep-learning-based methods for detecting and segmenting metastases with different lesion volume sizes on 3D brain MR images. Approach: A 2.5D U-Net and a 3D U-Net were selected. We also evaluated weak learner fusion of the prediction features generated by the 2.5D and the 3D networks. A 3D fully convolutional one-stage (FCOS) detector was selected as a representative of bounding-box regression-based detection methods. A total of 422 3D post-contrast T1-weighted scans from patients with brain metastases were used. Performances were analyzed based on lesion volume, total metastatic volume per patient, and number of lesions per patient. Results: The performance of detection of the 2.5D and 3D U-Net methods had recall of > 0.83 and precision of > 0.44 for lesion volume > 0.3 cm 3 but deteriorated as metastasis size decreased below 0.3 cm 3 to 0.58 to 0.74 in recall and 0.16 to 0.25 in precision. Compared the two U-Nets for detection capability, high precision was achieved by the 2.5D network, but high recall was achieved by the 3D network for all lesion sizes. The weak learner fusion achieved a balanced performance between the 2.5D and 3D U-Nets; particularly, it increased precision to 0.83 for lesion volumes of 0.1 to 0.3 cm 3 but decreased recall to 0.59. The 3D FCOS detector did not outperform the U-Net methods in detecting either the small or large metastases presumably because of the limited data size. Conclusions: Our study provides the performances of four deep learning methods in relationship to lesion size, total metastasis volume, and number of lesions per patient, providing insight into further development of the deep learning networks.

Quantifying and leveraging predictive uncertainty for medical image assessment.

The interpretation of medical images is a challenging task, often complicated by the presence of artifacts, occlusions, limited contrast and more. Most notable is the case of chest radiography, where there is a high inter-rater variability in the detection and classification of abnormalities. This is largely due to inconclusive evidence in the data or subjective definitions of disease appearance. An additional example is the classification of anatomical views based on 2D Ultrasound images. Often, the anatomical context captured in a frame is not sufficient to recognize the underlying anatomy. Current machine learning solutions for these problems are typically limited to providing probabilistic predictions, relying on the capacity of underlying models to adapt to limited information and the high degree of label noise. In practice, however, this leads to overconfident systems with poor generalization on unseen data. To account for this, we propose a system that learns not only the probabilistic estimate for classification, but also an explicit uncertainty measure which captures the confidence of the system in the predicted output. We argue that this approach is essential to account for the inherent ambiguity characteristic of medical images from different radiologic exams including computed radiography, ultrasonography and magnetic resonance imaging. In our experiments we demonstrate that sample rejection based on the predicted uncertainty can significantly improve the ROC-AUC for various tasks, e.g., by 8% to 0.91 with an expected rejection rate of under 25% for the classification of different abnormalities in chest radiographs. In addition, we show that using uncertainty-driven bootstrapping to filter the training data, one can achieve a significant increase in robustness and accuracy. Finally, we present a multi-reader study showing that the predictive uncertainty is indicative of reader errors.

Stereotactic Transcranial Focused Ultrasound Targeting System for Murine Brain Models.

An inexpensive, accurate focused ultrasound stereotactic targeting method guided by pretreatment magnetic resonance imaging (MRI) images for murine brain models is presented. An uncertainty of each sub-component of the stereotactic system was analyzed. The entire system was calibrated using clot phantoms. The targeting accuracy of the system was demonstrated with an in vivo mouse glioblastoma (GBM) model. The accuracy was quantified by the absolute distance difference between the prescribed and ablated points visible on the pre treatment and posttreatment MR images, respectively. A precalibration phantom study ( N = 6 ) resulted in an error of 0.32 ± 0.31, 0.72 ± 0.16, and 1.06 ± 0.38 mm in axial, lateral, and elevational axes, respectively. A postcalibration phantom study ( N = 8 ) demonstrated a residual error of 0.09 ± 0.01, 0.15 ± 0.09, and 0.47 ± 0.18 mm in axial, lateral, and elevational axes, respectively. The calibrated system showed significantly reduced ( ) error of 0.20 ± 0.21, 0.34 ± 0.24, and 0.28 ± 0.21 mm in axial, lateral, and elevational axes, respectively, in the in vivo GBM tumor-bearing mice ( N = 10 ).

Investigating the SPECT Dose-Function Metrics Associated With Radiation-Induced Lung Toxicity Risk in Patients With Non-small Cell Lung Cancer Undergoing Radiation Therapy.

PURPOSE: Dose to normal lung has commonly been linked with radiation-induced lung toxicity (RILT) risk, but incorporating functional lung metrics in treatment planning may help further optimize dose delivery and reduce RILT incidence. The purpose of this study was to investigate the impact of the dose delivered to functional lung regions by analyzing perfusion (Q), ventilation (V), and combined V/Q single-photon-emission computed tomography (SPECT) dose-function metrics with regard to RILT risk in patients with non-small cell lung cancer (NSCLC) patients who received radiation therapy (RT). METHODS AND MATERIALS: SPECT images acquired from 88 patients with locally advanced NSCLC before undergoing conventionally fractionated RT were retrospectively analyzed. Dose was converted to the nominal dose equivalent per 2 Gy fraction, and SPECT intensities were normalized. Regional lung segments were defined, and the average dose delivered to each lung region was quantified. Three functional categorizations were defined to represent low-, normal-, and high-functioning lungs. The percent of functional lung category receiving ≥20 Gy and mean functional intensity receiving ≥20 Gy (iV20) were calculated. RILT was defined as grade 2+ radiation pneumonitis and/or clinical radiation fibrosis. A logistic regression was used to evaluate the association between dose-function metrics and risk of RILT. RESULTS: By analyzing V/Q normalized intensities and functional distributions across the population, a wide range in functional capability (especially in the ipsilateral lung) was observed in patients with NSCLC before RT. Through multivariable regression models, global lung average dose to the lower lung was found to be significantly associated with RILT, and Q and V iV20 were correlated with RILT when using ipsilateral lung metrics. Through a receiver operating characteristic analysis, combined V/Q low-function receiving ≥20 Gy (low-functioning V/Q20) in the ipsilateral lung was found to be the best predictor (area under the curce: 0.79) of RILT risk. CONCLUSIONS: Irradiation of the inferior lung appears to be a locational sensitivity for RILT risk. The multivariable correlation between ipsilateral lung iV20 and RILT, as well as the association of low-functioning V/Q20 and RILT, suggest that irradiating low-functioning regions in the lung may lead to higher toxicity rates.

Gender Differences in Work-Life Integration Among Medical Physicists.

PURPOSE: To generate an understanding of the primary concerns facing medical physicists regarding integration of a demanding technical career with their personal lives. METHODS AND MATERIALS: In 2019, we recruited 32 medical physics residents, faculty, and staff via emails to US medical physics residency program directors to participate in a 1-hour, semistructured interview that elicited their thoughts on several topics, including work-life integration. Standard techniques of qualitative thematic analysis were used to generate the research findings. RESULTS: Of the participants, 50% were women and 69% were non-Hispanic White individuals, with a mean (SD) age of 37.5 (7.4) years. They were evenly split between residents and faculty or staff. Participant responses centered around 5 primary themes: the gendered distribution of household responsibilities, the effect of career or work on home and family life, the effect of family on career or work, support and strategies for reconciling work-life conflicts, and the role of professional societies in addressing work-life integration. Participants expressed concern about the effect of heavy workloads on home life, with female respondents more likely to report carrying the majority of the household burden. CONCLUSIONS: Medical physicists experience challenges in managing work-life conflict amid a diverse array of personal and professional responsibilities. Further investigations are needed to quantitatively assess the division of work and household labor by gender in medical physics, particularly after the outbreak of the COVID-19 pandemic, but this study's qualitative findings suggest that the profession should consider ways to address root causes of work-life conflict to promote the future success and well-being of all medical physicists, and perhaps women in particular.