Nonlinear tube-fitting for the analysis of anatomical and functional structures

We are concerned with the estimation of the exterior surface and interior summaries of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT) and to a white-matter tract image produced using diffusion tensor imaging (DTI).Comment: Published in at http://dx.doi.org/10.1214/10-AOAS384 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Colon centreline calculation for CT colonography using optimised 3D opological thinning

CT colonography is an emerging technique for colorectal cancer screening. This technique facilitates noninvasive imaging of the colon interior by generating virtual reality models of the colon lumen. Manual navigation through these models is a slow and tedious process. It is possible to automate navigation by calculating the centreline of the colon lumen. There are numerous well documented approaches for centreline calculation. Many of these techniques have been developed as alternatives to 3D topological thinning which has been discounted by others due to its computationally intensive nature. This paper describes a fully automated, optimised version of 3D topological thinning that has been specifically developed for calculating the centreline of the human colon

Fast colon centreline calculation using optimised 3D topological thinning

Topological thinning can be used to accurately identify the central path through a computer model of the colon generated using computed tomography colonography. The central path can subsequently be used to simplify the task of navigation within the colon model. Unfortunately standard topological thinning is an extremely inefficient process. We present an optimised version of topological thinning that significantly improves the performance of centreline calculation without compromising the accuracy of the result. This is achieved by using lookup tables to reduce the computational burden associated with the thinning process

NONLINEAR TUBE-FITTING FOR THE ANALYSIS OF ANATOMICAL AND FUNCTIONAL STRUCTURES

We are concerned with the estimation of the exterior surface of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT)and to a white-matter tract image produced using diffusion tensor `imaging (DTI)

A CASE STUDY IN PHARMACOLOGIC IMAGING USING PRINCIPAL CURVES IN SINGLE PHOTON EMISSION COMPUTED TOMOGRAPHY

In this manuscript we are concerned with functional imaging of the colon to assess the kinetics of a microbicide lubricant. The overarching goal is to understand the distribution of the lubricant in the colon. Such information is crucial for understanding the potential impact of the microbicide on HIV viral transmission. The experiment was conducted by imaging a radiolabeled lubricant distributed in the subject’s colon. The tracer imaging was conducted via single photon emission computed tomography (SPECT), a non-invasive, in-vivo functional imaging technique. We develop a novel principal curve algorithm to construct a three dimensional curve through the colon images. The developed algorithm is tested and debugged on several difficult two dimensional images of familiar curves where the original principal curve algorithm does not apply. The final curve fit to the colon data is compared with experimental sigmoidoscope collection

Surgical Subtask Automation for Intraluminal Procedures using Deep Reinforcement Learning

Intraluminal procedures have opened up a new sub-field of minimally invasive surgery that use flexible instruments to navigate through complex luminal structures of the body, resulting in reduced invasiveness and improved patient benefits. One of the major challenges in this field is the accurate and precise control of the instrument inside the human body. Robotics has emerged as a promising solution to this problem. However, to achieve successful robotic intraluminal interventions, the control of the instrument needs to be automated to a large extent. The thesis first examines the state-of-the-art in intraluminal surgical robotics and identifies the key challenges in this field, which include the need for safe and effective tool manipulation, and the ability to adapt to unexpected changes in the luminal environment. To address these challenges, the thesis proposes several levels of autonomy that enable the robotic system to perform individual subtasks autonomously, while still allowing the surgeon to retain overall control of the procedure. The approach facilitates the development of specialized algorithms such as Deep Reinforcement Learning (DRL) for subtasks like navigation and tissue manipulation to produce robust surgical gestures. Additionally, the thesis proposes a safety framework that provides formal guarantees to prevent risky actions. The presented approaches are evaluated through a series of experiments using simulation and robotic platforms. The experiments demonstrate that subtask automation can improve the accuracy and efficiency of tool positioning and tissue manipulation, while also reducing the cognitive load on the surgeon. The results of this research have the potential to improve the reliability and safety of intraluminal surgical interventions, ultimately leading to better outcomes for patients and surgeons

Developing advanced mathematical models for detecting abnormalities in 2D/3D medical structures.

Detecting abnormalities in two-dimensional (2D) and three-dimensional (3D) medical structures is among the most interesting and challenging research areas in the medical imaging field. Obtaining the desired accurate automated quantification of abnormalities in medical structures is still very challenging. This is due to a large and constantly growing number of different objects of interest and associated abnormalities, large variations of their appearances and shapes in images, different medical imaging modalities, and associated changes of signal homogeneity and noise for each object. The main objective of this dissertation is to address these problems and to provide proper mathematical models and techniques that are capable of analyzing low and high resolution medical data and providing an accurate, automated analysis of the abnormalities in medical structures in terms of their area/volume, shape, and associated abnormal functionality. This dissertation presents different preliminary mathematical models and techniques that are applied in three case studies: (i) detecting abnormal tissue in the left ventricle (LV) wall of the heart from delayed contrast-enhanced cardiac magnetic resonance images (MRI), (ii) detecting local cardiac diseases based on estimating the functional strain metric from cardiac cine MRI, and (iii) identifying the abnormalities in the corpus callosum (CC) brain structure—the largest fiber bundle that connects the two hemispheres in the brain—for subjects that suffer from developmental brain disorders. For detecting the abnormal tissue in the heart, a graph-cut mathematical optimization model with a cost function that accounts for the object’s visual appearance and shape is used to segment the the inner cavity. The model is further integrated with a geometric model (i.e., a fast marching level set model) to segment the outer border of the myocardial wall (the LV). Then the abnormal tissue in the myocardium wall (also called dead tissue, pathological tissue, or infarct area) is identified based on a joint Markov-Gibbs random field (MGRF) model of the image and its region (segmentation) map that accounts for the pixel intensities and the spatial interactions between the pixels. Experiments with real in-vivo data and comparative results with ground truth (identified by a radiologist) and other approaches showed that the proposed framework can accurately detect the pathological tissue and can provide useful metrics for radiologists and clinicians. To estimate the strain from cardiac cine MRI, a novel method based on tracking the LV wall geometry is proposed. To achieve this goal, a partial differential equation (PDE) method is applied to track the LV wall points by solving the Laplace equation between the LV contours of each two successive image frames over the cardiac cycle. The main advantage of the proposed tracking method over traditional texture-based methods is its ability to track the movement and rotation of the LV wall based on tracking the geometric features of the inner, mid-, and outer walls of the LV. This overcomes noise sources that come from scanner and heart motion. To identify the abnormalities in the CC from brain MRI, the CCs are aligned using a rigid registration model and are segmented using a shape-appearance model. Then, they are mapped to a simple unified space for analysis. This work introduces a novel cylindrical mapping model, which is conformal (i.e., one to one transformation and bijective), that enables accurate 3D shape analysis of the CC in the cylindrical domain. The framework can detect abnormalities in all divisions of the CC (i.e., splenium, rostrum, genu and body). In addition, it offers a whole 3D analysis of the CC abnormalities instead of only area-based analysis as done by previous groups. The initial classification results based on the centerline length and CC thickness suggest that the proposed CC shape analysis is a promising supplement to the current techniques for diagnosing dyslexia. The proposed techniques in this dissertation have been successfully tested on complex synthetic and MR images and can be used to advantage in many of today’s clinical applications of computer-assisted medical diagnostics and intervention

Automatic extraction of bronchus and centerline determination from CT images for three dimensional virtual bronchoscopy.

Law Tsui Ying.Thesis (M.Phil.)--Chinese University of Hong Kong, 2000.Includes bibliographical references (leaves 64-70).Abstracts in English and Chinese.Acknowledgments --- p.iiChapter 1 --- Introduction --- p.1Chapter 1.1 --- Structure of Bronchus --- p.3Chapter 1.2 --- Existing Systems --- p.4Chapter 1.2.1 --- Virtual Endoscope System (VES) --- p.4Chapter 1.2.2 --- Virtual Reality Surgical Simulator --- p.4Chapter 1.2.3 --- Automated Virtual Colonoscopy (AVC) --- p.5Chapter 1.2.4 --- QUICKSEE --- p.5Chapter 1.3 --- Organization of Thesis --- p.6Chapter 2 --- Three Dimensional Visualization in Medicine --- p.7Chapter 2.1 --- Acquisition --- p.8Chapter 2.1.1 --- Computed Tomography --- p.8Chapter 2.2 --- Resampling --- p.9Chapter 2.3 --- Segmentation and Classification --- p.9Chapter 2.3.1 --- Segmentation by Thresholding --- p.10Chapter 2.3.2 --- Segmentation by Texture Analysis --- p.10Chapter 2.3.3 --- Segmentation by Region Growing --- p.10Chapter 2.3.4 --- Segmentation by Edge Detection --- p.11Chapter 2.4 --- Rendering --- p.12Chapter 2.5 --- Display --- p.13Chapter 2.6 --- Hazards of Visualization --- p.13Chapter 2.6.1 --- Adding Visual Richness and Obscuring Important Detail --- p.14Chapter 2.6.2 --- Enhancing Details Incorrectly --- p.14Chapter 2.6.3 --- The Picture is not the Patient --- p.14Chapter 2.6.4 --- Pictures-'R'-Us --- p.14Chapter 3 --- Overview of Advanced Segmentation Methodologies --- p.15Chapter 3.1 --- Mathematical Morphology --- p.15Chapter 3.2 --- Recursive Region Search --- p.16Chapter 3.3 --- Active Region Models --- p.17Chapter 4 --- Overview of Centerline Methodologies --- p.18Chapter 4.1 --- Thinning Approach --- p.18Chapter 4.2 --- Volume Growing Approach --- p.21Chapter 4.3 --- Combination of Mathematical Morphology and Region Growing Schemes --- p.22Chapter 4.4 --- Simultaneous Borders Identification Approach --- p.23Chapter 4.5 --- Tracking Approach --- p.24Chapter 4.6 --- Distance Transform Approach --- p.25Chapter 5 --- Automated Extraction of Bronchus Area --- p.27Chapter 5.1 --- Basic Idea --- p.27Chapter 5.2 --- Outline of the Automated Extraction Algorithm --- p.28Chapter 5.2.1 --- Selection of a Start Point --- p.28Chapter 5.2.2 --- Three Dimensional Region Growing Method --- p.29Chapter 5.2.3 --- Optimization of the Threshold Value --- p.29Chapter 5.3 --- Retrieval of Start Point Algorithm Using Genetic Algorithm --- p.29Chapter 5.3.1 --- Introduction to Genetic Algorithm --- p.30Chapter 5.3.2 --- Problem Modeling --- p.31Chapter 5.3.3 --- Algorithm for Determining a Start Point --- p.33Chapter 5.3.4 --- Genetic Operators --- p.33Chapter 5.4 --- Three Dimensional Painting Algorithm --- p.34Chapter 5.4.1 --- Outline of the Three Dimensional Painting Algorithm --- p.34Chapter 5.5 --- Optimization of the Threshold Value --- p.36Chapter 6 --- Automatic Centerline Determination Algorithm --- p.38Chapter 6.1 --- Distance Transformations --- p.38Chapter 6.2 --- End Points Retrieval --- p.41Chapter 6.3 --- Graph Based Centerline Algorithm --- p.44Chapter 7 --- Experiments and Discussion --- p.48Chapter 7.1 --- Experiment of Automated Determination of Bronchus Algorithm --- p.48Chapter 7.2 --- Experiment of Automatic Centerline Determination Algorithm --- p.54Chapter 8 --- Conclusion --- p.62Bibliography --- p.6