3-D lung deformation and function from respiratory-gated 4-D x-ray CT images : application to radiation treatment planning.

Many lung diseases or injuries can cause biomechanical or material property changes that can alter lung function. While the mechanical changes associated with the change of the material properties originate at a regional level, they remain largely asymptomatic and are invisible to global measures of lung function until they have advanced significantly and have aggregated. In the realm of external beam radiation therapy of patients suffering from lung cancer, determination of patterns of pre- and post-treatment motion, and measures of regional and global lung elasticity and function are clinically relevant. In this dissertation, we demonstrate that 4-D CT derived ventilation images, including mechanical strain, provide an accurate and physiologically relevant assessment of regional pulmonary function which may be incorporated into the treatment planning process. Our contributions are as follows: (i) A new volumetric deformable image registration technique based on 3-D optical flow (MOFID) has been designed and implemented which permits the possibility of enforcing physical constraints on the numerical solutions for computing motion field from respiratory-gated 4-D CT thoracic images. The proposed optical flow framework is an accurate motion model for the thoracic CT registration problem. (ii) A large displacement landmark-base elastic registration method has been devised for thoracic CT volumetric image sets containing large deformations or changes, as encountered for example in registration of pre-treatment and post-treatment images or multi-modality registration. (iii) Based on deformation maps from MOFIO, a novel framework for regional quantification of mechanical strain as an index of lung functionality has been formulated for measurement of regional pulmonary function. (iv) In a cohort consisting of seven patients with non-small cell lung cancer, validation of physiologic accuracy of the 4-0 CT derived quantitative images including Jacobian metric of ventilation, Vjac, and principal strains, (V?1, V?2, V?3, has been performed through correlation of the derived measures with SPECT ventilation and perfusion scans. The statistical correlations with SPECT have shown that the maximum principal strain pulmonary function map derived from MOFIO, outperforms all previously established ventilation metrics from 40-CT. It is hypothesized that use of CT -derived ventilation images in the treatment planning process will help predict and prevent pulmonary toxicity due to radiation treatment. It is also hypothesized that measures of regional and global lung elasticity and function obtained during the course of treatment may be used to adapt radiation treatment. Having objective methods with which to assess pre-treatment global and regional lung function and biomechanical properties, the radiation treatment dose can potentially be escalated to improve tumor response and local control

Investigation of time-resolved volumetric MRI to enhance MR-guided radiotherapy of moving lung tumors

In photon radiotherapy of lung cancer, respiratory-induced motion introduces systematic and statistical uncertainties in treatment planning and dose delivery. By integrating magnetic resonance imaging (MRI) in the treatment planning process in MR-guided radiotherapy (MRgRT), uncertainties in target volume definition can be reduced with respect to state-of-the-art X-ray-based approaches. Furthermore, MR-guided linear accelerators (MR-Linacs) offer dose delivery with enhanced accuracy and precision through daily treatment plan adaptation and gated beam delivery based on real-time MRI. Today, the potential of MRgRT of moving targets is, however, not fully exploited due to the lack of time-resolved four-dimensional MRI (4D-MRI) in clinical practice. Therefore, the aim of this thesis was to develop and experimentally validate new methods for motion characterization and estimation with 4D-MRI for MRgRT of lung cancer. Different concepts were investigated for all phases of the clinical workflow - treatment planning, beam delivery, and post-treatment analysis. Firstly, a novel internal target volume (ITV) definition method based on the probability-of-presence of moving tumors derived from real-time 4D-MRI was developed. The ability of the ITVs to prospectively account for changes occurring over the course of several weeks was assessed in retrospective geometric analyses of lung cancer patient data. Higher robustness of the probabilistic 4D-MRI-based ITVs against interfractional changes was observed compared to conventional target volumes defined with four-dimensional computed tomography (4D-CT). The study demonstrated that motion characterization over extended times enabled by real-time 4D-MRI can reduce systematic and statistical uncertainties associated with today’s standard workflow. Secondly, experimental validation of a published motion estimation method - the propagation method - was conducted with a porcine lung phantom under realistic patient-like conditions. Estimated 4D-MRIs with a temporal resolution of 3.65 Hz were created based on orthogonal 2D cine MRI acquired at the scanner unit of an MR-Linac. A comparison of these datasets with ground truth respiratory-correlated 4D-MRIs in geometric analyses showed that the propagation method can generate geometrically accurate estimated 4D-MRIs. These could decrease target localization errors and enable 3D motion monitoring during beam delivery at the MR-Linac in the future. Lastly, the propagation method was extended to create continuous time-resolved estimated synthetic CTs (tresCTs). The proposed method was experimentally tested with the porcine lung phantom, successively imaged at a CT scanner and an MR-Linac. A high agreement of the images and corresponding dose distributions of the tresCTs and measured ground truth 4D-CTs was found in geometric and dosimetric analyses. The tresCTs could be used for post-treatment time-resolved reconstruction of the delivered dose to guide treatment adaptations in the future. These studies represent important steps towards a clinical application of time-resolved 4D-MRI methods for enhanced MRgRT of lung tumors in the near future

Motion Compensation for Free-Breathing Abdominal Diffusion-Weighted Imaging (MoCo DWI)

Diffusion-weighted imaging (DWI) is a common technique in medical diagnostics. One challenge of thoracic and abdominal DWI is respiratory motion which can result in motion artifacts. To eliminate these artifacts, a new kind of retrospective, respiratory motion compensation for DWI was developed and tested. This new technique — MoCo DWI — is the first in DWI which provides fully-deformable motion compensation. To enable this, despite the low image quality of DWI, two free-breathing sequences were used: (1) a gradient echo sequence (GRE) with a configuration for optimal respiratory motion estimation and (2) a DWI in a configuration of clinical interest. The DWI acquisition was gated into 10 motion phases. Each motion phase was then co-aligned with the motion estimation. The implementation was tested with eleven volunteers. The results showed that MoCo DWI can reduce motion blurring in single b-value images, especially at the liver-lung interface. The improvement of ADC-maps was even more prominent. Individual slices showed motion induced artifacts which could be reduced or even eliminated by MoCo DWI. This was also reflected by expected more homogeneous ADC values in the liver in all data sets. These results promise to reduce measurements with limited diagnostic value while keeping or increasing patient comfort

How extensive of a 4D dataset is needed to estimate cumulative dose distribution plan evaluation metrics in conformal lung therapy?a)

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134757/1/mp0624.pd

Doctor of Philosophy

dissertationX-ray computed tomography (CT) is a widely popular medical imaging technique that allows for viewing of in vivo anatomy and physiology. In order to produce high-quality images and provide reliable treatment, CT imaging requires the precise knowledge of t

Registration based respiratory motion models for use in lung radiotherapy.

Respiratory motion is a major factor contributing to errors and uncertainties in Radiotherapy (RT) treatment of lung tumours. Knowledge of this motion may improve the planning and delivery of RT treatment for lung cancer patients. This thesis develops and evaluates methods of building patient specific respiratory motion models. These relate the internal motion to respiratory parameters derived from an external surrogate signal that can be measured during data acquisition and treatment delivery. The models offer a number of advantages over current methods of imaging and analysing respiratory motion, in particular their ability to account for variations in the respiratory motion. Computer Tomography (CT) data is acquired over several respiratory cycles to sample some of the variation in the respiratory motion. B-spline registrations are used to recover the motion and deformation from the CT data. The models are then constructed by fitting functions that relate the registration results to the respiratory parameters. This thesis describes the CT data and respiratory parameters that have been used to construct the motion models. It details the registrations protocols used and evaluates their results. The initial models presented in the thesis relate the registration results to a single parameter, the phase of the respiratory cycle, and average out any variation in the respiratory motion. The later models relate the registration results to two respiratory parameters, with the intention of modelling some of the variation. A number of different functions are assessed for both the single and two parameter models. The results show that the models can predict the respiratory motion in the CT data very accurately (mean error < 1.4 mm). This thesis also discusses some of the uses of the motion models in RT and, in particular, explores the use of the motion models for 'tracking' respiratory motion while delivering intensity modulated RT.

IMAGE-BASED RESPIRATORY MOTION EXTRACTION AND RESPIRATION-CORRELATED CONE BEAM CT (4D-CBCT) RECONSTRUCTION

Accounting for respiration motion during imaging helps improve targeting precision in radiation therapy. Respiratory motion can be a major source of error in determining the position of thoracic and upper abdominal tumor targets during radiotherapy. Thus, extracting respiratory motion is a key task in radiation therapy planning. Respiration-correlated or four-dimensional CT (4DCT) imaging techniques have been recently integrated into imaging systems for verifying tumor position during treatment and managing respiration-induced tissue motion. The quality of the 4D reconstructed volumes is highly affected by the respiratory signal extracted and the phase sorting method used. This thesis is divided into two parts. In the first part, two image-based respiratory signal extraction methods are proposed and evaluated. Those methods are able to extract the respiratory signals from CBCT images without using external sources, implanted markers or even dependence on any structure in the images such as the diaphragm. The first method, called Local Intensity Feature Tracking (LIFT), extracts the respiratory signal depending on feature points extracted and tracked through the sequence of projections. The second method, called Intensity Flow Dimensionality Reduction (IFDR), detects the respiration signal by computing the optical flow motion of every pixel in each pair of adjacent projections. Then, the motion variance in the optical flow dataset is extracted using linear and non-linear dimensionality reduction techniques to represent a respiratory signal. Experiments conducted on clinical datasets showed that the respiratory signal was successfully extracted using both proposed methods and it correlates well with standard respiratory signals such as diaphragm position and the internal markers’ signal. In the second part of this thesis, 4D-CBCT reconstruction based on different phase sorting techniques is studied. The quality of the 4D reconstructed images is evaluated and compared for different phase sorting methods such as internal markers, external markers and image-based methods (LIFT and IFDR). Also, a method for generating additional projections to be used in 4D-CBCT reconstruction is proposed to reduce the artifacts that result when reconstructing from an insufficient number of projections. Experimental results showed that the feasibility of the proposed method in recovering the edges and reducing the streak artifacts

Characterisation and correction of respiratory-motion artefacts in cardiac PET-CT

Respiratory motion during cardiac Positron Emission Tomography (PET) Computed Tomography (CT) imaging results in blurring of the PET data and can induce mismatches between the PET and CT datasets, leading to attenuation-correction artefacts. The aim of this project was to develop a method of motion-correction to overcome both of these problems. The approach implemented was to transform a single CT to match the frames of a gated PET study, to facilitate respiratory-matched attenuation-correction, without the need for a gated CT. This is benecial for lowering the radiation dose to the patient and in reducing PETCT mismatches, which can arise even in gated studies. The heart and diaphragm were identied through phantom studies as the structures responsible for generating attenuation-correction artefacts in the heart and their motions therefore needed to be considered in transforming the CT. Estimating heart motion was straight-forward, due to its high contrast in PET, however the poor diaphragm contrast meant that additional information was required to track its position. Therefore a diaphragm shape model was constructed using segmented diaphragm surfaces, enabling complete diaphragm surfaces to be produced from incomplete and noisy initial estimates. These complete surfaces, in combination with the estimated heart motions were used to transform the CT. The PET frames were then attenuation-corrected with the transformed CT, reconstructed, aligned and summed, to produce motion-free images. It was found that motion-blurring was reduced through alignment, although benets were marginal in the presence of small respiratory motions. Quantitative accuracy was improved from use of the transformed CT for attenuation-correction (compared with no CT transformation), which was attributed to both the heart and the diaphragm transformations. In comparison to a gated CT, a substantial dose saving and a reduced dependence on gating techniques were achieved, indicating the potential value of the technique in routine clinical procedures

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

Segmentation, tracking, and kinematics of lung parenchyma and lung tumors from 4D CT with application to radiation treatment planning.

This thesis is concerned with development of techniques for efficient computerized analysis of 4-D CT data. The goal is to have a highly automated approach to segmentation of the lung boundary and lung nodules inside the lung. The determination of exact lung tumor location over space and time by image segmentation is an essential step to track thoracic malignancies. Accurate image segmentation helps clinical experts examine the anatomy and structure and determine the disease progress. Since 4-D CT provides structural and anatomical information during tidal breathing, we use the same data to also measure mechanical properties related to deformation of the lung tissue including Jacobian and strain at high resolutions and as a function of time. Radiation Treatment of patients with lung cancer can benefit from knowledge of these measures of regional ventilation. Graph-cuts techniques have been popular for image segmentation since they are able to treat highly textured data via robust global optimization, avoiding local minima in graph based optimization. The graph-cuts methods have been used to extract globally optimal boundaries from images by s/t cut, with energy function based on model-specific visual cues, and useful topological constraints. The method makes N-dimensional globally optimal segmentation possible with good computational efficiency. Even though the graph-cuts method can extract objects where there is a clear intensity difference, segmentation of organs or tumors pose a challenge. For organ segmentation, many segmentation methods using a shape prior have been proposed. However, in the case of lung tumors, the shape varies from patient to patient, and with location. In this thesis, we use a shape prior for tumors through a training step and PCA analysis based on the Active Shape Model (ASM). The method has been tested on real patient data from the Brown Cancer Center at the University of Louisville. We performed temporal B-spline deformable registration of the 4-D CT data - this yielded 3-D deformation fields between successive respiratory phases from which measures of regional lung function were determined. During the respiratory cycle, the lung volume changes and five different lobes of the lung (two in the left and three in the right lung) show different deformation yielding different strain and Jacobian maps. In this thesis, we determine the regional lung mechanics in the Lagrangian frame of reference through different respiratory phases, for example, Phase10 to 20, Phase10 to 30, Phase10 to 40, and Phase10 to 50. Single photon emission computed tomography (SPECT) lung imaging using radioactive tracers with SPECT ventilation and SPECT perfusion imaging also provides functional information. As part of an IRB-approved study therefore, we registered the max-inhale CT volume to both VSPECT and QSPECT data sets using the Demon\u27s non-rigid registration algorithm in patient subjects. Subsequently, statistical correlation between CT ventilation images (Jacobian and strain values), with both VSPECT and QSPECT was undertaken. Through statistical analysis with the Spearman\u27s rank correlation coefficient, we found that Jacobian values have the highest correlation with both VSPECT and QSPECT