Image databases in medical applications

The number of medical images acquired yearly in hospitals increases all the time. These imaging data contain lots of information on the characteristics of anatomical structures and on their variations. This information can be utilized in numerous medical applications. In deformable model-based segmentation and registration methods, the information in the image databases can be used to give a priori information on the shape of the object studied and the gray-level values in the image, and on their variations. On the other hand, by studying the variations of the object of interest in different populations, the effects of, for example, aging, gender, and diseases on anatomical structures can be detected. In the work described in this Thesis, methods that utilize image databases in medical applications were studied. Methods were developed and compared for deformable model-based segmentation and registration. Model selection procedure, mean models, and combination of classifiers were studied for the construction of a good a priori model. Statistical and probabilistic shape models were generated to constrain the deformations in segmentation and registration so that only the shapes typical to the object studied were accepted. In the shape analysis of the striatum, both volume and local shape changes were studied. The effects of aging and gender, and also the asymmetries were examined. The results proved that the segmentation and registration accuracy of deformable model-based methods can be improved by utilizing the information in image databases. The databases used were relatively small. Therefore, the statistical and probabilistic methods were not able to model all the population-specific variation. On the other hand, the simpler methods, the model selection procedure, mean models, and combination of classifiers, gave good results also with the small image databases. Two main applications were the reconstruction of 3-D geometry from incomplete data and the segmentation of heart ventricles and atria from short- and long-axis magnetic resonance images. In both applications, the methods studied provided promising results. The shape analysis of the striatum showed that the volume of the striatum decreases in aging. Also, the shape of the striatum changes locally. Asymmetries in the shape were found, too, but any gender-related local shape differences were not found.reviewe

Statistics of populations of images and its embedded objects: driving applications in neuroimaging

journal articleWork in progress towards modeling shape statistics of multi-object complexes is presented. Constraints defined by the set of objects such as a compact representation of object shape relationships and correlation of shape changes might have advantages for automatic segmentation and group discrimination. We present a concept for statistical multi-object modeling and discuss the major challenges which are a reduction to a small set of descriptive features, calculation of mean and variability via curved statistics, the choice of aligning sets of multiple objects, and the problem of describing the statistics of object pose and object shape and their interrelationship. Shape modeling and analysis is demonstrated with an application to a longitudinal autism study, with shape modeling of sets of 10 subcortical structures in a population of 20 subjects

Development of an Atlas-Based Segmentation of Cranial Nerves Using Shape-Aware Discrete Deformable Models for Neurosurgical Planning and Simulation

Twelve pairs of cranial nerves arise from the brain or brainstem and control our sensory functions such as vision, hearing, smell and taste as well as several motor functions to the head and neck including facial expressions and eye movement. Often, these cranial nerves are difficult to detect in MRI data, and thus represent problems in neurosurgery planning and simulation, due to their thin anatomical structure, in the face of low imaging resolution as well as image artifacts. As a result, they may be at risk in neurosurgical procedures around the skull base, which might have dire consequences such as the loss of eyesight or hearing and facial paralysis. Consequently, it is of great importance to clearly delineate cranial nerves in medical images for avoidance in the planning of neurosurgical procedures and for targeting in the treatment of cranial nerve disorders. In this research, we propose to develop a digital atlas methodology that will be used to segment the cranial nerves from patient image data. The atlas will be created from high-resolution MRI data based on a discrete deformable contour model called 1-Simplex mesh. Each of the cranial nerves will be modeled using its centerline and radius information where the centerline is estimated in a semi-automatic approach by finding a shortest path between two user-defined end points. The cranial nerve atlas is then made more robust by integrating a Statistical Shape Model so that the atlas can identify and segment nerves from images characterized by artifacts or low resolution. To the best of our knowledge, no such digital atlas methodology exists for segmenting nerves cranial nerves from MRI data. Therefore, our proposed system has important benefits to the neurosurgical community

Fitting Skeletal Object Models Using Spherical Harmonics Based Template Warping

We present a scheme that propagates a reference skeletal model (s-rep) into a particular case of an object, thereby propagating the initial shape-related layout of the skeleton-to-boundary vectors, called spokes. The scheme represents the surfaces of the template as well as the target objects by spherical harmonics and computes a warp between these via a thin plate spline. To form the propagated s-rep, it applies the warp to the spokes of the template s-rep and then statistically refines. This automatic approach promises to make s-rep fitting robust for complicated objects, which allows s-rep based statistics to be available to all. The improvement in fitting and statistics is significant compared with the previous methods and in statistics compared with a state-of-the-art boundary based method

Smart Cage Active Contours and their application to brain image segmentation

In this work we present a new segmentation method named Smart Cage Active Contours (SCAC) that combines a parametrized active contour framework named Cage Active Contours (CAC), based on a ne trans- formations, with Active Shape Models (ASM). Our method e ectively restricts the shapes the evolving contours can take without the need of the training images to be manually landmarked. We apply our method to segment the caudate nuclei subcortical structure of a set of 40 subjects in magnetic resonance brain images, with promising results

Computational Anatomy for Multi-Organ Analysis in Medical Imaging: A Review

The medical image analysis field has traditionally been focused on the development of organ-, and disease-specific methods. Recently, the interest in the development of more 20 comprehensive computational anatomical models has grown, leading to the creation of multi-organ models. Multi-organ approaches, unlike traditional organ-specific strategies, incorporate inter-organ relations into the model, thus leading to a more accurate representation of the complex human anatomy. Inter-organ relations are not only spatial, but also functional and physiological. Over the years, the strategies 25 proposed to efficiently model multi-organ structures have evolved from the simple global modeling, to more sophisticated approaches such as sequential, hierarchical, or machine learning-based models. In this paper, we present a review of the state of the art on multi-organ analysis and associated computation anatomy methodology. The manuscript follows a methodology-based classification of the different techniques 30 available for the analysis of multi-organs and multi-anatomical structures, from techniques using point distribution models to the most recent deep learning-based approaches. With more than 300 papers included in this review, we reflect on the trends and challenges of the field of computational anatomy, the particularities of each anatomical region, and the potential of multi-organ analysis to increase the impact of 35 medical imaging applications on the future of healthcare.Comment: Paper under revie

Construction of boundary element models in bioelectromagnetism

Multisensor electro- and magnetoencephalographic (EEG and MEG) as well as electro- and magnetocardiographic (ECG and MCG) recordings have been proved useful in noninvasively extracting information on bioelectric excitation. The anatomy of the patient needs to be taken into account, when excitation sites are localized by solving the inverse problem. In this work, a methodology has been developed to construct patient specific boundary element models for bioelectromagnetic inverse problems from magnetic resonance (MR) data volumes as well as from two orthogonal X-ray projections. The process consists of three main steps: reconstruction of 3-D geometry, triangulation of reconstructed geometry, and registration of the model with a bioelectromagnetic measurement system. The 3-D geometry is reconstructed from MR data by matching a 3-D deformable boundary element template to images. The deformation is accomplished as an energy minimization process consisting of image and model based terms. The robustness of the matching is improved by multi-resolution and global-to-local approaches as well as using oriented distance maps. A boundary element template is also used when 3-D geometry is reconstructed from X-ray projections. The deformation is first accomplished in 2-D for the contours of simulated, built from the template, and real X-ray projections. The produced 2-D vector field is back-projected and interpolated on the 3-D template surface. A marching cube triangulation is computed for the reconstructed 3-D geometry. Thereafter, a non-iterative mesh-simplification method is applied. The method is based on the Voronoi-Delaunay duality on a 3-D surface with discrete distance measures. Finally, the triangulated surfaces are registered with a bioelectromagnetic measurement utilizing markers. More than fifty boundary element models have been successfully constructed from MR images using the methods developed in this work. A simulation demonstrated the feasibility of X-ray reconstruction; some practical problems of X-ray imaging need to be solved to begin tests with real data.reviewe

Statistical Shape Modelling and Segmentation of the Respiratory Airway

The human respiratory airway consists of the upper (nasal cavity, pharynx) and the lower (trachea, bronchi) respiratory tracts. Accurate segmentation of these two airway tracts can lead to better diagnosis and interpretation of airway-specific diseases, and lead to improvement in the localization of abnormal metabolic or pathological sites found within and/or surrounding the respiratory regions. Due to the complexity and the variability displayed in the anatomical structure of the upper respiratory airway along with the challenges in distinguishing the nasal cavity from non-respiratory regions such as the paranasal sinuses, it is difficult for existing algorithms to accurately segment the upper airway without manual intervention. This thesis presents an implicit non-parametric framework for constructing a statistical shape model (SSM) of the upper and lower respiratory tract, capable of distinct shape generation and be adapted for segmentation. An SSM of the nasal cavity was successfully constructed using 50 nasal CT scans. The performance of the SSM was evaluated for compactness, specificity and generality. An averaged distance error of 1.47 mm was measured for the generality assessment. The constructed SSM was further adapted with a modified locally constrained random walk algorithm to segment the nasal cavity. The proposed algorithm was evaluated on 30 CT images and outperformed comparative state-of-the-art and conventional algorithms. For the lower airway, a separate algorithm was proposed to automatically segment the trachea and bronchi, and was designed to tolerate the image characteristics inherent in low-contrast CT images. The algorithm was evaluated on 20 clinical low-contrast CT from PET-CT patient studies and demonstrated better performance (87.1±2.8 DSC and distance error of 0.37±0.08 mm) in segmentation results against comparative state-of-the-art algorithms

A Heuristic Image Search Algorithm for Active Shape Model Segmentation of the Caudate Nucleus and Hippocampus in Brain MR Images of Children with FASD

Magnetic Resonance Imaging provides a non-invasive means to study the neural correlates of Fetal Alcohol Spectrum Disorder (FASD) - the most common form of preventable mental retardation worldwide. One approach aims to detect brain abnormalities through an assessment of volume and shape of two sub-cortical structures, the caudate nucleus and hippocampus. We present a method for automatically segmenting these structures from high-resolution MR images captured as part of an ongoing study into the neural correlates of FASD. Our method incorporates an Active Shape Model, which is used to learn shape variation from manually segmented training data. A modified discrete Geometrically Deformable Model is used to generate point correspondence between training models. An ASM is then created from the landmark points. Experiments were conducted on the image search phase of ASM segmentation, in order to find the technique best suited to segmentation of the hippocampus and caudate nucleus. Various popular image search techniques were tested, including an edge detection method and a method based on grey profile Mahalanobis distance measurement. A novel heuristic image search method was also developed and tested. This heuristic method improves image segmentation by taking advantage of characteristics specific to the target data, such as a relatively homogeneous tissue colour in target structures. Results show that ASMs that use the heuristic image search technique produce the most accurate segmentations. An ASM constructed using this technique will enable researchers to quickly, reliably, and automatically segment test data for use in the FASD study