EE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging. Martin Uecker
|
|
- Hillary Andrews
- 5 years ago
- Views:
Transcription
1 EE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging Martin Uecker
2 Tentative Syllabus 01: Jan 27 Introduction 02: Feb 03 Parallel Imaging as Inverse Problem 03: Feb 10 Iterative Reconstruction Algorithms : Feb 17 (holiday) 04: Feb 24 Non-Cartesian MRI : Mar 03 (cancelled) 05: Mar 10 GRAPPA/SPIRiT 06: Mar 17 Nonlinear Inverse Reconstruction : Mar 24 (spring recess) 08: Mar 31 SAKE/ESPIRiT 09: Apr 07 Model-based Reconstruction 10: Apr 14 Compressed Sensing 11: Apr 21 Compressed Sensing 12: Apr 28 TBA
3 Parallel MRI Goal: Reduction of measurement time I Subsampling of k-space I Simultaneous acquisition with multiple receive coils I Coil sensitivities provide spatial information I Compensation for missing k-space data
4 Parallel MRI: Undersampling Undersampling Aliasing k phase k read k partition k phase
5 Image Reconstruction as Inverse Problem Forward problem: y = Fx + n x image (and more), F (nonlinear) operator, n noise, y data Regularized solution: Advantages: x = argmin Fx y 2 2 }{{} + αr(x) }{{} data consistency regularization Simple extension to non-cartesian trajectories Modelling of physical effects Prior knowledge via suitable regularization terms
6 Parallel MRI as Inverse Problem Signal from multiple coils (image ρ, sensitivities c j ): s j (t) = d x ρ( x)c j ( x)e i x k(t) Assumption: known sensitivities c j linear relation between image ρ and data y Image reconstruction is a linear inverse problem: V Aρ = y Ra and Rim, Magn Reson Med, 30: (1993) Pruessmann et al., Magn Reson Med, 4: (1999)
7 Autocalibration: Motivation Estimation of the coil sensitivities during the measurement Motivation: Object influences sensitivities (dielectric) Problems with consistency due to movement (e.g. breathing) Example: Movement (R = 2x2) (Sensitivities calibrated in one position and used in the other)
8 Autocalibration: Conventional Methods Method: Complete acquisition of the k-space center Reconstruction of the fully sampled center Division by RSS image (removal of the image content) Postprocessing: smoothing, extrapolation,... k y k x (GRAPPA: Estimation of weights)
9 Autocalibration: Conventional Methods Problems: Artifacts due to erraneous calibration Limitation to small acceleration factors Need for a high number of reference lines SENSE/auto GRAPPA Acceleration factor: 4=2x2, reference lines: 16x16
10 Purpose Problems of previous methods: Calibration only uses part of the data Inaccurate removal of image content (image not known) Goals: Optimal use of all available data Consistent estimation of image and sensitivities New approach necessary
11 Parallel MRI as Nonlinear Inverse Problem New approach: Joint estimation of image and sensitivities from all data Need to solve signal equation for ρ and c j at the same time: s j (t) = d x ρ( x)c j ( x) e i x k(t) V }{{} unknown Ying L, Sheng J. Magn Reson Med, 57: , 2007 Uecker M, Hohage T, Block KT, Frahm J, Magn Reson Med, 60: , 2008
12 Parallel MRI as Nonlinear Inverse Problem Challenges: 1. Inversion of a nonlinear equation F ( x) + n = y 2. Equations are underdetermined d x c j( x) a( x)ρ( x) e i x k(t) a( x) V
13 Newton-Method How does one solve a non-linear equation? Solve: f (x) y = 0 Idea: linearize f (x + x) df (x) x + f (x) Algorithm: until f (x n ) y < ɛ df (x n ) x = y f (x n ) x n+1 x n + x x 0 y
14 Newton-Method How does one solve a non-linear equation? Solve: f (x) y = 0 Idea: linearize f (x + x) df (x) x + f (x) Algorithm: until f (x n ) y < ɛ df (x n ) x = y f (x n ) x n+1 x n + x x 0 y
15 Newton-Method How does one solve a non-linear equation? Solve: f (x) y = 0 Idea: linearize f (x + x) df (x) x + f (x) Algorithm: until f (x n ) y < ɛ df (x n ) x = y f (x n ) x n+1 x n + x x 0 x 1 y
16 Newton-Method How does one solve a non-linear equation? Solve: f (x) y = 0 Idea: linearize f (x + x) df (x) x + f (x) Algorithm: until f (x n ) y < ɛ df (x n ) x = y f (x n ) x n+1 x n + x x 0 x 1 y
17 Newton-Method How does one solve a non-linear equation? Solve: f (x) y = 0 Idea: linearize f (x + x) df (x) x + f (x) Algorithm: until f (x n ) y < ɛ df (x n ) x = y f (x n ) x n+1 x n + x x 0 x 1 x 2 y
18 Newton-Method (Quadratic) convergence: Initial guess close to solution Continuously differentiable Non-zero derivative (Second derivative exists)
19 Newton-Method Nonlinear system of equations: F : R n R m Taylor series: (DF : Jacobian) Fx = y Newton method: y = Fx n+1 = Fx n + DF (x n+1 x n ) + DF (x n+1 x n ) = y Fx n Banach spaces: Fréchet derivative Convergence: Kantorovich theorem
20 Newton for Optimization Minimization of non-linear least-squares problem: φ(x) = 1 2 Fx y 2 2 = 1 2 (Fx y)t (Fx y) Optimality at the solution: φ = 0 DF T (Fx y) = 0 Newton Method: DF T (Fx y) + H (x n+1 x n ) = 0 Hessian: H = DF T DF + terms with higher derivatives (Note: Approximated Hessian is positive semi-definite)
21 Gauss-Newton Method: Landweber Method: DF H DF (x n+1 x n ) = DF H (y Fx n ) x n+1 x n = µdf H (y Fx n ) Levenberg-Marquardt Algorithm: ( ) DF H DF + αi (x n+1 x n ) = DF H (y Fx n ) Linear optimization problem: x n+1 x n = argmin x DF x + Fx n y α x 2 2
22 Iteratively Regularized Gauss-Newton Method (IRGNM) Update rule: x n+1 x n = argmin x DF x + Fx n y α n x + x n 2 2 Similar to Levenberg-Marquardt algorithm Optimization of: Fx y α x 2 2 Iterative reduction of regularization: ( ) 2 n α n = α 0 3
23 Iteratively Regularized Gauss-Newton Method (IRGNM) Solve: F (x) y = 0 F (x + x) DF (x) x + F (x) Algorithm Newton Normal equations Regularization (Number of iterations determines regularization) while F (x n ) y 2 > ɛ Solve (with CG): ( DF H (x n )DF (x n ) + α n I ) x = DF H (x n ) (y F (x n )) + α n (x n x 0 ) x n+1 x n + x ( ) n 2 α n = α 0 3
24 Iteratively Regularized Gauss-Newton Method (IRGNM) Solve: F (x) y = 0 F (x + x) DF (x) x + F (x) Algorithm Newton Normal equations Regularization (Number of iterations determines regularization) while F (x n ) y 2 > ɛ Solve (with CG): ( DF H (x n )DF (x n ) + α n I ) x = DF H (x n ) (y F (x n )) + α n (x n x 0 ) x n+1 x n + x ( ) n 2 α n = α 0 3
25 Iteratively Regularized Gauss-Newton Method (IRGNM) Solve: F (x) y = 0 F (x + x) DF (x) x + F (x) Algorithm Newton Normal equations Regularization (Number of iterations determines regularization) while F (x n ) y 2 > ɛ Solve (with CG): ( DF H (x n )DF (x n ) + α n I ) x = DF H (x n ) (y F (x n )) + α n (x n x 0 ) x n+1 x n + x ( ) n 2 α n = α 0 3
26 Iteratively Regularized Gauss-Newton Method (IRGNM) Solve: F (x) y = 0 F (x + x) DF (x) x + F (x) Algorithm Newton Normal equations Regularization (Number of iterations determines regularization) while F (x n ) y 2 > ɛ Solve (with CG): ( DF H (x n )DF (x n ) + α n I ) x = DF H (x n ) (y F (x n )) + α n (x n x 0 ) x n+1 x n + x ( ) n 2 α n = α 0 3
27 Parallel Imaging Direct application to parallel imaging not possible! I System of equations F (ρ, cj ) = ~y is underdetermined. c If (ρ, cj ) solution, then also (ρ a, aj ): Z F (ρ/a, cj a) = I d~z cj (~z ) ~ a(~z )ρ(~z ) e i k ~z = F (ρ, ci ) a(~z ) Every image is solution with certain sensitivities! A constraint to reasonable solutions is necessary Prior knowledge: Only very smooth sensitivities are possible
28 Adaptation of the Regularization So far: Regularization used to suppress noise New: Regularization suppresses variation in the sensitivities Mathematical realization: Multiplication with weight matrix W Solution of transformed system: F = F W Choice of weight matrix: Image component of W : identity Coil components of W : penalization of high derivatives (e.g. differential operator (I ) l/2 )
29 Nonlinear Inversion The signal equation for unknown image ρ and unknown coil sensitivities c j is a nonlinear equation Fx = y. Forward operator: F : H l ([0, 1] 3, C N ) L 2 ([0, 1] 3, C) L 2 (range( k), C N ) (c j, ρ) y j = d x c j ( x)ρ( x) e i k(t) x Reconstruction: Iteratively regularized Gauss-Newton method (IRGNM) Smoothness penalty for the coil sensitivities: ρ (1 + s k 2 ) l FTc j 2 2
30 Nonlinear Inversion Algorithm: I Initialization: ρ = 1 and cj = 0 I Update rule (IRGNM): (DF (xn )H DF (xn ) + αn I )δx = DF (xn )H (y F (xn )) + αn (x0 xn ) (solved with the conjugate gradient algorithm) Regularization: I αn = q n α0, e.g. q = 1/2 I αn kρk22 + αn k(1 + s ~k 2 )l FTcj k22 (smoothness of the sensitivities) I Implementation: multiplication with a weighting matrix: F 0 = F W Uecker et al., Magn Reson Med 60: (2008)
31 Nonlinear Inversion Example: Siemens Tim Trio 3 T, 12-channel head coil 3D-FLASH, acceleration R = 2 2 iterative reconstruction of image and coil sensitivities Uecker et al., Magn Reson Med 60: (2008)
32 Nonlinear Inversion with Non-Quadratic Regularization Iteratively Regularized Gauss Newton Method (IRGNM) x n+1 x n = argmin δx DF H (x n )δx + F (x n ) y α n R(δx + x n ) Previously: Image regularized with L 2 -norm R(x) = ρ (1 + s k 2 ) l FTc j 2 2 Now: Different regularization terms R(x) = R(ρ) + (1 + s k 2 ) l FTc j 2 2 Knoll F, Clason C, Bredies K, Uecker M, Stollberger R, Magn Reson Med, 67:34-41 (2012).
33 Nonlinear Inversion Quality of the reconstructed images can be improved Acceleration: 3 x 2 L1-Wavelet: Cohen-Daubechies-Feauveau 9/7
34 Nonlinear Inverse Reconstruction with Variational Penalties Experiments: Siemens Tim Trio 3 T, 12-channel head coil 3D FLASH, acceleration: R = 4 (pseudorandom sampling) Knoll F, Clason C, Bredies K, Uecker M, Stollberger R, Magn Reson Med, 67:34-41 (2012).
35 Nonlinear Inversion for Real-Time MRI Advantages: Can be used directly with non-cartesian data No fully-sampled calibration region needed Dynamic update of coil sensitivities Prior knowledge using regularization Disadvantage: High computational demand!
36 Real-Time MRI α Echo RF G z G x G y TR Sequence for fast low-angle shot (FLASH) MRI and interleaved radial k-space scheme. Advantages of radial sampling: Robustness to motion Tolerance to undersampling Continuous updating of image data Self calibration for parallel imaging Zhang et al., J Magn Reson Imaging, 31: , 2010.
37 Gridding Reconstruction 75 spokes, 150 ms Problems: low SNR, streaking artefacts, low temporal resolution
38 Real-Time MRI Reconstruction: Autocalibrated parallel imaging based on nonlinear inversion Algorithm extended to non-cartesian (radial) sampling Further improvements for real-time imaging gridding nonlinear inversion real-time version Figure: short-axis view of a human heart, 15 spokes (30 ms)
39 Image Reconstruction on Graphical Processing Units 30 seconds of a movie: 1500 frames, 1GB data 8x TITAN GPUs Integrated with MRI scanner (TCP/IP) Array compression: 8-12 virtual channels Reconstruction time: 42ms (24 fps) (optimal quality) Uecker et al., Magn Reson Med 63: (2010) Schätz and Uecker, In: Lecture Notes in Computer Science, 7439: (2012)
40 Numerical Implementation Forward operator: F = P k FT M FOV C W W C M FOV FT P k Weighting matrix Multiplication of image and coil sensitivities Restriction to the field of view (FOV) Fourier transformation Sampling operator At the beginning: interpolation onto Cartesian grid: M FOV FT H P H k y Each iteration: convolution with point spread function (PSF) M FOV FT H Pk H }{{ P kft } convolution M FOV Uecker et al., Magn Reson Med 63: (2010)
41 Image Reconstruction on Graphical Processing Units CPU: Preprocessing measurement data whitening and array compression interpolation onto Cartesian grid GPU: Iterative reconstruction iteratively regularized Gauss-Newton method initial guess convolutionbased conjugate gradient algorithm image and sensitivities initialization with previous frame
42 Real-Time MRI Basic nonlinear inverse reconstruction: Short-axis view of the human heart
43 Real-Time MRI β = 0 β = 0.8 Figure: short-axis view of a human heart, 15 spokes (30 ms) Improved regularization: Previous frame as prior knowledge: x βx prev 2 2 Damping factor β to avoid accumulation of errors Enhanced recovery of high frequencies
44 Real-Time MRI unfiltered filtered Figure: short-axis view of a human heart, 15 spokes (30 ms) Median Filter: Applied in the temporal domain Removes streaking artefacts Preserves sharp transitions t
45 Real-Time MRI unfiltered filtered Figure: short-axis view of a human heart, 15 spokes (30 ms) Median Filter: Applied in the temporal domain Removes streaking artefacts Preserves sharp transitions t
46 Real-Time MRI Why does the median filter remove streaking artefacts? Interleaved k-space sampling scheme median Median: invariant to reordering Removes flickering artefacts for static image content Median as L 1 minimization: x { N } x = argmin x x k x 1 x 4 x 5 x 2 x 1 x 3 k=1
47 Real-Time MRI before after Figure: short-axis view of a human heart, 15 spokes (30 ms) Image filter: Edge enhancement Denoising
48 Real-Time MRI Figure: short-axis view of a human heart, 15 spokes (30 ms) Reconstruction steps: 1. Parallel imaging with nonlinear inverse reconstruction 2. Improved regularization 3. Median filter 4. Further image enhancement
49 Real-Time MRI Experiments: Siemens Tim Trio 3 T 32 channel cardiac coil (array compression to 12 virtual channels) RF-spoiled radial FLASH Healthy volunteers Free breathing, no synchronisation to ECG Image reconstruction with eight GTX 580 GPUs (Nvidia)
50 Real-Time MRI: Movies of the Human Heart Acquisition time 50 ms (25 spokes) Spatial resolution 2.0x2.0x8 mm 3
51 Real-Time MRI: Movies of the Human Heart Acquisition time 33 ms (15 spokes) Spatial resolution 1.5x1.5x8 mm 3 Acquisition time 22 ms (11 spokes) Spatial resolution 2.0x2.0x8 mm 3
52 Real-Time MRI: Swallowing and Speaking Acquisition time 50 ms (25 spokes) Spatial resolution 1.7x1.7x10 mm 3 Acquisition time 50 ms (25 spokes) Spatial resolution 1.7x1.7x10 mm 3 Niebergall A, Zhang S, Kunay E, Keydana G, Job M, Uecker M, Frahm J. Magn Reson Med, submitted 2011.
53 Project 2: Non-Cartesian MRI Project: Implement non-cartesian reconstruction Tools: Matlab, reconstruction toolbox, python,... Deadline: Mar 21 Hand in: Working code and plots/results with description. See website for data and instructions.
54 Project 2: Non-Cartesian MRI Step 1: Computation of non-cartesian samples Compute non-uniform DFT matrix Implement nufft and its adjoint 1D and 2D versions Hints: Oversampling 2 Gaussian kernel (not optimal)
55 Project 2: Non-Cartesian MRI Step 2: Implement iterative gridding Reconstruction from non-cartesian data Use Landweber and/or CG Step 3: Implement non-cartesian SENSE Implement non-cartesian SENSE 1. JD O Sullivan. A fast sinc function gridding algorithm for Fourier inversion in computer tomography. IEEE Trans Med Imaging, 4: , JI Jackson, CH Meyer, DG Nishimura, and A Macovski. Selection of a convolution function for Fourier inversion using gridding. IEEE Trans Med Imaging, 3: , KP Pruessmann, M Weiger, P Börnert, and P Boesiger. Advances in sensitivity encoding with arbitrary k-space trajectories. Magn Reson Med, 46: , 2001.
56 Final Project Project:? Tools: Matlab, reconstruction toolbox, python,... Deadline: Apr 28 Hand in: Working code and plots/results with description.
From Image to Video: Real-time Medical Imaging with MRI
From Image to Video: Real-time Medical Imaging with MRI Sebastian Schaetz, Martin Uecker BiomedNMR Forschungs GmbH at the MPI for biophysical Chemistry, Goettingen, Germany Electrical Engineering and Computer
More informationEE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging. Martin Uecker
EE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging Martin Uecker Tentative Syllabus 01: Jan 27 Introduction 02: Feb 03 Parallel Imaging as Inverse Problem 03: Feb 10 Iterative Reconstruction
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 informationImage reconstruction using compressed sensing for individual and collective coil methods.
Biomedical Research 2016; Special Issue: S287-S292 ISSN 0970-938X www.biomedres.info Image reconstruction using compressed sensing for individual and collective coil methods. Mahmood Qureshi *, Muhammad
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 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 informationNuts & Bolts of Advanced Imaging. Image Reconstruction Parallel Imaging
Nuts & Bolts of Advanced Imaging Image Reconstruction Parallel Imaging Michael S. Hansen, PhD Magnetic Resonance Technology Program National Institutes of Health, NHLBI Declaration of Financial Interests
More informationTHE STEP BY STEP INTERACTIVE GUIDE
COMSATS Institute of Information Technology, Islamabad PAKISTAN A MATLAB BASED INTERACTIVE GRAPHICAL USER INTERFACE FOR ADVANCE IMAGE RECONSTRUCTION ALGORITHMS IN MRI Medical Image Processing Research
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 informationFOV. ] are the gradient waveforms. The reconstruction of this signal proceeds by an inverse Fourier Transform as:. [2] ( ) ( )
Gridding Procedures for Non-Cartesian K-space Trajectories Douglas C. Noll and Bradley P. Sutton Dept. of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA 1. Introduction The data collected
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 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 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 informationReal-time MRI at a resolution of 20 ms
Research Article Received: 10 February 2010, Revised: 2 June 2010, Accepted: 3 June 2010, Published online in Wiley Online Library: 27 August 2010 (wileyonlinelibrary.com) DOI:10.1002/nbm.1585 Real-time
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 informationCollaborative Sparsity and Compressive MRI
Modeling and Computation Seminar February 14, 2013 Table of Contents 1 T2 Estimation 2 Undersampling in MRI 3 Compressed Sensing 4 Model-Based Approach 5 From L1 to L0 6 Spatially Adaptive Sparsity MRI
More informationZigzag Sampling for Improved Parallel Imaging
Magnetic Resonance in Medicine 60:474 478 (2008) Zigzag Sampling for Improved Parallel Imaging Felix A. Breuer, 1 * Hisamoto Moriguchi, 2 Nicole Seiberlich, 3 Martin Blaimer, 1 Peter M. Jakob, 1,3 Jeffrey
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 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 informationSeparate Magnitude and Phase Regularization via Compressed Sensing
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 9, SEPTEMBER 2012 1713 Separate Magnitude and Phase Regularization via Compressed Sensing Feng Zhao*, Douglas C. Noll, Senior Member, IEEE, Jon-Fredrik
More informationInstitute of Cardiovascular Science, UCL Centre for Cardiovascular Imaging, London, United Kingdom, 2
Grzegorz Tomasz Kowalik 1, Jennifer Anne Steeden 1, Bejal Pandya 1, David Atkinson 2, Andrew Taylor 1, and Vivek Muthurangu 1 1 Institute of Cardiovascular Science, UCL Centre for Cardiovascular Imaging,
More informationCompressed Sensing for Rapid MR Imaging
Compressed Sensing for Rapid Imaging Michael Lustig1, Juan Santos1, David Donoho2 and John Pauly1 1 Electrical Engineering Department, Stanford University 2 Statistics Department, Stanford University rapid
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 informationConstrained Reconstruction of Sparse Cardiac MR DTI Data
Constrained Reconstruction of Sparse Cardiac MR DTI Data Ganesh Adluru 1,3, Edward Hsu, and Edward V.R. DiBella,3 1 Electrical and Computer Engineering department, 50 S. Central Campus Dr., MEB, University
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 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 informationCompressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI
1 Compressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI Kevin T. Looby klooby@stanford.edu ABSTRACT The temporal resolution necessary for dynamic contrast enhanced (DCE) magnetic resonance
More informationTradeoffs and complexities in model-based MR image reconstruction
Tradeoffs and complexities in model-based MR image reconstruction Jeffrey A. Fessler EECS Dept., The University of Michigan fessler@umich.edu December 6, 2007 Abstract Model-based methods for image reconstruction
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 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 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 informationOn the Temporal Fidelity of Nonlinear Inverse Reconstructions for Real- Time MRI The Motion Challenge
Send Orders for Reprints to reprints@benthamscience.org The Open Medical Imaging Journal, 2014, 8, 1-7 1 Open Access On the Temporal Fidelity of Nonlinear Inverse Reconstructions for Real- Time MRI The
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 informationCS 395T Lecture 12: Feature Matching and Bundle Adjustment. Qixing Huang October 10 st 2018
CS 395T Lecture 12: Feature Matching and Bundle Adjustment Qixing Huang October 10 st 2018 Lecture Overview Dense Feature Correspondences Bundle Adjustment in Structure-from-Motion Image Matching Algorithm
More informationEfficient MR Image Reconstruction for Compressed MR Imaging
Efficient MR Image Reconstruction for Compressed MR Imaging Junzhou Huang, Shaoting Zhang, and Dimitris Metaxas Division of Computer and Information Sciences, Rutgers University, NJ, USA 08854 Abstract.
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 informationarxiv: v2 [physics.med-ph] 3 Aug 2017
arxiv:1705.04135v2 [physics.med-ph] 3 Aug 2017 Simultaneous Multi-Slice MRI using Cartesian and Radial FLASH and Regularized Nonlinear Inversion: SMS-NLINV Sebastian Rosenzweig 1, H. Christian M. Holme
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 informationarxiv: v2 [cs.cv] 5 Jan 2016
Estimating Absolute-Phase Maps Using ESPIRiT and Virtual Conjugate Coils arxiv:1509.03557v2 [cs.cv] 5 Jan 2016 Martin Uecker 1 and Michael Lustig 2 1 Institute for Diagnostic and Interventional Radiology,
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 informationSpiral keyhole imaging for MR fingerprinting
Spiral keyhole imaging for MR fingerprinting Guido Buonincontri 1, Laura Biagi 1,2, Pedro A Gómez 3,4, Rolf F Schulte 4, Michela Tosetti 1,2 1 IMAGO7 Research Center, Pisa, Italy 2 IRCCS Stella Maris,
More informationk-t BLAST Reconstruction From Non-Cartesian k-t Space Sampling
Magnetic Resonance in Medicine 55:85 91 (006) k-t BLAST Reconstruction From Non-Cartesian k-t Space Sampling Michael S. Hansen, 1, * Christof Baltes, 1 Jeffrey Tsao, 3 Sebastian Kozerke, 1 Klaas P. Pruessmann,
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 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 informationSingle Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction
Single Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction Felix Lugauer 1,3, Jens Wetzl 1, Christoph Forman 2, Manuel Schneider 1, Berthold Kiefer
More informationSupplementary methods
Supplementary methods This section provides additional technical details on the sample, the applied imaging and analysis steps and methods. Structural imaging Trained radiographers placed all participants
More informationNIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.
NIH Public Access Author Manuscript Published in final edited form as: Med Phys. 2008 February ; 35(2): 660 663. Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic
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 informationA Rapid and Robust Numerical Algorithm for Sensitivity Encoding with Sparsity Constraints: Self-Feeding Sparse SENSE
A Rapid and Robust Numerical Algorithm for Sensitivity Encoding with Sparsity Constraints: Self-Feeding Sparse SENSE Feng Huang, 1 * Yunmei Chen, Wotao Yin, 3 Wei Lin, 1 Xiaojing Ye, Weihong Guo, 4 and
More informationVD-AUTO-SMASH Imaging
Magnetic Resonance in Medicine 45:1066 1074 (2001) VD-AUTO-SMASH Imaging Robin M. Heidemann, Mark A. Griswold, Axel Haase, and Peter M. Jakob* Recently a self-calibrating SMASH technique, AUTO-SMASH, was
More informationEvaluations of k-space Trajectories for Fast MR Imaging for project of the course EE591, Fall 2004
Evaluations of k-space Trajectories for Fast MR Imaging for project of the course EE591, Fall 24 1 Alec Chi-Wah Wong Department of Electrical Engineering University of Southern California 374 McClintock
More informationSpread Spectrum Using Chirp Modulated RF Pulses for Incoherent Sampling Compressive Sensing MRI
Spread Spectrum Using Chirp Modulated RF Pulses for Incoherent Sampling Compressive Sensing MRI Sulaiman A. AL Hasani Department of ECSE, Monash University, Melbourne, Australia Email: sulaiman.alhasani@monash.edu
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 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 informationDYNAMIC magnetic resonance imaging (MRI), which captures
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL 26, NO 7, JULY 2007 917 Accelerating Dynamic Spiral MRI by Algebraic Reconstruction From Undersampled k t Space Taehoon Shin*, Student Member, IEEE, Jon-Fredrik
More informationGPUs Open New Avenues in Medical MRI
GPUs Open New Avenues in Medical MRI Chris A. Cocosco D. Gallichan, F. Testud, M. Zaitsev, and J. Hennig Dept. of Radiology, Medical Physics, UNIVERSITY MEDICAL CENTER FREIBURG 1 Our research group: Biomedical
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 informationarxiv: v1 [physics.med-ph] 29 Feb 2016
arxiv:1603.00040v1 [physics.med-ph] 29 Feb 2016 Model-Based Iterative Reconstruction for Radial Fast Spin-Echo MRI Kai Tobias Block, Martin Uecker, and Jens Frahm Abstract In radial fast spin-echo MRI,
More informationSupplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging
Supplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging Paper ID: 1195 No Institute Given 1 More Experiment Results 1.1 Visual Comparisons We apply all methods on four 2D
More informationRandomized sampling strategies
Randomized sampling strategies Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers & Acquisition costs impediments Full-waveform inversion
More informationInformation about presenter
Information about presenter 2013-now Engineer R&D ithera Medical GmbH 2011-2013 M.Sc. in Biomedical Computing (TU München) Thesis title: A General Reconstruction Framework for Constrained Optimisation
More informationGRAPPA Operator for Wider Radial Bands (GROWL) with Optimally Regularized Self-Calibration
GRAPPA Operator for Wider Radial Bands (GROWL) with Optimally Regularized Self-Calibration Wei Lin,* Feng Huang, Yu Li, and Arne Reykowski Magnetic Resonance in Medicine 64:757 766 (2010) A self-calibrated
More informationParallel Magnetic Resonance Imaging (pmri): How Does it Work, and What is it Good For?
Parallel Magnetic Resonance Imaging (pmri): How Does it Work, and What is it Good For? Nathan Yanasak, Ph.D. Chair, AAPM TG118 Department of Radiology Georgia Regents University Overview Phased-array coils
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 informationHartley - Zisserman reading club. Part I: Hartley and Zisserman Appendix 6: Part II: Zhengyou Zhang: Presented by Daniel Fontijne
Hartley - Zisserman reading club Part I: Hartley and Zisserman Appendix 6: Iterative estimation methods Part II: Zhengyou Zhang: A Flexible New Technique for Camera Calibration Presented by Daniel Fontijne
More informationFast Implementation of Iterative Image Reconstruction
Fast Implementation of Iterative Image Reconstruction Tobias Schaeffter 1, Michael Schacht Hansen 2, Thomas Sangild Sørensen 3 1 Division of Imaging Sciences and Biomedical Engineering, Kings College London
More informationA Convex Optimization Approach to pmri Reconstruction
1 A Convex Optimization Approach to pmri Reconstruction Cishen Zhang and Ifat Al Baqee arxiv:1311.366v1 [physics.med-ph] 11 Nov 013 Abstract In parallel magnetic resonance imaging (pmri) reconstruction
More informationMidterm Review
Midterm Review - 2017 EE369B Concepts Noise Simulations with Bloch Matrices, EPG Gradient Echo Imaging 1 About the Midterm Monday Oct 30, 2017. CCSR 4107 Up to end of C2 1. Write your name legibly on this
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 MR Image Reconstruction Using Augmented Lagrangian Methods Sathish Ramani*, Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE
694 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 30, NO. 3, MARCH 2011 Parallel MR Image Reconstruction Using Augmented Lagrangian Methods Sathish Ramani*, Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE
More informationIMAGE reconstruction in conventional magnetic resonance
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 3, MARCH 2005 325 Conjugate Phase MRI Reconstruction With Spatially Variant Sample Density Correction Douglas C. Noll*, Member, IEEE, Jeffrey A. Fessler,
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationRedundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity
Redundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity Vimal Singh and Ahmed H. Tewfik Electrical and Computer Engineering Dept., The University of Texas at Austin, USA Abstract. For
More informationHigher Degree Total Variation for 3-D Image Recovery
Higher Degree Total Variation for 3-D Image Recovery Greg Ongie*, Yue Hu, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa ISBI 2014 Beijing, China Motivation: Compressed
More informationThe cost of parallel imaging in functional MRI of the human brain
Magnetic Resonance Imaging 24 (2006) 1 5 Original contributions The cost of parallel imaging in functional MRI of the human brain Henry Lqtcke4, Klaus-Dietmar Merboldt, Jens Frahm Biomedizinische NMR Forschungs
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 informationFast reconstruction for multichannel compressed sensing using a hierarchically semiseparable solver
Fast reconstruction for multichannel compressed sensing using a hierarchically semiseparable solver The MIT Faculty has made this article openly available. Please share how this access benefits you. Your
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 informationSecond Order Total Generalized Variation (TGV) for MRI
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz Second Order Total Generalized Variation (TGV) for MRI F. Knoll K. Bredies T. Pock
More informationComputational Acceleration for MR Image Reconstruction in Partially Parallel Imaging Xiaojing Ye*, Yunmei Chen, and Feng Huang
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 30, NO. 5, MAY 2011 1055 Computational Acceleration for MR Image Reconstruction in Partially Parallel Imaging Xiaojing Ye*, Yunmei Chen, and Feng Huang Abstract
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationAnimation Lecture 10 Slide Fall 2003
Animation Lecture 10 Slide 1 6.837 Fall 2003 Conventional Animation Draw each frame of the animation great control tedious Reduce burden with cel animation layer keyframe inbetween cel panoramas (Disney
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationIntroduction to Topics in Machine Learning
Introduction to Topics in Machine Learning Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani 1/ 27 Compressed Sensing / Sparse Recovery: Given y :=
More informationMR-Encephalography (MREG)
MR-Encephalography (MREG) C2P package Version 2.0 For SIEMENS Magnetom Verio/TRIO/Skyra/Prisma Installation and User s Guide NUMARIS/4 VB17, VD11, VD13 Jakob Assländer Pierre LeVan, Ph.D. 10.09.2014 Magnetic
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 informationFast, Iterative Image Reconstruction for MRI in the Presence of Field Inhomogeneities
178 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 2, FEBRUARY 2003 Fast, Iterative Image Reconstruction for MRI in the Presence of Field Inhomogeneities Bradley P. Sutton*, Student Member, IEEE, Douglas
More informationMRI reconstruction from partial k-space data by iterative stationary wavelet transform thresholding
MRI reconstruction from partial k-space data by iterative stationary wavelet transform thresholding Mohammad H. Kayvanrad 1,2, Charles A. McKenzie 1,2,3, Terry M. Peters 1,2,3 1 Robarts Research Institute,
More informationA practical acceleration algorithm for real-time imaging
A practical acceleration algorithm for real-time imaging The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher
More informationBlind Sparse Motion MRI with Linear Subpixel Interpolation
Blind Sparse Motion MRI with Linear Subpixel Interpolation Anita Möller, Marco Maaß, Alfred Mertins Institute for Signal Processing, University of Lübeck moeller@isip.uni-luebeck.de Abstract. Vital and
More informationSparse Reconstruction / Compressive Sensing
Sparse Reconstruction / Compressive Sensing Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani Sparse Reconstruction / Compressive Sensing 1/ 20 The
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 informationarxiv: v3 [math.na] 14 Jun 2018
LAP: A LINEARIZE AND PROJECT METHOD FOR SOLVING INVERSE PROBLEMS WITH COUPLED VARIABLES JAMES L. HERRING, JAMES G. NAGY, AND LARS RUTHOTTO arxiv:1705.09992v3 [math.na] 14 Jun 2018 Abstract. Many inverse
More informationAdaptive Kernel Regression for Image Processing and Reconstruction
Adaptive Kernel Regression for Image Processing and Reconstruction Peyman Milanfar* EE Department University of California, Santa Cruz *Joint work with Sina Farsiu, Hiro Takeda AFOSR Sensing Program Review,
More information2D spatially selective excitation pulse design and the artifact evaluation
EE 591 Project 2D spatially selective excitation pulse design and the artifact evaluation 12/08/2004 Zungho Zun Two-dimensional spatially selective excitation is used to excite a volume such as pencil
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 informationSparse MRI: The Application of Compressed Sensing for Rapid MR Imaging
Sparse MRI: The Application of Compressed Sensing for Rapid MR Imaging Michael Lustig, David Donoho and John M. Pauly Magnetic Resonance Systems Research Laboratory, Department of Electrical Engineering,
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationThe Benefit of Tree Sparsity in Accelerated MRI
The Benefit of Tree Sparsity in Accelerated MRI Chen Chen and Junzhou Huang Department of Computer Science and Engineering, The University of Texas at Arlington, TX, USA 76019 Abstract. The wavelet coefficients
More information