Assessment of intrafractional prostate motion and its dosimetric impact in MRI-guided online adaptive radiotherapy with gating.

PURPOSE: This study aimed to evaluate the intrafractional prostate motion captured during gated magnetic resonance imaging (MRI)-guided online adaptive radiotherapy for prostate cancer and analyze its impact on the delivered dose as well as the effect of gating. METHODS: Sagittal 2D cine-MRI scans were acquired at 4 Hz during treatment at a ViewRay MRIdian (ViewRay Inc., Oakwood Village, OH, USA) MR linac. Prostate shifts in anterior-posterior (AP) and superior-inferior (SI) directions were extracted separately. Using the static dose cloud approximation, the planned fractional dose was shifted according to the 2D gated motion (residual motion in gating window) to estimate the delivered dose by superimposing and averaging the shifted dose volumes. The dose of a hypothetical non-gated delivery was reconstructed similarly using the non-gated motion. For the clinical target volume (CTV), rectum, and bladder, dose-volume histogram parameters of the planned and reconstructed doses were compared. RESULTS: In total, 174 fractions (15.7 h of cine-MRI) from 10 patients were evaluated. The average (±1 σ) non-gated prostate motion was 0.6 ± 1.0 mm in the AP and 0.0 ± 0.6 mm in the SI direction with respect to the centroid position of the gating boundary. 95% of the shifts were within [-3.5, 2.7] mm in the AP and [-2.9, 3.2] mm in the SI direction. For the gated treatment and averaged over all fractions, CTV D98% decreased by less than 2% for all patients. The rectum and the bladder D2% increased by less than 3% and 0.5%, respectively. Doses reconstructed for gated and non-gated delivery were similar for most fractions. CONCLUSION: A pipeline for extraction of prostate motion during gated MRI-guided radiotherapy based on 2D cine-MRI was implemented. The 2D motion data enabled an approximate estimation of the delivered dose. For the majority of fractions, the benefit of gating was negligible, and clinical dosimetric constraints were met, indicating safety of the currently adopted gated MRI-guided treatment workflow.

A convolutional neural network with self-attention for fully automated metabolic tumor volume delineation of head and neck cancer in [Formula: see text]F]FDG PET/CT.

PURPOSE: PET-derived metabolic tumor volume (MTV) and total lesion glycolysis of the primary tumor are known to be prognostic of clinical outcome in head and neck cancer (HNC). Including evaluation of lymph node metastases can further increase the prognostic value of PET but accurate manual delineation and classification of all lesions is time-consuming and prone to interobserver variability. Our goal, therefore, was development and evaluation of an automated tool for MTV delineation/classification of primary tumor and lymph node metastases in PET/CT investigations of HNC patients. METHODS: Automated lesion delineation was performed with a residual 3D U-Net convolutional neural network (CNN) incorporating a multi-head self-attention block. 698 [Formula: see text]F]FDG PET/CT scans from 3 different sites and 5 public databases were used for network training and testing. An external dataset of 181 [Formula: see text]F]FDG PET/CT scans from 2 additional sites was employed to assess the generalizability of the network. In these data, primary tumor and lymph node (LN) metastases were interactively delineated and labeled by two experienced physicians. Performance of the trained network models was assessed by 5-fold cross-validation in the main dataset and by pooling results from the 5 developed models in the external dataset. The Dice similarity coefficient (DSC) for individual delineation tasks and the primary tumor/metastasis classification accuracy were used as evaluation metrics. Additionally, a survival analysis using univariate Cox regression was performed comparing achieved group separation for manual and automated delineation, respectively. RESULTS: In the cross-validation experiment, delineation of all malignant lesions with the trained U-Net models achieves DSC of 0.885, 0.805, and 0.870 for primary tumor, LN metastases, and the union of both, respectively. In external testing, the DSC reaches 0.850, 0.724, and 0.823 for primary tumor, LN metastases, and the union of both, respectively. The voxel classification accuracy was 98.0% and 97.9% in cross-validation and external data, respectively. Univariate Cox analysis in the cross-validation and the external testing reveals that manually and automatically derived total MTVs are both highly prognostic with respect to overall survival, yielding essentially identical hazard ratios (HR) ([Formula: see text]; [Formula: see text] vs. [Formula: see text]; [Formula: see text] in cross-validation and [Formula: see text]; [Formula: see text] vs. [Formula: see text]; [Formula: see text] in external testing). CONCLUSION: To the best of our knowledge, this work presents the first CNN model for successful MTV delineation and lesion classification in HNC. In the vast majority of patients, the network performs satisfactory delineation and classification of primary tumor and lymph node metastases and only rarely requires more than minimal manual correction. It is thus able to massively facilitate study data evaluation in large patient groups and also does have clear potential for supervised clinical application.

ExacTrac Dynamic workflow evaluation: Combined surface optical/thermal imaging and X-ray positioning.

In modern radiotherapy (RT), especially for stereotactic radiotherapy or stereotactic radiosurgery treatments, image guidance is essential. Recently, the ExacTrac Dynamic (EXTD) system, a new combined surface-guided RT and image-guided RT (IGRT) system for patient positioning, monitoring, and tumor targeting, was introduced in clinical practice. The purpose of this study was to provide more information about the geometric accuracy of EXTD and its workflow in a clinical environment. The surface optical/thermal- and the stereoscopic X-ray imaging positioning systems of EXTD was evaluated and compared to cone-beam computed tomography (CBCT). Additionally, the congruence with the radiation isocenter was tested. A Winston Lutz test was executed several times over 1 year, and repeated end-to-end positioning tests were performed. The magnitude of the displacements between all systems, CBCT, stereoscopic X-ray, optical-surface imaging, and MV portal imaging was within the submillimeter range, suggesting that the image guidance provided by EXTD is accurate at any couch angle. Additionally, results from the evaluation of 14 patients with intracranial tumors treated with open-face masks are reported, and limited differences with a maximum of 0.02 mm between optical/thermal- and stereoscopic X-ray imaging were found. As the optical/thermal positioning system showed a comparable accuracy to other IGRT systems, and due to its constant monitoring capability, it can be an efficient tool for detecting intra-fractional motion and for real-time tracking of the surface position during RT.

Evaluation of proton and photon dose distributions recalculated on 2D and 3D Unet-generated pseudoCTs from T1-weighted MR head scans.

Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.

Comparing cone-beam CT intensity correction methods for dose recalculation in adaptive intensity-modulated photon and proton therapy for head and neck cancer.

BACKGROUND: Adaptive intensity-modulated photon and proton radiotherapy (IMRT and IMPT) of head and neck (H&N) cancer requires frequent three-dimensional (3D) dose calculation. We compared two approaches for dose recalculation on the basis of intensity-corrected cone-beam (CB) x-ray computed tomography (CT) images. MATERIAL AND METHODS: For nine H&N tumor patients, virtual CTs (vCT) were generated by deformable image registration of the planning CT (pCT) to the CBCT. The second intensity correction approach used population-based lookup tables for scaling CBCT intensities to the pCT HU range (CBCTLUT). IMRT and IMPT plans were generated with a commercial treatment planning system. Dose recalculations on vCT and CBCTLUT were analyzed using a (3%, 3 mm) gamma-index analysis and comparison of normal tissue and tumor dose/volume parameters. A replanning CT (rpCT) acquired within three days of the CBCT served as reference. Single field uniform dose (SFUD) proton plans were created and recalculated on vCT and CBCTLUT for proton range comparison. RESULTS: Dose/volume parameters showed minor differences between rpCT, vCT and CBCTLUT in IMRT, but clinically relevant deviations between CBCTLUT and rpCT in the spinal cord for IMPT. Gamma-index pass-rates were found increased for vCT with respect to CBCTLUT in IMPT (by up to 21 percentage points) and IMRT (by up to 9 percentage points) for most cases. The SFUD-based proton range assessment showed improved agreement of vCT and rpCT, with 88-99% of the depth dose profiles in beam's eye view agreeing within 3 mm. For CBCTLUT, only 80-94% of the profiles fulfilled this criterion. CONCLUSION: vCT and CBCTLUT are suitable options for dose recalculation in adaptive IMRT. In the scope of IMPT, the vCT approach is preferable.

A simulation study on proton computed tomography (CT) stopping power accuracy using dual energy CT scans as benchmark.

BACKGROUND: Accurate stopping power estimation is crucial for treatment planning in proton therapy, and the uncertainties in stopping power are currently the largest contributor to the employed dose margins. Dual energy x-ray computed tomography (CT) (clinically available) and proton CT (in development) have both been proposed as methods for obtaining patient stopping power maps. The purpose of this work was to assess the accuracy of proton CT using dual energy CT scans of phantoms to establish reference accuracy levels. MATERIAL AND METHODS: A CT calibration phantom and an abdomen cross section phantom containing inserts were scanned with dual energy and single energy CT with a state-of-the-art dual energy CT scanner. Proton CT scans were simulated using Monte Carlo methods. The simulations followed the setup used in current prototype proton CT scanners and included realistic modeling of detectors and the corresponding noise characteristics. Stopping power maps were calculated for all three scans, and compared with the ground truth stopping power from the phantoms. RESULTS: Proton CT gave slightly better stopping power estimates than the dual energy CT method, with root mean square errors of 0.2% and 0.5% (for each phantom) compared to 0.5% and 0.9%. Single energy CT root mean square errors were 2.7% and 1.6%. Maximal errors for proton, dual energy and single energy CT were 0.51%, 1.7% and 7.4%, respectively. CONCLUSION: Better stopping power estimates could significantly reduce the range errors in proton therapy, but requires a large improvement in current methods which may be achievable with proton CT.

Intra-frame motion deterioration effects and deep-learning-based compensation in MR-guided radiotherapy.

BACKGROUND: Current commercially available hybrid magnetic resonance linear accelerators (MR-Linac) use 2D+t cine MR imaging to provide intra-fractional motion monitoring. However, given the limited temporal resolution of cine MR imaging, target intra-frame motion deterioration effects, resulting in effective time latency and motion artifacts in the image domain, can be appreciable, especially in the case of fast breathing. PURPOSE: The aim of this work is to investigate intra-frame motion deterioration effects in MR-guided radiotherapy (MRgRT) by simulating the motion-corrupted image acquisition, and to explore the feasibility of deep-learning-based compensation approaches, relying on the intra-frame motion information which is spatially and temporally encoded in the raw data (k-space). METHODS: An intra-frame motion model was defined to simulate motion-corrupted MR images, with 4D anthropomorphic digital phantoms being exploited to provide ground truth 2D+t cine MR sequences. A total number of 10 digital phantoms were generated for lung cancer patients, with randomly selected eight patients for training or validation and the remaining two for testing. The simulation code served as the data generator, and a dedicated motion pattern perturbation scheme was proposed to build the intra-frame motion database, where three degrees of freedom were designed to guarantee the diversity of intra-frame motion trajectories, enabling a thorough exploration in the domain of the potential anatomical structure positions. U-Nets with three types of loss functions: L1 or L2 loss defined in image or Fourier domain, referred to as NNImgLoss-L1 , NNFloss-L1 and NNL2-Loss were trained to extract information from the motion-corrupted image and used to estimate the ground truth final-position image, corresponding to the end of the acquisition. Images before and after compensation were evaluated in terms of (i) image mean-squared error (MSE) and mean absolute error (MAE), and (ii) accuracy of gross tumor volume (GTV) contouring, based on optical-flow image registration. RESULTS: Image degradation caused by intra-frame motion was observed: for a linearly and fully acquired Cartesian readout k-space trajectory, intra-frame motion resulted in an imaging latency of approximately 50% of the acquisition time; in comparison, the motion artifacts exhibited only a negligible contribution to the overall geometric errors. All three compensation models led to a decrease in image MSE/MAE and GTV position offset compared to the motion-corrupted image. In the investigated testing dataset for GTV contouring, the average dice similarity coefficients (DSC) improved from 88% to 96%, and the 95th percentile Hausdorff distance (HD95 ) dropped from 4.8 mm to 2.1 mm. Different models showed slight performance variations across different intra-frame motion amplitude categories: NNImgLoss-L1 excelled for small/medium amplitudes, whereas NNFloss-L1 demonstrated higher DSC median values at larger amplitudes. The saliency maps of the motion-corrupted image highlighted the major contribution of the later acquired k-space data, as well as the edges of the moving anatomical structures at their final positions, during the model inference stage. CONCLUSIONS: Our results demonstrate the deep-learning-based approaches have the potential to compensate for intra-frame motion by utilizing the later acquired data to drive the convergence of the earlier acquired k-space components.

Repeatability quantification of brain diffusion-weighted imaging for future clinical implementation at a low-field MR-linac.

BACKGROUND: Longitudinal assessments of apparent diffusion coefficients (ADCs) derived from diffusion-weighted imaging (DWI) during intracranial radiotherapy at magnetic resonance imaging-guided linear accelerators (MR-linacs) could enable early response assessment by tracking tumor diffusivity changes. However, DWI pulse sequences are currently unavailable in clinical practice at low-field MR-linacs. Quantifying the in vivo repeatability of ADC measurements is a crucial step towards clinical implementation of DWI sequences but has not yet been reported on for low-field MR-linacs. This study assessed ADC measurement repeatability in a phantom and in vivo at a 0.35 T MR-linac. METHODS: Eleven volunteers and a diffusion phantom were imaged on a 0.35 T MR-linac. Two echo-planar imaging DWI sequence variants, emphasizing high spatial resolution ("highRes") and signal-to-noise ratio ("highSNR"), were investigated. A test-retest study with an intermediate outside-scanner-break was performed to assess repeatability in the phantom and volunteers' brains. Mean ADCs within phantom vials, cerebrospinal fluid (CSF), and four brain tissue regions were compared to literature values. Absolute relative differences of mean ADCs in pre- and post-break scans were calculated for the diffusion phantom, and repeatability coefficients (RC) and relative RC (relRC) with 95% confidence intervals were determined for each region-of-interest (ROI) in volunteers. RESULTS: Both DWI sequence variants demonstrated high repeatability, with absolute relative deviations below 1% for water, dimethyl sulfoxide, and polyethylene glycol in the diffusion phantom. RelRCs were 7% [5%, 12%] (CSF; highRes), 12% [9%, 22%] (CSF; highSNR), 9% [8%, 12%] (brain tissue ROIs; highRes), and 6% [5%, 7%] (brain tissue ROIs; highSNR), respectively. ADCs measured with the highSNR variant were consistent with literature values for volunteers, while smaller mean values were measured for the diffusion phantom. Conversely, the highRes variant underestimated ADCs compared to literature values, indicating systematic deviations. CONCLUSIONS: High repeatability of ADC measurements in a diffusion phantom and volunteers' brains were measured at a low-field MR-linac. The highSNR variant outperformed the highRes variant in accuracy and repeatability, at the expense of an approximately doubled voxel volume. The observed high in vivo repeatability confirms the potential utility of DWI at low-field MR-linacs for early treatment response assessment.

Simultaneous object detection and segmentation for patient-specific markerless lung tumor tracking in simulated radiographs with deep learning.

BACKGROUND: Real-time tumor tracking is one motion management method to address motion-induced uncertainty. To date, fiducial markers are often required to reliably track lung tumors with X-ray imaging, which carries risks of complications and leads to prolonged treatment time. A markerless tracking approach is thus desirable. Deep learning-based approaches have shown promise for markerless tracking, but systematic evaluation and procedures to investigate applicability in individual cases are missing. Moreover, few efforts have been made to provide bounding box prediction and mask segmentation simultaneously, which could allow either rigid or deformable multi-leaf collimator tracking. PURPOSE: The purpose of this study was to implement a deep learning-based markerless lung tumor tracking model exploiting patient-specific training which outputs both a bounding box and a mask segmentation simultaneously. We also aimed to compare the two kinds of predictions and to implement a specific procedure to understand the feasibility of markerless tracking on individual cases. METHODS: We first trained a Retina U-Net baseline model on digitally reconstructed radiographs (DRRs) generated from a public dataset containing 875 CT scans and corresponding lung nodule annotations. Afterwards, we used an independent cohort of 97 lung patients to develop a patient-specific refinement procedure. In order to determine the optimal hyperparameters for automatic patient-specific training, we selected 13 patients for validation where the baseline model predicted a bounding box on planning CT (PCT)-DRR with intersection over union (IoU) with the ground-truth higher than 0.7. The final test set contained the remaining 84 patients with varying PCT-DRR IoU. For each testing patient, the baseline model was refined on the PCT-DRR to generate a patient-specific model, which was then tested on a separate 10-phase 4DCT-DRR to mimic the intrafraction motion during treatment. A template matching algorithm served as benchmark model. The testing results were evaluated by four metrics: the center of mass (COM) error and the Dice similarity coefficient (DSC) for segmentation masks, and the center of box (COB) error and the DSC for bounding box detections. Performance was compared to the benchmark model including statistical testing for significance. RESULTS: A PCT-DRR IoU value of 0.2 was shown to be the threshold dividing inconsistent (68%) and consistent (100%) success (defined as mean bounding box DSC > 0.6) of PS models on 4DCT-DRRs. Thirty-seven out of the eighty-four testing cases had a PCT-DRR IoU above 0.2. For these 37 cases, the mean COM error was 2.6 mm, the mean segmentation DSC was 0.78, the mean COB error was 2.7 mm, and the mean box DSC was 0.83. Including the validation cases, the model was applicable to 50 out of 97 patients when using the PCT-DRR IoU threshold of 0.2. The inference time per frame was 170 ms. The model outperformed the benchmark model on all metrics, and the comparison was significant (p < 0.001) over the 37 PCT-DRR IoU > 0.2 cases, but not over the undifferentiated 84 testing cases. CONCLUSIONS: The implemented patient-specific refinement approach based on a pre-trained baseline model was shown to be applicable to markerless tumor tracking in simulated radiographs for lung cases.

Comparison of deep learning networks for fully automated head and neck tumor delineation on multi-centric PET/CT images.

OBJECTIVES: Deep learning-based auto-segmentation of head and neck cancer (HNC) tumors is expected to have better reproducibility than manual delineation. Positron emission tomography (PET) and computed tomography (CT) are commonly used in tumor segmentation. However, current methods still face challenges in handling whole-body scans where a manual selection of a bounding box may be required. Moreover, different institutions might still apply different guidelines for tumor delineation. This study aimed at exploring the auto-localization and segmentation of HNC tumors from entire PET/CT scans and investigating the transferability of trained baseline models to external real world cohorts. METHODS: We employed 2D Retina Unet to find HNC tumors from whole-body PET/CT and utilized a regular Unet to segment the union of the tumor and involved lymph nodes. In comparison, 2D/3D Retina Unets were also implemented to localize and segment the same target in an end-to-end manner. The segmentation performance was evaluated via Dice similarity coefficient (DSC) and Hausdorff distance 95th percentile (HD95). Delineated PET/CT scans from the HECKTOR challenge were used to train the baseline models by 5-fold cross-validation. Another 271 delineated PET/CTs from three different institutions (MAASTRO, CRO, BERLIN) were used for external testing. Finally, facility-specific transfer learning was applied to investigate the improvement of segmentation performance against baseline models. RESULTS: Encouraging localization results were observed, achieving a maximum omnidirectional tumor center difference lower than 6.8 cm for external testing. The three baseline models yielded similar averaged cross-validation (CV) results with a DSC in a range of 0.71-0.75, while the averaged CV HD95 was 8.6, 10.7 and 9.8 mm for the regular Unet, 2D and 3D Retina Unets, respectively. More than a 10% drop in DSC and a 40% increase in HD95 were observed if the baseline models were tested on the three external cohorts directly. After the facility-specific training, an improvement in external testing was observed for all models. The regular Unet had the best DSC (0.70) for the MAASTRO cohort, and the best HD95 (7.8 and 7.9 mm) in the MAASTRO and CRO cohorts. The 2D Retina Unet had the best DSC (0.76 and 0.67) for the CRO and BERLIN cohorts, and the best HD95 (12.4 mm) for the BERLIN cohort. CONCLUSION: The regular Unet outperformed the other two baseline models in CV and most external testing cohorts. Facility-specific transfer learning can potentially improve HNC segmentation performance for individual institutions, where the 2D Retina Unets could achieve comparable or even better results than the regular Unet.

Intrafractional motion detection for spine SBRT via X-ray imaging using ExacTrac Dynamic.

Purpose: Due to its close vicinity to critical structures, especially the spinal cord, standards for safety for spine stereotactic body radiotherapy (SBRT) should be high. This study was conducted, to evaluate intrafractional motion during spine SBRT for patients without individualized immobilization (e.g., vacuum cushions) using high accuracy patient monitoring via orthogonal X-ray imaging. Methods: Intrafractional X-ray data were collected from 29 patients receiving 79 fractions of spine SBRT. No individualized immobilization devices were used during the treatment. Intrafractional motion was monitored using the ExacTrac Dynamic (ETD) System (Brainlab AG, Munich, Germany). Deviations were detected in six degrees of freedom (6 DOF). Tolerances for repositioning were 0.7 mm for translational and 0.5° for rotational deviations. Patients were repositioned when the tolerance levels were exceeded. Results: Out of the 925 pairs of stereoscopic X-ray images examined, 138 (15 %) showed at least one deviation exceeding the predefined tolerance values. In all 6 DOF together, a total of 191 deviations out of tolerance were recorded. The frequency of deviations exceeding the tolerance levels varied among patients but occurred in all but one patient. Deviations out of tolerance could be seen in all 6 DOF. Maximum translational deviations were 2.6 mm, 2.3 mm and 2.8 mm in the lateral, longitudinal and vertical direction. Maximum rotational deviations were 1.8°, 2.6° and 1.6° for pitch, roll and yaw, respectively. Translational deviations were more frequent than rotational ones, and frequency and magnitude of deviations showed an inverse correlation. Conclusion: Intrafractional motion detection and patient repositioning during spine SBRT using X-ray imaging via the ETD System can lead to improved safety during the application of high BED in critical locations. When using intrafractional imaging with low thresholds for re-positioning individualized immobilization devices (e.g. vacuum cushions) may be omitted.

Impact of daily plan adaptation on accumulated doses in ultra-hypofractionated magnetic resonance-guided radiation therapy of prostate cancer.

Background and purpose: Ultra-hypofractionated online adaptive magnetic resonance-guided radiotherapy (MRgRT) is promising for prostate cancer. However, the impact of online adaptation on target coverage and organ-at-risk (OAR) sparing at the level of accumulated dose has not yet been reported. Using deformable image registration (DIR)-based accumulation, we compared the delivered adapted dose with the simulated non-adapted dose. Materials and methods: Twenty-three prostate cancer patients treated at two clinics with 0.35 T magnetic resonance-guided linear accelerator (MR-linac) following the same treatment protocol (5 × 7.5 Gy with urethral sparing and daily adaptation) were included. The fraction MR images were deformably registered to the planning MR image. Both non-adapted and adapted fraction doses were accumulated with the corresponding vector fields. Two DIR approaches were implemented. PTV* (planning target volume minus urethra+2mm) D95%, CTV* (clinical target volume minus urethra) D98%, and OARs (urethra+2mm, bladder, and rectum) D0.2cc, were evaluated. Statistical significance was inferred from a two-tailed Wilcoxon signed-rank test (p < 0.05). Results: Normalized to the baseline, the accumulated PTV* D95% increased significantly by 2.7 % ([1.5, 4.3]%) through adaptation, and the CTV* D98% by 1.2 % ([0.1, 1.7]%). For the OARs after adaptation, accumulated bladder D0.2cc decreased by 0.4 % ([-1.2, 0.4]%), urethra+2mmD0.2cc by 0.8 % ([-1.6, -0.1]%), while rectum D0.2cc increased by 2.6 % ([1.2, 4.9]%). For all patients, rectum D0.2cc was still below the clinical constraint. Results of both DIR approaches differed on average by less than 0.2 %. Conclusions: Online adaptation in MRgRT improved target coverage and OARs sparing at the level of accumulated dose.

Reduction of cone-beam CT artifacts in a robotic CBCT device using saddle trajectories with integrated infrared tracking.

BACKGROUND: Cone beam computed tomography (CBCT) is widely used in many medical fields. However, conventional CBCT circular scans suffer from cone beam (CB) artifacts that limit the quality and reliability of the reconstructed images due to incomplete data. PURPOSE: Saddle trajectories in theory might be able to improve the CBCT image quality by providing a larger region with complete data. Therefore, we investigated the feasibility and performance of saddle trajectory CBCT scans and compared them to circular trajectory scans. METHODS: We performed circular and saddle trajectory scans using a novel robotic CBCT scanner (Mobile ImagingRing (IRm); medPhoton, Salzburg, Austria). For the saddle trajectory, the gantry executed yaw motion up to ± 10 ∘ $\pm 10^{\circ }$ using motorized wheels driving on the floor. An infrared (IR) tracking device with reflective markers was used for online geometric calibration correction (mainly floor unevenness). All images were reconstructed using penalized least-squares minimization with the conjugate gradient algorithm from RTK with 0.5 × 0.5 × 0.5 mm 3 $0.5 \times 0.5\times 0.5 \text{ mm}^3$ voxel size. A disk phantom and an Alderson phantom were scanned to assess the image quality. Results were correlated with the local incompleteness value represented by tan ( ψ ) $\tan (\psi)$ , which was calculated at each voxel as a function of the source trajectory and the voxel's 3D coordinates. We assessed the magnitude of CB artifacts using the full width half maximum (FWHM) of each disk profile in the axial center of the reconstructed images. Spatial resolution was also quantified by the modulation transfer function at 10% (MTF10). RESULTS: When using the saddle trajectory, the region without CB artifacts was increased from 43 to 190 mm in the SI direction compared to the circular trajectory. This region coincided with low values for tan ( ψ ) $\tan (\psi)$ . When tan ( ψ ) $\tan (\psi)$ was larger than 0.02, we found there was a linear relationship between the FWHM and tan ( ψ ) $\tan (\psi)$ . For the saddle, IR tracking allowed the increase of MTF10 from 0.37 to 0.98 lp/mm. CONCLUSIONS: We achieved saddle trajectory CBCT scans with a novel CBCT system combined with IR tracking. The results show that the saddle trajectory provides a larger region with reliable reconstruction compared to the circular trajectory. The proposed method can be used to evaluate other non-circular trajectories.

Real-time motion management in MRI-guided radiotherapy: Current status and AI-enabled prospects.

MRI-guided radiotherapy (MRIgRT) is a highly complex treatment modality, allowing adaptation to anatomical changes occurring from one treatment day to the other (inter-fractional), but also to motion occurring during a treatment fraction (intra-fractional). In this vision paper, we describe the different steps of intra-fractional motion management during MRIgRT, from imaging to beam adaptation, and the solutions currently available both clinically and at a research level. Furthermore, considering the latest developments in the literature, a workflow is foreseen in which motion-induced over- and/or under-dosage is compensated in 3D, with minimal impact to the radiotherapy treatment time. Considering the time constraints of real-time adaptation, a particular focus is put on artificial intelligence (AI) solutions as a fast and accurate alternative to conventional algorithms.

The role of online MR-guided multi-fraction stereotactic ablative radiotherapy in lung tumours.

Background: The aim of this prospective observational study was to evaluate the dosimetry benefits, changes in pulmonary function, and clinical outcome of online adaptive MR-guided SBRT. Methods: From 11/2020-07/2022, 45 consecutive patients with 59 lesions underwent multi-fraction SBRT (3-8 fractions) at our institution. Patients were eligible if they had biopsy-proven NSCLC or lung cancer/metastases diagnosed via clinical imaging. Endpoints were local control (LC) and overall survival (OS). We evaluated PTV/GTV dose coverage, organs at risk exposure, and changes in pulmonary function (PF). Acute toxicity was classified per the National Cancer Institute-Common Terminology Criteria for Adverse Events version 5.0. Results: The median PTV was 14.4 cm3 (range: 3.4 - 96.5 cm3). In total 195/215 (91%) plans were reoptimised. In the reoptimised vs. predicted plans, PTV coverage by the prescribed dose increased in 94.6% of all fractions with a median increase in PTV VPD of 5.6% (range: -1.8 - 44.6%, p < 0.001), increasing the number of fractions with PTV VPD ≥ 95% from 33% to 98%. The PTV D95% and D98% (BED10) increased in 93% and 95% of all fractions with a median increase of 7.7% (p < 0.001) and 10.6% (p < 0.001). The PTV D95% (BED10) increased by a mean of 9.6 Gy (SD: 10.3 Gy, p < 0.001). At a median follow-up of 21.4 months (95% CI: 12.3-27.0 months), 1- and 2-year LC rates were 94.8% (95% CI: 87.6 - 100.0%) and 91.1% (95% CI: 81.3 - 100%); 1- and 2-year OS rates were 85.6% (95% CI: 75.0 - 96.3%) and 67.1 % (95% CI: 50.3 - 83.8%). One grade ≥ 3 toxicity and no significant reduction in short-term PF parameters were recorded. Conclusions: Online adaptive MR-guided SBRT is an effective, safe and generally well tolerated treatment option for lung tumours achieving encouraging local control rates with significantly improved target volume coverage.

Auto-segmentation of pelvic organs at risk on 0.35T MRI using 2D and 3D Generative Adversarial Network models.

PURPOSE: Manual recontouring of targets and Organs At Risk (OARs) is a time-consuming and operator-dependent task. We explored the potential of Generative Adversarial Networks (GAN) to auto-segment the rectum, bladder and femoral heads on 0.35T MRIs to accelerate the online MRI-guided-Radiotherapy (MRIgRT) workflow. METHODS: 3D planning MRIs from 60 prostate cancer patients treated with 0.35T MR-Linac were collected. A 3D GAN architecture and its equivalent 2D version were trained, validated and tested on 40, 10 and 10 patients respectively. The volumetric Dice Similarity Coefficient (DSC) and 95th percentile Hausdorff Distance (HD95th) were computed against expert drawn ground-truth delineations. The networks were also validated on an independent external dataset of 16 patients. RESULTS: In the internal test set, the 3D and 2D GANs showed DSC/HD95th of 0.83/9.72 mm and 0.81/10.65 mm for the rectum, 0.92/5.91 mm and 0.85/15.72 mm for the bladder, and 0.94/3.62 mm and 0.90/9.49 mm for the femoral heads. In the external test set, the performance was 0.74/31.13 mm and 0.72/25.07 mm for the rectum, 0.92/9.46 mm and 0.88/11.28 mm for the bladder, and 0.89/7.00 mm and 0.88/10.06 mm for the femoral heads. The 3D and 2D GANs required on average 1.44 s and 6.59 s respectively to generate the OARs' volumetric segmentation for a single patient. CONCLUSIONS: The proposed 3D GAN auto-segments pelvic OARs with high accuracy on 0.35T, in both the internal and the external test sets, outperforming its 2D equivalent in both segmentation robustness and volume generation time.

Impact of bias field correction on 0.35 T pelvic MR images: evaluation on generative adversarial network-based OARs' auto-segmentation and visual grading assessment.

Purpose: Magnetic resonance imaging (MRI)-guided radiotherapy enables adaptive treatment plans based on daily anatomical changes and accurate organ visualization. However, the bias field artifact can compromise image quality, affecting diagnostic accuracy and quantitative analyses. This study aims to assess the impact of bias field correction on 0.35 T pelvis MRIs by evaluating clinical anatomy visualization and generative adversarial network (GAN) auto-segmentation performance. Materials and methods: 3D simulation MRIs from 60 prostate cancer patients treated on MR-Linac (0.35 T) were collected and preprocessed with the N4ITK algorithm for bias field correction. A 3D GAN architecture was trained, validated, and tested on 40, 10, and 10 patients, respectively, to auto-segment the organs at risk (OARs) rectum and bladder. The GAN was trained and evaluated either with the original or the bias-corrected MRIs. The Dice similarity coefficient (DSC) and 95th percentile Hausdorff distance (HD95th) were computed for the segmented volumes of each patient. The Wilcoxon signed-rank test assessed the statistical difference of the metrics within OARs, both with and without bias field correction. Five radiation oncologists blindly scored 22 randomly chosen patients in terms of overall image quality and visibility of boundaries (prostate, rectum, bladder, seminal vesicles) of the original and bias-corrected MRIs. Bennett's S score and Fleiss' kappa were used to assess the pairwise interrater agreement and the interrater agreement among all the observers, respectively. Results: In the test set, the GAN trained and evaluated on original and bias-corrected MRIs showed DSC/HD95th of 0.92/5.63 mm and 0.92/5.91 mm for the bladder and 0.84/10.61 mm and 0.83/9.71 mm for the rectum. No statistical differences in the distribution of the evaluation metrics were found neither for the bladder (DSC: p = 0.07; HD95th: p = 0.35) nor for the rectum (DSC: p = 0.32; HD95th: p = 0.63). From the clinical visual grading assessment, the bias-corrected MRI resulted mostly in either no change or an improvement of the image quality and visualization of the organs' boundaries compared with the original MRI. Conclusion: The bias field correction did not improve the anatomy visualization from a clinical point of view and the OARs' auto-segmentation outputs generated by the GAN.

A joint ESTRO and AAPM guideline for development, clinical validation and reporting of artificial intelligence models in radiation therapy.

BACKGROUND AND PURPOSE: Artificial Intelligence (AI) models in radiation therapy are being developed with increasing pace. Despite this, the radiation therapy community has not widely adopted these models in clinical practice. A cohesive guideline on how to develop, report and clinically validate AI algorithms might help bridge this gap. METHODS AND MATERIALS: A Delphi process with all co-authors was followed to determine which topics should be addressed in this comprehensive guideline. Separate sections of the guideline, including Statements, were written by subgroups of the authors and discussed with the whole group at several meetings. Statements were formulated and scored as highly recommended or recommended. RESULTS: The following topics were found most relevant: Decision making, image analysis, volume segmentation, treatment planning, patient specific quality assurance of treatment delivery, adaptive treatment, outcome prediction, training, validation and testing of AI model parameters, model availability for others to verify, model quality assurance/updates and upgrades, ethics. Key references were given together with an outlook on current hurdles and possibilities to overcome these. 19 Statements were formulated. CONCLUSION: A cohesive guideline has been written which addresses main topics regarding AI in radiation therapy. It will help to guide development, as well as transparent and consistent reporting and validation of new AI tools and facilitate adoption.

CT-scan prediction of thyroid cartilage invasion for early laryngeal squamous cell carcinoma.

Treatment choice for laryngeal cancer may be influenced by the diagnosis of thyroid cartilage invasion on preoperative computed tomography (CT). Our objective was to determine the predictive value of CT for thyroid cartilage invasion in early- to mid-stage laryngeal cancer. Retrospective study (1992-2008) of laryngeal squamous cell carcinoma treated with open partial laryngectomy and resection of at least part of the thyroid cartilage. Previous laser surgery, radiation therapy, chemotherapy and second primaries were excluded. CT prediction of thyroid cartilage invasion was determined by specialized radiologists. Tumor characteristics and pathologic thyroid cartilage invasion were compared to the radiologic assessment. 236 patients were treated by vertical (20 %), supracricoid (67 %) or supraglottic partial laryngectomy (13 %) for tumors staged cT1 (26 %), cT2 (55 %), and cT3 (19 %). The thyroid cartilage was invaded on pathology in 19 cases (8 %). CT's sensitivity was 10.5 %, specificity 94 %, positive predictive value 13 %, and negative predictive value 92 %. CT correctly predicted thyroid cartilage invasion in only two cases for an overall accuracy of 87 %. Among the false-positive CT's, tumors involving the anterior commissure were significantly over-represented (61.5 % vs. 27 %, p = .004). Tumors with decreased vocal fold (VF) mobility were significantly over-represented in the group of false-negatives (41 vs. 13 %, p = .0035). Preoperative CT was not effective in predicting thyroid cartilage invasion in these early- to mid-stage lesions, overestimating cartilage invasion for AC lesions and underestimating invasion for lesions with decreased VF mobility.

The role of artificial intelligence in radiotherapy clinical practice.

This review article visits the current state of artificial intelligence (AI) in radiotherapy clinical practice. We will discuss how AI has a place in the modern radiotherapy workflow at the level of automatic segmentation and planning, two applications which have seen real-work implementation. A special emphasis will be placed on the role AI can play in online adaptive radiotherapy, such as performed at MR-linacs, where online plan adaptation is a procedure which could benefit from automation to reduce on-couch time for patients. Pseudo-CT generation and AI for motion tracking will be introduced in the scope of online adaptive radiotherapy as well. We further discuss the use of AI for decision-making and response assessment, for example for personalized prescription and treatment selection, risk stratification for outcomes and toxicities, and AI for quantitative imaging and response assessment. Finally, the challenges of generalizability and ethical aspects will be covered. With this, we provide a comprehensive overview of the current and future applications of AI in radiotherapy.