New full-counts phase-matched data-driven gated (DDG) PET/CT

BACKGROUND

Data-driven gated (DDG) PET has gained clinical acceptance and has been shown to match or outperform external-device gated (EDG) PET. However, in most clinical applications, DDG PET is matched with helical CT acquired in free breathing (FB) at a random respiratory phase, leaving registration, and optimal attenuation correction (AC) to chance. Furthermore, DDG PET requires additional scan time to reduce image noise as it only preserves 35%-50% of the PET data at or near the end-expiratory phase of the breathing cycle.

PURPOSE

A new full-counts, phase-matched (FCPM) DDG PET/CT was developed based on a low-dose cine CT to improve registration between DDG PET and DDG CT, to reduce image noise, and to avoid increasing acquisition times in DDG PET.

METHODS

A new DDG CT was developed for three respiratory phases of CT images from a low dose cine CT acquisition of 1.35 mSv for a coverage of about 15.4 cm: end-inspiration (EI), average (AVG), and end-expiration (EE) to match with the three corresponding phases of DDG PET data: -10% to 15%; 15% to 30%, and 80% to 90%; and 30% to 80%, respectively. The EI and EE phases of DDG CT were selected based on the physiological changes in lung density and body outlines reflected in the dynamic cine CT images. The AVG phase was derived from averaging of all phases of the cine CT images. The cine CT was acquired over the lower lungs and/or upper abdomen for correction of misregistration between PET and FB CT as well as DDG PET and FB CT. The three phases of DDG CT were used for AC of the corresponding phases of PET. After phase-matched AC of each PET dataset, the EI and AVG PET data were registered to the EE PET data with deformable image registration. The final result was FCPM DDG PET/CT which accounts for all PET data registered at the EE phase. We applied this approach to 14 18F-FDG lung cancer patient studies acquired at 2 min/bed position on the GE Discovery MI (25-cm axial FOV) and evaluated its efficacy in improved quantification and noise reduction.

RESULTS

Relative to static PET/CT, the SUVmax increases for the EI, AVG, EE, and FCPM DDG PET/CT were 1.67 ± 0.40, 1.50 ± 0.28, 1.64 ± 0.36, and 1.49 ± 0.28, respectively. There were 10.8% and 9.1% average decreases in SUVmax from EI and EE to FCPM DDG PET/CT, respectively. EI, AVG, and EE DDG PET/CT all maintained increased image noise relative to static PET/CT. However, the noise levels of FCPM and static PET were statistically equivalent, suggesting the inclusion of all counts was able to decrease the image noise relative to EI and EE DDG PET/CT.

CONCLUSIONS

A new FCPM DDG PET/CT has been developed to account for 100% of collected PET data in DDG PET applications. Image noise in FCPM is comparable to static PET, while small decreases in SUVmax were also observed in FCPM when compared to either EI or EE DDG PET/CT.

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

BACKGROUND

Volumetric reconstruction of magnetic resonance imaging (MRI) from sparse samples is desirable for 3D motion tracking and promises to improve magnetic resonance (MR)-guided radiation treatment precision. Data-driven sparse MRI reconstruction, however, requires large-scale training datasets for prior learning, which is time-consuming and challenging to acquire in clinical settings.

PURPOSE

To investigate volumetric reconstruction of MRI from sparse samples of two orthogonal slices aided by sparse priors of two static 3D MRI through implicit neural representation (NeRP) learning, in support of 3D motion tracking during MR-guided radiotherapy.

METHODS

A multi-layer perceptron network was trained to parameterize the NeRP model of a patient-specific MRI dataset, where the network takes 4D data coordinates of voxel locations and motion states as inputs and outputs corresponding voxel intensities. By first training the network to learn the NeRP of two static 3D MRI with different breathing motion states, prior information of patient breathing motion was embedded into network weights through optimization. The prior information was then augmented from two motion states to 31 motion states by querying the optimized network at interpolated and extrapolated motion state coordinates. Starting from the prior-augmented NeRP model as an initialization point, we further trained the network to fit sparse samples of two orthogonal MRI slices and the final volumetric reconstruction was obtained by querying the trained network at 3D spatial locations. We evaluated the proposed method using 5-min volumetric MRI time series with 340 ms temporal resolution for seven abdominal patients with hepatocellular carcinoma, acquired using golden-angle radial MRI sequence and reconstructed through retrospective sorting. Two volumetric MRI with inhale and exhale states respectively were selected from the first 30 s of the time series for prior embedding and augmentation. The remaining 4.5-min time series was used for volumetric reconstruction evaluation, where we retrospectively subsampled each MRI to two orthogonal slices and compared model-reconstructed images to ground truth images in terms of image quality and the capability of supporting 3D target motion tracking.

RESULTS

Across the seven patients evaluated, the peak signal-to-noise-ratio between model-reconstructed and ground truth MR images was 38.02 ± 2.60 dB and the structure similarity index measure was 0.98 ± 0.01. Throughout the 4.5-min time period, gross tumor volume (GTV) motion estimated by deforming a reference state MRI to model-reconstructed and ground truth MRI showed good consistency. The 95-percentile Hausdorff distance between GTV contours was 2.41 ± 0.77 mm, which is less than the voxel dimension. The mean GTV centroid position difference between ground truth and model estimation was less than 1 mm in all three orthogonal directions.

CONCLUSION

A prior-augmented NeRP model has been developed to reconstruct volumetric MRI from sparse samples of orthogonal cine slices. Only one exhale and one inhale 3D MRI were needed to train the model to learn prior information of patient breathing motion for sparse image reconstruction. The proposed model has the potential of supporting 3D motion tracking during MR-guided radiotherapy for improved treatment precision and promises a major simplification of the workflow by eliminating the need for large-scale training datasets.

Resolving signal drift in the wall-mounted camera of the RGSC system

PURPOSE

To present a modified calibration method to reduce signal drift due to table sagging in Respiratory Gating for Scanner (RGSC) systems with a wall-mounted camera.

MATERIALS AND METHODS

Approximately 70 kg of solid water phantoms were evenly distributed on the CT couch, mimicking the patient's weight. New calibration measurements were performed at 9 points at the combination of three lateral positions, the CT isocenter and ±10 cm laterally from the isocenter, and three longitudinal locations, the CT isocenter and ±30 cm or ±40 cm from the isocenter. The new calibration was tested in two hospitals.

RESULTS

Implementing the new weighed calibration method at the extended distance yielded improved results during the DIBH scan, reducing the drift to within 1 from 3 mm. The extended calibration positions exhibited similarly reduced drift in both hospitals, reinforcing the method's robustness and its potential applicability across all centers.

CONCLUSION

This proposed solution aims to minimize the systematic error in radiation delivery for patients undergoing motion management with wall-mounted camera RGSC systems, especially in conjunction with a bariatric CT couchtop.

Correcting CT misregistration in data-driven gated (DDG) PET with PET self-gating and deformable image registration

BACKGROUND

Misregistration between CT and PET data can result in mis-localization and inaccurate quantification of functional uptake in whole body PET/CT imaging. This problem is exacerbated when an abnormal inspiration occurs during the free-breathing helical CT (FB CT) used for attenuation correction of PET data. In data-driven gated (DDG) PET, the data selected for reconstruction is typically derived from the end-expiration (EE) phase of the breathing cycle, making this potential issue worse.

PURPOSE

The objective of this study is to develop a deformable image registration (DIR)-based respiratory motion model to improve the registration and quantification between misregistered FB CT and PET.

METHODS

Twenty-two whole-body 18 F-FDG PET/CT scans encompassing 48 lesions in misregistered regions were analyzed in this study. End-inspiration (EI) and EE PET data were derived from -10% to 15% and 30% to 80% of the breathing cycle, respectively. DIR was used to estimate a motion model from the EE to EI phase of the PET data. The model was then used to generate PET images at any phase of up to four times the amplitude of motion between EE and EI for correlation with the misregistered FB CT. Once a matched phase of the FB CT was determined, FB CT was deformed to a pseudo CT at the EE phase (DIR CT). DIR CT was compared with the ground truth DDG CT for AC and localization of the DDG PET.

RESULTS

Between DDG PET/FB CT and DDG PET/DIR CT, a significant increase in ∆%SUV was observed (p < 0.01), with median values elevating from 26.7% to 42.4%. This new method was most effective for lesions ≤3 cm proximal to the diaphragm (p < 0.001) but showed decreasing efficacy as the distance increased. When FB CT was severely misregistered with DDG PET (>3 cm), DDG PET/DIR CT outperformed DDG PET/FB CT alone (p < 0.05). Even when patients showed varied breathing patterns during the PET/CT scan, DDG PET/DIR CT still surpassed the efficiency of DDG PET/FB CT (p < 0.01). Though DDG PET/DIR CT couldn't match the performance of the DDG PET/CT ground truth (42.4% vs. 53.6%, p < 0.01), it reached 84% of its quantification, demonstrating good agreement and a strong overall correlation (regression coefficient of 0.94, p < 0.0001). In some cases, anatomical distortion and blurring, and misregistration error were observed in DIR CT, rendering it still unable to correct inaccurate localization near the boundaries of two organs.

CONCLUSIONS

Based on the motion model derived from gated PET data, DIR CT can significantly improve the quantification and localization of DDG PET. This approach can achieve a performance level of about 84% of the ground truth established by DDG PET/CT. These results show that self-gated PET and DIR CT may offer an alternative clinical solution to DDG PET and FB CT for quantification without the need for additional cine-CT imaging. DIR CT was at times inferior to DDG CT due to some distortion and blurring of anatomy and misregistration error.

Reconstruction of multi-phase parametric maps in 4D-magnetic resonance fingerprinting (4D-MRF) by optimization of local T1 and T2 sensitivities

BACKGROUND

Time-resolved magnetic resonance fingerprinting (MRF), or 4D-MRF, has been demonstrated its feasibility in motion management in radiotherapy (RT). However, the prohibitive long acquisition time is one of challenges of the clinical implementation of 4D-MRF. The shortening of acquisition time causes data insufficiency in each respiratory phase, leading to poor accuracies and consistencies of the predicted tissues' properties of each phase.

PURPOSE

To develop a technique for the reconstruction of multi-phase parametric maps in four-dimensional magnetic resonance fingerprinting (4D-MRF) through the optimization of local T1 and T2 sensitivities.

METHODS

The proposed technique employed an iterative optimization to tailor the data arrangement of each phase by manipulation of inter-phase frames, such that the T1 and T2 sensitivities, which were quantified by the modified Minkowski distance, of the truncated signal evolution curve was maximized. The multi-phase signal evolution curves were modified by sliding window reconstruction and inter-phase frame sharing (SWIFS). Motion correction (MC) and dot product matching were sequentially performed on the modified signal evolution and dictionary to reconstruct the multi-parametric maps. The proposed technique was evaluated by numerical simulations using the extended cardiac-torso (XCAT) phantom with regular and irregular breathing patterns, and by in vivo MRF data of three health volunteers and six liver cancer patients acquired at a 3.0 T scanner.

RESULTS

In simulation study, the proposed SWIFS approach achieved the overall mean absolute percentage error (MAPE) of 8.62% ± 1.59% and 16.2% ± 3.88% for the eight-phases T1 and T2 maps, respectively, in the sagittal view with irregular breathing patterns. In contrast, the overall MAPE of T1 and T2 maps generated by the conventional approach with multiple MRF repetitions were 22.1% ± 11.0% and 30.8% ± 14.9%, respectively. For in-vivo study, the predicted mean T1 and T2 of liver by the proposed SWIFS approach were 795 ms ± 38.9 ms and 58.3 ms ± 11.7 ms, respectively.

CONCLUSIONS

Both simulation and in vivo results showed that the approach empowered by T1 and T2 sensitivities optimization and sliding window under the shortened acquisition of MRF had superior performance in the estimation of multi-phase T1 and T2 maps as compared to the conventional approach with oversampling of MRF data.

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.

And Yet It Moves: Clinical Outcomes and Motion Management in Stereotactic Body Radiation Therapy (SBRT) of Centrally Located Non-Small Cell Lung Cancer (NSCLC): Shedding Light on the Internal Organ at Risk Volume (IRV) Concept

The internal organ at risk volume (IRV) concept might improve toxicity profiles in stereotactic body radiation therapy (SBRT) for non-small cell lung cancer (NSCLC). We studied (1) clinical aspects in central vs. peripheral tumors, (2) the IRV concept in central tumors, (3) organ motion, and (4) associated normal tissue complication probabilities (NTCPs). We analyzed patients who received SBRT for NSCLC (clinical aspects, n = 78; motion management, n = 35). We found lower biologically effective doses, larger planning target volume sizes, higher lung doses, and worse locoregional control for central vs. peripheral tumors. Organ motion was greater in males and tall patients (bronchial tree), whereas volume changes were lower in patients with a high body mass index (BMI) (esophagus). Applying the IRV concept (retrospectively, without new optimization), we found an absolute increase of >10% in NTCPs for the bronchial tree in three patients. This study emphasizes the need to optimize methods to balance dose escalation with toxicities in central tumors. There is evidence that organ motion/volume changes could be more pronounced in males and tall patients, and less pronounced in patients with higher BMI. Since recent studies have made efforts to further subclassify central tumors to refine treatment, the IRV concept should be considered for optimal risk assessment.

Shallow kinetics induced by a metronome (SKIM): A novel contactless respiratory motion management

OBJECTIVES

As an alternative to conventional compression amidst the COVID-19 pandemic, we developed a contactless motion management strategy. By increasing the patient's breathing rate to induce shallow breathing with the aid of a metronome, our hypothesis is that the motion magnitude of the target may be minimized without physical contact or compression.

METHODS

Fourteen lung stereotactic body radiation therapy (SBRT) patients treated under fast shallow-breathing (FSB) were selected for inclusion in this retrospective study. Our proposed method is called shallow kinetics induced by a metronome (SKIM). We induce FSB by setting the beats-per-minute (BPM) high (typically in the range of 50-60). This corresponded to a patient breathing rate of 25-30 (breathing) cycles per minute. The magnitude of target motion in 3D under SKIM was evaluated using 4DCT datasets. Comparison with free breathing (FB) 4DCT was also made for a subset for which FB data available.

RESULTS

The overall effectiveness of SKIM was evaluated with 18 targets (14 patients). Direct comparison with FB was performed with 12 targets (10 patients). The vector norm mean ± SD value of motion magnitude under SKIM for 18 targets was 8.2 ± 4.1 mm. The mean ± SD metronome BPM was 54.9 ± 4.0 in this group. The vector norm means ± SD values of target motion for FB and SKIM in the 12 target sub-group were 14.6 ± 8.5 mm and 9.3 ± 3.7 mm, respectively. The mean ± SD metronome BPM for this sub-group was 56.3 ± 2.5.

CONCLUSION

Compared with FB, SKIM can significantly reduce respiratory motion magnitude of thoracic targets. The difference in maximum motion reduction in the overall vector norm, S-I, and A-P directions was significant (p = 0.033, 0.042, 0.011). Our proposed method can be an excellent practical alternative to conventional compression due to its flexibility and ease of implementation.

Technical note: Feasibility of gating for dynamic trajectory radiotherapy - Mechanical accuracy and dosimetric performance

BACKGROUND

Dynamic trajectory radiotherapy (DTRT) extends state-of-the-art volumetric modulated arc therapy (VMAT) by dynamic table and collimator rotations during beam-on. The effects of intrafraction motion during DTRT delivery are unknown, especially regarding the possible interplay between patient and machine motion with additional dynamic axes.

PURPOSE

To experimentally assess the technical feasibility and quantify the mechanical and dosimetric accuracy of respiratory gating during DTRT delivery.

METHODS

A DTRT and VMAT plan are created for a clinically motivated lung cancer case and delivered to a dosimetric motion phantom (MP) placed on the table of a TrueBeam system using Developer Mode. The MP reproduces four different 3D motion traces. Gating is triggered using an external marker block, placed on the MP. Mechanical accuracy and delivery time of the VMAT and DTRT deliveries with and without gating are extracted from the logfiles. Dosimetric performance is assessed by means of gamma evaluation (3% global/2 mm, 10% threshold).

RESULTS

The DTRT and VMAT plans are successfully delivered with and without gating for all motion traces. Mechanical accuracy is similar for all experiments with deviations <0.14° (gantry angle), <0.15° (table angle), <0.09° (collimator angle) and <0.08 mm (MLC leaf positions). For DTRT (VMAT), delivery times are 1.6-2.3 (1.6- 2.5) times longer with than without gating for all motion traces except one, where DTRT (VMAT) delivery is 5.0 (3.6) times longer due to a substantial uncorrected baseline drift affecting only DTRT delivery. Gamma passing rates with (without) gating for DTRT/VMAT were ≥96.7%/98.5% (≤88.3%/84.8%). For one VMAT arc without gating it was 99.6%.

CONCLUSION

Gating is successfully applied during DTRT delivery on a TrueBeam system for the first time. Mechanical accuracy is similar for VMAT and DTRT deliveries with and without gating. Gating substantially improved dosimetric performance for DTRT and VMAT.

Carbon-fiber alternative to the commercial gating surrogate for the Varian Truebeam™

BACKGROUND

In this study we present the Tracking Accessory 3 (TA3) as an alternative to the commercial gating block (GB) surrogate for the Varian Truebeam™ gating system (TGS). The TGS requires three visible reflectors to track the surrogate, presenting an opportunity for a surrogate to be made with less material and thus smaller dosimetric footprint than the commercial four reflector model.

MATERIALS AND METHODS

Relative dose and depth dose profiles below the TA3 and the GB were measured with radiosensitive film. Accuracy and reproducibility of the detected motion amplitude for three TA3s and one GB were determined using a respiratory phantom with surrogate to determine the camera's tracking volume. Clinical performance was evaluated prospectively in 10 breast cancer patients treated with deep inspiration breath hold monitored with TA3 and compared to previously published results. Non-parametric statistics were applied to test for significance.

RESULTS AND CONCLUSIONS

Surface doses were increased up to 94% and 187% for the TA3 and GB, respectively, compared to no surrogate. The surface area influenced by at least 25% increase in dose was 12 cm2 and 105 cm2 for the TA3 and GB, respectively. The water equivalent thickness of the surrogates was found to be 1 mm for the TA3 and 3 mm for GB. The difference in measured amplitude were <0.2 mm for TA3 compared to the GB. The TA3s and GB were detected at all extremes of the clinically relevant tracking volume of the TGS. Clinical performance showed no significant differences. The TA3 caused less surface dose increase compared to the commercial GB. In the tested range all surrogates measured motion amplitude within 0.2 mm of reference value, which is not a clinically relevant difference. The TA3 showed no significant differences in clinical performance to similarly positioned surrogates.

Planning CT Identifies Patients at Risk of High Prostate Intrafraction Motion

Prostate motion (standard deviation, range of motion, and diffusion coefficient) was calculated from 4D ultrasound data of 1791 fractions of radiation therapy in N = 100 patients. The inner diameter of the lesser pelvis was obtained from transversal slices through the pubic symphysis in planning CTs. On the lateral and craniocaudal axes, motility increases significantly (t-test, p < 0.005) with the inner diameter of the lesser pelvis. A diameter of >106 mm (ca. 6th decile) is a good predictor for high prostate intrafraction motion (ca. 9th decile). The corresponding area under the receiver operator curve (AUROC) is 80% in the lateral direction, 68% to 80% in the craniocaudal direction, and 62% to 70% in the vertical direction. On the lateral x-axis, the proposed test is 100% sensitive and has a 100% negative predictive value for all three characteristics (standard deviation, range of motion, and diffusion coefficient). On the craniocaudal z-axis, the proposed test is 79% to 100% sensitive and reaches 95% to 100% negative predictive value. On the vertical axis, the proposed test still delivers 98% negative predictive value but is not particularly sensitive. Overall, the proposed predictor is able to help identify patients at risk of high prostate motion based on a single planning CT.

A patient-specific deep learning framework for 3D motion estimation and volumetric imaging during lung cancer radiotherapy

Objective. Respiration introduces a constant source of irregular motion that poses a significant challenge for the precise irradiation of thoracic and abdominal cancers. Current real-time motion management strategies require dedicated systems that are not available in most radiotherapy centers. We sought to develop a system that estimates and visualises the impact of respiratory motion in 3D given the 2D images acquired on a standard linear accelerator.Approach. In this paper we introduceVoxelmap, a patient-specific deep learning framework that achieves 3D motion estimation and volumetric imaging using the data and resources available in standard clinical settings. Here we perform a simulation study of this framework using imaging data from two lung cancer patients.Main results. Using 2D images as input and 3D-3DElastixregistrations as ground-truth,Voxelmapwas able to continuously predict 3D tumor motion with mean errors of 0.1 ± 0.5, -0.6 ± 0.8, and 0.0 ± 0.2 mm along the left-right, superior-inferior, and anterior-posterior axes respectively.Voxelmapalso predicted 3D thoracoabdominal motion with mean errors of -0.1 ± 0.3, -0.1 ± 0.6, and -0.2 ± 0.2 mm respectively. Moreover, volumetric imaging was achieved with mean average error 0.0003, root-mean-squared error 0.0007, structural similarity 1.0 and peak-signal-to-noise ratio 65.8.Significance. The results of this study demonstrate the possibility of achieving 3D motion estimation and volumetric imaging during lung cancer treatments on a standard linear accelerator.

Real-time myocardial landmark tracking for MRI-guided cardiac radio-ablation using Gaussian Processes

The high speed of cardiorespiratory motion introduces a unique challenge for cardiac stereotactic radio-ablation (STAR) treatments with the MR-linac. Such treatments require tracking myocardial landmarks with a maximum latency of 100 ms, which includes the acquisition of the required data. The aim of this study is to present a new method that enables tracking myocardial landmarks from few readouts of MRI data, thereby achieving a latency sufficient for STAR treatments. We present a tracking framework that requires few readouts of k-space data as input, which can be acquired at least an order of magnitude faster than MR-images. Combined with the real-time tracking speed of a probabilistic machine learning framework called Gaussian Processes, this allows to track myocardial landmarks with a sufficiently low latency for cardiac STAR guidance. This includes both the acquisition of required data, and the tracking inference. The framework is demonstrated in 2D on a motion phantom, and in vivo on volunteers and a ventricular tachycardia (arrhythmia) patient. Moreover, the feasibility of an extension to 3D was demonstrated by in silico 3D experiments with a digital motion phantom. The framework was compared with template matching - a reference, image-based, method - and linear regression methods. Results indicate an order of magnitude lower total latency (<10 ms) for the proposed framework in comparison with alternative methods. The root-mean-square-distances and mean end-point-distance with the reference tracking method was less than 0.8 mm for all experiments, showing excellent (sub-voxel) agreement. The high accuracy in combination with a total latency of less than 10 ms - including data acquisition and processing - make the proposed method a suitable candidate for tracking during STAR treatments. Additionally, the probabilistic nature of the Gaussian Processes also gives access to real-time prediction uncertainties, which could prove useful for real-time quality assurance during treatments.

Surface guided electron FLASH radiotherapy for canine cancer patients

BACKGROUND

During recent years FLASH radiotherapy (FLASH-RT) has shown promising results in radiation oncology, with the potential to spare normal tissue while maintaining the antitumor effects. The high speed of the FLASH-RT delivery increases the need for fast and precise motion monitoring to avoid underdosing the target. Surface guided radiotherapy (SGRT) uses surface imaging (SI) to render a 3D surface of the patient. SI provides real-time motion monitoring and has a large scanning field of view, covering off-isocentric positions. However, SI has so far only been used for human patients with conventional setup and treatment.

PURPOSE

The aim of this study was to investigate the performance of SI as a motion management tool during electron FLASH-RT of canine cancer patients.

METHODS

To evaluate the SI system's ability to render surfaces of fur, three fur-like blankets in white, grey, and black were used to imitate the surface of canine patients and the camera settings were optimized for each blanket. Phantom measurements using the fur blankets were carried out, simulating respiratory motion and sudden shift. Respiratory motion was simulated using the QUASAR Respiratory Motion Phantom with the fur blankets placed on the phantom platform, which moved 10 mm vertically with a simulated respiratory period of 4 s. Sudden motion was simulated with an in-house developed phantom, consisting of a platform which was moved vertically in a stepwise motion at a chosen frequency. For sudden measurements, 1, 2, 3, 4, 5, 6, 7, and 10 Hz were measured. All measurements were both carried out at the conventional source-to-surface distance (SSD) of 100 cm, and in the locally used FLASH-RT setup at SSD = 70 cm. The capability of the SI system to reproduce the simulated motion and the sampling time were evaluated. As an initial step towards clinical implementation, the feasibility of SI for surface guided FLASH-RT was evaluated for 11 canine cancer patients.

RESULTS

The SI camera was capable of rendering surfaces for all blankets. The deviation between simulated and measured mean peak-to-peak breathing amplitude was within 0.6 mm for all blankets. The sampling time was generally higher for the black fur than for the white and grey fur, for the measurement of both respiratory and sudden motion. The SI system could measure sudden motion within 62.5 ms and detect motion with a frequency of 10 Hz. The feasibility study of the canine patients showed that the SI system could be an important tool to ensure patient safety. By using this system we could ensure and document that 10 out of 11 canine patients had a total vector offset from the reference setup position <2 mm immediately before and after irradiation.

CONCLUSIONS

We have shown that SI can be used for surface guided FLASH-RT of canine patients. The SI system is currently not fast enough to interrupt a FLASH-RT beam while irradiating but with the short sampling time sudden motion can be detected. The beam can therefore be held just prior to irradiation, preventing treatment errors such as underdosing the target.

Deep learning enables MV-based real-time image guided radiation therapy for prostate cancer patients

Using MV images for real-time image guided radiation therapy (IGRT) is ideal as it does not require additional imaging equipment, adds no additional imaging dose and provides motion data in the treatment beam frame of reference. However, accurate tracking using MV images is challenging due to low contrast and modulated fields. Here, a novel real-time marker tracking system based on a convolutional neural network (CNN) classifier was developed and evaluated on retrospectively acquired patient data for MV-based IGRT for prostate cancer patients. &#xD;&#xD;MV images, acquired from 29 VMAT prostate cancer patients treated in a multi-institutional clinical trial, were used to train and evaluate a CNN-based marker tracking system. The CNN was trained using labelled MV images from 9 prostate cancer patients (35 fractions) with implanted markers. CNN performance was evaluated on an independent cohort of unseen MV images from 20 patients (78 fractions), using a Precision-Recall curve (PRC), area under the PRC plot (AUC) and sensitivity and specificity. The accuracy of the tracking system was evaluated on the same unseen dataset and quantified by calculating mean absolute (± 1 SD) and [1st, 99th] percentiles of the geometric tracking error in treatment beam co-ordinates using manual identification as the ground truth. &#xD;&#xD;The CNN had an AUC of 0.99, sensitivity of 98.31% and specificity of 99.87%. The mean absolute geometric tracking error was 0.30 ± 0.27 and 0.35 ± 0.31 mm in the lateral and superior-inferior directions of the MV images, respectively. The [1st, 99th] percentiles of the error were [-1.03, 0.90] and [-1.12, 1.12] mm in the lateral and SI directions, respectively.&#xD;&#xD;The high classification performance on unseen MV images demonstrates the CNN can successfully identify implanted prostate markers. Furthermore, the sub-millimetre accuracy and precision of the marker tracking system demonstrates potential for adaptation to real-time applications.

Cardiac motion and its dosimetric impact during radioablation for refractory ventricular tachycardia

INTRODUCTION

Cardiac radioablation (CR) is a noninvasive treatment option for patients with refractory ventricular tachycardia (VT) during which high doses of radiation, typically 25 Gy, are delivered to myocardial scar. In this study, we investigate motion from cardiac cycle and evaluate the dosimetric impact in a cohort of patients treated with CR.

METHODS

This retrospective study included eight patients treated at our institution who had respiratory-correlated and ECG-gated 4DCT scans acquired within 2 weeks of CR. Deformable image registration was applied between maximum systole (SYS) and diastole (DIAS) CTs to assess cardiac motion. The average respiratory-correlated CT (AVG_resp ) was deformably registered to the average cardiac (AVG_cardiac ), SYS, and DIAS CTs, and contours were propagated using the deformation vector fields (DVFs). Finally, the original treatment plan was recalculated on the deformed AVG_resp CT for dosimetric assessment.

RESULTS

Motion magnitudes were measured as the mean (SD) value over the DVFs within each structure. Displacement during the cardiac cycle for all chambers was 1.4 (0.9) mm medially/laterally (ML), 1.6 (1.0) mm anteriorly/posteriorly (AP), and 3.0 (2.8) mm superiorly/inferiorly (SI). Displacement for the 12 distinct clinical target volumes (CTVs) was 1.7 (1.5) mm ML, 2.4 (1.1) mm AP, and 2.1 (1.5) SI. Displacements between the AVG_resp and AVG_cardiac scans were 4.2 (2.0) mm SI and 5.8 (1.4) mm total. Dose recalculations showed that cardiac motion may impact dosimetry, with dose to 95% of the CTV dropping from 27.0 (1.3) Gy on the AVG_resp to 20.5 (7.1) Gy as estimated on the AVG_cardiac .

CONCLUSIONS

Cardiac CTV motion in this patient cohort is on average below 3 mm, location-dependent, and when not accounted for in treatment planning may impact target coverage. Further study is needed to assess the impact of cardiac motion on clinical outcomes.

Initial clinical experience with magnetic resonance-guided radiotherapy in pediatric patients: Lessons learned from a single institution with proton therapy

Purpose/Objectives

Magnetic resonance-guided radiotherapy (MRgRT) is increasingly used in a variety of adult cancers. To date, published experience regarding the use of MRgRT in pediatric patients is limited to two case reports. We report on the use of MRgRT for pediatric patients at our institution during a four-year period and describe important considerations in the selection and application of this technology in children.

Materials/Methods

All patients treated with MRgRT since inception at our institution between 4/2018 and 4/2022 were retrospectively reviewed. We also evaluated all pediatric patients treated at our institution during the same period who received either imaging or treatment using our magnetic resonance-guided linear accelerator (MR Linac). We summarize four clinical cases where MRgRT was selected for treatment in our clinic, including disease outcomes and toxicities and describe our experience using the MR Linac for imaging before and during treatment for image fusion and tumor assessments.

Results

Between 4/2018 and 4/2022, 535 patients received MRgRT at our center, including 405 (75.7%) with stereotactic ablative radiotherapy (SABR). During this period, 347 distinct radiotherapy courses were delivered to pediatric patients, including 217 (62.5%) with proton therapy. Four pediatric patients received MRgRT. One received SABR for lung metastasis with daily adaptive replanning and a second was treated for liver metastasis using a non-adaptive workflow. Two patients received fractionated MRgRT for an ALK-rearranged non-small cell lung cancer and neuroblastoma. No Grade 2 or higher toxicities were observed or reported during MRgRT or subsequent follow-up. Twelve patients underwent MR imaging without contrast during treatment for brain tumors to assess for tumor/cystic changes. Two patients treated with other modalities underwent MR simulation for target volume delineation and organ at risk sparing due to anatomic changes during treatment or unexpected delays in obtaining diagnostic MR appointments.

Conclusions

In four pediatric patients treated with MRgRT, treatment was well tolerated with no severe acute effects. At our center, most pediatric patients are treated with proton therapy, but the cases selected for MRgRT demonstrated significant organ at risk sparing compared to alternative modalities. In particular, MRgRT may provide advantages for thoracic/abdominal/pelvic targets using gated delivery and adaptive replanning, but selected patients treated with fractionated radiotherapy may also benefit MRgRT through superior organ at risk sparing.

Institutional experience of using active breathing control for paediatric and teenage patients receiving thoraco-abdominal radiotherapy

Introduction

Active Breathing Control (ABC) is a motion management strategy that facilitates reproducible breath-hold for thoracic radiotherapy (RT), which may reduce radiation dose to organs at risk (OARs). Reduction of radiation-induced toxicity is of high importance in younger patients. However, there is little published literature on the feasibility of ABC in this group. The purpose of this study was to report our experience of using ABC for paediatric and teenage patients.

Methods

Patients ≤18 years referred for thoracic RT using ABC at our centre from 2013-2021 were identified. Electronic records were retrospectively reviewed to obtain information on diagnosis, RT dose and technique, OAR dosimetry, tolerability of ABC, post-treatment imaging and early toxicity rates.

Results

12 patients completed RT and were able to comply with ABC during planning and for the duration of RT. Median age was 15.5 years (10-18 years). Diagnoses were: Hodgkin lymphoma (n = 5), mediastinal B-cell lymphoma (n = 1), Ewing sarcoma (n = 5) and rhabdomyosarcoma (n = 1). For mediastinal RT cases (n = 6), median dose delivered was 30.6Gy(19.8-40Gy), median mean heart dose was 11.4Gy(4.8-19.4Gy), median mean lung dose was 9.9Gy(5.7-14.5Gy) and mean lung V20 was 10.9%. For ipsilateral RT cases, (n = 6), median hemithorax and total doses to primary tumour were 18Gy(15-20Gy) and 52.2Gy(36-60Gy) respectively. Median mean heart dose was 19.5Gy(10.6-33.2Gy) and median mean lung dose was 17.7Gy(16.3-30.5Gy). Mean bilateral lung V20 was 39.6%. Median mean contralateral lung dose was 5.2Gy(3.5-11.6Gy) and mean contralateral lung V20 was 1.5%. At a median follow-up of 36 months, only 1 patient had symptomatic radiation pneumonitis having received further thoracic RT following relapse.

Conclusions

ABC is feasible and well tolerated in younger patients receiving RT. Children as young as 10 years are able to comply. Use of ABC results in OAR dosimetry which is comparable to similar data in adults and can facilitate RT for extensive thoracic sarcoma.

Intrafraction tumor motion monitoring and dose reconstruction for liver pencil beam scanning proton therapy

Background

Pencil beam scanning (PBS) proton therapy can provide highly conformal target dose distributions and healthy tissue sparing. However, proton therapy of hepatocellular carcinoma (HCC) is prone to dosimetrical uncertainties induced by respiratory motion. This study aims to develop intra-treatment tumor motion monitoring during respiratory gated proton therapy and combine it with motion-including dose reconstruction to estimate the delivered tumor doses for individual HCC treatment fractions.

Methods

Three HCC-patients were planned to receive 58 GyRBE (n=2) or 67.5 GyRBE (n=1) of exhale respiratory gated PBS proton therapy in 15 fractions. The treatment planning was based on the exhale phase of a 4-dimensional CT scan. Daily setup was based on cone-beam CT (CBCT) imaging of three implanted fiducial markers. An external marker block (RPM) on the patient's abdomen was used for exhale gating in free breathing. This study was based on 5 fractions (patient 1), 1 fraction (patient 2) and 6 fractions (patient 3) where a post-treatment control CBCT was available. After treatment, segmented 2D marker positions in the post-treatment CBCT projections provided the estimated 3D motion trajectory during the CBCT by a probability-based method. An external-internal correlation model (ECM) that estimated the tumor motion from the RPM motion was built from the synchronized RPM signal and marker motion in the CBCT. The ECM was then used to estimate intra-treatment tumor motion. Finally, the motion-including CTV dose was estimated using a dose reconstruction method that emulates tumor motion in beam's eye view as lateral spot shifts and in-depth motion as changes in the proton beam energy. The CTV homogeneity index (HI) The CTV homogeneity index (HI) was calculated as D 2 %  -  D 98 % D 50 %   × 100 % .

Results

The tumor position during spot delivery had a root-mean-square error of 1.3 mm in left-right, 2.8 mm in cranio-caudal and 1.7 mm in anterior-posterior directions compared to the planned position. On average, the CTV HI was larger than planned by 3.7%-points (range: 1.0-6.6%-points) for individual fractions and by 0.7%-points (range: 0.3-1.1%-points) for the average dose of 5 or 6 fractions.

Conclusions

A method to estimate internal tumor motion and reconstruct the motion-including fraction dose for PBS proton therapy of HCC was developed and demonstrated successfully clinically.

Proton liver stereotactic body radiation therapy: Treatment techniques and dosimetry feasibility from a single institution

Purpose

To assess the resulting dosimetry characteristics of simulation and planning techniques for proton stereotactic body radiation therapy (SBRT) of primary and secondary liver tumors.

Methods

Consecutive patients treated under volumetric daily image guidance with liver proton SBRT between September 2019 and March 2022 at Emory Proton Therapy Center were included in this study. Prescriptions ranged from 40 Gy to 60 Gy in 3- or 5-fraction regimens, and motion management techniques were used when target motion exceeded 5 mm. 4D robust optimization was used when necessary. Dosimetry evaluation was conducted for ITV V100, D99, Dmax, and liver-ITV mean dose and D700cc. Statistical analysis was performed using independent-samples Mann-Whitney U tests.

Results

Thirty-six tumors from 29 patients were treated. Proton therapy for primary and secondary liver tumors using motion management techniques and robust optimization resulted in high target coverage and low doses to critical organs. The median ITV V100% was 100.0%, and the median ITV D99% was 111.3%. The median liver-ITV mean dose and D700cc were 499 cGy and 5.7 cGy, respectively. The median conformity index (CI) was 1.03, and the median R50 was 2.56. Except for ITV D99% (primary 118.1% vs. secondary 107.2%, p = 0.005), there were no significant differences in age, ITV volume, ITV V100%, ITV maximum dose, liver-ITV mean dose, or D700cc between primary and secondary tumor groups.

Conclusion

The study demonstrated that proton therapy with motion management techniques and robust optimization achieves excellent target coverage with low normal liver doses for primary and secondary liver tumors. The results showed high target coverage, high conformality, and low doses to the liver.

Advances in treatment planning and management for the safety and accuracy of lung stereotactic body radiation therapy using proton pencil beam scanning: Simulation, planning, quality assurance, and delivery recommendations

This study presents the clinical experiences of the New York Proton Center in employing proton pencil beam scanning (PBS) for the treatment of lung stereotactic body radiation therapy. It encompasses a comprehensive examination of multiple facets, including patient simulation, delineation of target volumes and organs at risk, treatment planning, plan evaluation, quality assurance, and motion management strategies. By sharing the approaches of the New York Proton Center and providing recommendations across simulation, treatment planning, and treatment delivery, it is anticipated that the valuable experience will be provided to a broader proton therapy community, serving as a useful reference for future clinical practice and research endeavors in the field of stereotactic body proton therapy for lung tumors.

Monte Carlo Modeling of Dynamic Tumor Tracking on a Gimbaled Linear Accelerator

Purpose and Aim

The Vero4DRT (Brainlab AG) linear accelerator is capable of dynamic tumor tracking (DTT) by panning/tilting the radiation beam to follow respiratory-induced tumor motion in real time. In this study, the panning/tilting motion is modeled in Monte Carlo (MC) for quality assurance (QA) of four-dimensional (4D) dose distributions created within the treatment planning system (TPS).

Materials and Methods

Step-and-shoot intensity-modulated radiation therapy plans were optimized for 10 previously treated liver patients. These plans were recalculated on multiple phases of a 4D computed tomography (4DCT) scan using MC while modeling panning/tilting. The dose distributions on each phase were accumulated to create a respiratory-weighted 4D dose distribution. Differences between the TPS and MC modeled doses were examined.

Results

On average, 4D dose calculations in MC showed the maximum dose of an organ at risk (OAR) to be 10% greater than the TPS' three-dimensional dose calculation (collapsed cone [CC] convolution algorithm) predicted. MC's 4D dose calculations showed that 6 out of 24 OARs could exceed their specified dose limits, and calculated their maximum dose to be 4% higher on average (up to 13%) than the TPS' 4D dose calculations. Dose differences between MC and the TPS were greatest in the beam penumbra region.

Conclusion

Modeling panning/tilting for DTT has been successfully modeled with MC and is a useful tool to QA respiratory-correlated 4D dose distributions. The dose differences between the TPS and MC calculations highlight the importance of using 4D MC to confirm the safety of OAR doses before DTT treatments.

The dosimetric error due to uncorrected tumor rotation during real-time adaptive prostate stereotactic body radiation therapy

BACKGROUND

During prostate stereotactic body radiation therapy (SBRT), prostate tumor translational motion may deteriorate the planned dose distribution. Most of the major advances in motion management to date have focused on correcting this one aspect of the tumor motion, translation. However, large prostate rotation up to 30° has been measured. As the technological innovation evolves toward delivering increasingly precise radiotherapy, it is important to quantify the clinical benefit of translational and rotational motion correction over translational motion correction alone.

PURPOSE

The purpose of this work was to quantify the dosimetric impact of intrafractional dynamic rotation of the prostate measured with a six degrees-of-freedom tumor motion monitoring technology.

METHODS

The delivered dose was reconstructed including (a) translational and rotational motion and (b) only translational motion of the tumor for 32 prostate cancer patients recruited on a 5-fraction prostate SBRT clinical trial. Patients on the trial received 7.25 Gy in a treatment fraction. A 5 mm clinical target volume (CTV) to planning target volume (PTV) margin was applied in all directions except the posterior direction where a 3 mm expansion was used. Prostate intrafractional translational motion was managed using a gating strategy, and any translation above the gating threshold was corrected by applying an equivalent couch shift. The residual translational motion is denoted as T r e s $T_{res}$ . Prostate intrafractional rotational motion R u n c o r r $R_{uncorr}$ was recorded but not corrected. The dose differences from the planned dose due to T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ , ΔD( T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ ) and due to T r e s $T_{res}$ alone, ΔD( T r e s $T_{res}$ ), were then determined for CTV D98, PTV D95, bladder V6Gy, and rectum V6Gy. The residual dose error due to uncorrected rotation, R u n c o r r $R_{uncorr}$ was then quantified: Δ D R e s i d u a l $\Delta D_{Residual}$ = ΔD( T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ ) - ΔD( T res ${T}_{\textit{res}}$ ).

RESULTS

Fractional data analysis shows that the dose differences from the plan (both ΔD( T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ ) and ΔD( T r e s $T_{res}$ )) for CTV D98 was less than 5% in all treatment fractions. ΔD( T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ ) was larger than 5% in one fraction for PTV D95, in one fraction for bladder V6Gy, and in five fractions for rectum V6Gy. Uncorrected rotation, R u n c o r r $R_{uncorr}$ induced residual dose error, Δ D R e s i d u a l $\Delta D_{Residual}$ , resulted in less dose to CTV and PTV in 43% and 59% treatment fractions, respectively, and more dose to bladder and rectum in 51% and 53% treatment fractions, respectively. The cumulative dose over five fractions, ∑D( T r e s $T_{res}$ + R u n c o r r $R_{uncorr}$ ) and ∑D( T r e s $T_{res}$ ), was always within 5% of the planned dose for all four structures for every patient.

CONCLUSIONS

The dosimetric impact of tumor rotation on a large prostate cancer patient cohort was quantified in this study. These results suggest that the standard 3-5 mm CTV-PTV margin was sufficient to account for the intrafraction prostate rotation observed for this cohort of patients, provided an appropriate gating threshold was applied to correct for translational motion. Residual dose errors due to uncorrected prostate rotation were small in magnitude, which may be corrected using different treatment adaptation strategies to further improve the dosimetric accuracy.

Surface guided 3DCRT in deep-inspiration breath-hold for left sided breast cancer radiotherapy: implementation and first clinical experience in Iran

Background

The aim of the study is to evaluate the overall accuracy of the surface-guided radiotherapy (SGRT) workflow through a comprehensive commissioning and quality assurance procedures and assess the potential benefits of deep-inspiration breath-hold (DIBH) radiotherapy as a cardiac and lung dose reduction approach for left-sided breast cancer irradiation.

Materials and methods

Accuracy and reproducibility of the optical surface scanner used for DIBH treatment were evaluated using different phantoms. Patient positioning accuracy and reproducibility of DIBH treatment were evaluated. Twenty patients were studied for treatment plan quality in target dose coverage and healthy organ sparing for the two different treatment techniques.

Results

Reproducibility tests for the surface scanner showed good stability within 1 mm in all directions. The maximum position variation between applied shifts on the couch and the scanner measured offsets is 1 mm in all directions. The clinical study of 200 fractions showed good agreement between the surface scanner and portal imaging with the isocenter position deviation of less than 3 mm in each lateral, longitudinal, and vertical direction. The standard deviation of the DIBH level showed a value of < 2 mm during all evaluated DIBHs. Compared to the free breathing (FB) technique, DIBH showed significant reduction of 48% for heart mean dose, 43% for heart V25, and 20% for ipsilateral lung V20.

Conclusion

Surface-guided radiotherapy can be regarded as an accurate tool for patient positioning and monitoring in breast radiotherapy. DIBH treatment are considered to be effective techniques in heart and ipsilateral lung dose reductions for left breast radiotherapy.

Experience With Normal Breathhold Planning Scans for Radiosurgery of Moving Targets With Live Tracking

PURPOSE

Utilization of breathhold scans with live tracking has a long track record of good published outcomes for stereotactic body radiation therapy (SBRT) and is recommended by the manufacturer of the Synchrony tracking system. However, the popularity of four-dimensional computed tomography (4DCT) scans challenges the validity of the breathhold scan with live tracking technique. Although this study is not intended to prove the superiority of either method, we demonstrate the feasibility of using the breathhold scans with a phantom test and clinical examples.

METHODS

A 4DCT of a perfect sphere was scanned at 20 breaths per minute and compared to a 4DCT of a small lung tumor in one patient and a 4DCT of a larger renal tumor in another patient, as well as to fiducial matching in a patient with pancreatic cancer. Normal exhale and normal inhale breathhold CT scans were performed for the pancreatic cancer patient, combined with Synchrony tracking on CyberKnife (Sunnyvale, CA: Accuray) for treatment.

RESULTS

The 4DCT scan of the phantom exhibited considerable apparent deformation, which must be entirely due to imaging artifact since the perfect sphere in the phantom is known to be completely rigid. The 4DCT of the lung and renal tumors in patients had similar apparent deformation. Usually in patients, from 4DCT alone, it is difficult to determine how much was due to deformation and how much was due to artifact. Fiducial positions in the final normal exhale and normal inhale breathhold scans for Synchrony matched each other within 1mm for the pancreatic cancer patient.

CONCLUSION

We demonstrated the feasibility of breathhold scans with Synchrony live tracking, as recommended by the manufacturer. More studies will be needed to determine whether this method is better than using a 4DCT.

Measurement of the temporal latency of a respiratory gating system using two distinct approaches

Purpose

To develop a methodology that can be used to measure the temporal latency of a respiratory gating system.

Methods

The gating system was composed of an automatic gating interface (Response) and an in‐house respiratory motion monitoring system featuring an optically tracked surface marker. Two approaches were used to measure gating latencies. A modular approach involved measuring separately the latency of the gating system's complementary metal–oxide–semiconductor tracking camera, tracking software, and a gating latency of the LINAC. Additionally, an end‐to‐end approach was used to measure the total latency of the gating system. End‐to‐end latencies were measured using the displacement of a radiographic target moving at known velocities during the gating process.

Results

Summing together the latencies of each of the modular components investigated yielded a total beam‐on latency of 1.55 s and a total beam‐off latency of 0.49 s. End‐to‐end beam‐on and beam‐off latency was found to be 1.49 and 0.34 s, respectively. In each case, no statistically significant differences were found between the end‐to‐end latency of the gating system and the summation of the individually measured components.

Conclusion

Two distinct approaches to quantify gating latencies were presented. Measuring the end‐to‐end latency of the gating system provided an independent means of validating the modular approach. It is expected that the beam‐on latencies reported in this work could be reduced by altering the control system configuration of the LINAC. The modular approach can be used to decouple the individual latencies of the gating system, but future improvements in the temporal resolution of the service graphing feature are needed to reduce the uncertainty of LINAC‐related gating latencies measured using this approach. Both approaches are generalizable and can be used together when designing a quality assurance program for a respiratory gating system.

Efficacy and Optimal Pressure of Continuous Positive Airway Pressure in Intensity-Modulated Radiotherapy for Locally Advanced Lung Cancer

We aimed to determine the optimal pressure of continuous positive airway pressure (CPAP) for radiotherapy (RT) through changes in the dosimetric parameters and lung volume according to pressure. Patients with locally advanced lung cancer, who underwent CPAP during computed tomography (CT) simulation, were included. The air pressure was raised in five steps of 4, 7, 10, 14, and 17 cmH2O and a CT scan was performed at the baseline and at each pressure step, accompanied by contouring and RT planning. Paired t- and Wilcoxon signed rank tests were used to compare the volumetric and dosimetric parameters according to pressure and interpressure. A total of 29 patients were selected, and 158 CT datasets were obtained. The lung volume increased significantly at all pressures (p < 0.01). The Dmean of the lung decreased significantly from 7 cmH2O (p < 0.01), the V5, V10, V15, and V20 of the lung decreased significantly from 7 cmH2O with increasing pressure, and the Dmean and V5 of the heart decreased significantly from 14 cmH2O with increasing pressure. The V50 showed no significant differences at any pressure. We recommend the use of at least 7 cmH2O with 14 cmH2O as the optimal pressure to achieve the effect of heart preservation.

An integrated ultrasound imaging and abdominal compression device for respiratory motion management in radiation therapy

BACKGROUND

Radiotherapy to tumors in the abdomen is challenging because of the significant organ movement and tissue deformation caused by respiration.

PURPOSE

A motion management strategy that integrated ultrasound (US) imaging with abdominal compression was developed and evaluated, where US was used to real-time monitor organ motion after abdominal compression.

METHODS

A device that combined a US imaging system and an abdominal compression plate (ACP) was developed. Twenty-one healthy volunteers were involved to evaluate the motion management efficacy. Each volunteer was immobilized on a flat bench by the device. Abdominal US data were successively collected with and without ACP compression and experiments were repeated three times to verify the imaging reproducibility. A template matching algorithm based on normalized cross correlation (NCC) was implemented to track the targets (vessels in the liver, pancreas and stomach) automatically. The matching algorithm was validated by comparing with the manual references. Automatic tracking was judged as failed if the center of mass difference from manual tracking was beyond a failure threshold. Based on the locations obtained through the template matching algorithm, the motion correlation between liver and pancreas/stomach was investigated using Pearson correlation test. Paired Student's t-test was used to analyze the difference between the results without and with ACP compression.

RESULTS

The liver motion amplitude over all 21 volunteers was significantly (p<0.001) reduced from 14.9 ± 5.5/3.4 ± 1.8 mm in superior-inferior (SI)/anterior-posterior (AP) direction before ACP compression to 7.3 ± 1.5/1.6 ± 0.7 mm after ACP compression. The mean liver motion standard deviation before compression was on average 2.8/1.4 mm in SI/AP direction, and was significantly (p<0.001) reduced to 0.9/0.4 mm after compression. The failure rates of automatic tracking for liver, pancreas and stomach were reduced for failure thresholds of 1-5 mm after applying ACP. The Pearson correlation coefficients between liver and pancreas/stomach were 0.98/0.97 without ACP and 0.96/0.94 with ACP in SI direction, and were 0.68/0.68 and 0.43/0.42 in AP direction. The motion prediction errors for pancreas/stomach with ACP have significantly (p<0.001) reduced to 0.45 ± 0.36/0.52 ± 0.43 mm from 0.69 ± 0.56/0.71 ± 0.66 mm without ACP in SI direction, and to 0.38 ± 0.33/0.39 ± 0.27 mm from 0.44 ± 0.35/0.61 ± 0.59 mm in AP direction.

CONCLUSIONS

The proposed strategy that combines real-time US imaging and abdominal compression has the potential to reduce the abdominal organ motion while improving both target tracking reliability and motion reproducibility. Furthermore, the observed correlation between liver and pancreas/stomach motion indicates the possibility of indirect pancreas/stomach tracking using liver markers as tracking surrogates. The strategy is expected to provide an alternative for respiratory motion management in the radiation treatment of abdominal tumors. This article is protected by copyright. All rights reserved.

Dosimetric evaluation of respiratory gating on a 0.35-T magnetic resonance-guided radiotherapy linac

PURPOSE

The commercial 0.35-T magnetic resonance imaging (MRI)-guided radiotherapy vendor ViewRay recently introduced upgraded real-time imaging frame rates based on compressed sensing techniques. Furthermore, additional motion tracking algorithms were made available. Compressed sensing allows for increased image frame rates but may compromise image quality. To assess the impact of this upgrade on respiratory gating accuracy, we evaluated gated dose distributions pre- and post-upgrade using a motion phantom and radiochromic film.

METHODS

Seven motion waveforms (four artificial, two patient-derived free-breathing, and one breath-holding) were used to drive an MRI-compatible motion phantom. A treatment plan was developed to deliver a 3-cm diameter spherical dose distribution typical of a stereotactic body radiotherapy plan. Gating was performed using 4-frames per second (fps) imaging pre-upgrade on the "default" tracking algorithm and 8-fps post-upgrade using the "small mobile targets" (SMT) and "large deforming targets" (LDT) tracking algorithms. Radiochromic film was placed in a moving insert within the phantom to measure dose. The planned and delivered dose distributions were compared using the gamma index with 3%/3-mm criteria. Dose-area histograms were produced to calculate the dose to 95% (D95) of the sphere planning target volume (PTV) and two simulated gross tumor volumes formed by contracting the PTV by 3 and 5 mm, respectively.

RESULTS

Gamma pass rates ranged from 18% to 93% over the 21 combinations of breathing trace and gating conditions examined. D95 ranged from 206 to 514 cGy. On average, the LDT algorithm yielded lower gamma and D95 values than the default and SMT algorithms.

CONCLUSION

Respiratory gating at 8 fps with the new tracking algorithms provides similar gating performance to the original algorithm with 4 fps, although the LDT algorithm had lower accuracy for our non-deformable target. This indicates that the choice of deformable image registration algorithm should be chosen deliberately based on whether the target is rigid or deforming.

Radiotherapy respiratory motion management in hepatobiliary and pancreatic malignancies: a systematic review of patient factors influencing effectiveness of motion reduction with abdominal compression

BACKGROUND

The effectiveness of abdominal compression for motion management in hepatobiliary-pancreatic (HPB) radiotherapy has not been systematically evaluated.

METHODS & MATERIALS

A systematic review was carried out using PubMed/Medline, Cochrane Library, Web of Science, and CINAHL databases up to 1 July 2021. No date restrictions were applied. Additional searches were carried out using the University of Manchester digital library, Google Scholar and of retrieved papers' reference lists. Studies conducted evaluating respiratory motion utilising imaging with and without abdominal compression in the same patients available in English were included. Studies conducted in healthy volunteers or majority non-HPB sites, not providing descriptive motion statistics or patient characteristics before and after compression in the same patients or published without peer-review were excluded. A narrative synthesis was employed by tabulating retrieved studies and organising chronologically by abdominal compression device type to help identify patterns in the evidence.

RESULTS

The inclusion criteria were met by 6 studies with a total of 152 patients. Designs were a mix of retrospective and prospective quantitative designs with chronological, non-randomised recruitment. Abdominal compression reduced craniocaudal respiratory motion in the majority of patients, although in four studies there were increases seen in at least one direction. The influence of patient comorbidities on effectiveness of compression, and/or comfort with compression was not evaluated in any study.

CONCLUSION

Abdominal compression may not be appropriate for all patients, and benefit should be weighed with potential increase in motion or discomfort in patients with small initial motion (<5 mm). Patient factors including male sex, and high body mass index (BMI) were found to impact the effectiveness of compression, however with limited evidence. High-quality studies are warranted to fully assess the clinical impact of abdominal compression on treatment outcomes and toxicity prospective in comparison to other motion management strategies.

A standardized workflow for respiratory-gated motion management decision-making

Purpose

Motion management of tumors within the lung and abdomen is challenging because it requires balancing tissue sparing with accuracy of hitting the target, while considering treatment delivery efficiency. Physicists can play an important role in analyzing four‐dimensional computed tomography (4DCT) data to recommend the optimal respiratory gating parameters for a patient. The goal of this work was to develop a standardized procedure for making recommendations regarding gating parameters and planning margins for lung and gastrointestinal stereotactic body radiotherapy (SBRT) treatments. In doing so, we hoped to simplify decision‐making and analysis, and provide a tool for troubleshooting complex cases.

Methods

Factors that impact gating decisions and planning target volume (PTV) margins were identified. The gating options included gating on exhale with approximately a 50% duty cycle (Gate3070), exhale gating with a reduced duty cycle (Gate4060), and treating for most of respiration, excluding only extreme inhales and exhales (Gate100). A standard operating procedure was developed, as well as a physics consult document to communicate motion management recommendations to other members of the treatment team. This procedure was implemented clinically for 1 year and results are reported below.

Results

Identified factors that impact motion management included the magnitude of motion observed on 4DCT, the regularity of breathing and quality of 4DCT data, and ability to observe the target on fluoroscopy. These were collated into two decision tables—one specific to lung tumors and another for gastrointestinal tumors—such that a physicist could answer a series of questions to determine the optimal gating and PTV margin. The procedure was used clinically for 252 sites from 213 patients treated with respiratory‐gated SBRT and standardized practice across our 12‐member physics team.

Conclusion

Implementation of a standardized procedure for respiratory gating had a positive impact in our clinic, improving efficiency and ease of 4DCT analysis and standardizing gating decision‐making amongst physicists.

Management of Motion and Anatomical Variations in Charged Particle Therapy: Past, Present, and Into the Future

The major aim of radiation therapy is to provide curative or palliative treatment to cancerous malignancies while minimizing damage to healthy tissues. Charged particle radiotherapy utilizing carbon ions or protons is uniquely suited for this task due to its ability to achieve highly conformal dose distributions around the tumor volume. For these treatment modalities, uncertainties in the localization of patient anatomy due to inter- and intra-fractional motion present a heightened risk of undesired dose delivery. A diverse range of mitigation strategies have been developed and clinically implemented in various disease sites to monitor and correct for patient motion, but much work remains. This review provides an overview of current clinical practices for inter and intra-fractional motion management in charged particle therapy, including motion control, current imaging and motion tracking modalities, as well as treatment planning and delivery techniques. We also cover progress to date on emerging technologies including particle-based radiography imaging, novel treatment delivery methods such as tumor tracking and FLASH, and artificial intelligence and discuss their potential impact towards improving or increasing the challenge of motion mitigation in charged particle therapy.

Dosimetry impact of gating latency in cine magnetic resonance image guided breath-hold pancreatic cancer radiotherapy

Objective.We investigated dosimetry effect of gating latency in cine magnetic resonance image (cine MRI) guided breath-hold pancreatic cancer radiotherapy.Approach.The gating latency was calculated based on cine MRI obtained from 17 patients who received MRI guided radiotherapy. Because of the cine MRI-related latency, beam overshoot occurs when beam remains on while the tracking target already moves out of the target boundary. The number of beam on/off events was calculated from the cine MRI data. We generated both IMRT and VMAT plans for all 17 patients using 33 Gy prescription, and created motion plans by applying isocenter shift that corresponds to motion-induced tumor displacement. The GTV and PTV coverage and dose to nearby critical structures were compared between the motion and original plan to evaluate the dosimetry change caused by cine MRI latency.Main results.The time ratio of cine MRI imaging latency over the treatment duration is 6.6 ± 3.1%, the mean and median percentage of beam-on events <4 s are 67.0 ± 14.3% and 66.6%. When a gating boundary of 4 mm and a target-out threshold of 5% is used, there is no significant difference for GTV V33Gy between the motion and original plan (p = 0.861 and 0.397 for IMRT and VMAT planning techniques, respectively). However, the PTV V33Gy and stomach Dmax for the motion plans are significantly lower; duodenum V12.5 Gy and V18Gy are significantly higher when compared with the original plans, for both IMRT and VMAT planning techniques.Significance.The cine MRI gating latency can significantly decrease the dose delivered to the PTV, and increase the dose to the nearby critical structures. However, no significant difference is observed for the GTV coverage. The dosimetry impact can be mitigated by implementing additional beam-on control techniques which reduces unnecessary beam on events and/or by using faster cine MRI sequences which reduces the latency period.

A Dual-supervised Deformation Estimation Model (DDEM) for constructing ultra-quality 4D-MRI based on a commercial low-quality 4D-MRI for liver cancer radiation therapy

BACKGROUND

Most available 4D-MRI techniques are limited by insufficient image quality and long acquisition times or require specially designed sequences or hardware that are not available in the clinic. These limitations have greatly hindered the clinical implementation of 4D-MRI.

PURPOSE

This study aims to develop a fast ultra-quality (UQ) 4D-MRI reconstruction method using a commercially available 4D-MRI sequence and dual-supervised deformation estimation model (DDEM).

METHODS

Thirty-nine patients receiving radiotherapy for liver tumors were included. Each patient was scanned using a TWIST-VIBE MRI sequence to acquire 4D-MR images. They also received 3D T1-/T2-weighted MRI scans as prior images and UQ 4D-MRI at any instant was considered a deformation of them. A DDEM was developed to obtain a 4D deformable vector field (DVF) from 4D-MRI data, and the prior images were deformed using this 4D-DVF to generate UQ 4D-MR images. The registration accuracies of the DDEM, VoxelMorph (normalized cross-correlation (NCC) supervised), VoxelMorph (end-to-end point error (EPE) supervised), and the parametric total variation (pTV) algorithm were compared. Tumor motion on UQ 4D-MRI was evaluated quantitatively using region-of-interest (ROI) tracking errors, while image quality was evaluated using the contrast-to-noise ratio (CNR), lung-liver edge sharpness, and perceptual blur metric (PBM).

RESULTS

The registration accuracy of the DDEM was significantly better than those of VoxelMorph (NCC supervised), VoxelMorph (EPE supervised) and the pTV algorithm (all, p < 0.001), with an inference time of 69.3 ± 5.9 ms. UQ 4D-MRI yielded ROI tracking errors of 0.79 ± 0.65, 0.50 ± 0.55, and 0.51 ± 0.58 mm in the superior-inferior, anterior-posterior, and mid-lateral directions, respectively. From the original 4D-MRI to UQ 4D-MRI, the CNR increased from 7.25 ± 4.89 to 18.86 ± 15.81; the lung-liver edge full-width-at-half-maximum decreased from 8.22 ± 3.17 to 3.65 ± 1.66 mm in the in-plane direction and from 8.79 ± 2.78 to 5.04 ± 1.67 mm in the cross-plane direction, and the PBM decreased from 0.68 ± 0.07 to 0.38 ± 0.01.

CONCLUSION

This novel DDEM method successfully generated UQ 4D-MR images based on a commercial 4D-MRI sequence. It shows great promise for improving liver tumor motion management during radiation therapy. This article is protected by copyright. All rights reserved.

A transit portal dosimetry method for respiratory gating quality assurance with a dynamic 3D printed tumor phantom

Backgrounds

Respiratory gating is one of the motion management techniques that is used to deliver radiation dose to a tumor at a specific position under free breathing. However, due to the dynamic feedback process of this approach, regular equipment quality assurance (QA) and patient‐specific QA checks need to be performed. This work proposes a new QA methodology using electronic portal imaging detector (EPID) to determine the target localization accuracy of phase gating.

Methods

QA tools comprising 3D printed spherical tumor phantoms, programmable stages, and an EPID detector are characterized and assembled. Algorithms for predicting portal dose (PD) through moving phantoms are developed and verified using gamma analysis for two spherical tumor phantoms (2 cm and 4 cm), two different 6 MV volumetric modulated arc therapy plans, and two different gating windows (30%–70% and 40%–60%). Comparison between the two gating windows is then performed using the Wilcoxon signed‐rank test. An optimizer routine, which is used to determine the optimal window, based on maximal gamma passing rate (GPR), was applied to an actual breathing curve and breathing plan. This was done to ascertain if our method yielded a similar result with the actual gating window.

Results

High GPRs of more than 97% and 91% were observed when comparing the predicted PD with the measured PD in moving phantom at 2 mm/2% and 1 mm/1% levels, respectively. Analysis of gamma heatmaps shows an excellent agreement with the tumor phantom.

The GPR of 40%–60% PD was significantly lower than that of the 30%–70% PD at the 1 mm/1% level (p = 0.0064). At the 2 mm/2% level, no significant differences were observed. The optimizer routine could accurately predict the center of the gating window to within a 10% range.

Conclusion

We have successfully performed and verified a new method for QA with the use of a moving phantom with EPID for phase gating with real‐time position management.

Report of AAPM Task Group 290: Respiratory motion management for particle therapy

Dose uncertainty induced by respiratory motion remains a major concern for treating thoracic and abdominal lesions using particle beams. This Task Group report reviews the impact of tumor motion and dosimetric considerations in particle radiotherapy, current motion‐management techniques, and limitations for different particle‐beam delivery modes (i.e., passive scattering, uniform scanning, and pencil‐beam scanning). Furthermore, the report provides guidance and risk analysis for quality assurance of the motion‐management procedures to ensure consistency and accuracy, and discusses future development and emerging motion‐management strategies. This report supplements previously published AAPM report TG76, and considers aspects of motion management that are crucial to the accurate and safe delivery of particle‐beam therapy. To that end, this report produces general recommendations for commissioning and facility‐specific dosimetric characterization, motion assessment, treatment planning, active and passive motion‐management techniques, image guidance and related decision‐making, monitoring throughout therapy, and recommendations for vendors. Key among these recommendations are that: (1) facilities should perform thorough planning studies (using retrospective data) and develop standard operating procedures that address all aspects of therapy for any treatment site involving respiratory motion; (2) a risk‐based methodology should be adopted for quality management and ongoing process improvement.

An Evaluation of Total Internal Motions of Locally Advanced Pancreatic Cancer during SABR Using Calypso® Extracranial Tracking, and Its Possible Clinical Impact on Motion Management

(1) Background: the aims of this study were to determine the total extent of pancreatic cancer’s internal motions, using Calypso^® extracranial tracking, and to indicate possible clinical advantages of continuous intrafractional fiducial-based tumor motion tracking during SABR. (2) Methods: thirty-four patients were treated with SABR for LAPC using Calypso^® for motion management. Planning MSCTs in FB and DBH, and 4D-CTs were performed. Using data from Calypso^® and 4D-CTs, the movements of the lesions in the CC, AP and LR directions, as well as the volumes of the 4D-CT-based ITV and the volumes of the Calypso^®-based ITV were compared. (3) Results: significantly larger medians of tumor excursions were found with Calypso^® than with 4D-CT: CC: 29 mm (p < 0.001); AP: 14 mm (p < 0.001) and LR: 11 mm (p < 0.039). The median volume of the Calypso^®-based ITV was significantly larger than that of the 4D-CT based ITV (p < 0.001). (4) Conclusion: beside known respiratory-induced internal motions, pancreatic cancer seems to have significant additional motions which should be considered during respiratory motion management. Only direct and continuous intrafractional fiducial-based motion tracking seems to provide complete coverage of the target lesion with the prescribed isodose, which could allow for safe tumor dose escalation.

Application of an automatic, uncertainty model-guided, target-generating algorithm to lung stereotactic body radiotherapy

PURPOSE

This work evaluated a new radiotherapy target-generating framework (the αTarget algorithm) for creating internal target volumes for lung SBRT.

METHODS

Nineteen patients previously treated with definitive intent SBRT to the lung were identified from a clinical database. For each patient's 4DCT simulation scan, deformable image registration was used between phases of the scan in order to generate voxelized models of motion for 35 individual gross tumor volumes. These motion models were then used with a new implementation of a previously described target-generating algorithm to create new internal target volumes (αITVs). The resulting αITVs were analyzed with respect to their volume and the coverage they provided each tumor voxel per that voxel's motion model. The clinically used ITVs were similarly analyzed, and were then compared to the αITVs using paired Student's t-tests. In addition, isotropic margins were added to the αITVs in order to determine the largest margin magnitude that could be added without exceeding the volume of the clinical ITVs.

RESULTS

The αITVs increased the target coverage provided to each tumor's 5th-percentile-most-covered-voxel an average of 50.3% compared to the clinical ITVs (p < 0.0001). At the same time, the αITVs had volumes that were, on average, 31.4% smaller (p < 0.0001). The differences in volume were large enough that, on average, an extra 2 mm isotropic margin could be added to the αITV before it had a volume greater than the clinical ITV.

CONCLUSIONS

The αTarget algorithm can generate more effective lung SBRT internal target volumes that provide greater coverage with smaller volumes. In combination with numerous other advantages of the framework, this effectiveness makes the αTarget algorithm a powerful new method for advanced IGRT or adaptive radiotherapy techniques.

Personalized treatment gating thresholds in frameless stereotactic radiosurgery using predictions of dosimetric fidelity and treatment interruption

PURPOSE

Gamma Knife Icon (GKI) enables a user-defined gating threshold for intrafraction motion during stereotactic radiosurgery (SRS). An optimal threshold would ensure dosimetric fidelity of the planned distribution and minimize treatment time extension by gating. A prediction of motion characteristics for a patient based on a retrospective database of motion traces could be beneficial to evaluating the choice of gating threshold. A short acquisition of motion may help to define a personalized threshold that balances dosimetric accuracy and treatment length. This study aims to evaluate the performance of a prediction of motion and the resultant dosimetric consequences for a range of motion gating thresholds.

METHODS

A database of 2552 motion traces (776 patients) was analyzed using previously published methods to characterize patient intrafraction motion on the GKI. For a selection of six fractionated SRS patient cases (two patients with single brain metastasis, four vestibular schwannomas), a 10-min sample of motion was used to classify motion and identify traces in the database with similar metrics. The similar motion traces were used to perform a predictive reconstruction of the selected patient's delivered dose for a range of motion thresholds. The remaining fractions were reconstructed and compared to that predicted. From the six cases, 26 fractions were used to predict the number of interruptions (n = 26), change in target coverage (n = 26), and change in brainstem maximum dose (vestibular cases only, n = 20). The difference between mean predicted and reconstructed values was compared for accuracy.

RESULTS

The difference between mean prediction and reconstructed values was 0.32 ± 0.38% in target coverage, 2.36 ± 5.06 interruptions, and 0.15 ± 0.24 Gy for the brainstem maximum dose. Sixty-seven of the 72 predictions (26 coverage, 26 interruptions, and 20 brainstem maximum dose) were within one standard deviation of the predicted mean.

CONCLUSIONS

Large databases of motion traces were used to characterize patient performance and predict motion performance. Dosimetric deterioration due to motion and extension of treatment duration can be predicted in some cases using only a short acquisition of motion and the treatment plan. This reconstruction may provide benefit in generating a patient-specific motion threshold.

Potential Morbidity Reduction for Lung Stereotactic Body Radiation Therapy Using Respiratory Gating

Simple Summary

Lung stereotactic body radiotherapy (SBRT) is the standard of care for early-stage lung cancer and oligometastases. For SBRT, motion has to be considered to avoid misdosage. Respiratory phase gating, meaning to irradiate the target volume only in a predefined gating motion phase window, can be applied to mitigate motion-induced effects. The aim of this study was to exploit the clinical benefit of gating for lung SBRT. For the majority of 14 lung tumor patients and various gating windows, we could prove a reduced dose to normal tissue by gating simulation. A normal tissue complication probability (NTCP) model analysis revealed a major reduction of normal tissue toxicity for moderate gating window sizes. The most beneficial effect of gating was found for those patients with the highest prior toxicity risk. The presented results are useful for personalized risk assessment prior to treatment and may help to select patients and optimal gating windows.

Abstract

We investigated the potential of respiratory gating to mitigate the motion-caused misdosage in lung stereotactic body radiotherapy (SBRT). For fourteen patients with lung tumors, we investigated treatment plans for a gating window (GW) including three breathing phases around the maximum exhalation phase, GW40–60. For a subset of six patients, we also assessed a preceding three-phase GW20–40 and six-phase GW20–70. We analyzed the target volume, lung, esophagus, and heart doses. Using normal tissue complication probability (NTCP) models, we estimated radiation pneumonitis and esophagitis risks. Compared to plans without gating, GW40–60 significantly reduced doses to organs at risk without impairing the tumor doses. On average, the mean lung dose decreased by 0.6 Gy (p < 0.001), treated lung V20Gy by 2.4% (p = 0.003), esophageal dose to 5cc by 2.0 Gy (p = 0.003), and maximum heart dose by 3.2 Gy (p = 0.009). The model-estimated mean risks of 11% for pneumonitis and 12% for esophagitis without gating decreased upon GW40–60 to 7% and 9%, respectively. For the highest-risk patient, gating reduced the pneumonitis risk from 43% to 32%. Gating is most beneficial for patients with high-toxicity risks. Pre-treatment toxicity risk assessment may help optimize patient selection for gating, as well as GW selection for individual patients.

Initial analysis of the dosimetric benefit and clinical resource cost of CBCT-based online adaptive radiotherapy for patients with cancers of the cervix or rectum

Purpose

This provides a benchmark of dosimetric benefit and clinical cost of cone‐beam CT‐based online adaptive radiotherapy (ART) technology for cervical and rectal cancer patients.

Methods

An emulator of a CBCT‐based online ART system was used to simulate more than 300 treatments for 13 cervical and 15 rectal cancer patients. CBCT images were used to generate adaptive replans. To measure clinical resource cost, the six phases of the workflow were timed. To evaluate the dosimetric benefit, changes in dosimetric values were assessed. These included minimum dose (Dmin) and volume receiving 95% of prescription (V95%) for the planning target volume (PTV) and the clinical target volume (CTV), and maximum 2 cc's (D2cc) of the bladder, bowel, rectum, and sigmoid colon.

Results

The average duration of the workflow was 24.4 and 9.2 min for cervical and rectal cancer patients, respectively. A large proportion of time was dedicated to editing target contours (13.1 and 2.7 min, respectively). For cervical cancer patients, the replan changed the Dmin to the PTVs and CTVs for each fraction 0.25 and 0.25 Gy, respectively. The replan changed the V95% by 9.2 and 7.9%. The D2cc to the bladder, bowel, rectum, and sigmoid colon for each fraction changed −0.02, −0.08, −0.07, and −0.04 Gy, respectively. For rectal cancer patients, the replan changed the Dmin to the PTVs and CTVs for each fraction of 0.20 and 0.24 Gy, respectively. The replan changed the V95% by 4.1 and 1.5%. The D2cc to the bladder and bowel for each fraction changed 0.02 and −0.02 Gy, respectively.

Conclusions

Dosimetric benefits can be achieved with CBCT‐based online ART that is amenable to conventional appointment slots. The clinical significance of these benefits remains to be determined. Managing contours was the primary factor affecting the total duration and is imperative for safe and effective adaptive radiotherapy.

Technical Note: Four-dimensional deformable digital phantom for MRI sequence development

PURPOSE

MR-guided radiotherapy has different requirements for the images than diagnostic radiology, thus requiring development of novel imaging sequences. MRI simulation is an excellent tool for optimizing these new sequences; however, currently available software does not provide all the necessary features. In this paper, we present a digital framework for testing MRI sequences that incorporates anatomical structure, respiratory motion, and realistic presentation of MR physics.

METHODS

The extended Cardiac-Torso (XCAT) software was used to create T1 , T2 , and proton density maps that formed the anatomical structure of the phantom. Respiratory motion model was based on the XCAT deformation vector fields, modified to create a motion model driven by a respiration signal. MRI simulation was carried out with JEMRIS, an open source Bloch simulator. We developed an extension for JEMRIS, which calculates the motion of each spin independently, allowing for deformable motion.

RESULTS

The performance of the framework was demonstrated through simulating the acquisition of a two-dimensional (2D) cine and demonstrating expected motion ghosts from T2 weighted spin echo acquisitions with different respiratory patterns. All simulations were consistent with behavior previously described in literature. Simulations with deformable motion were not more time consuming than with rigid motion.

CONCLUSIONS

We present a deformable four-dimensional (4D) digital phantom framework for MR sequence development. The framework incorporates anatomical structure, realistic breathing patterns, deformable motion, and Bloch simulation to achieve accurate simulation of MRI. This method is particularly relevant for testing novel imaging sequences for the purpose of MR-guided radiotherapy in lungs and abdomen.

Validation of proton dose calculation on scatter corrected 4D cone beam computed tomography using a porcine lung phantom

Proton therapy treatment for lungs remains challenging as images enabling the detection of inter- and intra-fractional motion, which could be used for proton dose adaptation, are not readily available. 4D computed tomography (4DCT) provides high image quality but is rarely available in-room, while in-room 4D cone beam computed tomography (4DCBCT) suffers from image quality limitations stemming mostly from scatter detection. This study investigated the feasibility of using virtual 4D computed tomography (4DvCT) as a prior for a phase-per-phase scatter correction algorithm yielding a 4D scatter corrected cone beam computed tomography image (4DCBCT_cor), which can be used for proton dose calculation. 4DCT and 4DCBCT scans of a porcine lung phantom, which generated reproducible ventilation, were acquired with matching breathing patterns. Diffeomorphic Morphons, a deformable image registration algorithm, was used to register the mid-position 4DCT to the mid-position 4DCBCT and yield a 4DvCT. The 4DCBCT was reconstructed using motion-aware reconstruction based on spatial and temporal regularization (MA-ROOSTER). Successively for each phase, digitally reconstructed radiographs of the 4DvCT, simulated without scatter, were exploited to correct scatter in the corresponding CBCT projections. The 4DCBCT_corwas then reconstructed with MA-ROOSTER using the corrected CBCT projections and the same settings and deformation vector fields as those already used for reconstructing the 4DCBCT. The 4DCBCT_corand the 4DvCT were evaluated phase-by-phase, performing proton dose calculations and comparison to those of a ground truth 4DCT by means of dose-volume-histograms (DVH) and gamma pass-rates (PR). For accumulated doses, DVH parameters deviated by at most 1.7% in the 4DvCT and 2.0% in the 4DCBCT_corcase. The gamma PR for a (2%, 2 mm) criterion with 10% threshold were at least 93.2% (4DvCT) and 94.2% (4DCBCT_cor), respectively. The 4DCBCT_cortechnique enabled accurate proton dose calculation, which indicates the potential for applicability to clinical 4DCBCT scans.

Surface guided frameless positioning for lung stereotactic body radiation therapy

Background and purpose

When treating lung tumors with stereotactic body radiation therapy (SBRT), patient immobilization is of outmost importance. In this study, the intra‐fractional shifts of the patient (based on bony anatomy) and the tumor (based on the visible target volume) are quantified, and the associated impact on the delivered dose is estimated for a frameless immobilization approach in combination with surface guided radiation therapy (SGRT) monitoring.

Methods

Cone beam computed tomographies (CBCT) were collected in free breathing prior and after each treatment for 25 patients with lung tumors, in total 137 fractions. The CBCT collected after each treatment was registered to the CBCT collected before each treatment with focus on bony anatomy to determine the shift of the patient, and with focus on the visible target volume to determine the shift of the tumor. Rigid registrations with 6 degrees of freedom were used. The patients were positioned in frameless immobilizations with their position and respiration continuously monitored by a commercial SGRT system. The patients were breathing freely within a preset gating window during treatment delivery. The beam was automatically interrupted if isocenter shifts >4 mm or breathing amplitudes outside the gating window were detected by the SGRT system. The time between the acquisition of the CBCTs was registered for each fraction to examine correlations between treatment time and patient shift. The impact of the observed shifts on the dose to organs at risk (OAR) and the gross tumor volume (GTV) was assessed.

Results

The shift of the patient in the CBCTs was ≤2 mm for 132/137 fractions in the vertical (vrt) and lateral (lat) directions, and 134/137 fractions in the longitudinal (lng) direction and ≤4 mm in 134/137 (vrt) and 137/137 (lat, lng) of the fractions. The shift of the tumor was ≤2 mm in 116/137 (vrt), 123/137 (lat) and 115/137 (lng) fractions and ≤4 mm in 136/137 (vrt), 137/137 (lat), and 135/137 (lng) fractions. The maximal observed shift in the evaluated CBCT data was 4.6 mm for the patient and 7.2 mm for the tumor. Rotations were ≤3.3ᵒ for all fractions and the mean/standard deviation were 0.2/1.0ᵒ (roll), 0.1/0.8ᵒ (yaw), and 0.3/1.0ᵒ (pitch). The SGRT system interrupted the beam due to intra‐fractional isocenter shifts >4 mm for 21% of the fractions, but the patients always returned within tolerance without the need of repositioning. The maximal observed isocenter shift by the SGRT system during the beam holds was 8 mm. For the respiration monitoring, the beam was interrupted at least one time for 54% of the fractions.

The visual tumor was within the planned internal target volume (ITV) for 136/137 fractions in the evaluated CBCT data collected at the end of each fraction. For the fraction where the tumor was outside the ITV, the D_98% for the GTV decreased with 0.4 Gy. For the OARs, the difference between planned and estimated dose from the CBCT data (D_2% or D_mean) was ≤2.6% of the prescribed PTV dose. No correlation was found between treatment time and the magnitude of the patient shift.

Conclusions

Using SGRT for motion management and respiration monitoring in combination with a frameless immobilization is a feasible approach for lung SBRT.

Transmission study of the Abdominal Compression plate (BodyFIX Diaphragm Control) for abdominal and stereotactic body radiotherapy

Purpose

Abdominal Compression is one of the methods available to minimize breathing motion during stereotactic body radiotherapy (SBRT), particularly for abdominal malignancies. It might be necessary to treat some tumors with radiation entering through the compression device. One clinically available compression plate device (Elekta BodyFIX Diaphragm Control) was evaluated to understand its impact on dosimetry during clinical treatments.

Methods

The BodyFIX compression device was CT scanned following departmental stereo scanning protocols. Treatment planning system (TPS) calculations were used to determine attenuation ratios through each section of the compression device: the outer frame, compression plate, and higher density couch fixation points and compression screw. TPS calculated skin doses where the compression plate will come in contact with the skin were recorded. All attenuation ratio fields were measured on an Elekta Versa HD linear accelerator. Where differences in attenuation were observed, TPS density overrides were found to bring calculated doses into agreement with measurement.

Results

The compression plate and frame showed low dose attenuation (3%–4%). Only minor density overrides for the frame were required due to artefacts from the limited CT field‐of‐view. The high‐density materials in the couch fixation points resulted in higher attenuation (14%–20%). Similarly, the compression screw recorded very high attenuation (44%–65%), depending on the length of screw used. Skin doses assessed from the TPS calculations showed dose build‐up under the compression plate that would result in skin receiving the maximum dose.

Conclusion

Compression devices can cause significant dose attenuation. Density overrides for TPS calculations are recommended for correcting attenuation in some sections of the device. High‐density structures like the fixation screw and frame fixation points create high levels of dosimetric uncertainty, and beam entry through those areas has been disallowed.

Real-time CT image generation based on voxel-by-voxel modeling of internal deformation by utilizing the displacement of fiducial markers

PURPOSE

To show the feasibility of real-time CT image generation technique utilizing internal fiducial markers that facilitate the evaluation of internal deformation.

METHODS

In the proposed method, a linear regression model that can derive internal deformation from the displacement of fiducial markers is built for each voxel in the training process before the treatment session. Marker displacement and internal deformation are derived from the four-dimensional computed tomography (4DCT) dataset. In the treatment session, the three-dimensional deformation vector field is derived according to the marker displacement, which is monitored by the real-time imaging system. The whole CT image can be synthesized by deforming the reference CT image with a deformation vector field in real-time. To show the feasibility of the technique, image synthesis accuracy and tumor localization accuracy were evaluated using the dataset generated by extended NURBS-Based Cardiac-Torso (XCAT) phantom and clinical 4DCT datasets from six patients, containing 10 CT datasets each. In the validation with XCAT phantom, motion range of the tumor in training data and validation data were about 10 and 15 mm, respectively, so as to simulate motion variation between 4DCT acquisition and treatment session. In the validation with patient 4DCT dataset, eight CT datasets from the 4DCT dataset were used in the training process. Two excluded inhale CT datasets can be regarded as the datasets with large deformations more than training dataset. CT images were generated for each respiratory phase using the corresponding marker displacement. Root mean squared error (RMSE), normalized RMSE (NRMSE), and structural similarity index measure (SSIM) between the original CT images and the synthesized CT images were evaluated as the quantitative indices of the accuracy of image synthesis. The accuracy of tumor localization was also evaluated.

RESULTS

In the validation with XCAT phantom, the mean NRMSE, SSIM, and three-dimensional tumor localization error were 7.5 ± 1.1%, 0.95 ± 0.02, and 0.4 ± 0.3 mm, respectively. In the validation with patient 4DCT dataset, the mean RMSE, NRMSE, SSIM, and three-dimensional tumor localization error in six patients were 73.7 ± 19.6 HU, 9.2 ± 2.6%, 0.88 ± 0.04, and 0.8 ± 0.6 mm, respectively. These results suggest that the accuracy of the proposed technique is adequate when the respiratory motion is within the range of the training dataset. In the evaluation with a marker displacement larger than that of the training dataset, the mean RMSE, NRMSE, and tumor localization error were about 100 HU, 13%, and <2.0 mm, respectively, except for one case having large motion variation. The performance of the proposed method was similar to those of previous studies. Processing time to generate the volumetric image was <100 ms.

CONCLUSION

We have shown the feasibility of the real-time CT image generation technique for volumetric imaging.

Validation of dose distribution for liver tumors treated with real-time-image gated spot-scanning proton therapy by log data based dose reconstruction

In spot scanning proton therapy (SSPT), the spot position relative to the target may fluctuate through tumor motion even when gating the radiation by utilizing a fiducial marker. We have established a procedure that evaluates the delivered dose distribution by utilizing log data on tumor motion and spot information. The purpose of this study is to show the reliability of the dose distributions for liver tumors treated with real-time-image gated SSPT (RGPT). In the evaluation procedure, the delivered spot information and the marker position are synchronized on the basis of log data on the timing of the spot irradiation and fluoroscopic X-ray irradiation. Then a treatment planning system reconstructs the delivered dose distribution. Dose distributions accumulated for all fractions were reconstructed for eight liver cases. The log data were acquired in all 168 fractions for all eight cases. The evaluation was performed for the values of maximum dose, minimum dose, D99, and D5-D95 for the clinical target volumes (CTVs) and mean liver dose (MLD) scaled by the prescribed dose. These dosimetric parameters were statistically compared between the planned dose distribution and the reconstructed dose distribution. The mean difference of the maximum dose was 1.3% (95% confidence interval [CI]: 0.6%-2.1%). Regarding the minimum dose, the mean difference was 0.1% (95% CI: -0.5%-0.7%). The mean differences of D99, D5-D95 and MLD were below 1%. The reliability of dose distributions for liver tumors treated with RGPT-SSPT was shown by the evaluation of the accumulated dose distributions.

Design of an anthropomorphic PET phantom with elastic lungs and respiration modeling

PURPOSE

Respiratory motion during positron emission tomography (PET) scans can be a major detriment to image quality in oncological imaging. The impact of motion on lesion quantification and detectability can be assessed using phantoms with realistic anatomy representation and motion modeling. In this work, we develop an anthropomorphic phantom for PET imaging that combines anatomic fidelity and a realistic breathing mechanism with deformable lungs.

METHODS

We start from a previously developed anatomically accurate but static phantom of a human torso, and add elastic lungs with a highly controllable actuation mechanism which replicates the physics of breathing. The space outside the lungs is filled with a radioactive water solution. To maintain anatomical accuracy and realistic gamma ray attenuation in the torso, all motion mechanisms and actuators are positioned outside of the phantom compartment. The actuation mechanism can produce custom respiratory waveforms with breathing rates up to 25 breaths per minute and tidal volumes up to 1200 mL.

RESULTS

Several tests were performed to validate the performance of the phantom assembly, in which the phantom was filled with water and given respiratory waveforms to execute. All parts demonstrated expected performance. Force requirements were not exceeded and no leaks were detected, although continued use of the phantom is required to evaluate wear. The motion of the lungs was determined to be within a reasonable realistic range.

CONCLUSIONS

The full mechanical design is described in this paper, as well as a software application with graphical user interface which was developed to plan and visualize respiratory patterns. Both are available online as open source files. The developed phantom will facilitate future work in evaluating the impact of respiratory motion on lesion quantification and detectability in clinical practice.

On the evaluation of mobile target trajectory between four-dimensional computer tomography and four-dimensional cone-beam computer tomography

Purpose

For mobile lung tumors, four‐dimensional computer tomography (4D CT) is often used for simulation and treatment planning. Localization accuracy remains a challenge in lung stereotactic body radiation therapy (SBRT) treatments. An attractive image guidance method to increase localization accuracy is 4D cone‐beam CT (CBCT) as it allows for visualization of tumor motion with reduced motion artifacts. However, acquisition and reconstruction of 4D CBCT differ from that of 4D CT. This study evaluates the discrepancies between the reconstructed motion of 4D CBCT and 4D CT imaging over a wide range of sine target motion parameters and patient waveforms.

Methods

A thorax motion phantom was used to examine 24 sine motions with varying amplitudes and cycle times and seven patient waveforms. Each programmed motion was imaged using 4D CT and 4D CBCT. The images were processed to auto segment the target. For sine motion, the target centroid at each phase was fitted to a sinusoidal curve to evaluate equivalence in amplitude between the two imaging modalities. The patient waveform motion was evaluated based on the average 4D data sets.

Results

The mean difference and root‐mean‐square‐error between the two modalities for sine motion were −0.35 ± 0.22 and 0.60 mm, respectively, with 4D CBCT slightly overestimating amplitude compared with 4D CT. The two imaging methods were determined to be significantly equivalent within ±1 mm based on two one‐sided t tests (p < 0.001). For patient‐specific motion, the mean difference was 1.5 ± 2.1 (0.8 ± 0.6 without outlier), 0.4 ± 0.3, and 0.8 ± 0.6 mm for superior/inferior (SI), anterior/posterior (AP), and left/right (LR), respectively.

Conclusion

In cases where 4D CT is used to image mobile tumors, 4D CBCT is an attractive localization method due to its assessment of motion with respect to 4D CT, particularly for lung SBRT treatments where accuracy is paramount.