G Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing
|
|
- Luke Summers
- 6 years ago
- Views:
Transcription
1 G Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine Compressed Sensing Ricardo Otazo, PhD
2 Compressed sensing: The big picture Exploit image sparsity/compressibility to reconstruct undersampled data Original Gradient (sparse) Which image would require fewer samples? Nyquist: same FOV, same number of samples Common sense: fewer non-zero pixels, fewer samples
3 Image compression Essential tool for modern data storage and transmission Exploit pixel correlations to reduce number of bits First reconstruct, then compress Fully-sampled acquisition (Nyquist rate) Sparsifying transform (e.g. wavelets) Recover image from sparse coefficients 3 MB 300 kb Store or transmit non-zero coefficients only 10-fold compression
4 Image compression is inefficient Question for image compression Why do we need to acquire samples at the Nyquist rate if we are going to throw away most of them? Answer We don t. Do compressed sensing instead Build data compression in the acquisition First compress, then reconstruct Candès E,Romberg J,Tao T. IEEE Trans Inf Theory 2006; 52(2): Donoho D. IEEE Trans Inf Theory 2006; 52(4):
5 Compressed sensing components Sparsity Represent images with a few coefficients Transform: wavelets, gradient, etc. Incoherence Sparsity is preserved Noise-like aliasing artifacts Non-linear reconstruction Denoise aliasing artifacts
6 Simple compressed sensing example k-space No undersampling Image space Regular undersampling Coherent aliasing Random undersampling Incoherent aliasing Courtesy of Miki Lustig
7 Simple compressed sensing example Undersampled signal Threshold - Threshold Final sparse solution Courtesy of Miki Lustig
8 Sparsifying transforms Finite differences y ( n) = x( n) x( n 1) Inmatrix form: y = Dx 1 1 D = Total variation TV N ( x) = x( n) x( n 1 ) mintv ( x) = min Dx 1 n= 2
9 Sparsifying transforms Wavelets Space-frequency localization More efficient than sines and cosines to represent discontinuities
10 Sparsifying transforms Wavelets Multiresolution image representation Recursive application of the wavelet function modified by the scaling function (each resolution is twice of that of the previous scale) Brain image Daubechies 4-tap wavelet transform
11 Incoherence Incoherence property: the acquisition domain and the sparse domain must have low correlation E: acquisition matrix W: sparsifying transform matrix <a i,w j > is small for all basis functions
12 Transform Point Spread Function (TPSF) Encoding model: s = Em Representation model: m = Wp (p is sparse) s = EWp Transform point spread function for position r TPSF( r) = W H E H EW ( r )
13 Transform Point Spread Function (TPSF) Tool to check incoherence in the sparse domain Computation Apply inverse FFT to sampling mask (1=sampled, 0=otherwise) Apply sparsifying transform Incoherence = ratio of the main peak to the std of the pseudo-noise (incoherent artifacts) Sampling mask PSF TPSF (Wavelet)
14 Incoherence of k-space trajectories Lustig M, Donoho DL, Santos JM, Pauly JM. Compressed Sensing. MRI, IEEE Signal Processing Magazine, 2008; 25(2): 72-82
15 Incoherent sampling patterns What are the best samples?? Sampling pattern Image domain Sparse domain Random undersampling (Cartesian) Variable-density random undersampling (Cartesian) Regular undersampling (radial)
16 Dimensionality and incoherence Sampling pattern (R=4) PSF 64x64 space I=71 I= x128 space
17 Compressed sensing reconstruction Acquisition model: s = Em s : acquired data E : undersampled Fourier transform m : image to reconstruct Sparsifying transform: W min Wm 1 m subject to Em = s Data consistency Sparsity N x = (l 1 - norm of x) 1 n= 1 x i
18 Compressed sensing reconstruction Unconstrained optimization (in practice) min m Em s λ Wm 1 Regularization parameter λ Trade-off between data fidelity and removal of aliasing artifacts High λ: artifact removal and denoising at the expense of image corruption (blurring, ringing, blocking, etc) Low λ : no image corruption, but residual aliasing
19 Compressed sensing reconstruction Gradient descent Cost function: Iterations: m C( m) = Em s + λ Wm n+ 1 = mn α C( mn 2 2 ) 1 C( m n ) = E H ( ) H 1 Emn s + λw M Wmn M is a diagonal matrix: Classwork: compute the gradient of C(m) * Hint: x = x x + µ
20 Compressed sensing reconstruction Proximal gradient descent (iterative soft-thresholding) Soft-thresholding operation 4 3 S( x, λ) = 0, if x x x λ ( x λ), if x > λ S(x, 1 ) x Our gradient gradient-descent algorithm becomes m [ ( [ H S W m E ( Em )],λ 1 n+ 1 = W n n s
21 Iterative soft-thresholding Sparse signal denoising Original With noise
22 Iterative soft-thresholding Sparse signal denoising Original After 5 iterations
23 Iterative soft-thresholding Sparse signal denoising Original After 15 iterations
24 Iterative soft-thresholding in Matlab function m=cs_ista(s,lambda,nite) % W and Wi are the forward and inverse sparsifying % transforms m = ifft2_mri(s); % initial solution sm = s~=0; % sampling mask for ite=1:nite, % enforce sparsity m = Wi (SoftT(W(m),lambda)); % enforce data consistency m = m ifft2_mri( (fft2_mri(m).*sm s).*sm ); end % soft-thresholding function function y=softt(x,p) y = (abs(x)-p).*x./abs(x).*(abs(x)>p); end
25 Iterative soft-thresholding Initial solution Iterations Inverse FT of the zero-filled k-space Sparse representation S(, λ) w i After soft-thresholding After data consistency k-space representation
26 Iterative soft-thresholding (R=3) Initial solution Inverse FT of the zero-filled k-space After 30 iterations Iterative soft-thresholding
27 How many samples do we need in practice? ~4K samples for a K-sparse image N = 6.5x10 4 K = 6.5x10 3 Compression ratio C = 10 R = 2.5 (safe undersampling factor) Rule of thumb: about 4 incoherent samples for each sparse coefficient
28 Image quality in compressed sensing SNR is not a good metric Loss of small coefficients in the sparse domain Loss of contrast Blurring Blockiness Ringing Images look more synthetic R=2 R=4
29 Combination of compressed sensing and parallel imaging
30 Why would CS & PI make sense? Image sparsity and coil-sensitivity encoding are complementary sources of information Compressed sensing can regularize the inverse problem in parallel imaging Parallel imaging can reduce the incoherent aliasing artifacts
31 Why would CS & PI make no sense? CS requires irregular k-space sampling while PI requires regular k-space sampling Non-linear reconstruction in CS can prevent high g- factors due to irregular sampling Incoherence along the coil dimension is very limited because the sampling pattern for all coils is the same It does not matter, because we are not exploiting incoherence along the coil dimension
32 Approaches for CS&PI CS with SENSE parallel imaging model Muticoil radial imaging with TV minimization Block T et al. MRM 2007 SPARSE-SENSE Liang D et al. MRM 2009 Otazo R et al. MRM 2010 CS with GRAPPA parallel imaging model l 1 -SPIRiT Lustig M et al. MRM 2010
33 CS with SENSE parallel imaging model Enforce joint sparsity on the multicoil system coil sensitivities E sparsifying transform W Coil 1 s 1 Coil 2 s 2 s CS & PI m Coil N s N min m { 2 Em - s + λ Wm } 2 1 SENSE data consistency Joint sparsity m represents the contribution from all coils
34 How does PI help CS? Reduces incoherent aliasing artifacts Sampling pattern (R=4) Point spread function (PSF) 1 coil 8 coils k y k x y y
35 How does PI help CS? Reduces incoherent aliasing artifacts 1 coil 8 coils Which one would be easier to reconstruct?
36 Limits of Performance Compressed sensing alone ~4K samples for a K-sparse image Compressed sensing & parallel imaging Close to K samples for large number of coils Noise-free simulation K = 32 N c = number of coils Otazo et al. ISMRM 2009; 378 UISNR coil geometry
37 CS & PI for 2D imaging Siemens 3T Tim Trio 12-channel matrix coil array 4-fold acceleration GRAPPA Coil-by-coil CS Joint CS & PI
38 Application to musculoskeletal imaging Conventional clinical protocol for knee imaging Multislice 2D-FSE (axial) Multislice 2D-FSE (coronal) minutes Multislice2D-FSE (sagittal) New compressed sensing protocol (Sparse-SPACE) 3D-FSE 4-5 minutes NYU-Siemens collaboration to be presented at RSNA 2014
39 Application to musculoskeletal imaging First clinical study 50 patients >90% agreement between conventional 2D-FSE and Sparse- SPACE 2D-FSE (10:56 min) Sparse-SPACE (4:36 min) NYU-Siemens collaboration to be presented at RSNA 2014
40 Summary Compressed sensing New sampling theorem Information rate (sparsity) rather than pixel rate (bandwidth) Ingredients Sparsity Incoherence Non-linear reconstruction Fast imaging tool for MRI MR images are naturally compressible Data acquisition in k-space facilitates incoherence Can be combined with parallel imaging
Sparse 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 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 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 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 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 informationOptimal Sampling Geometries for TV-Norm Reconstruction of fmri Data
Optimal Sampling Geometries for TV-Norm Reconstruction of fmri Data Oliver M. Jeromin, Student Member, IEEE, Vince D. Calhoun, Senior Member, IEEE, and Marios S. Pattichis, Senior Member, IEEE Abstract
More information29 th NATIONAL RADIO SCIENCE CONFERENCE (NRSC 2012) April 10 12, 2012, Faculty of Engineering/Cairo University, Egypt
K1. High Performance Compressed Sensing MRI Image Reconstruction Ahmed Abdel Salam, Fadwa Fawzy, Norhan Shaker, Yasser M.Kadah Biomedical Engineering, Cairo University, Cairo, Egypt, ymk@k-space.org Computer
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationSparse MRI: The Application of Compressed Sensing for Rapid MR Imaging
Magnetic Resonance in Medicine 58:1182 1195 (2007) Sparse MRI: The Application of Compressed Sensing for Rapid MR Imaging Michael Lustig, 1 David Donoho, 2 and John M. Pauly 1 The sparsity which is implicit
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 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 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 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 informationEE369C: Assignment 6
EE369C Fall 2017-18 Medical Image Reconstruction 1 EE369C: Assignment 6 Due Wednesday Nov 15 In this assignment we will explore some of the basic elements of Compressed Sensing: Sparsity, Incoherency and
More informationIncoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Incoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY The separation of signal and noise is a key issue in seismic data processing. By
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 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 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 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 informationCompressive Sensing Based Image Reconstruction using Wavelet Transform
Compressive Sensing Based Image Reconstruction using Wavelet Transform Sherin C Abraham #1, Ketki Pathak *2, Jigna J Patel #3 # Electronics & Communication department, Gujarat Technological University
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 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 informationBlind Compressed Sensing Using Sparsifying Transforms
Blind Compressed Sensing Using Sparsifying Transforms Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and Coordinated Science Laboratory University of Illinois
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 informationWeighted-CS for reconstruction of highly under-sampled dynamic MRI sequences
Weighted- for reconstruction of highly under-sampled dynamic MRI sequences Dornoosh Zonoobi and Ashraf A. Kassim Dept. Electrical and Computer Engineering National University of Singapore, Singapore E-mail:
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 informationComputerLab: compressive sensing and application to MRI
Compressive Sensing, 207-8 ComputerLab: compressive sensing and application to MRI Aline Roumy This computer lab addresses the implementation and analysis of reconstruction algorithms for compressive sensing.
More informationMeasurements and Bits: Compressed Sensing meets Information Theory. Dror Baron ECE Department Rice University dsp.rice.edu/cs
Measurements and Bits: Compressed Sensing meets Information Theory Dror Baron ECE Department Rice University dsp.rice.edu/cs Sensing by Sampling Sample data at Nyquist rate Compress data using model (e.g.,
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 informationTERM PAPER ON The Compressive Sensing Based on Biorthogonal Wavelet Basis
TERM PAPER ON The Compressive Sensing Based on Biorthogonal Wavelet Basis Submitted By: Amrita Mishra 11104163 Manoj C 11104059 Under the Guidance of Dr. Sumana Gupta Professor Department of Electrical
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 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 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 informationImage Reconstruction from Multiple Sparse Representations
Image Reconstruction from Multiple Sparse Representations Robert Crandall Advisor: Professor Ali Bilgin University of Arizona Program in Applied Mathematics 617 N. Santa Rita, Tucson, AZ 85719 Abstract
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 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 informationCompressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction
Mathematical Modelling and Applications 2017; 2(6): 75-80 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170206.14 ISSN: 2575-1786 (Print); ISSN: 2575-1794 (Online) Compressed Sensing
More informationPERFORMANCE ANALYSIS BETWEEN TWO SPARSITY-CONSTRAINED MRI METHODS: HIGHLY CONSTRAINED BACKPROJECTION (HYPR) AND
PERFORMANCE ANALYSIS BETWEEN TWO SPARSITY-CONSTRAINED MRI METHODS: HIGHLY CONSTRAINED BACKPROJECTION (HYPR) AND COMPRESSED SENSING (CS) FOR DYNAMIC IMAGING A Thesis by NIBAL ARZOUNI Submitted to the Office
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 24 Compressed Sensing III M. Lustig, EECS UC Berkeley RADIOS https://inst.eecs.berkeley.edu/~ee123/ sp15/radio.html Interfaces and radios on Wednesday -- please
More informationLearning best wavelet packet bases for compressed sensing of classes of images: application to brain MR imaging
Learning best wavelet packet bases for compressed sensing of classes of images: application to brain MR imaging Michal P. Romaniuk, 1 Anil W. Rao, 1 Robin Wolz, 1 Joseph V. Hajnal 2 and Daniel Rueckert
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 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 informationCompressed Sensing MRI with Multichannel Data Using Multicore Processors
Compressed Sensing MRI with Multichannel Data Using Multicore Processors Ching-Hua Chang and Jim Ji* Magnetic Resonance in Medicine 64:1135 1139 (2010) Compressed sensing (CS) is a promising method to
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 informationVARIABLE DENSITY COMPRESSED IMAGE SAMPLING
VARIABLE DENSITY COMPRESSED IMAGE SAMPLING Zhongmin Wang, Gonzalo R. Arce and Jose L. Paredes Dep. of Electrical and Computer Engineering, University of Delaware, Newark, DE, 19716 USA. e-mail:{zhongmin,arce}@ee.udel.edu
More informationImproving the quality of compressed sensing MRI that exploits adjacent slice similarity
Improving the quality of compressed sensing MRI that exploits adjacent slice similarity Norihito Inamuro and Akira Hirabayashi College of Information Science and Engineering, Ritsumeikan University, Kusatsu,
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 informationLow-Rank and Sparse Matrix Decomposition for Accelerated Dynamic MRI with Separation of Background and Dynamic Components
Low-Rank and Sparse Matrix Decomposition for Accelerated Dynamic MRI with Separation of Background and Dynamic Components Ricardo Otazo 1, Emmanuel Candès 2, Daniel K. Sodickson 1 1 Department of Radiology,
More informationSairam Geethanath, Ph.D. Medical Imaging Research Centre Dayananda Sagar Institutions, Bangalore
Sairam Geethanath, Ph.D. Medical Imaging Research Centre Dayananda Sagar Institutions, Bangalore Contrast SNR MRI Speed Data provided by Baek Number of non-zero coefficients in a data vector Importance
More informationSignal Reconstruction from Sparse Representations: An Introdu. Sensing
Signal Reconstruction from Sparse Representations: An Introduction to Compressed Sensing December 18, 2009 Digital Data Acquisition Suppose we want to acquire some real world signal digitally. Applications
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More informationEE369C: Assignment 4 Solutions
EE369C Fall 2017-18 Medical Image Reconstruction 1 EE369C: Assignment 4 Solutions Due Wednesday, Oct. 25 The problems this week will be concerned with 2DFT SENSE reconstruction. The data is an axial brain
More informationMODEL-BASED FREE-BREATHING CARDIAC MRI RECONSTRUCTION USING DEEP LEARNED & STORM PRIORS: MODL-STORM
MODEL-BASED FREE-BREATHING CARDIAC MRI RECONSTRUCTION USING DEEP LEARNED & STORM PRIORS: MODL-STORM Sampurna Biswas, Hemant K. Aggarwal, Sunrita Poddar, and Mathews Jacob Department of Electrical and Computer
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 informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
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 informationA model-based method with joint sparsity constant for direct diffusion tensor estimation. Zhu, Y; Wu, Y; Zheng, YJ; Wu, EX; Ying, L; Liang, D
Title A model-based method with oint sparsity constant for direct diffusion tensor estimation Author(s) Zhu, Y; Wu, Y; Zheng, YJ; Wu, EX; Ying, L; Liang, D Citation The 9th IEEE International Symposium
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 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 informationCompressed Sensing MRI. [A look at how CS can improve on current imaging techniques] Digital Object Identifier /MSP.2007.
DIGITAL VISION Compressed Sensing MRI [A look at how CS can improve on current imaging techniques] [ Michael Lustig, David L. Donoho, Juan M. Santos, and John M. Pauly ] Compressed sensing (CS) aims to
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 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 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 informationMATHEMATICAL MODELING AND FAST NUMERICAL ALGORITHMS IN MEDICAL IMAGING APPLICATIONS
MATHEMATICAL MODELING AND FAST NUMERICAL ALGORITHMS IN MEDICAL IMAGING APPLICATIONS By XIAOJING YE A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF
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 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 informationCompressed sensing and undersampling k-space
Boston University OpenBU Theses & Dissertations http://open.bu.edu Boston University Theses & Dissertations 2015 Compressed sensing and undersampling k-space Rahimpour, Yashar https://hdl.handle.net/2144/16257
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationDeep Learning of Compressed Sensing Operators with Structural Similarity Loss
Deep Learning of Compressed Sensing Operators with Structural Similarity Loss Y. Zur and A. Adler Abstract Compressed sensing CS is a signal processing framework for efficiently reconstructing a signal
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 informationIntroduction. Wavelets, Curvelets [4], Surfacelets [5].
Introduction Signal reconstruction from the smallest possible Fourier measurements has been a key motivation in the compressed sensing (CS) research [1]. Accurate reconstruction from partial Fourier data
More informationSUMMARY. In combination with compressive sensing, a successful reconstruction
Higher dimensional blue-noise sampling schemes for curvelet-based seismic data recovery Gang Tang, Tsinghua University & UBC-Seismic Laboratory for Imaging and Modeling (UBC-SLIM), Reza Shahidi, UBC-SLIM,
More informationFast Dynamic MRI for Radiotherapy
1 Fast Dynamic MRI for Radiotherapy KE SHENG, PH.D., FAAPM DEPARTMENT OF RADIATION ONCOLOGY UNIVERSITY OF CALIFORNIA, LOS ANGELES Disclosures I receive research grants from Varian Medical Systems I am
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 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 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 informationComparison of MRI Under-Sampling Techniques for Compressed Sensing with Translation Invariant Wavelets Using FastTestCS: A Flexible Simulation Tool
Journal of Signal and Information Processing, 2016, 7, 252-271 http://www.scirp.org/journal/jsip ISSN Online: 2159-4481 ISSN Print: 2159-4465 Comparison of MRI Under-Sampling Techniques for Compressed
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 informationMotion Artifacts and Suppression in MRI At a Glance
Motion Artifacts and Suppression in MRI At a Glance Xiaodong Zhong, PhD MR R&D Collaborations Siemens Healthcare MRI Motion Artifacts and Suppression At a Glance Outline Background Physics Common Motion
More informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
More informationCompressed Sensing and Applications by using Dictionaries in Image Processing
Advances in Computational Sciences and Technology ISSN 0973-6107 Volume 10, Number 2 (2017) pp. 165-170 Research India Publications http://www.ripublication.com Compressed Sensing and Applications by using
More informationAn Iterative Approach for Reconstruction of Arbitrary Sparsely Sampled Magnetic Resonance Images
An Iterative Approach for Reconstruction of Arbitrary Sparsely Sampled Magnetic Resonance Images Hamed Pirsiavash¹, Mohammad Soleymani², Gholam-Ali Hossein-Zadeh³ ¹Department of electrical engineering,
More informationELEG Compressive Sensing and Sparse Signal Representations
ELEG 867 - Compressive Sensing and Sparse Signal Representations Gonzalo R. Arce Depart. of Electrical and Computer Engineering University of Delaware Fall 211 Compressive Sensing G. Arce Fall, 211 1 /
More informationMODEL-based recovery of images from noisy and sparse. MoDL: Model Based Deep Learning Architecture for Inverse Problems
1 MoDL: Model Based Deep Learning Architecture for Inverse Problems Hemant K. Aggarwal, Member, IEEE, Merry P. Mani, and Mathews Jacob, Senior Member, IEEE arxiv:1712.02862v3 [cs.cv] 10 Aug 2018 Abstract
More informationCardiac diffusion tensor imaging based on compressed sensing using joint sparsity and low-rank approximation
Technology and Health Care 24 2016) S593 S599 DOI 10.3233/THC-161186 IOS Press S593 Cardiac diffusion tensor imaging based on compressed sensing using joint sparsity and low-rank approximation Jianping
More informationAutomated Aperture GenerationQuantitative Evaluation of Ape
. Objectives Approaches Randomness Reconstructing Images from K-space References Automated Aperture Generation Quantitative Evaluation of Aperture April 5, 011 . Objectives Approaches Randomness Reconstructing
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
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 informationP-LORAKS: Low-Rank Modeling of Local k-space Neighborhoods with Parallel Imaging Data
P-LORAKS: Low-Rank Modeling of Local k-space Neighborhoods with Parallel Imaging Data Justin P. Haldar 1, Jingwei Zhuo 2 1 Electrical Engineering, University of Southern California, Los Angeles, CA, USA
More informationModule 5: Dynamic Imaging and Phase Sharing. (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review. Improving Temporal Resolution.
MRES 7005 - Fast Imaging Techniques Module 5: Dynamic Imaging and Phase Sharing (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review Improving Temporal Resolution True-FISP (I) True-FISP (II) Keyhole
More informationCompressive Sensing MRI with Wavelet Tree Sparsity
Compressive Sensing MRI with Wavelet Tree Sparsity Chen Chen and Junzhou Huang Department of Computer Science and Engineering University of Texas at Arlington cchen@mavs.uta.edu jzhuang@uta.edu Abstract
More informationCompressive Sensing Applications and Demonstrations: Synthetic Aperture Radar
Compressive Sensing Applications and Demonstrations: Synthetic Aperture Radar Shaun I. Kelly The University of Edinburgh 1 Outline 1 SAR Basics 2 Compressed Sensing SAR 3 Other Applications of Sparsity
More informationClustered Compressive Sensing: Application on Medical Imaging
Clustered Compressive Sensing: Application on Medical Imaging Solomon A. Tesfamicael, IACSIT Member and Faraz Barzideh Abstract This paper provides clustered compressive sensing (CCS) based image processing
More informationMain Menu. Summary. sampled) f has a sparse representation transform domain S with. in certain. f S x, the relation becomes
Preliminary study on Dreamlet based compressive sensing data recovery Ru-Shan Wu*, Yu Geng 1 and Lingling Ye, Modeling and Imaging Lab, Earth & Planetary Sciences/IGPP, University of California, Santa
More informationImage Denoising Based on Hybrid Fourier and Neighborhood Wavelet Coefficients Jun Cheng, Songli Lei
Image Denoising Based on Hybrid Fourier and Neighborhood Wavelet Coefficients Jun Cheng, Songli Lei College of Physical and Information Science, Hunan Normal University, Changsha, China Hunan Art Professional
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 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 informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationCompressive Sensing Algorithms for Fast and Accurate Imaging
Compressive Sensing Algorithms for Fast and Accurate Imaging Wotao Yin Department of Computational and Applied Mathematics, Rice University SCIMM 10 ASU, Tempe, AZ Acknowledgements: results come in part
More informationImaging and Deconvolution
Imaging and Deconvolution Urvashi Rau National Radio Astronomy Observatory, Socorro, NM, USA The van-cittert Zernike theorem Ei E V ij u, v = I l, m e sky j 2 i ul vm dldm 2D Fourier transform : Image
More information