Biomechanical Modeling for Lung Tumor Motion Prediction during Brachytherapy and Radiotherapy

A novel technique is proposed to develop a biomechanical model for estimating lung’s tumor position as a function of respiration cycle time. Continuous tumor motion is a major challenge in lung cancer treatment techniques where the tumor needs to be targeted; e.g. in external beam radiotherapy and brachytherapy. If not accounted for, this motion leads to areas of radiation over and/or under dosage for normal tissue and tumors. In this thesis, biomechanical models were developed for lung tumor motion predication in two distinct cases of lung brachytherapy and lung external beam radiotherapy. The lung and other relevant surrounding organs geometries, loading, boundary conditions and mechanical properties were considered and incorporated properly for each case. While using material model with constant incompressibility is sufficient to model the lung tissue in the brachytherapy case, in external beam radiation therapy the tissue incompressibility varies significantly due to normal breathing. One of the main issues tackled in this research is characterizing lung tissue incompressibility variations and measuring its corresponding parameters as a function of respiration cycle time. Results obtained from an ex-vivo porcine deflated lung indicated feasibility and reliability of using the developed biomechanical model to predict tumor motion during brachytherapy. For external beam radiotherapy, in-silico studies indicated very significant impact of considering the lung tissue incompressibility on the accuracy of predicting tumor motion. Furthermore, ex-vivo porcine lung experiments demonstrated the capability and reliability of the proposed approach for predicting tumor motion as a function of cyclic time. As such, the proposed models have a good potential to be incorporated effectively in computer assisted lung radiotherapy treatment systems

A biomechanical approach for real-time tracking of lung tumors during External Beam Radiation Therapy (EBRT)

Lung cancer is the most common cause of cancer related death in both men and women. Radiation therapy is widely used for lung cancer treatment. However, this method can be challenging due to respiratory motion. Motion modeling is a popular method for respiratory motion compensation, while biomechanics-based motion models are believed to be more robust and accurate as they are based on the physics of motion. In this study, we aim to develop a biomechanics-based lung tumor tracking algorithm which can be used during External Beam Radiation Therapy (EBRT). An accelerated lung biomechanical model can be used during EBRT only if its boundary conditions (BCs) are defined in a way that they can be updated in real-time. As such, we have developed a lung finite element (FE) model in conjunction with a Neural Networks (NNs) based method for predicting the BCs of the lung model from chest surface motion data. To develop the lung FE model for tumor motion prediction, thoracic 4D CT images of lung cancer patients were processed to capture the lung and diaphragm geometry, trans-pulmonary pressure, and diaphragm motion. Next, the chest surface motion was obtained through tracking the motion of the ribcage in 4D CT images. This was performed to simulate surface motion data that can be acquired using optical tracking systems. Finally, two feedforward NNs were developed, one for estimating the trans-pulmonary pressure and another for estimating the diaphragm motion from chest surface motion data. The algorithm development consists of four steps of: 1) Automatic segmentation of the lungs and diaphragm, 2) diaphragm motion modelling using Principal Component Analysis (PCA), 3) Developing the lung FE model, and 4) Using two NNs to estimate the trans-pulmonary pressure values and diaphragm motion from chest surface motion data. The results indicate that the Dice similarity coefficient between actual and simulated tumor volumes ranges from 0.76±0.04 to 0.91±0.01, which is favorable. As such, real-time lung tumor tracking during EBRT using the proposed algorithm is feasible. Hence, further clinical studies involving lung cancer patients to assess the algorithm performance are justified

A Free-Breathing Lung Motion Model

Lung cancer has been the leading cause of cancer deaths for decades in the United States. Although radiotherapy is one of the most effective treatments, side effects from error in delivery of radiation due to organ motion during breathing remain a significant issue. To compensate the breathing motion during the treatment, a free breathing lung motion model, x= x0+αv+βf, was developed and discussed, where x is the position of a piece of tissue located at reference position x0. α is a parameter which characterizes the motion due to local air filling: motion as a function of tidal volume) and β is the parameter that accounts for the motion due to the imbalance of dynamical stress distributions during inspiration and exhalation which cause lung motion hysteresis: motion as a function of airflow). The parameters α and β together provide a quantitative characterization of breathing motion that inherently includes the complex hysteresis interplay. The theoretical foundation of the model was built by investigating the stress distribution inside of a lung and the biomechanical properties of the lung tissues. Accuracy of the model was investigated by using 49 free-breathing patient data sets. Applications of the model in localizing lung cancer, monitoring radiation damage and suppressing artifacts in free-breathing PET images were also discussed. This work supported in part by NIHR01CA096679 and NIHR01CA11671

2012 Activity Report of the Regional Research Programme on Hadrontherapy for the ETOILE Center

2012 is the penultimate year of financial support by the CPER 2007-2013 for ETOILE's research program, sustained by the PRRH at the University Claude Bernard. As with each edition we make the annual review of the research in this group, so active for over 12 years now. Over the difficulties in the decision-making process for the implementation of the ETOILE Center, towards which all our efforts are focussed, some "themes" (work packages) were strengthened, others have progressed, or have been dropped. This is the case of the eighth theme (technological developments), centered around the technology for rotative beam distribution heads (gantries) and, after being synchronized with the developments of ULICE's WP6, remained so by ceasing its activities, coinciding also with the retirement of its historic leader at IPNL, Marcel Bajard. Topic number 5 ("In silico simulations") has suffered the departure of its leader, Benjamin Ribba, although the work has still been provided by Branka Bernard, a former postdoctoral fellow in Lyon Sud, and now back home in Croatia, still in contract with UCBL for the ULICE project. Aside from these two issues (and the fact that the theme "Medico-economical simulations" is now directly linked to the first one ("Medical Project"), the rest of the teams are growing, as evidenced by the publication statistics at the beginning of this report. This is obviously due to the financial support of our always faithful regional institutions, but also to the synergy that the previous years, the European projects, the arrival of the PRIMES LabEx, and the national France Hadron infrastructure have managed to impulse. The Rhone-Alpes hadron team, which naturally includes the researchers of LPC at Clermont, should also see its influence result in a strong presence in France Hadron's regional node, which is being organized. The future of this regional research is not yet fully guaranteed, especially in the still uncertain context of ETOILE, but the tracks are beginning to emerge to allow past and present efforts translate into a long future that we all want to see established. Each of the researchers in PRRH is aware that 2013 will be (and already is) the year of great challenge : for ETOILE, for the PRRH, for hadron therapy in France, for French hadrontherapy in Europe (after the opening and beginning of treatments in the German [HIT Heidelberg, Marburg], Italian [CNAO, Pavia] and Austrian [MedAustron, Wien Neuerstadt]) centers. Let us meet again in early 2014 for a comprehensive review of the past and a perspective for the future ..

A Composite Material-based Computational Model for Diaphragm Muscle Biomechanical Simulation

Lung cancer is the most common cause of cancer related death among both men and women. Radiation therapy is the most widely used treatment for this disease. Motion compensation for tumor movement is often clinically important and biomechanics-based motion models may provide the most robust method as they are based on the physics of motion. In this study, we aim to develop a patient specific biomechanical model that predicts the deformation field of the diaphragm muscle during respiration. The first part of the project involved developing an accurate and adaptable micro-to-macro mechanical approach for skeletal muscle tissue modelling for application in a FE solver. The next objective was to develop the FE-based mechanical model of the diaphragm muscle based on patient specific 4D-CT data. The model shows adaptability to pathologies and may have the potential to be incorporated into respiratory models for the aid in treatment and diagnosis of diseases

A Heterogeneous Patient-Specific Biomechanical Model of the Lung for Tumor Motion Compensation and Effective Lung Radiation Therapy Planning

Radiation therapy is a main component of treatment for many lung cancer patients. However, the respiratory motion can cause inaccuracies in radiation delivery that can lead to treatment complications. In addition, the radiation-induced damage to healthy tissue limits the effectiveness of radiation treatment. Motion management methods have been developed to increase the accuracy of radiation delivery, and functional avoidance treatment planning has emerged to help reduce the chances of radiation-induced toxicity. In this work, we have developed biomechanical model-based techniques for tumor motion estimation, as well as lung functional imaging. The proposed biomechanical model accurately estimates lung and tumor motion/deformation by mimicking the physiology of respiration, while accounting for heterogeneous changes in the lung mechanics caused by COPD, a common lung cancer comorbidity. A biomechanics-based image registration algorithm is developed and is combined with an air segmentation algorithm to develop a 4DCT-based ventilation imaging technique, with potential applications in functional avoidance therapies

Inverse-Consistent Determination of Young\u27s Modulus of Human Lung

Human lung undergoes respiration-induced deformation due to sequential inhalation and exhalation. Accurate determination of lung deformation is crucial for tumor localization and targeted radiotherapy in patients with lung cancer. Numerical modeling of human lung dynamics based on underlying physics and physiology enables simulation and virtual visualization of lung deformation. Dynamical modeling is numerically complicated by the lack of information on lung elastic behavior, structural heterogeneity as well as boundary constrains. This study integrates physics-based modeling and image-based data acquisition to develop the patient-specific biomechanical model and consequently establish the first consistent Young\u27s modulus (YM) of human lung. This dissertation has four major components: (i) develop biomechanical model for computation of the flow and deformation characteristics that can utilize subject-specific, spatially-dependent lung material property; (ii) develop a fusion algorithm to integrate deformation results from a deformable image registration (DIR) and physics-based modeling using the theory of Tikhonov regularization; (iii) utilize fusion algorithm to establish unique and consistent patient specific Young\u27s modulus and; (iv) validate biomechanical model utilizing established patient-specific elastic property with imaging data. The simulation is performed on three dimensional lung geometry reconstructed from four-dimensional computed tomography (4DCT) dataset of human subjects. The heterogeneous Young\u27s modulus is estimated from a linear elastic deformation model with the same lung geometry and 4D lung DIR. The biomechanical model adequately predicts the spatio-temporal lung deformation, consistent with data obtained from imaging. The accuracy of the numerical solution is enhanced through fusion with the imaging data beyond the classical comparison of the two sets of data. Finally, the fused displacement results are used to establish unique and consistent patient-specific elastic property of the lung

Applications of a Biomechanical Patient Model for Adaptive Radiation Therapy

Biomechanical patient modeling incorporates physical knowledge of the human anatomy into the image processing that is required for tracking anatomical deformations during adaptive radiation therapy, especially particle therapy. In contrast to standard image registration, this enforces bio-fidelic image transformation. In this thesis, the potential of a kinematic skeleton model and soft tissue motion propagation are investigated for crucial image analysis steps in adaptive radiation therapy. The first application is the integration of the kinematic model in a deformable image registration process (KinematicDIR). For monomodal CT scan pairs, the median target registration error based on skeleton landmarks, is smaller than (1.6 ± 0.2) mm. In addition, the successful transferability of this concept to otherwise challenging multimodal registration between CT and CBCT as well as CT and MRI scan pairs is shown to result in median target registration error in the order of 2 mm. This meets the accuracy requirement for adaptive radiation therapy and is especially interesting for MR-guided approaches. Another aspect, emerging in radiotherapy, is the utilization of deep-learning-based organ segmentation. As radiotherapy-specific labeled data is scarce, the training of such methods relies heavily on augmentation techniques. In this work, the generation of synthetically but realistically deformed scans used as Bionic Augmentation in the training phase improved the predicted segmentations by up to 15% in the Dice similarity coefficient, depending on the training strategy. Finally, it is shown that the biomechanical model can be built-up from automatic segmentations without deterioration of the KinematicDIR application. This is essential for use in a clinical workflow