The use of the Kalman filter in the automated segmentation of EIT lung images

In this paper, we present a new pipeline for the fast and accurate segmentation of impedance images of the lungs using electrical impedance tomography (EIT). EIT is an emerging, promising, non-invasive imaging modality that produces real-time, low spatial but high temporal resolution images of impedance inside a body. Recovering impedance itself constitutes a nonlinear ill-posed inverse problem, therefore the problem is usually linearized, which produces impedance-change images, rather than static impedance ones. Such images are highly blurry and fuzzy along object boundaries. We provide a mathematical reasoning behind the high suitability of the Kalman filter when it comes to segmenting and tracking conductivity changes in EIT lung images. Next, we use a two-fold approach to tackle the segmentation problem. First, we construct a global lung shape to restrict the search region of the Kalman filter. Next, we proceed with augmenting the Kalman filter by incorporating an adaptive foreground detection system to provide the boundary contours for the Kalman filter to carry out the tracking of the conductivity changes as the lungs undergo deformation in a respiratory cycle. The proposed method has been validated by using performance statistics such as misclassified area, and false positive rate, and compared to previous approaches. The results show that the proposed automated method can be a fast and reliable segmentation tool for EIT imaging

Automatic Bootstrapping and Tracking of Object Contours

This work introduces a new fully automatic object tracking and segmentation framework. The framework consists of a motion based bootstrapping algorithm concurrent to a shape based active contour. The shape based active contour uses a finite shape memory that is automatically and continuously built from both the bootstrap process and the active contour object tracker. A scheme is proposed to ensure the finite shape memory is continuously updated but forgets unnecessary information. Two new ways of automatically extracting shape information from image data given a region of interest are also proposed. Results demonstrate that the bootstrapping stage provides important motion and shape information to the object tracker

A method and software for segmentation of anatomic object ensembles by deformable m-reps: Deformable M-Reps

Deformable shape models (DSMs) comprise a general approach that shows great promise for automatic image segmentation. Published studies by others and our own research results strongly suggest that segmentation of a normal or near-normal object from 3D medical images will be most successful when the DSM approach uses 1) knowledge of the geometry of not only the target anatomic object but also the ensemble of objects providing context for the target object and 2) knowledge of the image intensities to be expected relative to the geometry of the target and contextual objects. The segmentation will be most efficient when the deformation operates at multiple object-related scales and uses deformations that include not just local translations but the biologically important transformations of bending and twisting, i.e., local rotation, and local magnification. In computer vision an important class of DSM methods uses explicit geometric models in a Bayesian statistical framework to provide a priori information used in posterior optimization to match the DSM against a target image. In this approach a DSM of the object to be segmented is placed in the target image data and undergoes a series of rigid and non-rigid transformations that deform the model to closely match the target object. The deformation process is driven by optimizing an objective function that has terms for the geometric typicality and model-to-image match for each instance of the deformed model. The success of this approach depends strongly on the object representation, i.e., the structural details and parameter set for the DSM, which in turn determines the analytic form of the objective function. This paper describes a form of DSM called m-reps that has or allows these properties, and a method of segmentation consisting of large to small scale posterior optimization of m-reps. Segmentation by deformable m-reps, together with the appropriate data representations, visualizations, and user interface, has been implemented in software that accomplishes 3D segmentations in a few minutes. Software for building and training models has also been developed. The methods underlying this software and its abilities are the subject of this paper

Automated voxel-wise brain DTI analysis of fitness and aging

Diffusion Tensor Imaging (DTI) has become a widely used MR modality to investigate white matter integrity in the brain. This paper presents the application of an automated method for voxel-wise group comparisons of DTI images in a study of fitness and aging. The automated processing method consists of 3 steps: 1) preprocessing including image format converting, image quality control, eddy-current and motion artifact correction, skull stripping and tensor image estimation, 2) study-specific unbiased DTI atlas computation via diffeomorphic fluid-based and demons deformable registration and 3) voxel-wise statistical analysis via heterogeneous linear regression and a wild bootstrap technique for correcting for multiple comparisons. Our results show that this fully automated method is suitable for voxel-wise group DTI analysis. Furthermore, in older adults, the results suggest a strong link between reduced fractional anisotropy (FA) values, fitness and aging

Three Dimensional Nonlinear Statistical Modeling Framework for Morphological Analysis

This dissertation describes a novel three-dimensional (3D) morphometric analysis framework for building statistical shape models and identifying shape differences between populations. This research generalizes the use of anatomical atlases on more complex anatomy as in case of irregular, flat bones, and bones with deformity and irregular bone growth. The foundations for this framework are: 1) Anatomical atlases which allow the creation of homologues anatomical models across populations; 2) Statistical representation for output models in a compact form to capture both local and global shape variation across populations; 3) Shape Analysis using automated 3D landmarking and surface matching. The proposed framework has various applications in clinical, forensic and physical anthropology fields. Extensive research has been published in peer-reviewed image processing, forensic anthropology, physical anthropology, biomedical engineering, and clinical orthopedics conferences and journals. The forthcoming discussion of existing methods for morphometric analysis, including manual and semi-automatic methods, addresses the need for automation of morphometric analysis and statistical atlases. Explanations of these existing methods for the construction of statistical shape models, including benefits and limitations of each method, provide evidence of the necessity for such a novel algorithm. A novel approach was taken to achieve accurate point correspondence in case of irregular and deformed anatomy. This was achieved using a scale space approach to detect prominent scale invariant features. These features were then matched and registered using a novel multi-scale method, utilizing both coordinate data as well as shape descriptors, followed by an overall surface deformation using a new constrained free-form deformation. Applications of output statistical atlases are discussed, including forensic applications for the skull sexing, as well as physical anthropology applications, such as asymmetry in clavicles. Clinical applications in pelvis reconstruction and studying of lumbar kinematics and studying thickness of bone and soft tissue are also discussed

Analysis and Manipulation of Repetitive Structures of Varying Shape

Self-similarity and repetitions are ubiquitous in man-made and natural objects. Such structural regularities often relate to form, function, aesthetics, and design considerations. Discovering structural redundancies along with their dominant variations from 3D geometry not only allows us to better understand the underlying objects, but is also beneficial for several geometry processing tasks including compact representation, shape completion, and intuitive shape manipulation. To identify these repetitions, we present a novel detection algorithm based on analyzing a graph of surface features. We combine general feature detection schemes with a RANSAC-based randomized subgraph searching algorithm in order to reliably detect recurring patterns of locally unique structures. A subsequent segmentation step based on a simultaneous region growing is applied to verify that the actual data supports the patterns detected in the feature graphs. We introduce our graph based detection algorithm on the example of rigid repetitive structure detection. Then we extend the approach to allow more general deformations between the detected parts. We introduce subspace symmetries whereby we characterize similarity by requiring the set of repeating structures to form a low dimensional shape space. We discover these structures based on detecting linearly correlated correspondences among graphs of invariant features. The found symmetries along with the modeled variations are useful for a variety of applications including non-local and non-rigid denoising. Employing subspace symmetries for shape editing, we introduce a morphable part model for smart shape manipulation. The input geometry is converted to an assembly of deformable parts with appropriate boundary conditions. Our method uses self-similarities from a single model or corresponding parts of shape collections as training input and allows the user also to reassemble the identified parts in new configurations, thus exploiting both the discrete and continuous learned variations while ensuring appropriate boundary conditions across part boundaries. We obtain an interactive yet intuitive shape deformation framework producing realistic deformations on classes of objects that are difficult to edit using repetition-unaware deformation techniques

Taking shape: The data science of elastic shape analysis with practical applications

This thesis was submitted for the award of Doctor of Philosophy and was awarded by Brunel University London.A mathematical curve can represent many different objects, both physical and abstract, from the outline curve of an artefact in an image to the weight of growing animal to the set of frequencies used in a sound. Regardless of these variations, the curves can almost always vary non-linearly. One way to study shapes and their potential variations is elastic shape analysis, a rich theory of which has developed over the past twenty years. However, methods of elastic shape analysis are seldom utilized in practical applications on real-world data, especially outside of the mathematical shape analysis community. Our aim in this thesis is to explore some practical applications of elastic shape analysis. To do this, we work with various types of shape data, the majority of which are based on image datasets. As our focus is on two-dimensional curves, it is important to be able to robustly extract contours from images, before we can apply elastic shape analysis tools. In order to analyse the shapes in a dataset, we turn to methods of machine learning, to investigate the applications of elastic shape analysis in classification. In this thesis, we introduce an anthology of projects, in order to emphasise and under- stand the potential of elastic shape analysis in practical applications. There are four main projects in this thesis: (i) Classification of objects using outlines and the comparisons between methods of elastic shape analysis, geometric morphometrics, and human experts, with a focus on ancient Greek vases, (ii) Mussel species identification and a demonstra- tion that shape may not be enough in some applications, (iii) A novel tool to monitor the development of k ̄ak ̄ap ̄o chicks, and (iv) Classifying individual kiwi based on acoustic data from their calls. By combining tools from computer vision and machine learning with methods of elastic shape analysis, we introduce a practical framework for the application of elastic shape analysis, through a data science lens

Doctor of Philosophy in Computing

dissertationStatistical shape analysis has emerged as an important tool for the quantitative analysis of anatomy in many medical imaging applications. The correspondence based approach to evaluate shape variability is a popular method, based on comparing configurations of carefully placed landmarks on each shape. In recent years, methods for automatic placement of landmarks have enhanced the ability of this approach to capture statistical properties of shape populations. However, biomedical shapes continue to present considerable difficulties in automatic correspondence optimization due to inherent geometric complexity and the need to correlate shape change with underlying biological parameters. This dissertation addresses these technical difficulties and presents improved shape correspondence models. In particular, this dissertation builds on the particle-based modeling (PBM) framework described by Joshua Cates' 2010 Ph.D. dissertation. In the PBM framework, correspondences are modeled as a set of dynamic points or a particle system, positioned automatically on shape surfaces by optimizing entropy contained in the model, with the idea of balancing model simplicity against accuracy of the particle system representation of shapes. This dissertation is a collection of four papers that extend the PBM framework to include shape regression and longitudinal analysis and also adds new methods to improve modeling of complex shapes. It also includes a summary of two applications from the field of orthopaedics. Technical details of the PBM framework are provided in Chapter 2, after which the first topic related to the study of shape change over time is addressed (Chapters 3 and 4). In analyses of normative growth or disease progression, shape regression models allow characterization of the underlying biological process while also facilitating comparison of a sample against a normative model. The first paper introduces a shape regression model into the PBM framework to characterize shape variability due to an underlying biological parameter. It further confirms the statistical significance of this relationship via systematic permutation testing. Simple regression models are, however, not sufficient to leverage information provided by longitudinal studies. Longitudinal studies collect data at multiple time points for each participant and have the potential to provide a rich picture of the anatomical changes occurring during development, disease progression, or recovery. The second paper presents a linear-mixed-effects (LME) shape model in order to fully leverage the high-dimensional, complex features provided by longitudinal data. The parameters of the LME shape model are estimated in a hierarchical manner within the PBM framework. The topic of geometric complexity present in certain biological shapes is addressed next (Chapters 5 and 6). Certain biological shapes are inherently complex and highly variable, inhibiting correspondence based methods from producing a faithful representation of the average shape. In the PBM framework, use of Euclidean distances leads to incorrect particle system interactions while a position-only representation leads to incorrect correspondences around sharp features across shapes. The third paper extends the PBM framework to use efficiently computed geodesic distances and also adds an entropy term based on the surface normal. The fourth paper further replaces the position-only representation with a more robust distance-from-landmark feature in the PBM framework to obtain isometry invariant correspondences. Finally, the above methods are applied to two applications from the field of orthopaedics. The first application uses correspondences across an ensemble of human femurs to characterize morphological shape differences due to femoroacetabular impingement. The second application involves an investigation of the short bone phenotype apparent in mouse models of multiple osteochondromas. Metaphyseal volume deviations are correlated with deviations in length to quantify the effect of cancer toward the apparent shortening of long bones (femur, tibia-fibula) in mouse models