Non-rigid Image Registration using Electric Current Flow
|
|
- Eleanore Jewel Boone
- 5 years ago
- Views:
Transcription
1 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, The Hong Kong University of Science and Technology, Hong Kong Abstract In this paper, a novel electric current flow (ECF) based model is proposed to perform feature based non-rigid brain image registration. The ECF features simultaneously capture both voxel intensity and inter-voxel distance information. In the proposed ECF framework, each voxel is regarded as exhibiting electric potential proportional to voxel intensity. Voxels are connected by conductive wires in a pairwise manner. Each conductive wire has resistance, in which the resistance value is proportional to the length of the wire. The electric potential difference among connected pixels induces electric current passing through their connected wire. The amount of the electric current is the ratio between the voxel potential difference and the wire resistance. The potential difference and resistance are respectively proportional to the voxel intensity difference and the inter-voxel distance. By analyzing the electric current induced by the connection between a reference voxel and its counterparts in a given range, the ECF algorithm searches for the most salient connection to construct the ECF features. The ECF features are incorporated in the Markov random field labeling framework for non-rigid image registration. The registration quality of the proposed method has been evaluated intensively on both BrainWeb and IBSR databases. It is compared with four related approaches. Experimental results illustrate that the proposed method consistently achieves the highest registration accuracy among all the compared methods on both databases. 1. Introduction A variety of non-rigid image registration approaches have been developed for medical image analysis. Intensity based approaches have widely been used for registering template and subject images based on inter-image intensity similarity estimated from image joint intensity distributions. However, intensity similarity does not necessarily equivalent to anatomical similarity. As pointed out in [9], the mere use of intensity information is possibly ineffective to distinguish between different brain anatomical structures. On the other hand, feature based registration methods use feature vector as anatomical signature for each voxel, where the registration problem is formulated as a feature matching problem. Observing intensity or gradient similarity between subject and template images are common criteria for finding correspondences [4]. Apart from voxel intensity and local image gradient, as will be demonstrated in this paper, it is beneficial to consider inter-voxel relation between distant voxels rather than relation between adjacent voxels. It is because the distant voxel relation offers additional intensity and geometric information within a predefined region of a voxel. A new feature, namely Electric Current Flow (ECF), is proposed for representing both the geometric information (inter-voxel distance) and intensity information (voxel intensity difference) at various scales. The ECF feature is integrated with the Markov random field deformation model to perform non-rigid image registration. It is evaluated on two publicly available databases - BrainWeb and IBSR. The proposed method is compared with four state-of-the-art registration approaches. 2. The electric current flow (ECF) feature Each pair of voxels, denoted as v 1 and v 2, is assumed to be connected by a conductive wire. At each voxel position, it is assumed that its electric potential is represented by the its intensity, denoted as G( v 1 ) and G( v 2 ), respectively. The relationship between v 1 and v 2 is established based on the Ohm s law. The voltage between v 1 and v 2 is therefore the absolute value of the potential difference between them, which is expressed as V ( v 1, v 2 ) = G( v 1 ) G( v 2 ). The resistance between v 1 and v 2 is modeled by their voxel distance, which is given by R( v 1, v 2 ) = v 1 v 2. It is therefore assumed that the resistance increases proportional to the distance between two voxels. Such assumption reflects the physical phenomenon that the longer the wire, the larger the resistance. Then, the amount of current passing through v 1 and v 2 is obtained by dividing the voltage value 169
2 by the resistance according to the Ohm s law which is given as I( v 1, v 2 ) = V ( v 1, v 2 ) R( v 1, v 2 ) = G( v 1) G( v 2 ). (1) v 1 v 2 In this paper, it is assumed that the larger the amount of current passing through v 1 and v 2, the more salient the anatomical information encoded in the inter-voxel relationship between v 1 and v 2. For a reference voxel v, voxel v is called the ECF source point of v if the amount of current calculated by Equation 1 passing through v and v is maximized with respect to v. The ECF source point v of v can be expressed by Equation 2, v = arg max I( v, v), (2) v Ω where Ω represents the image domain. The amount of current passing through v and v is used as the voxel signature for v, which is defined by Equation 3 ECF ( v) = I( v, v), (3) where v is the ECF source point of v calculated by Equation 2. As pointed out in [10], feature saliency is closely related to the scale from which it is extracted. Extracting features from only a single scale of interest generally is insufficient to capture all the anatomical information around the reference voxel. To extract the ECF features from different scales, the search domains for finding the ECF source points in Equation 2 are confined in spheres centered at v with different values of radius r instead of the whole image domain, denoted as v r, where r reflects the scale of interest. Therefore, now the ECF signature of v with scale of interest r is defined as, ECF ( v, r) = I( ˆv r, v) = V ( ˆv r, v) R( ˆv r, v), (4) where ˆv r is the ECF source point of v with scale of interest r, it is defined as: ˆvr = arg max I( v, v). (5) v v r The ECF feature signature of each reference voxel is obtained by concatenating the ECF features calculated by Equation 4 from different scales r. Algorithm 1 summarizes the procedure of calculating the ECF feature signatured image F with k scales of interest. The ECF feature has several distinctive properties. Property 1: The ECF feature is 3D rotation invariant. It is given that the ECF feature is the maximum amount of current passing through the reference voxel and Algorithm 1 Computation of the ECF feature signature Input: Input image G, k scales of interest: r 1,...,r k. Output: The ECF feature signatured image F of G. 1. FOR each voxel v G 2. Initialize an empty feature vector Y. 3. FOR i (1...k) 4. Calculate the ECF feature of voxel v at scale r i by Equation 4, and the calculated ECF feature is denoted as f i. 5. Y = [ Y ; fi]. 6. END FOR 7. F ( v) = Y. 8. END FOR 9. Return F. its corresponding ECF source point within a local spherical search domain, as stated in Eqns 4 and 5. Although the image is rotated, the ECF source point remains the same within such spherical search domain because the maximum current is still the same between the reference voxel and its original ECF source point. Thus the ECF feature remains the same, despite of image rotation. Property 2: The conventional gradient magnitude feature is a special case of the ECF feature. It should be noted that the ECF source point is calculated based on Equation 5, where the ECF source point is found within a spherical search domain centered at the reference voxel with a specified scale of interest r (i.e., the radius of the spherical search domain). When r equals one, it means that the ECF source point is now only found within the immediate neighboring voxels with respect to the reference voxel. In this case the denominator of Equation 4 is equal to one, and the numerator of Equation 4 now specifies the largest intensity changes between the reference voxel and its immediate neighboring voxels. This is exactly the physical meaning of the conventional gradient magnitude feature. Property 3: The maximum distance gradient magnitude (MDGM) feature [2] is a special case of the ECF feature. The MDGM feature is defined based on the distance gradient operator d between two arbitrary different voxels v 1 and v 2, which is defined as [2]: v d G( v 1, v 2 ) = (G( v 1 ) G( v 2 )) 1 v 2 v 1 v 2, where G 2 denotes the input image, G( v 1 ) and G( v 2 ) denote the input image intensity values at positions v 1 and v 2, respectively. Then the MDGM feature for voxel v is defined as, MDGM( v) = d G( v, v), where v is the MDGM source point of v defined as: v = arg max v Ω dg( v, v), where Ω is the image domain of G. Therefore, the ECF feature becomes exactly the MDGM feature when the spherical search domain in Equation 5 to find the ECF source point is large enough to cover the whole image domain. 170
3 (a) (b) (c) (d) (e) (f) (g) Figure 1. (a) An image obtained from BrainWeb, the features of the voxel highlighted with the green cross is compared with the features of other voxels; (b) The difference map of image shown in (a) using ECF, blue indicates high similarity; (c) The difference map obtained by comparing the ECF feature of the highlighted voxel in (a) with all the voxels in the image shown in (d); (d) An image obtained from BrainWeb with a different patient from (a). The green cross in (d) denotes the voxel in the image (d) with the most similar ECF feature compared to the highlighted voxel in (a); (e) and (f) are the intra-patient and inter-patient difference maps using MDGM [2] respectively. The red cross in (g) denotes the voxel in the image (d) with the most similar conventional MDGM feature compared to the highlighted voxel in (a). Please refer to the electronic version of this paper for better visualization. We now demonstrate the superior discriminant power of the ECF feature signature over the MDGM feature. Figure 1a is an MR image obtained from BrainWeb, and the green cross in Figure 1a denotes the referencing voxel v ref. The ECF feature signature of this referencing voxel is compared with the ECF feature signatures of other voxels. Figure 1b is the color-coded difference map of the image shown in Figure 1a with the ECF feature signatures calculated from three different scales r 1 = 3, r 2 = 5, r 3 = 8 for each voxel. Figure 1b shows that (from the blue voxels) the ECF features of v ref in Fig. 1a are only similar to a small number of voxels with similar anatomical properties with v ref (i.e., near the ventricle corner). It reflects the discriminant power of ECF feature signature to identify different anatomical structures. Moreover, voxels of similar anatomical positions across different patients should have similar feature signatures for establishing reliable correspondences. The difference map between the ECF feature signatures of v ref in Figure 1a and all the voxels in the image shown in Figure 1d is also calculated and shown in Figure 1c. The MR image shown in Figure 1d is also obtained from BrainWeb but with a different patient from the one shown in Figure 1a. It is observed that the ECF feature signatures of v ref in Figure 1a are only similar to the voxels in Figure 1c which are also located near the ventricle corner. For the image shown in Figure 1d, the voxel with the most similar ECF feature signature compared to v ref in Figure 1a is also highlighted with a green cross. It is also located near the ventricle corner. Thus, the ECF feature signatures can also provide reliable correspondences between voxels with similar anatomical properties across two patients. The difference map generated by the MDGM feature [2] is also included for comparison. Figures 1e and f are the intra-patient and inter-patient difference maps obtained by the MDGM feature with the same referencing pixel in Figure 1a respectively, and lots of false similarity matches are identified. The voxel with the most similar MDGM feature of the image shown in Figure 1d compared to v ref in Figure 1a is highlighted with a red cross shown in Figure 1g, which is obviously located at a different anatomical position. Therefore, the higher discriminant power of ECF over MDGM is implied in this example. 171
4 3. Deformation Model for ECF In this paper, the ECF signatured image is integrated with the Markov random field (MRF) labeling framework [3, 7] to conduct registration. The general form of the MRF energy function is defined as, E f = E data + E smoothness = D p (l p ) + V p,q (l p, l q ), (6) p Ω (p,q) N where Ω is the image domain, and N is the neighborhood system. The 6-connected neighborhood system is used in this paper. D p (l p ) denotes the cost function of the data term. It represents the cost of assigning label l p to voxel p. V p,q (l p, l q ) is the cost function of the smoothness term, and it penalizes the cost of label discrepancy between two neighboring voxels. The deformation space is quantized into a discrete set of labels L {l 1, l 2,..., l n }. Each label l i (1 i n) corresponds to a displacement vector d i. Labeling the voxel p with label l p denotes moving p to a new position by the displacement vector dlp. The quantization step described in [7] is adopted in this paper, where each voxel can move from the original position to a new position bounded by a discretized window Ψ = {0, ±s, ±2s,..., ±ws} d with dimension d. In this paper, we set s = 1. The cost function D p (l p ) of the data term is defined based on the ECF features: D p (l p ) = D p (G template ( p), G subject ( p + d lp ) = D p ( F template ( p), F subject ( p + d lp )) = F template ( p) F subject ( p + d lp ) L2,(7) where G template and G subject denote the template and subject images, F template and F subject are the ECF feature signature images of the template and subject images respectively calculated from different scales of interest by Algorithm 1 at the current iteration. L2 denotes the L 2 distance. The piece-wise truncated absolute distance is used as the smoothness cost function given by V p,q (l p, l q ) = min(λ, d lp d lq ), where λ is a constant representing the maximum penalty, and it was set to 20 in this paper. The α-expansion algorithm [11] is used to minimize the energy function defined in Equation 6 in this paper and it guarantees that the solution is within a known constant range of the global minimum. 4. Experimental Results We have evaluated the proposed method on both the simulated and real 3D databases obtained from BrainWeb where A and B denote the regions of a specific tissue in two images. The mean values and SDs of P for WM, GM and CSF are listed in Table 1 for various methods. It is observed that the proposed ECF feature achieves the highest value of P among all the compared methods. The average computation times for FFD, Demons, MDGM, HAMMER and ECF are 4.8hr, 2.6hr, 18.5hr, 5.4hr and 5.2hr, respectively. Then, we evaluated the proposed method on the IBSR dataset with all twelve brain MR images. Each image has size of voxels. The expert segmentation results of WM, GM and CSF are provided by IBSR. The ex- and IBSR 2, respectively. The proposed method was compared with three widely used methods: FFD [5], Demons [8] and HAMMER [6]. For comparison purpose, the results obtained by using the MDGM feature [2] with the same MRF based deformation model are included. In all the experiments, five different scales were used for the ECF features: r 1 = 3, r 2 = 5, r 3 = 8, r 4 = 10, r 5 = 12. The 3D displacement window used in this paper was Ψ = {0, ±1, ±2,..., ±12} 3. The control point spacing for FFD [5] was set to 2.5mm as suggested in [1]. Each image was preprocessed by a skull stripping process using the Brainsuite software 3 as it is a required step for the HAMMER [6] method compared in this paper. Affine registrations were first applied to register the subject images to the template image before conducting the non-rigid image registration. All the experiments were run on a 3GHz duo core CPU, 2GB RAM computer. The FFD and Demons algorithms were implemented based on ITK 4. For HAMMER, its software provided by the authors of HAMMER was adopted 5. The proposed method was also implemented based on ITK. For the BrainWeb dataset, all twenty brain MR images were obtained. Each image has volume size of voxels and the segmentation results of each image, including white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF), are provided by BrainWeb. One of the images was selected as the template image, and others were used as subject images. To visualize the registration accuracy, the same cross section of the average brain images reconstructed from the warped subject images obtained by using different approaches are shown in Figure 2. The same cross section of the template image is also included in Figure 2a for reference. It is observed from Figure 2 that the average brain image obtained by using the proposed method is most similar to the template image than those obtained via other methods. The tissue overlap measure proposed in [1] was also adopted to quantitatively evaluate the registration accuracy of different methods. It is defined as P = N(A B) N(A B)
5 (a) Template (b) Our Method (c) MDGM (d) HAMMER (e) FFD (f) Demons Figure 2. (a) The template image; From (b) to (f): The average brain images obtained by using the proposed method, MDGM [2], HAMMER [6], FFD [5] and Demons [8] respectively for the BrainWeb database. Regions with significant differences are highlighted with green circles. Note that the skulls of input images had already been removed by using the Brainsuite software. Tissue BR FFD Demons MDGM HAMMER ECF Gray 0.402± ± ± ± ± ±0.04 White 0.462± ± ± ± ± ±0.06 CSF 0.393± ± ± ± ± ±0.01 Table 1. The mean values of P and SDs of GM, WM and CSF with different methods on the simulated 3D database. BR denotes before registration. periment setting is similar to the setting in the experiments on BrainWeb datasets. The means and SDs of P are listed in Table 2 for various methods. It is observed that the proposed method consistently achieves the highest value of P among all the compared methods, and the robustness of the proposed method is strongly implied. The average computation times for FFD, Demons, MDGM, HAMMER and ECF are 3.1hr, 1.5hr, 13.7hr, 4.6hr and 3.9hr, respectively. It should be noted that the MDGM feature [2] does not give very good registration accuracies on both the BrainWeb and IBSR datasets, which reflects the importance of extracting anatomical features from input images with multiple scales. 5. Conclusion A feature, namely Electric Current Flow (ECF), is proposed for feature based non-rigid image registration. The ECF feature aims at modeling the inter-voxel relationship as the electric current flow. It is demonstrated that both the conventional gradient magnitude and the recently proposed maximum distance gradient magnitude feature are the spe- 173
6 Tissue BR FFD Demons MDGM HAMMER ECF Gray 0.427± ± ± ± ± ±0.06 White 0.464± ± ± ± ± ±0.05 CSF 0.403± ± ± ± ± ±0.06 Table 2. The mean values of P and SDs of GM, WM and CSF with different methods on the real 3D database. BR denotes before registration. cial cases of the proposed ECF feature. The ECF features are integrated with the Markov random field based deformation model to drive the non-rigid registration process. The proposed method has been compared with four related registration approaches (FFD, Demons, MDGM and HAM- MER) on both the simulated and real 3D image databases. It is shown that the proposed method consistently achieves the highest registration accuracies among all the compared methods and its computation efficiency is comparable. [11] B. Yuri, V. Olga, and Z. Ramin. Fast approximate energy minimization via graph cuts. PAMI, 23: , Acknowledgment The authors would like to acknowledge the generous support from the K. S. Lo Foundation and Hong Kong GRF grant References [1] W. Crum, D. Rueckert, M. Jenkinson, D. Kennedy, and M. Smith. A framework for detailed objective comparison of non-rigid registration algorithms in neuroimaging. In MIC- CAI, pages , [2] R. Gan, A. Chung, and S. Liao. Maximum distance-gradient for robust image registration. MedIA, 12: , [3] B. Glocker, N. Komodakis, N. Paragios, G. Tziritas, and N. Navab. Inter and intra-modal deformable registration: Continuous deformations meet efficient optimal linear programming. In IPMI, pages , [4] J. Pluim, J. Maintz, and M. Viergever. Image registration by maximization of combined mutual information and gradient information. TMI, 19: , [5] D. Rueckert, L. Sonoda, C. Hayes, D. Hill, M. Leach, and D. Hawkes. Nonrigid registration using free-form deformations: application to breast mr images. TMI, 18: , [6] D. Shen and C. Davatzikos. Hammer: hierarchical attribute matching mechanism for elastic registration. TMI, 21: , [7] W. Tang and A. Chung. Non-rigid image registration using graph-cuts. In MICCAI, pages , [8] J. Thirion. Image matching as a diffusion process: an analogy with maxwell s demons. Med. Image Anal., 3: , [9] Z. Tu, K. Narr, P. Dollar, I. Dinov, P. Thompson, and A. Toga. Brain anatomical structure segmentation by hybrid discriminative/generative models. TMI, 27: , [10] G. Wu, F. Qi, and D. Shen. Learning-based deformable registration of mr brain images. TMI, 25: ,
Attribute Similarity and Mutual-Saliency Weighting for Registration and Label Fusion
Attribute Similarity and Mutual-Saliency Weighting for Registration and Label Fusion Yangming Ou, Jimit Doshi, Guray Erus, and Christos Davatzikos Section of Biomedical Image Analysis (SBIA) Department
More informationNonrigid 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 informationDRAMMS: Deformable Registration via Attribute Matching and Mutual-Saliency weighting
DRAMMS: Deformable Registration via Attribute Matching and Mutual-Saliency weighting Yangming Ou, Christos Davatzikos Section of Biomedical Image Analysis (SBIA) University of Pennsylvania Outline 1. Background
More informationNon-rigid Registration using Discrete MRFs: Application to Thoracic CT Images
Non-rigid Registration using Discrete MRFs: Application to Thoracic CT Images Ben Glocker 1, Nikos Komodakis 2, Nikos Paragios 3,4, and Nassir Navab 1 1 Computer Aided Medical Procedures (CAMP), TU Mï
More informationMulti-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 informationCorrespondence 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 informationADAPTIVE 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 informationRegional Manifold Learning for Deformable Registration of Brain MR Images
Regional Manifold Learning for Deformable Registration of Brain MR Images Dong Hye Ye, Jihun Hamm, Dongjin Kwon, Christos Davatzikos, and Kilian M. Pohl Department of Radiology, University of Pennsylvania,
More informationEnsemble registration: Combining groupwise registration and segmentation
PURWANI, COOTES, TWINING: ENSEMBLE REGISTRATION 1 Ensemble registration: Combining groupwise registration and segmentation Sri Purwani 1,2 sri.purwani@postgrad.manchester.ac.uk Tim Cootes 1 t.cootes@manchester.ac.uk
More informationAutomatic 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 informationGLIRT: Groupwise and Longitudinal Image Registration Toolbox
Software Release (1.0.1) Last updated: March. 30, 2011. GLIRT: Groupwise and Longitudinal Image Registration Toolbox Guorong Wu 1, Qian Wang 1,2, Hongjun Jia 1, and Dinggang Shen 1 1 Image Display, Enhancement,
More informationMath in image processing
Math in image processing Math in image processing Nyquist theorem Math in image processing Discrete Fourier Transformation Math in image processing Image enhancement: scaling Math in image processing Image
More informationNon-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation
M. HEINRICH et al.: MULTIMODAL REGISTRATION USING GRADIENT ORIENTATION 1 Non-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation Mattias P. Heinrich 1 mattias.heinrich@eng.ox.ac.uk
More informationHybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information
Hybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information Andreas Biesdorf 1, Stefan Wörz 1, Hans-Jürgen Kaiser 2, Karl Rohr 1 1 University of Heidelberg, BIOQUANT, IPMB,
More informationComparison 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 informationMARS: Multiple Atlases Robust Segmentation
Software Release (1.0.1) Last updated: April 30, 2014. MARS: Multiple Atlases Robust Segmentation Guorong Wu, Minjeong Kim, Gerard Sanroma, and Dinggang Shen {grwu, mjkim, gerard_sanroma, dgshen}@med.unc.edu
More informationInter and Intra-Modal Deformable Registration:
Inter and Intra-Modal Deformable Registration: Continuous Deformations Meet Efficient Optimal Linear Programming Ben Glocker 1,2, Nikos Komodakis 1,3, Nikos Paragios 1, Georgios Tziritas 3, Nassir Navab
More informationTransformation Model and Constraints Cause Bias in Statistics on Deformation Fields
Transformation Model and Constraints Cause Bias in Statistics on Deformation Fields Torsten Rohlfing Neuroscience Program, SRI International, Menlo Park, CA, USA torsten@synapse.sri.com http://www.stanford.edu/
More informationMulti-modal Image Registration Using the Generalized Survival Exponential Entropy
Multi-modal Image Registration Using the Generalized Survival Exponential Entropy Shu Liao and Albert C.S. Chung Lo Kwee-Seong Medical Image Analysis Laboratory, Department of Computer Science and Engineering,
More informationImage Comparison on the Base of a Combinatorial Matching Algorithm
Image Comparison on the Base of a Combinatorial Matching Algorithm Benjamin Drayer Department of Computer Science, University of Freiburg Abstract. In this paper we compare images based on the constellation
More informationTopology Correction for Brain Atlas Segmentation using a Multiscale Algorithm
Topology Correction for Brain Atlas Segmentation using a Multiscale Algorithm Lin Chen and Gudrun Wagenknecht Central Institute for Electronics, Research Center Jülich, Jülich, Germany Email: l.chen@fz-juelich.de
More informationUsing K-means Clustering and MI for Non-rigid Registration of MRI and CT
Using K-means Clustering and MI for Non-rigid Registration of MRI and CT Yixun Liu 1,2 and Nikos Chrisochoides 2 1 Department of Computer Science, College of William and Mary, enjoywm@cs.wm.edu 2 Department
More informationMEDICAL 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 informationGeometric Registration for Deformable Shapes 3.3 Advanced Global Matching
Geometric Registration for Deformable Shapes 3.3 Advanced Global Matching Correlated Correspondences [ASP*04] A Complete Registration System [HAW*08] In this session Advanced Global Matching Some practical
More informationSemantic Context Forests for Learning- Based Knee Cartilage Segmentation in 3D MR Images
Semantic Context Forests for Learning- Based Knee Cartilage Segmentation in 3D MR Images MICCAI 2013: Workshop on Medical Computer Vision Authors: Quan Wang, Dijia Wu, Le Lu, Meizhu Liu, Kim L. Boyer,
More informationGraph-based Deformable Image Registration
Graph-based Deformable Image Registration Aristeidis Sotiras, Yangming Ou, Nikos Paragios and Christos Davatzikos Abstract Deformable image registration is a field that has received considerable attention
More informationSubvoxel 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 informationMEDICAL 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 informationFree-Form B-spline Deformation Model for Groupwise Registration
Free-Form B-spline Deformation Model for Groupwise Registration The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Balci
More informationCHAPTER 6 MODIFIED FUZZY TECHNIQUES BASED IMAGE SEGMENTATION
CHAPTER 6 MODIFIED FUZZY TECHNIQUES BASED IMAGE SEGMENTATION 6.1 INTRODUCTION Fuzzy logic based computational techniques are becoming increasingly important in the medical image analysis arena. The significant
More informationMethods for data preprocessing
Methods for data preprocessing John Ashburner Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK. Overview Voxel-Based Morphometry Morphometry in general Volumetrics VBM preprocessing
More informationAtlas 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 informationAn 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 informationDRAMMS: deformable registration via attribute matching and mutual-saliency weighting
University of Pennsylvania ScholarlyCommons Departmental Papers (BE) Department of Bioengineering 7-2009 DRAMMS: deformable registration via attribute matching and mutual-saliency weighting Yangming Ou
More informationDense Image Registration through MRFs and Efficient Linear Programming
Revised Manuscript Dense Image Registration through MRFs and Efficient Linear Programming Ben Glocker a,b,, Nikos Komodakis a,c, Georgios Tziritas c, Nassir Navab b, Nikos Paragios a a GALEN Group, Laboratoire
More informationThe Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy
The Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy Sokratis K. Makrogiannis, PhD From post-doctoral research at SBIA lab, Department of Radiology,
More informationFree-Form B-spline Deformation Model for Groupwise Registration
Free-Form B-spline Deformation Model for Groupwise Registration Serdar K. Balci 1, Polina Golland 1, Martha Shenton 2, and William M. Wells 2 1 CSAIL, MIT, Cambridge, MA, USA, 2 Brigham & Women s Hospital,
More informationNon-Rigid Image Registration III
Non-Rigid Image Registration III CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (CS6240) Non-Rigid Image Registration
More informationAutomatic 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 informationBiomedical Image Analysis Using Markov Random Fields & Efficient Linear Programing
Biomedical Image Analysis Using Markov Random Fields & Efficient Linear Programing Nikos Komodakis Ahmed Besbes Ben Glocker Nikos Paragios Abstract Computer-aided diagnosis through biomedical image analysis
More informationMR 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 informationRigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach
Rigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach Julien Jomier and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab The University of North Carolina at Chapel
More informationRobust Linear Registration of CT images using Random Regression Forests
Robust Linear Registration of CT images using Random Regression Forests Ender Konukoglu a Antonio Criminisi a Sayan Pathak b Duncan Robertson a Steve White b David Haynor c Khan Siddiqui b a Microsoft
More informationIschemic 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 informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationPreprocessing 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 informationMeasuring 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 informationAutomatic Vascular Tree Formation Using the Mahalanobis Distance
Automatic Vascular Tree Formation Using the Mahalanobis Distance Julien Jomier, Vincent LeDigarcher, and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab, Department of Radiology The University
More informationReliability-Driven, Spatially-Adaptive Regularization for Deformable Registration
Reliability-Driven, Spatially-Adaptive Regularization for Deformable Registration Lisa Tang 1, Ghassan Hamarneh 1, and Rafeef Abugharbieh 2 1 Medical Image Analysis Lab., School Computing Science, Simon
More informationFast Image Registration via Joint Gradient Maximization: Application to Multi-Modal Data
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Fast Image Registration via Joint Gradient Maximization: Application to Multi-Modal Data Xue Mei, Fatih Porikli TR-19 September Abstract We
More informationPerformance 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 informationRIGID IMAGE REGISTRATION
RIGID IMAGE REGISTRATION Duygu Tosun-Turgut, Ph.D. Center for Imaging of Neurodegenerative Diseases Department of Radiology and Biomedical Imaging duygu.tosun@ucsf.edu What is registration? Image registration
More informationIntroduction to Medical Image Registration
Introduction to Medical Image Registration Sailesh Conjeti Computer Aided Medical Procedures (CAMP), Technische Universität München, Germany sailesh.conjeti@tum.de Partially adapted from slides by: 1.
More informationTranslation Symmetry Detection: A Repetitive Pattern Analysis Approach
2013 IEEE Conference on Computer Vision and Pattern Recognition Workshops Translation Symmetry Detection: A Repetitive Pattern Analysis Approach Yunliang Cai and George Baciu GAMA Lab, Department of Computing
More informationBLUT : Fast and Low Memory B-spline Image Interpolation
BLUT : Fast and Low Memory B-spline Image Interpolation David Sarrut a,b,c, Jef Vandemeulebroucke a,b,c a Université de Lyon, F-69622 Lyon, France. b Creatis, CNRS UMR 5220, F-69622, Villeurbanne, France.
More informationNonrigid Registration with Adaptive, Content-Based Filtering of the Deformation Field
Nonrigid Registration with Adaptive, Content-Based Filtering of the Deformation Field Marius Staring*, Stefan Klein and Josien P.W. Pluim Image Sciences Institute, University Medical Center Utrecht, P.O.
More informationAuxiliary Anatomical Labels for Joint Segmentation and Atlas Registration
Auxiliary Anatomical Labels for Joint Segmentation and Atlas Registration Tobias Gass, Gabor Szekely and Orcun Goksel Computer Vision Lab, ETH Zurich, Switzerland. {gasst, szekely, ogoksel}@vision.ee.ethz.ch
More informationMulti-Atlas Brain MRI Segmentation with Multiway Cut
Multi-Atlas Brain MRI Segmentation with Multiway Cut Duygu Sarikaya, Liang Zhao, and Jason J. Corso SUNY Buffalo, Computer Science and Engineering Department, 338 Davis Hall- Buffalo, New York, USA 14260-2500
More informationAutomatic 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 informationA Multiple-Layer Flexible Mesh Template Matching Method for Nonrigid Registration between a Pelvis Model and CT Images
A Multiple-Layer Flexible Mesh Template Matching Method for Nonrigid Registration between a Pelvis Model and CT Images Jianhua Yao 1, Russell Taylor 2 1. Diagnostic Radiology Department, Clinical Center,
More informationDeformable Medical Image Registration: Methods
Deformable Medical Image Registration 1 Deformable Medical Image Registration: Setting The State Of The Art With Discrete Methods Ben Glocker Microsoft Research, Cambridge, UK Aristeidis Sotiras Ecole
More informationMARS: Multiple Atlases Robust Segmentation
Software Release (1.0.1) Last updated: April 30, 2014. MARS: Multiple Atlases Robust Segmentation Guorong Wu, Minjeong Kim, Gerard Sanroma, and Dinggang Shen {grwu, mjkim, gerard_sanroma, dgshen}@med.unc.edu
More informationNeuroimaging 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 informationHST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007
MIT OpenCourseWare http://ocw.mit.edu HST.582J / 6.555J / 16.456J Biomedical Signal and Image Processing Spring 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationSegmenting Glioma in Multi-Modal Images using a Generative Model for Brain Lesion Segmentation
Segmenting Glioma in Multi-Modal Images using a Generative Model for Brain Lesion Segmentation Bjoern H. Menze 1,2, Koen Van Leemput 3, Danial Lashkari 4 Marc-André Weber 5, Nicholas Ayache 2, and Polina
More informationAutomatic Construction of 3D Statistical Deformation Models Using Non-rigid Registration
Automatic Construction of 3D Statistical Deformation Models Using Non-rigid Registration D. Rueckert 1, A.F. Frangi 2,3, and J.A. Schnabel 4 1 Visual Information Processing, Department of Computing, Imperial
More informationEMSegmenter Tutorial (Advanced Mode)
EMSegmenter Tutorial (Advanced Mode) Dominique Belhachemi Section of Biomedical Image Analysis Department of Radiology University of Pennsylvania 1/65 Overview The goal of this tutorial is to apply the
More informationSlicer3 Training Tutorial Using EM Segmenter with Non- Human Primate Images
Slicer3 Training Compendium Slicer3 Training Tutorial Using EM Segmenter with Non- Human Primate Images Vidya Rajagopalan Christopher Wyatt BioImaging Systems Lab Dept. of Electrical Engineering Virginia
More informationAutomatic 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 informationRegistration Using Sparse Free-Form Deformations
Registration Using Sparse Free-Form Deformations Wenzhe Shi 1, Xiahai Zhuang 2, Luis Pizarro 1,WenjiaBai 1, Haiyan Wang 1, Kai-Pin Tung 1, Philip Edwards 1, and Daniel Rueckert 1 1 Biomedical Image Analysis
More informationSubcortical Structure Segmentation using Probabilistic Atlas Priors
Subcortical Structure Segmentation using Probabilistic Atlas Priors Sylvain Gouttard 1, Martin Styner 1,2, Sarang Joshi 3, Brad Davis 2, Rachel G. Smith 1, Heather Cody Hazlett 1, Guido Gerig 1,2 1 Department
More informationTopology Preserving Brain Tissue Segmentation Using Graph Cuts
Topology Preserving Brain Tissue Segmentation Using Graph Cuts Xinyang Liu 1, Pierre-Louis Bazin 2, Aaron Carass 3, and Jerry Prince 3 1 Brigham and Women s Hospital, Boston, MA 1 xinyang@bwh.harvard.edu
More informationThe Insight Toolkit. Image Registration Algorithms & Frameworks
The Insight Toolkit Image Registration Algorithms & Frameworks Registration in ITK Image Registration Framework Multi Resolution Registration Framework Components PDE Based Registration FEM Based Registration
More informationCurvature guided surface registration using level sets
Curvature guided surface registration using level sets Marcel Lüthi, Thomas Albrecht, Thomas Vetter Department of Computer Science, University of Basel, Switzerland Abstract. We present a new approach
More informationCoupled Bayesian Framework for Dual Energy Image Registration
in Proc. IEEE Conf. on Computer Vision and Pattern Recognition (CVPR), Vol. 2, pp 2475-2482, 2006 Coupled Bayesian Framework for Dual Energy Image Registration Hao Wu University of Maryland College Park,
More informationAtlas of Classifiers for Brain MRI Segmentation
Atlas of Classifiers for Brain MRI Segmentation B. Kodner 1,2, S. H. Gordon 1,2, J. Goldberger 3 and T. Riklin Raviv 1,2 1 Department of Electrical and Computer Engineering, 2 The Zlotowski Center for
More informationLilla Zöllei A.A. Martinos Center, MGH; Boston, MA
Lilla Zöllei lzollei@nmr.mgh.harvard.edu A.A. Martinos Center, MGH; Boston, MA Bruce Fischl Gheorghe Postelnicu Jean Augustinack Anastasia Yendiki Allison Stevens Kristen Huber Sita Kakonoori + the FreeSurfer
More informationNorbert 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 informationAutomatic 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 informationDeformable Registration of Cortical Structures via Hybrid Volumetric and Surface Warping
Deformable Registration of Cortical Structures via Hybrid Volumetric and Surface Warping Tianming Liu, Dinggang Shen, and Christos Davatzikos Section of Biomedical Image Analysis, Department of Radiology,
More informationNorbert 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Überatlas: Robust Speed-Up of Feature-Based Registration and Multi-Atlas Based Segmentation. Jennifer Alvén
Überatlas: Robust Speed-Up of Feature-Based Registration and Multi-Atlas Based Segmentation Jennifer Alvén January 2015 Abstract Registration is a key component in multi-atlas approaches to medical image
More informationA Generic Framework for Non-rigid Registration Based on Non-uniform Multi-level Free-Form Deformations
A Generic Framework for Non-rigid Registration Based on Non-uniform Multi-level Free-Form Deformations Julia A. Schnabel 1, Daniel Rueckert 2, Marcel Quist 3, Jane M. Blackall 1, Andy D. Castellano-Smith
More informationRegistration Techniques
EMBO Practical Course on Light Sheet Microscopy Junior-Prof. Dr. Olaf Ronneberger Computer Science Department and BIOSS Centre for Biological Signalling Studies University of Freiburg Germany O. Ronneberger,
More informationKnowledge-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 informationDeep Similarity Learning for Multimodal Medical Images
Deep Similarity Learning for Multimodal Medical Images Xi Cheng, Li Zhang, and Yefeng Zheng Siemens Corporation, Corporate Technology, Princeton, NJ, USA Abstract. An effective similarity measure for multi-modal
More information3D Registration based on Normalized Mutual Information
3D Registration based on Normalized Mutual Information Performance of CPU vs. GPU Implementation Florian Jung, Stefan Wesarg Interactive Graphics Systems Group (GRIS), TU Darmstadt, Germany stefan.wesarg@gris.tu-darmstadt.de
More informationAdaptive Local Multi-Atlas Segmentation: Application to Heart Segmentation in Chest CT Scans
Adaptive Local Multi-Atlas Segmentation: Application to Heart Segmentation in Chest CT Scans Eva M. van Rikxoort, Ivana Isgum, Marius Staring, Stefan Klein and Bram van Ginneken Image Sciences Institute,
More informationScene-Based Segmentation of Multiple Muscles from MRI in MITK
Scene-Based Segmentation of Multiple Muscles from MRI in MITK Yan Geng 1, Sebastian Ullrich 2, Oliver Grottke 3, Rolf Rossaint 3, Torsten Kuhlen 2, Thomas M. Deserno 1 1 Department of Medical Informatics,
More informationChapter 3 Image Registration. Chapter 3 Image Registration
Chapter 3 Image Registration Distributed Algorithms for Introduction (1) Definition: Image Registration Input: 2 images of the same scene but taken from different perspectives Goal: Identify transformation
More informationWhite Matter Lesion Segmentation (WMLS) Manual
White Matter Lesion Segmentation (WMLS) Manual 1. Introduction White matter lesions (WMLs) are brain abnormalities that appear in different brain diseases, such as multiple sclerosis (MS), head injury,
More informationFast Free-Form Deformation using the Normalised Mutual Information gradient and Graphics Processing Units
Fast Free-Form Deformation using the Normalised Mutual Information gradient and Graphics Processing Units Marc Modat 1, Zeike A. Taylor 1, Josephine Barnes 2, David J. Hawkes 1, Nick C. Fox 2, and Sébastien
More informationLow-Rank to the Rescue Atlas-based Analyses in the Presence of Pathologies
Low-Rank to the Rescue Atlas-based Analyses in the Presence of Pathologies Xiaoxiao Liu 1, Marc Niethammer 2, Roland Kwitt 3, Matthew McCormick 1, and Stephen Aylward 1 1 Kitware Inc., USA 2 University
More informationBildverarbeitung für die Medizin 2007
Bildverarbeitung für die Medizin 2007 Image Registration with Local Rigidity Constraints Jan Modersitzki Institute of Mathematics, University of Lübeck, Wallstraße 40, D-23560 Lübeck 1 Summary Registration
More informationDeformable MRI-Ultrasound Registration via Attribute Matching and Mutual-Saliency Weighting for Image-Guided Neurosurgery
Deformable MRI-Ultrasound Registration via Attribute Matching and Mutual-Saliency Weighting for Image-Guided Neurosurgery Inês Machado 1,2(&), Matthew Toews 3, Jie Luo 1,4, Prashin Unadkat 5, Walid Essayed
More informationDistance Transforms in Multi Channel MR Image Registration
Distance Transforms in Multi Channel MR Image Registration Min Chen 1, Aaron Carass 1, John Bogovic 1, Pierre-Louis Bazin 2 and Jerry L. Prince 1 1 Image Analysis and Communications Laboratory, 2 The Laboratory
More informationNon-Rigid Registration of Medical Images: Theory, Methods and Applications
Non-Rigid Registration of Medical Images: Theory, Methods and Applications Daniel Rueckert Paul Aljabar Medical image registration [1] plays an increasingly important role in many clinical applications
More informationNon-Rigid Registration of Medical Images: Theory, Methods and Applications
Non-Rigid Registration of Medical Images: Theory, Methods and Applications Daniel Rueckert Paul Aljabar Medical mage registration [1] plays an increasingly important role in many clinical applications
More informationMRI 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 informationMarkov Networks in Computer Vision
Markov Networks in Computer Vision Sargur Srihari srihari@cedar.buffalo.edu 1 Markov Networks for Computer Vision Some applications: 1. Image segmentation 2. Removal of blur/noise 3. Stereo reconstruction
More information