averaged) 2. This makes it necessary to employ higher level knowledge about distribution of grey/white matter in the brain and lesion features.

Size: px
Start display at page:

Download "averaged) 2. This makes it necessary to employ higher level knowledge about distribution of grey/white matter in the brain and lesion features."

Transcription

1 Fusing Markov Random Fields with Anatomical Knowledge and Shape Based Analysis to Segment Multiple Sclerosis White Matter Lesions in Magnetic Resonance Images of the Brain Stephan Al-Zubi a, c, Klaus Toennies a, c, Nils Bodammer b, d, Hermann Hinrichs b, e a Institute for Simulation and Graphics, b Department of Neurology II Otto-von-Guericke University of Magdeburg c {stephan, klaus}@isg.cs.uni-magdeburg.de, d bodammer@neuro.med.uni-magdeburg.de e hermann.hinrichs@medizin.uni-magdeburg.de ABSTRACT This paper proposes an image analysis system to segment multiple sclerosis lesions of magnetic resonance (MR) brain volumes consisting of 3 mm thick slices using three channels (images showing T-, T- and PD -weighted contrast). The method uses the statistical model of Markov Random Fields (MRF) both at low and high levels. The neighborhood system used in this MRF is defined in three types: () Voxel to voxel: a low-level heterogeneous neighborhood system is used to restore noisy images. () Voxel to segment: a fuzzy atlas, which indicates the probability distribution of each tissue type in the brain, is registered elastically with the MRF. It is used by the MRF as a-priori knowledge to correct miss-classified voxels. (3) Segment to segment: Remaining lesion candidates are processed by a feature based classifier that looks at unary and neighborhood information to eliminate more false positives. An expert s manual segmentation was compared with the algorithm. Keywords: Markov Random Field, Registration, Multiple Sclerosis, Classification, Shape based analysis.. INTRODUCTION Multiple Sclerosis is a disease of the central nervous system. It appears as the myelin sheathes protecting nerve axons break down causing plaques. As a result 90% - 95% of lesions occur within white matter tissue. Certain parts of the nervous system are more affected by these MS-plaques: The anterior angles of the lateral ventricles, the corpus callosum and periventriclar areas. Lesion sizes vary considerably in size from a few millimeters to lesions involving the entire Centrum semiovale. MSplaques in white matter are ellipsoidal in shape. Their intensity is characterized by a bright center with decreasing intensity as we move away from the center. Figure. Left: example of a slice showing PD contrast. Middle: lesion mask obtained by intensity based classification alone and the actual lesion mask. Right: Manual segmentation. As seen much grey matter is falsely classified as lesions. Segmentation based on intensity alone will not succeed because the intensity histograms of grey matter and lesions overlap in MRI images (See Figure ). This means that any classifier based on voxel intensity alone will classify some gray matter as lesions and lesions as gray matter. This effect is evident where partial volume effect exists (e.g. where cerebrospinal fluid (CSF) and grey matter are

2 averaged). This makes it necessary to employ higher level knowledge about distribution of grey/white matter in the brain and lesion features. There are a number of different approaches that employ different models for segmentation. Some researchers proposed matching a brain volume with an anatomical atlas to create grey/white mask constraints used to label lesions correctly. Warfield proposes segmenting the cortex by region growing and constraining its boundary by an elastically registered anatomical atlas. This is used to make the white matter mask containing lesions and white matter that can be separated by a simple intensity based classifier. In Warfield 4 a statistical k-nn classifier for gray levels is used to classify tissues based on a set of manually selected prototypes. A spatial mask obtained by elastic matching adds spatial features to the k-nn classifier resolving the ambiguity of overlapping histograms. The elastic model is re-aligned with the volume using the improved segmentation. The algorithm iteratively applies k-nn with elastic matching until convergence. Similar work can be found for segmenting tumors in Kaus et al. 3 and Kamber et al. 5. They propose the use of a brain tissue probability model. This model gives the probability distribution of each tissue type within a standardized Talairach space. A set of segmented brain volumes were registered into Talairach space and averaged together to construct the model. This model was used as a mask to confine search within white matter and as a way to provide geometric features used to classify lesions (e.g. probability of ventricular CSF is useful in detecting periventriclar lesions). Another way to classify lesions is using a feature space collected from possible candidate lesions to sort out false positives. Ardizzone 6 uses the fuzzy c-means algorithm by first obtaining a set of oversegmented regions followed by a re-clustering phase. The re-clustering uses shape and intensity features to label or split unknown clusters. After this phase three masks corresponding to WM, WM+GM and WM+GM+CSF are built. The holes in these masks correspond to candidate lesions. A feature vector for each candidate is passed through a neural classifier to decide whether it s a lesion or not. Those features include: Contact with WM and GM and CSF, mean intensity, shape measures like compactness and elongation and position of the lesions like the distance from the ventricular area. Fuzzy connectedness 7 is used to detect connected fuzzy objects representing white, grey and CSF objects. The holes remaining represent locations of potential lesions. False positives are eliminated by operator assistance. Fuzzy connectedness is defined by a measure that computes a connected component such that the strength of connection between any two pixels in the segment is above a certain threshold. The strength of connection between pixels is determined by the distance and intensity homogeneity. Johnson 8 applies Iterated Conditional Mode (ICM) algorithm used to segment noisy images with morphological operators. This method consists of initially segmenting the volume by ICM using all labels. The segmentation is post processed to get the white matter + lesion mask. Partial volume results in classifying the areas around lesions as grey matter. These gaps in the mask are closed by binary morphology. The resulting mask is processed again with ICM but this time classifying only two tissue types. Small lesions are eliminated from the final output. Kirshnan 9 uses the ratio image of PD/T to eliminate magnetic field inhomogeneity. He then applies a series of thresholds to extract the WM, CSF mask. The lesions are identified by intensity thresholding as the brightest spots and then applying those masks to remove false lesions outside white matter. To summarize, many authors aim to find the white matter mask containing the lesions and then segmenting the lesions afterwards. Some apply registration with an anatomical model to determine where grey matter is. Others apply feature based classifiers to eliminate false positive. The method represented in this paper combines alignment with an anatomical model, feature based classification and image restoration using Markov Random Fields. This paper is structured as follows: section () introduces the Markov Random Fields which will be used to restore noisy data at a voxel level and later to eliminate false positives for candidate lesion segments. Section (3) will describe the process pipeline to classify MS lesions and surveys some literature pertaining to each step. Section (4) shows the experimental results followed by discussion of future directions.

3 . MARKOV RANDOM FIELDS (MRF) Let Y= {y i,k,m i N x, j N y, k N z, m M} be a multi-modal input image of dimension N x N y N z and M channels. Let X= {x i,k,m : i N x, j N y, k N z } be the corresponding segmentation of Y with each x i,k,m taking a value form a set L={, l} of labels. Let x Ω be a specific segmentation of X. Ω is the space of all segmentations. Let x be an instance of X and similarly y an instance of Y. The objective is to maximize the probability of P(X Y). Given Bayes Theorem P(X Y) P(Y X)P(X) we get x * = arg max P( y x) P( x) () x Ω The term P(y x) in () is modeled as a multivariate Gaussian additive noise which is independent of location. Let µ p, p L be the vector representing the mean intensity values for a label (p) and Σ p be the corresponding covariance matrix. Then we can define P(y x) as l ( ) ( ) Σ x exp y i k x x ( i k x ) i j k, j, µ Σ y, i, j, k i, j, k, j, i, j, k P( x) = µ, i, k y π () The advantage of using the multi-variate normal () is that it can better compensate for correlated tissue variations due to partial volume effect 6. P(x) in () is the regularization term in which each voxel interacts with its neighbors to restore noisy images by inferring the correct value from neighbors. The homogeneous multi-level logistic (MLL) distribution will be used to define P(x). Specifically, for each voxel a first order neighborhood system consisting of its six adjacent voxels is defined: { x, x x } η (3) x =,,,,, i, k i± j k i j± k i, j, k± A potential function is defined on cliques of one voxel and cliques of adjacent voxel pairs, specifically, given two adjacent sites u, u : β u = u V ( u, u ) = (4) β otherwise And for each label type (k) we define an a-priori expectation α k and so that a single site potential is Using (4) and (5), P(x) is defined as V ( x i, k ) α x i, k Z is the normalizing constant called the partition function. = (5) exp V ( u, u ) + V ( u) adjacent u u u x (, ) P(x ) = (6) Z The local optimization method of iterated conditional modes ICM by Besag 0 works by iteratively updating one site at a time in raster scanning fashion until label changes per image iteration falls below a threshold. In the case of optimizing a single site (x i,k =u) we only have to maximize () for a single site which simplifies as the minimization of: U ( ) = y i, k µ u Σu ( y i, k µ u ) + ln( Σu ) V ( u, u ) α u (7) ( u) u η u

4 3. METHOD The algorithm pipeline goes through the following steps (See Figure ): () Initial gray level segmentation. () Heterogeneous clique computation for the MRF. (3) Magnetic field inhomogeneity correction. (4) Image restoration by the ICM algorithm. (5) Registration with a tissue probability distribution atlas. (6) Using the aligned atlas to correct misclassified voxels. (7) Eliminating more false positive lesions by applying an MRF classifier on lesion segments. It looks at shape features of those segments and their neighborhood. The following are the individual steps in detail: initial T, T, PD data () initial gray level segmentation () compute heterogeneous clique weights (3) magnetic field correction (4) ICM (6) apply atlas on areas marked as lesions and regions that surround them using ICM (5) register atlas with segmented data (7) apply shape based ICM to eliminate false positives (8) add resulting segmentation as an instance to the atlas 3.. Initial Grey Level Segmentation Figure. Algorithm pipeline The first step in the pipeline is to perform use a multivariate Gaussian grey level classifier where the sample mean and standard deviation are acquired by manual segmentations of some brain slices (see Figure ). 3.. Compute Heterogeneous Clique Weights The interaction between neighbors can be further improved if we include edge information into P(x). The voxel-voxel weight β in V (see eq. 4) is redefined as a local field {β u,v }between voxels u and v to better adjust to likely edges by a function similar to eq. (3): β β u, v = (8) + ( u ) / λ Where u is the direction between u and v. and y σ is the gradient of a Gaussian smoothed image. y This method represents a low cost approximation to the heterogeneous field which enables a better behavior for the ICM algorithm (See Figure 3). Many authors include edge information into the Markov Random Field, 3, 4. Equation (8) is an adaptation of anisotropic diffusion filters found in Weickert Intensity Inhomogeneity Correction The image Y is intensity corrected for magnetic field inhomogeneities by using the white matter mask found in step () and modeling the inhomogeneity as a multiplicative field (y corrected (i,k)=y i,k *m i,k ). A grid of 3D points is defined where at each point a local estimate of the mean gray value of white matter is calculated from a minimum neighborhood size of Q pixels. The multiplicative field is defined as the ratio between the local and the global average (m i,k =µ global / µ local (i,k)). Bilinear interpolation provides the estimation for the remaining voxels. This approach is similar to the one used in Rajapakse et al.. σ

5 Figure 3. Heterogeneous MRF neighborhood weights in X, Y, Z directions 3.4. Applying the ICM Algorithm to Restore Miss-classified Voxels The image is restored using the method described in section () Registering Atlas with Segmented Data X In this paper the atlas P atlas ={ p i,k (l) } is a probability distribution of segments where p i,k (l) represents the probability that the voxel i,k has a label l. The transformation (T =T affine +T B-Spline ) is used to warp the atlas to the input image y or its segmentation x where: T affine: Represents the global extrinsic transform of 9 parameters defined by the affine space. T B-Spline : B-Splines were used in to model elastic intrinsic deformation because of their: () built in smoothness () locality of control that enables fast registration (3) a small parameter space that includes only the xyz distances between control points (for more details see Rueckert et al. 7 and Studholme et al. 9 ). The normalized entropy measure 8 {H(A,B)=(H(A)+H(B))/H(A,B)} enables the registration between x and y directly because of its ability to compare images of different modalities. A better similarity measure between P atlas and x is defined by S( P atlas, x) = pi, k ( xi, k ) (9) n i, k It is possible to build an atlas by averaging several brain volumes by registration to build a probability distribution model for example see Chen et al. 0 and Studhome et al. 9. In this paper an initial atlas was obtained by segmenting a simulated volume in brainweb then blurring the distribution by a Gaussian kernel. The blurring improves registration by minimizing chances of local minima in the fitting process. As more segmented images are obtained we can refine the atlas by registering the segmentations with P atlas and adding the new instance to form a new distribution. n+ npatlas + x P atlas = (0) n + The atlas has the ability to include the a-priori probability distributions of lesions in the brain where pi,k(l lesion ) represents the expectancy of a lesion at location i,k. Initially we set p to some threshold (αpi,k(l white-matter )) which gets refined as more segmentation instances are added Lesion Detection Using the registered atlas in section 3.5, we can differentiate between true and false lesions. This is done by reapplying ICM at segments labeled as lesions and their neighboring voxels with an additional term in equation (7) using P atlas. n

6 Figure 4. Improvement after applying the anatomical atlas. Left: the lesions of Figure. Middle: lesions reclassified by applying the anatomical atlas. Right: the expert s manual segmentation. As we see the atlas was able to eliminate a great deal of false positives. The remaining false positives are handled by shape based MRF. ( ) U ( u) = y i, k µ u Σu ( y i, k µ u ) + ln( Σu ) V ( u, u ) u + ln( pi, k ( u)) α γ () This has the effect of adding special information to each voxel resulting in a better classification (see Figure 4). Registration errors may still misclassify certain voxels therefore feature based classification is needed to improve results Applying Shape Based MRF False positives resulting from section 3.6 are further eliminated by defining individual candidate lesions as shape units in a Markov Random Field. Each shape unit can be reassigned a state {MS lesion, grey matter} based on shape and neighborhood features. An iterative process similar to ICM is used to re-label those units until convergence. The following is the description of the steps in detail:. The watershed transform is used to isolate most individual plaques in each XY slice using the property that generally the intensity is bright at the center and decreases towards the edges. The connected components thus formed within white matter generally have an elliptical shape (unless they contact other structures like the ventricular system). Per slice connected components are used instead of 3D connected components because the Z aspect ratio of (3:) makes calculation of good shape features in 3D difficult (i.e. Lesions generally do not persist more than -4 slices and partial volume makes them stick together in the Z direction).. The shape units calculated in step () are considered random variables (sites) which may be relabeled as {lesion, grey matter}: S= {s s k }, L S : S {l grey_matter, l lesion } 3. Unary features are computed for each shape unit that measure () compactness () elongation (3) distance from the center of ventricular area (φ unary : S v unary ). 4. A neighborhood system (η S z : S S ) is defined for each shape unit where η S z (s i ) are all other shape units that have contact with s i in the Z direction. Similarly we define η S xy as the neighborhood that specifies contact in the XY direction. Binary features are computed for each shape unit that includes: () the area of contact with gray matter () the area of contact with other shape units in the neighborhood system {η S z, η S xy } (φ binary : S v binary ). 5. A binary classifier C: S η S z η S xy (P lesion, P grey_matter ) is defined. It assigns higher probability that a shape unit is a lesion whenever it is more oval or contacts more other lesions or has a certain location. It assigns it a higher probability that it is grey matter when it contacts more grey matter or has a certain location. u η u

7 Figure 5. Left: Lesion mask after applying the anatomical atlas. Middle: Lesion mask after applying shape based MRF. Right: Experts manual segmentation 6. Using the neighborhood system and shape features, an iterated conditional mode algorithm is defined where each site is assigned a state of lesion or gray matter based the classifier C. The Algorithm is iterated a few times until convergence (see Figure 5). The algorithm uses contact information between shape units to mimic what the expert does. For example, when the expert finds an isolated oval region in white matter and looking at a the slice before he discovers that this structure is connected to cortical gray matter so it must be a continuation of cortical matter, however, if it contacts another lesion then it is likely to be a continuation of that lesion structure. Another example are elongated structures embedded in gray matter are less likely to be lesions but classification errors. 4. RESULTS Experiments were conducted on three brain volumes with T, T, and PD weighted contrast. Each modality is of dimension where the voxel size is mm. Manual segmentations of the multiple-sclerosis lesions for the three volumes were obtained by a medical expert. Two measures were used to compare the manual and automatic segmentation:. Similarity index 8 : The similarity index between two segments A, A is a number between [0, ] defined by S = A I A A + A (). Automatic and Manual segmentation volumes which show the quantity of false positives eliminated after each stage of the algorithm. Figure 6 shows similarity index for a typical MS patient after each stage of the algorithm. The ICM algorithm restored the image and yielded a 5% improvement in similarity index for the whole volume. Using the atlas yields a 3% improvement in similarity. Using the shape based ICM yields a 3% improvement in similarity especially after slice -6 as shown in the figure. Figure 7 shows automatic versus the manual segmentation volumes for the same patient. Initially the overall correlation between manual and intensity based segmentation is 0.3. The correlation became 0.87 after applying the atlas because most false positives were eliminated at this stage. After applying the shape based ICM the correlation became 0.95 as can be seen between the closely matching curves. Experiments with the other two patients yielded a per-slice-correlation of 0.93 and 0.9 between the volumes of manual and automatic segmentation.

8 Intensity based Segmentation ICM segmentation Atlas segmentation Shape ICM Similarity Index Slice Number Figure 6. Similarity Index for after each stage of the algorithm pipeline for one patient volume 3500 Manual Segmentation Intensity based Segmentation ICM segementation Atlas segmentation Shape ICM Volume on MS lesions Slice Number Figure 7. Manual segmentation volume vs. Automatic segmentation volume for each stage of the algorithm. 5. CONCLUSION AND FUTURE WORK A method to segment multiple sclerosis was presented. It employs a three stage algorithm which first segments and restores the image at low level and then an anatomical atlas is used to disambiguate between lesions and gray matter. A shape based MRF in which shape units consisting of lesion slices and a neighborhood system representing contact with other lesions and gray matter is used to eliminate false positives. Future work will concentrate on refining the shape based classifier. For example: () defining a spatial probability map for contact information. () Improving the feature space to include more lesion properties (3) Improving of classification rules to be more close to that of the expert.

9 6. REFERENCES. C. Poster, An Atlas of Multiple Sclerosis, The Parthenon Publishing Group, S. Warfield, J. Dengler, J. Zaers, R. Guttman, W. Wells, G. Ettinger, J. Hiller, R. Kikinis, Automatic identification of grey matter structures from MRI to improve the segmentation of white matter lesions. J Image Guid Surg, pp , M. Kaus, S. Warfield, F. Jolesz, R. Kikinis, Adaptive Template Moderated Brain Tumor Segmentation in MRI. Bildverbeitung für die Medizin, 0-06, S. Warfield, M. Kaus, F. Jolesz, R. Kikinis, Adaptive, Template Mpderated, Spatially varying Statistical Classification, Medical Image Analysis 4(), pp , M. Kamber, R. Shinghal, L. Collins, G. Francis, A. Evans, Model-Based 3D Segmentation of Multiple Sclerosis lesions in Magnetic Resonance Brain Images, IEEE Trans. On medical Images 4(3), pp , E. Ardizzone, R. Pirrone, An Architecture for the Recognition and Classification of Multiple Sclerosis Lesions in MR Images, IDAMAP, J. Udupa, L. Wei, S. Samarasekera, Y. Miki, M. Vanbuchem, R. Grossman, Multiple sclerosis lesion quantification using fuzzy- connectedness principles, IEEE Trans. on medical imaging 6(), pp , B. Johnston, M. Atkins, B. Mackiewich, M. Anderson, Segmentation of Multiple Sclerosis Lesions in Intensity Corrected Multispectral MRI, IEEE Trans. On Medical Imaging 5(), K. Kirshnan, M. Atkins, Segmenttion of Multiple Sclerosis in MRI- An Image Analysis Approach, SPIE Conf. On Image Processing, pp. 06-6, San Diego, J. Besag, On the statistical analysis of dirty pictures, J. Roy. Statist. Soc. B 48(3), pp , S. Lakshmanan, D. Haluk, Simultaneous Parameter Estimation and Segmentation of Gibbs Random Fields Using Simulated Annealing, IEEE Trans. On Pattern Analysis and Machine Intelligence (8), pp , S. Geman, D. Geman, Stochastic Relaxtation, Gibbs Distributions, and the Bayesian Restoration of Images, IEEE Trans. On Pattern Analysis and Machine Intelligence 6(6), S. Nadabar, A. Jain, Parameter Estimation in Markov Random Field Contextual Models Using Geometric Models of Objects, IEEE Trans. On Pattern Analysis and machine Intelligence 8(3), R. Aykroyd, Bayesian Estimation for Homogeneous and Inhomogeneous Gaussian Random Fields, IEEE Trans. On Pattern Analysis and machine Intelligence 0(5), J. Weickert, Anisotropic Diffusion in Image Processing, B.G Teubner Stuttgart, M. Desco, J. Gispert, S. Reig, A. Santos, J. Pascau, N. Malpica, P. Garcia-Barreno, Statistical Segmentation of multidimensional brain datasets, SPIE Medical Imaging, pp 84-93, D. Rueckert, L. Sonoda, C. Hayes, D. Hill, M. Leach, D. Hawks, Nonrigid Registration Using Free-Form Deformations: Application to Breast MR Images, IEEE Trans. on Medical Imaging 8(8), pp 7-7, C. Studholme, D. Hill, D. Hawks, An overlap invariant entropy measure of 3D medical image alignment, Pattern Recognition 3, pp 7-86, C. Studholme, V. Cardenas, M. Weiner, Multi scale image and multi scale deformation of brain anatomy for building average brain atlases, SPIE medical Imaging, pp , M. Chen, T. Kanade, D. Pomerleau, J. Schneider, Probabilistic Registration of 3-D Medical Images, tech. report CMU-RI-TR-99-6, Robotics Institute, Carnegie Mellon University, J. Rajapakse, J. giedd, J. Rapoport, Statistical Approach to Segmentation of Single-Channel Cerebral MR Images, IEEE Trans. on medical imaging 6( ), pp , Brainweb 3. C.A. Cocosco, V. Kollokian, R.K.-S. Kwan, A.C. Evans : "BrainWeb: Online Interface to a 3D MRI Simulated Brain Database", NeuroImage 5(4), part /4, S45, Proceedings of 3-rd International Conference on Functional Mapping of the Human Brain, Copenhagen, 997.

MR IMAGE SEGMENTATION

MR IMAGE SEGMENTATION MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification

More information

An Introduction To Automatic Tissue Classification Of Brain MRI. Colm Elliott Mar 2014

An Introduction To Automatic Tissue Classification Of Brain MRI. Colm Elliott Mar 2014 An Introduction To Automatic Tissue Classification Of Brain MRI Colm Elliott Mar 2014 Tissue Classification Tissue classification is part of many processing pipelines. We often want to classify each voxel

More information

Norbert Schuff VA Medical Center and UCSF

Norbert Schuff VA Medical Center and UCSF Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role

More information

Automatic segmentation of the cortical grey and white matter in MRI using a Region Growing approach based on anatomical knowledge

Automatic segmentation of the cortical grey and white matter in MRI using a Region Growing approach based on anatomical knowledge Automatic segmentation of the cortical grey and white matter in MRI using a Region Growing approach based on anatomical knowledge Christian Wasserthal 1, Karin Engel 1, Karsten Rink 1 und André Brechmann

More information

Preprocessing II: Between Subjects John Ashburner

Preprocessing II: Between Subjects John Ashburner Preprocessing II: Between Subjects John Ashburner Pre-processing Overview Statistics or whatever fmri time-series Anatomical MRI Template Smoothed Estimate Spatial Norm Motion Correct Smooth Coregister

More information

ADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION

ADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION ADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION Abstract: MIP Project Report Spring 2013 Gaurav Mittal 201232644 This is a detailed report about the course project, which was to implement

More information

Histograms. h(r k ) = n k. p(r k )= n k /NM. Histogram: number of times intensity level rk appears in the image

Histograms. h(r k ) = n k. p(r k )= n k /NM. Histogram: number of times intensity level rk appears in the image Histograms h(r k ) = n k Histogram: number of times intensity level rk appears in the image p(r k )= n k /NM normalized histogram also a probability of occurence 1 Histogram of Image Intensities Create

More information

Automatic Registration-Based Segmentation for Neonatal Brains Using ANTs and Atropos

Automatic Registration-Based Segmentation for Neonatal Brains Using ANTs and Atropos Automatic Registration-Based Segmentation for Neonatal Brains Using ANTs and Atropos Jue Wu and Brian Avants Penn Image Computing and Science Lab, University of Pennsylvania, Philadelphia, USA Abstract.

More information

Computational Neuroanatomy

Computational Neuroanatomy Computational Neuroanatomy John Ashburner john@fil.ion.ucl.ac.uk Smoothing Motion Correction Between Modality Co-registration Spatial Normalisation Segmentation Morphometry Overview fmri time-series kernel

More information

Automatic Generation of Training Data for Brain Tissue Classification from MRI

Automatic Generation of Training Data for Brain Tissue Classification from MRI MICCAI-2002 1 Automatic Generation of Training Data for Brain Tissue Classification from MRI Chris A. Cocosco, Alex P. Zijdenbos, and Alan C. Evans McConnell Brain Imaging Centre, Montreal Neurological

More information

Norbert Schuff Professor of Radiology VA Medical Center and UCSF

Norbert Schuff Professor of Radiology VA Medical Center and UCSF Norbert Schuff Professor of Radiology Medical Center and UCSF Norbert.schuff@ucsf.edu 2010, N.Schuff Slide 1/67 Overview Definitions Role of Segmentation Segmentation methods Intensity based Shape based

More information

Probabilistic Registration of 3-D Medical Images

Probabilistic Registration of 3-D Medical Images Probabilistic Registration of 3-D Medical Images Mei Chen, Takeo Kanade, Dean Pomerleau, Jeff Schneider CMU-RI-TR-99-16 The Robotics Institute Carnegie Mellon University Pittsburgh, Pennsylvania 1513 July,

More information

Automatic MS Lesion Segmentation by Outlier Detection and Information Theoretic Region Partitioning Release 0.00

Automatic MS Lesion Segmentation by Outlier Detection and Information Theoretic Region Partitioning Release 0.00 Automatic MS Lesion Segmentation by Outlier Detection and Information Theoretic Region Partitioning Release 0.00 Marcel Prastawa 1 and Guido Gerig 1 Abstract July 17, 2008 1 Scientific Computing and Imaging

More information

Introduction to Medical Image Processing

Introduction to Medical Image Processing Introduction to Medical Image Processing Δ Essential environments of a medical imaging system Subject Image Analysis Energy Imaging System Images Image Processing Feature Images Image processing may be

More information

Computer-Aided Tracking of MS lesions

Computer-Aided Tracking of MS lesions Computer-Aided Tracking of MS lesions Deborah Sturm* 1, Devorah Gurwitz Kletenik 2 Phillip Koshy 3 1 Department of Computer Science, College of Staten Island, City University of New York, 2800 Victory

More information

Norbert Schuff Professor of Radiology VA Medical Center and UCSF

Norbert Schuff Professor of Radiology VA Medical Center and UCSF Norbert Schuff Professor of Radiology Medical Center and UCSF Norbert.schuff@ucsf.edu Slide 1/67 Overview Definitions Role of Segmentation Segmentation methods Intensity based Shape based Texture based

More information

Non-rigid Image Registration using Electric Current Flow

Non-rigid Image Registration using Electric Current Flow Non-rigid Image Registration using Electric Current Flow Shu Liao, Max W. K. Law and Albert C. S. Chung Lo Kwee-Seong Medical Image Analysis Laboratory, Department of Computer Science and Engineering,

More information

Automatic Generation of Training Data for Brain Tissue Classification from MRI

Automatic Generation of Training Data for Brain Tissue Classification from MRI Automatic Generation of Training Data for Brain Tissue Classification from MRI Chris A. COCOSCO, Alex P. ZIJDENBOS, and Alan C. EVANS http://www.bic.mni.mcgill.ca/users/crisco/ McConnell Brain Imaging

More information

mritc: A Package for MRI Tissue Classification

mritc: A Package for MRI Tissue Classification mritc: A Package for MRI Tissue Classification Dai Feng 1 Luke Tierney 2 1 Merck Research Labratories 2 University of Iowa July 2010 Feng & Tierney (Merck & U of Iowa) MRI Tissue Classification July 2010

More information

Image Registration + Other Stuff

Image Registration + Other Stuff Image Registration + Other Stuff John Ashburner Pre-processing Overview fmri time-series Motion Correct Anatomical MRI Coregister m11 m 21 m 31 m12 m13 m14 m 22 m 23 m 24 m 32 m 33 m 34 1 Template Estimate

More information

Comparison Study of Clinical 3D MRI Brain Segmentation Evaluation

Comparison Study of Clinical 3D MRI Brain Segmentation Evaluation Comparison Study of Clinical 3D MRI Brain Segmentation Evaluation Ting Song 1, Elsa D. Angelini 2, Brett D. Mensh 3, Andrew Laine 1 1 Heffner Biomedical Imaging Laboratory Department of Biomedical Engineering,

More information

Spatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields

Spatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields Spatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields Wei Liu 1, Peihong Zhu 1, Jeffrey S. Anderson 2, Deborah Yurgelun-Todd 3, and P. Thomas Fletcher 1 1 Scientific

More information

Automatic Segmentation of Neonatal Brain MRI

Automatic Segmentation of Neonatal Brain MRI Automatic Segmentation of Neonatal Brain MRI 1 Marcel Prastawa, 2 John Gilmore, 3 Weili Lin, and 1,2 Guido Gerig Department of 1 Computer Science, 2 Psychiatry, and 3 Radiology University of North Carolina,

More information

Nonrigid Registration using Free-Form Deformations

Nonrigid Registration using Free-Form Deformations Nonrigid Registration using Free-Form Deformations Hongchang Peng April 20th Paper Presented: Rueckert et al., TMI 1999: Nonrigid registration using freeform deformations: Application to breast MR images

More information

Supplementary methods

Supplementary methods Supplementary methods This section provides additional technical details on the sample, the applied imaging and analysis steps and methods. Structural imaging Trained radiographers placed all participants

More information

The organization of the human cerebral cortex estimated by intrinsic functional connectivity

The organization of the human cerebral cortex estimated by intrinsic functional connectivity 1 The organization of the human cerebral cortex estimated by intrinsic functional connectivity Journal: Journal of Neurophysiology Author: B. T. Thomas Yeo, et al Link: https://www.ncbi.nlm.nih.gov/pubmed/21653723

More information

MRI Segmentation. MRI Bootcamp, 14 th of January J. Miguel Valverde

MRI Segmentation. MRI Bootcamp, 14 th of January J. Miguel Valverde MRI Segmentation MRI Bootcamp, 14 th of January 2019 Segmentation Segmentation Information Segmentation Algorithms Approach Types of Information Local 4 45 100 110 80 50 76 42 27 186 177 120 167 111 56

More information

Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit

Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit John Melonakos 1, Ramsey Al-Hakim 1, James Fallon 2 and Allen Tannenbaum 1 1 Georgia Institute of Technology, Atlanta GA 30332,

More information

1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12

1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12 Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....

More information

Correspondence Detection Using Wavelet-Based Attribute Vectors

Correspondence Detection Using Wavelet-Based Attribute Vectors Correspondence Detection Using Wavelet-Based Attribute Vectors Zhong Xue, Dinggang Shen, and Christos Davatzikos Section of Biomedical Image Analysis, Department of Radiology University of Pennsylvania,

More information

Neuroimaging and mathematical modelling Lesson 2: Voxel Based Morphometry

Neuroimaging and mathematical modelling Lesson 2: Voxel Based Morphometry Neuroimaging and mathematical modelling Lesson 2: Voxel Based Morphometry Nivedita Agarwal, MD Nivedita.agarwal@apss.tn.it Nivedita.agarwal@unitn.it Volume and surface morphometry Brain volume White matter

More information

Subvoxel Segmentation and Representation of Brain Cortex Using Fuzzy Clustering and Gradient Vector Diffusion

Subvoxel Segmentation and Representation of Brain Cortex Using Fuzzy Clustering and Gradient Vector Diffusion Subvoxel Segmentation and Representation of Brain Cortex Using Fuzzy Clustering and Gradient Vector Diffusion Ming-Ching Chang Xiaodong Tao GE Global Research Center {changm, taox} @ research.ge.com SPIE

More information

Region-based Segmentation

Region-based Segmentation Region-based Segmentation Image Segmentation Group similar components (such as, pixels in an image, image frames in a video) to obtain a compact representation. Applications: Finding tumors, veins, etc.

More information

Performance Evaluation of the TINA Medical Image Segmentation Algorithm on Brainweb Simulated Images

Performance Evaluation of the TINA Medical Image Segmentation Algorithm on Brainweb Simulated Images Tina Memo No. 2008-003 Internal Memo Performance Evaluation of the TINA Medical Image Segmentation Algorithm on Brainweb Simulated Images P. A. Bromiley Last updated 20 / 12 / 2007 Imaging Science and

More information

Image Segmentation and Registration

Image Segmentation and Registration Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation

More information

A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction

A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction Tina Memo No. 2007-003 Published in Proc. MIUA 2007 A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction P. A. Bromiley and N.A. Thacker Last updated 13 / 4 / 2007 Imaging Science and

More information

Adaptive Fuzzy Connectedness-Based Medical Image Segmentation

Adaptive Fuzzy Connectedness-Based Medical Image Segmentation Adaptive Fuzzy Connectedness-Based Medical Image Segmentation Amol Pednekar Ioannis A. Kakadiaris Uday Kurkure Visual Computing Lab, Dept. of Computer Science, Univ. of Houston, Houston, TX, USA apedneka@bayou.uh.edu

More information

Registration-Based Segmentation of Medical Images

Registration-Based Segmentation of Medical Images School of Computing National University of Singapore Graduate Research Paper Registration-Based Segmentation of Medical Images by Li Hao under guidance of A/Prof. Leow Wee Kheng July, 2006 Abstract Medical

More information

FROM IMAGE RECONSTRUCTION TO CONNECTIVITY ANALYSIS: A JOURNEY THROUGH THE BRAIN'S WIRING. Francesca Pizzorni Ferrarese

FROM IMAGE RECONSTRUCTION TO CONNECTIVITY ANALYSIS: A JOURNEY THROUGH THE BRAIN'S WIRING. Francesca Pizzorni Ferrarese FROM IMAGE RECONSTRUCTION TO CONNECTIVITY ANALYSIS: A JOURNEY THROUGH THE BRAIN'S WIRING Francesca Pizzorni Ferrarese Pipeline overview WM and GM Segmentation Registration Data reconstruction Tractography

More information

Where are we now? Structural MRI processing and analysis

Where are we now? Structural MRI processing and analysis Where are we now? Structural MRI processing and analysis Pierre-Louis Bazin bazin@cbs.mpg.de Leipzig, Germany Structural MRI processing: why bother? Just use the standards? SPM FreeSurfer FSL However:

More information

Image Processing with Nonparametric Neighborhood Statistics

Image Processing with Nonparametric Neighborhood Statistics Image Processing with Nonparametric Neighborhood Statistics Ross T. Whitaker Scientific Computing and Imaging Institute School of Computing University of Utah PhD Suyash P. Awate University of Pennsylvania,

More information

Validation of Image Segmentation and Expert Quality with an Expectation-Maximization Algorithm

Validation of Image Segmentation and Expert Quality with an Expectation-Maximization Algorithm Validation of Image Segmentation and Expert Quality with an Expectation-Maximization Algorithm Simon K. Warfield, Kelly H. Zou, and William M. Wells Computational Radiology Laboratory and Surgical Planning

More information

Learning-based Neuroimage Registration

Learning-based Neuroimage Registration Learning-based Neuroimage Registration Leonid Teverovskiy and Yanxi Liu 1 October 2004 CMU-CALD-04-108, CMU-RI-TR-04-59 School of Computer Science Carnegie Mellon University Pittsburgh, PA 15213 Abstract

More information

ECSE 626 Project Report Multimodality Image Registration by Maximization of Mutual Information

ECSE 626 Project Report Multimodality Image Registration by Maximization of Mutual Information ECSE 626 Project Report Multimodality Image Registration by Maximization of Mutual Information Emmanuel Piuze McGill University Montreal, Qc, Canada. epiuze@cim.mcgill.ca Abstract In 1997, Maes et al.

More information

Lesion Segmentation and Bias Correction in Breast Ultrasound B-mode Images Including Elastography Information

Lesion Segmentation and Bias Correction in Breast Ultrasound B-mode Images Including Elastography Information Lesion Segmentation and Bias Correction in Breast Ultrasound B-mode Images Including Elastography Information Gerard Pons a, Joan Martí a, Robert Martí a, Mariano Cabezas a, Andrew di Battista b, and J.

More information

Ulrik Söderström 16 Feb Image Processing. Segmentation

Ulrik Söderström 16 Feb Image Processing. Segmentation Ulrik Söderström ulrik.soderstrom@tfe.umu.se 16 Feb 2011 Image Processing Segmentation What is Image Segmentation? To be able to extract information from an image it is common to subdivide it into background

More information

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 4: Pre-Processing Medical Images (II)

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 4: Pre-Processing Medical Images (II) SPRING 2016 1 MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 4: Pre-Processing Medical Images (II) Dr. Ulas Bagci HEC 221, Center for Research in Computer Vision (CRCV), University of Central Florida (UCF),

More information

INTERPOLATION is a commonly used operation in image

INTERPOLATION is a commonly used operation in image IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 2, FEBRUARY 1999 137 A Task-Specific Evaluation of Three-Dimensional Image Interpolation Techniques George J. Grevera, Jayaram K. Udupa,* Senior Member,

More information

An Intensity Consistent Approach to the Cross Sectional Analysis of Deformation Tensor Derived Maps of Brain Shape

An Intensity Consistent Approach to the Cross Sectional Analysis of Deformation Tensor Derived Maps of Brain Shape An Intensity Consistent Approach to the Cross Sectional Analysis of Deformation Tensor Derived Maps of Brain Shape C. Studholme, V. Cardenas, A. Maudsley, and M. Weiner U.C.S.F., Dept of Radiology, VAMC

More information

Structural Segmentation

Structural Segmentation Structural Segmentation FAST tissue-type segmentation FIRST sub-cortical structure segmentation FSL-VBM voxelwise grey-matter density analysis SIENA atrophy analysis FAST FMRIB s Automated Segmentation

More information

Surface-based Analysis: Inter-subject Registration and Smoothing

Surface-based Analysis: Inter-subject Registration and Smoothing Surface-based Analysis: Inter-subject Registration and Smoothing Outline Exploratory Spatial Analysis Coordinate Systems 3D (Volumetric) 2D (Surface-based) Inter-subject registration Volume-based Surface-based

More information

Ischemic Stroke Lesion Segmentation Proceedings 5th October 2015 Munich, Germany

Ischemic Stroke Lesion Segmentation   Proceedings 5th October 2015 Munich, Germany 0111010001110001101000100101010111100111011100100011011101110101101012 Ischemic Stroke Lesion Segmentation www.isles-challenge.org Proceedings 5th October 2015 Munich, Germany Preface Stroke is the second

More information

A Binary Entropy Measure to Assess Nonrigid Registration Algorithms

A Binary Entropy Measure to Assess Nonrigid Registration Algorithms A Binary Entropy Measure to Assess Nonrigid Registration Algorithms Simon K. Warfield 1, Jan Rexilius 1, Petra S. Huppi 2, Terrie E. Inder 3, Erik G. Miller 1, William M. Wells III 1, Gary P. Zientara

More information

Application of fuzzy set theory in image analysis. Nataša Sladoje Centre for Image Analysis

Application of fuzzy set theory in image analysis. Nataša Sladoje Centre for Image Analysis Application of fuzzy set theory in image analysis Nataša Sladoje Centre for Image Analysis Our topics for today Crisp vs fuzzy Fuzzy sets and fuzzy membership functions Fuzzy set operators Approximate

More information

STIC AmSud Project. Graph cut based segmentation of cardiac ventricles in MRI: a shape-prior based approach

STIC AmSud Project. Graph cut based segmentation of cardiac ventricles in MRI: a shape-prior based approach STIC AmSud Project Graph cut based segmentation of cardiac ventricles in MRI: a shape-prior based approach Caroline Petitjean A joint work with Damien Grosgeorge, Pr Su Ruan, Pr JN Dacher, MD October 22,

More information

Simon K. Warfield, Michael Kaus, Ferenc A. Jolesz and Ron Kikinis

Simon K. Warfield, Michael Kaus, Ferenc A. Jolesz and Ron Kikinis Medical Image Analysis (2000) volume 4, number 1, pp 43 55 c Elsevier Science BV Adaptive, Template Moderated, Spatially Varying Statistical Classification Simon K. Warfield, Michael Kaus, Ferenc A. Jolesz

More information

LOCUS: LOcal Cooperative Unified Segmentation of MRI Brain Scans

LOCUS: LOcal Cooperative Unified Segmentation of MRI Brain Scans Author manuscript, published in "MICCAI07-10th International Conference on Medical Image Computing and Computer Assisted Intervention, Brisbane : Australia (2007)" DOI : 10.1007/978-3-540-75757-3_27 LOCUS:

More information

Whole Body MRI Intensity Standardization

Whole Body MRI Intensity Standardization Whole Body MRI Intensity Standardization Florian Jäger 1, László Nyúl 1, Bernd Frericks 2, Frank Wacker 2 and Joachim Hornegger 1 1 Institute of Pattern Recognition, University of Erlangen, {jaeger,nyul,hornegger}@informatik.uni-erlangen.de

More information

Multi-Sectional Views Textural Based SVM for MS Lesion Segmentation in Multi-Channels MRIs

Multi-Sectional Views Textural Based SVM for MS Lesion Segmentation in Multi-Channels MRIs 56 The Open Biomedical Engineering Journal, 2012, 6, 56-72 Open Access Multi-Sectional Views Textural Based SVM for MS Lesion Segmentation in Multi-Channels MRIs Bassem A. Abdullah*, Akmal A. Younis and

More information

Discriminative, Semantic Segmentation of Brain Tissue in MR Images

Discriminative, Semantic Segmentation of Brain Tissue in MR Images Discriminative, Semantic Segmentation of Brain Tissue in MR Images Zhao Yi 1, Antonio Criminisi 2, Jamie Shotton 2, and Andrew Blake 2 1 University of California, Los Angeles, USA. zyi@ucla.edu. 2 Microsoft

More information

Measuring longitudinal brain changes in humans and small animal models. Christos Davatzikos

Measuring longitudinal brain changes in humans and small animal models. Christos Davatzikos Measuring longitudinal brain changes in humans and small animal models Christos Davatzikos Section of Biomedical Image Analysis University of Pennsylvania (Radiology) http://www.rad.upenn.edu/sbia Computational

More information

Segmentation of Cerebral MRI Scans Using a Partial Volume Model, Shading Correction, and an Anatomical Prior

Segmentation of Cerebral MRI Scans Using a Partial Volume Model, Shading Correction, and an Anatomical Prior Segmentation of Cerebral MRI Scans Using a Partial Volume Model, Shading Correction, and an Anatomical Prior Aljaž Noe a, Stanislav Kovačič a, James C. Gee b a Faculty of Electrical Engineering, University

More information

Neuroimage Processing

Neuroimage Processing Neuroimage Processing Instructor: Moo K. Chung mkchung@wisc.edu Lecture 2-3. General Linear Models (GLM) Voxel-based Morphometry (VBM) September 11, 2009 What is GLM The general linear model (GLM) is a

More information

Image Segmentation Based on Watershed and Edge Detection Techniques

Image Segmentation Based on Watershed and Edge Detection Techniques 0 The International Arab Journal of Information Technology, Vol., No., April 00 Image Segmentation Based on Watershed and Edge Detection Techniques Nassir Salman Computer Science Department, Zarqa Private

More information

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem SPRING 07 MEDICAL IMAGE COMPUTING (CAP 97) LECTURE 0: Medical Image Segmentation as an Energy Minimization Problem Dr. Ulas Bagci HEC, Center for Research in Computer Vision (CRCV), University of Central

More information

A Design Toolbox to Generate Complex Phantoms for the Evaluation of Medical Image Processing Algorithms

A Design Toolbox to Generate Complex Phantoms for the Evaluation of Medical Image Processing Algorithms A Design Toolbox to Generate Complex Phantoms for the Evaluation of Medical Image Processing Algorithms Omar Hamo, Georg Nelles, Gudrun Wagenknecht Central Institute for Electronics, Research Center Juelich,

More information

Automatic White Matter Lesion Segmentation in FLAIR MRI

Automatic White Matter Lesion Segmentation in FLAIR MRI Automatic White Matter Lesion Segmentation in FLAIR MRI Sai Prasanth.S, Shinu Rajan Mathew 2, Mobin Varghese Mathew 3 & Rahul.B.R 4 UG Student, Department of ECE, Rajas Engineering College, Vadakkangulam,

More information

NIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2014 October 07.

NIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2014 October 07. NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2014 March 21; 9034: 903442. doi:10.1117/12.2042915. MRI Brain Tumor Segmentation and Necrosis Detection

More information

Atlas Based Segmentation of the prostate in MR images

Atlas Based Segmentation of the prostate in MR images Atlas Based Segmentation of the prostate in MR images Albert Gubern-Merida and Robert Marti Universitat de Girona, Computer Vision and Robotics Group, Girona, Spain {agubern,marly}@eia.udg.edu Abstract.

More information

Segmenting Lesions in Multiple Sclerosis Patients James Chen, Jason Su

Segmenting Lesions in Multiple Sclerosis Patients James Chen, Jason Su Segmenting Lesions in Multiple Sclerosis Patients James Chen, Jason Su Radiologists and researchers spend countless hours tediously segmenting white matter lesions to diagnose and study brain diseases.

More information

STATISTICAL ATLAS-BASED SUB-VOXEL SEGMENTATION OF 3D BRAIN MRI

STATISTICAL ATLAS-BASED SUB-VOXEL SEGMENTATION OF 3D BRAIN MRI STATISTICA ATAS-BASED SUB-VOXE SEGMENTATION OF 3D BRAIN MRI Marcel Bosc 1,2, Fabrice Heitz 1, Jean-Paul Armspach 2 (1) SIIT UMR-7005 CNRS / Strasbourg I University, 67400 Illkirch, France (2) IPB UMR-7004

More information

BAYESIAN SEGMENTATION OF THREE DIMENSIONAL IMAGES USING THE EM/MPM ALGORITHM. A Thesis. Submitted to the Faculty.

BAYESIAN SEGMENTATION OF THREE DIMENSIONAL IMAGES USING THE EM/MPM ALGORITHM. A Thesis. Submitted to the Faculty. BAYESIAN SEGMENTATION OF THREE DIMENSIONAL IMAGES USING THE EM/MPM ALGORITHM A Thesis Submitted to the Faculty of Purdue University by Lauren Christopher In Partial Fulfillment of the Requirements for

More information

Image Analysis Image Segmentation (Basic Methods)

Image Analysis Image Segmentation (Basic Methods) Image Analysis Image Segmentation (Basic Methods) Christophoros Nikou cnikou@cs.uoi.gr Images taken from: R. Gonzalez and R. Woods. Digital Image Processing, Prentice Hall, 2008. Computer Vision course

More information

Fmri Spatial Processing

Fmri Spatial Processing Educational Course: Fmri Spatial Processing Ray Razlighi Jun. 8, 2014 Spatial Processing Spatial Re-alignment Geometric distortion correction Spatial Normalization Smoothing Why, When, How, Which Why is

More information

Medical Image Registration by Maximization of Mutual Information

Medical Image Registration by Maximization of Mutual Information Medical Image Registration by Maximization of Mutual Information EE 591 Introduction to Information Theory Instructor Dr. Donald Adjeroh Submitted by Senthil.P.Ramamurthy Damodaraswamy, Umamaheswari Introduction

More information

Modified Expectation Maximization Method for Automatic Segmentation of MR Brain Images

Modified Expectation Maximization Method for Automatic Segmentation of MR Brain Images Modified Expectation Maximization Method for Automatic Segmentation of MR Brain Images R.Meena Prakash, R.Shantha Selva Kumari 1 P.S.R.Engineering College, Sivakasi, Tamil Nadu, India 2 Mepco Schlenk Engineering

More information

Graphical Models, Bayesian Method, Sampling, and Variational Inference

Graphical Models, Bayesian Method, Sampling, and Variational Inference Graphical Models, Bayesian Method, Sampling, and Variational Inference With Application in Function MRI Analysis and Other Imaging Problems Wei Liu Scientific Computing and Imaging Institute University

More information

Structural Segmentation

Structural Segmentation Structural Segmentation FAST tissue-type segmentation FIRST sub-cortical structure segmentation FSL-VBM voxelwise grey-matter density analysis SIENA atrophy analysis FAST FMRIB s Automated Segmentation

More information

Integrated Approaches to Non-Rigid Registration in Medical Images

Integrated Approaches to Non-Rigid Registration in Medical Images Work. on Appl. of Comp. Vision, pg 102-108. 1 Integrated Approaches to Non-Rigid Registration in Medical Images Yongmei Wang and Lawrence H. Staib + Departments of Electrical Engineering and Diagnostic

More information

Mine Boundary Detection Using Markov Random Field Models

Mine Boundary Detection Using Markov Random Field Models Mine Boundary Detection Using Markov Random Field Models Xia Hua and Jennifer Davidson Department of Electrical & Computer Engineering Iowa State University Noel Cressie Department of Statistics Iowa State

More information

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem

MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem SPRING 06 MEDICAL IMAGE COMPUTING (CAP 97) LECTURE 0: Medical Image Segmentation as an Energy Minimization Problem Dr. Ulas Bagci HEC, Center for Research in Computer Vision (CRCV), University of Central

More information

Level Set Evolution with Region Competition: Automatic 3-D Segmentation of Brain Tumors

Level Set Evolution with Region Competition: Automatic 3-D Segmentation of Brain Tumors 1 Level Set Evolution with Region Competition: Automatic 3-D Segmentation of Brain Tumors 1 Sean Ho, 2 Elizabeth Bullitt, and 1,3 Guido Gerig 1 Department of Computer Science, 2 Department of Surgery,

More information

Automated Brain-Tissue Segmentation by Multi-Feature SVM Classification

Automated Brain-Tissue Segmentation by Multi-Feature SVM Classification Automated Brain-Tissue Segmentation by Multi-Feature SVM Classification Annegreet van Opbroek 1, Fedde van der Lijn 1 and Marleen de Bruijne 1,2 1 Biomedical Imaging Group Rotterdam, Departments of Medical

More information

Markov Networks in Computer Vision. Sargur Srihari

Markov Networks in Computer Vision. Sargur Srihari Markov Networks in Computer Vision Sargur srihari@cedar.buffalo.edu 1 Markov Networks for Computer Vision Important application area for MNs 1. Image segmentation 2. Removal of blur/noise 3. Stereo reconstruction

More information

Dept of CSE, CIT Gubbi, Tumkur, Mysore, India

Dept of CSE, CIT Gubbi, Tumkur, Mysore, India Volume 5, Issue 6, June 2015 ISSN: 2277 128X International Journal of Advanced Research in Computer Science and Software Engineering Research Paper Available online at: www.ijarcsse.com MRI Tissue Segmentation

More information

Multi-Atlas Segmentation of the Cardiac MR Right Ventricle

Multi-Atlas Segmentation of the Cardiac MR Right Ventricle Multi-Atlas Segmentation of the Cardiac MR Right Ventricle Yangming Ou, Jimit Doshi, Guray Erus, and Christos Davatzikos Section of Biomedical Image Analysis (SBIA) Department of Radiology, University

More information

Short Survey on Static Hand Gesture Recognition

Short Survey on Static Hand Gesture Recognition Short Survey on Static Hand Gesture Recognition Huu-Hung Huynh University of Science and Technology The University of Danang, Vietnam Duc-Hoang Vo University of Science and Technology The University of

More information

Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging

Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging 1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant

More information

Automatic Subthalamic Nucleus Targeting for Deep Brain Stimulation. A Validation Study

Automatic Subthalamic Nucleus Targeting for Deep Brain Stimulation. A Validation Study Automatic Subthalamic Nucleus Targeting for Deep Brain Stimulation. A Validation Study F. Javier Sánchez Castro a, Claudio Pollo a,b, Jean-Guy Villemure b, Jean-Philippe Thiran a a École Polytechnique

More information

A Review on Label Image Constrained Multiatlas Selection

A Review on Label Image Constrained Multiatlas Selection A Review on Label Image Constrained Multiatlas Selection Ms. VAIBHAVI NANDKUMAR JAGTAP 1, Mr. SANTOSH D. KALE 2 1PG Scholar, Department of Electronics and Telecommunication, SVPM College of Engineering,

More information

A Novel Nonrigid Registration Algorithm and Applications

A Novel Nonrigid Registration Algorithm and Applications A Novel Nonrigid Registration Algorithm and Applications J. Rexilius 1, S.K. Warfield 1, C.R.G. Guttmann 1, X. Wei 1, R. Benson 2, L. Wolfson 2, M. Shenton 1, H. Handels 3, and R. Kikinis 1 1 Surgical

More information

IMAGE SEGMENTATION BY FUZZY C-MEANS CLUSTERING ALGORITHM WITH A NOVEL PENALTY TERM

IMAGE SEGMENTATION BY FUZZY C-MEANS CLUSTERING ALGORITHM WITH A NOVEL PENALTY TERM Computing and Informatics, Vol. 26, 2007, 17 31 IMAGE SEGMENTATION BY FUZZY C-MEANS CLUSTERING ALGORITHM WITH A NOVEL PENALTY TERM Yong Yang School of Information Management Jiangxi University of Finance

More information

Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques

Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October

More information

Segmentation of 3-D medical image data sets with a combination of region based initial segmentation and active surfaces

Segmentation of 3-D medical image data sets with a combination of region based initial segmentation and active surfaces Header for SPIE use Segmentation of 3-D medical image data sets with a combination of region based initial segmentation and active surfaces Regina Pohle, Thomas Behlau, Klaus D. Toennies Otto-von-Guericke

More information

Non-rigid Image Registration

Non-rigid Image Registration Overview Non-rigid Image Registration Introduction to image registration - he goal of image registration - Motivation for medical image registration - Classification of image registration - Nonrigid registration

More information

IMAGE SEGMENTATION. Václav Hlaváč

IMAGE SEGMENTATION. Václav Hlaváč IMAGE SEGMENTATION Václav Hlaváč Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/ hlavac, hlavac@fel.cvut.cz

More information

Efficient population registration of 3D data

Efficient population registration of 3D data Efficient population registration of 3D data Lilla Zöllei 1, Erik Learned-Miller 2, Eric Grimson 1, William Wells 1,3 1 Computer Science and Artificial Intelligence Lab, MIT; 2 Dept. of Computer Science,

More information

Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit

Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit Knowledge-Based Segmentation of Brain MRI Scans Using the Insight Toolkit John Melonakos 1, Ramsey Al-Hakim 1, James Fallon 2 and Allen Tannenbaum 1 1 Georgia Institute of Technology, Atlanta GA 30332,

More information

MEDICAL IMAGE NOISE REDUCTION AND REGION CONTRAST ENHANCEMENT USING PARTIAL DIFFERENTIAL EQUATIONS

MEDICAL IMAGE NOISE REDUCTION AND REGION CONTRAST ENHANCEMENT USING PARTIAL DIFFERENTIAL EQUATIONS MEDICAL IMAGE NOISE REDUCTION AND REGION CONTRAST ENHANCEMENT USING PARTIAL DIFFERENTIAL EQUATIONS Miguel Alemán-Flores, Luis Álvarez-León Departamento de Informática y Sistemas, Universidad de Las Palmas

More information

Texture Modeling using MRF and Parameters Estimation

Texture Modeling using MRF and Parameters Estimation Texture Modeling using MRF and Parameters Estimation Ms. H. P. Lone 1, Prof. G. R. Gidveer 2 1 Postgraduate Student E & TC Department MGM J.N.E.C,Aurangabad 2 Professor E & TC Department MGM J.N.E.C,Aurangabad

More information