An adaptive-focus statistical shape model for segmentation and shape modeling of 3-D brain structures

Targeted prostate biopsy using statistical image analysis

In this paper, a method for maximizing the probability of prostate cancer detection via biopsy is presented, by combining image analysis and optimization techniques. This method consists of three major steps. First, a statistical atlas of the spatial distribution of prostate cancer is constructed from histological images obtained from radical prostatectomy specimen. Second, a probabilistic optimization framework is employed to optimize the biopsy strategy, so that the probability of cancer detection is maximized under needle placement uncertainties. Finally, the optimized biopsy strategy generated in the atlas space is mapped to a specific patient space using an automated segmentation and elastic registration method. Cross-validation experiments showed that the predictive power of the optimized biopsy strategy for cancer detection reached the 94%-96% levels for 6-7 biopsy cores, which is significantly better than standard random-systematic biopsy protocols, thereby encouraging further investigation of optimized biopsy strategies in prospective clinical studies.

Registering histological and MR images of prostate for image-based cancer detection

This paper presents a deformable registration method to co-register histological images with MR images of the same prostate. By considering various distortion and cutting artifacts in histological images and also fundamentally different nature of histological and MR images, our registration method is thus guided by two types of landmark points that can be reliably detected in both histological and MR images, i.e., prostate boundary points, and internal salient points that can be identified by a scale-space analysis method. The similarity between these automatically detected landmarks in histological and MR images are defined by geometric features and normalized mutual information, respectively. By optimizing a function, which integrates the similarities between landmarks with spatial constraints, the correspondences between the landmarks as well as the deformable transformation between histological and MR images can be simultaneously obtained. The performance of our proposed registration algorithm has been evaluated by various designed experiments. This work is part of a larger effort to develop statistical atlases of prostate cancer using both imaging and histological information, and to use these atlases for optimal biopsy and therapy planning.

Brain functional localization: a survey of image registration techniques

Functional localization is a concept which involves the application of a sequence of geometrical and statistical image processing operations in order to define the location of brain activity or to produce functional/parametric maps with respect to the brain structure or anatomy. Considering that functional brain images do not normally convey detailed structural information and, thus, do not present an anatomically specific localization of functional activity, various image registration techniques are introduced in the literature for the purpose of mapping functional activity into an anatomical image or a brain atlas. The problems addressed by these techniques differ depending on the application and the type of analysis, i.e., single-subject versus group analysis. Functional to anatomical brain image registration is the core part of functional localization in most applications and is accompanied by intersubject and subject-to-atlas registration for group analysis studies. Cortical surface registration and automatic brain labeling are some of the other tools towards establishing a fully automatic functional localization procedure. While several previous survey papers have reviewed and classified general-purpose medical image registration techniques, this paper provides an overview of brain functional localization along with a survey and classification of the image registration techniques related to this problem.

Comparing MR image intensity standardization against tissue characterizability of magnetization transfer ratio imaging

PURPOSE

To evaluate existing methods of standardization by exploiting the well-known tissue characterizing property of magnetization transfer ratio (MTR) values obtained from MT imaging, and compare the tissue characterizability of standardized T2, proton density (PD), and T1 images against the MTR images.

MATERIALS AND METHODS

Image intensity standardization is a postprocessing method that was designed to correct for acquisition-to-acquisition signal intensity variations (nonstandardness) inherent in magnetic resonance (MR) images. The main idea of this technique is to deform the volume image histogram of each study to match a standard histogram, and to utilize the resulting transformations to map the image intensities into a standard scale. The method has been shown to produce a significant gain in similarity of resulting images and to achieve numeric tissue characterization. In this work we compared PD-, T2-, and T1-weighted images before and after standardization with the corresponding MT images for 10 patient MRI studies of the brain, in terms of the normalized median values on the corresponding image histograms.

RESULTS

No statistically significant difference was observed between the standardized PD-, T2-, and T1-weighted images and the corresponding MTR images. However, a statistically significant difference was found between the pre- and poststandardized PD-, T2-, and T1-weighted images, and between the prestandardized PD-, T2-, and T1-weighted images and the corresponding MTR images.

CONCLUSION

These results suggest that standardized T2, PD, and T1 images and their tissue-specific intensity signatures may be useful for characterizing disease.

Localized maximum entropy shape modelling

A core part of many medical image segmentation techniques is the point distribution model, i.e., the landmark-based statistical shape model which describes the type of shapes under consideration. To build a proper model, that is flexible and generalizes well, one typically needs a large amount of landmarked training data, which can be hard to obtain. This problem becomes worse with increasing shape complexity and dimensionality. This work presents a novel methodology applicable to principal component-based shape model building and similar techniques. The main idea of the method is to make regular PCA shape modelling more flexible by using merely covariances between neighboring landmarks. The remaining unknown second order moments are determined using the maximum entropy principle based on which the full covariance matrix--as employed in the PCA--is determined using matrix completion. The method presented can be applied in a variety of situations and in conjunction with other technique facilitating model building. The experiments on point distributions demonstrate that improved shape models can be obtained using this localized maximum entropy modelling.

Computer technology in detection and staging of prostate carcinoma: A review

After two decades of increasing interest and research activity, computer-assisted diagnostic approaches are reaching the stage where more routine deployment in clinical practice is becoming a possibility [Kruppinski, E.A., 2004. Computer-aided detection in clinical environment: Benefits and challenges for radiologists. Radiology 231, 7-9]. This is particularly the case in the analysis of mammographic images [Helvie, M.A., Hadjiiski, L., Makariou, E., Chan, H.P., Petrick, N., Sahiner, B., Lo, S.C., Freedman, M., Adler, D., Bailey, J., Blane, C., Hoff, D., Hunt, K., Joynt, L., Klein, K., Paramagul, C., Patterson, S.K., Roubidoux, M.A., 2004. Sensitivity of noncommercial computer-aided detection system for mammographic breast cancer detection: pilot clinical trial. Radiology 231, 208-214] and in the detection of pulmonary nodules [Reeves, A.P., Kostis, W.J., 2000. Computer-aided diagnosis for lung cancer. Radiol. Clin. North Am. 38, 497-509]. However, similar approaches can be applied more widely with the promise of increasing clinical utility in other areas. We review how computer-aided approaches may be applied in the diagnosis and staging of prostatic cancer. The current status of computer technology is reviewed, covering artificial neural networks for detection and staging, computerised biopsy simulation and computer-assisted analysis of ultrasound and magnetic resonance images.

Deformable segmentation of 3-D ultrasound prostate images using statistical texture matching method

This paper presents a novel deformable model for automatic segmentation of prostates from three-dimensional ultrasound images, by statistical matching of both shape and texture. A set of Gabor-support vector machines (G-SVMs) are positioned on different patches of the model surface, and trained to adaptively capture texture priors of ultrasound images for differentiation of prostate and nonprostate tissues in different zones around prostate boundary. Each G-SVM consists of a Gabor filter bank for extraction of rotation-invariant texture features and a kernel support vector machine for robust differentiation of textures. In the deformable segmentation procedure, these pretrained G-SVMs are used to tentatively label voxels around the surface of deformable model as prostate or nonprostate tissues by a statistical texture matching. Subsequently, the surface of deformable model is driven to the boundary between the tentatively labeled prostate and non-prostate tissues. Since the step of tissue labeling and the step of label-based surface deformation are dependent on each other, these two steps are repeated until they converge. Experimental results by using both synthesized and real data show the good performance of the proposed model in segmenting prostates from ultrasound images.

Segmenting magnetic resonance images via hierarchical mixture modelling

We present a statistically innovative as well as scientifically and practically relevant method for automatically segmenting magnetic resonance images using hierarchical mixture models. Our method is a general tool for automated cortical analysis which promises to contribute substantially to the science of neuropsychiatry. We demonstrate that our method has advantages over competing approaches on a magnetic resonance brain imagery segmentation task.

A coarse-to-fine approach to prostate boundary segmentation in ultrasound images

Background

In this paper a novel method for prostate segmentation in transrectal ultrasound images is presented.

Methods

A segmentation procedure consisting of four main stages is proposed. In the first stage, a locally adaptive contrast enhancement method is used to generate a well-contrasted image. In the second stage, this enhanced image is thresholded to extract an area containing the prostate (or large portions of it). Morphological operators are then applied to obtain a point inside of this area. Afterwards, a Kalman estimator is employed to distinguish the boundary from irrelevant parts (usually caused by shadow) and generate a coarsely segmented version of the prostate. In the third stage, dilation and erosion operators are applied to extract outer and inner boundaries from the coarsely estimated version. Consequently, fuzzy membership functions describing regional and gray-level information are employed to selectively enhance the contrast within the prostate region. In the last stage, the prostate boundary is extracted using strong edges obtained from selectively enhanced image and information from the vicinity of the coarse estimation.

Results

A total average similarity of 98.76%(± 0.68) with gold standards was achieved.

Conclusion

The proposed approach represents a robust and accurate approach to prostate segmentation.

Asymmetry analysis of deformable hippocampal model using the principal component in schizophrenia

The hippocampus is thought to play an important role in learning and memory processing, and impairments in memory, attention, and decision making are found commonly in schizophrenia. Although many studies have reported decreases in hippocampal volume in the left hemisphere in schizophrenia, regionally specific hippocampal volume loss has not been revealed consistently using volume analysis. Recently, many studies have analyzed shape asymmetry using 3-D models; however, inconsistent results have been reported, mainly due to methodologic differences. We therefore used an active, flexible, deformable shape model for surface parameterization, and compared shape asymmetry based on principal component analysis (PCA) in the hippocampi of schizophrenic patients with those of the normal controls. Although the overall pattern of the statistical results did not change according to the number of principal components, the reconstructed results based on six major components were much more distinguishable. Although the left hemispheric hippocampal volume was larger than the right hemispheric was in this study, the difference was not significant. In shape asymmetry analysis, the right hemisphere hippocampus was bilaterally larger than the left hemisphere hippocampus was in the head of the superior CA1 and smaller in the tail and head of the inferior CA1. The asymmetry in the schizophrenia group was statistically smaller than that in the control group through reduction of the left hemisphere hippocampus volume.

Automatic extraction of brain surface and mid-sagittal plane from PET images applying deformable models

In this study, we propose and evaluate new methods for automatic extraction of the brain surface and the mid-sagittal plane from functional positron emission tomography (PET) images. Designing methods for these segmentation tasks is challenging because the spatial distribution of intensity values in a PET image depends on the applied radiopharmaceutical and the contrast to noise ratio in a PET image is typically low. We extracted the brain surface with a deformable model which is based on a global optimization algorithm. The global optimization allows reliable automation of the extraction task. Based on the extracted brain surface, the mid-sagittal plane was determined. The method was tested with the image of the Hoffman brain phantom (FDG) and the images from the brain studies with the FDG (17 images) and the C11-Raclopride tracers (4 images). In addition to the brain surfaces, we applied the deformable model for extraction of the coarse cortical structure based on the tracer uptake from FDG-PET brain images. The proposed segmentation methods provide a promising direction for automatic processing and analysis of PET brain images.

Interplay between intensity standardization and inhomogeneity correction in MR image processing

Image intensity standardization is a postprocessing method designed for correcting acquisition-to-acquisition signal intensity variations (nonstandardness) inherent in magnetic resonance (MR) images. Inhomogeneity correction is a process used to suppress the low frequency background nonuniformities (inhomogeneities) of the image domain that exist in MR images. Both these procedures have important implications in MR image analysis. The effects of these postprocessing operations on improvement of image quality in isolation has been well documented. However, the combined effects of these two processes on MR images and how the processes influence each other have not been studied thus far. In this paper, we evaluate the effect of inhomogeneity correction followed by standardization and vice-versa on MR images in order to determine the best sequence to follow for enhancing image quality. We conducted experiments on several clinical and phantom data sets (nearly 4000 three-dimensional MR images were analyzed) corresponding to four different MRI protocols. Different levels of artificial nonstandardness, and different models and levels of artificial background inhomogeneity were used in these experiments. Our results indicate that improved standardization can be achieved by preceding it with inhomogeneity correction. There is no statistically significant difference in image quality obtained between the results of standardization followed by correction and that of correction followed by standardization from the perspective of inhomogeneity correction. The correction operation is found to bias the effect of standardization. We demonstrate this bias both qualitatively and quantitatively by using two different methods of inhomogeneity correction. We also show that this bias in standardization is independent of the specific inhomogeneity correction method used. The effect of this bias due to correction was also seen in magnetization transfer ratio (MTR) images, which are naturally endowed with the standardness property. Standardization, on the other hand, does not seem to influence the correction operation. It is also found that longer sequences of repeated correction and standardization operations do not considerably improve image quality. These results were found to hold for the clinical and the phantom data sets, for different MRI protocols, for different levels of artificial nonstandardness, for different models and levels of artificial inhomogeneity, for different correction methods, and for images that were endowed with inherent standardness as well as for those that were standardized by using the intensity standardization method. Overall, we conclude that inhomogeneity correction followed by intensity standardization is the best sequence to follow from the perspective of both image quality and computational efficiency.

Discovery and integrative neuroscience

Hypothesis driven research has been shown to be an excellent model for pursuing investigations in neuroscience. The Human Genome Project demonstrated the added value of discovery research, especially in areas where large amounts of data are produced. Neuroscience has become a data rich field, and one that would be enhanced by incorporating the discovery approach. Databases, as well as analytical, modeling and simulation tools, will have to be developed, and they will need to be interoperable and federated. This paper presents an overview of the development of the field of neuroscience databases and associate tools: Neuroinformatics. The primary focus is on the impact of NIH funding of this process. The important issues of data sharing, as viewed from the perspective of the scientist and private and public funding organizations, are discussed. Neuroinformatics will provide more than just a sophisticated array of information technologies to help scientists understand and integrate nervous system data. It will make available powerful models of neural functions and facilitate discovery, hypothesis formulation and electronic collaboration.

Non-rigid image registration: theory and practice

Image registration is an important enabling technology in medical image analysis. The current emphasis is on development and validation of application-specific non-rigid techniques, but there is already a plethora of techniques and terminology in use. In this paper we discuss the current state of the art of non-rigid registration to put on-going research in context and to highlight current and future clinical applications that might benefit from this technology. The philosophy and motivation underlying non-rigid registration is discussed and a guide to common terminology is presented. The core components of registration systems are described and outstanding issues of validity and validation are confronted.

Quantitative MR imaging of the neocortex

This article provides an overview of novel MR image analysis methods applied to the quantitative assessment of the neocortex in various forms of epilepsy. Postacquisition processing methods, such as voxel-based morphometry and texture analysis, involve the use of computer software to manipulate, enhance, and classify image information in a digital format. These techniques have the potential to demonstrate subtle abnormalities that are not identified by eye because of anatomic variability. Information provided by quantitative MR imaging of the neocortex may be important for the identification of accurate predictors of surgical outcome and may refine the selection of surgical candidates, particularly those with "nonlesional" neocortical epilepsy.

A relational-tubular (ReTu) deformable model for vasculature quantification of zebrafish embryo from microangiography image series

Embryonic cardiovascular system plays a vital role in embryonic development of human and animal. In this work, we introduce a novel deformable model, which we called Relational-tubular (ReTu) deformable model for segmenting and quantifying the embryonic vasculature of zebrafish embryo from microangiography image series. Particularly, to incorporate additional constraints on the spatial relationships among vessel branches, we introduce a new energy term called relation energy into the model energy function. This energy item acts as a repulsion force between neighboring vessels during the deformation to encourage them to move towards their respective volume data. Using the ReTu deformable model, the deformation process is an iterative two-stage procedure: vascular axis deformation and vascular surface deformation. The efficiency and robustness of this approach are demonstrated by experiments which show that satisfactory quantifications of the vasculature can be obtained after 3-4 iterations.

Determining correspondence in 3-D MR brain images using attribute vectors as morphological signatures of voxels

Finding point correspondence in anatomical images is a key step in shape analysis and deformable registration. This paper proposes an automatic correspondence detection algorithm for intramodality MR brain images of different subjects using wavelet-based attribute vectors (WAVs) defined on every image voxel. The attribute vector (AV) is extracted from the wavelet subimages and reflects the image structure in a large neighborhood around the respective voxel in a multiscale fashion. It plays the role of a morphological signature for each voxel, and our goal is, therefore, to make it distinctive of the respective voxel. Correspondence is then determined from similarities of AVs. By incorporating the prior knowledge of the spatial relationship among voxels, the ability of the proposed algorithm to find anatomical correspondence is further improved. Experiments with MR images of human brains show that the algorithm performs similarly to experts, even for complex cortical structures.

Deformable registration of cortical structures via hybrid volumetric and surface warping

Registration of cortical structures across individuals is a very important step for quantitative analysis of the human brain cortex. This paper presents a method for deformable registration of cortical structures across individuals, using hybrid volumetric and surface warping. In the first step, a feature-based volumetric registration algorithm is used to warp a model cortical surface to the individual's space. This step greatly reduces the variation between the model and individual, thus providing a good initialization for the next step of surface warping. In the second step, a surface registration method, based on matching geometric attributes, warps the model surface to the individual. Point correspondences are also established at this step. The attribute vector, as the morphological signature of surface, was designed to be as distinctive as possible, so that each vertex on the model surface can find its correspondence on the individual surface. Experimental results on both synthesized and real brain data demonstrate the performance of the proposed method in the registration of cortical structures across individuals.

Optimized prostate biopsy via a statistical atlas of cancer spatial distribution

A methodology is presented for constructing a statistical atlas of spatial distribution of prostate cancer from a large patient cohort, and it is used for optimizing needle biopsy. An adaptive-focus deformable model is used for the spatial normalization and registration of 100 prostate histological samples, which were provided by the Center for Prostate Disease Research of the US Department of Defense, resulting in a statistical atlas of spatial distribution of prostate cancer. Based on this atlas, a statistical predictive model was developed to optimize the needle biopsy sites, by maximizing the probability of detecting cancer. Experimental results using cross-validation show that the proposed method can detect cancer with a 99% success rate using seven needles, in these samples.

Estimating topology preserving and smooth displacement fields

We propose a method for enforcing topology preservation and smoothness onto a given displacement field. We first analyze the conditions for topology preservation on two- and three-dimensional displacement fields over a discrete rectangular grid. We then pose the problem of finding the closest topology preserving displacement field in terms of its complete set of gradients, which we later solve using a cyclic projections framework. Adaptive smoothing of a displacement field is then formulated as an extension of topology preservation, via constraints imposed on the Jacobian of the displacement field. The simulation results indicate that this technique is a fast and reliable method to estimate a topology preserving displacement field from a noisy observation that does not necessarily preserve topology. They also show that the proposed smoothing method can render morphometric analysis methods that are based on displacement field of shape transformations more robust to noise without removing important morphologic characteristics.

A Bayesian morphometry algorithm

Most methods for structure-function analysis of the brain in medical images are usually based on voxel-wise statistical tests performed on registered magnetic resonance (MR) images across subjects. A major drawback of such methods is the inability to accurately locate regions that manifest nonlinear associations with clinical variables. In this paper, we propose Bayesian morphological analysis methods, based on a Bayesian-network representation, for the analysis of MR brain images. First, we describe how Bayesian networks (BNs) can represent probabilistic associations among voxels and clinical (function) variables. Second, we present a model-selection framework, which generates a BN that captures structure-function relationships from MR brain images and function variables. We demonstrate our methods in the context of determining associations between regional brain atrophy (as demonstrated on MR images of the brain), and functional deficits. We employ two data sets for this evaluation: the first contains MR images of 11 subjects, where associations between regional atrophy and a functional deficit are almost linear; the second data set contains MR images of the ventricles of 84 subjects, where the structure-function association is nonlinear. Our methods successfully identify voxel-wise morphological changes that are associated with functional deficits in both data sets, whereas standard statistical analysis (i.e., t-test and paired t-test) fails in the nonlinear-association case.

A new algorithm of brain volume contours segmentation

This paper explores brain CT slices segmentation technique and some related problems, including contours segmentation algorithms, edge detector, algorithm evaluation and experimental results. This article describes a method for contour-baed segmentation of anatomical structures in 3D medical data sets. With this method, the user manually traces one or more 2D contours of an anatomical structure of interest on parallel planes arbitrarily cutting the data set. The experimental results show the segmentation based on 3D brain volume and 2D CT slices. The main creative contributions in this paper are: (1) contours segmentation algorithm; (2) edge detector; (3) algorithm evaluation.

Segmentation of prostate boundaries from ultrasound images using statistical shape model

This paper presents a statistical shape model for the automatic prostate segmentation in transrectal ultrasound images. A Gabor filter bank is first used to characterize the prostate boundaries in ultrasound images in both multiple scales and multiple orientations. The Gabor features are further reconstructed to be invariant to the rotation of the ultrasound probe and incorporated in the prostate model as image attributes for guiding the deformable segmentation. A hierarchical deformation strategy is then employed, in which the model adaptively focuses on the similarity of different Gabor features at different deformation stages using a multiresolution technique, i.e., coarse features first and fine features later. A number of successful experiments validate the algorithm.

Hierarchical active shape models, using the wavelet transform

Active shape models (ASMs) are often limited by the inability of relatively few eigenvectors to capture the full range of biological shape variability. This paper presents a method that overcomes this limitation, by using a hierarchical formulation of active shape models, using the wavelet transform. The statistical properties of the wavelet transform of a deformable contour are analyzed via principal component analysis, and used as priors in the contour's deformation. Some of these priors reflect relatively global shape characteristics of the object boundaries, whereas, some of them capture local and high-frequency shape characteristics and, thus, serve as local smoothness constraints. This formulation achieves two objectives. First, it is robust when only a limited number of training samples is available. Second, by using local statistics as smoothness constraints, it eliminates the need for adopting ad hoc physical models, such as elasticity or other smoothness models, which do not necessarily reflect true biological variability. Examples on magnetic resonance images of the corpus callosum and hand contours demonstrate that good and fully automated segmentations can be achieved, even with as few as five training samples.

A framework for callosal fiber distribution analysis

This paper presents a framework for analyzing the spatial distribution of neural fibers in the brain, with emphasis on interhemispheric fiber bundles crossing through the corpus callosum. The proposed approach combines methodologies for fiber tracking and spatial normalization and is applied on diffusion tensor images and standard magnetic resonance images.

HAMMER: hierarchical attribute matching mechanism for elastic registration

A new approach is presented for elastic registration of medical images, and is applied to magnetic resonance images of the brain. Experimental results demonstrate very high accuracy in superposition of images from different subjects. There are two major novelties in the proposed algorithm. First, it uses an attribute vector, i.e., a set of geometric moment invariants (GMIs) that are defined on each voxel in an image and are calculated from the tissue maps, to reflect the underlying anatomy at different scales. The attribute vector, if rich enough, can distinguish between different parts of an image, which helps establish anatomical correspondences in the deformation procedure; it also helps reduce local minima, by reducing ambiguity in potential matches. This is a fundamental deviation of our method, referred to as the hierarchical attribute matching mechanism for elastic registration (HAMMER), from other volumetric deformation methods, which are typically based on maximizing image similarity. Second, in order to avoid being trapped by local minima, i.e., suboptimal poor matches, HAMMER uses a successive approximation of the energy function being optimized by lower dimensional smooth energy functions, which are constructed to have significantly fewer local minima. This is achieved by hierarchically selecting the driving features that have distinct attribute vectors, thus, drastically reducing ambiguity in finding correspondence. A number of experiments demonstrate that the proposed algorithm results in accurate superposition of image data from individuals with significant anatomical differences.

On the metrics and euler-lagrange equations of computational anatomy

This paper reviews literature, current concepts and approaches in computational anatomy (CA). The model of CA is a Grenander deformable template, an orbit generated from a template under groups of diffeomorphisms. The metric space of all anatomical images is constructed from the geodesic connecting one anatomical structure to another in the orbit. The variational problems specifying these metrics are reviewed along with their associated Euler-Lagrange equations. The Euler equations of motion derived by Arnold for the geodesics in the group of divergence-free volume-preserving diffeomorphisms of incompressible fluids are generalized for the larger group of diffeomorphisms used in CA with nonconstant Jacobians. Metrics that accommodate photometric variation are described extending the anatomical model to incorporate the construction of neoplasm. Metrics on landmarked shapes are reviewed as well as Joshi's diffeomorphism metrics, Bookstein's thin-plate spline approximate-metrics, and Kendall's affine invariant metrics. We conclude by showing recent experimental results from the Toga & Thompson group in growth, the Van Essen group in macaque and human cortex mapping, and the Csernansky group in hippocampus mapping for neuropsychiatric studies in aging and schizophrenia.

Spatial normalization of spine MR images for statistical correlation of lesions with clinical symptoms

An image analysis method was developed for spatial normalization of spine magnetic resonance images. A deformable shape model of the spine is first constructed, and it is subsequently used by an automated algorithm to find a shape transformation that places patient data into a stereotactic space. Very good agreement with manual segmentations was observed. The main application of this method is in lesion-deficit analysis for determining associations between structural damage and clinical symptoms.

3-D active appearance models: segmentation of cardiac MR and ultrasound images

A model-based method for three-dimensional image segmentation was developed and its performance assessed in segmentation of volumetric cardiac magnetic resonance (MR) images and echocardiographic temporal image sequences. Comprehensive design of a three-dimensional (3-D) active appearance model (AAM) is reported for the first time as an involved extension of the AAM framework introduced by Cootes et al. The model's behavior is learned from manually traced segmentation examples during an automated training stage. Information about shape and image appearance of the cardiac structures is contained in a single model. This ensures a spatially and/or temporally consistent segmentation of three-dimensional cardiac images. The clinical potential of the 3-D AAM is demonstrated in short-axis cardiac MR images and four-chamber echocardiographic sequences. The method's performance was assessed by comparison with manually identified independent standards in 56 clinical MR and 64 clinical echo image sequences. The AAM method showed good agreement with the independent standard using quantitative indexes of border positioning errors, endo- and epicardial volumes, and left ventricular mass. In MR, the endocardial volumes, epicardial volumes, and left ventricular wall mass correlation coefficients between manual and AAM were R2 = 0.94, 0.97, 0.82, respectively. For echocardiographic analysis, the area correlation was R2 = 0.79. The AAM method shows high promise for successful application to MR and echocardiographic image analysis in a clinical setting.

Measuring size and shape of the hippocampus in MR images using a deformable shape model

A method for segmentation and quantification of the shape and size of the hippocampus is proposed, based on an automated image analysis algorithm. The algorithm uses a deformable shape model to locate the hippocampus in magnetic resonance images and to determine a geometric representation of its boundary. The deformable model combines three types of information. First, it employs information about the geometric properties of the hippocampal boundary, from a local and relatively finer scale to a more global and relatively coarser scale. Second, the model includes a statistical characterization of normal shape variation across individuals, serving as prior knowledge to the algorithm. Third, the algorithm utilizes a number of manually defined boundary points, which can help guide the model deformation to the appropriate boundaries, wherever these boundaries are weak or not clearly defined in MR images. Excellent agreement is demonstrated between the algorithm and manual segmentations by well-trained raters, with a correlation coefficient equal to 0.97 and algorithm/rater differences statistically equivalent to interrater differences for manual definitions.

A framework for predictive modeling of anatomical deformations

A framework for modeling and predicting anatomical deformations is presented, and tested on simulated images. Although a variety of deformations can be modeled in this framework, emphasis is placed on surgical planning, and particularly on modeling and predicting changes of anatomy between preoperative and intraoperative positions, as well as on deformations induced by tumor growth. Two methods are examined. The first is purely shape-based and utilizes the principal modes of co-variation between anatomy and deformation in order to statistically represent deformability. When a patient's anatomy is available, it is used in conjunction with the statistical model to predict the way in which the anatomy will/can deform. The second method is related, and it uses the statistical model in conjunction with a biomechanical model of anatomical deformation. It examines the principal modes of co-variation between shape and forces, with the latter driving the biomechanical model, and thus predicting deformation. Results are shown on simulated images, demonstrating that systematic deformations, such as those resulting from change in position or from tumor growth, can be estimated very well using these models. Estimation accuracy will depend on the application, and particularly on how systematic a deformation of interest is.