Anatomical Image Analysis John Ashburner Functional Imaging Lab, 12 Queen Square, London, UK. A Growing Trend Larger and more complex models are being produced to explain brain imaging data. Bigger and better computers
allow more powerful models to be used More experience among software developers Older and wiser More engineers - rather than e.g. psychiatrists & biochemists This presentation is about combining various preprocessing procedures for anatomical images into a single generative model.
Traditional View of Pre-processing Brain image processing is often thought of as a pipeline procedure. One tool applied before another etc... For example Original Image Skull Strip
Classify Brain Tissues Non-uniformity Correct Extract Brain Surfaces Another example (for VBM) Bias Correction helps Registration
MRI images are corrupted by a smooth intensity non-uniformity (bias). Image intensity non-uniformity artefact has a negative impact on most registration approaches. Much better if this artefact is corrected.
Image with bias artefact Corrected image Bias Correction helps Segmentation Similar tissues no longer have similar intensities.
Artefact should be corrected to enable intensity-based tissue classification. Registration helps Segmentation SPM99 and SPM2 require tissue probability maps to be overlaid prior to segmentation. Segmentation helps Bias Correction
Bias correction should not eliminate differences between tissue classes. Can be done by make all white matter about the same intensity make all grey matter about the same intensity etc Currently fairly standard practice to combine bias correction and tissue classification
Segmentation helps Registration A convoluted method using SPM2 Template Spatially Normalised MRI Original MRI Affine register Affine Transform
Tissue probability maps Segmen t Grey Matter Spatial Normalisation - estimation
Spatial Normalisation - writing Deformation Unified Segmentation The solution to this circularity is to put everything in the same generative model. A
solution is found by repeatedly alternating among classification, bias correction and registration steps. The generative model involves: Mixture of Gaussians (MOG) Bias Correction Component
Warping (Non-linear Registration) Component Gaussian Probability Density If intensities are assumed to be Gaussian of mean k and variance 2k, then the probability of a value yi is: Non-Gaussian Probability Distribution
A non-Gaussian probability density function can be modelled by a Mixture of Gaussians (MOG): Mixing proportion - positive and sums to one Belonging Probabilities Belonging probabilities are assigned by normalising
to one. Mixing Proportions The mixing proportion k represents the prior probability of a voxel being drawn from class k - irrespective of its intensity. So: Non-Gaussian Intensity
Distributions Multiple Gaussians per tissue class allow nonGaussian intensity distributions to be modelled. Probability of Whole Dataset If the voxels are assumed to be independent, then the probability of the whole image is the product of the probabilities of each voxel: It
is often easier to work with negative logprobabilities: Modelling a Bias Field A bias field is included, such that the required scaling at voxel i, parameterised by , is i(). Replace the means by / () Replace the variances by ( / ()) k i
k i 2 Modelling a Bias Field After rearranging: y ()
y () Tissue Probability Maps Tissue probability maps (TPMs) are used instead of the proportion of voxels in each Gaussian as the prior. ICBM Tissue Probabilistic Atlases.
These tissue probability maps are kindly provided by the International Consortium for Brain Mapping, John C. Mazziotta and Arthur W. Toga. Mixing Proportions Tissue probability maps for each class are included. The probability of obtaining class k at voxel i, given weights is then: Deforming the Tissue Probability
Maps Tissue probability images are deformed according to parameters . The probability of obtaining class k at voxel i, given weights and parameters is then:
The Extended Model By combining the modified P(ci=k|) and P(yi|ci=k,), the overall objective function (E) becomes: The Objective Function Optimisation The best parameters are those that
minimise this objective function. Optimisation involves finding them. Begin with starting estimates, and repeatedly change them so that the objective function decreases each time. Steepest Descent Start Optimum Alternate between optimising different
groups of parameters Schematic of optimisation Repeat until convergence... Hold , , 2 and constant, and minimise E w.r.t. - Levenberg-Marquardt strategy, using dE/d and d2E/d2 Hold , , 2 and constant, and minimise E w.r.t. - Levenberg-Marquardt strategy, using dE/d and d2E/d2 Hold and constant, and minimise E w.r.t. , and 2 -Use an Expectation Maximisation (EM) strategy.
end Levenberg-Marquardt Optimisation LM optimisation is used for nonlinear registration () and bias correction (). Requires first and second derivatives of the objective function (E). Parameters and are updated by Increase
to improve stability (at expense of decreasing speed of convergence). Expectation Maximisation is used to update , 2 and For iteration (n), alternate between: E-step: Estimate belonging probabilities by: M-step: Set (n+1))
to values that reduce: Linear Regularisation Some bias fields and distortions are more probable (a priori) than others. Encoded using Bayes rule: Prior probability distributions can be modelled by a multivariate normal distribution. Mean vector and
Covariance matrix and -log[P()] = (- ( + const T
-1) Initial Affine Registration The procedure begins with
a Mutual Information affine registration of the image with the tissue probability maps. MI is computed from a 4x256 joint probability histogram. Joint Probability Histogram See D'Agostino, Maes, Vandermeulen & P. Suetens. Non-rigid Atlas-toImage Registration by Minimization of Class-Conditional Image Entropy. Proc. MICCAI 2004. LNCS 321)6, 2004.
Background voxels Pages 745-753. excluded Background Voxels are Excluded An intensity threshold is found by fitting image intensities to a mixture of two Gaussians. This threshold is used to exclude most of the voxels containing only air. Spatially
normalised BrainWeb phantoms (T1), T2 and PD) Tissue probability maps of GM and WM Cocosco, Kollokian, Kwan & Evans. BrainWeb: Online Interface to a 3D MRI Simulated Brain Database. NeuroImage 5(4):S425 (1)997)
Further Reading Ashburner & Friston. Unified Segmentation. To appear in NeuroImage. SPM Web Pages Look out for SPM5 http://www.fil.ion.ucl.ac.uk/spm/ Koen Van Leemputs page
contains some nice slides on tissue classification http://users.tkk.fi/~vanleemp/ A View of Science Science is about building models that can make predictions about the world. If its not predictive, then its not science. Biological
sciences are messy and kind of fuzzy. Need to work probabilistically. The only consistent system for working with probabilities is Bayesian. Dutch Book arguments.
Bayes Rule y P(|y) P(y|) P() P(y)
P(,y) - the data - a theory, model, or set of parameters - probability probability probability probability probability of
of of of of given y (posterior probability) y given (likelihood) (prior probability) y (evidence) and y (joint probability)
WRITING WORKSHOP REVIEW OF WRITING TERMS & SCHAFFER FORMAT WRITING TERMS Attention-Grabber (aka the HOOK) Statement that grabs reader's attention 1st sentence in intro paragraph Quote, statistic, definition, short anecdote What should you NEVER use as an AG statement?
Enroll in a flexible spending account - a health care and/or dependent care account, under the FSAFEDS Program. Employees MUST re-enroll in FSAFEDS each year to participate. Enroll in, change, or cancel an existing enrollment in a dental and/or a...
Thats Johns book; yours is on the back table with Carlos things. Jesus and Jerrys project was the last submitted in 15. I couldve sworn that there were two ss in buses, but my sister says its not true. "Wash...
Political cartooning is an art form used to express opinion. These opinions could be about anything, but normally address political, economic or social issues. Political cartoonists use certain techniques (or tools) to get their point across. They don't use all...
Carl Bender Physics Department Washington University in collaboration with Daniel Hook Theoretical Physics Imperial College Pup Tent Pup Tent Probability contours in the complex plane More interesting contours... First excited state - one node Second excited state - two nodes...
Library holdings can be described at a very granular level, including details about each library's copy of a work. In a context like Passage, how much of that granularity is appropriate? Is there still a need for a local system...
Aryan India was unified by a common Hindu culture but the arrival of Islamic influences led India to become divided into Hindu and Muslim areas. Aztecs and Yaquis of Mexico - it is believed that they used to share the...