1 / 53

MRI Image Formation

MRI Image Formation. Karla Miller FMRIB Physics Group. Image Formation. Gradients and spatial encoding Sampling k -space Trajectories and acquisition strategies Fast imaging Acquiring multiple slices Image reconstruction and artifacts. field strength. field offset.

idallas
Télécharger la présentation

MRI Image Formation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. MRI Image Formation Karla Miller FMRIB Physics Group

  2. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  3. field strength field offset MR imaging is based on precession z y x [courtesy William Overall] Spins precess at the Larmor rate: • = g (B0 + DB)

  4. High field B0 Low field Magnetic Gradients Gradient: Additional magnetic field which varies over space • Gradient adds to B0, so field depends on position • Precessional frequency varies with position! • “Pulse sequence” modulates size of gradient

  5. High field B0 Low field Magnetic Gradients • Spins at each position sing at different frequency • RF coil hears all of the spins at once • Differentiate material at a given position by selectively listening to that frequency Fast precession Slow precession

  6. Simple “imaging” experiment (1D) increasing field

  7. Simple “imaging” experiment (1D) Signal Fourier transform “Image” position Fourier Transform: determines amount of material at a given location by selectively “listening” to the corresponding frequency

  8. y ky 1D Signal 1D “Image” x kx 2DFT 1DFT 2D Signal 2D Image 2D Imaging via 2D Fourier Transform

  9. Analogy: Weather Mapping

  10. y ky x kx 2D Fourier Transform 2DFT Measured signal (frequency-, or k-space) Reconstructed image FT can be applied in any number of dimensions MRI: signal acquired in 2D frequency space (k-space) (Usually) reconstruct image with 2DFT

  11. Gradients and image acquisition • Magnetic field gradients encode spatial position in precession frequency • Signal is acquired in the frequency domain (k-space) • To get an image, acquire spatial frequencies along both x and y • Image is recovered from k-space data using a Fourier transform

  12. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  13. x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x FT Sampling k-space Perfect reconstruction of an object would require measurement of all locations in k-space (infinite!) Data is acquired point-by-point in k-space (sampling)

  14. kx ky kx 2 kxmax Sampling k-space • What is the highest frequency we need to sample in k-space (kmax)? • How close should the samples be in k-space (Dk)?

  15. Frequency spectrum What is the maximum frequency we need to measure? Or, what is the maximum k-space value we must sample (kmax)? FT -kmax kmax

  16. Frequency spectrum

  17. Frequency spectrum

  18. Frequency spectrum

  19. Frequency spectrum

  20. Frequency spectrum

  21. Frequency spectrum Higher frequencies make the reconstruction look more like the original object! Large kmax increases resolution (allows us to distinguish smaller features)

  22. 2D Extension increasing kmax kymax kxmax kxmax kymax 2 kxmax kmaxdetermines image resolution Large kmaxmeans high resolution !

  23. kx ky kx 2 kxmax Sampling k-space • What is the highest frequency we need to sample in k-space (kmax)? • How close should the samples be in k-space (Dk)?

  24. Nyquist Sampling Theorem A given frequency must be sampled at least twice per cycle in order to reproduce it accurately 1 samp/cyc 2 samp/cyc Upper waveform is resolved! Cannot distinguish between waveforms

  25. increasing field Nyquist Sampling Theorem Insufficient sampling forces us to interpret that both samples are at the same location: aliasing

  26. Aliasing (ghosting): inability to differentiate between 2 frequencies makes them appear to be at same location x x Aliased image Applied FOV max ive frequency max ive frequency

  27. kx ky FOV = 1/kx kx x= 1/(2*kxmax) 2 kxmax k-space relations:FOV and Resolution

  28. kx ky xmax = 1/kx kx 2 kxmax = 1/x 2 kxmax k-space relations:FOV and Resolution k-space and image-space are inversely related: resolution in one domain determines extent in other

  29. k-space Image Full-FOV, high-res Full sampling 2DFT Full-FOV, low-res: blurred Reduce kmax Low-FOV, high-res: may be aliased Increase k

  30. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  31. Visualizing k-space trajectories kx(t) =  Gx(t) dt ky(t) =  Gy(t) dt k-space location is proportional to accumulated area under gradient waveforms Gradients move us along a trajectory through k-space !

  32. Raster-scan (2DFT) Acquisition Acquire k-space line-by-line (usually called “2DFT”) Gx causes frequency shift along x: “frequency encode” axis Gy causes phase shift along y: “phase ecode” axis

  33. Echo-planar Imaging (EPI) Acquisition Single-shot (snap-shot): acquire all data at once

  34. Many possible trajectories through k-space…

  35. SNR   Ttotal =  N  Tread Trajectory considerations • Longer readout = more image artifacts • Single-shot (EPI & spiral) warping or blurring • PR & 2DFT have very short readouts and few artifacts • Cartesian (2DFT, EPI) vs radial (PR, spiral) • 2DFT & EPI = ghosting & warping artifacts • PR & spiral = blurring artifacts • SNR for N shots with time per shot Tread :

  36. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  37. Partial k-space If object is entirely real, quadrants of k-space contain redundant information 2 1 c+id a+ib aib cid ky 3 4 kx

  38. Partial k-space Idea: just acquire half of k-space and “fill in” missing data Symmetry isn’t perfect, so must get slightly more than half 1 c+id a+ib measured data aib cid missing data ky kx

  39. ky ky kx kx Multiple approaches Acquire half of each frequency encode Reduced phase encode steps

  40. Parallel imaging(SENSE, SMASH, GRAPPA, iPAT, etc) Surface coils Object in 8-channel array Single coil sensitivity Multi-channel coils: Array of RF receive coils Each coil is sensitive to a subset of the object

  41. Parallel imaging(SENSE, SMASH, GRAPPA, iPAT, etc) Surface coils Object in 8-channel array Single coil sensitivity Coil sensitivity to encode additional information Can “leave out” large parts of k-space (more than 1/2!) Similar uses to partial k-space (faster imaging, reduced distortion, etc), but can go farther

  42. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  43. Slice Selection RF frequency 0 gradient Gz excited slice

  44. t1 t2 t3 t4 t5 t6 2D Multi-slice Imaging excited slice All slices excited and acquired sequentially (separately) Most scans acquired this way (including FMRI, DTI)

  45. excited volume “True” 3D imaging excited volume Repeatedly excite all slices simultaneously, k-space acquisition extended from 2D to 3D Higher SNR than multi-slice, but may take longer Typically used in structural scans

  46. Image Formation • Gradients and spatial encoding • Sampling k-space • Trajectories and acquisition strategies • Fast imaging • Acquiring multiple slices • Image reconstruction and artifacts

  47. Motion Artifacts Motion causes inconsistencies between readouts in multi-shot data (structurals) Usually looks like replication of object edges along phase encode direction PE

  48. Gibbs Ringing (Truncation) Abruptly truncating signal in k-space introduces “ringing” to the image

  49. EPI distortion (warping) field offset image distortion Field map EPI image (uncorrected) Magnetization precesses at a different rate than expected Reconstruction places the signal at the wrong location

  50. EPI unwarping (FUGUE) field map uncorrected corrected Field map tells us where there are problems Estimate distortion from field map and remove it

More Related