Development of Phase Contrast Sequence
|
|
- Christopher Neal
- 5 years ago
- Views:
Transcription
1 École Polytechnique de Montréal Development of Phase Contrast Sequence Author : Dominique Mathieu Matricule Supervisor : Rick Hoge, PhD September 1 st, 2009
2 Contents 1 Acknowledgements 2 2 Project description 2 3 Working environment IDEA T Scanner Neurolens Theoretical concepts Effect of velocity on spin phase Flow compensation Phase contrast angiography Implementation FLASH.cpp sample Implementation of PCA sequence Preliminary results First scan session : Untriggered PCA acquisition Second scan session : PCA acquisition with ECG triggering Conclusion 13 8 References 14 1
3 1 Acknowledgements I would like to express gratitude to Professor Rick Hoge for this project offer and his integration to the lab team. Getting familiar with pulse programing has helped me understand practical aspects of MRI technique. Also, special thanks to Felipe Tancredi, master student at the Neurovascular Imaging Lab, for his constant support on both technical and theoretical aspects and his knowledge sharing in pulse sequence programming. 2 Project description L Unité de participation et d initiation à la recherche (UPIR) is an academic program from Polytechnique Montreal enabling undergraduate students to get familiar with research activities. Students are invited to collaborate a half day s work for 25 weeks with a Polytechnique teacher or anyone from the biomedical or nuclear institute. During this project, I had the opportunity to work with Rick Hoge s lab team. Professor Rick Hoge is the director of the Neurovascular Imaging Lab 1, a lab developing technology to image blood flow in the brain in order to study physiological regulatory processes during neural activation. The original goal of the project was to analyze the diffusional attenuation during global increases in cerebral blood flow. The objectives were to : 1. Use phase contrast angiography (PCA) sequence available to measure the blood flow in arteries and veins associated with primary motor cortex at rest and during execution of a motor task or during hypercapnia. 2. Perform numerical simulations with Matlab to predict effects of diffusion attenuation gradients on signal in such vessels. At the beginning of the project, we realized that the angiography package required to perform velocity measurements was not available on the scanner. The project goal was redefined and would then be to implement a PCA sequence using a FLASH template. Hence, the new project objectives became to : 1. Get familiar with MRI pulse sequence and object-oriented C++ development. 2. Implement PCA sequence from FLASH template. 3. Perform velocity measurements on the scanner. Table 1 presents in detail the time repartition of the project. Task description Hours Theoretical review 10 Familiarization with IDEA environment 25 Sequence programming 60 Scan test 5 Data analysis 5 Report writing 20 Total 125 Table 1: Time repartition 1 Neurovascular Imaging Lab Rick Hoge Ph.D. 2
4 3 Working environment The working environment provided the equipment required for sequence implementation, data acquisition and image processing. 3.1 IDEA IDEA stands for Integrated Development Environment for Applications. This framework was used for the development of MRI pulse sequence. It contains a shell environment used for sequence programming called the sequence development environment (SDE). This environment enables mutiple tasks such as 1. Sequence testing. 2. Protocol editing. 3. Timing table display. 4. etc. Several source code templates, for instance the FLASH sequence, were included in the software for development. These were useful as a starting point for the implementation of the PCA sequence T Scanner The images were acquired on a Siemens Trio system (3T). Scan tests were performed under the supervison of Carollyn Hurst, radiology technologist of the UNF. 3.3 Neurolens Neurolens is a visualization and analysis platform for quantitative physiological neuroimaging that provides fast image processing. Neurolens was used for flow quantification image analysis. 4 Theoretical concepts This section establishes the theoretical basis which is essential for a complete understanding of PCA sequence implementation. It describes the effects observed when imaging moving objects such as blood flow, and addresses the physical concepts under phase correction using the flow compensation technique. Finally, a discussion on PCA, which uses phase behavior of moving spin to encode velocities, is proposed. 4.1 Effect of velocity on spin phase The phase behavior of spins with constant velocity (as those found in a plug flow) exposed to a gradient magnetic field can be derived from the Newton s laws of motion. Let s define a spin with initial position, x 0, moving in the read encoding direction (ˆx) with constant velocity v x, exposed to a typical readout gradient waveform where 0 t 3τ as shown in Figure 1. x (t) = x 0 + v x t (1) During the application of a gradient magnetic field, Larmor precession becomes time dependent ω(t) = γb 0 + γgz (t). As a result, the phase accumulated by moving spin φ(t) = ω (t) dt varies as a quadratic function of time and doesn t come to zero at the echo time (t = 2τ). 3
5 This phase behavior during the rephasing lobe (τ t 3τ ) can be expressed in terms of stationary and velocity components (Figure 1) φ(t ) = φ s (t ) + φ v (t ) (2) where φ s (t ) and φ v (t ) are respectively position and velocity dependent φ s (t ) = γgx 0 t (3) with t = t 2τ and τ t τ. φ v (t ) = 1 2 γgv x ( t 2 + 4τt + 2τ 2) (4) Figure 1: Phase behavior of static and mobile spins during typical readout gradient application Solving Equation 2 at the echo time (t = 0) indicates that the phase acquired by moving spin is proportional to the velocity φ(t ) = φ v (0) = γgv x τ 2 (5) t =0 This demonstration can be extended to a laminar flow which can be decomposed over a short period of time and a tiny dot of space (pixel) as the sum of multiple plug flows at different velocities. To avoid non-zero phase at the echo, it is possible to add flow compensation gradients to the typical gradient waveform as discussed in the next section. 4.2 Flow compensation Flow compensation is a technique used to correct phase dispersion caused by motion or to reduce flow artifacts such as ghosting due to periodic flow. It is based on the principle of even-echo rephasing and consists in the addition of extra gradient lobes along the motion direction. Particular flow compensation waveform can be used to correct both constant velocity or accelerated spin. 4
6 One solution is here proposed to correct phase acquired by spins with constant velocity moving in the read encoding direction, as discussed in the previous section. For this scheme, two preparatory lobes G 1 and G 2 must be added prior to the readout gradient G. To facilitate the resolution, the duration τ of each three lobes is set (Figure 2). The third lobe added to the traditional readout gradients must ensure that phase is zero at echo time for both stationary and moving spins. For the stationary spins, this condition leads to the first equation on the zeroth moment M 0 of the gradient waveform G (t) which leads to M 0 (te) = ˆ te 0 G 1 t τ + G 2t 0 2τ τ G (t) dt = 0 (6) + Gt 3τ = 0 (7) 2τ G 1 τ + G 2 τ + Gτ = 0 (8) Zero phase at the echo for moving spins leads to the second condition on the first moment M 1 of G (t) leading to ˆ te M 1 (te) = v x G (t) dt = 0 (9) t 2 τ G G t τ τ + G t2 2 3τ 2τ = 0 (10) 1 2 G 1τ G ( 2 4τ 2 τ 2) G ( 9τ 2 4τ 2) = 0 (11) Solving this two-equation-two-unknown system guides to the appropriate gradient waveform { G 1 + G 2 + G = 0 (12) G 1 + 3G 2 + 5G = 0 which has for solution { G 1 = G G 2 = 2G Figure 2 presents the phase behavior as a function of time of both stationary and constant velocity spins in a velocity compensated readout gradient waveform. (13) 5
7 Figure 2: Phase behavior of static and mobile spins during flow compensated readout gradient application As predicted, phase comes to zero at the echo time for both stationary and moving spins, independently of the speed of flow. Sometimes, the second gradient amplitude ( 2G) cannot be performed, in which case the time application must be increased. Usually, velocity compensated sequences are designed to minimize the echo time to avoid T2 signal loss or flow acceleration effects. This demonstration addressed velocity compensation in the read direction, but the same rational reasoning can be used to obtain the gradient waveform for compensation in the slice selection direction. In this case, the zeroth and first moment of the slice selection gradient must be zero at the end of the slice select gradient, with t = 0 fixed at the center of the radio frequency pulse. Further development leads to a (G, 2G, G) gradient scheme. 4.3 Phase contrast angiography PCA uses the effect of motion on spin phase in order to encode velocities. Square bipolar pulses are added prior to the data collection to the flow compensated gradient scheme found in Figure 2. With an analysis similar to the one found in Section 4.2, it is seen that spin moving in the gradient axe develops a phase proportional to their velocity. At the end of the bipolar pulses, the phase φ acquired by moving spins becomes where τ stands for the duration of each two bipolar gradient lobes. φ = γgv z τ 2 (14) Note that there are other sources of phase alteration due to field inhomogeneities. To suppress these background phases, two images are collected in a standard PCA Kernel with an inversion of the polarity of the bipolar gradient lobes. The corresponding phase images are then subtracted and velocity can be extracted using φ = 2γGv z τ 2 (15) 6
8 In order to get proper velocity measurements, the range of phase must stand between π and π to avoid aliasing effects. This condition must be transposed to Equation 14 γgv z τ 2 < π In a PCA sequence, velocity encoding bipolar lobes are determined from a maximum velocity venc parameter set by the user venc = π /γgτ 2 For accurate velocity measurements, experimentator must set venc parameter nearest to the maximum velocity expected. Also it is essential to ensure that flow arises strictly perpendicular to the encoding plane to avoid partial volume effects or to underevaluate velocity. Velocity maps obtained with PCA sequences have found multiple application such as evaluation of blood volume in human vessel with 2D acquisition. Velocity maps can be used straightforwardly to determine the actual flow through a vessel ˆ F = v z da Rigourous blood flow quantification requires cardiac triggering and can lead to a good approximation of blood volume. A series of transverse images are generally collected during cardiac cycle across great vessels like the aorta or the vena cava. 7
9 5 Implementation This section describes the approach used for the implementation of the PCA sequence. 5.1 FLASH.cpp sample All implementations were achieved using the FLASH sample (Figure 3) provided with IDEA installation. This sample proposes a typical gradient echo sequence with multiple acquisition options (preparation pulses, reconstruction options, physiological triggering, etc.). GSlice GSliceSpoil GS GSliReph GPhase GPhaseRew GP GRead GReadSpoil GR GReadDeph Figure 3: Flash sample kernel diagram which consists in a standard 2D gradient echo with spoiler gradients added in the slice and read encoding directions The sample s sequence library provided four mandatory points implemented as C++ member functions which were slightly modified to perform PCA : 1. fseqinit : Set the hard limits of the sequence (maximum gradient amplitude, minimum rising time, etc.). 2. fseqprep : Verify if protocol is consistent, perform calculations before measurement. 3. fseqrun : Handle sequence timing, provide Image Calculation Environment (ICE) programs with information. 4. fseqcheck : Support the gradient overflow check. Next section describes succinctly the changes made in the sample source code for the PCA sequence. 8
10 5.2 Implementation of PCA sequence Before starting the implementation, the desired sequence kernel (Figure 4) was drawn to get a general view of the problem. GSlice GSFQOuter GSlice -GSFQInner GS GSFQInner -GSFQOuter GPhase GPhaseRew GPhase GPhaseRew GP GRFCOuter GRFCOuter GRead GReadSpoil GRead GReadSpoil GR GRFCInner GRFCInner Figure 4: PCA kernel diagram Comparing Figure 3 and Figure 4 lead to the following modifications in the FLASH sample source code : 1. In fseqinit : Allow two images calculation (There are two different acquisitions for each phase step). 2. In fseqprep : Create and prepare flow quantification/compensation gradients in the slice selection and read encoding directions. 3. In fseqprep : Adapt Echo Time (TE) calculation to incorporate new gradient objects. 4. In fseqrun : Adapt sequence timing. To allow two images calculation, sequence contrast parameter was set to two as in a Double Echo sequence. This procedure ensured Image Calcul Environment (ICE) programs to expect two different set of measurement data headers (Mdh) and reconstruct two images. Flow quantification was performed in the slice selection direction and read encoding direction was flow compensated to avoid phase alteration caused by non-perpendicular field of view (FOV) positioning. As discussed in the theoretical section, flow compensation requires a (G, 2G, G) gradient scheme. Since gradient calculation is generally restricted by amplitude and minimum timing, it can be arduous 9
11 to get the perfect compromise between maximum amplitude and minimum TE. Implementation of such optimization from scratch could have been delicate, but fortunately flow quantification and compensation functions were supported in the sequence utility library. The IDEA library contained two functions for both flow quantification and compensation : 1. fsugslflowcalc : Function returning the appropriate bipolar gradient pair characteristics to fill all the time available before TE (Figure 5). 2. fsucalcshortestflowcomp : Calculates a gradient pair in the shortest possible time. These functions require general information about the gradient to be compensated and return the pair of compensating gradients. Functions could be called with multiple variable input or with pointers to the gradient pulses. Function fsugslflowcalc was used in the implementation. One of the key input parameters of fsugslflowcalc was the maximum velocity encoded, venc. For full flow compensation, venc had to be set to 0. For flow quantification gradients calculation, this parameter had to be modified in the protocol by the user. There were several ways to enable venc parameter adjustement in the user interface. For learning purposes, the sequence special card, which is specially designed for sequence programmers, was used. A special venc parameter could be find in this sub-card of the sequence card for the expected maximum velocity in the slice direction. This venc could admit values from 0 to 5000 mm /s with increment of 100 mm /s. F0S R0S R1 F1 ΔG GS G0 G1 G2 ΔG R2 F2 Time available Flow quantification bipolar Gradient ΔG ΔG Figure 5: Approach behind flow quantificaton gradient calculation with fsugslflowcalc Flow sensitive gradient pair G are added to the flow compensated gradients G 1 and G 2. A greater venc parameter leads to a smaller G. The user has to be aware that there is a minimum 10
12 venc available in the scanner restricted by the maximum gradient amplitude. If flow quantification gradient exceeds these hard limits all gradients are automatically set to 0 for safety purposes. All flow quantification and compensation gradient objects were incorporated in the timing table to obtain the desired sequence kernel. The new gradient durations were also included in the time available calculation to calculate the proper TE. All the lines added in the original sample code can be found in the PCAngio.cpp source code and are commented with \\DOMINIQUE to facilitate research. 6 Preliminary results Two scan sessions were organized on the MRI platform during the project to confirm that the implementation could produce results. 6.1 First scan session : Untriggered PCA acquisition This session was performed the 12 th March 2009 on a 22 year-old male subject. The test purpose was to ensure that the PCA sequence could generate two different phase images with proper bipolar gradient inversion. Plane of study was positioned around the carotid junction level using a 2D multislice FLASH acquisition. Magnitude and phase reconstruction images were obtained with this set of parameters : T E = 10ms, T R = 35ms, F lipangle = 35, SliceT hickness = 5mm, M atrixsize = Maximum velocity was set to venc = 1500 mm/s to avoid aliasing. Table 2 presents the results. Table 2: PCA images with V max set to 1500 mm /s (Left : Modulus ; Center : Positive bipolar acquisition ; Right : negative bipolar acquisition) The sequence implemented could be used to reconstruct two different phase images with opposite flow encoding bipolar gradients. Both phase images revealed that there was an obvious phase pattern caused by the magnetic field inhomogeneities. This field pattern was not affected by the gradient inversion. The phase images were subtracted after the scan session using Neurolens to spot phase modifications in the blood flowing regions. In Table 3, the image to the left shows the basic subtraction, the middle one shows the subtraction with adapted intensity level to mask noise, and the right one shows the superposition of the second image and the modulus image. 11
13 Table 3: PCA phase images subtraction with V max set to 1500 mm /s With appropriate intensity adjustments as shown in the third image, at least two blood vessels with opposite flowing directions were vizualised. On the other hand, there was a lot of flow artifacts caused by the blood flow pulsativity and possibly the plane orientation which could have not been strictly perpendicular to the vessels. These artifacts caused periodic repetition of vessel regions along phase encoding direction corrupting dramatically the velocity maps. 6.2 Second scan session : PCA acquisition with ECG triggering To remove flow artifacts, an other scan test was performed the 29 th May 2009 on a 21 year-old female subject. The test purpose was to certify the PCA sequence implemented could handle ECG triggering. Plane of study was positioned around the carotid junction level using a 2D multislice FLASH acquisition. Magnitude and phase reconstruction images were obtained with the same set of parameters defined in Section 6.1. Maximum velocity was fixed to venc = 1500 mm/s and venc = 1000 mm/s to observe how phase varied with bipolar gradient amplitude. Scan session was performed with prospective ECG triggering. Acquisition window was reduced in order to acquire one k-space line per R wave. Table 4 shows the results for venc = 1500 mm/s. Table 4: PCA images with venc set to 1500 mm /s (Left : Modulus ; Center : Basic subtraction ; Right : Subtraction with adapted intensity) With the triggered method, phase variations were much more condensed in the expected vessel regions. Also, ghosting along the phase encoding direction due to flow pulsativity was reduced. To verify the dependence between venc and phase accumulation, another set of images was obtained 12
14 with venc = 1000 mm/s. Table 5 displays the maps obtained with the same intensity adjustments as those used in the previous Table. Table 5: PCA images with venc set to 1000 mm /s As predicted by the theory, subtracted phase images revealed a greater phase accumulation in the vessel regions with venc = 1000 mm/s. The relation among flow encoding gradient amplitude, phase accumulation and velocity was observed practically. 7 Conclusion To conclude, this challenging project was an excellent opportunity to get more familiar with sequence programming and to see how physics principles could be integrated in a source code. The final version of the PCA sequence is not ready to perform accurate velocity measurements, but it can still be used to observe practically how moving spin exposed to bipolar gradients can accumulate a phase term. Optimization of the sequence could include implementation of flow compensation gradient lobes in the phase encoding direction to avoid mismeasurements due to inacurate plane of study positioning. Calculation of such gradients was not supported in the IDEA library and would have to be implemented from scratch. Also, ICE programs could be modified to perform the basic subtraction of phase images and the velocity map normalization. 13
15 8 References BERNSTEIN M. A. et al. Handbook of MRI pulse sequences. Elsevier, p. HAACKE E. M. et al. Magnetic Resonance Imaging. New-York : Wiley-Liss, p. HOA D. et al. L IRM pas à pas. Paris : Campus Medica, p. MATHIEU D. Rapport de stage. Grenoble : Institut des Neurosciences, 42 p. POPE J. M. et YAO S. Quantitative NMR Imaging of Flow. Conceptis in Magnetic Resonance, 1993, vol. 5, issue 4, pp IDEA Manual. Integrated Development Environment for Applications p. 14
Correction Gradients. Nov7, Reference: Handbook of pulse sequence
Correction Gradients Nov7, 2005 Reference: Handbook of pulse sequence Correction Gradients 1. Concomitant-Field Correction Gradients 2. Crusher Gradients 3. Eddy-Current Compensation 4. Spoiler Gradients
More informationChapter 24 MRA and Flow quantification. Yongquan Ye, Ph.D. Assist. Prof. Radiology, SOM Wayne State University
Chapter 24 MRA and Flow quantification Yongquan Ye, Ph.D. Assist. Prof. Radiology, SOM Wayne State University Previous classes Flow and flow compensation (Chap. 23) Steady state signal (Cha. 18) Today
More informationM R I Physics Course. Jerry Allison Ph.D., Chris Wright B.S., Tom Lavin B.S., Nathan Yanasak Ph.D. Department of Radiology Medical College of Georgia
M R I Physics Course Jerry Allison Ph.D., Chris Wright B.S., Tom Lavin B.S., Nathan Yanasak Ph.D. Department of Radiology Medical College of Georgia M R I Physics Course Spin Echo Imaging Hahn Spin Echo
More informationMR Advance Techniques. Flow Phenomena. Class I
MR Advance Techniques Flow Phenomena Class I Flow Phenomena In this class we will explore different phenomenona produced from nuclei that move during the acquisition of data. Flowing nuclei exhibit different
More informationNMR/MRI examination (8N080 / 3F240)
NMR/MRI examination (8N080 / 3F240) Remarks: 1. This test consists of 3 problems with at total of 26 sub-questions. 2. Questions are in English. You are allowed to answer them in English or Dutch. 3. Please
More informationA Study of Flow Effects on the Gradient Echo Sequence
-MR Flow Imaging- A Study of Flow Effects on the Gradient Echo Sequence Cylinder filled with doped water α pulse α pulse Flowing water Plastic pipes Slice Phase Read a TE b Signal sampling TR Thesis for
More informationIntroduction to Biomedical Imaging
Alejandro Frangi, PhD Computational Imaging Lab Department of Information & Communication Technology Pompeu Fabra University www.cilab.upf.edu MRI advantages Superior soft-tissue contrast Depends on among
More informationEL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam
EL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam (closed book, 1 sheets of notes double sided allowed, no calculator or other electronic devices allowed) 1. Ultrasound Physics (15 pt) A) (9
More informationContrast Mechanisms in MRI. Michael Jay Schillaci
Contrast Mechanisms in MRI Michael Jay Schillaci Overview Image Acquisition Basic Pulse Sequences Unwrapping K-Space Image Optimization Contrast Mechanisms Static and Motion Contrasts T1 & T2 Weighting,
More informationOutlines: (June 11, 1996) Instructor:
Magnetic Resonance Imaging (June 11, 1996) Instructor: Tai-huang Huang Institute of Biomedical Sciences Academia Sinica Tel. (02) 2652-3036; Fax. (02) 2788-7641 E. mail: bmthh@ibms.sinica.edu.tw Reference:
More informationK-space. Spin-Warp Pulse Sequence. At each point in time, the received signal is the Fourier transform of the object s(t) = M( k x
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2015 MRI Lecture 4 k (t) = γ 2π k y (t) = γ 2π K-space At each point in time, the received signal is the Fourier transform of the object
More informationNMR and MRI : an introduction
Intensive Programme 2011 Design, Synthesis and Validation of Imaging Probes NMR and MRI : an introduction Walter Dastrù Università di Torino walter.dastru@unito.it \ Introduction Magnetic Resonance Imaging
More informationNavigator Echoes. BioE 594 Advanced Topics in MRI Mauli. M. Modi. BioE /18/ What are Navigator Echoes?
Navigator Echoes BioE 594 Advanced Topics in MRI Mauli. M. Modi. 1 What are Navigator Echoes? In order to correct the motional artifacts in Diffusion weighted MR images, a modified pulse sequence is proposed
More informationApodization. Gibbs Artifact. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2013 MRI Lecture 5. rect(k x )
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2013 MRI Lecture 5 GE Medical Systems 2003 Gibbs Artifact Apodization rect(k ) Hanning Window h(k )=1/2(1+cos(2πk ) 256256 image 256128
More informationRelaxation times in nuclear magnetic resonance
Relaxation times in TEP Related topics Nuclear spins, atomic nuclei with a magnetic moment, precession movement of the nuclear spins, Landau-Lifshitz equation, Bloch equation, magnetisation, resonance
More informationIntroduction to MRI Acquisition
Introduction to MRI Acquisition James Meakin FMRIB Physics Group FSL Course, Bristol, September 2012 1 What are we trying to achieve? 2 What are we trying to achieve? Informed decision making: Protocols
More informationPhysics of MR Image Acquisition
Physics of MR Image Acquisition HST-583, Fall 2002 Review: -MRI: Overview - MRI: Spatial Encoding MRI Contrast: Basic sequences - Gradient Echo - Spin Echo - Inversion Recovery : Functional Magnetic Resonance
More informationIntroduction to Magnetic Resonance Imaging (MRI) Pietro Gori
Introduction to Magnetic Resonance Imaging (MRI) Pietro Gori Enseignant-chercheur Equipe IMAGES - Télécom ParisTech pietro.gori@telecom-paristech.fr September 20, 2017 P. Gori BIOMED 20/09/2017 1 / 76
More informationBNG/ECE 487 FINAL (W16)
BNG/ECE 487 FINAL (W16) NAME: 4 Problems for 100 pts This exam is closed-everything (no notes, books, etc.). Calculators are permitted. Possibly useful formulas and tables are provided on this page. Fourier
More informationROCHESTER INSTITUTE OF TECHNOLOGY COURSE OUTLINE FORM COLLEGE OF SCIENCE. Chester F. Carlson Center for Imaging Science
ROCHESTER INSTITUTE OF TECHNOLOGY COURSE OUTLINE FORM COLLEGE OF SCIENCE Chester F. Carlson Center for Imaging Science NEW COURSE: COS-IMGS-730 Magnetic Resonance Imaging 1.0 Course Designations and Approvals
More informationPulse Sequences: RARE and Simulations
Pulse Sequences: RARE and Simulations M229 Advanced Topics in MRI Holden H. Wu, Ph.D. 2018.04.19 Department of Radiological Sciences David Geffen School of Medicine at UCLA Class Business Final project
More informationEE225E/BIOE265 Spring 2016 Principles of MRI. Assignment 4. Due Friday Feb 19st, 2016, Self Grading Due Monday Feb 22nd, 2016
EE225E/BIOE265 Spring 2016 Principles of MRI Miki Lustig Assignment 4 Due Friday Feb 19st, 2016, Self Grading Due Monday Feb 22nd, 2016 1. Finish reading Nishimura Ch.4 and Ch. 5. 2. The following pulse
More informationDiffusion Tensor Imaging (DTI): An overview of key concepts
Diffusion Tensor Imaging (DTI): An overview of key concepts (Supplemental material for presentation) Prepared by: Nadia Barakat BMB 601 Chris Conklin Thursday, April 8 th 2010 Diffusion Concept [1,2]:
More informationBasic MRI physics and Functional MRI
Basic MRI physics and Functional MRI Gregory R. Lee, Ph.D Assistant Professor, Department of Radiology June 24, 2013 Pediatric Neuroimaging Research Consortium Objectives Neuroimaging Overview MR Physics
More informationLab 2: Magnetic Resonance Imaging
EE225E/BIOE265 Spring 2013 Principles of MRI Miki Lustig Developed by: Galen Reed and Miki Lustig Lab 2: Magnetic Resonance Imaging Introduction In this lab, we will get some hands-on experience with an
More informationMRI in Practice. Catherine Westbrook MSc, DCRR, CTC Senior Lecturer Anglia Polytechnic University Cambridge UK. John Talbot MSc, DCRR
MRI in Practice Third edition Catherine Westbrook MSc, DCRR, CTC Senior Lecturer Anglia Polytechnic University Cambridge UK and Carolyn Kaut RothRT(R) (MR) (CT) (M) (CV) Fellow SMRT (Section for Magnetic
More informationField trip: Tuesday, Feb 5th
Pulse Sequences Field trip: Tuesday, Feb 5th Hardware tour of VUIIIS Philips 3T Meet here at regular class time (11.15) Complete MRI screening form! Chuck Nockowski Philips Service Engineer Reminder: Project/Presentation
More informationVelocity Images. Phase Contrast Technique. G. Reiter 1,2, U. Reiter 1, R. Rienmüller 1
Velocity Images - the MR Phase Contrast Technique G. Reiter 1,2, U. Reiter 1, R. Rienmüller 1 SSIP 2004 12 th Summer School in Image Processing, Graz, Austria 1 Interdisciplinary Cardiac Imaging Center,
More informationRochester Institute of Technology Rochester, New York. COLLEGE of Science Department of Chemistry. NEW (or REVISED) COURSE:
Rochester Institute of Technology Rochester, New York COLLEGE of Science Department of Chemistry NEW (or REVISED) COURSE: 1014-730 1.0 Title: Magnetic Resonance Imaging (MRI) Date: July 2006 Credit Hours:
More informationThe NMR Inverse Imaging Problem
The NMR Inverse Imaging Problem Nuclear Magnetic Resonance Protons and Neutrons have intrinsic angular momentum Atoms with an odd number of proton and/or odd number of neutrons have a net magnetic moment=>
More informationIntroduction to MRI. Spin & Magnetic Moments. Relaxation (T1, T2) Spin Echoes. 2DFT Imaging. K-space & Spatial Resolution.
Introduction to MRI Spin & Magnetic Moments Relaxation (T1, T2) Spin Echoes 2DFT Imaging Selective excitation, phase & frequency encoding K-space & Spatial Resolution Contrast (T1, T2) Acknowledgement:
More informationChapter 15:Magnetic Resonance Imaging
Chapter 15:Magnetic Resonance Imaging Slide set of 242 slides based on the chapter authored by Martin O. Leach of the publication (ISBN 978-92-0-131010-1): Diagnostic Radiology Physics: A Handbook for
More informationTissue Parametric Mapping:
Tissue Parametric Mapping: Contrast Mechanisms Using SSFP Sequences Jongho Lee Department of Radiology University of Pennsylvania Tissue Parametric Mapping: Contrast Mechanisms Using bssfp Sequences Jongho
More informationMRI Physics II: Gradients, Imaging. Douglas C. Noll, Ph.D. Dept. of Biomedical Engineering University of Michigan, Ann Arbor
MRI Physics II: Gradients, Imaging Douglas C., Ph.D. Dept. of Biomedical Engineering University of Michigan, Ann Arbor Magnetic Fields in MRI B 0 The main magnetic field. Always on (0.5-7 T) Magnetizes
More informationRad Tech 4912 MRI Registry Review. Outline of the Registry Exam: Certification Fees
Rad Tech 4912 MRI Registry Review Outline of the Registry Exam: Category: # of questions: A. Patient Care 30 B. Imaging Procedures 62 C. Data Acquisition and Processing 65 D. Physical Principles of Image
More informationFundamental MRI Principles Module 2 N. Nuclear Magnetic Resonance. X-ray. MRI Hydrogen Protons. Page 1. Electrons
Fundamental MRI Principles Module 2 N S 1 Nuclear Magnetic Resonance There are three main subatomic particles: protons positively charged neutrons no significant charge electrons negatively charged Protons
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 informationTechnical University of Denmark
Technical University of Denmark Page 1 of 10 pages Written test, 12 December 2012 Course name: Introduction to medical imaging Course no. 31540 Aids allowed: None. Pocket calculator not allowed "Weighting":
More informationAdvanced Topics and Diffusion MRI
Advanced Topics and Diffusion MRI Slides originally by Karla Miller, FMRIB Centre Modified by Mark Chiew (mark.chiew@ndcn.ox.ac.uk) Slides available at: http://users.fmrib.ox.ac.uk/~mchiew/teaching/ MRI
More informationFundamentals of MR Imaging
Fundamentals of MR Imaging Shantanu Sinha. Department of Radiology UCSD School of Medicine, San Diego, CA-92103. E-mail: shsinha@ucsd.edu Background References: R.B.Lufkin, The MRI Manual (2nd Edition).
More informationApplications of Spin Echo and Gradient Echo: Diffusion and Susceptibility Contrast
Applications of Spin Echo and Gradient Echo: Diffusion and Susceptibility Contrast Chunlei Liu, PhD Department of Electrical Engineering & Computer Sciences and Helen Wills Neuroscience Institute University
More informationIntroduction to functional MRI in humans. Michael Hallquist University of Pittsburgh
Introduction to functional MRI in humans Michael Hallquist University of Pittsburgh Goals of human neuroimaging Localization of brain function (mapping) Understanding large-scale functional integration
More informationPart III: Sequences and Contrast
Part III: Sequences and Contrast Contents T1 and T2/T2* Relaxation Contrast of Imaging Sequences T1 weighting T2/T2* weighting Contrast Agents Saturation Inversion Recovery JUST WATER? (i.e., proton density
More informationBasic Pulse Sequences I Saturation & Inversion Recovery UCLA. Radiology
Basic Pulse Sequences I Saturation & Inversion Recovery Lecture #5 Learning Objectives Explain what the most important equations of motion are for describing spin systems for MRI. Understand the assumptions
More informationEE225E/BIOE265 Spring 2013 Principles of MRI. Assignment 9 Solutions. Due April 29th, 2013
EE5E/BIOE65 Spring 013 Principles of MRI Miki Lustig This is the last homework in class. Enjoy it. Assignment 9 Solutions Due April 9th, 013 1) In class when we presented the spin-echo saturation recovery
More informationMRI Physics I: Spins, Excitation, Relaxation
MRI Physics I: Spins, Excitation, Relaxation Douglas C. Noll Biomedical Engineering University of Michigan Michigan Functional MRI Laboratory Outline Introduction to Nuclear Magnetic Resonance Imaging
More informationFundamental MRI Principles Module Two
Fundamental MRI Principles Module Two 1 Nuclear Magnetic Resonance There are three main subatomic particles: protons neutrons electrons positively charged no significant charge negatively charged Protons
More informationChapter 14:Physics of Magnetic Resonance
Chapter 14:Physics of Magnetic Resonance Slide set of 141 slides based on the chapter authored by Hee Kwon Song of the publication (ISBN 978-92-0-131010-1): Diagnostic Radiology Physics: A Handbook for
More informationNMR Spectroscopy Laboratory Experiment Introduction. 2. Theory
1. Introduction 64-311 Laboratory Experiment 11 NMR Spectroscopy Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful and theoretically complex analytical tool. This experiment will introduce to
More informationSuppression of Static Magnetic Field in Diffusion Measurements of Heterogeneous Materials
PIERS ONLINE, VOL. 5, NO. 1, 2009 81 Suppression of Static Magnetic Field in Diffusion Measurements of Heterogeneous Materials Eva Gescheidtova 1 and Karel Bartusek 2 1 Faculty of Electrical Engineering
More informationSequence Overview. Gradient Echo Spin Echo Magnetization Preparation Sampling and Trajectories Parallel Imaging. B.Hargreaves - RAD 229
Sequence Overview Gradient Echo Spin Echo Magnetization Preparation Sampling and Trajectories Parallel Imaging 75 Pulse Sequences and k-space RF k y G z k x G x 3D k-space G y k y k z Acq. k x 76 Gradient
More informationLecture #7 In Vivo Water
Lecture #7 In Vivo Water Topics Hydration layers Tissue relaxation times Magic angle effects Magnetization Transfer Contrast (MTC) CEST Handouts and Reading assignments Mathur-De Vre, R., The NMR studies
More informationBMB 601 MRI. Ari Borthakur, PhD. Assistant Professor, Department of Radiology Associate Director, Center for Magnetic Resonance & Optical Imaging
BMB 601 MRI Ari Borthakur, PhD Assistant Professor, Department of Radiology Associate Director, Center for Magnetic Resonance & Optical Imaging University of Pennsylvania School of Medicine A brief history
More informationSynchronous Multi-Directional Motion Encoding. in Magnetic Resonance Elastography. DAVID A. BURNS B.S., University of Illinois, Urbana-Champaign, 2010
Synchronous Multi-Directional Motion Encoding in Magnetic Resonance Elastography BY DAVID A. BURNS B.S., University of Illinois, Urbana-Champaign, 2010 THESIS Submitted as partial fulfillment of the requirements
More informationDoppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.
1 Doppler echocardiography & Magnetic Resonance Imaging History: - Langevin developed sonar. - 1940s development of pulse-echo. - 1950s development of mode A and B. - 1957 development of continuous wave
More informationMagnetization Preparation Sequences
Magnetization Preparation Sequences Acquisition method may not give desired contrast Prep block adds contrast (and/or encoding) MP-RAGE = Magnetization prepared rapid acquisition with gradient echo (Mugler,
More informationLab 1: Intro to NMR. March 10, 2014
Lab 1: Intro to NMR March 10, 2014 Christine Leon-Swisher (GSI), Miki Lustig (Prof) 1 Preliminaries Never bring anything metal into the room with an MRI (i.e. keys, metallic jewelry, chairs) Do not enter
More informationNuclear Magnetic Resonance Imaging
Nuclear Magnetic Resonance Imaging Jeffrey A. Fessler EECS Department The University of Michigan NSS-MIC: Fundamentals of Medical Imaging Oct. 20, 2003 NMR-0 Background Basic physics 4 magnetic fields
More informationExam 8N080 - Introduction to MRI
Exam 8N080 - Introduction to MRI Friday April 10 2015, 18.00-21.00 h For this exam you may use an ordinary calculator (not a graphical one). In total there are 5 assignments and a total of 50 points can
More informationIntroductory MRI Physics
C HAPR 18 Introductory MRI Physics Aaron Sodickson EXRNAL MAGNETIC FIELD, PROTONS AND EQUILIBRIUM MAGNETIZATION Much of the bulk of the magnetic resonance imaging (MRI) scanner apparatus is dedicated to
More informationNMR Imaging in porous media
NMR Imaging in porous media What does NMR give us. Chemical structure. Molecular structure. Interactions between atoms and molecules. Incoherent dynamics (fluctuation, rotation, diffusion). Coherent flow
More informationVelocity k-space analysis of Flow Effects in Echo-Planar, Spiral and Projection Reconstruction Imaging. Sangeetha Somayajula
Velocity k-space analysis of Flow Effects in Echo-Planar, Spiral and Projection Reconstruction Imaging Sangeetha Somayajula Introduction Flow and motion in the object causes distortion in the MRI signal
More informationNuclear Magnetic Resonance Imaging
Nuclear Magnetic Resonance Imaging Simon Lacoste-Julien Electromagnetic Theory Project 198-562B Department of Physics McGill University April 21 2003 Abstract This paper gives an elementary introduction
More informationBiomedical Imaging Magnetic Resonance Imaging
Biomedical Imaging Magnetic Resonance Imaging Charles A. DiMarzio & Eric Kercher EECE 4649 Northeastern University May 2018 Background and History Measurement of Nuclear Spins Widely used in physics/chemistry
More informationSpin Echo Imaging Sequence
1 MRI In Stereotactic Procedures Edward F. Jackson, Ph.D. The University of Texas M.D. Anderson Cancer Center Houston, Texas 2 RF G slice G phase G freq Signal k-space Spin Echo Imaging Sequence TE 1st
More informationLecture k-space. k-space illustrations. Zeugmatography 3/7/2011. Use of gradients to make an image echo. K-space Intro to k-space sampling
Lecture 21-3-16 K-space Intro to k-space sampling (chap 3) Frequenc encoding and Discrete sampling (chap 2) Point Spread Function K-space properties K-space sampling principles (chap 3) Basic Contrast
More informationEffect of Bulk Tissue Motion on Quantitative Perfusion and Diffusion Magnetic Resonance Imaging *
MAGNETIC RESONANCE IN MEDICINE 19,261-265 (1991) Effect of Bulk Tissue Motion on Quantitative Perfusion and Diffusion Magnetic Resonance Imaging * THOMAS L. CHENEVERT AND JAMES G. PIPE University of Michigan
More informationChapter 1 Pulsed Field Gradient NMR Sequences
Chapter 1 Pulsed Field Gradient NMR Sequences Abstract The mechanism via which an NMR signal is generated and how diffusion can be measured from solving the equation of motion of the nuclear spins is described.
More informationProfessor Stuart Bunt 217
Professor Stuart Bunt 217 Traditional Anatomy Phrenology, the study of bumps on the skull. Measuring brain weights and size (still being done..see the fuss about Einstein s brain). Little link between
More informationPart II: Magnetic Resonance Imaging (MRI)
Part II: Magnetic Resonance Imaging (MRI) Contents Magnetic Field Gradients Selective Excitation Spatially Resolved Reception k-space Gradient Echo Sequence Spin Echo Sequence Magnetic Resonance Imaging
More informationA Lecture on Selective RF-pulses in MRI
A Lecture on Selective RF-pulses in MRI Charles L. Epstein February, 23 Introduction We describe the notion of selective excitation and explain how the Bloch equation is used to design selective RF-pulses.
More informationMRI in Review: Simple Steps to Cutting Edge Part I
MRI in Review: Simple Steps to Cutting Edge Part I DWI is now 2 years old... Mike Moseley Radiology Stanford DWI, b = 1413 T2wt, 28/16 ASN 21 San Francisco + Disclosures: Funding NINDS, NCRR, NCI 45 minutes
More informationThe Basics of Magnetic Resonance Imaging
The Basics of Magnetic Resonance Imaging Nathalie JUST, PhD nathalie.just@epfl.ch CIBM-AIT, EPFL Course 2013-2014-Chemistry 1 Course 2013-2014-Chemistry 2 MRI: Many different contrasts Proton density T1
More informationSpatial encoding in Magnetic Resonance Imaging. Jean-Marie BONNY
Spatial encoding in Magnetic Resonance Imaging Jean-Marie BONNY What s Qu est an image ce qu une? image? «a reproduction of a material object by a camera or a related technique» Multi-dimensional signal
More informationFREQUENCY SELECTIVE EXCITATION
PULSE SEQUENCES FREQUENCY SELECTIVE EXCITATION RF Grad 0 Sir Peter Mansfield A 1D IMAGE Field Strength / Frequency Position FOURIER PROJECTIONS MR Image Raw Data FFT of Raw Data BACK PROJECTION Image Domain
More informationLecture 21. Nuclear magnetic resonance
Lecture 21 Nuclear magnetic resonance A very brief history Stern and Gerlach atomic beam experiments Isidor Rabi molecular beam exp.; nuclear magnetic moments (angular momentum) Felix Bloch & Edward Purcell
More informationThe physics US and MRI. Prof. Peter Bogner
The physics US and MRI Prof. Peter Bogner Sound waves mechanical disturbance, a pressure wave moves along longitudinal wave compression rarefaction zones c = nl, (c: velocity, n: frequency, l: wavelength
More informationHow does this work? How does this method differ from ordinary MRI?
361-Lec41 Tue 18nov14 How does this work? How does this method differ from ordinary MRI? NEW kinds of MRI (magnetic resononance imaging (MRI) Diffusion Magnetic Resonance Imaging Tractographic reconstruction
More informationMRI based quantification of outflow boundary conditions for wall shear stress calculations in stenosed human carotid bifurcations
MRI based quantification of outflow boundary conditions for wall shear stress calculations in stenosed human carotid bifurcations Lenette Simons BME 09.34 Master Thesis Medical Engineering (8Z170) 13 th
More informationHST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS
HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT
More informationOn Signal to Noise Ratio Tradeoffs in fmri
On Signal to Noise Ratio Tradeoffs in fmri G. H. Glover April 11, 1999 This monograph addresses the question of signal to noise ratio (SNR) in fmri scanning, when parameters are changed under conditions
More informationMRI Physics (Phys 352A)
MRI Physics (Phys 352A) Manus J. Donahue: mj.donahue@vanderbilt.edu Department of Radiology, Neurology, Physics, and Psychiatry Office: Vanderbilt University Institute of Imaging Science (VUIIS) AAA-3115
More informationBackground II. Signal-to-Noise Ratio (SNR) Pulse Sequences Sampling and Trajectories Parallel Imaging. B.Hargreaves - RAD 229.
Background II Signal-to-Noise Ratio (SNR) Pulse Sequences Sampling and Trajectories Parallel Imaging 1 SNR: Signal-to-Noise Ratio Signal: Desired voltage in coil Noise: Thermal, electronic Noise Thermal
More informationRECOGNITION ALGORITHM FOR DEVELOPING A BRAIN- COMPUTER INTERFACE USING FUNCTIONAL NEAR INFRARED SPECTROSCOPY
The 212 International Conference on Green Technology and Sustainable Development (GTSD212) RECOGNITION ALGORITHM FOR DEVELOPING A BRAIN- COMPUTER INTERFACE USING FUNCTIONAL NEAR INFRARED SPECTROSCOPY Cuong
More informationPhysical fundamentals of magnetic resonance imaging
Physical fundamentals of magnetic resonance imaging Stepan Sereda University of Bonn 1 / 26 Why? Figure 1 : Full body MRI scan (Source: [4]) 2 / 26 Overview Spin angular momentum Rotating frame and interaction
More informationRF Pulse Design. Multi-dimensional Excitation I. M229 Advanced Topics in MRI Kyung Sung, Ph.D Class Business
RF Pulse Design Multi-dimensional Excitation I M229 Advanced Topics in MRI Kyung Sung, Ph.D. 2018.04.10 Class Business Office hours - Instructors: Fri 10-12pm TAs: Xinran Zhong and Zhaohuan Zhang (time:
More informationMagnetic resonance imaging MRI
Magnetic resonance imaging MRI Introduction What is MRI MRI is an imaging technique used primarily in medical settings that uses a strong magnetic field and radio waves to produce very clear and detailed
More informationThe physics of medical imaging US, CT, MRI. Prof. Peter Bogner
The physics of medical imaging US, CT, MRI Prof. Peter Bogner Clinical radiology curriculum blocks of lectures and clinical practice (7x2) Physics of medical imaging Neuroradiology Head and neck I. Head
More informationThe Theory of Nuclear Magnetic Resonance Behind Magnetic Resonance Imaging. Catherine Wasko Physics 304 Physics of the Human Body May 3, 2005
The Theory of Nuclear Magnetic Resonance Behind Magnetic Resonance Imaging Catherine Wasko Physics 304 Physics of the Human Body May 3, 2005 Magnetic resonance imaging (MRI) is a tool utilized in the medical
More informationThe ASL signal. Parenchy mal signal. Venous signal. Arterial signal. Input Function (Label) Dispersion: (t e -kt ) Relaxation: (e -t/t1a )
Lecture Goals Other non-bold techniques (T2 weighted, Mn contrast agents, SSFP, Dynamic Diffusion, ASL) Understand Basic Principles in Spin labeling : spin inversion, flow vs. perfusion ASL variations
More informationBME I5000: Biomedical Imaging
BME I5000: Biomedical Imaging Lecture 9 Magnetic Resonance Imaging (imaging) Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/ 1 Schedule 1. Introduction, Spatial Resolution,
More informationREFAAT E. GABR, PHD Fannin Street, MSE R102D, Houston, Texas 77030
NAME: Refaat Elsayed Gabr REFAAT E. GABR, PHD 3-Jul-13 5 pages PRESENT TITLE: ADDRESS: BIRTHDATE: CITIZENSHIP: Assistant Professor of Radiology Department of Diagnostic and Interventional Imaging University
More information} B 1 } Coil } Gradients } FFT
Introduction to MRI Daniel B. Ennis, Ph.D. Requirements for MRI UCLA DCVI Requirements for MRI Dipoles to Images MR Active uclei e.g. 1 H in H20 Cryogen Liquid He and 2 Magnetic Field (B0) Polarizer ystem
More informationRAD229: Final Exam 2014/ SOLUTIONS You will have 3 hours to complete this Exam
RAD229: Final Exam 2014/2015 - SOLUTIONS You will have 3 hours to complete this Exam Solutions are given in Blue. In some cases, different interpretations may have led to different, but reasonable answers,
More informationFunctional Magnetic Resonance Imaging (FMRI) is an imaging technique for
Chapter 2 Principles of FMRI Functional Magnetic Resonance Imaging (FMRI) is an imaging technique for examining brain function. Since its first appearance in 1991 (Belliveau et al.[8]) the use of FMRI
More informationMagnetic Resonance Imaging. Pål Erik Goa Associate Professor in Medical Imaging Dept. of Physics
Magnetic Resonance Imaging Pål Erik Goa Associate Professor in Medical Imaging Dept. of Physics pal.e.goa@ntnu.no 1 Why MRI? X-ray/CT: Great for bone structures and high spatial resolution Not so great
More informationCambridge University Press MRI from A to Z: A Definitive Guide for Medical Professionals Gary Liney Excerpt More information
Main glossary Aa AB systems Referring to molecules exhibiting multiply split MRS peaks due to spin-spin interactions. In an AB system, the chemical shift between the spins is of similar magnitude to the
More informationMagnetic Resonance Imaging
Magnetic Resonance Imaging History Nuclear magnetic resonance was first described by Isidor Rabi in 1938 - Columbia University, New York City, (Nobel Prize Nobel Prize in Physics 1944) 1946 - Edward Mills
More informationSpatial encoding in Magnetic Resonance Imaging. Jean-Marie BONNY
Spatial encoding in Magnetic Resonance Imaging Jean-Marie BONNY What s Qu est an image ce qu une? image? «a reproduction of a material object by a camera or a related technique» Multi-dimensional signal
More informationCHEM / BCMB 4190/6190/8189. Introductory NMR. Lecture 10
CHEM / BCMB 490/690/889 Introductory NMR Lecture 0 - - CHEM 490/690 Spin-Echo The spin-echo pulse sequence: 90 - τ - 80 - τ(echo) Spins echoes are widely used as part of larger pulse sequence to refocus
More information