Statistical Analysis of Image Reconstructed Fully-Sampled and Sub-Sampled fmri Data
|
|
- Claude Moody
- 5 years ago
- Views:
Transcription
1 Statistical Analysis of Image Reconstructed Fully-Sampled and Sub-Sampled fmri Data Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette University February 1, 2016 NIHR21 NS
2 3 Outline Introduction Varied methods to produce fmri images with varied properties. Reconstruction Voxels are not directly measured (k-space). Reconstructed! Processing Images are processed for enhancement & artifact reduction. Implications Effects of image reconstruction & processing? Mean, Var, Corr? Discussion How was our data was produced and what was done to it? 2
3 Introduction In fmri and fcmri, there has been an amazing amount of advanced analysis and interpretations presented, but little attention has been paid to what the data truly are. Fox et al, PNAS,
4 Introduction In general, reconstructed GRE EPI images look like below. How do we get from the below to the previous activation? And the below isn t even our original measurements mm FOV 2.5 mm 2 In-Plane R Are we ahead of the data with our analyses and interpretations? I 4
5 Reconstruction In fmri and MRI, the measurements taken by the machine are an array of complex-valued spatial frequencies. This array of complex-valued spatial frequencies need to be reconstructed into an image for us to see, analyze, and interpret. The array of complex-valued spatial frequencies are reconstructed into an image via the inverse Fourier transform. So lets briefly remind ourselves what the FT and IFT are. 5
6 Coil wraps around so uniform sensitivity. Reconstruction Single Coil Acquisition Coil measures k-space. Coil ky kx Δky Δkx Unaccelerated Acquisition (A=1). Kumar et al: JMR, 18(1):69-83,
7 Reconstruction (FOV=240 mm) (n x =n y =96, Δx=Δ y=2.5 mm) We inverse Fourier transform spatial freqs to generate image i spatial frequencies 7
8 Reconstruction (FOV=240 mm) (n x =n y =96, Δx=Δ y=2.5 mm) We inverse Fourier transform spatial freqs to generate image. + i +i +i IFT matrix spatial frequencies IFT matrix 8
9 Reconstruction (FOV=240 mm) (n x =n y =96, Δx=Δ y=2.5 mm) We inverse Fourier transform spatial freqs to generate image. + i +i +i = +i IFT matrix spatial frequencies IFT matrix image 9
10 Reconstruction (FOV=240 mm) (n x =n y =96, Δx=Δ y=2.5 mm) We inverse Fourier transform spatial freqs to generate image. + i +i +i exp(i ) IFT matrix spatial frequencies IFT matrix Phase discarded. Ask for the other ½ of YOUR data. image 10
11 Reconstruction The machine Fourier encodes the image. Measure spatial freq i +i +i = +i 0 FT matrix image FT matrix spatial frequencies 11
12 R Reconstruction We can stack freq. rows of reals over rows of imaginaries, I 12
13 R Reconstruction We can stack freq. rows of reals over rows of imaginaries, make one IFT reconstruction matrix from the two, I f R f I 13
14 R Reconstruction We can stack freq. rows of reals over rows of imaginaries, make one IFT reconstruction matrix from the two, to get the rows of reals over rows of imaginaries. I y R f R y I f I 14
15 Processing Many processing operations are performed by the scanner, by physicists, and by engineers before statistical analysis. + i k-space Processing Nyquist Ghost Correction Static B0 Field Correction Zero Fill Interpolation Non-Cartesian Interpolation Ramp Sampling Interpolation Homodyne Interpolation Apodization And many more Image Reconstruction 2D inverse Fourier transform In-Plane SENSE/GRAPPA Through-Plane SENSE Image Processing Image Smoothing Global Normalization Motion Correction And many more Time Series Processing Filtering Smoothing Dynamic B0 Correction Slice Timing And many more Show ones in blue. 15
16 Reconstruction We can stack freq. rows of reals over rows of imaginaries, make one IFT reconstruction matrix from the two, to get the rows of reals over rows of imaginaries. y R f R O I O k y I f I 16
17 Processing f k-space 17
18 Processing y O I OR O k f 1 1 S Ω A Z = Processed Image Smoothing IFT Reconstruction Apodization Zero-Filling k-space y=of 18
19 Processing We measure an array of complex-valued numbers, perform complex-valued image reconstruction to this array, to generate complex-valued images in real and imaginary, along the way, there is complex-valued image processing. What are the implications of what was done to the data? 19
20 Implications In statistics, we know the rule that says: If a vector f has a mean δ, and a covariance Γ, Then y=of has a mean μ=oδ, and a covariance Σ=OΓO T. Then Σ can converted into a correlation matrix R=D -1/2 ΣD -1/2. 1/2 Where D 1/ diag( ). Assume k-space measurements independent so Γ=I. 20
21 Implications Operators, O. S Ω A Z
22 P M BIRS Neuroimaging Data Analysis Implications Mean, µ=of. O=Ω O=ΩZ O=ΩA O=ΩAZ O=SΩAZ 22
23 Implications Correlation matrix and correlation image x5 image cor(y) = 25x25 correlation matrix x5 correlation image 23
24 Implications Correlation, R=D -1/2 ΣD -1/2. R R R RR IR R R RI II +1-1 TH=
25 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition Coil 1 k y k x Coil 2 Δky Δkx Each coil measures k-space. Coil 3 N C =4, A=1 Hyde et al.: JMR, 70: , Pruessmann et al.: Proc. ISRM, 579,
26 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition S 1 S 4 Coil 1 v S 2 Coil 2 Coil 3 Each coil measures k-space. S 3 N C =4, A=1 26
27 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition a 4 =S 4 v S 1 a 1 =S 1 v Coil 1 S 4 v S 2 Coil 2 Coil 3 Each coil measures k-space. a 3 =S 3 v S 3 a 2 =S 2 v N C =4, A=1 27
28 Reconstruction SENSE Measured Coil Images Estimated Sensitivities Combined Image ^ v & N C =4, A=1 28
29 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition S 11 S 12 S 13 Coil 1 v 2 S 41 v 1 S 21 S 42 v 2 Coil 2 S 22 S 43 Coil 3 v 3 S 23 Each coil measures k-space. S 31 S 32 N C =4, A=3 S 33 29
30 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition S 11 S 12 a 4 =S 41 v 1 +S 42 v 2 +S 43 v 3 a 1 =S 11 v 1 +S 12 v 2 +S 13 v 3 S 13 Coil 1 v 2 S 41 v 1 S 21 S 42 v 2 Coil 2 S 22 S 43 Coil 3 v 3 S 23 Each coil measures k-space. S 31 S 32 N C =4, A=3 a 3 =S 31 v 1 +S 43 v 2 +S 33 v 3 a 2 =S 21 v 1 +S 22 v 2 +S 23 v 3 S 33 30
31 Reconstruction SENSE Measured Coil Images Estimated Sensitivities Separated Combined Image ^ v 1 & ^v 2 ^ v 3 N C =4, A=3 31
32 Reconstruction/Processing SENSE Image vector Coil 4 Coil 1 Coil 3 Coil 2 u OP I U P CS Ok 0 u p 0 permute to by folded voxel reconstruct N c =4 images k-space vector of N c images 32
33 Implications In statistics, we know the rule that says: If a vector f has a mean δ, and a covariance Γ, Then y=of has a mean μ=oδ, and a covariance Σ=OΓO T. Then Σ can converted into a correlation matrix R=D -1/2 ΣD -1/2. 1/2 Where D 1/ diag( ). Assume k-space measurements independent so Γ=I. 33
34 Implications SENSE induces long-range in-plane correlation. fold Theoretical Results SENSE A=3 smoothed R-R I-I R-I M 2 +1 Center voxel fold -1 Basically multiplying voxel values a t by same 3 numbers over time t to lead to correlated voxels. a b c 34
35 Implications SENSE Reconstruction induces long-range correlation. Experimental Results SENSE A=3 smoothed +1 R-R I-I R-I M 2 R-R I-I R-I M 2 Karaman,Nencka,Bruce,Rowe: Brain Connect, 4: ,
36 Implications GRAPPA reconstruction induces long-range correlation. Experimental Results GRAPPA A=3 smoothed +1 R-R I-I R-I M 2 R-R I-I R-I M 2 Bruce: PhD Dissertation,
37 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition Simultaneous Multi-Slice (SMS) Coil 1 v 2 v 3 v 1 Coil 2 Coil 3 Each coil measures k-space. v j N C =4, A=3 j=1:a Larkman et al: JMRI, 13: ,
38 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition Simultaneous Multi-Slice (SMS) S 11 S 12 S 13 Coil 1 S 43 v 3 S 23 S 42 S 41 v 2 v 1 Coil 2 S 21 S 22 Coil 3 Each coil measures k-space. S 31 S 32 S 33 N C =4, A=3 v j,s kj j=1:a, k=1:n c 38
39 Coil 4 BIRS Neuroimaging Data Analysis Coil local so non uniform sensitivity. Reconstruction Multi-Coil Acquisition Simultaneous Multi-Slice (SMS) S 11 S 12 S 13 a 4 =S 41 v 1 +S 42 v 2 +S 43 v 3 a 1 =S 11 v 1 +S 12 v 2 +S 13 v 3 Coil 1 S 41 S 42 S 43 v 1 v 2 v 3 Coil 2 S 21 S 22 S 23 Coil 3 Each coil measures k-space. a 3 =S 31 v 1 +S 43 v 2 +S 33 v 3 a 2 =S 21 v 1 +S 22 v 2 +S 23 v 3 S 31 S 32 S 33 N C =4, A=3 v j,s kj,a k j=1:a, k=1:n c 39
40 Reconstruction SENSE SMS Measured Coil Images Estimated Sensitivities Separated Combined Images ^ v 1 & ^ v 2 ^ v 3 N C =4, A=3 40
41 Reconstruction/Processing SENSE SMS Image vector 4 u OP I U P CS Ok 0 u p 0 permute to by folded voxel reconstruct N c =4 images k-space vector of N c images 41
42 Implications In statistics, we know the rule that says: If a vector f has a mean δ, and a covariance Γ, Then y=of has a mean μ=oδ, and a covariance Σ=OΓO T. Then Σ can converted into a correlation matrix R=D -1/2 ΣD -1/2. 1/2 Where D 1/ diag( ). Assume k-space measurements independent so Γ=I. 42
43 Coil 4 Slice 3 Coil 3 Slice 2 Coil 2 Slice 1 Coil 1 BIRS Neuroimaging Data Analysis Implications SENSE induces long-range through-plane correlation. Coil Images Separated Images R1 R1 I1 R2 I2 R3 I3 I1 R2 smooth I2 R3 I3 43
44 Coil 4 Slice 3 Coil 3 Slice 2 Coil 2 Slice 1 Coil 1 BIRS Neuroimaging Data Analysis Implications SENSE induces long-range in-plane correlation. Coil Images Separated Images R1 R1 I1 R2 I2 R3 I3 I1 R2 smooth I2 R3 I3 44
45 Coil 4 Slice 3 Coil 3 Slice 2 Coil 2 Slice 1 Coil 1 BIRS Neuroimaging Data Analysis Implications SENSE induces long-range through-plane correlation. Coil Images Separated Images 2 M1 M1 M2 M smooth 2 M2 2 M3 45
46 Coil 4 Slice 3 Coil 3 Slice 2 Coil 2 Slice 1 Coil 1 BIRS Neuroimaging Data Analysis Implications SENSE induces long-range in-plane correlation. Coil Images Separated Images 2 M1 M1 M2 M smooth 2 M2 2 M3 46
47 Discussion Care needs to be taken when we obtain data. We should get data in its originally measured form. We should do any required processing ourselves. We should develop models that incorporate processing. We should use all of the data (magnitude and phase). 47
48 Thank You! This work is joint with former & current PhD students: Dr. Andrew S. Nencka, Medical College of Wisconsin Dr. Andrew D. Hahn, University of Wisconsin-Madison Dr. Iain P. Bruce, Duke University Dr. M. Muge Karaman, University if Illinois-Chicago Ms. Mary C. Kociuba, Marquette University Mr. Kevin K. Liu, Marquette University 48
An Introduction to Image Reconstruction, Processing, and their Effects in FMRI
An Introduction to Image Reconstruction, Processing, and their Effects in FMRI Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette University
More informationA Bayesian Approach to SENSE Image Reconstruction in FMRI
A Bayesian Approach to SENSE Image Reconstruction in FMRI Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette University July 31, 017
More informationSignal and Noise in Complex-Valued SENSE MR Image Reconstruction
Signal and Noise in omplex-valued SENSE M mage econstruction Daniel B. owe, Ph.D. (Joint with ain P. Bruce) Associate Professor omputational Sciences Program Department of Mathematics, Statistics, and
More informationA TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING. Mary C. Kociuba
A TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING by Mary C. Kociuba A Thesis Submitted to the Faculty of the Graduate School, Marquette University,
More informationMarquette University Iain P. Bruce Marquette University, M. Muge Karaman Marquette University
Marquette University e-publications@marquette Mathematics, Statistics and omputer Science Faculty Research and Publications Mathematics, Statistics and omputer Science, Department of 10-1-2012 The SENSE-Isomorphism
More informationA Resampling of Calibration Images for a Multi-Coil Separation of Parallel Encoded Complex-valued Slices in fmri
*Manuscript Click here to view linked References 1 A Resampling of Calibration Images for a Multi-Coil Separation of Parallel Encoded Complex-valued Slices in fmri Mary C. Kociuba a, Andrew S. Nencka b,
More informationQuantifying the Statistical Impact of GRAPPA in fcmri Data With a Real-Valued Isomorphism
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 33, NO. 2, FEBRUARY 2014 495 Quantifying the Statistical Impact of GRAPPA in fcmri Data With a Real-Valued Isomorphism Iain P. Bruce and Daniel B. Rowe* Abstract
More informationMRI Physics II: Gradients, Imaging
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 informationAccelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing
Accelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing Peng Hu, Ph.D. Associate Professor Department of Radiological Sciences PengHu@mednet.ucla.edu 310-267-6838 MRI... MRI has low
More informationSparse sampling in MRI: From basic theory to clinical application. R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology
Sparse sampling in MRI: From basic theory to clinical application R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology Objective Provide an intuitive overview of compressed sensing
More informationSampling, Ordering, Interleaving
Sampling, Ordering, Interleaving Sampling patterns and PSFs View ordering Modulation due to transients Temporal modulations Slice interleaving Sequential, Odd/even, bit-reversed Arbitrary Other considerations:
More informationFMRI statistical brain activation from k-space data
FMRI statistical brain activation from k-space data Daniel B. Rowe, Department of Biophysics, Division of Bistatistics, Medical College of Wisconsin, Milwaukee, WI USA KEY WORDS: image reconstruction,
More informationRole of Parallel Imaging in High Field Functional MRI
Role of Parallel Imaging in High Field Functional MRI Douglas C. Noll & Bradley P. Sutton Department of Biomedical Engineering, University of Michigan Supported by NIH Grant DA15410 & The Whitaker Foundation
More informationHead motion in diffusion MRI
Head motion in diffusion MRI Anastasia Yendiki HMS/MGH/MIT Athinoula A. Martinos Center for Biomedical Imaging 11/06/13 Head motion in diffusion MRI 0/33 Diffusion contrast Basic principle of diffusion
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 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 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 informationSteen Moeller Center for Magnetic Resonance research University of Minnesota
Steen Moeller Center for Magnetic Resonance research University of Minnesota moeller@cmrr.umn.edu Lot of material is from a talk by Douglas C. Noll Department of Biomedical Engineering Functional MRI Laboratory
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 informationXI Signal-to-Noise (SNR)
XI Signal-to-Noise (SNR) Lecture notes by Assaf Tal n(t) t. Noise. Characterizing Noise Noise is a random signal that gets added to all of our measurements. In D it looks like this: while in D
More informationComputational Aspects of MRI
David Atkinson Philip Batchelor David Larkman Programme 09:30 11:00 Fourier, sampling, gridding, interpolation. Matrices and Linear Algebra 11:30 13:00 MRI Lunch (not provided) 14:00 15:30 SVD, eigenvalues.
More informationCombination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T
Combination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T DEDALE Workshop Nice Loubna EL GUEDDARI (NeuroSPin) Joint work with: Carole LAZARUS, Alexandre VIGNAUD and Philippe
More informationMRI image formation 8/3/2016. Disclosure. Outlines. Chen Lin, PhD DABR 3. Indiana University School of Medicine and Indiana University Health
MRI image formation Indiana University School of Medicine and Indiana University Health Disclosure No conflict of interest for this presentation 2 Outlines Data acquisition Spatial (Slice/Slab) selection
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 informationSampling, Ordering, Interleaving
Sampling, Ordering, Interleaving Sampling patterns and PSFs View ordering Modulation due to transients Temporal modulations Timing: cine, gating, triggering Slice interleaving Sequential, Odd/even, bit-reversed
More informationParallel Imaging. Marcin.
Parallel Imaging Marcin m.jankiewicz@gmail.com Parallel Imaging initial thoughts Over the last 15 years, great progress in the development of pmri methods has taken place, thereby producing a multitude
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 informationMotion Correction in fmri by Mapping Slice-to-Volume with Concurrent Field-Inhomogeneity Correction
Motion Correction in fmri by Mapping Slice-to-Volume with Concurrent Field-Inhomogeneity Correction Desmond T.B. Yeo 1,2, Jeffery A. Fessler 2, and Boklye Kim 1 1 Department of Radiology, University of
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 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 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 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 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 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 informationModule 4. K-Space Symmetry. Review. K-Space Review. K-Space Symmetry. Partial or Fractional Echo. Half or Partial Fourier HASTE
MRES 7005 - Fast Imaging Techniques Module 4 K-Space Symmetry Review K-Space Review K-Space Symmetry Partial or Fractional Echo Half or Partial Fourier HASTE Conditions for successful reconstruction Interpolation
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 informationBayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri
Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri Galin L. Jones 1 School of Statistics University of Minnesota March 2015 1 Joint with Martin Bezener and John Hughes Experiment
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 informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
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 informationPartial k-space Recconstruction
Partial k-space Recconstruction John Pauly September 29, 2005 1 Motivation for Partial k-space Reconstruction a) Magnitude b) Phase In theory, most MRI images depict the spin density as a function of position,
More informationsurface Image reconstruction: 2D Fourier Transform
2/1/217 Chapter 2-3 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 informationG Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing
G16.4428 Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine Compressed Sensing Ricardo Otazo, PhD ricardo.otazo@nyumc.org Compressed
More informationDynamic Autocalibrated Parallel Imaging Using Temporal GRAPPA (TGRAPPA)
Magnetic Resonance in Medicine 53:981 985 (2005) Dynamic Autocalibrated Parallel Imaging Using Temporal GRAPPA (TGRAPPA) Felix A. Breuer, 1 * Peter Kellman, 2 Mark A. Griswold, 1 and Peter M. Jakob 1 Current
More informationFast Imaging UCLA. Class Business. Class Business. Daniel B. Ennis, Ph.D. Magnetic Resonance Research Labs. Tuesday (3/7) from 6-9pm HW #1 HW #2
Fast Imaging Daniel B. Ennis, Ph.D. Magnetic Resonance Research Labs Class Business Tuesday (3/7) from 6-9pm 6:00-7:30pm Groups Avanto Sara Said, Yara Azar, April Pan Skyra Timothy Marcum, Diana Lopez,
More informationK-Space Trajectories and Spiral Scan
K-Space and Spiral Scan Presented by: Novena Rangwala nrangw2@uic.edu 1 Outline K-space Gridding Reconstruction Features of Spiral Sampling Pulse Sequences Mathematical Basis of Spiral Scanning Variations
More informationFast Imaging Trajectories: Non-Cartesian Sampling (1)
Fast Imaging Trajectories: Non-Cartesian Sampling (1) M229 Advanced Topics in MRI Holden H. Wu, Ph.D. 2018.05.03 Department of Radiological Sciences David Geffen School of Medicine at UCLA Class Business
More informationAdvanced Imaging Trajectories
Advanced Imaging Trajectories Cartesian EPI Spiral Radial Projection 1 Radial and Projection Imaging Sample spokes Radial out : from k=0 to kmax Projection: from -kmax to kmax Trajectory design considerations
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 informationPartial k-space Reconstruction
Chapter 2 Partial k-space Reconstruction 2.1 Motivation for Partial k- Space Reconstruction a) Magnitude b) Phase In theory, most MRI images depict the spin density as a function of position, and hence
More informationHigh Spatial Resolution EPI Using an Odd Number of Interleaves
Magnetic Resonance in Medicine 41:1199 1205 (1999) High Spatial Resolution EPI Using an Odd Number of Interleaves Michael H. Buonocore* and David C. Zhu Ghost artifacts in echoplanar imaging (EPI) arise
More informationarxiv: v2 [physics.med-ph] 22 Jul 2014
Multichannel Compressive Sensing MRI Using Noiselet Encoding arxiv:1407.5536v2 [physics.med-ph] 22 Jul 2014 Kamlesh Pawar 1,2,3, Gary Egan 4, and Jingxin Zhang 1,5,* 1 Department of Electrical and Computer
More informationNUFFT for Medical and Subsurface Image Reconstruction
NUFFT for Medical and Subsurface Image Reconstruction Qing H. Liu Department of Electrical and Computer Engineering Duke University Duke Frontiers 2006 May 16, 2006 Acknowledgment Jiayu Song main contributor
More informationSIEMENS MAGNETOM Skyra syngo MR D13
Page 1 of 8 SIEMENS MAGNETOM Skyra syngo MR D13 \\USER\CIND\StudyProtocols\PTSA\*dm_ep2d_mono70_b0_p2_iso2.0 TA:1:05 PAT:2 Voxel size:2.0 2.0 2.0 mm Rel. SNR:1.00 :epse Properties Routine Prio Recon Load
More informationFunctional magnetic resonance imaging brain activation directly from k-space
Available online at www.sciencedirect.com Magnetic Resonance Imaging 7 (009) 1370 1381 Functional magnetic resonance imaging brain activation directly from k-space Daniel B. Rowe a,b,, Andrew D. Hahn a,
More informationThis Time. fmri Data analysis
This Time Reslice example Spatial Normalization Noise in fmri Methods for estimating and correcting for physiologic noise SPM Example Spatial Normalization: Remind ourselves what a typical functional image
More informationSIEMENS MAGNETOM Verio syngo MR B17
\\USER\Dr. Behrmann\routine\Ilan\ep2d_bold_PMU_resting TA: 8:06 PAT: Voxel size: 3.03.03.0 mm Rel. SNR: 1.00 USER: ep2d_bold_pmu Properties Special sat. Prio Recon System Before measurement Body After
More informationCOBRE Scan Information
COBRE Scan Information Below is more information on the directory structure for the COBRE imaging data. Also below are the imaging parameters for each series. Directory structure: var/www/html/dropbox/1139_anonymized/human:
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 informationMotion Robust Magnetic Susceptibility and Field Inhomogeneity Estimation Using Regularized Image Restoration Techniques for fmri
Motion Robust Magnetic Susceptibility and Field Inhomogeneity Estimation Using Regularized Image Restoration Techniques for fmri Desmond Tec Beng Yeo 1,, Jeffrey A. Fessler 1,, and Bolye Kim 1 1 Department
More informationMRI Imaging Options. Frank R. Korosec, Ph.D. Departments of Radiology and Medical Physics University of Wisconsin Madison
MRI Imaging Options Frank R. Korosec, Ph.D. Departments of Radiolog and Medical Phsics Universit of Wisconsin Madison f.korosec@hosp.wisc.edu As MR imaging becomes more developed, more imaging options
More informationClinical Importance. Aortic Stenosis. Aortic Regurgitation. Ultrasound vs. MRI. Carotid Artery Stenosis
Clinical Importance Rapid cardiovascular flow quantitation using sliceselective Fourier velocity encoding with spiral readouts Valve disease affects 10% of patients with heart disease in the U.S. Most
More informationControlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA)
www.siemens.com/magnetom-world Controlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA) Felix Breuer; Martin Blaimer; Mark Griswold; Peter Jakob Answers for life. Controlled
More informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationImaging Notes, Part IV
BME 483 MRI Notes 34 page 1 Imaging Notes, Part IV Slice Selective Excitation The most common approach for dealing with the 3 rd (z) dimension is to use slice selective excitation. This is done by applying
More informationIMPROVING FMRI ANALYSIS AND MR RECONSTRUCTION WITH THE INCORPORATION OF MR RELAXIVITIES AND CORRELATION EFFECT EXAMINATION. M.
IMPROVING FMRI ANALYSIS AND MR RECONSTRUCTION WITH THE INCORPORATION OF MR RELAXIVITIES AND CORRELATION EFFECT EXAMINATION by M. Muge Karaman A Dissertation submitted to the Faculty of the Graduate School,
More informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More informationPhilips MRI Protocol Dump Created on Comment Software Stream
Page 1 of 5 Philips MRI Protocol Dump Created on 2/17/2011 4:11:01 PM Comment Created by ExamCard_to_XML with inputs: "J:\ADNI GO - ADNI 2 Phantom5.ExamCard" on system (BU SCHOOL OF MEDICINE :: 192.168.71.10)
More informationUse of Multicoil Arrays for Separation of Signal from Multiple Slices Simultaneously Excited
JOURNAL OF MAGNETIC RESONANCE IMAGING 13:313 317 (2001) Technical Note Use of Multicoil Arrays for Separation of Signal from Multiple Slices Simultaneously Excited David J. Larkman, PhD, 1 * Joseph V.
More informationEnhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr.
Enhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr. Greg Zaharchuk, Associate Professor in Radiology, Stanford
More informationSIEMENS MAGNETOM TrioTim syngo MR B17
\\USER\KNARRGROUP\MultiBand\LavretskyMultiBand\trufi localizer 3-plane TA: 5.1 s PAT: Voxel size: 1.2 1.2 5. Rel. SNR: 1.00 SIEMENS: trufi Load to stamp Slice group 1 Slices 1 Dist. factor 20 % Phase enc.
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 information(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
EE225E/BIOE265 Spring 2013 Principles of MRI Miki Lustig Assignment 8 Solutions 1. Nishimura 7.1 P)jc%z 5 ivcr3. 1. I Due Wednesday April 10th, 2013 (a Scrhon5 R2iwd b 0 ZOms Xn,s r cx > qs 4-4 8ni6 4
More informationPhase Difference Reconstruction. Outline
Advanced MRI Phase Difference Reconstruction Faik Can MERAL Outline Introduction Quantitative Description Arctangent operation ATAN2 Phased-Array Multiple Coil Data Correction of Predictable Phase Errors
More informationM R I Physics Course
M R I Physics Course Multichannel Technology & Parallel Imaging Nathan Yanasak, Ph.D. Jerry Allison Ph.D. Tom Lavin, B.S. Department of Radiology Medical College of Georgia References: 1) The Physics of
More informationA GPU Accelerated Interactive Interface for Exploratory Functional Connectivity Analysis of fmri Data
A GPU Accelerated Interactive Interface for Exploratory Functional Connectivity Analysis of fmri Data Anders Eklund, Ola Friman, Mats Andersson and Hans Knutsson Linköping University Post Print N.B.: When
More informationNon-Cartesian Parallel Magnetic Resonance Imaging
Non-Cartesian Parallel Magnetic Resonance Imaging Dissertation zur Erlangung des naturwissenschaftlichen Doktorgrades der Bayerischen Julius-Maximilians-Universität Würzburg vorgelegt von Robin Heidemann
More informationImproved 3D image plane parallel magnetic resonance imaging (pmri) method
B. Wu, R. P. Millane, R. Watts, P. J. Bones, Improved 3D Image Plane Parallel Magnetic Resonance Imaging (pmri) Method, Proceedings of Image and Vision Computing New Zealand 2007, pp. 311 316, Hamilton,
More informationDynamic Contrast enhanced MRA
Dynamic Contrast enhanced MRA Speaker: Yung-Chieh Chang Date : 106.07.22 Department of Radiology, Taichung Veterans General Hospital, Taichung, Taiwan 1 Outline Basic and advanced principles of Diffusion
More informationAssignment 2. Due Feb 3, 2012
EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Due Feb 3, 2012 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir
More informationControlled Aliasing in Volumetric Parallel Imaging (2D CAIPIRINHA)
Magnetic Resonance in Medicine 55:549 556 (2006) Controlled Aliasing in Volumetric Parallel Imaging (2D CAIPIRINHA) Felix A. Breuer,* Martin Blaimer, Matthias F. Mueller, Nicole Seiberlich, Robin M. Heidemann,
More informationCognitive States Detection in fmri Data Analysis using incremental PCA
Department of Computer Engineering Cognitive States Detection in fmri Data Analysis using incremental PCA Hoang Trong Minh Tuan, Yonggwan Won*, Hyung-Jeong Yang International Conference on Computational
More informationSlide 1. Technical Aspects of Quality Control in Magnetic Resonance Imaging. Slide 2. Annual Compliance Testing. of MRI Systems.
Slide 1 Technical Aspects of Quality Control in Magnetic Resonance Imaging Slide 2 Compliance Testing of MRI Systems, Ph.D. Department of Radiology Henry Ford Hospital, Detroit, MI Slide 3 Compliance Testing
More informationADVANCED RECONSTRUCTION TECHNIQUES IN MRI - 2
ADVANCED RECONSTRUCTION TECHNIQUES IN MRI - 2 Presented by Rahil Kothari PARTIAL FOURIER RECONSTRUCTION WHAT IS PARTIAL FOURIER RECONSTRUCTION? In Partial Fourier Reconstruction data is not collected symmetrically
More informationCompressed Sensing And Joint Acquisition Techniques In Mri
Wayne State University Wayne State University Theses 1-1-2013 Compressed Sensing And Joint Acquisition Techniques In Mri Rouhollah Hamtaei Wayne State University, Follow this and additional works at: http://digitalcommons.wayne.edu/oa_theses
More informationAdvanced methods for image reconstruction in fmri
Advanced methods for image reconstruction in fmri Jeffrey A. Fessler EECS Department The University of Michigan Regional Symposium on MRI Sep. 28, 27 Acknowledgements: Doug Noll, Brad Sutton, Outline MR
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 informationEffect of age and dementia on topology of brain functional networks. Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand
Effect of age and dementia on topology of brain functional networks Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand 1 Structural changes in aging brain Age-related changes in
More informationFollowing on from the two previous chapters, which considered the model of the
Chapter 5 Simulator validation Following on from the two previous chapters, which considered the model of the simulation process and how this model was implemented in software, this chapter is concerned
More informationAcknowledgments and financial disclosure
AAPM 2012 Annual Meeting Digital breast tomosynthesis: basic understanding of physics principles James T. Dobbins III, Ph.D., FAAPM Director, Medical Physics Graduate Program Ravin Advanced Imaging Laboratories
More information3D MAGNETIC RESONANCE IMAGING OF THE HUMAN BRAIN NOVEL RADIAL SAMPLING, FILTERING AND RECONSTRUCTION
3D MAGNETIC RESONANCE IMAGING OF THE HUMAN BRAIN NOVEL RADIAL SAMPLING, FILTERING AND RECONSTRUCTION Maria Magnusson 1,2,4, Olof Dahlqvist Leinhard 2,4, Patrik Brynolfsson 2,4, Per Thyr 2,4 and Peter Lundberg
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 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 informationA Study of Nonlinear Approaches to Parallel Magnetic Resonance Imaging
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations December 2012 A Study of Nonlinear Approaches to Parallel Magnetic Resonance Imaging Yuchou Chang University of Wisconsin-Milwaukee
More informationField Maps. 1 Field Map Acquisition. John Pauly. October 5, 2005
Field Maps John Pauly October 5, 25 The acquisition and reconstruction of frequency, or field, maps is important for both the acquisition of MRI data, and for its reconstruction. Many of the imaging methods
More informationDIGITAL TERRAIN MODELS
DIGITAL TERRAIN MODELS 1 Digital Terrain Models Dr. Mohsen Mostafa Hassan Badawy Remote Sensing Center GENERAL: A Digital Terrain Models (DTM) is defined as the digital representation of the spatial distribution
More informationSIEMENS MAGNETOM Verio syngo MR B15V
\\USER\ZAHID_RESEARCH\MS\No Name\3D SWI TA: 6:39 PAT: 2 Voxel size: 1.0 0.5 2.0 mm Rel. SNR: 1.00 SIEMENS: gre Properties Prio Recon Before measurement After measurement Load to viewer Inline movie Auto
More informationDevelopment of fast imaging techniques in MRI From the principle to the recent development
980-8575 2-1 2012 10 13 Development of fast imaging techniques in MRI From the principle to the recent development Yoshio MACHIDA and Issei MORI Health Sciences, Tohoku University Graduate School of Medicine
More informationA Support-Based Reconstruction for SENSE MRI
Sensors 013, 13, 409-4040; doi:10.3390/s13040409 Article OPEN ACCESS sensors ISSN 144-80 www.mdpi.com/journal/sensors A Support-Based Reconstruction for SENSE MRI Yudong Zhang, Bradley S. Peterson and
More informationJoint SENSE Reconstruction for Faster Multi-Contrast Wave Encoding
Joint SENSE Reconstruction for Faster Multi-Contrast Wave Encoding Berkin Bilgic1, Stephen F Cauley1, Lawrence L Wald1, and Kawin Setsompop1 1Martinos Center for Biomedical Imaging, Charlestown, MA, United
More information