Intra-fractional lung tumor motion monitoring using arbitrary gantry angles during radiotherapy treatment.

Intensity-based 2D/3D registration using kilo-voltage (kV) and mega-voltage (MV) on-board imaging is a promising approach for real-time tumor motion tracking. So far, the performance of the kV images as well as kV-MV image pairs for 2D/3D registration using only one gantry angle (in anterior-posterior (AP) direction) has been investigated on patient data. In stereotactic body radiation therapy (SBRT), however, various gantry angles are typically used. This study attempts to answer the question of whether automatic 2D/3D registration is possible using kV images as well as kV-MV image pairs for gantry angles other than the AP direction. We also investigated the effect of additional portal MV images paired with kV images to improve 2D/3D registration in extracting cranio-caudal (CC) and AP displacement at arbitrary gantry angles and different fractions. The kV and MV image sequences as well as 3D volume data from five patients suffering from non-small cell lung cancer undergoing SBRT were used. Diaphragm motion served as the reference signal. The CC and AP displacements resulting from the registration results were compared with the corresponding reference motion signal. Pearson correlation coefficients (R value) was used to calculate the similarity measure between reference signal and the extracted displacements resulting from the registration. Signals we found that using 2D/3D registration tumor motion in 5 degrees of freedom (DOF) with kV images and in 6 degrees of freedom with kV-MV image pairs can be extracted for most gantry angles in all patients. Furthermore, our results have shown that the use of kV-MV image pairs increases the overall chance of tumor visibility and therefore leads to more successful extraction of CC as well as AP displacements for almost all gantry angles in all patients. We observed an improvement in registration of at least 0.29% more gantry angle for all patients when we used kV-MV images compared to kV images alone. In addition, an improvement in the R-value was observed in up to 16 fractions in various patients.

Real-time 2D/3D registration using kV-MV image pairs for tumor motion tracking in image guided radiotherapy.

Intra-fractional respiratory motion during radiotherapy leads to a larger planning target volume (PTV). Real-time tumor motion tracking by two-dimensional (2D)/3D registration using on-board kilo-voltage (kV) imaging can allow for a reduction of the PTV though motion along the imaging beam axis cannot be resolved using only one projection image. We present a retrospective patient study investigating the impact of paired portal mega-voltage (MV) and kV images on registration accuracy. Material and methods. We used data from 10 patients suffering from non-small cell lung cancer (NSCLC) undergoing stereotactic body radiation therapy (SBRT) lung treatment. For each patient we acquired a planning computed tomography (CT) and sequences of kV and MV images during treatment. We compared the accuracy of motion tracking in six degrees-of-freedom (DOF) using the anterior-posterior (AP) kV sequence or the sequence of kV-MV image pairs. Results. Motion along cranial-caudal direction could accurately be extracted when using only the kV sequence but in AP direction we obtained large errors. When using kV-MV pairs, the average error was reduced from 2.9 mm to 1.5 mm and the motion along AP was successfully extracted. Mean registration time was 188 ms. Conclusion. Our evaluation shows that using kV-MV image pairs leads to improved motion extraction in six DOF and is suitable for real-time tumor motion tracking with a conventional LINAC.

Endoclamp balloon visualization and automatic placement system.

OBJECTIVES: Aortic occlusion is one of the most important open discussions in minimally invasive cardiac surgery. Different techniques can be employed, and all have benefits and drawbacks. The objective of our work is to improve the safety of internal aortic occlusion with the Port Access technique, which employs an endoclamp balloon catheter. We propose a combined information and positioning system based on augmented reality technology and robotics in which the position of the balloon can be seen at all times and can be automatically controlled by a robotic actuator. METHODS: The system was designed by a multidisciplinary team of engineers, medical doctors, and human factor specialists in a human-centered design approach. We measure the balloon position in real time with a magnetic tracking system. This position is superimposed on a 3-dimensional scan of the patient's thorax, with the balloon in the artery shown at all times. The position measurement is also used to control the robotic catheter inserter that places and maintains the balloon position at a specified target. The system was evaluated in 2 user studies that compared it with other visual aids. RESULTS: The user tests have shown that the system effectively supports the surgeon in the placement task, with an increase in placement accuracy and a reduction in time compared with the current visualization technique. The users also rated the system as supporting them well. CONCLUSIONS: The clinical feasibility of the system was proved. The system provides better visualization and position control and can effectively increase the safety of the procedure. This system has the potential of making Port Access a more attractive technique.

Investigating conditional GAN performance with different generator architectures, an ensemble model, and different MR scanners for MR-sCT conversion.

Recent developments in magnetic resonance (MR) to synthetic computed tomography (sCT) conversion have shown that treatment planning is possible without an initial planning CT. Promising conversion results have been demonstrated recently using conditional generative adversarial networks (cGANs). However, the performance is generally only tested on images from one MR scanner, which neglects the potential of neural networks to find general high-level abstract features. In this study, we explored the generalizability of the generator models, trained on a single field strength scanner, to data acquired with higher field strengths. T2-weighted 0.35T MRIs and CTs from 51 patients treated for prostate (40) and cervical cancer (11) were included. 25 of them were used to train four different generators (SE-ResNet, DenseNet, U-Net, and Embedded Net). Further, an ensemble model was created from the four network outputs. The models were validated on 16 patients from a 0.35T MR scanner. Further, the trained models were tested on the Gold Atlas dataset, containing T2-weighted MR scans of different field strengths; 1.5T(7) and 3T(12), and 10 patients from the 0.35T scanner. The sCTs were dosimetrically compared using clinical VMAT plans for all test patients. For the same scanner (0.35T), the results from the different models were comparable on the test set, with only minor differences in the mean absolute error (MAE) (35-51HU body). Similar results were obtained for conversions of 3T GE Signa and the 3T GE Discovery images (40-62HU MAE) for three of the models. However, larger differences were observed for the 1.5T images (48-65HU MAE). The overall best model was found to be the ensemble model. All dose differences were below 1%. This study shows that it is possible to generalize models trained on images of one scanner to other scanners and different field strengths. The best metric results were achieved by the combination of all networks.

The first prospective implementation of markerless lung target tracking in an experimental quality assurance procedure on a standard linear accelerator.

The ability to track tumour motion without implanted markers on a standard linear accelerator (linac) could enable wide access to real-time adaptive radiotherapy for cancer patients. We previously have retrospectively validated a method for 3D markerless target tracking using intra-fractional kilovoltage (kV) projections acquired on a standard linac. This paper presents the first prospective implementation of markerless lung target tracking on a standard linac and its quality assurance (QA) procedure. The workflow and the algorithm developed to track the 3D target position during volumetric modulated arc therapy treatment delivery were optimised. The linac was operated in clinical QA mode, while kV projections were streamed to a dedicated computer using a frame-grabber software. The markerless target tracking accuracy and precision were measured in a lung phantom experiment under the following conditions: static localisation of seven distinct positions, dynamic localisation of five patient-measured motion traces, and dynamic localisation with treatment interruption. The QA guidelines were developed following the AAPM Task Group 147 report with the requirement that the tracking margin components, the margins required to account for tracking errors, did not exceed 5 mm in any direction. The mean tracking error ranged from 0.0 to 0.9 mm (left-right), -0.6 to -0.1 mm (superior-inferior) and -0.7 to 0.1 mm (anterior-posterior) over the three tests. Larger errors were found in cases with large left-right or anterior-posterior and small superior-inferior motion. The tracking margin components did not exceed 5 mm in any direction and ranged from 0.4 to 3.2 mm (left-right), 0.7 to 1.6 mm (superior-inferior) and 0.8 to 1.5 mm (anterior-posterior). This study presents the first prospective implementation of markerless lung target tracking on a standard linac and provides a QA procedure for its safe clinical implementation, potentially enabling real-time adaptive radiotherapy for a large population of lung cancer patients.

Optimization for customized trajectories in cone beam computed tomography.

PURPOSE: We developed a target-based cone beam computed tomography (CBCT) imaging framework for optimizing an unconstrained three dimensional (3D) source-detector trajectory by incorporating prior image information. Our main aim is to enable a CBCT system to provide topical information about the target using a limited angle noncircular scan orbit with a minimal number of projections. Such a customized trajectory should include enough information to sufficiently reconstruct a particular volume of interest (VOI) under kinematic constraints, which may result from the patient size or additional surgical or radiation therapy-related equipment. METHODS: A patient-specific model from a prior diagnostic computed tomography (CT) volume is used as a digital phantom for CBCT trajectory simulations. Selection of the best projection views is accomplished through maximizing an objective function fed by the imaging quality provided by different x-ray positions on the digital phantom data. The final optimized trajectory includes a limited angular range and a minimal number of projections which can be applied to a C-arm device capable of general source-detector positioning. The performance of the proposed framework is investigated in experiments involving an in-house-built box phantom including spherical targets as well as an Alderson-Rando head phantom. In order to quantify the image quality of the reconstructed image, we use the average full-width-half-maximum (FWHMavg ) for the spherical target and feature similarity index (FSIM), universal quality index (UQI), and contrast-to-noise ratio (CNR) for an anatomical target. RESULTS: Our experiments based on both the box and head phantom showed that optimized trajectories could achieve a comparable image quality in the VOI with respect to the standard C-arm circular CBCT while using approximately one quarter of projections. We achieved a relative deviation <7% for FWHMavg between the reconstructed images from the optimized trajectories and the standard C-arm CBCT for all spherical targets. Furthermore, for the anatomical target, the relative deviation of FSIM, UQI, and CNR between the reconstructed image related to the proposed trajectory and the standard C-arm circular CBCT was found to be 5.06%, 6.89%, and 8.64%, respectively. We also compared our proposed trajectories to circular trajectories with equivalent angular sampling as the optimized trajectories. Our results show that optimized trajectories can outperform simple partial circular trajectories in the VOI in term of image quality. Typically, an angular range between 116° and 152° was used for the optimized trajectories. CONCLUSION: We demonstrated that applying limited angle noncircular trajectories with optimized orientations in 3D space can provide a suitable image quality for particular image targets and has a potential for limited angle and low-dose CBCT-based interventions under strong spatial constraints.

Clinical implementations of 4D pencil beam scanned particle therapy: Report on the 4D treatment planning workshop 2016 and 2017.

In 2016 and 2017, the 8th and 9th 4D treatment planning workshop took place in Groningen (the Netherlands) and Vienna (Austria), respectively. This annual workshop brings together international experts to discuss research, advances in clinical implementation as well as problems and challenges in 4D treatment planning, mainly in spot scanned proton therapy. In the last two years several aspects like treatment planning, beam delivery, Monte Carlo simulations, motion modeling and monitoring, QA phantoms as well as 4D imaging were thoroughly discussed. This report provides an overview of discussed topics, recent findings and literature review from the last two years. Its main focus is to highlight translation of 4D research into clinical practice and to discuss remaining challenges and pitfalls that still need to be addressed and to be overcome.

Advanced Radiation DOSimetry phantom (ARDOS): a versatile breathing phantom for 4D radiation therapy and medical imaging.

A novel breathing phantom was designed for being used in conventional and ion-beam radiotherapy as well as for medical imaging. Accurate dose delivery and patient safety are aimed to be verified for four-dimensional (4D) treatment techniques compensating for breathing-induced tumor motion. The phantom includes anthropomorphic components representing an average human thorax. It consists of real tissue equivalent materials to fulfill the requirements for dosimetric experiments and imaging purposes. The different parts of the torso (lungs, chest wall, and ribs) and the tumor can move independently. Simple regular movements, as well as more advanced patient-specific breathing cycles are feasible while a reproducible setup can be guaranteed. The phantom provides the flexibility to use different types of dosimetric devices and was designed in a way that it is robust, transportable and easy to handle. Tolerance levels and the reliability of the phantom setup were determined in combination with tests on motion accuracy and reproducibility by using infrared optical tracking technology. Different imaging was performed including positron emission tomography imaging, 4D computed tomography as well as real-time in-room imaging. The initial dosimetric benchmarking studies were performed in a photon beam where dose parameters are predictable and the dosimetric procedures well established.

Clinical Assessment of 2D/3D Registration Accuracy in 4 Major Anatomic Sites Using On-Board 2D Kilovoltage Images for 6D Patient Setup.

To provide a comprehensive assessment of patient setup accuracy in 6 degrees of freedom (DOFs) using 2-dimensional/3-dimensional (2D/3D) image registration with on-board 2-dimensional kilovoltage (OB-2 DkV) radiographic images, we evaluated cranial, head and neck (HN), and thoracic and abdominal sites under clinical conditions. A fast 2D/3D image registration method using graphics processing unit GPU was modified for registration between OB-2 DkV and 3D simulation computed tomography (simCT) images, with 3D/3D registration as the gold standard for 6 DOF alignment. In 2D/3D registration, body roll rotation was obtained solely by matching orthogonal OB-2 DkV images with a series of digitally reconstructed radiographs (DRRs) from simCT with a small rotational increment along the gantry rotation axis. The window/level adjustments for optimal visualization of the bone in OB-2 DkV and DRRs were performed prior to registration. Ideal patient alignment at the isocenter was calculated and used as an initial registration position. In 3D/3D registration, cone-beam CT (CBCT) was aligned to simCT on bony structures using a bone density filter in 6DOF. Included in this retrospective study were 37 patients treated in 55 fractions with frameless stereotactic radiosurgery or stereotactic body radiotherapy for cranial and paraspinal cancer. A cranial phantom was used to serve as a control. In all cases, CBCT images were acquired for patient setup with subsequent OB-2 DkV verification. It was found that the accuracy of the 2D/3D registration was 0.0 ± 0.5 mm and 0.1° ± 0.4° in phantom. In patient, it is site dependent due to deformation of the anatomy: 0.2 ± 1.6 mm and -0.4° ± 1.2° on average for each dimension for the cranial site, 0.7 ± 1.6 mm and 0.3° ± 1.3° for HN, 0.7 ± 2.0 mm and -0.7° ± 1.1° for the thorax, and 1.1 ± 2.6 mm and -0.5° ± 1.9° for the abdomen. Anatomical deformation and presence of soft tissue in 2D/3D registration affect the consistency with 3D/3D registration in 6 DOF: the discrepancy increases in superior to inferior direction.

A quantitative comparison of the performance of three deformable registration algorithms in radiotherapy.

We present an evaluation of various non-rigid registration algorithms for the purpose of compensating interfractional motion of the target volume and organs at risk areas when acquiring CBCT image data prior to irradiation. Three different deformable registration (DR) methods were used: the Demons algorithm implemented in the iPlan Software (BrainLAB AG, Feldkirchen, Germany) and two custom-developed piecewise methods using either a Normalized Correlation or a Mutual Information metric (featureletNC and featureletMI). These methods were tested on data acquired using a novel purpose-built phantom for deformable registration and clinical CT/CBCT data of prostate and lung cancer patients. The Dice similarity coefficient (DSC) between manually drawn contours and the contours generated by a derived deformation field of the structures in question was compared to the result obtained with rigid registration (RR). For the phantom, the piecewise methods were slightly superior, the featureletNC for the intramodality and the featureletMI for the intermodality registrations. For the prostate cases in less than 50% of the images studied the DSC was improved over RR. Deformable registration methods improved the outcome over a rigid registration for lung cases and in the phantom study, but not in a significant way for the prostate study. A significantly superior deformation method could not be identified.

High-performance GPU-based rendering for real-time, rigid 2D/3D-image registration and motion prediction in radiation oncology.

A common problem in image-guided radiation therapy (IGRT) of lung cancer as well as other malignant diseases is the compensation of periodic and aperiodic motion during dose delivery. Modern systems for image-guided radiation oncology allow for the acquisition of cone-beam computed tomography data in the treatment room as well as the acquisition of planar radiographs during the treatment. A mid-term research goal is the compensation of tumor target volume motion by 2D/3D Registration. In 2D/3D registration, spatial information on organ location is derived by an iterative comparison of perspective volume renderings, so-called digitally rendered radiographs (DRR) from computed tomography volume data, and planar reference x-rays. Currently, this rendering process is very time consuming, and real-time registration, which should at least provide data on organ position in less than a second, has not come into existence. We present two GPU-based rendering algorithms which generate a DRR of 512×512 pixels size from a CT dataset of 53 MB size at a pace of almost 100 Hz. This rendering rate is feasible by applying a number of algorithmic simplifications which range from alternative volume-driven rendering approaches - namely so-called wobbled splatting - to sub-sampling of the DRR-image by means of specialized raycasting techniques. Furthermore, general purpose graphics processing unit (GPGPU) programming paradigms were consequently utilized. Rendering quality and performance as well as the influence on the quality and performance of the overall registration process were measured and analyzed in detail. The results show that both methods are competitive and pave the way for fast motion compensation by rigid and possibly even non-rigid 2D/3D registration and, beyond that, adaptive filtering of motion models in IGRT.

Monitoring tumor motion by real time 2D/3D registration during radiotherapy.

BACKGROUND AND PURPOSE: In this paper, we investigate the possibility to use X-ray based real time 2D/3D registration for non-invasive tumor motion monitoring during radiotherapy. MATERIALS AND METHODS: The 2D/3D registration scheme is implemented using general purpose computation on graphics hardware (GPGPU) programming techniques and several algorithmic refinements in the registration process. Validation is conducted off-line using a phantom and five clinical patient data sets. The registration is performed on a region of interest (ROI) centered around the planned target volume (PTV). RESULTS: The phantom motion is measured with an rms error of 2.56 mm. For the patient data sets, a sinusoidal movement that clearly correlates to the breathing cycle is shown. Videos show a good match between X-ray and digitally reconstructed radiographs (DRR) displacement. Mean registration time is 0.5 s. CONCLUSIONS: We have demonstrated that real-time organ motion monitoring using image based markerless registration is feasible.