Difference between revisions of "Projects:MultimodalAtlas"

From NAMIC Wiki
Jump to: navigation, search
Line 1: Line 1:
 
  Back to [[NA-MIC_Internal_Collaborations:StructuralImageAnalysis|NA-MIC Collaborations]], [[Algorithm:MIT|MIT Algorithms]], [[DBP2:Harvard|Harvard DBP2]]
 
  Back to [[NA-MIC_Internal_Collaborations:StructuralImageAnalysis|NA-MIC Collaborations]], [[Algorithm:MIT|MIT Algorithms]], [[DBP2:Harvard|Harvard DBP2]]
 
__NOTOC__
 
__NOTOC__
−
In many medical image analysis applications we have a large collection of data sets available to us. In such applications, we are faced with at least two different challenges: first, most analyses we are interested in carrying out on the data-sets require point-wise spatial correspondence across the images. Co-registration of the image data set directly or indirectly attempts to solve this problem by spatially warping the images into a common, so-called “atlas space.” Another important challenge we face is coming up with a summary of the collection – this summary is usually called a template or an atlas, which is usually an average image. The famous MNI template brain, commonly used in neuroimaging studies is one example. Templates can be used to “normalize” new images, that is bring them into the common atlas space: for localizing regions of interests in the new subject or comparing groups of interest.  
+
Today, computational anatomy studies are mainly hypothesis-driven, aiming to identify and characterize structural or functional differences between, for instance a group of patients with a specific disorder and control subjects. This type of approach has two premises: clinical classification of the subjects and spatial correspondence across the images. In practice, achieving either can be challenging. First, the complex spectrum of symptoms of neuro-degenerative disorders like schizophrenia and overlapping symptoms across different types of dementia like Alzheimer's disease, delirium and depression make a diagnosis based on standardized clinical tests like the mental status examination difficult. Second, across-subject correspondence in the images is a particularly hard problem that requires different approaches in various contexts. A popular technique is to normalize all subjects into a standard space, such as the Talairach space, by registering each image with a single, universal template image that represents an average brain. However, the quality of such an approach is limited by the accuracy with which the universal template represents the population in the study.
  
−
Recent work has acknowledged the fact that in some cases a single template may not suffice to summarize the whole population. For example, in Blezek and Miller MICCAI '06, the authors propose an atlas stratification algorithm that computes multiple templates. In this work, we propose and investigate an algorithm that jointly co-registers a collection of images while computing multiple templates. The algorithm, called '''iCluster''' for Image Clustering, is based on the following idea: given the templates, the co-registration problem becomes simple, reducing to a number of pairwise registration instances. On the other hand, if we are given a collection of images that have been co-registered, i.e. in perfect alignment, then an off-the shelf clustering or averaging algorithm can be used to compute the templates.  
+
With the increasing availability of medical images, data-driven algorithms offer the ability to probe a population and potentially discover sub-groups that may differ in unexpected ways. In this paper, we propose and demonstrate an efficient probabilistic clustering algorithm, called '''iCluster''', that:
 +
 
 +
* computes a small number of templates that summarize a given population of images,
 +
* simultaneously co-registers all the images using a nonlinear transformation model,
 +
* assigns each input image to a template.
 +
 
 +
The templates are guaranteed to live in an affine-normalized space, i.e., they are spatially aligned with respect to an affine transformation model.
  
 
= Description =
 
= Description =

Revision as of 20:09, 16 May 2008

Home < Projects:MultimodalAtlas
Back to NA-MIC Collaborations, MIT Algorithms, Harvard DBP2

Today, computational anatomy studies are mainly hypothesis-driven, aiming to identify and characterize structural or functional differences between, for instance a group of patients with a specific disorder and control subjects. This type of approach has two premises: clinical classification of the subjects and spatial correspondence across the images. In practice, achieving either can be challenging. First, the complex spectrum of symptoms of neuro-degenerative disorders like schizophrenia and overlapping symptoms across different types of dementia like Alzheimer's disease, delirium and depression make a diagnosis based on standardized clinical tests like the mental status examination difficult. Second, across-subject correspondence in the images is a particularly hard problem that requires different approaches in various contexts. A popular technique is to normalize all subjects into a standard space, such as the Talairach space, by registering each image with a single, universal template image that represents an average brain. However, the quality of such an approach is limited by the accuracy with which the universal template represents the population in the study.

With the increasing availability of medical images, data-driven algorithms offer the ability to probe a population and potentially discover sub-groups that may differ in unexpected ways. In this paper, we propose and demonstrate an efficient probabilistic clustering algorithm, called iCluster, that:

  • computes a small number of templates that summarize a given population of images,
  • simultaneously co-registers all the images using a nonlinear transformation model,
  • assigns each input image to a template.

The templates are guaranteed to live in an affine-normalized space, i.e., they are spatially aligned with respect to an affine transformation model.

Description

Generative Model used in iCluster.

iCluster is derived from a simple generative model. We assume that there are a fixed and known number of template images. Then the process that generates an observed image is as follows: a template is randomly drawn – note that the probability that governs this process doesn’t have to be uniform. Next, the chosen template is warped with a random transformation and i.i.d Gaussian noise is added to this warped image to generate an observed image. This process is repeated multiple times to generate a collection of images.

We formulate the problem as a maximum likelihood solution. We employ a Generalized Maximum Likelihood algorithm to solve the problem.

E-step.

In the E-step, the algorithm computes membership probabilities.

T-step.

We divide the M-step into two sub-steps: the first one, called the T (template) step is where we have closed-form expressions for the optimization. The algorithm updates the template images as weighted averages of the images, where the weights are the memberships. The template priors are updated as normalized sums of memberships.

R-step.

In the second part of the M-step, called R (registration) step, the algorithm performs a collection of independent pairwise registration instances. Each image is registered with a unique, "effective" template, which is a weighted average of the current templates and the weights are the corresponding memberships.

Results

We ran iCluster with k = 2 on a 50 subject MR data set. The following shows the two templates (modes) discovered by the algorithm.


Two templates in a 50 subject MRI.

'Software'

The algorithm is currently implemented in the Insight ToolKit (ITK) and will be made publicly available. We also plan to integrate it into Slicer.

Key Investigators

MIT Algorithms: Mert R. Sabuncu and Polina Golland

Harvard DBP 2: M.E. Shenton, M. Kubicki and S. Bouix

Publications

In Print

Retrieved from "https://www.na-mic.org/w/index.php?title=Projects:MultimodalAtlas&oldid=24680"