An interdimensional correlation framework for real-time estimation of six degree of freedom target motion using a single x-ray imager during radiotherapy

© 2017 Institute of Physics and Engineering in Medicine. Increasing evidence suggests that intrafraction tumour motion monitoring needs to include both 3D translations and 3D rotations. Presently, methods to estimate the rotation motion require the 3D translation of the target to be known first. However, ideally, translation and rotation should be estimated concurrently. We present the first method to directly estimate six-degree-of-freedom (6DoF) motion from the target's projection on a single rotating x-ray imager in real-time. This novel method is based on the linear correlations between the superior-inferior translations and the motion in the other five degrees-of-freedom. The accuracy of the method was evaluated in silico with 81 liver tumour motion traces from 19 patients with three implanted markers. The ground-truth motion was estimated using the current gold standard method where each marker's 3D position was first estimated using a Gaussian probability method, and the 6DoF motion was then estimated from the 3D positions using an iterative method. The 3D position of each marker was projected onto a gantry-mounted imager with an imaging rate of 11 Hz. After an initial 110° gantry rotation (200 images), a correlation model between the superior-inferior translations and the five other DoFs was built using a least square method. The correlation model was then updated after each subsequent frame to estimate 6DoF motion in real-time. The proposed algorithm had an accuracy (±precision) of -0.03 ± 0.32 mm, -0.01 ± 0.13 mm and 0.03 ± 0.52 mm for translations in the left-right (LR), superior-inferior (SI) and anterior-posterior (AP) directions respectively; and, 0.07 ± 1.18°, 0.07 ± 1.00° and 0.06 ± 1.32° for rotations around the LR, SI and AP axes respectively on the dataset. The first method to directly estimate real-time 6DoF target motion from segmented marker positions on a 2D imager was devised. The algorithm was evaluated using 81 motion traces from 19 liver patients and was found to have sub-mm and sub-degree accuracy

Real-time intrafraction motion monitoring in external beam radiotherapy

© 2019 Institute of Physics and Engineering in Medicine. Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT

An augmented correlation framework for the estimation of tumour translational and rotational motion during external beam radiotherapy treatments using intermittent monoscopic x-ray imaging and an external respiratory signal

© 2018 Institute of Physics and Engineering in Medicine. Increasing evidence shows that intrafraction tumour motion monitoring must include both six degrees of freedom (6DoF): 3D translations and 3D rotations. Existing real-time algorithms for 6DoF target motion estimation require continuous intrafraction fluoroscopic imaging at high frequency, thereby exposing patients to additional high imaging dose. This paper presents the first method capable of 6DoF motion monitoring using intermittent 2D kV imaging and a continuous external respiratory signal. Our approach is to optimise a state-augmented linear correlation model between an external signal and internal 6DoF motion. In standard treatments, the model can be built using information obtained during pre-treatment cone beam CT (CBCT). Real-time 6DoF tumor motion can then be estimated using just the external signal. Intermittent intrafraction kV images are used to update the model parameters, accounting for changes in correlation and baseline shifts. The method was evaluated in silico using data from 6 lung SABR patients, with the internal tumour motion recorded with electromagnetic beacons and the external signal from a bellows belt. Projection images from CBCT (10 Hz) and intermittent kV images were simulated by projecting the 3D Calypso beacon positions onto an imager. IMRT and VMAT treatments were simulated with increasing imaging update intervals: 0.1 s, 1 s, 3 s, 10 s and 30 s. For all the tested clinical scenarios, translational motion estimates with our method had sub-mm accuracy (mean) and precision (standard deviation) while rotational motion estimates were accurate to < and precise to . Motion estimation errors increased as the imaging update interval increased. With the largest imaging update interval (30 s), the errors were mm, mm and mm for translation in the left-right, superior-inferior and anterior-posterior directions, respectively, and , and for rotation around the aforementioned axes for both VMAT and IMRT treatments. In conclusion, we developed and evaluated a novel method for highly accurate real-time 6DoF motion monitoring on a standard linear accelerator without requiring continuous kV imaging. The proposed method achieved sub-mm and sub-degree accuracy on a lung cancer patient dataset

A comparison of gantry-mounted x-ray-based real-time target tracking methods.

PURPOSE: Most modern radiotherapy machines are built with a 2D kV imaging system. Combining this imaging system with a 2D-3D inference method would allow for a ready-made option for real-time 3D tumor tracking. This work investigates and compares the accuracy of four existing 2D-3D inference methods using both motion traces inferred from external surrogates and measured internally from implanted beacons. METHOD: Tumor motion data from 160 fractions (46 thoracic/abdominal patients) of Synchrony traces (inferred traces), and 28 fractions (7 lung patients) of Calypso traces (internal traces) from the LIGHT SABR trial (NCT02514512) were used in this study. The motion traces were used as the ground truth. The ground truth trajectories were used in silico to generate 2D positions projected on the kV detector. These 2D traces were then passed to the 2D-3D inference methods: interdimensional correlation, Gaussian probability density function (PDF), arbitrary-shape PDF, and the Kalman filter. The inferred 3D positions were compared with the ground truth to determine tracking errors. The relationships between tracking error and motion magnitude, interdimensional correlation, and breathing periodicity index (BPI) were also investigated. RESULTS: Larger tracking errors were observed from the Calypso traces, with RMS and 95th percentile 3D errors of 0.84-1.25 mm and 1.72-2.64 mm, compared to 0.45-0.68 mm and 0.74-1.13 mm from the Synchrony traces. The Gaussian PDF method was found to be the most accurate, followed by the Kalman filter, the interdimensional correlation method, and the arbitrary-shape PDF method. Tracking error was found to strongly and positively correlate with motion magnitude for both the Synchrony and Calypso traces and for all four methods. Interdimensional correlation and BPI were found to negatively correlate with tracking error only for the Synchrony traces. The Synchrony traces exhibited higher interdimensional correlation than the Calypso traces especially in the anterior-posterior direction. CONCLUSION: Inferred traces often exhibit higher interdimensional correlation, which are not true representation of thoracic/abdominal motion and may underestimate kV-based tracking errors. The use of internal traces acquired from systems such as Calypso is advised for future kV-based tracking studies. The Gaussian PDF method is the most accurate 2D-3D inference method for tracking thoracic/abdominal targets. Motion magnitude has significant impact on 2D-3D inference error, and should be considered when estimating kV-based tracking error

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

PURPOSE: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. METHODS: Because the superior-inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left-right and anterior-posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors' simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. RESULTS: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. CONCLUSIONS: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.NHMRC, National Research Foundation of Kore

Development and clinical evaluation of automatic fiducial detection for tumor tracking in cine megavoltage images during volumetric modulated arc therapy

PURPOSE: Real-time tracking of implanted fiducials in cine megavoltage (MV) imaging during volumetric modulated arc therapy (VMAT) delivery is complicated due to the inherent low contrast of MV images and potential blockage of dynamic leaves configurations. The purpose of this work is to develop a clinically practical autodetection algorithm for motion management during VMAT. METHODS: The expected field-specific segments and the planned fiducial position from the Eclipse (Varian Medical Systems, Palo Alto, CA) treatment planning system were projected onto the MV images. The fiducials were enhanced by applying a Laplacian of Gaussian filter in the spatial domain for each image, with a blob-shaped object as the impulse response. The search of implanted fiducials was then performed on a region of interest centered on the projection of the fiducial when it was within an open field including the case when it was close to the field edge or partially occluded by the leaves. A universal template formula was proposed for template matching and normalized cross correlation was employed for its simplicity and computational efficiency. The search region for every image was adaptively updated through a prediction model that employed the 3D position of the fiducial estimated from the localized positions in previous images. This prediction model allowed the actual fiducial position to be tracked dynamically and was used to initialize the search region. The artifacts caused by electronic interference during the acquisition were effectively removed. A score map was computed by combining both morphological information and image intensity. The pixel location with the highest score was selected as the detected fiducial position. The sets of cine MV images taken during treatment were analyzed with in-house developed software written in MATLAB (The Mathworks, Inc., Natick, MA). Five prostate patients were analyzed to assess the algorithm performance by measuring their positioning accuracy during treatment. RESULTS: The algorithm was able to accurately localize the fiducial position on MV images with success rates of more than 90% per case. The percentage of images in which each fiducial was localized in the studied cases varied between 23% and 65%, with at least one fiducial having been localized between 40% and 95% of the images. This depended mainly on the modulation of the plan and fiducial blockage. The prostate movement in the presented cases varied between 0.8 and 3.5 mm (mean values). The maximum displacement detected among all patients was of 5.7 mm. CONCLUSIONS: An algorithm for automatic detection of fiducial markers in cine MV images has been developed and tested with five clinical cases. Despite the challenges posed by complex beam aperture shapes, fiducial localization close to the field edge, partial occlusion of fiducials, fast leaf and gantry movement, and inherently low MV image quality, good localization results were achieved in patient images. This work provides a technique for enabling real-time accurate fiducial detection and tumor tracking during VMAT treatments without the use of extra imaging dose

IGRT and motion management during lung SBRT delivery.

Patient motion can cause misalignment of the tumour and toxicities to the healthy lung tissue during lung stereotactic body radiation therapy (SBRT). Any deviations from the reference setup can miss the target and have acute toxic effects on the patient with consequences onto its quality of life and survival outcomes. Correction for motion, either immediately prior to treatment or intra-treatment, can be realized with image-guided radiation therapy (IGRT) and motion management devices. The use of these techniques has demonstrated the feasibility of integrating complex technology with clinical linear accelerator to provide a higher standard of care for the patients and increase their quality of life

Adaptive Radiation Therapy (ART) Strategies and Technical Considerations: A State of the ART Review From NRG Oncology

The integration of adaptive radiation therapy (ART), or modifying the treatment plan during the treatment course, is becoming more widely available in clinical practice. ART offers strong potential for minimizing treatment-related toxicity while escalating or de-escalating target doses based on the dose to organs at risk. Yet, ART workflows add complexity into the radiation therapy planning and delivery process that may introduce additional uncertainties. This work sought to review presently available ART workflows and technological considerations such as image quality, deformable image registration, and dose accumulation. Quality assurance considerations for ART components and minimum recommendations are described. Personnel and workflow efficiency recommendations are provided, as is a summary of currently available clinical evidence supporting the implementation of ART. Finally, to guide future clinical trial protocols, an example ART physician directive and a physics template following standard NRG Oncology protocol is provided