250 likes | 402 Views
Unified Framework for Automatic Segmentation, Probabilistic Atlas Construction, Registration and Clustering of Brain MR Images. Annemie Ribbens Jeroen Hermans, Frederik Maes, Dirk Vandermeulen, Paul Suetens. Introduction. Computer–aided diagnosis. Introduction. Segmentation. Introduction.
E N D
Unified Framework for Automatic Segmentation, Probabilistic Atlas Construction, Registration and Clustering of Brain MR Images Annemie Ribbens Jeroen Hermans, Frederik Maes, Dirk Vandermeulen, Paul Suetens
Introduction Computer–aided diagnosis
Introduction Segmentation
Introduction Φ Atlas & Atlas-to-image registration
Introduction Population Specific Atlases
Introduction Images I Registrations Atlas previous iteration Atlas Construction
Images I Deformed images New Atlas Averaging
Introduction Computer aided-diagnosis Segmentation Prob. Atlases Registration
Framework • Aspects: • Segmentation • Clustering (i.e. computer-aided diagnosis) (+ Localization of cluster specific morphological differences) • Groupwise registration (nonrigid probabilistic atlasesper cluster) • Atlas-to-image registration • Advantages: • Less prior information necessary • Cooperation • Statistical framework convergence
Framework Segmentation Atlas-to-image registration Atlas formation & Clustering
Framework: model K = tissue classes number of Gaussians Y = intensities Image i
Framework: model Atlas t (Gray matter map) Image i (Gray matter map)
Framework: model Uniform prior for all voxels in an image
Framework: model G1 Deformations G2
Framework: MAP MAP: Jensen’s inequality Expectation maximization framework
Framework: EM algorithm: E-step Gaussian mixture model Uniform prior on the cluster memberships i = images j = voxels k = tissue classes t = clusters Gaussian prior on the deformations of each cluster Per cluster: atlas deformed towards image
Framework: EM Posterior Posterior = (clustering) * (segmentation using the atlas of the same cluster) Clustering = probability that voxel j of image i belongs to cluster t = sum over all tissue classes of the posterior = (prior of clustering) * (atlas is sharp & close to intensity model) * (subject specific registration close to groupwise) Segmentation= probability that a voxel belongs to a certain tissue class = sum over all clusters of the posterior = weighted sum of the segmentations using a specific atlas
Framework: EM algorithm: M-step • Maximum likelihood • Q-function parameters • All solutions close form (except registration) • Solutions (e.g. atlas) ~ literature
Framework: EM algorithm: M-step • Atlas • Prior cluster memberships • Groupwise registration • Atlas-to-image registration • No closed form solution • Spatial regularization • Viscous fluid model on derivative Weighting terms per voxel Gaussian mixture parameters:
w8 w8 w1 w1
Experiments Dice = • Brainweb data: • 20 simulated normal images • One cluster Segmentation & Atlas:
Experiments 8 brain MR images of healthy persons (normals) 8 brain MR images of Huntington disease patients (HD) Cluster memberships: all correctly classified
Conclusion • Statistical framework combining: • Segmentation • Clustering • Atlas construction per cluster (weighted) • Registration • Convergence & cooperation & less prior information needed • Validation promising • Cluster specific morphological differences are found • Easily extendable to incorporate clinical/spatial prior knowledge