Correction of Partial Volume Effects in Arterial Spin Labeling MRI
|
|
- Maurice Matthews
- 5 years ago
- Views:
Transcription
1 Correction of Partial Volume Effects in Arterial Spin Labeling MRI By: Tracy Ssali Supervisors: Dr. Keith St. Lawrence and Udunna Anazodo Medical Biophysics 3970Z Six Week Project April 13 th 2012
2 Introduction The amount of blood being delivered to brain tissue, also known as brain perfusion, is indicative of function within a given region. By measuring the amount of perfusion to each area, it is possible to map out which areas of the brain correspond to various actions. Using these observations, we can discover how the brain responds to different stimuli. In addition, we can gain insight on which areas of the brain are under or over stimulated in some disease states, by comparing the results of healthy subjects and patients. There are a number of different techniques which can be used to measure perfusion in the brain. In this experiment, perfusion will be measured using Arterial Spin Labeling MRI (ASL-MRI). The brain consists of grey matter (GM), white matter (WM) and cerebral spinal fluids (CSF). Since the majority of the diseases studied show irregular flow in grey matter (e.g. Alzheimer s and Parkinson s disease), most ASL-MRI studies are interested in grey matter perfusion [1]. Due to the resolution of ASL imaging, it is difficult to separate the signals coming from the three types of tissue. To give a better approximation of the perfusion signal coming from the grey matter tissue, a partial volume effects (PVE) correction is applied to the images. The objective of this experiment was to determine the signal to noise (SNR) in whole brain signal before and after the application of an in-house implementation of the PVE correction.
3 Theory Arterial Spin Labeling MRI Perfusion in the brain can be measured using diffusible tracers such as xenon in Perfusion CT or gadolinium DTPA in Dynamic Contrast enhanced MRI. The use of these established techniques often require the administration of radioactive or exogenous tracers which cannot be used on certain patient populations, and require long clearance times. Arterial Spin Labeling MRI is a novel technique which uses the magnetization of water molecules in arterial blood as a tracer to measure tissue perfusion non-invasively. ASL MRI consists of two experiments. The first experiment creates a control image while the second experiment creates a tagged image. To create a tagged image, arterial blood water is magnetically labeled by directing radiofrequency pulses at an area proximal to the desired imaging region. For perfusion images of the brain, the tag pulses are applied to vessels in the neck. A delay time is given to allow the labeled blood to reach the brain. When the labeled arterial water interacts with the magnetic field, the longitudinal magnetization is affected which in turn affects the signal being produced. The control image is taken to account for the static magnetization that is produced by the water molecules already present in the brain. Since we are interested in the change of flow, the difference of the signal produced from the labeled and the control image is taken. The difference of the control and the label image is known as the ASL signal change.
4 Figure 1: Schematic of how the ASL image is created. The first 3 images represent an imaged voxel. The difference of the control and the label image is taken to get the ASL signal, shown in the third image. Areas where the signal changes between the control and the label are highlighted in red. The rightmost image is a CBF image. It is proportional to the signals produced in the ASL image. [2] The fractional ASL signal change is proportional to the amount of perfusion within a tissue. This is expressed as a ratio between the difference of the control and the labeled signal in the tissue (ΔM), scaled by the fully relaxed tissue signal (M 0 ) i.e. (ΔM/M 0 ). The images are acquired in quick succession using Echo Planar Imaging (EPI). EPI has the fastest acquisition of images at the expense of having a lower spatial resolution, hence, the heterogeneity of the signals in each voxel. The partial volume corrected fractional signal change can be expressed as: [a] ( ) ( ) ( ) Where; P GM, P WM, and P CSF represent the probability maps of GM, WM and CSF respectively. ΔM GM, ΔM WM and ΔM CSF represent the ASL perfusion weighted signal of the gray matter, white matter and CSF. M 0_GM, M 0_WM, M 0_CSF represent the fully relaxed EPI magnetization in gray matter, white matter and CSF. The perfusion weighted contribution
5 from CSF is set to zero since the labeled and the control magnetizations in CSF are essentially the same [1]. The SNR can be estimated by dividing the mean ASL signal by the standard deviation or variance of the signal. Since the primary focus is on grey matter, the SNR is expressed as: [b] Where; µm GM is the mean ASL signal and σ GM is the standard deviation of the grey matter signal. In an image, the mean signal is represented by the mean value of all the voxels at a particular point. The standard deviation is the square root of the variance in these voxels. Partial Volume Effects The intensity of signals acquired from ASL MRI are inherently low. Consequently, multiple tag and control pairs are collected and averaged to increase the SNR. In the case of where the brain is atrophic, the resolution of ASL-EPI (which is approximately 3x3x3mm per voxel) is not small enough to distinguish the different tissue types. Additionally, the problem of point spread blurring is also often run into due to the rapid image acquisition. As a result of these two tendencies, the contributions from GM WM and CSF are combined. Since perfusion in the three types of tissues are not the same, a calculated difference in the flow between two voxels may be caused by the fact that the
6 voxels are heterogeneous rather than an actual difference in the flow. This is known as the partial volume effect. Prior to quantifying blood perfusion, ASL difference images are usually smoothened with a Gaussian filter to reduce the amount of noise by decreasing the variance. GM flow can then be estimated by setting a threshold value such that brain tissues that fall within the specified threshold value are classified accordingly. Generally, a threshold value of 0.8 is set for GM matter tissue [1]. All voxels which have values greater than 0.8 are used in the estimation of GM flow. This however, does not work well in brain areas that have voxels with high heterogeneity as it is prone to signal cross contamination (the combining of the different signals). By using this thresholding method instead of applying a true PVE correction, GM signal can be underestimated by up to 24% [2]. Kernel Regression Algorithm Instead of using the threshold to classify the tissue, a kernel regression algorithm can be used to account for partial volume effect. For this project, a PVE correction was implemented with, a kernel regression algorithm outlined by Iris Asllaniet. al [1]. The regression algorithm estimates the partial signal contribution of each tissue type in the ASL image based on MRI tissue contrast information from high resolution anatomical (MPRAGE) MRI image volume. The Kernel regression algorithm estimates the value of the GM, WM and CSF ASL signal contribution from voxels within the radius of a predefined kernel size. Increasing the kernel size will increase the radius of voxel values included as nearest neighbor points in the calculation. The regression algorithm uses the assumption that the GM values (or WM or CSF) are homogeneous over the kernel size.
7 As a result, the kernel size or k-fliter needs to be chosen precisely. Large kernel sizes tend to smoothen out the image too much whereas small kernel sizes are sensitive to noise. Methods Image Acquisition 5 Chronic Regional Pain syndrome (CRPS) patients were scanned on a Siemens 3.0T system. For each subject a high-resolution T1-weighted 3D MPRAGE anatomical image volume, 8 2D EPI-M0 images, and 64 label/control 2D EPI-ASL image pairs were acquired. Image Preprocessing For each subject, the skull was removed from the raw MPRAGE image volume. Using SPM8 image analysis software (FIL, UCL, London, UK), the MPRAGE image volume was segmented into GM, WM and CSF probability maps. The raw EPI images were preprocessed in SPM8 and in-house code written in MATLAB (The MathWorks, Natick, MA, USA). To remove motion errors, all the EPI images were aligned to the first acquired image. Average perfusion-weighted images were generated by time averaging the EPI-ASL images and taking the difference of the image pairs. Since the spatial resolution of the ASL MRI is coarser than the anatomical MRI images, the anatomical MRI images are down sampled using a coregistration algorithm to match the ASL MRI images.
8 PVE Correction The PVE correction algorithm was applied to the motion corrected EPI image time series of 5 CRPS subjects, to estimate partial ΔM and M 0 signal contribution for GM, WM and CSF. This experiment was carried out for a kernel size of 5 and 9. The SNR was estimated for ΔM in gray matter using formula [b] and compared to ASL data pre and post-pve correction to see whether there is a significant improvement in the signal to noise after the PVE correction is applied. Using an independent two sample t-test the results from the SNR calculation pre and post-correction were compared once for kernel size 5 and again for kernel size 9. K-filter The kernel size used in the regression model was adjusted to equal a voxel size of 5, 7, 9, 11 and 15. Using the different kernel sizes, the PVE correction was applied on the ASL images of an arbitrarily chosen subject. For each kernel size, the SNR was calculated. Results Over all, the uncorrected SNR calculated from the data is within the range reported by the literature. In an experiment conducted by Y.Chen et al, the reported average raw pcasl SNR was 1.58 [3]. The SNR for Subjects 1 through 5 were 1.27, 1.12, 1.33, 1.52 and 1.40 respectively. This yields an average of 1.33.
9 Signal to Noise Ratio PVE correction at k=5 for 5 CRPS patients Following the PVE correction of the images using a kernel size of 5, there was a significant decrease in SNR (P<0.05). The signal was consistently reduced by over 80% for each of the patients. These results are depicted in figure 2. 2 SNR before and after PV correction k = Uncorrected SNR PV corrected SNR 0 Subject 1 Subject 2 Subject 3 Subject 4 Subject 5 Complex Regional Pain Syndrome Patient Figure 2: Graphical representation of SNR before and after the PV correction was applied for the 5 CRPS patients using a kernel size of 5. The dotted bars represent the uncorrected SNR and the line patterned bars represent the PVE corrected SNR. PVE correction at k=9 for 5 CRPS patients In the second experiment, the kernel size was changed from 5 to 9. As seen in figure 3, the SNR after the PV correction yielded inconsistent results for each of the patients. The SNR of the ASL MRI images increased for subject 5, marginally changed for subjects 1 and 2 and considerably decreased for subjects 3 and 4. The results of the student t-test indicate that there was not a significant change in the SNR (P<0.05) after the application of the PVE correction. If the data for subject 5 is removed, there is a significant decrease in the SNR (P<0.1).
10 Signal to noise Ratio Signal to Noise Ratio SNR before and after PV correction k = Subject 1 Subject 2 Subject 3 Subject 4 Subject 5 Complex Regional Pain Syndrome Patient Uncorrected SNR PV corrected SNR Figure 3: Graphical representation of the SNR before and after the PV correction for 5 CPRS patients using a kernel size of 9. The dotted bars represent the uncorrected SNR and the line patterned bars represent the PV corrected SNR. The error bars represent the standard error of the data set. Kernel size and SNR of Subject 5 In the third experiment the k-filter was set equal to a voxel size of 5, 7, 9, 11, and 15. The results, shown in figure 4, show that as the kernel size increased, the SNR increased. The uncorrected SNR remains the same since there was no filters or corrections applied to it. Kernel Size vs SNR Kernel Size Uncorrected SNR PV Corrected SNR Figure 4: Graphical representation of the change in the SNR as the Kernel size changes. The dotted bars represent the Uncorrected SNR and the line patterned bars represent the PV corrected SNR.
11 Discussion PVE Correction at K=5 The results of the PVE correction experiment yielded unexpected results. Using the kernel size of 5, it was found that the SNR decreased considerably for each of the patients. A possible explanation for this trend is that the ASL MRI images had a low resolution and were noisy. The kernel algorithm works by reassigning a value to the pixel at the center based on the pixels surrounding it. If there is a lot of noise in the image, the filter will corrupt the data by assigning incorrect values to the voxels and possibly amplifying the variation of the voxels. Since the SNR is the signal divided by the noise, and the noise is the square root of the variance, an increase in variance will cause the SNR to decrease. In addition, if there is a substantial amount of noise in the image, the voxels which are being used in the regression algorithm cannot be considered homogenous. This is one of the assumptions that are made in using the regression algorithm. The SNR after the PV correction was so low that using it as an image processing technique is not beneficial. In this scenario, it may be practical to use the crude thresholding technique that was mentioned in the theory. PVE Correction at K=9 In the second experiment, the kernel size was increased to 9 and applied to all 5 subjects to test if a larger kernel size would provide a more accurate estimate the partial
12 signal contribution. The idea was that increasing the kernel filter will include more voxels in the algorithm and resultantly minimize the effect of noise. The results from this investigation varied from subject to subject; some of the SNR s increased, while the remaining ones either decreased or marginally changed. In terms of fitting the data, a k-filter of 9 is more stable. While this was an improvement from the k- filter of 5, the inconsistency of the results suggests that there is an underlying problem. One of the possibilities could be that there are bugs within the in-house PVE algorithm code that need to be investigated. Initially the code was intended to find the SNR voxel wise for each of the horizontal slices as well as the SNR for the entire brain. After reviewing the code, it appeared that it was calculating the SNR voxel wise vertically for each slice. With this in mind, there may have been other issues within the code that need to be ironed out. The student t-test of this data yielded insignificant results. This suggests that there is no difference from the mean signal before the PVE correction for this patient population. By looking at Figure 4, subject 5 had a very high SNR. This produced a relatively large increase in the P value which affected its significance. Kernel Filter and SNR Since the k-filter dictates the number of pixel elements which were used to determine the value of the signal, a change in the kernel size had a substantial impact on the SNR.
13 PVE corrected images which used small kernel sizes had lower SNRs while larger kernel sizes yielded higher SNRs. The low SNR could be attributed to the fact that the algorithm is unstable at lower kernel sizes. Since the algorithm is looking at a smaller radius, the variance within that area is higher. This especially happens at the edges between the background and the feature. Smaller kernel sizes tend to have lower SNR, take less time to compute and produce less smoothened images. The high SNR of the large kernel sizes can be explained by the larger region of voxels the algorithm uses. The algorithm works with a larger region and thus makes the entire image more uniform. After the PVE correction, the images are smoothened. As a result the SNRs tend to be greatly over estimated. In comparison to the small kernel sizes, the larger ones take less time to compute. As previously mentioned, it is important to use a kernel which doesn t smooth the data too much or break down as a result of variation due to noise. From analyzing the experimental data in figures 2, 3 and 4 it is difficult to decide on a kernel size which will optimize the SNR. A kernel size of 5 tends to result in very low SNRs. A kernel size of 9 yields inconsistent results or a significant decrease in the SNR (when patient subject 5 is removed). Above a kernel size of 9, the images are smoothed. Further tests are required to get a better idea of which kernel size would be most suitable.
14 Conclusion Based on the results from this experiment it would be difficult to justify the use of the PVE correction we implemented. In ASL MRI imaging it is common to underestimate the cerebral blood flow (CBF) as a result of partial volume effects. From the experimental data collected, the use of this PVE correction at a kernel size of 5, results in a decreased SNR (P<0.05) which would lead to an underestimation of the CBF. This could lead to misdiagnosis of both healthy and diseased patients. Using a kernel size of 9 would also lead to detrimental effects because the results were inconsistent across the patient population. As a result of the variance in the SNRs, there was no significant change in the SNR (P>0.05). By increasing the kernel sizes, we discovered that the SNR and the smoothness of the images increase. Resultantly, the kernel size has to be selected with this in mind. References [1] Asllani, Iris, AjnaBorogovac, and Truman R. Brown. "Regression Algorithm Correcting for Partial Volume Effects in Arterial Spin Labeling MRI." Magnetic Resonance in Medicine 60.6 (2008): [2] Asllani, Iris, AjnaBorogovac. "Arterial Spin Labeling (ASL) FMRI: Advantages, Theoretical Constrains and Experimental Challenges in Neurosciences." International Journal of Biomedical Imaging 2012 (2011). [3] Chen, Y., M. Korczykowski, and J. Wang. "Comparison of Reproducibility between Continuous, Pulsed, and Pseudo-continuous Arterial Spin Labeling." Magnetic Resonance in Medicine (2008): 2008.
15 [4] Liu, Thomas T., and Gregory G. Brown. "Measurement of Cerebral Perfusion with Arterial Spin Labeling: Part 1. Methods." Journal of the International Neuropsychological Society (2007). [5] Wolf, R., and J. Detre. "Clinical Neuroimaging Using Arterial Spin-Labeled Perfusion Magnetic Resonance Imaging." Neurotherapeutics 4.3 (2007):
MB-EPI PCASL. Release Notes for Version February 2015
MB-EPI PCASL Release Notes for Version 1.0 20 February 2015 1 Background High-resolution arterial spin labeling (ASL) imaging is highly desirable in both neuroscience research and clinical applications
More informationASAP_2.0 (Automatic Software for ASL Processing) USER S MANUAL
ASAP_2.0 (Automatic Software for ASL Processing) USER S MANUAL ASAP was developed as part of the COST Action "Arterial Spin Labelling Initiative in Dementia (AID)" by: Department of Neuroimaging, Institute
More informationA Spatio-temporal Denoising Approach based on Total Variation Regularization for Arterial Spin Labeling
A Spatio-temporal Denoising Approach based on Total Variation Regularization for Arterial Spin Labeling Cagdas Ulas 1,2, Stephan Kaczmarz 3, Christine Preibisch 3, Jonathan I Sperl 2, Marion I Menzel 2,
More informationA 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 informationSupplementary 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 informationClassification 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 informationPre-processing of ASL data T CT
Wed October 2, 2013 Image Processing Pre-processing: motion correction, denoising, outlier detection Alessandra Bertoldo Pre-processing of ASL data T CT C T C Single TI ASL T T T T C CCC average Pre-processing
More informationSupplementary Figure 1
Supplementary Figure 1 BOLD and CBV functional maps showing EPI versus line-scanning FLASH fmri. A. Colored BOLD and CBV functional maps are shown in the highlighted window (green frame) of the raw EPI
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 informationIntroduction to Neuroimaging Janaina Mourao-Miranda
Introduction to Neuroimaging Janaina Mourao-Miranda Neuroimaging techniques have changed the way neuroscientists address questions about functional anatomy, especially in relation to behavior and clinical
More informationChapter 3 Set Redundancy in Magnetic Resonance Brain Images
16 Chapter 3 Set Redundancy in Magnetic Resonance Brain Images 3.1 MRI (magnetic resonance imaging) MRI is a technique of measuring physical structure within the human anatomy. Our proposed research focuses
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 informationFmri 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 informationIntroduction to fmri. Pre-processing
Introduction to fmri Pre-processing Tibor Auer Department of Psychology Research Fellow in MRI Data Types Anatomical data: T 1 -weighted, 3D, 1/subject or session - (ME)MPRAGE/FLASH sequence, undistorted
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
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 informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationSPM8 for Basic and Clinical Investigators. Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
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 informationFunctional MRI data preprocessing. Cyril Pernet, PhD
Functional MRI data preprocessing Cyril Pernet, PhD Data have been acquired, what s s next? time No matter the design, multiple volumes (made from multiple slices) have been acquired in time. Before getting
More informationXI Conference "Medical Informatics & Technologies" VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD
XI Conference "Medical Informatics & Technologies" - 2006 medical imaging, MRI, brain perfusion Bartosz KARCZEWSKI 1, Jacek RUMIŃSKI 1 VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD Brain perfusion imaging
More informationThe 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 informationDeepASL: Kinetic Model Incorporated Loss for Denoising Arterial Spin Labeled MRI via Deep Residual Learning
DeepASL: Kinetic Model Incorporated Loss for Denoising Arterial Spin Labeled MRI via Deep Residual Learning Cagdas Ulas 1, Giles Tetteh 1, Stephan Kaczmarz 2, Christine Preibisch 2, and Bjoern H. Menze
More informationWhite Pixel Artifact. Caused by a noise spike during acquisition Spike in K-space <--> sinusoid in image space
White Pixel Artifact Caused by a noise spike during acquisition Spike in K-space sinusoid in image space Susceptibility Artifacts Off-resonance artifacts caused by adjacent regions with different
More informationDeepASL: Kinetic Model Incorporated Loss for Denoising Arterial Spin Labeled MRI via Deep Residual Learning
DeepASL: Kinetic Model Incorporated Loss for Denoising Arterial Spin Labeled MRI via Deep Residual Learning Cagdas Ulas 1, Giles Tetteh 1, Stephan Kaczmarz 2, Christine Preibisch 2, and Bjoern H. Menze
More informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Proceeding] Motion artifact correction in ASL images: automated procedure an improved Original Citation: Di Cataldo S., Ficarra E., Acquaviva A., Macii
More informationn o r d i c B r a i n E x Tutorial DSC Module
m a k i n g f u n c t i o n a l M R I e a s y n o r d i c B r a i n E x Tutorial DSC Module Please note that this tutorial is for the latest released nordicbrainex. If you are using an older version please
More informationRETROSPECTIVE TECHNIQUES FOR SEGMENTATION OF STRUCTURAL AND FUNCTIONAL MR BRAIN IMAGES
RETROSPECTIVE TECHNIQUES FOR SEGMENTATION OF STRUCTURAL AND FUNCTIONAL MR BRAIN IMAGES SURESH ANAND SADANANTHAN School of Computer Engineering A thesis submitted to the Nanyang Technological University
More informationAdvanced MRI Techniques (and Applications)
Advanced MRI Techniques (and Applications) Jeffry R. Alger, PhD Department of Neurology Ahmanson-Lovelace Brain Mapping Center Brain Research Institute Jonsson Comprehensive Cancer Center University of
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 informationAutomatic 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 informationComputational 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 informationStatistical Analysis of Neuroimaging Data. Phebe Kemmer BIOS 516 Sept 24, 2015
Statistical Analysis of Neuroimaging Data Phebe Kemmer BIOS 516 Sept 24, 2015 Review from last time Structural Imaging modalities MRI, CAT, DTI (diffusion tensor imaging) Functional Imaging modalities
More informationStatistical Analysis of MRI Data
Statistical Analysis of MRI Data Shelby Cummings August 1, 2012 Abstract Every day, numerous people around the country go under medical testing with the use of MRI technology. Developed in the late twentieth
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 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 informationfmri Image Preprocessing
fmri Image Preprocessing Rick Hoge, Ph.D. Laboratoire de neuroimagerie vasculaire (LINeV) Centre de recherche de l institut universitaire de gériatrie de Montréal, Université de Montréal Outline Motion
More informationCHAPTER 3 TUMOR DETECTION BASED ON NEURO-FUZZY TECHNIQUE
32 CHAPTER 3 TUMOR DETECTION BASED ON NEURO-FUZZY TECHNIQUE 3.1 INTRODUCTION In this chapter we present the real time implementation of an artificial neural network based on fuzzy segmentation process
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationSegmenting 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 informationSpatial 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 informationMRI. When to use What sequences. Outline 2012/09/19. Sequence: Definition. Basic Principles: Step 2. Basic Principles: Step 1. Govind Chavhan, MD
MRI When to use What sequences Govind Chavhan, MD Assistant Professor and Staff Radiologist The Hospital For Sick Children, Toronto Planning Acquisition Post processing Interpretation Patient history and
More informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More informationDiffusion MRI Acquisition. Karla Miller FMRIB Centre, University of Oxford
Diffusion MRI Acquisition Karla Miller FMRIB Centre, University of Oxford karla@fmrib.ox.ac.uk Diffusion Imaging How is diffusion weighting achieved? How is the image acquired? What are the limitations,
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 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 informationCHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd
CHAPTER 2 Morphometry on rodent brains A.E.H. Scheenstra J. Dijkstra L. van der Weerd This chapter was adapted from: Volumetry and other quantitative measurements to assess the rodent brain, In vivo NMR
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 informationDimension Reduction for Big Data Analysis. Dan Shen. Department of Mathematics & Statistics University of South Florida.
Dimension Reduction for Big Data Analysis Dan Shen Department of Mathematics & Statistics University of South Florida danshen@usf.edu October 24, 2014 1 Outline Multiscale weighted PCA for Image Analysis
More informationNew Technology Allows Multiple Image Contrasts in a Single Scan
These images were acquired with an investigational device. PD T2 T2 FLAIR T1 MAP T1 FLAIR PSIR T1 New Technology Allows Multiple Image Contrasts in a Single Scan MR exams can be time consuming. A typical
More informationMeasuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI
Measuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI Revision date: 09/15/2008 Overview This document describes the procedure for measuring baseline whole-brain perfusion
More informationCP Generalize Concepts in Abstract Multi-dimensional Image Model Component Semantics. David Clunie.
CP-1390 - Generalize Concepts in Abstract Multi-dimensional Image Model Semantics Page 1 STATUS Date of Last Update Person Assigned Submitter Name Submission Date Assigned 2014/06/09 David Clunie mailto:dclunie@dclunie.com
More informationWhole 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 informationEmpirical analyses of null-hypothesis perfusion FMRI data at 1.5 and 4 T
NeuroImage 19 (2003) 1449 1462 www.elsevier.com/locate/ynimg Empirical analyses of null-hypothesis perfusion FMRI data at 1.5 and 4 T Jiongjiong Wang, a,b, * Geoffrey K. Aguirre, a,c Daniel Y. Kimberg,
More informationParameter estimation in arterial spin labeling MRI: comparing the four phase model and the buxton model with fourier transform
Original Article Parameter estimation in arterial spin labeling MRI: comparing the four phase model and the buxton model with fourier transform Jim Ji 1, Vincent Pham 1, Xiao-Ping Zhu 2, Ka-loh Li 2 1
More informationINDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS
INDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS V. Calhoun 1,2, T. Adali, 2 and G. Pearlson 1 1 Johns Hopkins University Division of Psychiatric Neuro-Imaging,
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 informationHistograms. 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 information6 credits. BMSC-GA Practical Magnetic Resonance Imaging II
BMSC-GA 4428 - Practical Magnetic Resonance Imaging II 6 credits Course director: Ricardo Otazo, PhD Course description: This course is a practical introduction to image reconstruction, image analysis
More informationCHAPTER 9: Magnetic Susceptibility Effects in High Field MRI
Figure 1. In the brain, the gray matter has substantially more blood vessels and capillaries than white matter. The magnified image on the right displays the rich vasculature in gray matter forming porous,
More informationStructural 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 informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationfmri pre-processing Juergen Dukart
fmri pre-processing Juergen Dukart Outline Why do we need pre-processing? fmri pre-processing Slice time correction Realignment Unwarping Coregistration Spatial normalisation Smoothing Overview fmri time-series
More informationAnalysis of fmri data within Brainvisa Example with the Saccades database
Analysis of fmri data within Brainvisa Example with the Saccades database 18/11/2009 Note : All the sentences in italic correspond to informations relative to the specific dataset under study TP participants
More informationEMSegment Tutorial. How to Define and Fine-Tune Automatic Brain Compartment Segmentation and the Detection of White Matter Hyperintensities
EMSegment Tutorial How to Define and Fine-Tune Automatic Brain Compartment Segmentation and the Detection of White Matter Hyperintensities This documentation serves as a tutorial to learn to customize
More informationMeasuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI
Measuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI Revision date: 11/20/2006 Overview This document describes the procedure for measuring baseline whole-brain perfusion
More informationTumor Detection and classification of Medical MRI UsingAdvance ROIPropANN Algorithm
International Journal of Engineering Research and Advanced Technology (IJERAT) DOI:http://dx.doi.org/10.31695/IJERAT.2018.3273 E-ISSN : 2454-6135 Volume.4, Issue 6 June -2018 Tumor Detection and classification
More informationMEDICAL 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 informationQUANTITATION OF THE PREMATURE INFANT BRAIN VOLUME FROM MR IMAGES USING WATERSHED TRANSFORM AND BAYESIAN SEGMENTATION
QUANTITATION OF THE PREMATURE INFANT BRAIN VOLUME FROM MR IMAGES USING WATERSHED TRANSFORM AND BAYESIAN SEGMENTATION Merisaari Harri 1;2, Teräs Mika 2, Alhoniemi Esa 1, Parkkola Riitta 2;3, Nevalainen
More informationGLM for fmri data analysis Lab Exercise 1
GLM for fmri data analysis Lab Exercise 1 March 15, 2013 Medical Image Processing Lab Medical Image Processing Lab GLM for fmri data analysis Outline 1 Getting Started 2 AUDIO 1 st level Preprocessing
More informationTUMOR DETECTION IN MRI IMAGES
TUMOR DETECTION IN MRI IMAGES Prof. Pravin P. Adivarekar, 2 Priyanka P. Khatate, 3 Punam N. Pawar Prof. Pravin P. Adivarekar, 2 Priyanka P. Khatate, 3 Punam N. Pawar Asst. Professor, 2,3 BE Student,,2,3
More informationLearning to Identify Fuzzy Regions in Magnetic Resonance Images
Learning to Identify Fuzzy Regions in Magnetic Resonance Images Sarah E. Crane and Lawrence O. Hall Department of Computer Science and Engineering, ENB 118 University of South Florida 4202 E. Fowler Ave.
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
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 informationGlobal Thresholding Techniques to Classify Dead Cells in Diffusion Weighted Magnetic Resonant Images
Global Thresholding Techniques to Classify Dead Cells in Diffusion Weighted Magnetic Resonant Images Ravi S 1, A. M. Khan 2 1 Research Student, Department of Electronics, Mangalore University, Karnataka
More informationmritc: 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 informationUser Manual Parameter Estimation for a Compartmental Tracer Kinetic Model Applied to PET data
User Manual Parameter Estimation for a Compartmental Tracer Kinetic Model Applied to PET data Guadalupe Ayala 1 Christina Negoita 2 Rosemary A. Renaut 3 May 12, 2005 1 This work was supported in part by
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
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 informationAppendix E1. Supplementary Methods. MR Image Acquisition. MR Image Analysis
RSNA, 2015 10.1148/radiol.2015150532 Appendix E1 Supplementary Methods MR Image Acquisition By using a 1.5-T system (Avanto, Siemens Medical, Erlangen, Germany) under a program of regular maintenance (no
More informationDeviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies
g Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies Presented by Adam Kesner, Ph.D., DABR Assistant Professor, Division of Radiological Sciences,
More informationImproved Quantification of Arterial Spin Labelling Images using Partial Volume Correction Techniques. Ruth Abigail Oliver
Improved Quantification of Arterial Spin Labelling Images using Partial Volume Correction Techniques Ruth Abigail Oliver A thesis submitted to University College London for the degree of Doctor of Philosophy
More informationJournal of Articles in Support of The Null Hypothesis
Data Preprocessing Martin M. Monti, PhD UCLA Psychology NITP 2016 Typical (task-based) fmri analysis sequence Image Pre-processing Single Subject Analysis Group Analysis Journal of Articles in Support
More informationBasic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis
Basic principles of MR image analysis Basic principles of MR image analysis Julien Milles Leiden University Medical Center Terminology of fmri Brain extraction Registration Linear registration Non-linear
More informationA Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data
A Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data Alexey Samsonov, Julia Velikina Departments of Radiology and Medical
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 informationMEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 9: Medical Image Segmentation (III) (Fuzzy Connected Image Segmentation)
SPRING 2017 1 MEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 9: Medical Image Segmentation (III) (Fuzzy Connected Image Segmentation) Dr. Ulas Bagci HEC 221, Center for Research in Computer Vision (CRCV),
More informationCocozza S., et al. : ALTERATIONS OF FUNCTIONAL CONNECTIVITY OF THE MOTOR CORTEX IN FABRY'S DISEASE: AN RS-FMRI STUDY
ALTERATIONS OF FUNCTIONAL CONNECTIVITY OF THE MOTOR CORTEX IN FABRY'S DISEASE: AN RS-FMRI STUDY SUPPLEMENTARY MATERIALS Sirio Cocozza, MD 1*, Antonio Pisani, MD, PhD 2, Gaia Olivo, MD 1, Francesco Saccà,
More informationExisting Packages. Overview. Siemens. Michael A. Chappell
Existing Packages Michael A. Chappell michael.chappell@eng.ox.ac.uk www.ibme.ox.ac.uk/qubic Institute of Biomedical Engineering & Oxford Centre for Functional MRI of the Brain University of Oxford. Scanner
More informationarxiv: v2 [q-bio.qm] 16 Oct 2017
Gulban this is page 1 The relation between color spaces and compositional data analysis demonstrated with magnetic resonance image processing applications O.F. Gulban Maastricht University, Maastricht,
More informationStructural 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 informationQuantitative MRI of the Brain: Investigation of Cerebral Gray and White Matter Diseases
Quantities Measured by MR - Quantitative MRI of the Brain: Investigation of Cerebral Gray and White Matter Diseases Static parameters (influenced by molecular environment): T, T* (transverse relaxation)
More informationCS/NEUR125 Brains, Minds, and Machines. Due: Wednesday, April 5
CS/NEUR125 Brains, Minds, and Machines Lab 8: Using fmri to Discover Language Areas in the Brain Due: Wednesday, April 5 In this lab, you will analyze fmri data from an experiment that was designed to
More informationPatch-Based Super-Resolution of Arterial Spin Labeling Magnetic Resonance Images
Patch-Based Super-Resolution of Arterial Spin Labeling Magnetic Resonance Images Cédric Meurée, Pierre Maurel, Jean-Christophe Ferré, Christian Barillot To cite this version: Cédric Meurée, Pierre Maurel,
More informationI.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger:
Supplementary Materials I. Evidence of sex differences on eating behavior in children I.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger: Table 2. Parent Report for Child
More informationSurface Projection Method for Visualizing Volumetric Data
Surface Projection Method for Visualizing Volumetric Data by Peter Lincoln A senior thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Science With Departmental Honors
More informationA novel dementia diagnosis strategy on arterial spin labeling magnetic resonance images via pixel-wise partial volume correction and ranking
Multimed Tools Appl (216) 75:267 29 DOI 1.17/s1142-14-2395-2 A novel dementia diagnosis strategy on arterial spin labeling magnetic resonance images via pixel-wise partial volume correction and ranking
More informationRADIOMICS: potential role in the clinics and challenges
27 giugno 2018 Dipartimento di Fisica Università degli Studi di Milano RADIOMICS: potential role in the clinics and challenges Dr. Francesca Botta Medical Physicist Istituto Europeo di Oncologia (Milano)
More informationNormalization for clinical data
Normalization for clinical data Christopher Rorden, Leonardo Bonilha, Julius Fridriksson, Benjamin Bender, Hans-Otto Karnath (2012) Agespecific CT and MRI templates for spatial normalization. NeuroImage
More information