Deformable vector fields warping for modelling of irregular breathing.

Purpose 4D computed tomography (4DCT) is the clinical standard to image organ motion in radiotherapy, although it is limited in imaging breathing variability. We propose a method to transfer breathing motion across longitudinal imaging datasets to include intra-patient variability and verify its performance in lung cancer patients. Methods Five repeated control 4DCTs for 6 non-small cell lung cancer patients were combined into multi-breath datasets (m4DCT) by merging stages of deformable image registration to isolate respiratory motion. The displacement of the centre of mass of the primary tumour and its volume changes were evaluated to quantify intra-patient differences. Internal target volumes defined on the m4DCT were compared with those conventionally drawn on the 4DCT. Results Motion analysis suggests no discontinuity at the junction between successive breaths, confirming the method's ability to merge repeated imaging into a continuum. Motion (variability) is primarily in superior-inferior direction and goes from 14.4 mm (8.7 mm) down to 0.1 mm (0.6 mm), respectively for tumours located in the lower lobes or most apical ones. On average, up to 65% and 74% of the tumour volume was subject to expansion or contraction in the inhalation and exhalation phases. These variations lead to an enlargement of the ITV up to 8% of its volume in our dataset. Conclusion 4DCT can be extended to model variable breathing motion by adding synthetic phases from multiple time-resolved images. The inclusion of this improved knowledge of patients' breathing allows better definition of treatment volumes and their margins for radiation therapy. .

Retrospective reconstruction of four-dimensional magnetic resonance from interleaved cine imaging - A comparative study with four-dimensional computed tomography in the lung.

Background and purpose: Imaging of respiration-induced anatomical changes is essential to ensure high accuracy in radiotherapy of lung cancer. We expanded here on methods for retrospective reconstruction of time-resolved volumetric magnetic resonance (4DMR) of the thoracic region and benchmarked the results against 4D computed tomography (4DCT). Materials and method: MR data of six lung cancer patients were collected by interleaving cine-navigator images with 2D data frame images, acquired across the thorax. The data frame images have been stacked in volumes based on a similarity metric that considers the anatomical deformation of lungs, while addressing ambiguities in respiratory phase detection and interpolation of missing data. The resulting images were validated against cine-navigator images and compared to paired 4DCTs in terms of amplitude and period of motion, assessing differences in internal target volume (ITV) margin definition. Results: 4DMR-based motion amplitude was on average within 1.8 mm of that measured in the corresponding 2D cine-navigator images. In our dataset, the 4DCT motion and the 4DMR median amplitude were always within 3.8 mm. The median period was generally close to CT references, although deviations up to 24 % have been observed. These changes were reflected in the ITV, which was generally larger for MRI than for 4DCT (up to 39.7 %). Conclusions: The proposed algorithm for retrospective reconstruction of time-resolved volumetric MR provided quality anatomical images with high temporal resolution for motion modelling and treatment planning. The potential for imaging organ motion variability makes 4DMR a valuable complement to standard 4DCT imaging.

Exploring beamline momentum acceptance for tracking respiratory variability in lung cancer proton therapy: a simulation study.

Objective. Investigating the aspects of proton beam delivery to track organ motion with pencil beam scanning therapy. Considering current systems as a reference, specify requirements for next-generation units aiming at real-time image-guided treatments.Approach. Proton treatments for six non-small cell lung cancer (NSCLC) patients were simulated using repeated 4DCTs to model respiratory motion variability. Energy corrections required for this treatment site were evaluated for different approaches to tumour tracking, focusing on the potential for energy adjustment within beamline momentum acceptance (dp/p). A respiration-synchronised tracking, taking into account realistic machine delivery limits, was compared to ideal tracking scenarios, in which unconstrained energy corrections are possible. Rescanning and the use of multiple fields to mitigate residual interplay effects and dose degradation have also been investigated.Main results. Energy correction requirements increased with motion amplitudes, for all patients and tracking scenarios. Higher dose degradation was found for larger motion amplitudes, rescanning has beneficial effects and helped to improve dosimetry metrics for the investigated limited dp/pof 1.2% (realistic) and 2.4%. The median differences between ideal and respiratory-synchronised tracking show minimal discrepancies, 1% and 5% respectively for dose coverage (CTV V95) and homogeneity (D5-D95). Multiple-field planning improves D5-D95 up to 50% in the most extreme cases while it does not show a significant effect on V95.Significance. This work shows the potential of implementing tumour tracking in current proton therapy units and outlines design requirements for future developments. Energy regulation within momentum acceptance was investigated to tracking tumour motion with respiratory-synchronisation, achieving results in line with the performance of ideal tracking scenarios. ±5% Δp/p would allow to compensate for all range offsets in our NSCLC patient cohort, including breathing variability. However, the realistic momentum of 1.2% dp/prepresentative of existing medical units limitations, has been shown to preserve plan quality.

The impact of organ motion and the appliance of mitigation strategies on the effectiveness of hypoxia-guided proton therapy for non-small cell lung cancer.

BACKGROUND AND PURPOSE: To investigate the impact of organ motion on hypoxia-guided proton therapy treatments for non-small cell lung cancer (NSCLC) patients. MATERIALS AND METHODS: Hypoxia PET and 4D imaging data of six NSCLC patients were used to simulate hypoxia-guided proton therapy with different motion mitigation strategies including rescanning, breath-hold, respiratory gating and tumour tracking. Motion-induced dose degradation was estimated for treatment plans with dose painting of hypoxic tumour sub-volumes at escalated dose levels. Tumour control probability (TCP) and dosimetry indices were assessed to weigh the clinical benefit of dose escalation and motion mitigation. In addition, the difference in normal tissue complication probability (NTCP) between escalated proton and photon VMAT treatments has been assessed. RESULTS: Motion-induced dose degradation was found for target coverage (CTV V95% up to -4%) and quality of the dose-escalation-by-contour (QRMS up to 6%) as a function of motion amplitude and amount of dose escalation. The TCP benefit coming from dose escalation (+4-13%) outweighs the motion-induced losses (<2%). Significant average NTCP reductions of dose-escalated proton plans were found for lungs (-14%), oesophagus (-10%) and heart (-16%) compared to conventional VMAT plans. The best plan dosimetry was obtained with breath hold and respiratory gating with rescanning. CONCLUSION: NSCLC affected by hypoxia appears to be a prime target for proton therapy which, by dose-escalation, allows to mitigate hypoxia-induced radio-resistance despite the sensitivity to organ motion. Furthermore, substantial reduction in normal tissue toxicity can be expected compared to conventional VMAT. Accessibility and standardization of hypoxia imaging and clinical trials are necessary to confirm these findings in a clinical setting.

Increase of the transmission and emittance acceptance through a cyclotron-based proton therapy gantry.

PURPOSE: In proton therapy, the gantry, as the final part of the beamline, has a major effect on beam intensity and beam size at the isocenter. Most of the conventional beam optics of cyclotron-based proton gantries have been designed with an imaging factor between 1 and 2 from the coupling point (CP) at the gantry entrance to the isocenter (patient location) meaning that to achieve a clinically desirable (small) beam size at isocenter, a small beam size is also required at the CP. Here we will show that such imaging factors are limiting the emittance which can be transported through the gantry. We, therefore, propose the use of large beam size and low divergence beam at the CP along with an imaging factor of 0.5 (2:1) in a new design of gantry beam optics to achieve substantial improvements in transmission and thus increase beam intensity at the isocenter. METHODS: The beam optics of our gantry have been re-designed to transport higher emittance without the need of any mechanical modifications to the gantry beamline. The beam optics has been designed using TRANSPORT, with the resulting transmissions being calculated using Monte Carlo simulations (BDSIM code). Finally, the new beam optics have been tested with measurements performed on our Gantry 2 at PSI. RESULTS: With the new beam optics, we could maximize transmission through the gantry for a fixed emittance value. Additionally, we could transport almost four times higher emittance through the gantry compared to conventional optics, whilst achieving good transmissions through the gantry (>50%) with no increased losses in the gantry. As such, the overall transmission (cyclotron to isocenter) can be increased by almost a factor of 6 for low energies. Additionally, the point-to-point imaging inherent to the optics allows adjustment of the beam size at the isocenter by simply changing the beam size at the CP. CONCLUSION: We have developed a new gantry beam optics which, by selecting a large beam size and low divergence at the gantry entrance and using an imaging factor of 0.5 (2:1), increases the emittance acceptance of the gantry, leading to a substantial increase in beam intensity at low energies. We expect that this approach could easily be adapted for most types of existing gantries.

Beam properties within the momentum acceptance of a clinical gantry beamline for proton therapy.

PURPOSE: Energy changes in pencil beam scanning proton therapy can be a limiting factor in delivery time, hence, limiting patient throughput and the effectiveness of motion mitigation techniques requiring fast irradiation. In this study, we investigate the feasibility of performing fast and continuous energy modulation within the momentum acceptance of a clinical beamline for proton therapy. METHODS: The alternative use of a local beam degrader at the gantry coupling point has been compared with a more common upstream regulation. Focusing on clinically relevant parameters, a complete beam properties characterization has been carried out. In particular, the acquired empirical data allowed to model and parametrize the errors in range and beam current to deliver clinical treatment plans. RESULTS: For both options, the local and upstream degrader, depth-dose curves measured in water for off-momentum beams were only marginally distorted (γ(1%, 1 mm) > 90%) and the errors in the spot position were within the clinical tolerance, even though increasing at the boundaries of the investigated scan range. The impact on the beam size was limited for the upstream degrader, while dedicated strategies could be required to tackle the beam broadening through the local degrader. Range correction models were investigated for the upstream regulation. The impaired beam transport required a dedicated strategy for fine range control and compensation of beam intensity losses. Our current parameterization based on empirical data allowed energy modulation within acceptance with range errors (median 0.05 mm) and transmission (median -14%) compatible with clinical operation and remarkably low average 27 ms dead time for small energy changes. The technique, tested for the delivery of a skull glioma treatment, resulted in high gamma pass rates at 1%, 1 mm compared to conventional deliveries in experimental measurements with about 45% reduction of the energy switching time when regulation could be performed within acceptance. CONCLUSIONS: Fast energy modulation within beamline acceptance has potential for clinical applications and, when realized with an upstream degrader, does not require modification in the beamline hardware, therefore, being potentially applicable in any running facility. Centers with slow energy switching time can particularly profit from such a technique for reducing dead time during treatment delivery.