The uncertainties due to respiratory motion present significant challenges to accurate characterization of cancerous tissues both in terms of imaging and treatment. Currently available clinical lung imaging techniques are subject to inferior image quality and incorrect motion estimation, with consequences that can systematically impact the downstream treatment delivery and outcome. The main objective of this thesis is the development of the techniques of fast helical computed tomography (CT) imaging and deformable image registration for the radiotherapy applications in accurate breathing motion modeling, lung tissue density modeling and ventilation imaging.
Fast helical CT scanning was performed on 64-slice CT scanner using the shortest available gantry rotation time and largest pitch value such that scanning of the thorax region amounts to just two seconds, which is less than typical breathing cycle in humans. The scanning was conducted under free breathing condition. Any portion of the lung anatomy undergoing such scanning protocol would be irradiated for only a quarter second, effectively removing any motion induced image artifacts. The resulting CT data were pristine volumetric images that record the lung tissue position and density in a fraction of the breathing cycle. Following our developed protocol, multiple fast helical CT scans were acquired to sample the tissue positions in different breathing states. To measure the tissue displacement, deformable image registration was performed that registers the non-reference images to the reference one.
In modeling breathing motion, external breathing surrogate signal was recorded synchronously with the CT image slices. This allowed for the tissue-specific displacement to be modeled as parametrization of the recorded breathing signal using the 5D lung motion model. To assess the accuracy of the motion model in describing tissue position change, the model was used to simulate the original high-pitch helical CT scan geometries, employed as ground truth data. Image similarity between the simulated and ground truth scans was evaluated. The model validation experiments were conducted in a patient cohort of seventeen patients to assess the model robustness and inter-patient variation. The model error averaged over multiple tracked positions from several breathing cycles was found to be on the order of one millimeter.
In modeling the density change under free breathing condition, the determinant of Jacobian matrix from the registration-derived deformation vector field yielded volume change information of the lung tissues. Correlation of the Jacobian values to the corresponding voxel Housfield units (HU) reveals that the density variation for the majority of lung tissues can be very well described by mass conservation relationship. Different tissue types were identified and separately modeled. Large trials of validation experiments were performed. The averaged deviation between the modeled and the reference lung density was 30 HU, which was estimated to be the background CT noise level.
In characterizing the lung ventilation function, a novel method was developed to determine the extent of lung tissue volume change. Information on volume change was derived from the deformable image registration of the fast helical CT images in terms of Jacobian values with respect to a reference image. Assuming the multiple volume change measurements are independently and identically distributed, statistical formulation was derived to model ventilation distribution of each lung voxels and empirical minimum and maximum probability distribution of the Jacobian values was computed. Ventilation characteristic was evaluated as the difference of the expectation value from these extremal distributions. The resulting ventilation map was compared with an independently obtained ventilation image derived directly from the lung intensities and good correlation was found using statistical test. In addition, dynamic ventilation characterization was investigated by estimating the voxel-specific ventilation distribution. Ventilation maps were generated at different percentile levels using the tissue volume expansion metrics.