Publications

2007
Niethammer M, Bouix S, Aja-Fernández S, Westin C-F, Shenton ME. Outlier rejection for diffusion weighted imaging. Med Image Comput Comput Assist Interv. 2007;10 (Pt 1) :161-8.Abstract
This paper introduces an outlier rejection and signal reconstruction method for high angular resolution diffusion weighted imaging. The approach is based on the thresholding of Laplacian measurements over the sphere of the apparent diffusion coefficient profiles defined for a given set of gradient directions. Exemplary results are presented.
Ou W, Golland P, Hämäläinen M. Sources of variability in MEG. Med Image Comput Comput Assist Interv. 2007;10 (Pt 2) :751-9.Abstract
This paper investigates and characterizes sources of variability in MEG signals in multi-site, multi-subject studies. Understanding these sources will help to develop efficient strategies for comparing and pooling data across repetitions of an experiment, across subjects, and across sites. In this work, we investigated somatosensory MEG data collected at three different sites and applied variance component analysis and nonparametric KL divergence analysis in order to characterize the sources of variability. Our analysis showed that inter-subject differences are the biggest factor in the signal variability. We demonstrated that the timing of the deflections is very consistent in the early somatosensory response, which justifies a direct comparison of deflection peak times acquired from different visits, subjects, and systems. Compared with deflection peak times, deflection magnitudes have larger variation across sites; modeling of this variability is necessary for data pooling.
San José Estépar R, Stylopoulos N, Ellis R, Samset E, Westin C-F, Thompson C, Vosburgh K. Towards scarless surgery: an endoscopic ultrasound navigation system for transgastric access procedures. Comput Aided Surg. 2007;12 (6) :311-24.Abstract
OBJECTIVE: Scarless surgery is an innovative and promising technique that may herald a new era in surgical procedures. We have created a navigation system, named IRGUS, for endoscopic and transgastric access interventions and have validated it in in vivo pilot studies. Our hypothesis is that endoscopic ultrasound procedures will be performed more easily and efficiently if the operator is provided with approximately registered 3D and 2D processed CT images in real time that correspond to the probe position and ultrasound image. MATERIALS AND METHODS: The system provides augmented visual feedback and additional contextual information to assist the operator. It establishes correspondence between the real-time endoscopic ultrasound image and a preoperative CT volume registered using electromagnetic tracking of the endoscopic ultrasound probe position. Based on this positional information, the CT volume is reformatted in approximately the same coordinate frame as the ultrasound image and displayed to the operator. RESULTS: The system reduces the mental burden of probe navigation and enhances the operator's ability to interpret the ultrasound image. Using an initial rigid body registration, we measured the mis-registration error between the ultrasound image and the reformatted CT plane to be less than 5 mm, which is sufficient to enable the performance of novice users of endoscopic systems to approach that of expert users. CONCLUSIONS: Our analysis shows that real-time display of data using rigid registration is sufficiently accurate to assist surgeons in performing endoscopic abdominal procedures. By using preoperative data to provide context and support for image interpretation and real-time imaging for targeting, it appears probable that both preoperative and intraoperative data may be used to improve operator performance.
Michailovich O, Tannenbaum A. Blind deconvolution of medical ultrasound images: a parametric inverse filtering approach. IEEE Trans Image Process. 2007;16 (12) :3005-19.Abstract
The problem of reconstruction of ultrasound images by means of blind deconvolution has long been recognized as one of the central problems in medical ultrasound imaging. In this paper, this problem is addressed via proposing a blind deconvolution method which is innovative in several ways. In particular, the method is based on parametric inverse filtering, whose parameters are optimized using two-stage processing. At the first stage, some partial information on the point spread function is recovered. Subsequently, this information is used to explicitly constrain the spectral shape of the inverse filter. From this perspective, the proposed methodology can be viewed as a "hybridization" of two standard strategies in blind deconvolution, which are based on either concurrent or successive estimation of the point spread function and the image of interest. Moreover, evidence is provided that the "hybrid" approach can outperform the standard ones in a number of important practical cases. Additionally, the present study introduces a different approach to parameterizing the inverse filter. Specifically, we propose to model the inverse transfer function as a member of a principal shift-invariant subspace. It is shown that such a parameterization results in considerably more stable reconstructions as compared to standard parameterization methods. Finally, it is shown how the inverse filters designed in this way can be used to deconvolve the images in a nonblind manner so as to further improve their quality. The usefulness and practicability of all the introduced innovations are proven in a series of both in silico and in vivo experiments. Finally, it is shown that the proposed deconvolutioh algorithms are capable of improving the resolution of ultrasound images by factors of 2.24 or 6.52 (as judged by the autocorrelation criterion) depending on the type of regularization method used.
Georgiou T, Michailovich O, Rathi Y, Malcolm J, Tannenbaum A. Distribution Metrics and Image Segmentation. Linear Algebra Appl. 2007;425 (2-3) :663-672.Abstract
The purpose of this paper is to describe certain alternative metrics for quantifying distances between distributions, and to explain their use and relevance in visual tracking. Besides the theoretical interest, such metrics may be used to design filters for image segmentation, that is for solving the key visual task of separating an object from the background in an image. The segmenting curve is represented as the zero level set of a signed distance function. Most existing methods in the geometric active contour framework perform segmentation by maximizing the separation of intensity moments between the interior and the exterior of an evolving contour. Here one can use the given distributional metric to determine a flow which minimizes changes in the distribution inside and outside the curve.
2006
MacFall JR, Taylor WD, Rex DE, Pieper S, Payne ME, McQuoid DR, Steffens DC, Kikinis R, Toga AW, Krishnan RRK. Lobar distribution of lesion volumes in late-life depression: the Biomedical Informatics Research Network (BIRN). Neuropsychopharmacology. 2006;31 (7) :1500-7.Abstract
White matter hyperintense lesions on T2-weighted images are associated with late-life depression. Little work has been carried out examining differences in lesion location between elderly individuals with and without depression. In contrast to previous studies examining total brain white matter lesion volume, this study examined lobar differences in white matter lesion volumes derived from brain magnetic resonance imaging. This study examined 49 subjects with a DSM-IV diagnosis of major depression and 50 comparison subjects without depression. All participants were age 60 years or older. White matter lesion volumes were measured in each hemisphere using a semiautomated segmentation process and localized to lobar regions using a lobar atlas created for this sample using the imaging tools provided by the Biomedical Informatics Research Network (BIRN). The lobar lesion volumes were compared against depression status. After controlling for age and hypertension, subjects with depression exhibited significantly greater total white matter lesion volume in both hemispheres and in both frontal lobes than did control subjects. Although a similar trend was observed in the parietal lobes, the difference did not reach a level of statistical significance. Models of the temporal and occipital lobes were not statistically significant. Older individuals with depression have greater white matter disease than healthy controls, predominantly in the frontal lobes. These changes are thought to disrupt neural circuits involved in mood regulation, thus increasing the risk of developing depression.
Pohl KM, Fisher J, Grimson EWL, Kikinis R, Wells III WM. A Bayesian Model for Joint Segmentation and Registration. Neuroimage. 2006;31 (1) :228-39.Abstract

A statistical model is presented that combines the registration of an atlas with the segmentation of magnetic resonance images. We use an Expectation Maximization-based algorithm to find a solution within the model, which simultaneously estimates image artifacts, anatomical labelmaps, and a structure-dependent hierarchical mapping from the atlas to the image space. The algorithm produces segmentations for brain tissues as well as their substructures. We demonstrate the approach on a set of 22 magnetic resonance images. On this set of images, the new approach performs significantly better than similar methods which sequentially apply registration and segmentation.

Guttmann CRG, Meier DS, Holland CM. Can MRI reveal phenotypes of multiple sclerosis?. Magn Reson Imaging. 2006;24 (4) :475-81.Abstract
The multicontrast capability of magnetic resonance imaging (MRI) is discussed in its role in the search for phenotypes of multiple sclerosis (MS). Aspects of MRI specificity, putative markers for pathogenetic components of disease and issues of spatial and temporal distribution are discussed. While particular reference is made to MS, the concepts apply to common pathological features of many neurologic diseases and to neurodegenerative disease in general. The assessment and dissociation of disease activity and disease severity, as well as the combination of varied metrics for the purposes of inferential and predictive disease modeling, are explored with respect to biomarkers and clinical outcomes. By virtue of its noninvasive nature and multicontrast capabilities depicting multiple facets of MS pathology, MRI lends itself to the systematic search of pathogenetically distinct subtypes of MS in large populations of patients. In conjunction with clinical, immunological, serological and genetic information, clusters of MS patients with distinct clinical prognosis and diverse response profiles to available and future treatments may be identified.
Learned-Miller EG. Data driven image models through continuous joint alignment. IEEE Trans Pattern Anal Mach Intell. 2006;28 (2) :236-50.Abstract
This paper presents a family of techniques that we call congealing for modeling image classes from data. The idea is to start with a set of images and make them appear as similar as possible by removing variability along the known axes of variation. This technique can be used to eliminate "nuisance" variables such as affine deformations from handwritten digits or unwanted bias fields from magnetic resonance images. In addition to separating and modeling the latent images-i.e., the images without the nuisance variables-we can model the nuisance variables themselves, leading to factorized generative image models. When nuisance variable distributions are shared between classes, one can share the knowledge learned in one task with another task, leading to efficient learning. We demonstrate this process by building a handwritten digit classifier from just a single example of each class. In addition to applications in handwritten character recognition, we describe in detail the application of bias removal from magnetic resonance images. Unlike previous methods, we use a separate, nonparametric model for the intensity values at each pixel. This allows us to leverage the data from the MR images of different patients to remove bias from each other. Only very weak assumptions are made about the distributions of intensity values in the images. In addition to the digit and MR applications, we discuss a number of other uses of congealing and describe experiments about the robustness and consistency of the method.
Michailovich OV, Tannenbaum A. Despeckling of medical ultrasound images. IEEE Trans Ultrason Ferroelectr Freq Control. 2006;53 (1) :64-78.Abstract
Speckle noise is an inherent property of medical ultrasound imaging, and it generally tends to reduce the image resolution and contrast, thereby reducing the diagnostic value of this imaging modality. As a result, speckle noise reduction is an important prerequisite, whenever ultrasound imaging is used for tissue characterization. Among the many methods that have been proposed to perform this task, there exists a class of approaches that use a multiplicative model of speckled image formation and take advantage of the logarithmical transformation in order to convert multiplicative speckle noise into additive noise. The common assumption made in a dominant number of such studies is that the samples of the additive noise are mutually uncorrelated and obey a Gaussian distribution. The present study shows conceptually and experimentally that this assumption is oversimplified and unnatural. Moreover, it may lead to inadequate performance of the speckle reduction methods. The study introduces a simple preprocessing procedure, which modifies the acquired radio-frequency images (without affecting the anatomical information they contain), so that the noise in the log-transformation domain becomes very close in its behavior to a white Gaussian noise. As a result, the preprocessing allows filtering methods based on assuming the noise to be white and Gaussian, to perform in nearly optimal conditions. The study evaluates performances of three different, nonlinear filters--wavelet denoising, total variation filtering, and anisotropic diffusion--and demonstrates that, in all these cases, the proposed preprocessing significantly improves the quality of resultant images. Our numerical tests include a series of computer-simulated and in vivo experiments.
Talos I-F, Mian AZ, Zou KH, Hsu L, Goldberg-Zimring D, Haker S, Bhagwat JG, Mulkern RV. Magnetic resonance and the human brain: anatomy, function and metabolism. Cell Mol Life Sci. 2006;63 (10) :1106-24.Abstract
The introduction and development, over the last three decades, of magnetic resonance (MR) imaging and MR spectroscopy technology for in vivo studies of the human brain represents a truly remarkable achievement, with enormous scientific and clinical ramifications. These effectively non-invasive techniques allow for studies of the anatomy, the function and the metabolism of the living human brain. They have allowed for new understandings of how the healthy brain works and have provided insights into the mechanisms underlying multiple disease processes which affect the brain. Different MR techniques have been developed for studying anatomy, function and metabolism. The primary focus of this review is to describe these different methodologies and to briefly review how they are being employed to more fully appreciate the intricacies associated with the organ, which most distinctly differentiates the human species from the other animal forms on earth.
Angenent S, Pichon E, Tannenbaum A. MATHEMATICAL METHODS IN MEDICAL IMAGE PROCESSING. Bull New Ser Am Math Soc. 2006;43 :365-396.Abstract
In this paper, we describe some central mathematical problems in medical imaging. The subject has been undergoing rapid changes driven by better hardware and software. Much of the software is based on novel methods utilizing geometric partial differential equations in conjunction with standard signal/image processing techniques as well as computer graphics facilitating man/machine interactions. As part of this enterprise, researchers have been trying to base biomedical engineering principles on rigorous mathematical foundations for the development of software methods to be integrated into complete therapy delivery systems. These systems support the more effective delivery of many image-guided procedures such as radiation therapy, biopsy, and minimally invasive surgery. We will show how mathematics may impact some of the main problems in this area, including image enhancement, registration, and segmentation.
Magnotta VA, Friedman L. Measurement of Signal-to-Noise and Contrast-to-Noise in the fBIRN Multicenter Imaging Study. J Digit Imaging. 2006;19 (2) :140-7.Abstract
The ability to analyze and merge data across sites, vendors, and field strengths depends on one's ability to acquire images with the same image quality including image smoothness, signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). SNR can be used to compare different magnetic resonance scanners as a measure of comparability between the systems. This study looks at the SNR and CNR ratios in structural fast spin-echo T2-weighted scans acquired in five individuals across ten sites that are part of Functional Imaging Research of Schizophrenia Testbed Biomedical Informatics Research Network (fBIRN). Different manufacturers, field strengths, gradient coils, and RF coils were used at these sites. The SNR of gray matter was fairly uniform (41.3-43.3) across scanners at 1.5 T. The higher field scanners produced images with significantly higher SNR values (44.5-108.7 at 3 T and 50.8 at 4 T). Similar results were obtained for CNR measurements between gray/white matter at 1.5 T (9.5-10.2), again increasing at higher fields (10.1-28.9 at 3 T and 10.9 at 4 T).
Dimaio SP, Kacher DF, Ellis RE, Fichtinger G, Hata N, Zientara GP, Panych LP, Kikinis R, Jolesz FA. Needle Artifact Localization in 3T MR Images. Stud Health Technol Inform. 2006;119 :120-5.Abstract

This work explores an image-based approach for localizing needles during MRI-guided interventions, for the purpose of tracking and navigation. Susceptibility artifacts for several needles of varying thickness were imaged, in phantoms, using a 3 tesla MRI system, under a variety of conditions. The relationship between the true needle positions and the locations of artifacts within the images, determined both by manual and automatic segmentation methods, have been quantified and are presented here.

Dambreville S, Rathi Y, Tannenbaum A. Shape-Based Approach to Robust Image Segmentation using Kernel PCA. Proc IEEE Comput Soc Conf Comput Vis Pattern Recognit. 2006 :977-984.Abstract
Segmentation involves separating an object from the background. In this work, we propose a novel segmentation method combining image information with prior shape knowledge, within the level-set framework. Following the work of Leventon et al., we revisit the use of principal component analysis (PCA) to introduce prior knowledge about shapes in a more robust manner. To this end, we utilize Kernel PCA and show that this method of learning shapes outperforms linear PCA, by allowing only shapes that are close enough to the training data. In the proposed segmentation algorithm, shape knowledge and image information are encoded into two energy functionals entirely described in terms of shapes. This consistent description allows to fully take advantage of the Kernel PCA methodology and leads to promising segmentation results. In particular, our shape-driven segmentation technique allows for the simultaneous encoding of multiple types of shapes, and offers a convincing level of robustness with respect to noise, clutter, partial occlusions, or smearing.
Talos I-F, Zou KH, Ohno-Machado L, Bhagwat JG, Kikinis R, Black PM, Jolesz FA. Supratentorial low-grade glioma resectability: statistical predictive analysis based on anatomic MR features and tumor characteristics. Radiology. 2006;239 (2) :506-13.Abstract
PURPOSE: To retrospectively assess the main variables that affect the complete magnetic resonance (MR) imaging-guided resection of supratentorial low-grade gliomas. MATERIALS AND METHODS: Institutional review board approval was obtained for this retrospective HIPAA-compliant study, with the requirement for informed consent waived. Data from 101 patients (61 men, 40 women; mean age, 39 years; age range, 18-72 years) who had nonenhancing supratentorial mass lesions that were histopathologically diagnosed as low-grade (World Health Organization grade II) gliomas and consecutively underwent surgery with intraoperative MR imaging guidance were analyzed. There were 21 low-grade astrocytomas, 64 oligodendrogliomas, and 16 mixed oligoastrocytomas. Initial and residual tumor volumes were measured on intraoperative T2-weighted MR images and three-dimensional spoiled gradient-echo MR images. The anatomic relationships between the tumor and eloquent cortical and/or subcortical regions and the influence of these relationships on the extent of resection were analyzed on the basis of preoperative MR imaging findings. Summary measures, univariate Fisher exact test and t test, and multivariate logistic regression analyses were performed. RESULTS: Tumor volume ranged from 2.7-231.0 mL. Univariate analyses revealed the following tumor characteristics to be significant predictive variables of incomplete tumor resection: diffuse tumor margin on T2-weighted MR images, oligodendroglioma or oligoastrocytoma histopathologic type, and large tumor volume (P < .05 for all). Tumor involvement of the following structures was associated with incomplete resection: corpus callosum, corticospinal tract, insular lobe, middle cerebral artery, motor cortex, optic radiation, visual cortex, and basal ganglia (P < .05 for all). Multivariate analyses revealed that incomplete tumor resection was due to tumor involvement of the corticospinal tract (P < .01), large tumor volume (P < .01), and oligodendroglioma histopathologic type (P = .02). CONCLUSION: The main variables associated with incomplete tumor resection in 101 patients were identified by using statistical predictive analyses.
Archip N, Rohling R, Dessenne V, Erard P-J, Nolte LP. Anatomical structure modeling from medical images. Comput Methods Programs Biomed. 2006;82 (3) :203-15.Abstract
Some clinical applications, such as surgical planning, require volumetric models of anatomical structures represented as a set of tetrahedra. A practical method of constructing anatomical models from medical images is presented. The method starts with a set of contours segmented from the medical images by a clinician and produces a model that has high fidelity with the contours. Unlike most modeling methods, the contours are not restricted to lie on parallel planes. The main steps are a 3D Delaunay tetrahedralization, culling of non-object tetrahedra, and refinement of the tetrahedral mesh. The result is a high-quality set of tetrahedra whose surface points are guaranteed to match the original contours. The key is to use the distance map and bit volume structures that were created along with the contours. The method is demonstrated on computed tomography, MRI and 3D ultrasound data. Models of 170,000 tetrahedra are constructed on a standard workstation in approximately 10s. A comparison with related methods is also provided.
Wu Y, Warfield SK, Tan LI, Wells WM, Meier DS, van Schijndel RA, Barkhof F, Guttmann CRG. Automated segmentation of multiple sclerosis lesion subtypes with multichannel MRI. Neuroimage. 2006;32 (3) :1205-15.Abstract
PURPOSE: To automatically segment multiple sclerosis (MS) lesions into three subtypes (i.e., enhancing lesions, T1 "black holes", T2 hyperintense lesions). MATERIALS AND METHODS: Proton density-, T2- and contrast-enhanced T1-weighted brain images of 12 MR scans were pre-processed through intracranial cavity (IC) extraction, inhomogeneity correction and intensity normalization. Intensity-based statistical k-nearest neighbor (k-NN) classification was combined with template-driven segmentation and partial volume artifact correction (TDS+) for segmentation of MS lesions subtypes and brain tissue compartments. Operator-supervised tissue sampling and parameter calibration were performed on 2 randomly selected scans and were applied automatically to the remaining 10 scans. Results from this three-channel TDS+ (3ch-TDS+) were compared to those from a previously validated two-channel TDS+ (2ch-TDS+) method. The results of both the 3ch-TDS+ and 2ch-TDS+ were also compared to manual segmentation performed by experts. RESULTS: Intra-class correlation coefficients (ICC) of 3ch-TDS+ for all three subtypes of lesions were higher (ICC between 0.95 and 0.96) than that of 2ch-TDS+ for T2 lesions (ICC = 0.82). The 3ch-TDS+ also identified the three lesion subtypes with high specificity (98.7-99.9%) and accuracy (98.5-99.9%). Sensitivity of 3ch-TDS+ for T2 lesions was 16% higher than with 2ch-TDS+. Enhancing lesions were segmented with the best sensitivity (81.9%). "Black holes" were segmented with the least sensitivity (62.3%). CONCLUSION: 3ch-TDS+ is a promising method for automated segmentation of MS lesion subtypes.
Friman O, Farnebäck G, Westin C-F. A Bayesian approach for stochastic white matter tractography. IEEE Trans Med Imaging. 2006;25 (8) :965-78.Abstract
White matter fiber bundles in the human brain can be located by tracing the local water diffusion in diffusion weighted magnetic resonance imaging (MRI) images. In this paper, a novel Bayesian modeling approach for white matter tractography is presented. The uncertainty associated with estimated white matter fiber paths is investigated, and a method for calculating the probability of a connection between two areas in the brain is introduced. The main merits of the presented methodology are its simple implementation and its ability to handle noise in a theoretically justified way. Theory for estimating global connectivity is also presented, as well as a theorem that facilitates the estimation of the parameters in a constrained tensor model of the local water diffusion profile.
Szymczak A, Stillman A, Tannenbaum A, Mischaikow K. Coronary vessel trees from 3D imagery: a topological approach. Med Image Anal. 2006;10 (4) :548-59.Abstract
We propose a simple method for reconstructing vascular trees from 3D images. Our algorithm extracts persistent maxima of the intensity on all axis-aligned 2D slices of the input image. The maxima concentrate along 1D intensity ridges, in particular along blood vessels. We build a forest connecting the persistent maxima with short edges. The forest tends to approximate the blood vessels present in the image, but also contains numerous spurious features and often fails to connect segments belonging to one vessel in low contrast areas. We improve the forest by applying simple geometric filters that trim short branches, fill gaps in blood vessels and remove spurious branches from the vascular tree to be extracted. Experiments show that our technique can be applied to extract coronary trees from heart CT scans.

Pages