Heterogeneous volumetric data mapping and its medical applications

With the advance of data acquisition techniques, massive solid geometries are being collected routinely in scientific tasks, these complex and unstructured data need to be effectively correlated for various processing and analysis. Volumetric mapping solves bijective low-distortion correspondence between/among 3D geometric data, and can serve as an important preprocessing step in many tasks in compute-aided design and analysis, industrial manufacturing, medical image analysis, to name a few. This dissertation studied two important volumetric mapping problems: the mapping of heterogeneous volumes (with nonuniform inner structures/layers) and the mapping of sequential dynamic volumes. To effectively handle heterogeneous volumes, first, we studied the feature-aligned harmonic volumetric mapping. Compared to previous harmonic mapping, it supports the point, curve, and iso-surface alignment, which are important low-dimensional structures in heterogeneous volumetric data. Second, we proposed a biharmonic model for volumetric mapping. Unlike the conventional harmonic volumetric mapping that only supports positional continuity on the boundary, this new model allows us to have higher order continuity $C^1$ along the boundary surface. This suggests a potential model to solve the volumetric mapping of complex and big geometries through divide-and-conquer. We also studied the medical applications of our volumetric mapping in lung tumor respiratory motion modeling. We were building an effective digital platform for lung tumor radiotherapy based on effective volumetric CT/MRI image matching and analysis. We developed and integrated in this platform a set of geometric/image processing techniques including advanced image segmentation, finite element meshing, volumetric registration and interpolation. The lung organ/tumor and surrounding tissues are treated as a heterogeneous region and a dynamic 4D registration framework is developed for lung tumor motion modeling and tracking. Compared to the previous 3D pairwise registration, our new 4D parameterization model leads to a significantly improved registration accuracy. The constructed deforming model can hence approximate the deformation of the tissues and tumor

Digital Image Processing Applications

Digital image processing can refer to a wide variety of techniques, concepts, and applications of different types of processing for different purposes. This book provides examples of digital image processing applications and presents recent research on processing concepts and techniques. Chapters cover such topics as image processing in medical physics, binarization, video processing, and more

Surrogate-driven respiratory motion models for MRI-guided lung radiotherapy treatments

An MR-Linac integrates an MR scanner with a radiotherapy delivery system, providing non-ionizing real-time imaging of the internal anatomy before, during and after radiotherapy treatments. Due to spatio-temporal limitations of MR imaging, only high-resolution 2D cine-MR images can be acquired in real-time during MRI-guided radiotherapy (MRIgRT) to monitor the respiratory-induced motion of lung tumours and organs-at-risk. However, temporally-resolved 3D anatomical information is essential for accurate MR guidance of beam delivery and dose estimation of the actually delivered dose. Surrogate-driven respiratory motion models can estimate the 3D motion of the internal anatomy from surrogate signals, producing the required information. The overall aim of this thesis was to tailor a generalized respiratory motion modelling framework for lung MRIgRT. This framework can fit the model directly to unsorted 2D MR images sampling the 3D motion, and to surrogate signals extracted from the 2D cine-MR images acquired on an MR-Linac. It can model breath-to-breath variability and produce a motion compensated super-resolution reconstruction (MCSR) 3D image that can be deformed using the estimated motion. In this work novel MRI-derived surrogate signals were generated from 2D cine-MR images to model respiratory motion for lung cancer patients, by applying principal component analysis to the control point displacements obtained from the registration of the cine-MR images. An MR multi-slice interleaved acquisition potentially suitable for the MR-Linac was developed to generate MRI-derived surrogate signals and build accurate respiratory motion models with the generalized framework for lung cancer patients. The developed models and the MCSR images were thoroughly evaluated for lung cancer patients scanned on an MR-Linac. The results showed that respiratory motion models built with the generalized framework and minimal training data generally produced median errors within the MCSR voxel size of 2 mm, throughout the whole 3D thoracic field-of-view and over the expected lung MRIgRT treatment times

Statistical deformation reconstruction using multi-organ shape features for pancreatic cancer localization

Respiratory motion and the associated deformations of abdominal organs and tumors are essential information in clinical applications. However, inter- and intra-patient multi-organ deformations are complex and have not been statistically formulated, whereas single organ deformations have been widely studied. In this paper, we introduce a multi-organ deformation library and its application to deformation reconstruction based on the shape features of multiple abdominal organs. Statistical multi-organ motion/deformation models of the stomach, liver, left and right kidneys, and duodenum were generated by shape matching their region labels defined on four-dimensional computed tomography images. A total of 250 volumes were measured from 25 pancreatic cancer patients. This paper also proposes a per-region-based deformation learning using the non-linear kernel model to predict the displacement of pancreatic cancer for adaptive radiotherapy. The experimental results show that the proposed concept estimates deformations better than general per-patient-based learning models and achieves a clinically acceptable estimation error with a mean distance of 1.2 ± 0.7 mm and a Hausdorff distance of 4.2 ± 2.3 mm throughout the respiratory motion

Adaptive Radiation Therapy (ART) Strategies and Technical Considerations: A State of the ART Review From NRG Oncology

The integration of adaptive radiation therapy (ART), or modifying the treatment plan during the treatment course, is becoming more widely available in clinical practice. ART offers strong potential for minimizing treatment-related toxicity while escalating or de-escalating target doses based on the dose to organs at risk. Yet, ART workflows add complexity into the radiation therapy planning and delivery process that may introduce additional uncertainties. This work sought to review presently available ART workflows and technological considerations such as image quality, deformable image registration, and dose accumulation. Quality assurance considerations for ART components and minimum recommendations are described. Personnel and workflow efficiency recommendations are provided, as is a summary of currently available clinical evidence supporting the implementation of ART. Finally, to guide future clinical trial protocols, an example ART physician directive and a physics template following standard NRG Oncology protocol is provided

Segmentation of skin lesions in 2D and 3D ultrasound images using a spatially coherent generalized Rayleigh mixture model

This paper addresses the problem of jointly estimating the statistical distribution and segmenting lesions in multiple-tissue high-frequency skin ultrasound images. The distribution of multiple-tissue images is modeled as a spatially coherent finite mixture of heavy-tailed Rayleigh distributions. Spatial coherence inherent to biological tissues is modeled by enforcing local dependence between the mixture components. An original Bayesian algorithm combined with a Markov chain Monte Carlo method is then proposed to jointly estimate the mixture parameters and a label-vector associating each voxel to a tissue. More precisely, a hybrid Metropolis-within-Gibbs sampler is used to draw samples that are asymptotically distributed according to the posterior distribution of the Bayesian model. The Bayesian estimators of the model parameters are then computed from the generated samples. Simulation results are conducted on synthetic data to illustrate the performance of the proposed estimation strategy. The method is then successfully applied to the segmentation of in vivo skin tumors in high-frequency 2-D and 3-D ultrasound images

On the investigation of a novel x-ray imaging techniques in radiation oncology

Radiation therapy is indicated for nearly 50% of cancer patients in Australia. Radiation therapy requires accurate delivery of ionising radiation to the neoplastic tissue and pre-treatment in situ x-ray imaging plays an important role in meeting treatment accuracy requirements. Four dimensional cone-beam computed tomography (4D CBCT) is one such pre-treatment imaging technique that can help to visualise tumour target motion due to breathing at the time of radiation treatment delivery. Measuring and characterising the target motion can help to ensure highly accurate therapeutic x-ray beam delivery. In this thesis, a novel pre-treatment x-ray imaging technique, called Respiratory Triggered 4D cone-beam Computed Tomography (RT 4D CBCT), is conceived and investigated. Specifically, the aim of this work is to progress the 4D CBCT imaging technology by investigating the use of a patient’s breathing signal to improve and optimise the use of imaging radiation in 4D CBCT to facilitate the accurate delivery of radiation therapy. These investigations are presented in three main studies: 1. Introduction to the concept of respiratory triggered four dimensional conebeam computed tomography. 2. A simulation study exploring the behaviour of RT 4D CBCT using patientmeasured respiratory data. 3. The experimental realisation of RT 4D CBCT working in a real-time acquisitions setting. The major finding from this work is that RT 4D CBCT can provide target motion information with a 50% reduction in the x-ray imaging dose applied to the patient