Publications by Year: 2005

OBJECTIVE: The purpose of this study was to investigate abnormalities in the surface complexity of the prefrontal cortex and in the hemispheric asymmetry of cortical complexity in first-episode patients with schizophrenia.

We present a fast, robust and automatic method for computing centerline paths through tubular structures for application to virtual endoscopy. The key idea is to utilize a skeletonization algorithm which exploits properties of the average outward flux of the gradient vector field of a Euclidean distance function from the boundary of the structure. The algorithm is modified to yield a collection of 3D curves, each of which is locally centered. The approach requires no user interaction, is virtually parameter free and has low computational complexity. We validate the method quantitatively on a number of synthetic data sets with known centerlines and qualitatively on colon and vessel data segmented from CT and CRA images.

In magnetic resonance imaging (MRI) research, significant attention has been paid to the analysis of the hippocampus (HC) within the medial temporal lobe because of its importance in memory and learning, and its role in neurodegenerative diseases. Manual segmentation protocols have established a volume decline in the HC in conjunction with Alzheimer’s disease, epilepsy, post-traumatic stress disorder, and depression. Furthermore, recent studies have investigated age-related changes of HC volume which show an interaction with gender; in early adulthood, volume reduction of the HC is found in men but not in women. In this paper, we investigated gender differences in normal subjects in young adulthood by employing a shape analysis of the HC using medial surfaces. For each subject, the most prominent medial manifold of the HC was extracted and flattened. The flattened sheets were then registered using both a rigid and a non-rigid alignment technique, and the medial surface radius was expressed as a height function over them. This allowed for an investigation of the association between subject variables and the local width of the HC. With regard to the effects of age and gender, it could be shown that the previously observed gender differences were mostly due to volume loss in males in the lateral areas of the HC head and tail. We suggest that the analysis of HC shape using medial surfaces might thus serve as a complimentary technique to investigate group differences to the established segmentation protocols for volume quantification in MRI.

In this paper, we present an evaluation of seven automatic brain tissue classifiers based on level of agreements. A number of agreement measures are explained, and we show how they can be used to compare different segmentation techniques. We use the Simultaneous Truth and Performance Level Estimation (STAPLE) of Warfield et al. but also introduce a novel evaluation technique based on the Williams’ index. The methods are evaluated using these two techniques on a population of forty subjects, each having an SPGR scan and a co-registered T2 weighted scan. We provide an interpretation of the results and show how similar the output of the STAPLE analysis and Williams’ index are. When no ground truth is required, we recommend the use of Williams’ index as it is easy and fast to compute.

This paper details the derivation of rotationally invariant scalar measures from higher-rank diffusion tensors (DTs) and functions defined on a unit sphere. This was accomplished with the use of an expression that generalizes the evaluation of the trace operator to tensors of arbitrary rank, and even to functions whose domains are the unit sphere. It is shown that the mean diffusivity is invariant to the selection of tensor rank for the model used. However, this rank invariance does not apply to the anisotropy measures. Therefore, a variance-based, general anisotropy measure is proposed. Also an information theoretical parametrization of anisotropy is introduced that is frequently more consistent with the meaning attributed to anisotropy. We accomplished this by associating anisotropy with the amount of orientational information present in the data, regardless of the imaging technique used. Using a simplified model of fibrous tissue, we simulated anisotropy values with varying orientational complexity and tensor models. Simulations suggested that a lower-rank tensor model may produce artificially low anisotropy values in voxels with complex structure. This was confirmed with a spin-echo experiment performed on an excised rat brain.

This paper introduces a new, accurate and fast method for fiber orientation mapping using high angular resolution diffusion imaging (HARDI) data. The approach utilizes the Fourier relationship between the water displacement probabilities and diffusion attenuated magnetic resonance (MR) signal expressed in spherical coordinates. The Laplace series coefficients of the water displacement probabilities are evaluated at a fixed distance away from the origin. The computations take under one minute for most three-dimensional datasets. We present orientation maps computed from excised rat optic chiasm, brain and spinal cord images. The developed method will improve the reliability of tractography schemes and make it possible to correctly identify the neural connections between functionally connected regions of the nervous system.

We present a novel method for finding white matter fiber correspondences and clusters across a population of brains. Our input is a collection of paths from tractography in every brain. Using spectral methods we embed each path as a vector in a high dimensional space. We create the embedding space so that it is common across all brains, consequently similar paths in all brains will map to points near each other in the space. By performing clustering in this space we are able to find matching fiber tract clusters in all brains. In addition, we automatically obtain correspondence of tractographic paths across brains: by selecting one or several paths of interest in one brain, the most similar paths in all brains are obtained as the nearest points in the high-dimensional space.

This paper presents an extension of the phase correlation image alignment method to N-dimensional data sets. By the Fourier shift theorem, the motion model for translational shifts between N-dimensional images can be represented as a rank-one tensor. Through use of a high-order singular value decomposition, the phase correlation between two N-dimensional data sets can be decomposed to independently identify translational displacements along each dimension with subpixel resolution. Using three-dimensional MRI data sets, we demonstrate the effectiveness of this approach relative to other N-dimensional image registration methods.

Emphysema is characterized by the destruction and over distension of lung tissue, which manifest on high resolution computer tomography (CT) images as regions of low attenuation. Typically, it is diagnosed by clinical symptoms, physical examination, pulmonary function tests, and X-ray and CT imaging. In this paper we discuss a quantitative imaging approach to analyze emphysema which employs low-level segmentations of CT images that partition the data into perceptually relevant regions. We constructed multi-dimensional histograms of feature values computed over the image segmentation. For each region in the segmentation, we derive a rich set of feature measurements. While we can use any combination of physical and geometric features, we found that limiting the scope to two features - the mean attenuation across a region and the region area - is effective. The subject histogram is compared to a set of canonical histograms representative of various stages of emphysema using the Earth Mover’s Distance metric. Disease severity is assigned based on which canonical histogram is most similar to the subject histogram. Experimental results with 81 cases of emphysema at different stages of disease progression show good agreement against the reading of an expert radiologist.

We investigated a method, motion compensated integration (MCI), for enhancing stent Contrast-to-Noise Ratio (CNR) such that stent deployment may be more easily assessed. MCI registers fluoroscopic frames on the basis of stent motion and performs pixel-wise integration to reduce noise. Registration is based on marker balls, high contrast interventional devices which guide the clinician in stent placement. It is assumed that stent motion is identical to that of the marker balls. Detecting marker balls and identifying their centroids with a high degree of accuracy is a non-trivial task. To address the required registration accuracy, in this work we examine MCI’s visualization benefit as a function of registration error. We employ adaptive forced choice experiments to quantify human discrimination fidelity. Perception results are contrasted with CNR measurements. For each level of registration inaccuracy investigated, MCI conferred a benefit (p 0.05) on stent deployment assessment suggesting the technique is tolerant of modest registration error. We also consider the blurring effect of cardiac motion during the x-ray pulse and select frames for integration as a function of cardiac phase imaged.