Purpose

To develop a technique that assesses the accuracy of the breathing phase-specific volume image generation process by patient-specific breathing motion model using the original free-breathing computed tomographic (CT) scans as ground truths.

Methods

Sixteen lung cancer patients underwent a previously published protocol in which 25 free-breathing fast helical CT scans were acquired with a simultaneous breathing surrogate. A patient-specific motion model was constructed based on the tissue displacements determined by a state-of-the-art deformable image registration. The first image was arbitrarily selected as the reference image. The motion model was used, along with the free-breathing phase information of the original 25 image datasets, to generate a set of deformation vector fields that mapped the reference image to the 24 nonreference images. The high-pitch helically acquired original scans served as ground truths because they captured the instantaneous tissue positions during free breathing. Image similarity between the simulated and the original scans was assessed using deformable registration that evaluated the pointwise discordance throughout the lungs.

Results

Qualitative comparisons using image overlays showed excellent agreement between the simulated images and the original images. Even large 2-cm diaphragm displacements were very well modeled, as was sliding motion across the lung-chest wall boundary. The mean error across the patient cohort was 1.15 ± 0.37 mm, and the mean 95th percentile error was 2.47 ± 0.78 mm.

Conclusion

The proposed ground truth-based technique provided voxel-by-voxel accuracy analysis that could identify organ-specific or tumor-specific motion modeling errors for treatment planning. Despite a large variety of breathing patterns and lung deformations during the free-breathing scanning session, the 5-dimensionl CT technique was able to accurately reproduce the original helical CT scans, suggesting its applicability to a wide range of patients.

Purpose

To develop and validate a technique for radiation therapy gating using slow (≤1 frame per second) magnetic resonance imaging (MRI) and a motion model. Proposed uses of the technique include radiation therapy gating using T2-weighted images and conducting additional imaging studies during gated treatments.

Methods and materials

The technique uses a physiologically guided breathing motion model to interpolate deformed target position between 2-dimensional (2D) MRI images acquired every 1 to 3 seconds. The model is parameterized by a 1-dimensional respiratory bellows surrogate and is continuously updated with the most recently acquired 2D images. A phantom and 8 volunteers were imaged with a 0.35T MRI-guided radiation therapy system. A balanced steady-state free precession sequence with a 2D frame rate of 3 frames per second was used to evaluate the technique. The accuracy and beam-on positive predictive value (PPV) of the model-based gating decisions were evaluated using the gating decisions derived from imaging as a ground truth. A T2-weighted gating offline proof-of-concept study using a half-Fourier, single-shot, turbo-spin echo sequence is reported.

Results

Model-interpolated gating accuracy, beam-on PPV, and median absolute distances between model and image-tracked target centroids were, on average, 98.3%, 98.4%, and 0.33 mm, respectively, in the balanced steady-state free precession phantom studies and 93.7%, 92.1%, and 0.86 mm, respectively, in the volunteer studies. T2 model-interpolated gating in 6 volunteers yielded an average accuracy and PPV of 94.3% and 92.5%, respectively, and the mean absolute median distance between modeled and imaged target centroids was 0.86 mm.

Conclusions

This work demonstrates the concept of model-interpolated gating for MRI-guided radiation therapy. The technique was found to be potentially sufficiently accurate for clinical use. Further development is needed to accommodate out-of-plane motion and the use of an internal MR-based respiratory surrogate.

Purpose

To determine if and by how much the commercial 4DCT protocols under- and overestimate tumor breathing motion.

Methods

1D simulations were conducted that modeled a 16-slice CT scanner and tumors moving proportionally to breathing amplitude. External breathing surrogate traces of at least 5-min duration for 50 patients were used. Breathing trace amplitudes were converted to motion by relating the nominal tumor motion to the 90th percentile breathing amplitude, reflecting motion defined by the more recent 5DCT approach. Based on clinical low-pitch helical CT acquisition, the CT detector moved according to its velocity while the tumor moved according to the breathing trace. When the CT scanner overlapped the tumor, the overlapping slices were identified as having imaged the tumor. This process was repeated starting at successive 0.1 s time bin in the breathing trace until there was insufficient breathing trace to complete the simulation. The tumor size was subtracted from the distance between the most superior and inferior tumor positions to determine the measured tumor motion for that specific simulation. The effect of the scanning parameter variation was evaluated using two commercial 4DCT protocols with different pitch values. Because clinical 4DCT scan sessions would yield a single tumor motion displacement measurement for each patient, errors in the tumor motion measurement were considered systematic. The mean of largest 5% and smallest 5% of the measured motions was selected to identify over- and underdetermined motion amplitudes, respectively. The process was repeated for tumor motions of 1-4 cm in 1 cm increments and for tumor sizes of 1-4 cm in 1 cm increments.

Results

In the examined patient cohort, simulation using pitch of 0.06 showed that 30% of the patients exhibited a 5% chance of mean breathing amplitude overestimations of 47%, while 30% showed a 5% chance of mean breathing amplitude underestimations of 36%; with a separate simulation using pitch of 0.1 showing, respectively, 37% overestimation and 61% underestimation.

Conclusions

The simulation indicates that commercial low-pitch helical 4DCT processes potentially yield large tumor motion measurement errors, both over- and underestimating the tumor motion.

Objective

A recently proposed "5DCT" protocol uses deformable registration of free-breathing fast-helical CT scans to generate a breathing motion model. In order to allow accurate registration, free-breathing images are required to be free of doubling-artefacts, which arise when tissue motion is greater than scan speed.

Methods

Using a unique set of digital phantoms based on patient data and verified with a motion phantom, this work identifies the minimum scanner parameters required to successfully generate free-breathing artefact-free fast-helical scans. A motion phantom and 5 patients were imaged 25 times under free-breathing conditions in alternating directions with a 64-slice CT scanner employing a low-dose fast-helical protocol. A series of high temporal resolution (0.1 s) 5DCT scan data sets was generated in each case. A simulated CT scanner was used to "image" each free-breathing data set. Various CT scanner detector widths and rotation times were simulated, and verified using the motion phantom results. Motion-induced artefacts were quantified in patient images using structural similarity maps to determine the similarity between axial slices.

Results

Increasing amounts of motion-induced artefacts were observed with increasing rotation times >0.2 s for 16 mm detector configuration.

Conclusion

The current generation of 16-slice CT scanners, which are present in the majority of Radiation Oncology departments, are not capable of generating free-breathing sorting artefact-free images required for 5DCT. Advances in knowledge: A recently proposed "5DCT" protocol uses deformable registration of free-breathing fast-helical CT scans to generate a breathing motion model. In order to allow accurate registration, free-breathing images are required to be free of doubling-artefacts, which arise when tissue motion is greater than scan speed. The results suggest that the current generation of 16-slice CT scanners, present in the majority of Radiation Oncology departments, are not capable of generating the free-breathing images required for 5DCT.

Purpose

Magnetic resonance image (MRI) guided radiotherapy enables gating directly on the target position. We present an evaluation of an MRI-guided radiotherapy system's gating performance using an MRI-compatible respiratory motion phantom and radiochromic film. Our evaluation is geared toward validation of our institution's clinical gating protocol which involves planning to a target volume formed by expanding 5 mm about the gross tumor volume (GTV) and gating based on a 3 mm window about the GTV.

Methods

The motion phantom consisted of a target rod containing high-contrast target inserts which moved in the superior-inferior direction inside a body structure containing background contrast material. The target rod was equipped with a radiochromic film insert. Treatment plans were generated for a 3 cm diameter spherical planning target volume, and delivered to the phantom at rest and in motion with and without gating. Both sinusoidal trajectories and tumor trajectories measured during MRI-guided treatments were used. Similarity of the gated dose distribution to the planned, motion-frozen, distribution was quantified using the gamma technique.

Results

Without gating, gamma pass rates using 4%/3 mm criteria were 22-59% depending on motion trajectory. Using our clinical standard of repeated breath holds and a gating window of 3 mm with 10% target allowed outside the gating boundary, the gamma pass rate was 97.8% with 3%/3 mm gamma criteria. Using a 3 mm window and 10% allowed excursion, all of the patient tumor motion trajectories at actual speed resulting in at least 95% gamma pass rate at 4%/3 mm.

Conclusions

Our results suggest that the device can be used to compensate respiratory motion using a 3 mm gating margin and 10% allowed excursion results in conjunction with repeated breath holds. Full clinical validation requires a comprehensive evaluation of tracking performance in actual patient images, outside the scope of this study.

Purpose

To quantitatively compare the accuracy of tumor volume segmentation in amplitude-based and phase-based respiratory gating algorithms in respiratory-correlated positron emission tomography (PET).

Methods and materials

List-mode fluorodeoxyglucose-PET data was acquired for 10 patients with a total of 12 fluorodeoxyglucose-avid tumors and 9 lymph nodes. Additionally, a phantom experiment was performed in which 4 plastic butyrate spheres with inner diameters ranging from 1 to 4 cm were imaged as they underwent 1-dimensional motion based on 2 measured patient breathing trajectories. PET list-mode data were gated into 8 bins using 2 amplitude-based (equal amplitude bins [A1] and equal counts per bin [A2]) and 2 temporal phase-based gating algorithms. Gated images were segmented using a commercially available gradient-based technique and a fixed 40% threshold of maximum uptake. Internal target volumes (ITVs) were generated by taking the union of all 8 contours per gated image. Segmented phantom ITVs were compared with their respective ground-truth ITVs, defined as the volume subtended by the tumor model positions covering 99% of breathing amplitude. Superior-inferior distances between sphere centroids in the end-inhale and end-exhale phases were also calculated.

Results

Tumor ITVs from amplitude-based methods were significantly larger than those from temporal-based techniques (P=.002). For lymph nodes, A2 resulted in ITVs that were significantly larger than either of the temporal-based techniques (P<.0323). A1 produced the largest and most accurate ITVs for spheres with diameters of ≥2 cm (P=.002). No significant difference was shown between algorithms in the 1-cm sphere data set. For phantom spheres, amplitude-based methods recovered an average of 9.5% more motion displacement than temporal-based methods under regular breathing conditions and an average of 45.7% more in the presence of baseline drift (P<.001).

Conclusions

Target volumes in images generated from amplitude-based gating are larger and more accurate, at levels that are potentially clinically significant, compared with those from temporal phase-based gating.

OBJECTIVE:To evaluate variations in intra- and interfractional tumour motion, and the effect on internal target volume (ITV) contour accuracy, using deformable image registration of real-time two-dimensional-sagittal cine-mode MRI acquired during lung stereotactic body radiation therapy (SBRT) treatments. METHODS:Five lung tumour patients underwent free-breathing SBRT treatments on the ViewRay system, with dose prescribed to a planning target volume (defined as a 3-6 mm expansion of the 4DCT-ITV). Sagittal slice cine-MR images (3.5 × 3.5 mm2 pixels) were acquired through the centre of the tumour at 4 frames per second throughout the treatments (3-4 fractions of 21-32 min). Tumour gross tumour volumes (GTVs) were contoured on the first frame of the MR cine and tracked for the first 20 min of each treatment using offline optical-flow based deformable registration implemented on a GPU cluster. A ground truth ITV (MR-ITV20 min) was formed by taking the union of tracked GTV contours. Pseudo-ITVs were generated from unions of the GTV contours tracked over 10 s segments of image data (MR-ITV10 s). RESULTS:Differences were observed in the magnitude of median tumour displacement between days of treatments. MR-ITV10 s areas were as small as 46% of the MR-ITV20 min. CONCLUSION:An ITV offers a "snapshot" of breathing motion for the brief period of time the tumour is imaged on a specific day. Real-time MRI over prolonged periods of time and over multiple treatment fractions shows that ITV size varies. Further work is required to investigate the dosimetric effect of these results. Advances in knowledge: Five lung tumour patients underwent free-breathing MRI-guided SBRT treatments, and their tumours tracked using deformable registration of cine-mode MRI. The results indicate that variability of both intra- and interfractional breathing amplitude should be taken into account during planning of lung radiotherapy.

Purpose

To demonstrate that a "5DCT" technique which utilizes fast helical acquisition yields the same respiratory-gated images as a commercial technique for regular, mechanically produced breathing cycles.

Methods

Respiratory-gated images of an anesthetized, mechanically ventilated pig were generated using a Siemens low-pitch helical protocol and 5DCT for a range of breathing rates and amplitudes and with standard and low dose imaging protocols. 5DCT reconstructions were independently evaluated by measuring the distances between tissue positions predicted by a 5D motion model and those measured using deformable registration, as well by reconstructing the originally acquired scans. Discrepancies between the 5DCT and commercial reconstructions were measured using landmark correspondences.

Results

The mean distance between model predicted tissue positions and deformably registered tissue positions over the nine datasets was 0.65 ± 0.28 mm. Reconstructions of the original scans were on average accurate to 0.78 ± 0.57 mm. Mean landmark displacement between the commercial and 5DCT images was 1.76 ± 1.25 mm while the maximum lung tissue motion over the breathing cycle had a mean value of 27.2 ± 4.6 mm. An image composed of the average of 30 deformably registered images acquired with a low dose protocol had 6 HU image noise (single standard deviation) in the heart versus 31 HU for the commercial images.

Conclusions

An end to end evaluation of the 5DCT technique was conducted through landmark based comparison to breathing gated images acquired with a commercial protocol under highly regular ventilation. The techniques were found to agree to within 2 mm for most respiratory phases and most points in the lung.

Onboard magnetic resonance imaging (MRI) guided radiotherapy is now clinically available in nine centers in the world. This technology has facilitated the clinical implementation of online adaptive radiotherapy (OART), or the ability to alter the daily treatment plan based on tumor and anatomical changes in real-time while the patient is on the treatment table. However, due to the time sensitive nature of OART, implementation in a large and busy clinic has many potential obstacles as well as patient-related safety considerations. In this work, we have described the implementation of this new process of care in the Department of Radiation Oncology at the University of California, Los Angeles (UCLA). We describe the rationale, the initial challenges such as treatment time considerations, technical issues during the process of re-contouring, re-optimization, quality assurance, as well as our current solutions to overcome these challenges. In addition, we describe the implementation of a coverage system with a physician of the day as well as online planners (physicists or dosimetrists) to oversee each OART treatment with patient-specific 'hand-off' directives from the patient's treating physician. The purpose of this effort is to streamline the process without compromising treatment quality and patient safety. As more MRI-guided radiotherapy programs come online, we hope that our experience can facilitate successful adoption of OART in a way that maximally benefits the patient.