Advances in 4D radiation therapy for managing respiration: part II - 4D treatment planning

Comparison of translation algorithms in determining maximum allowable CTV shifts for Real-Time Gated Proton Therapy (RGPT) robustness evaluation in prostate cancers

INTRODUCTION

Real-Time Gated Proton Therapy (RGPT) is an active motion management technique that utilizes treatment gating and tumor tracking via fiducial markers. When performing RGPT treatment for prostate cancer, it is essential to account for the CTV displacement relative to the body in the clinical workflow. The workflow at the National Cancer Centre Singapore (NCCS) includes bone matching via CT-CBCT images, followed by fiducial matching via pulsed fluoroscopy (soft tissue matching), and finally, a robustness evaluation procedure to determine if the difference is within an allowable tolerance. In this study, we compare two CTV translation methods for robustness evaluation: (1) an in-house translation algorithm and (2) the RayStation "simulate organ motion" Deformable image registration (DIR) algorithm.

METHODS

Nine RGPT prostate patient plans with CTV volumes ranging from 17.1 to 96.72 cm2 were included in this study. An in-house translation algorithm and "simulate organ motion" DIR RayStation algorithm were used to generate CTV shifts along R-L, I-S, and P-A axes between ± $ \pm $ 10 mm at 2 mm steps. At each step, dose metrics, which include CTV Dmax, CTV D95%, and CTV D98%, were extracted and used as comparative metrics for CTV target coverage and hot spot evaluation.

RESULTS

Across all axes, there were no statistically significant differences between the two algorithms for all three dose metrics: CTV Dmax (P = 0.92, P = 0.91, and P = 0.47), CTV D95% (P = 0.97, P = 0.22, and P = 0.33), and CTV D98% (P = 0.85, P = 0.33, and P = 0.36). Further, the in-house translation algorithm evaluation time was less than 10 s, two orders of magnitude faster than the DIR algorithm.

CONCLUSION

Our results demonstrate that the simpler in-house algorithm performs equivalently to the realistic DIR algorithm when simulating CTV motion in prostate cancers. Furthermore, the in-house algorithm completes the robustness evaluation two orders of magnitude faster than the DIR algorithm. This significant reduction in evaluation time is crucial especially when preparatory time efficiency is of paramount importance in a busy clinic.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

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.

CycN-Net: A Convolutional Neural Network Specialized for 4D CBCT Images Refinement

Four-dimensional cone-beam computed tomography (4D CBCT) has been developed to provide a sequence of phase-resolved reconstructions in image-guided radiation therapy. However, 4D CBCT images are degraded by severe streaking artifacts and noise because the phase-resolved image is an extremely sparse-view CT procedure wherein a few under-sampled projections are used for the reconstruction of each phase. Aiming at improving the overall quality of 4D CBCT images, we proposed two CNN models, named N-Net and CycN-Net, respectively, by fully excavating the inherent property of 4D CBCT. To be specific, the proposed N-Net incorporates the prior image reconstructed from entire projection data based on U-Net to boost the image quality for each phase-resolved image. Based on N-Net, a temporal correlation among the phase-resolved images is also considered by the proposed CycN-Net. Extensive experiments on both XCAT simulation data and real patient 4D CBCT datasets were carried out to verify the feasibility of the proposed CNNs. Both networks can effectively suppress streaking artifacts and noise while restoring the distinct features simultaneously, compared with the existing CNN models and two state-of-the-art iterative algorithms. Moreover, the proposed method is robust in handling complicated tasks of various patient datasets and imaging devices, which implies its excellent generalization ability.

Respiratory deformation registration in 4D-CT/cone beam CT using deep learning

BACKGROUND

To investigate the feasibility of using a supervised convolutional neural network (CNN) to register phase-to-phase deformable vector field of lung 4D-CT/4D-cone beam CT for 4D dose accumulation, contour propagation, motion modeling, or target verification.

METHODS

We built a CNN-based deep learning method to register the deformation field directly between phases of patients' 4D-CT or 4D-cone beam CT. The input consists of patch pairs of two phases, while the output is the corresponding deformation field that registers the patch pairs. The centers of the patch pairs were uniformly sampled across the lung, and the size of the patches was chosen to cover the range of the respiratory motion. The network was trained to generate deformation field that matches with the reference deformation field generated by VelocityAI (Varian). The network is structured with four convolutional layers, two average pooling layers, and two fully connected layers. Half mean squared error is applied to guide the study as loss function. Nine patients with eleven sets of 4D-CT/cone beam CT image volumes were used for training and testing. The performance of the network was validated with intra-patient and inter-patient setups.

RESULTS

Registered images were generated with Velocity deformation field and the CNN deformation field, respectively. Main anatomic features such as the main vessels and the diaphragm matched well between two deformed images. In the diaphragm region, the coefficients of cross-correlation, root mean squared error, and structural similarity index measure (SSIM) between deformed images registered by CNN and VelocityAI was calculated. The cross-correlation was above 0.9 for all the intra-patient cases.

CONCLUSIONS

Patch-based deep learning methods achieved comparable deformable registration accuracy as VelocityAI. Compared to VelocityAI, the deep learning method is fully automatic and faster without user dependency, which makes it more preferable in clinical applications.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors

PURPOSE

To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario.

MATERIALS AND METHODS

The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed.

RESULTS

For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D_95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D_95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D_95% and gastro-intestinal OARs D_2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs.

CONCLUSIONS

The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

Reliability of ITV approach to varying treatment fraction time: a retrospective analysis based on 2D cine MR images

BACKGROUND

Internal Target Volume (ITV) is one of the most common strategies to passively manage tumour motion in Radiotherapy (RT). The reliability of this approach is based on the assumption that the tumour motion estimated during pre-treatment 4D Computed Tomography (CT) acquisition is representative of the motion during the whole RT treatment. With the introduction of Magnetic Resonance-guided RT (MRgRT), it has become possible to monitor tumour motion during the treatment and verify this assumption. Aim of this study was to investigate the reliability of the ITV approach with respect to the treatment fraction time (TFT) in abdominal and thoracic lesions.

METHODS

A total of 12 thoracic and 15 abdominal lesions was analysed. Before treatment, a 10-phase 4DCT was acquired and ITV margins were estimated considering the envelope of the lesion contoured on the different 4DCT phases. All patients underwent MRgRT treatment in free-breathing, monitoring the tumour position on a sagittal plane with 4 frames per second (sec). ITV margins were projected on the tumour trajectory and the percentage of treatment time in which the tumour was inside the ITV (%TT) was measured to varying of TFT. The ITV approach was considered moderately reliable when %TT ≥ 90% and strongly reliable when %TT ≥ 95%. Additional ITV margins required to achieve %TT ≥ 95% were also calculated.

RESULTS

In the analysed cohort of patients, ITV strategy can be considered strongly reliable only for lung lesions with TFT ≤ 7 min (min). The ITV strategy can be considered only moderately reliable for abdominal lesions, and additional margins are required to obtain %TT ≥ 95%. Considering a TFT ≤ 4 min, additional margins of 2 mm in cranio-caudal (CC) and 1 mm in antero-posterior (AP) are suggested for pancreatic lesions, 3 mm in CC and 2 mm in AP for renal and liver ones.

CONCLUSIONS

On the basis of the analysed cases, the ITV approach appears to be reliable in the thorax, while it results more challenging in the abdomen, due to the higher uncertainty in ITV definition and to the observed larger intra and inter-fraction motion variability. The addition of extra margins based on the TFT may represent a valid tool to compensate such limitations.

Online daily adaptive proton therapy

It is recognized that the use of a single plan calculated on an image acquired some time before the treatment is generally insufficient to accurately represent the daily dose to the target and to the organs at risk. This is particularly true for protons, due to the physical finite range. Although this characteristic enables the generation of steep dose gradients, which is essential for highly conformal radiotherapy, it also tightens the dependency of the delivered dose to the range accuracy. In particular, the use of an outdated patient anatomy is one of the most significant sources of range inaccuracy, thus affecting the quality of the planned dose distribution. A plan should be ideally adapted as soon as anatomical variations occur, ideally online. In this review, we describe in detail the different steps of the adaptive workflow and discuss the challenges and corresponding state-of-the art developments in particular for an online adaptive strategy.

The noise navigator: a surrogate for respiratory-correlated 4D-MRI for motion characterization in radiotherapy

Respiratory-correlated 4D-MRI can characterize respiratory-induced motion of tumors and organs-at-risk for radiotherapy treatment planning and is a necessity for image guidance of moving tumors treated on an MRI-linac. Essential for 4D-MRI generation is a robust respiratory surrogate signal. We investigated the feasibility of the noise navigator as respiratory surrogate signal for 4D-MRI generation. The noise navigator is based on the respiratory-induced modulation of the thermal noise variance measured by the receive coils during MR acquisition and thus is inherently present and synchronized with MRI data acquisition. Additionally, the noise navigator can be combined with any rectilinear readout strategy (e.g. radial and cartesian) and is independent of MR image contrast and imaging orientation. For radiotherapy applications, the noise navigator provides a robust respiratory signal for patients scanned with an elevated coil setup. This is particularly attractive for widely used cartesian sequences where currently a non-interfering self-navigation means is lacking for MRI-based simulation and MRI-guided radiotherapy. The feasibility of 4D-MRI generation with the noise navigator as respiratory surrogate signal was demonstrated for both cartesian and radial readout strategies in radiotherapy setup on four healthy volunteers and two radiotherapy patients on a dedicated 1.5 T MRI scanner and two radiotherapy patients on a 1.5 T MRI-linac system. Moreover, the respiratory-correlated 4D-MR images showed liver motion comparable to a reference 2D cine MRI series for the volunteers. For 2D cartesian cine MRI acquisitions, both the noise navigator and respiratory bellows were benchmarked against an image navigator. Respiratory phase detection based on the noise navigator agreed 1.4 times better with the image navigator than the respiratory bellows did. For a 3D Stack-of-Stars acquisitions, the noise navigator was compared to radial self-navigation and a 1.7 times higher respiratory phase detection agreement was observed than for the respiratory bellows compared to radial self-navigation.

Accounting for respiratory motion in small serial structures during radiotherapy planning: proof of concept in virtual bronchoscopy-guided lung functional avoidance radiotherapy

Respiratory motion management techniques in radiotherapy (RT) planning are primarily focused on maintaining tumor target coverage. An inadequately addressed need is accounting for motion in dosimetric estimations in smaller serial structures. Accurate dose estimations in such structures are more sensitive to motion because respiration can cause them to move completely in or out of a high dose-gradient field. In this work, we study three motion management strategies (m1-m3) to find an accurate method to estimate the dosimetry in airways. To validate these methods, we generated a 'ground truth' digital breathing model based on a 4DCT scan from a lung stereotactic ablative radiotherapy (SAbR) patient. We simulated 225 breathing cycles with  ±10% perturbations in amplitude, respiratory period, and time per respiratory phase. A high-resolution breath-hold CT (BHCT) was also acquired and used with a research virtual bronchoscopy software to autosegment 239 airways. Contours for planning target volume (PTV) and organs at risk (OARs) were defined on the maximum intensity projection of the 4DCT (CT_MIP) and transferred to the average of the 10 4DCT phases (CT_AVG). To design the motion management methods, the RT plan was recreated using different images and structure definitions. Methods m1 and m2 recreated the plan using the CT_AVG image. In method m1, airways were deformed to the CT_AVG. In m2, airways were deformed to each of the 4DCT phases, and union structures were transferred onto the CT_AVG. In m3, the RT plan was recreated on each of the 10 phases, and the dose distribution from each phase was deformed to the BHCT and summed. Dose errors (mean [min, max]) in airways were: m1: 21% (0.001%, 93%); m2: 45% (0.1%, 179%); and m3: 4% (0.006%, 14%). Our work suggests that accurate dose estimation in moving small serial structures requires customized motion management techniques (like m3 in this work) rather than current clinical and investigational approaches.

Comparison of planned dose on different CT image sets to four-dimensional Monte Carlo dose recalculation using the patient's actual breathing trace for lung stereotactic body radiation therapy

PURPOSE

The need for four-dimensional (4D) treatment planning becomes indispensable when it comes to radiation therapy for moving tumors in the thoracic and abdominal regions. The primary purpose of this study is to combine the actual breathing trace during each individual treatment fraction with the Linac's log file information and Monte Carlo 4D dose calculations. We investigated this workflow on multiple computed tomography (CT) datasets in a clinical environment for stereotactic body radiation therapy (SBRT) treatment planning.

METHODS

We have developed a workflow, which allows us to recalculate absorbed dose to a 4DCT dataset using Monte Carlo calculation methods and accumulate all 4D doses in order to compare them to the planned dose using the Linac's log file, a 4DCT dataset, and the patient's actual breathing curve for each individual fraction. For five lung patients, three-dimensional-conformal radiation therapy (3D-CRT) and volumetric modulated arc treatment (VMAT) treatment plans were generated on four different CT image datasets: a native free-breathing 3DCT, an average intensity projection (AIP) and a maximum intensity projection (MIP) CT both obtained from a 4DCT, and a 3DCT with density overrides based on the 3DCT (DO). The Monte Carlo 4D dose has been calculated on each 4DCT phase using the Linac's log file and the patient's breathing trace as a surrogate for tumor motion and dose was accumulated to the gross tumor volume (GTV) at the 50% breathing phase (end of exhale) using deformable image registration.

RESULTS

Δ D 98 % and Δ D 2 % between 4D dose and planned dose differed largely for 3DCT-based planning and also for DO in three patients. Least dose differences between planned and recalculated dose have been found for AIP and MIP treatment planning which both tend to be superior to DO, but the results indicate a dependency on the breathing variability, tumor motion, and size. An interplay effect has not been observed in the small patient cohort.

CONCLUSIONS

We have developed a workflow which, to our best knowledge, is the first incorporation of the patient breathing trace over the course of all individual treatment fractions with the Linac's log file information and 4D Monte Carlo recalculations of the actual treated dose. Due to the small patient cohort, no clear recommendation on which CT can be used for SBRT treatment planning can be given, but the developed workflow, after adaption for clinical use, could be used to enhance a priori 4D Monte Carlo treatment planning in the future and help with the decision on which CT dataset treatment planning should be carried out.

Deformable Registration for Dose Accumulation

As deformable image registration makes its way into the clinical routine, the summation of doses from fractionated treatment regimens to evaluate cumulative doses to targets and healthy tissues is also becoming a frequently utilized tool in the context of image-guided adaptive radiotherapy. Accounting for daily geometric changes using deformable image registration and dose accumulation potentially enables a better understanding of dose-volume-effect relationships, with the goal of translation of this knowledge to personalization of treatment, to further enhance treatment outcomes. Treatment adaptation involving image deformation requires patient-specific quality assurance of the image registration and dose accumulation processes, to ensure that uncertainties in the 3D dose distributions are identified and appreciated from a clinical relevance perspective. While much research has been devoted to identifying and managing the uncertainties associated with deformable image registration and dose accumulation approaches, there are still many unanswered questions. Here, we provide a review of current deformable image registration and dose accumulation techniques, and related clinical application. We also discuss salient issues that need to be deliberated when applying deformable algorithms for dose mapping and accumulation in the context of adaptive radiotherapy and response assessment.

MRI-guidance for motion management in external beam radiotherapy: current status and future challenges

High precision conformal radiotherapy requires sophisticated imaging techniques to aid in target localisation for planning and treatment, particularly when organ motion due to respiration is involved. X-ray based imaging is a well-established standard for radiotherapy treatments. Over the last few years, the ability of magnetic resonance imaging (MRI) to provide radiation-free images with high-resolution and superb soft tissue contrast has highlighted the potential of this imaging modality for radiotherapy treatment planning and motion management. In addition, these advantageous properties motivated several recent developments towards combined MRI radiation therapy treatment units, enabling in-room MRI-guidance and treatment adaptation. The aim of this review is to provide an overview of the state-of-the-art in MRI-based image guidance for organ motion management in external beam radiotherapy. Methodological aspects of MRI for organ motion management are reviewed and their application in treatment planning, in-room guidance and adaptive radiotherapy described. Finally, a roadmap for an optimal use of MRI-guidance is highlighted and future challenges are discussed.

A Monte-Carlo-based and GPU-accelerated 4D-dose calculator for a pencil beam scanning proton therapy system

PURPOSE

The presence of respiratory motion during radiation treatment leads to degradation of the expected dose distribution, both for target coverage and healthy tissue sparing, particularly for techniques like pencil beam scanning proton therapy which have dynamic delivery systems. While tools exist to estimate this degraded four-dimensional (4D) dose, they typically have one or more deficiencies such as not including the particular effects from a dynamic delivery, using analytical dose calculations, and/or using nonphysical dose-accumulation methods. This work presents a clinically useful 4D-dose calculator that addresses each of these shortcomings.

METHODS

To quickly compute the 4D dose, the three main tasks of the calculator were run on graphics processing units (GPUs). These tasks were (a) simulating the delivery of the plan using measured delivery parameters to distribute the plan amongst 4DCT phases characterizing the patient breathing, (b) using an in-house Monte Carlo simulation (MC) dose calculator to determine the dose delivered to each breathing phase, and (c) accumulating the doses from the various breathing phases onto a single phase for evaluation. The accumulation was performed by individually transferring the energy and mass of dose-grid subvoxels, a technique that models the transfer of dose in a more physically realistic manner. The calculator was run on three test cases, with lung, esophagus, and liver targets, respectively, to assess the various uncertainties in the beam delivery simulation as well as to characterize the dose-accumulation technique.

RESULTS

Four-dimensional doses were successfully computed for the three test cases with computation times ranging from 4-6 min on a server with eight NVIDIA Titan X graphics cards; the most time-consuming component was the MC dose engine. The subvoxel-based dose-accumulation technique produced stable 4D-dose distributions at subvoxel scales of 0.5-1.0 mm without impairing the total computation time. The uncertainties in the beam delivery simulation led to moderate variations of the dose-volume histograms for these cases; the variations were reduced by implementing repainting or phase-gating motion mitigation techniques in the calculator.

CONCLUSIONS

A MC-based and GPU-accelerated 4D-dose calculator was developed to estimate the effects of respiratory motion on pencil beam scanning proton therapy treatments. After future validation, the calculator could be used to assess treatment plans and its quick runtime would make it easily usable in a future 4D-robust optimization system.

Impact of deformable image registration on dose accumulation applied electrocardiograph-gated 4DCT in the heart and left ventricular myocardium during esophageal cancer radiotherapy

Background

The deformable image registration (DIR) technique has the potential to realize the dose accumulation during radiotherapy. This study will analyze the feasibility of evaluating dose-volume parameters for the heart and left ventricular myocardium (LVM) by applying DIR.

Methods

The electrocardiograph-gated four-dimensional CT (ECG-gated 4DCT) data of 21 patients were analyzed retrospectively. The heart and LVM were contoured on 20 phases of 4DCT (0%, 5%,…,95%). The heart and LVM in the minimum volume/dice similarity coefficient (DSC) phase (Volume _min/DSC _min) were deformed to the maximum volume/DSC phase (Volume _max/ DSC _max), which used the intensity-based free-form DIR algorithm of MIM software. The dose was deformed according to the deformation vector. The variations in volume, mean dose (D_mean), V₂₀, V₃₀ and V₄₀ for the heart and LVM before and after DIR were compared, and the reference phase was the Volume _max/DSC _max phase.

Results

For the heart, the difference between the pre- and post-registration Volume _min and Volume _max were reduced from 13.87 to 1.72%; the DSC was increased from 0.899 to 0.950 between the pre- and post-registration DSC _min phase relative to the DSC _max phase. The post-registration D_mean, V₂₀, V₃₀ and V₄₀ of the heart were statistically significant compared to those in the Volume _max/DSC _max phase (p < 0.05). For the LVM, the difference between the pre- and post-registration Volume _min and Volume _max were only reduced from 18.77 to 17.38%; the DSC reached only 0.733 in the post-registration DSC _min phase relative to the DSC _max phase. The pre- and post-registration volume, D_mean, V₂₀, V₃₀ and V₄₀ of the LVM were all statistically significant compared to those in the Volume _max/DSC _max phase (p < 0.05).

Conclusions

There was no significant relationship between the variation in dose-volume parameters and the variation in the volume and morphology for the heart; however, the inconsistency of the variation in the volume and morphology for the LVM was a major factor that led to uncertainty in the dose-volume evaluation. In addition, the individualized local deformation registration technology should be applied in dose accumulation for the heart and LVM.

Influence of deformable image registration on 4D dose simulation for extracranial SBRT: A multi-registration framework study

BACKGROUND AND PURPOSE

To evaluate the influence of deformable image registration approaches on correspondence model-based 4D dose simulation in extracranial SBRT by means of open source deformable image registration (DIR) frameworks.

MATERIAL AND METHODS

Established DIR algorithms of six different open source DIR frameworks were considered and registration accuracy evaluated using freely available 4D image data. Furthermore, correspondence models (regression-based correlation of external breathing signal measurements and internal structure motion field) were built and model accuracy evaluated. Finally, the DIR algorithms were applied for motion field estimation in radiotherapy planning 4D CT data of five lung and five liver lesion patients, correspondence model formation, and model-based 4D dose simulation. Deviations between the original, statically planned and the 4D-simulated VMAT dose distributions were analyzed and correlated to DIR accuracy differences.

RESULTS

Registration errors varied among the DIR approaches, with lower DIR accuracy translating into lower correspondence modeling accuracy. Yet, for lung metastases, indices of 4D-simulated dose distributions widely agreed, irrespective of DIR accuracy differences. In contrast, liver metastases 4D dose simulation results strongly vary for the different DIR approaches.

CONCLUSIONS

Especially in treatment areas with low image contrast (e.g. the liver), DIR-based 4D dose simulation results strongly depend on the applied DIR algorithm, drawing resulting dose simulations and indices questionable.

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.

Reduction of breathing irregularity-related motion artifacts in low-pitch spiral 4D CT by optimized projection binning

Background

Respiration-correlated CT (4D CT) is the basis of radiotherapy treatment planning of thoracic and abdominal tumors. Current clinical 4D CT images suffer, however, from artifacts due to unfulfilled assumptions concerning breathing pattern regularity. We propose and evaluate modifications to existing low-pitch spiral 4D CT reconstruction protocols to counteract respective artifacts.

Methods

The proposed advanced reconstruction (AR) approach consists of two steps that build on each other: (1) statistical analysis of the breathing signal recorded during CT data acquisition and extraction of a patient-specific reference breathing cycle for projection binning; (2) incorporation of an artifact measure into the reconstruction. 4D CT data of 30 patients were reconstructed by standard phase- and local amplitude-based reconstruction (PB, LAB) and compared with images obtained by AR. The number of artifacts was evaluated and artifact statistics correlated to breathing curve characteristics.

Results

AR reduced the number of 4D CT artifacts by 31% and 27% compared to PB and LAB; the reduction was most pronounced for irregular breathing curves.

Conclusions

We described a two-step optimization of low-pitch spiral 4D CT reconstruction to reduce artifacts in the presence of breathing irregularity and illustrated that the modifications to existing reconstruction solutions are effective in terms of artifact reduction.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-017-0835-7) contains supplementary material, which is available to authorized users.

Evaluation of gantry speed on image quality and imaging dose for 4D cone-beam CT acquisition

Background

This study investigates the effect of gantry speed on 4DCBCT image quality and dose for the Varian On-Board Imager®.

Methods

A thoracic 4DCBCT protocol was designed using a 125 kVp spectrum. Image quality parameters were evaluated for 4DCBCT acquisition using Catphan® phantom with real-time position management™ system for gantry speeds varying between 1.0 to 6.0°/s. Superior-inferior motion of the phantom was executed using a sinusoidal waveform with five second period. Scans were retrospectively sorted into 4 phases (CBCT-4 ph) and 10 phases (CBCT-10 ph); average 4DCBCT (CBCT-ave), using all image data from the 4DCBCT acquisitions was also evaluated. The 4DCBCT images were evaluated using the following image quality metrics: spatial resolution, contrast-to-noise ratio (CNR), and uniformity index (UI). Additionally, Hounsfield unit (HU) sensitivity compared to a baseline CBCT and percent differences and RMS errors (RMSE) of excursion were also determined. Imaging dose was evaluated using an IBA CC13 ion chamber placed within CIRS Thorax phantom using the same sinusoidal motion and image acquisition settings as mentioned above.

Results

Spatial resolution decreased linearly from 5.93 to 3.82 lp/cm as gantry speed increased from 1.0 to 6.0°/s. CNR decreased linearly from 4.80 to 1.82 with gantry speed increasing from 1.0 to 6.0°/s, respectively. No noteworthy variations in UI, HU sensitivity, or excursion metrics were observed with changes in gantry speed. Ion chamber dose rates measured ranged from 2.30 (lung) to 5.18 (bone) E-3 cGy/mAs.

Conclusions

A quantitative analysis of the Varian OBI’s 4DCBCT capabilities was explored. Changing gantry speed changes the number of projections used for reconstruction, affecting both image quality and imaging dose if x-ray tube current is held constant. From the results of this study, a gantry speed between 2 and 3°/s was optimal when considering image quality, dose, and reconstruction time. The future of 4DCBCT clinical utility relies on further investigation of image acquisition and reconstruction optimization.

Real time tracking in liver SBRT: comparison of CyberKnife and Vero by planning structure-based γ-evaluation and dose-area-histograms

The purpose of this study was to evaluate and compare two clinical tracking systems for radiosurgery with regard to their dosimetric and geometrical accuracy in liver SBRT: the robot-based CyberKnife and the gimbal-based Vero. Both systems perform real-time tumour tracking by correlating internal tumour and external surrogate motion. CyberKnife treatment plans were delivered to a high resolution 2D detector array mounted on a 4D motion platform, with the platform simulating (a) tumour motion trajectories extracted from the corresponding CyberKnife predictor log files and (b) the tumour motion trajectories with superimposed baseline-drift. Static reference and tracked dose measurements were compared and dosimetric as well as geometrical uncertainties analyzed by a planning structure-based evaluation. For (a), γ-passing rates inside the CTV (γ-criteria of 1% / 1 mm) ranged from 95% to 100% (CyberKnife) and 98% to 100% (Vero). However, dosimetric accuracy decreases in the presence of the baseline-drift. γ-passing rates for (b) ranged from 26% to 92% and 94% to 99%, respectively; i.e. the effect was more pronounced for CyberKnife. In contrast, the Vero system led to maximum dose deviations in the OAR between  +1.5 Gy to +6.0 Gy (CyberKnife: +0.5 Gy to +3.5 Gy). Potential dose shifts were interpreted as motion-induced geometrical tracking errors. Maximum observed shift ranges were  -1.0 mm to  +0.7 mm (lateral) /-0.6 mm to +0.1 mm (superior-inferior) for CyberKnife and  -0.8 mm to +0.2 mm /-0.8 mm to +0.4 mm for Vero. These values illustrate that CyberKnife and Vero provide high precision tracking of regular breathing patterns. Even for the modified motion trajectory, the obtained dose distributions appear to be clinical acceptable with regard to literature QA γ-criteria of 3% / 3 mm.

Advances in 4D treatment planning for scanned particle beam therapy - report of dedicated workshops

We report on recent progress in the field of mobile tumor treatment with scanned particle beams, as discussed in the latest editions of the 4D treatment planning workshop. The workshop series started in 2009, with about 20 people from 4 research institutes involved, all actively working on particle therapy delivery and development. The first workshop resulted in a summary of recommendations for the treatment of mobile targets, along with a list of requirements to apply these guidelines clinically. The increased interest in the treatment of mobile tumors led to a continuously growing number of attendees: the 2012 edition counted more than 60 participants from 20 institutions and commercial vendors. The focus of research discussions among workshop participants progressively moved from 4D treatment planning to complete 4D treatments, aiming at effective and safe treatment delivery. Current research perspectives on 4D treatments include all critical aspects of time resolved delivery, such as in-room imaging, motion detection, beam application, and quality assurance techniques. This was motivated by the start of first clinical treatments of hepato cellular tumors with a scanned particle beam, relying on gating or abdominal compression for motion mitigation. Up to date research activities emphasize significant efforts in investigating advanced motion mitigation techniques, with a specific interest in the development of dedicated tools for experimental validation. Potential improvements will be made possible in the near future through 4D optimized treatment plans that require upgrades of the currently established therapy control systems for time resolved delivery. But since also these novel optimization techniques rely on the validity of the 4DCT, research focusing on alternative 4D imaging technique, such as MRI based 4DCT generation will continue.

Toward a planning scheme for emission guided radiation therapy (EGRT): FDG based tumor tracking in a metastatic breast cancer patient

PURPOSE

Emission guided radiation therapy (EGRT) is a new modality that uses PET emissions in real-time for direct tumor tracking during radiation delivery. Radiation beamlets are delivered along positron emission tomography (PET) lines of response (LORs) by a fast rotating ring therapy unit consisting of a linear accelerator (Linac) and PET detectors. The feasibility of tumor tracking and a primitive modulation method to compensate for attenuation have been demonstrated using a 4D digital phantom in our prior work. However, the essential capability of achieving dose modulation as in conventional intensity modulated radiation therapy (IMRT) treatments remains absent. In this work, the authors develop a planning scheme for EGRT to accomplish sophisticated intensity modulation based on an IMRT plan while preserving tumor tracking.

METHODS

The planning scheme utilizes a precomputed LOR response probability distribution to achieve desired IMRT planning modulation with effects of inhomogeneous attenuation and nonuniform background activity distribution accounted for. Evaluation studies are performed on a 4D digital patient with a simulated lung tumor and a clinical patient who has a moving breast cancer metastasis in the lung. The Linac dose delivery is simulated using a voxel-based Monte Carlo algorithm. The IMRT plan is optimized for a planning target volume (PTV) that encompasses the tumor motion using the MOSEK package and a Pinnacle3™ workstation (Philips Healthcare, Fitchburg, WI) for digital and clinical patients, respectively. To obtain the emission data for both patients, the Geant4 application for tomographic emission (GATE) package and a commercial PET scanner are used. As a comparison, 3D and helical IMRT treatments covering the same PTV based on the same IMRT plan are simulated.

RESULTS

3D and helical IMRT treatments show similar dose distribution. In the digital patient case, compared with the 3D IMRT treatment, EGRT achieves a 15.1% relative increase in dose to 95% of the gross tumor volume (GTV) and a 31.8% increase to 50% of the GTV. In the patient case, EGRT yields a 15.2% relative increase in dose to 95% of the GTV and a 20.7% increase to 50% of the GTV. The organs at risk (OARs) doses are kept similar or lower for EGRT in both cases. Tumor tracking is observed in the presence of planning modulation in all EGRT treatments.

CONCLUSIONS

As compared to conventional IMRT treatments, the proposed EGRT planning scheme allows an escalated target dose while keeping dose to the OARs within the same planning limits. With the capabilities of incorporating planning modulation and accurate tumor tracking, EGRT has the potential to greatly improve targeting in radiation therapy and enable a practical and effective implementation of 4D radiation therapy for planning and delivery.