Blind Compressed Sensing Using Sparsifying Transforms
|
|
- Virgil Leonard
- 5 years ago
- Views:
Transcription
1 Blind Compressed Sensing Using Sparsifying Transforms Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and Coordinated Science Laboratory University of Illinois at Urbana-Champaign May 29,
2 Key Topics of Talk Non-adaptive Compressed Sensing (CS) Synthesis dictionary learning-based blind compressed sensing Transform learning vs. Dictionary learning Transform learning-based blind compressed sensing Application to magnetic resonance imaging (MRI) Transform learning based MRI (TLMRI) Conclusions 2
3 Compressed Sensing (CS) CS enables accurate recovery of images from far fewer measurements than the number of unknowns Sparsity of image in transform domain or dictionary measurement procedure incoherent with transform Reconstruction non-linear and expensive Reconstruction problem (NP-hard) - min x Data Fidelity Regularizer {}}{{}}{ Ax y 2 2 +λ x C P : image as vector, y C m : measurements. Ψx 0 (1) A C m P : sensing matrix (m < P), Ψ : transform (Wavelets, Contourlets, Total Variation). l 0 norm counts non-zeros. Iterative algorithms for CS reconstruction are usually expensive. 3
4 Application: Compressed Sensing MRI (CSMRI) Data - samples in k-space of spatial Fourier transform of object, acquired sequentially. Acquisition rate limited by MR physics, physiological constraints on RF energy deposition. CS accelerates the data acquisition in MRI. CSMRI with non-adaptive transforms or dictionaries limited to fold undersampling [Ma et al. 08]. Two directions to improve CSMRI - better or adaptive sparse modeling better choice of sampling pattern (F u) [EMBC, 2011] Fig. from Lustig et al. 07 4
5 Synthesis Model for Sparse Representation Given a signal y R n, and dictionary D R n K, we assume y = Dx with x 0 K a union of subspaces model. Real world signals modeled as y = Dx +e, e is deviation term. Given D, sparsity level s, the synthesis sparse coding problem is This problem is NP-hard. ˆx = argmin x y Dx 2 2 s.t. x 0 s 5
6 Synthesis Dictionary Learning The DL problem (NP-hard) - min D,B N R j x Db j 2 2 s.t. d k 2 = 1 k, b j 0 s j (2) R j x C n - n n patch indexed by location in image. R j C n P extracts patch. D C n K - patch based dictionary. b j C K - sparse, x j Db j. s - sparsity, B = [b 1 b 2... b N ]. DL minimizes fit error of all patches using sparse representations w.r.t. D. 6
7 Synthesis-based Blind Compressed Sensing (BCS) (P0) min x,d,b Sparse Fitting Regularizer {}}{ N R j x Db j 2 2 +ν Data Fidelity {}}{ Ax y 2 2 s.t. d k 2 = 1 k, b j 0 s j. B C n N : matrix that has the sparse codes b j as its columns. (P0) learns D C n K, and reconstructs x, from only undersampled y dictionary adaptive to underlying image. (P0) is NP-hard, non-convex even if l 0 norm relaxed to l 1. DLMRI 1 solves (P0) for MRI and works better than non-adaptive CS. Synthesis BCS algorithms have no guarantees and are expensive. 1 [Ravishankar & Bresler 11] 7
8 2D Random Sampling - 6 fold undersampling LDP 2 reconstruction (22 db) LDP error magnitude DLMRI reconstruction (32 db) DLMRI error magnitude 0 MRI data from Miki Lustig. 2 [Lustig et al. 07] 8
9 Alternative: Sparsifying Transform Model Given a signal y R n, and transform W R m n, we model Wy = x +η with x 0 m and η - error term. Natural signals are approximately sparse in Wavelets, DCT. Given W, and sparsity s, transform sparse coding is ˆx = argmin x Wy x 2 2 s.t. x 0 s ˆx = H s(wy) computed by thresholding Wy to the s largest magnitude elements. Sparse coding is cheap. Signal recovered as W ˆx. Sparsifying transforms exploited for compression (JPEG2000), etc. 9
10 Alternative: Sparsifying Transform Learning Square Transform Models Unstructured transform learning [IEEE TSP, 2013 & 2015] Doubly sparse transform learning [IEEE TIP, 2013] Online learning for Big Data [IEEE JSTSP, 2015] Convex formulations for transform learning [ICASSP, 2014] Overcomplete Transform Models Unstructured overcomplete transform learning [ICASSP, 2013] Learning structured overcomplete transforms with block cosparsity (OCTOBOS) [IJCV, 2014] Applications: Sparse representation, Image & Video denoising, Classification, Blind compressed sensing (BCS) for imaging. 10
11 Square Transform Learning Formulation (P1) min W,B Sparsification Error {}}{ N WR j x b j 2 2 +λ s.t. b j 0 s j Regularizer {}}{ (0.5 W 2F log detw ) Sparsification error - measures deviation of data in transform domain from perfect sparsity. Regularizer enables complete control over conditioning & scaling of W. If (Ŵ, ˆB) such that the condition number κ(ŵ) = 1, ŴR j x = ˆb j, ˆb j 0 s j globally identifiable by solving (P1). (P1) favors both a low sparsification error and good conditioning. The solution to (P1) is unitary as λ. 11
12 Transform-based Blind Compressed Sensing (BCS) (P2) min x,w,b s.t. Sparsification Error {}}{ N Data Fidelity Regularizer {}}{{}}{ WR j x b j 2 2 +ν Ax y 2 2 +λ v(w) N b j 0 s, x 2 C. (P2) learns W C n n, and reconstructs x, from only undersampled y transform adaptive to underlying image. v(w) log detw +0.5 W 2 F controls scaling and κ of W. x 2 C is an energy/range constraint. C > 0. 12
13 Transform BCS: Identifiability & Uniqueness Proposition 1 Let x C p, and let y = Ax with A C m p. Suppose x 2 C W C n n is a unitary transform N WR jx 0 s Further, let B denote the matrix that has WR j x as its columns. Then, (x,w,b) is a global minimizer of Problem (P2), i.e., it is identifiable by solving (P2). Given minimizer (x,w,b) of (P2), (x,θw,θb) is another equivalent minimizer Θ s.t. Θ H Θ = I, j Θb j 0 s. The optimal x is invariant to such transformations of (W,B). When W is constrained to be doubly sparse and unitary, uniqueness can be guaranteed under additional (e.g., spark) conditions. 13
14 Alternative Transform BCS Formulations (P3) min x,w,b N WR j x b j 2 2 +ν Ax y 2 2 s.t. W H W = I, N b j 0 s, x 2 C. (P3) is also a unitary synthesis dictionary-based BCS problem, with W H the synthesis dictionary. (P4) min x,w,b N WR j x b j 2 2 +ν Ax y 2 2 +λv(w)+ N η2 b j 0 s.t. x 2 C. 14
15 Block Coordinate Descent (BCD) Algorithm for (P2) (P2) solved by alternating between updating W, B, and x. Alternate a few times between the W and B updates, before performing an image update. Sparse Coding Step solves (P2) for B with fixed x, W. min B N WR j x b j 2 2 s.t. N b j 0 s. (3) Cheap Solution: Let Z C n N be the matrix with WR jx as its columns. Solution ˆB = H s(z) computed exactly by zeroing out all but the s largest magnitude coefficients in Z. 15
16 BCD Algorithm for (P2) Transform Update Step solves (P2) for W with fixed x, B. min W N WR j x b j λ W 2 F λlog detw (4) Let X C n N be the matrix with R jx as its columns. Closed-form solution: Ŵ = 0.5R ( ( ) Σ+ Σ 2 )1 2 +2λI V H L 1 (5) where XX H +0.5λI = LL H, and L 1 XB H has a full SVD of VΣR H. Solution is unique if and only if XB H is non-singular. 16
17 BCD Algorithm for (P2) Image Update Step solves (P2) for x with fixed W, B. min x N WR j x b j 2 2 +ν Ax y 2 2 s.t. x 2 C. (6) Least squares problem with l 2 norm constraint. Solution is unique as long as the set of overlapping patches cover all image pixels. Solve Least squares Lagrangian formulation: min x N WR jx b j 2 2 +ν Ax y 2 2 +µ ( x 2 2 C ) The optimal multiplier ˆµ R + is the smallest real such that ˆx 2 C. ˆµ and ˆx can be found cheaply. (7) 17
18 BCS Convergence Guarantees - Notations Define the barrier function ψ s (B) as { N 0, ψ s (B) = b j 0 s +, else χ C (x) is the barrier function corresponding to x 2 C. (P2) is equivalent to the problem of minimizing the objective g(w,b,x) = N WR jx b j 2 2 +ν Ax y 2 2 +λv(w)+ψ(b)+χ(x) For H C p q, ρ j (H) is the magnitude of the j th largest element (magnitude-wise) of H. X C n N denotes a matrix with R j x, 1 j N, as its columns. 18
19 Transform BCS Convergence Guarantees Theorem 1 For the sequence {W t,b t,x t } generated by the BCD Algorithm with initial (W 0,B 0,x 0 ), we have {g (W t,b t,x t )} g = g (W 0,B 0,x 0 ). {W t,b t,x t } is bounded, and all its accumulation points are equivalent, i.e., they achieve the same value g of the objective. x t x t as t. Every accumulation point (W,B,x) is a critical point of g satisfying the following partial global optimality conditions W argmin W x argmin g (W,B, x) (8) x ( ) ( ) g W,B,x, B argmin g W, B,x (9) B 19
20 Transform BCS Convergence Guarantees Theorem 2 Each accumulation point (W,B,x) of {W t,b t,x t } also satisfies the following partial local optimality conditions g(w + W,B + B,x) g(w,b,x) = g (10) g(w,b + B,x + x) g(w,b,x) = g (11) The conditions each hold for all x C p, and all W C n n satisfying W F ǫ for some ǫ = ǫ(w) > 0, and all B C n N in R1 R2 R1. The half-space Re ( tr { (WX B) B H}) 0. R2. The local region defined by B < ρ s (WX). Furthermore, if WX 0 s, then B can be arbitrary. 20
21 Global Convergence Guarantees Proposition 2 For each initialization, the iterate sequence in the BCD algorithm converges to an equivalence class (same objective values) of critical points of the objective, that are also partial global/local minimizers. Proposition 3 The BCD algorithm is globally convergent to a subset of the set of critical points of the objective. The subset includes all (W,B,x) that are at least partial global and partial local minimizers. 21
22 Computational Advantages of Transform BCS Cost per iteration of transform BCS: O(p 4 NL) N overlapping patches of size p p, W C n n, n p 2. L : # inner alternations between transform update & sparse coding. Cost per iteration of Synthesis BCS method DLMRI 3 : O(p 6 NJ). D C n K, n p 2, K n, sparsity s n. J : # of inner iterations of dictionary learning using K-SVD 4. In practice, transform BCS converges quickly and is much cheaper for large p. In 3D or 4D imaging, n = p 3 or p 4, and the gain in computations is about a factor n in order. 3 [Ravishankar & Bresler 11] 4 [Aharon et al. 06], 22
23 TLMRI Convergence - 4x Undersampling (s = 3.4%) 9.7 x 105 Reference Sampling mask 10 1 Objective Function x t x t Iteration Number Objective Iteration Number (t) x t x t 1 2 vs. t 23
24 Convergence & Learning - 4x Undersampling (s = 3.4%) Zero-filling (28.94 db) Zero-filling Error 0 TLMRI (32.66 db) real (top), imaginary (bottom) parts of learnt W 24
25 Comparison (PSNR & Runtime) to Recent Methods Sampling Scheme Undersampling Zero-filling LDP 5 PBDWS 6 DLMRI 7 PANO 8 TLMRI 2D Random 4x x Cartesian 4x x Avg. Runtime (s) TLMRI is up to 5.5 db better than LDP, that uses Wavelets + TV. TLMRI provides up to 1 db improvement in PSNR over the PBDWS method that uses redundant Wavelets and trained patch-based geometric directions, and is up to 1.6 db better than the non-local PANO method. It is up to 0.35 db better than DLMRI, that learns 4x overcomplete dictionary. TLMRI is 10x faster than DLMRI, and 4x faster than the PBDWS method. TLMRI provides the best reconstructions, and is the fastest. 5 [Lustig et al. 07] 6 [Ning et al. 13] 7 [Ravishankar & Bresler 11] 8 [Qu et al. 14] 25
26 Example - 2D random 5x Undersampling Reference DLMRI (28.54 db) TLMRI (30.47 db) Sampling Mask DLMRI Error TLMRI Error
27 Conclusions We introduced a transform-based BCS framework Proposed BCS algorithms have a low computational cost. We provided novel convergence guarantees for the algorithms, that do not require any restrictive assumptions. For CSMRI, the proposed approach is better than leading image reconstruction methods, while being much faster. Future work: convergence of algorithm to global minimizer & convergence rate. 27
28 Thank you! Questions?? 28
29 Convergence Guarantees - Definitions Definition 1 Let φ : R q (,+ ] be a proper function and let z domφ. The Fréchet sub-differential of the function φ at z is the following set: ˆ φ(z) { h R q : liminf b z,b z 1 b z (φ(b) φ(z) b z,h ) 0} (12) If z / domφ, then ˆ φ(z) =. The sub-differential of φ at z is defined as φ(z) { h R q : z k z,φ(z k ) φ(z),h k ˆ φ(z k ) h }. (13) Lemma 1 A necessary condition for z R q to be a minimizer of the function φ : R q (,+ ] is that z is a critical point of φ, i.e., 0 φ(z). If φ is a convex function, this condition is also sufficient. 29
Blind Compressed Sensing Using Sparsifying Transforms
Blind Compressed Sensing Using Sparsifying Transforms Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and the Coordinated Science Laboratory University of Illinois
More informationarxiv: v2 [cs.lg] 23 Oct 2015
SIAM J. IMAGING SCIENCES Vol. xx, pp. x c 2014 Society for Industrial and Applied Mathematics x x Efficient Blind Compressed Sensing Using Sparsifying Transforms with Convergence Guarantees and Application
More informationWhen Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration
When Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration Bihan Wen, Yanjun Li and Yoram Bresler Department of Electrical and Computer Engineering Coordinated
More informationData-Driven Learning of a Union of Sparsifying Transforms Model for Blind Compressed Sensing
1 Data-Driven Learning of a Union of Sparsifying Transforms Model for Blind Compressed Sensing Saiprasad Ravishankar, Member, IEEE, Yoram Bresler, Fellow, IEEE Abstract Compressed sensing is a powerful
More informationarxiv: v4 [cs.lg] 16 Oct 2017
FRIST - Flipping and Rotation Invariant Sparsifying Transform Learning and Applications arxiv:1511.06359v4 [cs.lg] 16 Oct 2017 Bihan Wen 1, Saiprasad Ravishankar 2, and Yoram Bresler 1 1 Department of
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 informationLEARNING OVERCOMPLETE SPARSIFYING TRANSFORMS WITH BLOCK COSPARSITY
LEARNING OVERCOMPLETE SPARSIFYING TRANSFORMS WITH BLOCK COSPARSITY Bihan Wen, Saiprasad Ravishankar, and Yoram Bresler Department of Electrical and Computer Engineering and the Coordinated Science Laboratory,
More informationLEARNING FLIPPING AND ROTATION INVARIANT SPARSIFYING TRANSFORMS
LEARNING FLIPPING AND ROTATION INVARIANT SPARSIFYING TRANSFORMS Bihan Wen 1, Saiprasad Ravishankar, and Yoram Bresler 1 1 Department of Electrical and Computer Engineering and the Coordinated Science Laboratory
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 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 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 informationUnion of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction
Union of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction Xuehang Zheng 1, Saiprasad Ravishankar 2, Yong Long 1, Jeff Fessler 2 1 University of Michigan - Shanghai Jiao Tong University
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 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 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 informationCOMPRESSED SENSING MRI WITH BAYESIAN DICTIONARY LEARNING.
COMPRESSED SENSING MRI WITH BAYESIAN DICTIONARY LEARNING Xinghao Ding 1, John Paisley 2, Yue Huang 1, Xianbo Chen 1, Feng Huang 3 and Xiao-ping Zhang 1,4 1 Department of Communications Engineering, Xiamen
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 informationLow-Rank and Adaptive Sparse Signal (LASSI) Models for Highly Accelerated Dynamic Imaging
1116 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 36, NO. 5, MAY 2017 Low-Rank and Adaptive Sparse Signal (LASSI) Models for Highly Accelerated Dynamic Imaging Saiprasad Ravishankar, Member, IEEE, Brian
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 informationECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis
ECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis Yuejie Chi Departments of ECE and BMI The Ohio State University September 24, 2015 Time, location, and office hours Time: Tue/Thu
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 informationInverse problem regularization using adaptive signal models
Inverse problem regularization using adaptive signal models Jeffrey A. Fessler William L. Root Professor of EECS EECS Dept., BME Dept., Dept. of Radiology University of Michigan http://web.eecs.umich.edu/
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 informationBayesian Nonparametric Dictionary Learning for Compressed Sensing MRI
1 Bayesian Nonparametric Dictionary Learning for Compressed Sensing MRI Yue Huang, John Paisley, Qin Lin, Xinghao Ding, Xueyang Fu and Xiao-ping Zhang, Senior Member, IEEE arxiv:12.2712v3 [cs.cv] 26 Jul
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 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 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 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 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 informationarxiv: v1 [cs.cv] 12 Feb 2013
1 MR Image Reconstruction from Undersampled k-space with Bayesian Dictionary Learning Yue Huang 1, John Paisley 2, Xianbo Chen 1, Xinghao Ding 1, Feng Huang 3 and Xiao-ping Zhang 4 arxiv:12.2712v1 [cs.cv]
More informationDetecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference
Detecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference Minh Dao 1, Xiang Xiang 1, Bulent Ayhan 2, Chiman Kwan 2, Trac D. Tran 1 Johns Hopkins Univeristy, 3400
More informationDEEP LEARNING OF COMPRESSED SENSING OPERATORS WITH STRUCTURAL SIMILARITY (SSIM) LOSS
DEEP LEARNING OF COMPRESSED SENSING OPERATORS WITH STRUCTURAL SIMILARITY (SSIM) LOSS ABSTRACT Compressed sensing (CS) is a signal processing framework for efficiently reconstructing a signal from a small
More informationMAGNETIC resonance imaging (MRI) is a widely used
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 23, NO. 12, DECEMBER 2014 5007 Bayesian Nonparametric Dictionary Learning for Compressed Sensing MRI Yue Huang, John Paisley, Qin Lin, Xinghao Ding, Xueyang
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 informationOutline Introduction Problem Formulation Proposed Solution Applications Conclusion. Compressed Sensing. David L Donoho Presented by: Nitesh Shroff
Compressed Sensing David L Donoho Presented by: Nitesh Shroff University of Maryland Outline 1 Introduction Compressed Sensing 2 Problem Formulation Sparse Signal Problem Statement 3 Proposed Solution
More information1. Introduction. performance of numerical methods. complexity bounds. structural convex optimization. course goals and topics
1. Introduction EE 546, Univ of Washington, Spring 2016 performance of numerical methods complexity bounds structural convex optimization course goals and topics 1 1 Some course info Welcome to EE 546!
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 informationSection 5 Convex Optimisation 1. W. Dai (IC) EE4.66 Data Proc. Convex Optimisation page 5-1
Section 5 Convex Optimisation 1 W. Dai (IC) EE4.66 Data Proc. Convex Optimisation 1 2018 page 5-1 Convex Combination Denition 5.1 A convex combination is a linear combination of points where all coecients
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 informationImage Restoration with Deep Generative Models
Image Restoration with Deep Generative Models Raymond A. Yeh *, Teck-Yian Lim *, Chen Chen, Alexander G. Schwing, Mark Hasegawa-Johnson, Minh N. Do Department of Electrical and Computer Engineering, University
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 informationConvex Optimization and Machine Learning
Convex Optimization and Machine Learning Mengliu Zhao Machine Learning Reading Group School of Computing Science Simon Fraser University March 12, 2014 Mengliu Zhao SFU-MLRG March 12, 2014 1 / 25 Introduction
More informationBilevel Sparse Coding
Adobe Research 345 Park Ave, San Jose, CA Mar 15, 2013 Outline 1 2 The learning model The learning algorithm 3 4 Sparse Modeling Many types of sensory data, e.g., images and audio, are in high-dimensional
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 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 informationA primal-dual framework for mixtures of regularizers
A primal-dual framework for mixtures of regularizers Baran Gözcü baran.goezcue@epfl.ch Laboratory for Information and Inference Systems (LIONS) École Polytechnique Fédérale de Lausanne (EPFL) Switzerland
More informationConvex or non-convex: which is better?
Sparsity Amplified Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering New York University Brooklyn, New York March 7 / 4 Convex or non-convex: which is better? (for sparse-regularized
More informationImage Deblurring Using Adaptive Sparse Domain Selection and Adaptive Regularization
Volume 3, No. 3, May-June 2012 International Journal of Advanced Research in Computer Science RESEARCH PAPER Available Online at www.ijarcs.info ISSN No. 0976-5697 Image Deblurring Using Adaptive Sparse
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 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 informationUniqueness in Bilinear Inverse Problems with Applications to Subspace and Sparsity Models
Uniqueness in Bilinear Inverse Problems with Applications to Subspace and Sparsity Models Yanjun Li Joint work with Kiryung Lee and Yoram Bresler Coordinated Science Laboratory Department of Electrical
More informationImage Denoising Using Sparse Representations
Image Denoising Using Sparse Representations SeyyedMajid Valiollahzadeh 1,,HamedFirouzi 1, Massoud Babaie-Zadeh 1, and Christian Jutten 2 1 Department of Electrical Engineering, Sharif University of Technology,
More informationImage Restoration Using DNN
Image Restoration Using DNN Hila Levi & Eran Amar Images were taken from: http://people.tuebingen.mpg.de/burger/neural_denoising/ Agenda Domain Expertise vs. End-to-End optimization Image Denoising and
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
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 informationANALYSIS SPARSE CODING MODELS FOR IMAGE-BASED CLASSIFICATION. Sumit Shekhar, Vishal M. Patel and Rama Chellappa
ANALYSIS SPARSE CODING MODELS FOR IMAGE-BASED CLASSIFICATION Sumit Shekhar, Vishal M. Patel and Rama Chellappa Center for Automation Research, University of Maryland, College Park, MD 20742 {sshekha, pvishalm,
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 informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More informationTwo-Parameter Selection Techniques for Projection-based Regularization Methods: Application to Partial-Fourier pmri
Two-Parameter Selection Techniques for Projection-based Regularization Methods: Application to Partial-Fourier pmri Misha E. Kilmer 1 Scott Hoge 1 Dept. of Mathematics, Tufts University Medford, MA Dept.
More informationConvex Optimization MLSS 2015
Convex Optimization MLSS 2015 Constantine Caramanis The University of Texas at Austin The Optimization Problem minimize : f (x) subject to : x X. The Optimization Problem minimize : f (x) subject to :
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 Splines for Sparse Tomographic Reconstruction. Elham Sakhaee and Alireza Entezari University of Florida
Learning Splines for Sparse Tomographic Reconstruction Elham Sakhaee and Alireza Entezari University of Florida esakhaee@cise.ufl.edu 2 Tomographic Reconstruction Recover the image given X-ray measurements
More informationCompressive Sensing for Multimedia. Communications in Wireless Sensor Networks
Compressive Sensing for Multimedia 1 Communications in Wireless Sensor Networks Wael Barakat & Rabih Saliba MDDSP Project Final Report Prof. Brian L. Evans May 9, 2008 Abstract Compressive Sensing is an
More informationCOMPRESSIVE VIDEO SAMPLING
COMPRESSIVE VIDEO SAMPLING Vladimir Stanković and Lina Stanković Dept of Electronic and Electrical Engineering University of Strathclyde, Glasgow, UK phone: +44-141-548-2679 email: {vladimir,lina}.stankovic@eee.strath.ac.uk
More informationConvex Optimization / Homework 2, due Oct 3
Convex Optimization 0-725/36-725 Homework 2, due Oct 3 Instructions: You must complete Problems 3 and either Problem 4 or Problem 5 (your choice between the two) When you submit the homework, upload a
More informationTomographic reconstruction: the challenge of dark information. S. Roux
Tomographic reconstruction: the challenge of dark information S. Roux Meeting on Tomography and Applications, Politecnico di Milano, 20-22 April, 2015 Tomography A mature technique, providing an outstanding
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More informationOff-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples
Off-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples Greg Ongie PhD Candidate Department of Applied Math and Computational Sciences University of Iowa April
More informationAdaptive Reconstruction Methods for Low-Dose Computed Tomography
Adaptive Reconstruction Methods for Low-Dose Computed Tomography Joseph Shtok Ph.D. supervisors: Prof. Michael Elad, Dr. Michael Zibulevsky. Technion IIT, Israel, 011 Ph.D. Talk, Apr. 01 Contents of this
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 informationSupport Vector Machines.
Support Vector Machines srihari@buffalo.edu SVM Discussion Overview 1. Overview of SVMs 2. Margin Geometry 3. SVM Optimization 4. Overlapping Distributions 5. Relationship to Logistic Regression 6. Dealing
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 informationBSIK-SVD: A DICTIONARY-LEARNING ALGORITHM FOR BLOCK-SPARSE REPRESENTATIONS. Yongqin Zhang, Jiaying Liu, Mading Li, Zongming Guo
2014 IEEE International Conference on Acoustic, Speech and Signal Processing (ICASSP) BSIK-SVD: A DICTIONARY-LEARNING ALGORITHM FOR BLOCK-SPARSE REPRESENTATIONS Yongqin Zhang, Jiaying Liu, Mading Li, Zongming
More informationEE 584 MACHINE VISION
EE 584 MACHINE VISION Binary Images Analysis Geometrical & Topological Properties Connectedness Binary Algorithms Morphology Binary Images Binary (two-valued; black/white) images gives better efficiency
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 informationLearning a Spatiotemporal Dictionary for Magnetic Resonance Fingerprinting with Compressed Sensing
Learning a Spatiotemporal Dictionary for Magnetic Resonance Fingerprinting with Compressed Sensing Pedro A Gómez 1,2,3, Cagdas Ulas 1, Jonathan I Sperl 3, Tim Sprenger 1,2,3, Miguel Molina-Romero 1,2,3,
More informationPatch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques
Patch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques Syed Gilani Pasha Assistant Professor, Dept. of ECE, School of Engineering, Central University of Karnataka, Gulbarga,
More informationNon-Local Compressive Sampling Recovery
Non-Local Compressive Sampling Recovery Xianbiao Shu 1, Jianchao Yang 2 and Narendra Ahuja 1 1 University of Illinois at Urbana-Champaign, 2 Adobe Research 1 {xshu2, n-ahuja}@illinois.edu, 2 jiayang@adobe.com
More informationDALM-SVD: ACCELERATED SPARSE CODING THROUGH SINGULAR VALUE DECOMPOSITION OF THE DICTIONARY
DALM-SVD: ACCELERATED SPARSE CODING THROUGH SINGULAR VALUE DECOMPOSITION OF THE DICTIONARY Hugo Gonçalves 1, Miguel Correia Xin Li 1 Aswin Sankaranarayanan 1 Vitor Tavares Carnegie Mellon University 1
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 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 informationNon-Parametric Bayesian Dictionary Learning for Sparse Image Representations
Non-Parametric Bayesian Dictionary Learning for Sparse Image Representations Mingyuan Zhou, Haojun Chen, John Paisley, Lu Ren, 1 Guillermo Sapiro and Lawrence Carin Department of Electrical and Computer
More informationImage Restoration and Background Separation Using Sparse Representation Framework
Image Restoration and Background Separation Using Sparse Representation Framework Liu, Shikun Abstract In this paper, we introduce patch-based PCA denoising and k-svd dictionary learning method for the
More informationImage Inpainting Using Sparsity of the Transform Domain
Image Inpainting Using Sparsity of the Transform Domain H. Hosseini*, N.B. Marvasti, Student Member, IEEE, F. Marvasti, Senior Member, IEEE Advanced Communication Research Institute (ACRI) Department of
More informationRobust image recovery via total-variation minimization
Robust image recovery via total-variation minimization Rachel Ward University of Texas at Austin (Joint work with Deanna Needell, Claremont McKenna College) February 16, 2012 2 Images are compressible
More informationCombinatorial Selection and Least Absolute Shrinkage via The CLASH Operator
Combinatorial Selection and Least Absolute Shrinkage via The CLASH Operator Volkan Cevher Laboratory for Information and Inference Systems LIONS / EPFL http://lions.epfl.ch & Idiap Research Institute joint
More informationA Patch Prior for Dense 3D Reconstruction in Man-Made Environments
A Patch Prior for Dense 3D Reconstruction in Man-Made Environments Christian Häne 1, Christopher Zach 2, Bernhard Zeisl 1, Marc Pollefeys 1 1 ETH Zürich 2 MSR Cambridge October 14, 2012 A Patch Prior for
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 informationA fast algorithm for sparse reconstruction based on shrinkage, subspace optimization and continuation [Wen,Yin,Goldfarb,Zhang 2009]
A fast algorithm for sparse reconstruction based on shrinkage, subspace optimization and continuation [Wen,Yin,Goldfarb,Zhang 2009] Yongjia Song University of Wisconsin-Madison April 22, 2010 Yongjia Song
More informationGEOMETRIC MANIFOLD APPROXIMATION USING LOCALLY LINEAR APPROXIMATIONS
GEOMETRIC MANIFOLD APPROXIMATION USING LOCALLY LINEAR APPROXIMATIONS BY TALAL AHMED A thesis submitted to the Graduate School New Brunswick Rutgers, The State University of New Jersey in partial fulfillment
More informationFace Recognition via Sparse Representation
Face Recognition via Sparse Representation John Wright, Allen Y. Yang, Arvind, S. Shankar Sastry and Yi Ma IEEE Trans. PAMI, March 2008 Research About Face Face Detection Face Alignment Face Recognition
More informationBlind compressive sensing dynamic MRI
1 Blind compressive sensing dynamic MRI Sajan Goud Lingala, Student Member, IEEE, Mathews Jacob, Senior Member, IEEE Abstract We propose a novel blind compressive sensing (BCS) frame work to recover dynamic
More informationAn Optimized Pixel-Wise Weighting Approach For Patch-Based Image Denoising
An Optimized Pixel-Wise Weighting Approach For Patch-Based Image Denoising Dr. B. R.VIKRAM M.E.,Ph.D.,MIEEE.,LMISTE, Principal of Vijay Rural Engineering College, NIZAMABAD ( Dt.) G. Chaitanya M.Tech,
More informationNon-Differentiable Image Manifolds
The Multiscale Structure of Non-Differentiable Image Manifolds Michael Wakin Electrical l Engineering i Colorado School of Mines Joint work with Richard Baraniuk, Hyeokho Choi, David Donoho Models for
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 informationIMA Preprint Series # 2211
LEARNING TO SENSE SPARSE SIGNALS: SIMULTANEOUS SENSING MATRIX AND SPARSIFYING DICTIONARY OPTIMIZATION By Julio Martin Duarte-Carvajalino and Guillermo Sapiro IMA Preprint Series # 2211 ( May 2008 ) INSTITUTE
More informationCS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 36
CS 473: Algorithms Ruta Mehta University of Illinois, Urbana-Champaign Spring 2018 Ruta (UIUC) CS473 1 Spring 2018 1 / 36 CS 473: Algorithms, Spring 2018 LP Duality Lecture 20 April 3, 2018 Some of the
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 informationTRADITIONAL patch-based sparse coding has been
Bridge the Gap Between Group Sparse Coding and Rank Minimization via Adaptive Dictionary Learning Zhiyuan Zha, Xin Yuan, Senior Member, IEEE arxiv:709.03979v2 [cs.cv] 8 Nov 207 Abstract Both sparse coding
More informationLecture 17 Sparse Convex Optimization
Lecture 17 Sparse Convex Optimization Compressed sensing A short introduction to Compressed Sensing An imaging perspective 10 Mega Pixels Scene Image compression Picture Why do we compress images? Introduction
More information