CHAPTER 9: Magnetic Susceptibility Effects in High Field MRI

Similar documents
Functional MRI. Jerry Allison, Ph. D. Medical College of Georgia

Slide 1. Technical Aspects of Quality Control in Magnetic Resonance Imaging. Slide 2. Annual Compliance Testing. of MRI Systems.

Following on from the two previous chapters, which considered the model of the

MRI Physics II: Gradients, Imaging

White Pixel Artifact. Caused by a noise spike during acquisition Spike in K-space <--> sinusoid in image space

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008

Diffusion MRI Acquisition. Karla Miller FMRIB Centre, University of Oxford

(a Scrhon5 R2iwd b. P)jc%z 5. ivcr3. 1. I. ZOms Xn,s. 1E IDrAS boms. EE225E/BIOE265 Spring 2013 Principles of MRI. Assignment 8 Solutions

Role of Parallel Imaging in High Field Functional MRI

ACQUIRING AND PROCESSING SUSCEPTIBILITY WEIGHTED IMAGING (SWI) DATA ON GE 3.0T

COBRE Scan Information

CORRECTION OF IMAGE DISTORTION IN ECHO PLANAR IMAGE SERIES USING PHASE AND INTENSITY. Ning Xu. Dissertation. Submitted to the Faculty of the

MR Advance Techniques. Vascular Imaging. Class III

Lucy Phantom MR Grid Evaluation

Lab Location: MRI, B2, Cardinal Carter Wing, St. Michael s Hospital, 30 Bond Street

This exercise uses one anatomical data set (ANAT1) and two functional data sets (FUNC1 and FUNC2).

A more accurate account of the effect of k-space sampling and signal decay on the effective spatial resolution in functional MRI

surface Image reconstruction: 2D Fourier Transform

A Transmission Line Matrix Model for Shielding Effects in Stents

Supplementary Figure 1

MRI. When to use What sequences. Outline 2012/09/19. Sequence: Definition. Basic Principles: Step 2. Basic Principles: Step 1. Govind Chavhan, MD

Fmri Spatial Processing

MB-EPI PCASL. Release Notes for Version February 2015

Module 5: Dynamic Imaging and Phase Sharing. (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review. Improving Temporal Resolution.

SIEMENS MAGNETOM Verio syngo MR B17

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008

Scan Acceleration with Rapid Gradient-Echo

Parallel Imaging. Marcin.

Supplementary Information

Magnetic Resonance Elastography (MRE) of Liver Disease

Page 1 of 9. Protocol: adult_other_adni3basichumanprotocol25x_ _ _1. 3 Plane Localizer. 3 Plane Localizer PATIENT POSITION

Imaging Notes, Part IV

Dynamic Contrast enhanced MRA

Partial k-space Reconstruction

New Technology Allows Multiple Image Contrasts in a Single Scan

Partial k-space Recconstruction

Acquisition Methods for fmri at 7 Tesla

SIEMENS MAGNETOM Avanto syngo MR B15

SIEMENS MAGNETOM Skyra syngo MR D13

Sources of Distortion in Functional MRI Data

Image Acquisition Systems

Field Maps. 1 Field Map Acquisition. John Pauly. October 5, 2005

Basic fmri Design and Analysis. Preprocessing

M R I Physics Course

TOF-MRA Using Multi-Oblique-Stack Acquisition (MOSA)

XI Signal-to-Noise (SNR)

SIEMENS MAGNETOM Avanto syngo MR B15

3D z-shim Method for Reduction of Susceptibility Effects in BOLD fmri

Philips MRI Protocol Dump Created on Comment Software Stream

Head motion in diffusion MRI

Supplementary methods

Improved Spatial Localization in 3D MRSI with a Sequence Combining PSF-Choice, EPSI and a Resolution Enhancement Algorithm

Breast MRI Accreditation Program Clinical Image Quality Guide

Chapter 3 Set Redundancy in Magnetic Resonance Brain Images

Midterm Review

SPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing

Qualitative Comparison of Conventional and Oblique MRI for Detection of Herniated Spinal Discs

Transform Introduction page 96 Spatial Transforms page 97

Sampling, Ordering, Interleaving

Measuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI

Computational Medical Imaging Analysis

MRI image formation 8/3/2016. Disclosure. Outlines. Chen Lin, PhD DABR 3. Indiana University School of Medicine and Indiana University Health

Functional MRI in Clinical Research and Practice Preprocessing

Measuring baseline whole-brain perfusion on GE 3.0T using arterial spin labeling (ASL) MRI

IMAGE reconstruction in conventional magnetic resonance

ThE ultimate, INTuITIVE Mr INTErFAcE

The cost of parallel imaging in functional MRI of the human brain

SIEMENS MAGNETOM Verio syngo MR B15V

EPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing

2.1 Signal Production. RF_Coil. Scanner. Phantom. Image. Image Production

SPM8 for Basic and Clinical Investigators. Preprocessing

CSImage Tutorial v August 2000

A Virtual MR Scanner for Education

MRI Imaging Options. Frank R. Korosec, Ph.D. Departments of Radiology and Medical Physics University of Wisconsin Madison

Continuation Format Page

Unaliasing by Fourier-Encoding the Overlaps Using the Temporal Dimension (UNFOLD), Applied to Cardiac Imaging and fmri

Super-Resolution Reconstruction of Diffusion-Weighted Images from Distortion Compensated Orthogonal Anisotropic Acquisitions.

Basic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis

Compressed Sensing for Rapid MR Imaging

A Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data

Correction for EPI Distortions Using Multi-Echo Gradient-Echo Imaging

GE Healthcare CLINICAL GALLERY. Discovery * MR750w 3.0T. This brochure is intended for European healthcare professionals.

Fast Imaging Trajectories: Non-Cartesian Sampling (1)

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

Sampling, Ordering, Interleaving

Exam 8N080 - Introduction MRI

Blipped-Controlled Aliasing in Parallel Imaging for Simultaneous Multislice Echo Planar Imaging With Reduced g-factor Penalty

Advanced Imaging Trajectories

CS/NEUR125 Brains, Minds, and Machines. Due: Wednesday, April 5

Imaging Guide Pancreatic Imaging

ADNI, ADNI_QH, SURVEY. Geometry. connection

Functional MRI data preprocessing. Cyril Pernet, PhD

SIEMENS MAGNETOM TrioTim syngo MR B17

Three-dimensional numerical simulations of susceptibility-induced magnetic field inhomogeneities in the human head

45 µm polystyrene bead embedded in scattering tissue phantom. (a,b) raw images under oblique

Magnetic Susceptibility Artefact Correction of Spin-Echo and Gradient- Echo EPI Images

Steen Moeller Center for Magnetic Resonance research University of Minnesota

Advanced Image Reconstruction Methods for Photoacoustic Tomography

Transcription:

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, sponge-like structures. The high concentration of blood in the gray matter vasculature makes the gray matter highly sensitive to the blood oxygen leveldependent (BOLD) effect.

Figure 2. The major air-containing structures included in the model are the frontal sinus (turquoise), the sphenoid sinus (pink), the maxillary sinuses (blue), mastoid air cells (red), nasal and oral cavities, and trachea (yellow). Calculated (top) and measured (bottom) B 0 distributions demonstrate local gradients in the sagittal (a), coronal (b), and axial (c) planes. Magnetic field deviations are color coded in parts per million (ppm). Reprinted with permission from [30].

Figure 3. Calculated B 0 distributions inside the human head model demonstrate the remarkable distortion by the air-containing structures inside the human head. Illustration courtesy of Drs. Shizhe Li, C.M. Collins, Q.X. Yang, and M.B. Smith. Figure 4. Axial plane view of numerically calculated B 0 field magnitudes in ppm deviation from applied field. The gray scale is set to reveal the B 0 deviation due to the susceptibility difference between gray and white matter in the brain. Reprinted with permission from [35].

Figure 5. Timing diagram of 2D gradient echo sequence for a voxel. A, slice selection gradient; B, slice rewinding gradient. During TE the local gradient in the slice selection dimension, G z l, causes intravoxel phase dispersion, which results in signal loss in GE images. A compensation gradient, G c, or a tailored RF pulse, α[φ 0 (z)] can be applied to cancel the phase dispersion by G z l. Figure 6. The GESEPI image sets of the human brain obtained at 3.0 Tesla before (2D-FT data set) and after (3D-FT) a Fourier transform along the slice direction with respect to k z by the compensation gradient. The 2D-FT image set was obtained by sequentially incrementing G c to obtain N = 16 frames. Note the removal of the dark artifacts in the images from 4 6 sub-slices (ca. 1 mm). Reprinted with permission from [46]).

Figure 7. A coronal view of the phantom and images of the axial slice along the midline between the two spheres with the GE and GESEPI methods. The phantom consists of two air-filled spheres positioned inside a gelatin-filled cylindrical container with its longitudinal axis perpendicular to B 0. The voxels P, Q, and R indicated in the GE image are labeled P', Q', and R' in the corresponding GESEPI images. The phantom images were obtained with a flip angle of 25º, TR/TE = 80 ms/25 ms, FOV = 25 25 cm 2, image matrix = 128 128, slice thickness = 5 mm. Adapted with permission from [46]. Figure 8. Serial signal intensity plots from the representative voxels P, Q, and R (see Figure 7) in the 2D-FT data set, and from the corresponding voxels P', Q', and R' in the 3D-FT data set. The image number axis is shown above each plot, corresponding to either k z (2D- FT, top row) or position z along the slice direction (3D-FT, bottom row). The dashed lines in the top row plots indicate k = 0, where 2D images were acquired. Reprinted with permission from [46].

Figure 9. GESEPI effectively removed signal-loss artifacts in heavily T 2 *-weighted images at 3.0 and 7.0 Tesla. All the images have the same slice thickness 5 mm. The GESEPI images were obtained with 8 G c steps. The intensities of the globus pallidus (GP) and substantia nigra (SN) are obscured in the GE images and enhanced in the GESEPI images. Note that the signal nonuniformity in the GESEPI image at 7.0 Tesla is caused by the B 1 field distortion due to the wave effect. Figure 10. T 2 *-weighted images of a mouse brain at 14.0 Tesla. By removing the artifacts in the GE image, the GESEPI image improves contrast resolution and displays anatomic details not observable in spin-echo images. Reprinted with permission from [50].

Figure 11. GE (top row) and GESEPI (bottom row) images of an immature rat brain as a function of TE taken on a 9.4T system. TR = 80 ms, matrix = 128 128, FOV = 20 20 mm 2, slice thickness = 1 mm. Reprinted with permission from [ 46]. Figure 12. Average R 2 * map of the human brain with GESEPI.

Figure 13. The voxel signal intensity from the excited slice of the GESEPI data is plotted as a function of k c and acquisition time t along the trajectory (a) and its Fourier transform (b). The plots were calculated using Eqs. (13) and (14) with T 2 * = 30 ms, G z l = 0.5 ppm/cm, z 0 = 5 mm, and B 0 = 3.0 Tesla. The plots in (c) are standard 2D voxel signal intensity along the trajectory with k l = 0 (Plot I), k l 0 (Plot II), and GESEPI intensity with k l 0 (Plot III), respectively. The Plot II in (c), obtained by intersecting the two-dimensional plot in (a) with k z = 0, shows a sinc-function modulation of T 2 * decay. Plot III, obtained by intersecting the plot in (b) with z = 0, coincides with Plot I of normal T 2 * relaxation. Figure 14. Top row: the signal intensities acquired without phase-encoding in the left right direction. Center row: image intensity plot, representing point-spread function (PSF) profile in phase-encoding direction. Bottom row: images acquired with phaseencoding, demonstrating the reduction of blurring artifacts in EPI using the GESEPI method. The quasi-point source phantom for axial images consists of two 5-mm diameter tubes filled with gelatin, with a 1.8-mm distance between the tubes. The tubes are aligned with the B 0 field. The image intensities are adjusted to facilitate visualization of PSF. Adapted with permission from [71].

Figure 15. Spin-echo EPI images acquired without (a) and with SENSE with a reduction factor r = 3 (b) from an inferior brain slice. In spin echo images, the signal-loss artifact is minimal so that the geometric distortion due to the in-plane gradient can be demonstrated clearly. In order to highlight the brain geometry, a diffusion gradient was applied to reduce the CSF signals. Both images were acquired from a Philips Intera 3.0T MRI system with the same acquisition parameters. Reprinted with permission from [75].

Figure 16. The EPI images were acquired from the same human brain volume with SENSE (r = 3) using the GE (a), GESEPI (b), and SE (c) methods. The typical signal-loss artifact in the GE images (a) in the inferior frontal and temporal brain regions indicated by the arrows is reduced significantly with GESEPI (b). Reprinted with permission from [71].

Figure 17. Relationship of T 2 * decay with Forward and Reversed spiral acquisition. Illustration provided courtesy of Dr. V. Andrew Stenger. Figure 18. Simulated k-space signals of a point source with reversed and forward spiral acquisitions (a) and the corresponding image profiles (b). Illustration provided courtesy of Dr. V. Andrew Stenger.

Figure 19. Simulated (a) and experimental (b) image profiles with different data acquisition methods in the absence and presence of local gradients. Illustration provided courtesy of Dr. V. Andrew Stenger. Figure 20. Spiral images taken from the inferior brain region show remarkable signal-loss artifact. As indicated by the arrows, the artifact is mitigated by reversed spiral and removed when GESEPI is incorporated into the reversed spiral technique.

Figure 21. Olfactory fmri activation at orbitofrontal and olfactory tract with GESEPIreversed-spiral imaging.