Evaluation of dose prediction error and optimization convergence error in four-dimensional inverse planning of robotic stereotactic lung radiotherapy

Inverse optimization of robotic stereotactic lung radiotherapy is typically performed using relatively simple dose calculation algorithm on a single instance of breathing geometry. Variations of patient geometry and tissue density during respiration could reduce the dose accuracy of these 3D optimized plans. To quantify the potential benefits of direct four-dimensional (4D) optimization in robotic lung radiosurgery, 4D optimizations using 1) ray-tracing algorithm with equivalent path-length heterogeneity correction (4EPL(opt)), and 2) Monte Carlo (MC) algorithm (4MC(opt)), were performed in 25 patients. The 4EPL(opt) plans were recalculated using MC algorithm (4MC(recal)) to quantify the dose prediction errors (DPEs). Optimization convergence errors (OCEs) were evaluated by comparing the 4MC(recal) and 4MC(opt) dose results. The results were analyzed by dose-volume histogram indices for selected organs. Statistical equivalence tests were performed to determine the clinical significance of the DPEs and OCEs, compared with a 3% tolerance. Statistical equivalence tests indicated that the DPE and the OCE are significant predominately in GTV D98%. The DPEs in V20 of lung, and D2% of cord, trachea, and esophagus are within 1.2%, while the OCEs are within 10.4% in lung V20 and within 3.5% in trachea D2%. The marked DPE and OCE suggest that 4D MC optimization is important to improve the dosimetric accuracy in robotic-based stereotactic body radiotherapy, despite the longer computation time.published_or_final_versio

Optimization strategies for respiratory motion management in stereotactic body radiation therapy

Various challenges arise during the treatment of lung tumors with stereotactic body radiation therapy (SBRT), which is a form of hypofractionated high precision conformal radiation therapy delivered to small targets. The dose is applied in only a few fractions and respiratory organ and tumor motion is a source of uncertainty additional to interfractional set-up errors. Respiratory organ and tumor motion is highly patient-specific and it affects the whole radiotherapy treatment chain. In this thesis, motion management techniques for SBRT are evaluated and improved in a clinical setting. A clinical need for improvement has been present at the LMU university hospital for each issue addressed in this thesis: Initially, the usage of respiratory correlated computed tomography (4DCT), which is vital for SBRT treatment, was seen as impractical and prone to uncertainties in the data reconstruction in its current form. Therefore, the 4DCT reconstruction workflow has been improved to minimize these potential error sources. Secondly, treatment planning for tumors affected by respiratory motion was evaluated and subsequently improved. Finally, the treatment technique of respiratory gating was implemented at the clinic, which led to the need of evaluating the respiratory gating characteristics of the novel system configuration. At first, the 4DCT reconstruction workflow used in clinical practice was investigated, as in the presence of respiratory motion the knowledge of tumor position over time is essential in SBRT treatments. Using 4DCT, the full motion range of the individual tumor can be determined. However, certain 4DCT reconstruction methods can under- or overestimate tumor motion due to limitations in the data acquisition scheme and due to the incorrect sorting of certain X-ray computed tomography (CT) image slices into different respiratory phases. As the regular clinical workflow of cycle-based sorting (CBS) without maximum inspiration detection (and therefore no clear starting point for the individual breathing cycles) seemed to be affected by these potential errors, the usage of CBS with correct maximum detection and another sorting algorithm of the respiration states, so-called local amplitude-based sorting (LAS), both have been implemented for a reduction of image artifacts and improved 4DCT quality. The three phase binning algorithms have been investigated in a phantom study (using 10 different breathing waveforms) and in a patient study (with 10 different patients). The mis-representation of the tumor volume was reduced in both implemented sorting algorithms compared to the previously used CBS approach (without correct maximum detection) in the phantom and the patient study. The clinical recommendation was the use of CBS with improved maximum detection, as too many manual interventions would be needed for the LAS workflow. Secondly, a combination of the actual patient breathing trace during treatment, the log files generated by the linear accelerator (LINAC), and Monte Carlo (MC) four-dimensional (4D) dose calculations for each individual fraction was implemented as a 4D dose evaluation tool. This workflow was tested in a clinical environment for SBRT treatment planning on multiple CT datasets featuring: a native free-breathing 3DCT, an average intensity projection (AIP) as well as a maximum intensity projection (MIP), both obtained from the patient's 4DCT, and density overrides (DOs) in a 3DCT. This study has been carried out for 5 SBRT patients for three-dimensional conformal radiation therapy (3D-CRT) and volumetric modulated arc therapy (VMAT) treatment plans. The dose has been recalculated on each 4DCT breathing phase according the the patient's breathing waveform and accumulated to the gross tumor volume (GTV) at the end-of-exhale (EOE) breathing phase using deformable image registration. Even though the least differences in planned and recalculated dose were found for AIP and MIP treatment planning, the results indicate a strong dependency on individual tumor motion due to the variability of breathing motion in general, and on tumor size. The combination of the patient's individual breathing trace during each SBRT fraction with 4D MC dose calculation based on the LINAC log file information leads to a good approximation of actual dose delivery. Finally, in order to ensure precise and accurate treatment for respiratory gating techniques, the technical characteristics of the LINAC in combination with a breathing motion monitoring system as s surrogate for tumor motion have to be identified. The dose delivery accuracy and the latency of a surface imaging system in connection with a modern medical LINAC were investigated using a dynamic breathing motion phantom. The dosimetric evaluation has been carried out using a static 2D-diode array. The measurement of the dose difference between gated and ungated radiation delivery was found to be below 1% (for clinical relevant gating levels of about 30%). The beam-on latency, or time delay, determined using radiographic films was found to be up to 851 ms±100 ms. With these known parameters, an adjustment of the pre-selected gating level or the internal target volume (ITV) margins could be made. With the highly patient-specific character of respiratory motion, lung SBRT faces many additional challenges besides the specific issues addressed in this thesis. However, the findings of this thesis have improved clinical workflows at the Department of Radiation Oncology of the LMU University hospital. In a future perspective, a workflow using evaluation of the actual 4D dose in combination with accurate 4DCT image acquisition and specialized treatment delivery (such as respiratory gating) has the potential for a safe further reduction of treatment margins and increased sparing of organs-at-risk (OARs) in SBRT without compromising tumor dose targeting accuracy

Optimization strategies for respiratory motion management in stereotactic body radiation therapy

Various challenges arise during the treatment of lung tumors with stereotactic body radiation therapy (SBRT), which is a form of hypofractionated high precision conformal radiation therapy delivered to small targets. The dose is applied in only a few fractions and respiratory organ and tumor motion is a source of uncertainty additional to interfractional set-up errors. Respiratory organ and tumor motion is highly patient-specific and it affects the whole radiotherapy treatment chain. In this thesis, motion management techniques for SBRT are evaluated and improved in a clinical setting. A clinical need for improvement has been present at the LMU university hospital for each issue addressed in this thesis: Initially, the usage of respiratory correlated computed tomography (4DCT), which is vital for SBRT treatment, was seen as impractical and prone to uncertainties in the data reconstruction in its current form. Therefore, the 4DCT reconstruction workflow has been improved to minimize these potential error sources. Secondly, treatment planning for tumors affected by respiratory motion was evaluated and subsequently improved. Finally, the treatment technique of respiratory gating was implemented at the clinic, which led to the need of evaluating the respiratory gating characteristics of the novel system configuration. At first, the 4DCT reconstruction workflow used in clinical practice was investigated, as in the presence of respiratory motion the knowledge of tumor position over time is essential in SBRT treatments. Using 4DCT, the full motion range of the individual tumor can be determined. However, certain 4DCT reconstruction methods can under- or overestimate tumor motion due to limitations in the data acquisition scheme and due to the incorrect sorting of certain X-ray computed tomography (CT) image slices into different respiratory phases. As the regular clinical workflow of cycle-based sorting (CBS) without maximum inspiration detection (and therefore no clear starting point for the individual breathing cycles) seemed to be affected by these potential errors, the usage of CBS with correct maximum detection and another sorting algorithm of the respiration states, so-called local amplitude-based sorting (LAS), both have been implemented for a reduction of image artifacts and improved 4DCT quality. The three phase binning algorithms have been investigated in a phantom study (using 10 different breathing waveforms) and in a patient study (with 10 different patients). The mis-representation of the tumor volume was reduced in both implemented sorting algorithms compared to the previously used CBS approach (without correct maximum detection) in the phantom and the patient study. The clinical recommendation was the use of CBS with improved maximum detection, as too many manual interventions would be needed for the LAS workflow. Secondly, a combination of the actual patient breathing trace during treatment, the log files generated by the linear accelerator (LINAC), and Monte Carlo (MC) four-dimensional (4D) dose calculations for each individual fraction was implemented as a 4D dose evaluation tool. This workflow was tested in a clinical environment for SBRT treatment planning on multiple CT datasets featuring: a native free-breathing 3DCT, an average intensity projection (AIP) as well as a maximum intensity projection (MIP), both obtained from the patient's 4DCT, and density overrides (DOs) in a 3DCT. This study has been carried out for 5 SBRT patients for three-dimensional conformal radiation therapy (3D-CRT) and volumetric modulated arc therapy (VMAT) treatment plans. The dose has been recalculated on each 4DCT breathing phase according the the patient's breathing waveform and accumulated to the gross tumor volume (GTV) at the end-of-exhale (EOE) breathing phase using deformable image registration. Even though the least differences in planned and recalculated dose were found for AIP and MIP treatment planning, the results indicate a strong dependency on individual tumor motion due to the variability of breathing motion in general, and on tumor size. The combination of the patient's individual breathing trace during each SBRT fraction with 4D MC dose calculation based on the LINAC log file information leads to a good approximation of actual dose delivery. Finally, in order to ensure precise and accurate treatment for respiratory gating techniques, the technical characteristics of the LINAC in combination with a breathing motion monitoring system as s surrogate for tumor motion have to be identified. The dose delivery accuracy and the latency of a surface imaging system in connection with a modern medical LINAC were investigated using a dynamic breathing motion phantom. The dosimetric evaluation has been carried out using a static 2D-diode array. The measurement of the dose difference between gated and ungated radiation delivery was found to be below 1% (for clinical relevant gating levels of about 30%). The beam-on latency, or time delay, determined using radiographic films was found to be up to 851 ms±100 ms. With these known parameters, an adjustment of the pre-selected gating level or the internal target volume (ITV) margins could be made. With the highly patient-specific character of respiratory motion, lung SBRT faces many additional challenges besides the specific issues addressed in this thesis. However, the findings of this thesis have improved clinical workflows at the Department of Radiation Oncology of the LMU University hospital. In a future perspective, a workflow using evaluation of the actual 4D dose in combination with accurate 4DCT image acquisition and specialized treatment delivery (such as respiratory gating) has the potential for a safe further reduction of treatment margins and increased sparing of organs-at-risk (OARs) in SBRT without compromising tumor dose targeting accuracy

Medical physics challenges in clinical MR-guided radiotherapy

The integration of magnetic resonance imaging (MRI) for guidance in external beam radiotherapy has faced significant research and development efforts in recent years. The current availability of linear accelerators with an embedded MRI unit, providing volumetric imaging at excellent soft tissue contrast, is expected to provide novel possibilities in the implementation of image-guided adaptive radiotherapy (IGART) protocols. This study reviews open medical physics issues in MR-guided radiotherapy (MRgRT) implementation, with a focus on current approaches and on the potential for innovation in IGART.Daily imaging in MRgRT provides the ability to visualize the static anatomy, to capture internal tumor motion and to extract quantitative image features for treatment verification and monitoring. Those capabilities enable the use of treatment adaptation, with potential benefits in terms of personalized medicine. The use of online MRI requires dedicated efforts to perform accurate dose measurements and calculations, due to the presence of magnetic fields. Likewise, MRgRT requires dedicated quality assurance (QA) protocols for safe clinical implementation.Reaction to anatomical changes in MRgRT, as visualized on daily images, demands for treatment adaptation concepts, with stringent requirements in terms of fast and accurate validation before the treatment fraction can be delivered. This entails specific challenges in terms of treatment workflow optimization, QA, and verification of the expected delivered dose while the patient is in treatment position. Those challenges require specialized medical physics developments towards the aim of fully exploiting MRI capabilities. Conversely, the use of MRgRT allows for higher confidence in tumor targeting and organs-at-risk (OAR) sparing.The systematic use of MRgRT brings the possibility of leveraging IGART methods for the optimization of tumor targeting and quantitative treatment verification. Although several challenges exist, the intrinsic benefits of MRgRT will provide a deeper understanding of dose delivery effects on an individual basis, with the potential for further treatment personalization

4D-CT Lung Registration and its Application for Lung Radiation Therapy

Radiation therapy has been successful in treating lung cancer patients, but its efficacy is limited by the inability to account for the respiratory motion during treatment planning and radiation dose delivery. Physics-based lung deformation models facilitate the motion computation of both tumor and local lung tissue during radiation therapy. In this dissertation, a novel method is discussed to accurately register 3D lungs across the respiratory phases from 4D-CT datasets, which facilitates the estimation of the volumetric lung deformation models. This method uses multi-level and multi-resolution optical flow registration coupled with thin plate splines (TPS), to address registration issue of inconsistent intensity across respiratory phases. It achieves higher accuracy as compared to multi-resolution optical flow registration and other commonly used registration methods. Results of validation show that the lung registration is computed with 3 mm Target Registration Error (TRE) and approximately 3 mm Inverse Consistency Error (ICE). This registration method is further implemented in GPU based real time dose delivery simulation to assist radiation therapy planning

Improving Dose-Response Correlations for Locally Advanced NSCLC Patients Treated with IMRT or PSPT

The standard of care for locally advanced non-small cell lung cancer (NSCLC) is concurrent chemo-radiotherapy. Despite recent advancements in radiation delivery methods, the median survival time of NSCLC patients remains below 28 months. Higher tumor dose has been found to increase survival but also a higher rate of radiation pneumonitis (RP) that affects breathing capability. In fear of such toxicity, less-aggressive treatment plans are often clinically preferred, leading to metastasis and recurrence. Therefore, accurate RP prediction is crucial to ensure tumor coverage to improve treatment outcome. Current models have associated RP with increased dose but with limited accuracy as they lack spatial correlation between accurate dose representation and quantitative RP representation. These models represent lung tissue damage with radiation dose distribution planned pre-treatment, which assumes a fixed patient geometry and inevitably renders imprecise dose delivery due to intra-fractional breathing motion and inter-fractional anatomy response. Additionally, current models employ whole-lung dose metrics as the contributing factor to RP as a qualitative, binary outcome but these global dose metrics discard microscopic, voxel-(3D pixel)-level information and prevent spatial correlations with quantitative RP representation. To tackle these limitations, we developed advanced deformable image registration (DIR) techniques that registered corresponding anatomical voxels between images for tracking and accumulating dose throughout treatment. DIR also enabled voxel-level dose-response correlation when CT image density change (IDC) was used to quantify RP. We hypothesized that more accurate estimates of biologically effective dose distributions actually delivered, achieved through (a) dose accumulation using deformable registration of weekly 4DCT images acquired over the course or radiotherapy and (b) the incorporation of variable relative biological effectiveness (RBE), would lead to statistically and clinically significant improvement in the correlation of RP with biologically effective dose distributions. Our work resulted in a robust intra-4DCT and inter-4DCT DIR workflow, with the accuracy meeting AAPM TG-132 recommendations for clinical implementation of DIR. The automated DIR workflow allowed us to develop a fully automated 4DCT-based dose accumulation pipeline in RayStation (RaySearch Laboratories, Stockholm, Sweden). With a sample of 67 IMRT patients, our results showed that the accumulated dose was statistically different than the planned dose across the entire cohort with an average MLD increase of ~1 Gy and clinically different for individual patients where 16% resulted in difference in the score of the normal tissue complication probability (NTCP) using an established, clinically used model, which could qualify the patients for treatment planning re-evaluation. Lastly, we associated dose difference with accuracy difference by establishing and comparing voxel-level dose-IDC correlations and concluded that the accumulated dose better described the localized damage, thereby a closer representation of the delivered dose. Using the same dose-response correlation strategy, we plotted the dose-IDC relationships for both photon patients (N = 51) and proton patients (N = 67), we measured the variable proton RBE values to be 3.07–1.27 from 9–52 Gy proton voxels. With the measured RBE values, we fitted an established variable proton RBE model with pseudo-R2 of 0.98. Therefore, our results led to statistically and clinically significant improvement in the correlation of RP with accumulated and biologically effective dose distributions and demonstrated the potential of incorporating the effect of anatomical change and biological damage in RP prediction models

Effect of deformable registration uncertainty on lung SBRT dose accumulation

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134767/1/mp8412.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/134767/2/mp8412_am.pd