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 phase space model of a Versa HD linear accelerator for application to Monte Carlo dose calculation in a real-time adaptive workflow.

PURPOSE: This study aims to develop and validate a simple geometric model of the accelerator head, from which a particle phase space can be calculated for application to fast Monte Carlo dose calculation in real-time adaptive photon radiotherapy. With this objective in view, the study investigates whether the phase space model can facilitate dose calculations which are compatible with those of a commercial treatment planning system, for convenient interoperability. MATERIALS AND METHODS: A dual-source model of the head of a Versa HD accelerator (Elekta AB, Stockholm, Sweden) was created. The model used parameters chosen to be compatible with those of 6-MV flattened and 6-MV flattening filter-free photon beams in the RayStation treatment planning system (RaySearch Laboratories, Stockholm, Sweden). The phase space model was used to calculate a photon phase space for several treatment plans, and the resulting phase space was applied to the Dose Planning Method (DPM) Monte Carlo dose calculation algorithm. Simple fields and intensity-modulated radiation therapy (IMRT) treatment plans for prostate and lung were calculated for benchmarking purposes and compared with the convolution-superposition dose calculation within RayStation. RESULTS: For simple square fields in a water phantom, the calculated dose distribution agrees to within ±2% with that from the commercial treatment planning system, except in the buildup region, where the DPM code does not model the electron contamination. For IMRT plans of prostate and lung, agreements of ±2% and ±6%, respectively, are found, with slightly larger differences in the high dose gradients. CONCLUSIONS: The phase space model presented allows convenient calculation of a phase space for application to Monte Carlo dose calculation, with straightforward translation of beam parameters from the RayStation beam model. This provides a basis on which to develop dose calculation in a real-time adaptive setting.

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol.

PURPOSE: Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time. METHODS: Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys. RESULTS: Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively. CONCLUSIONS: We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

A comprehensive model for heat-induced radio-sensitisation.

Combined radiotherapy (RT) and hyperthermia (HT) treatments may improve treatment outcome by heat induced radio-sensitisation. We propose an empirical cell survival model (AlphaR model) to describe this multimodality therapy. The model is motivated by the observation that heat induced radio-sensitisation may be explained by a reduction in the DNA damage repair capacity of heated cells. We assume that this repair is only possible up to a threshold level above which survival will decrease exponentially with dose. Experimental cell survival data from two cell lines (HCT116, Cal27) were considered along with that taken from the literature (baby hamster kidney [BHK] and Chinese hamster ovary cells [CHO]) for HT and combined RT-HT. The AlphaR model was used to study the dependence of clonogenic survival on treatment temperature, and thermal dose R2 ≥ 0.95 for all fits). For HT survival curves (0-80 CEM43 at 43.5-57 °C), the number of free fit AlphaR model parameters could be reduced to two. Both parameters increased exponentially with temperature. We derived the relative biological effectiveness (RBE) or HT treatments at different temperatures, to provide an alternative description of thermal dose, based on our AlphaR model. For combined RT-HT, our analysis is restricted to the linear quadratic arm of the model. We show that, for the range used (20-80 CEM43, 0-12 Gy), thermal dose is a valid indicator of heat induced radio-sensitisation, and that the model parameters can be described as a function thereof. Overall, the proposed model provides a flexible framework for describing cell survival curves, and may contribute to better quantification of heat induced radio-sensitisation, and thermal dose in general.

Motion monitoring during a course of lung radiotherapy with anchored electromagnetic transponders : Quantification of inter- and intrafraction motion and variability of relative transponder positions.

PURPOSE: Anchored electromagnetic transponders for tumor motion monitoring during lung radiotherapy were clinically evaluated. First, intrafractional motion patterns were analyzed as well as their interfractional variations. Second, intra- and interfractional changes of the geometric transponder positions were investigated. MATERIALS AND METHODS: Intrafractional motion data from 7 patients with an upper or middle lobe tumor and three implanted transponders each was used to calculate breathing amplitudes, overall motion amount and motion midlines in three mutual perpendicular directions and three-dimensionally (3D) for 162 fractions. For 6 patients intra- and interfractional variations in transponder distances and in the size of the triangle defined by the transponder locations over the treatment course were determined. RESULTS: Mean 3D values of all fractions were up to 4.0, 4.6 and 3.4 mm per patient for amplitude, overall motion amount and midline deviation, respectively. Intrafractional transponder distances varied with standard deviations up to 3.2 mm, while a maximal triangle shrinkage of 36.5% over 39 days was observed. CONCLUSIONS: Electromagnetic real-time motion monitoring was feasible for all patients. Detected respiratory motion was on average modest in this small cohort without lower lobe tumors, but changes in motion midline were of the same size as the amplitudes and greater midline motion can be observed in some fractions. Intra- and interfractional variations of the geometric transponder positions can be large, so for reliable motion management correlation between transponder and tumor motion needs to be evaluated per patient.

Physically constrained voxel-based penalty adaptation for ultra-fast IMRT planning.

Conventional treatment planning in intensity-modulated radiation therapy (IMRT) is a trial-and-error process that usually involves tedious tweaking of optimization parameters. Here, we present an algorithm that automates part of this process, in particular the adaptation of voxel-based penalties within normal tissue. Thereby, the proposed algorithm explicitly considers a priori known physical limitations of photon irradiation. The efficacy of the developed algorithm is assessed during treatment planning studies comprising 16 prostate and 5 head and neck cases. We study the eradication of hot spots in the normal tissue, effects on target coverage and target conformity, as well as selected dose volume points for organs at risk. The potential of the proposed method to generate class solutions for the two indications is investigated. Run-times of the algorithms are reported. Physically constrained voxel-based penalty adaptation is an adequate means to automatically detect and eradicate hot-spots during IMRT planning while maintaining target coverage and conformity. Negative effects on organs at risk are comparably small and restricted to lower doses. Using physically constrained voxel-based penalty adaptation, it was possible to improve the generation of class solutions for both indications. Considering the reported run-times of less than 20 s, physically constrained voxel-based penalty adaptation has the potential to reduce the clinical workload during planning and automated treatment plan generation in the long run, facilitating adaptive radiation treatments.

Comparison of CT number calibration techniques for CBCT-based dose calculation.

PURPOSE: The aim of this work was to compare and validate various computed tomography (CT) number calibration techniques with respect to cone beam CT (CBCT) dose calculation accuracy. METHODS: CBCT dose calculation accuracy was assessed for pelvic, lung, and head and neck (H&N) treatment sites for two approaches: (1) physics-based scatter correction methods (CBCTr); (2) density override approaches including assigning water density to the entire CBCT (W), assignment of either water or bone density (WB), and assignment of either water or lung density (WL). Methods for CBCT density assignment within a commercially available treatment planning system (RSauto), where CBCT voxels are binned into six density levels, were assessed and validated. Dose-difference maps and dose-volume statistics were used to compare the CBCT dose distributions with the ground truth of a planning CT acquired the same day as the CBCT. RESULTS: For pelvic cases, all CTN calibration methods resulted in average dose-volume deviations below 1.5 %. RSauto provided larger than average errors for pelvic treatments for patients with large amounts of adipose tissue. For H&N cases, all CTN calibration methods resulted in average dose-volume differences below 1.0 % with CBCTr (0.5 %) and RSauto (0.6 %) performing best. For lung cases, WL and RSauto methods generated dose distributions most similar to the ground truth. CONCLUSION: The RSauto density override approach is an attractive option for CTN adjustments for a variety of anatomical sites. RSauto methods were validated, resulting in dose calculations that were consistent with those calculated on diagnostic-quality CT images, for CBCT images acquired of the lung, for patients receiving pelvic RT in cases without excess adipose tissue, and for H&N cases.

Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy.

BACKGROUND: Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking. PURPOSE: We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions. METHODS: A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation. RESULTS: On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively. CONCLUSIONS: Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Clinical acceptance and dosimetric impact of automatically delineated elective target and organs at risk for head and neck MR-Linac patients.

Background: MR-Linac allows for daily online treatment adaptation to the observed geometry of tumor targets and organs at risk (OARs). Manual delineation for head and neck cancer (HNC) patients takes 45-75 minutes, making it unsuitable for online adaptive radiotherapy. This study aims to clinically and dosimetrically validate an in-house developed algorithm which automatically delineates the elective target volume and OARs for HNC patients in under a minute. Methods: Auto-contours were generated by an in-house model with 2D U-Net architecture trained and tested on 52 MRI scans via leave-one-out cross-validation. A randomized selection of 684 automated and manual contours (split half-and-half) was presented to an oncologist to perform a blind test and determine the clinical acceptability. The dosimetric impact was investigated for 13 patients evaluating the differences in dosage for all structures. Results: Automated contours were generated in 8 seconds per MRI scan. The blind test concluded that 114 (33%) of auto-contours required adjustments with 85 only minor and 15 (4.4%) of manual contours required adjustments with 12 only minor. Dosimetric analysis showed negligible dosimetric differences between clinically acceptable structures and structures requiring minor changes. The Dice Similarity coefficients for the auto-contours ranged from 0.66 ± 0.11 to 0.88 ± 0.06 across all structures. Conclusion: Majority of auto-contours were clinically acceptable and could be used without any adjustments. Majority of structures requiring minor adjustments did not lead to significant dosimetric differences, hence manual adjustments were needed only for structures requiring major changes, which takes no longer than 10 minutes per patient.

Respiratory motion modelling for MR-guided lung cancer radiotherapy: model development and geometric accuracy evaluation.

Objective.Respiratory motion of lung tumours and adjacent structures is challenging for radiotherapy. Online MR-imaging cannot currently provide real-time volumetric information of the moving patient anatomy, therefore limiting precise dose delivery, delivered dose reconstruction, and downstream adaptation methods.Approach.We tailor a respiratory motion modelling framework towards an MR-Linac workflow to estimate the time-resolved 4D motion from real-time data. We develop a multi-slice acquisition scheme which acquires thick, overlapping 2D motion-slices in different locations and orientations, interleaved with 2D surrogate-slices from a fixed location. The framework fits a motion model directly to the input data without the need for sorting or binning to account for inter- and intra-cycle variation of the breathing motion. The framework alternates between model fitting and motion-compensated super-resolution image reconstruction to recover a high-quality motion-free image and a motion model. The fitted model can then estimate the 4D motion from 2D surrogate-slices. The framework is applied to four simulated anthropomorphic datasets and evaluated against known ground truth anatomy and motion. Clinical applicability is demonstrated by applying our framework to eight datasets acquired on an MR-Linac from four lung cancer patients.Main results.The framework accurately reconstructs high-quality motion-compensated 3D images with 2 mm3isotropic voxels. For the simulated case with the largest target motion, the motion model achieved a mean deformation field error of 1.13 mm. For the patient cases residual error registrations estimate the model error to be 1.07 mm (1.64 mm), 0.91 mm (1.32 mm), and 0.88 mm (1.33 mm) in superior-inferior, anterior-posterior, and left-right directions respectively for the building (application) data.Significance.The motion modelling framework estimates the patient motion with high accuracy and accurately reconstructs the anatomy. The image acquisition scheme can be flexibly integrated into an MR-Linac workflow whilst maintaining the capability of online motion-management strategies based on cine imaging such as target tracking and/or gating.

Stereotactic Body Radiotherapy (SBRT) to Localised Prostate Cancer in the Era of MRI-Guided Adaptive Radiotherapy: Doses Delivered in the HERMES Trial Comparing Two- and Five-Fraction Treatments.

HERMES is a phase II trial of MRI-guided daily-adaptive radiotherapy (MRIgART) randomising men with localised prostate cancer to either 2-fractions of SBRT with a boost to the tumour or 5-fraction SBRT. In the context of this highly innovative regime the dose delivered must be carefully considered. The first ten patients recruited to HERMES were analysed in order to establish the dose received by the targets and organs at risk (OARS) in the context of intrafraction motion. A regression analysis was performed to measure how the volume of air within the rectum might further impact rectal dose secondary to the electron return effect (ERE). One hundred percent of CTV target objectives were achieved on the MRI taken prior to beam-on-time. The post-delivery MRI showed that high-dose CTV coverage was achieved in 90% of sub-fractions (each fraction is delivered in two sub-fractions) in the 2-fraction cohort and in 88% of fractions the 5-fraction cohort. Rectal D1 cm3 was the most exceeded constraint; three patients exceeded the D1 cm3 < 20.8 Gy in the 2-fraction cohort and one patient exceeded the D1 cm3 < 36 Gy in the 5-fraction cohort. The volume of rectal gas within 1 cm of the prostate was directly proportional to the increase in rectal D1 cm3, with a strong (R = 0.69) and very strong (R = 0.90) correlation in the 2-fraction and 5-fraction cohort respectively. Dose delivery specified in HERMES is feasible, although for some patients delivered doses to both target and OARs may vary from those planned.

Magnetic resonance image-guided adaptive radiotherapy enables safe CTV-to-PTV margin reduction in prostate cancer: a cine MRI motion study.

Introduction: We aimed to establish if stereotactic body radiotherapy to the prostate can be delivered safely using reduced clinical target volume (CTV) to planning target volume (PTV) margins on the 1.5T MR-Linac (MRL) (Elekta, Stockholm, Sweden), in the absence of gating. Methods: Cine images taken in 3 orthogonal planes during the delivery of prostate SBRT with 36.25 Gray (Gy) in 5 fractions on the MRL were analysed. Using the data from 20 patients, the percentage of radiotherapy (RT) delivery time where the prostate position moved beyond 1, 2, 3, 4 and 5 mm in the left-right (LR), superior-inferior (SI), anterior-posterior (AP) and any direction was calculated. Results: The prostate moved less than 3 mm in any direction for 90% of the monitoring period in 95% of patients. On a per-fraction basis, 93% of fractions displayed motion in all directions within 3 mm for 90% of the fraction delivery time. Recurring motion patterns were observed showing that the prostate moved with shallow drift (most common), transient excursions and persistent excursions during treatment. Conclusion: A 3 mm CTV-PTV margin is safe to use for the treatment of 5 fraction prostate SBRT on the MRL, without gating. In the context of gating this work suggests that treatment time will not be extensively lengthened when an appropriate gating window is applied.

Interim Toxicity Analysis From the Randomized HERMES Trial of 2- and 5-Fraction Magnetic Resonance Imaging-Guided Adaptive Prostate Radiation Therapy.

PURPOSE: Ultrahypofractionated radiation therapy (UHRT) is an effective treatment for localized prostate cancer with an acceptable toxicity profile; boosting the visible intraprostatic tumor has been shown to improve biochemical disease-free survival with no significant effect on genitourinary (GU) and gastrointestinal (GI) toxicity. METHODS AND MATERIALS: HERMES is a single-center noncomparative randomized phase 2 trial in men with intermediate or lower high risk prostate cancer. Patients were allocated (1:1) to 36.25 Gy in 5 fractions over 2 weeks or 24 Gy in 2 fractions over 8 days with an integrated boost to the magnetic resonance imaging (MRI) visible tumor of 27 Gy in 2 fractions. A minimization algorithm with a random element with risk group as a balancing factor was used for participant randomization. Treatment was delivered on the Unity MR-Linac (Elekta AB) with daily online adaption. The primary endpoint was acute GU Common Terminology Criteria for Adverse Events version 5.0 toxicity with the aim of excluding a doubling of the rate of acute grade 2+ GU toxicity seen in PACE. Analysis was by treatment received and included all participants who received at least 1 fraction of study treatment. This interim analysis was prespecified (stage 1 of a 2-stage Simon design) for when 10 participants in each treatment group had completed the acute toxicity monitoring period (12 weeks after radiation therapy). RESULTS: Acute grade 2 GU toxicity was reported in 1 (10%) patient in the 5-fraction group and 2 (20%) patients in the 2-fraction group. No grade 3+ GU toxicities were reported. CONCLUSIONS: At this interim analysis, the rate of GU toxicity in the 2-fraction and 5-fraction treatment groups was found to be below the prespecified threshold (5/10 grade 2+) and continuation of the study to complete recruitment of 23 participants per group was recommended.

Accelerating 4D image reconstruction for magnetic resonance-guided radiotherapy.

Background and purpose: Physiological motion impacts the dose delivered to tumours and vital organs in external beam radiotherapy and particularly in particle therapy. The excellent soft-tissue demarcation of 4D magnetic resonance imaging (4D-MRI) could inform on intra-fractional motion, but long image reconstruction times hinder its use in online treatment adaptation. Here we employ techniques from high-performance computing to reduce 4D-MRI reconstruction times below two minutes to facilitate their use in MR-guided radiotherapy. Material and methods: Four patients with pancreatic adenocarcinoma were scanned with a radial stack-of-stars gradient echo sequence on a 1.5T MR-Linac. Fast parallelised open-source implementations of the extra-dimensional golden-angle radial sparse parallel algorithm were developed for central processing unit (CPU) and graphics processing unit (GPU) architectures. We assessed the impact of architecture, oversampling and respiratory binning strategy on 4D-MRI reconstruction time and compared images using the structural similarity (SSIM) index against a MATLAB reference implementation. Scaling and bottlenecks for the different architectures were studied using multi-GPU systems. Results: All reconstructed 4D-MRI were identical to the reference implementation (SSIM > 0.99). Images reconstructed with overlapping respiratory bins were sharper at the cost of longer reconstruction times. The CPU  + GPU implementation was over 17 times faster than the reference implementation, reconstructing images in 60 ± 1 s and hyper-scaled using multiple GPUs. Conclusion: Respiratory-resolved 4D-MRI reconstruction times can be reduced using high-performance computing methods for online workflows in MR-guided radiotherapy with potential applications in particle therapy.

Evaluation of non-vendor magnetic resonance imaging sequences for use in bladder cancer magnetic resonance image guided radiotherapy.

Hybrid systems that combine Magnetic Resonance Imaging (MRI) and linear accelerators are available clinically to guide and adapt radiotherapy. Vendor-approved MRI sequences are provided, however alternative sequences may offer advantages. The aim of this study was to develop a systematic approach for non-vendor sequence evaluation, to determine safety, accuracy and overall clinical application of two potential sequences for bladder cancer MRI guided radiotherapy. Non-vendor sequences underwent and passed clinical image qualitative review, phantom quality assurance, and radiotherapy planning assessments. Volunteer workflow tests showed the potential for one sequence to reduce workflow time by 27% compared to the standard vendor sequence.

Joint radial trajectory correction for accelerated T2 * mapping on an MR-Linac.

BACKGROUND: T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes. PURPOSE: The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac). MATERIALS AND METHODS: The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared. RESULTS: Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV). CONCLUSION: Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac.

Position detection accuracy of a novel linac-mounted intrafractional x-ray imaging system.

PURPOSE: The authors have developed a system that monitors intrafractional target motion perpendicular to the treatment beam with the aid of radioopaque markers by means of separating kV image and megavoltage (MV) treatment field on a single flat-panel detector. METHODS: They equipped a research Siemens Artiste linear accelerator (linac) with a 41 × 41 cm(2) a-Si flat-panel detector underneath the treatment head. The in-line geometry allows kV (imaging) and MV (treatment) beams to share closely aligned beam axes. The kV source, usually mounted directly across from the flat-panel imager, was retracted toward the gantry by 13 cm to intentionally misalign kV and MV beams, resulting in a geometric separation of MV treatment field and kV image on the detector. Two consecutive images acquired within 140 ms (the first with MV-only and the second with kV and MV signal) were subtracted to generate a kV-only image. The images were then analyzed "online" with an automated threshold-based marker detection algorithm. They employed a 3D and a 4D phantom equipped with either a single radioopaque marker or three Calypso beacons to mimic respiratory motion. Measured room positions were either cross-referenced with a phantom voltage signal (single marker) or the Calypso system. The accuracy of the back-projection (from detected marker positions into room coordinates) was verified by a simulation study. RESULTS: A phantom study has demonstrated that the imaging framework is capable of automatically detecting marker positions and sending this information to the tracking tool at an update rate of 7.14 Hz. The system latency is 86.9 ± 1.0 ms for single marker detection in the absence of MV radiation. In the presence of a circular MV field of 5 cm diameter, the latency is 87.1 ± 0.9 ms. The total RMS position detection accuracy is 0.20 mm (without MV radiation) and 0.23 mm (with MV). CONCLUSIONS: Based on the evaluated motion patterns and MV field size, the positional accuracy and system latency indicate that this system is suitable for real-time adaptive applications.

Comparison of a multileaf collimator tracking system and a robotic treatment couch tracking system for organ motion compensation during radiotherapy.

PURPOSE: One limitation of accurate dose delivery in radiotherapy is intrafractional movement of the tumor or the entire patient which may lead to an underdosage of the target tissue or an overdosage of adjacent organs at risk. In order to compensate for this movement, different techniques have been developed. In this study the tracking performances of a multileaf collimator (MLC) tracking system and a robotic treatment couch tracking system were compared under equal conditions. METHODS: MLC tracking was performed using a tracking system based on the Siemens 160 MLC. A HexaPOD robotic treatment couch tracking system was also installed at the same linac. A programmable 4D motion stage was used to reproduce motion trajectories with different target phantoms. Motion localization of the target was provided by the 4D tracking system of Calypso Medical Inc. The gained positional data served as input signal for the control systems of the MLC and HexaPOD tracking systems attempting to compensate for the target motion. The geometric and dosimetric accuracy for the tracking of eight different respiratory motion trajectories was investigated for both systems. The dosimetric accuracy of both systems was also evaluated for the tracking of five prostate motion trajectories. RESULTS: For the respiratory motion the average root mean square error of all trajectories in y direction was reduced from 4.1 to 2.0 mm for MLC tracking and to 2.2 mm for HexaPOD tracking. In x direction it was reduced from 1.9 to 0.9 mm (MLC) and to 1.0 mm (HexaPOD). The average 2%/2 mm gamma pass rate for the respiratory motion trajectories was increased from 76.4% for no tracking to 89.8% and 95.3% for the MLC and the HexaPOD tracking systems, respectively. For the prostate motion trajectories the average 2%/2 mm gamma pass rate was 60.1% when no tracking was applied and was improved to 85.0% for MLC tracking and 95.3% for the HexaPOD tracking system. CONCLUSIONS: Both systems clearly increased the geometric and dosimetric accuracy during tracking of respiratory motion trajectories. Thereby, the geometric accuracy was increased almost equally by both systems, whereas the dosimetric accuracy of the HexaPOD tracking system was slightly better for all considered respiratory motion trajectories. Substantial improvement of the dosimetric accuracy was also observed during tracking of prostate motion trajectories during an intensity-modulated radiotherapy plan. Thereby, the HexaPOD tracking system showed better results than the MLC tracking.

Online adaptive radiotherapy for head and neck cancers on the MR linear Accelerator: Introducing a novel modified Adapt-to-Shape approach.

INTRODUCTION: The Elekta Unity MR-Linac (MRL) has enabled adaptive radiotherapy (ART) for patients with head and neck cancers (HNC). Adapt-To-Shape-Lite (ATS-Lite) is a novel Adapt-to-Shape strategy that provides ART without requiring daily clinician presence to perform online target and organ at risk (OAR) delineation. In this study we compared the performance of our clinically-delivered ATS-Lite strategy against three Adapt-To-Position (ATP) variants: Adapt Segments (ATP-AS), Optimise Weights (ATP-OW), and Optimise Shapes (ATP-OS). METHODS: Two patients with HNC received radical-dose radiotherapy on the MRL. For each fraction, an ATS-Lite plan was generated online and delivered and additional plans were generated offline for each ATP variant. To assess the clinical acceptability of a plan for every fraction, twenty clinical goals for targets and OARs were assessed for all four plans. RESULTS: 53 fractions were analysed. ATS-Lite passed 99.9% of mandatory dose constraints. ATP-AS and ATP-OW each failed 7.6% of mandatory dose constraints. The Planning Target Volumes for 54 Gy (D95% and D98%) were the most frequently failing dose constraint targets for ATP. ATS-Lite median fraction times for Patient 1 and 2 were 40 mins 9 s (range 28 mins 16 s - 47 mins 20 s) and 32 mins 14 s (range 25 mins 33 s - 44 mins 27 s), respectively. CONCLUSIONS: Our early data show that the novel ATS-Lite strategy produced plans that fulfilled 99.9% of clinical dose constraints in a time frame that is tolerable for patients and comparable to ATP workflows. Therefore, ATS-Lite, which bridges the gap between ATP and full ATS, will be further utilised and developed within our institute and it is a workflow that should be considered for treating patients with HNC on the MRL.