Breaking the coherence barrier - A new theory for compressed sensing
|
|
- Albert Robbins
- 6 years ago
- Views:
Transcription
1 Breaking the coherence barrier - A new theory for compressed sensing Anders C. Hansen (Cambridge) Joint work with: B. Adcock (Simon Fraser University) C. Poon (Cambridge) B. Roman (Cambridge) UCL-Duke Workshop, September 4, / 42
2 Compressed Sensing in Inverse Problems Typical analog/infinite-dimensional inverse problem where compressed sensing is/can be used: (i) Magnetic Resonance Imaging (MRI) (ii) X-ray Computed Tomography (iii) Thermoacoustic and Photoacoustic Tomography (iv) Single Photon Emission Computerized Tomography (v) Electrical Impedance Tomography (vi) Electron Microscopy (vii) Reflection seismology (viii) Radio interferometry (ix) Fluorescence Microscopy 2 / 42
3 Compressed Sensing in Inverse Problems Most of these problems are modelled by the Fourier transform Ff (ω) = f (x)e 2πiω x dx, R d or the Radon transform Rf : S R C (where S denotes the circle) Rf (θ, p) = f (x) dm(x), x,θ =p where dm denotes Lebesgue measure on the hyperplane {x : x, θ = p}. Fourier slice theorem both problems can be viewed as the problem of reconstructing f from pointwise samples of its Fourier transform. g = Ff, f L 2 (R d ). (1) 3 / 42
4 Compressed Sensing Given the linear system If Solve Ux 0 = y. min z 1 subject to P Ω Uz = P Ω y, where P Ω is a projection and Ω {1,..., N} is subsampled with Ω = m. m C N µ(u) s log(ɛ 1 ) log (N). then P(z = x 0 ) 1 ɛ, where µ(u) = max U i,j 2 i,j is referred to as the incoherence parameter. 4 / 42
5 Pillars of Compressed Sensing Sparsity Incoherence Uniform Random Subsampling In addition: The Restricted Isometry Property + uniform recovery. Problem: These concepts are absent in virtually all the problems listed above. Moreover, uniform random subsampling gives highly suboptimal results. Compressed sensing is currently used with great success in many of these fields, however the current theory does not cover this. 5 / 42
6 Uniform Random Subsampling 1 U = Udft Vdwt. 5% subsamp-map Reconstruction Enlarged 6 / 42
7 Sparsity The classical idea of sparsity in compressed sensing is that there are s important coefficients in the vector x 0 that we want to recover. The location of these coefficients is arbitrary. 7 / 42
8 Sparsity and the Flip Test Let x = and y = U df x, A = P Ω U df V 1 dw, where P Ω is a projection and Ω {1,..., N} is subsampled with Ω = m. Solve min z 1 subject to Az = P Ω y. 8 / 42
9 Sparsity - The Flip Test 2 Truncated (max = ) Figure : x10 Wavelet coefficients and subsampling reconstructions from 10% of Fourier coefficients with distributions (1 + ω12 + ω22 ) 1 and (1 + ω12 + ω22 ) 3/2. If sparsity is the right model we should be able to flip the coefficients. Let 2 Truncated (max = ) 1.5 zf = x10 9 / 42
10 Sparsity - The Flip Test Let Solve to get ẑ f. min z 1 ỹ = U df V 1 dw z f subject to Az = P Ω ỹ Flip the coefficients of ẑ f back to get ẑ, and let ˆx = V 1 dw ẑ. If the ordering of the wavelet coefficients did not matter i.e. sparsity is the right model, then ˆx should be close to x. 10 / 42
11 Sparsity- The Flip Test: Results Figure : The reconstructions from the reversed coefficients. Conclusion: The ordering of the coefficients did matter. Moreover, this phenomenon happens with all wavelets, curvelets, contourlets and shearlets and any reasonable subsampling scheme. Question: Is sparsity really the right model? 11 / 42
12 Sparsity - The Flip Test 512, 20% U Had V 1 dwt Fluorescence Microscopy CS reconstr. CS reconstr, w/ flip Subsampling coeffs. pattern 1024, 12% U Had V 1 dwt Compressive Imaging, Hadamard Spectroscopy 12 / 42
13 Sparsity - The Flip Test (contd.) CS reconstr. CS reconstr, w/ flip coeffs. Subsampling pattern 1024, 20% 1 Udft Vdwt Magnetic Resonance Imaging 512, 12% 1 Udft Vdwt Tomography, Electron Microscopy 13 / 42
14 Sparsity - The Flip Test (contd.) 1024, 10% U dft V 1 dwt Radio interferometry CS reconstr. CS reconstr, w/ flip Subsampling coeffs. pattern 14 / 42
15 Lesson Learned! The optimal sampling strategy depends on the sparsity structure of the signal. Compressed sensing theorems must therefore show how the optimal sampling strategy will depend on the structure. Such theorems currently do not exists. 15 / 42
16 What about the RIP? Did any of the matrices used in the examples satisfy the RIP? 16 / 42
17 Images are not sparse, they are asymptotically sparse How to measure asymptotic sparsity: Suppose f = β j ϕ j. j=1 Let N = k N{M k 1 + 1,..., M k }, where 0 = M 0 < M 1 < M 2 <... and {M k 1 + 1,..., M k } is the set of indices corresponding to the k th scale. Let ɛ (0, 1] and let { s k := s k (ɛ) = min K : K β π(i) ϕ π(i) ɛ i=1 M k i=m k 1 +1 β j ϕ j }, in order words, s k is the effective sparsity at the k th scale. Here π : {1,..., M k M k 1 } {M k 1 + 1,..., M k } is a bijection such that β π(i) β π(i+1). 17 / 42
18 Images are not sparse, they are asymptotically sparse Sparsity, sk(ǫ)/(mk Mk 1) Level 1 Level 2 Level 3 Level 4 Level 5 Level 6 Level 7 Level 8 Worst sparsity Best sparsity Sparsity, sk(ǫ)/(mk Mk 1) Level 1 Level 2 Level 3 Level 4 Level 5 Level 6 Level 7 Level 8 Worst sparsity Best sparsity Relative threshold, ǫ Relative threshold, ǫ Figure : Relative sparsity of Daubechies 8 wavelet coefficients. 18 / 42
19 Analog inverse problems are coherent Let U n = U df V 1 dw Cn n where U df is the discrete Fourier transform and V dw is the discrete wavelet transform. Then µ(u n ) = 1 for all n and all Daubechies wavelets! 19 / 42
20 Analog inverse problems are coherent, why? Note that where Thus, we will always have WOT-lim n U dfv 1 dw = U, ϕ 1, ψ 1 ϕ 2, ψ 1 U = ϕ 1, ψ 2 ϕ 2, ψ µ(u df V 1 dw ) c. 20 / 42
21 Analog inverse problems are asymptotically incoherent Fourier to DB4 Fourier to Legendre Polynomials Figure : Plots of the absolute values of the entries of the matrix U 21 / 42
22 We need a new theory Such theory must incorporates asymptotic sparsity and asymptotic incoherence. It must explain the two intriguing phenomena observed in practice: The optimal sampling strategy is signal structure dependent The success of compressed sensing is resolution dependent The theory cannot be RIP based (at least not with the classical definition of the RIP) 22 / 42
23 New Pillars of Compressed Sensing Asymptotic Sparsity Asymptotic Incoherence Multi-level Subsampling 23 / 42
24 Sparsity in levels Definition For r N let M = (M 1,..., M r ) N r with 1 M 1 <... < M r and s = (s 1,..., s r ) N r, with s k M k M k 1, k = 1,..., r, where M 0 = 0. We say that β l 2 (N) is (s, M)-sparse if, for each k = 1,..., r, k := supp(β) {M k 1 + 1,..., M k }, satisfies k s k. We denote the set of (s, M)-sparse vectors by Σ s,m. 24 / 42
25 Sparsity in levels Definition Let f = j N β jϕ j H, where β = (β j ) j N l 1 (N). Let σ s,m (f ) := min β η l 1. (2) η Σ s,m 25 / 42
26 Multi-level sampling scheme Definition Let r N, N = (N 1,..., N r ) N r with 1 N 1 <... < N r, m = (m 1,..., m r ) N r, with m k N k N k 1, k = 1,..., r, and suppose that Ω k {N k 1 + 1,..., N k }, Ω k = m k, k = 1,..., r, are chosen uniformly at random, where N 0 = 0. We refer to the set Ω = Ω N,m := Ω 1... Ω r. as an (N, m)-multilevel sampling scheme. 26 / 42
27 Local coherence Definition Let U C N N. If N = (N 1,..., N r ) N r and M = (M 1,..., M r ) N r with 1 N 1 <... N r and 1 M 1 <... < M r we define the (k, l) th local coherence of U with respect to N and M by µ N,M (k, l) = where N 0 = M 0 = 0. µ(p N k 1 N k UP M l 1 M l ) µ(p N k 1 N k U), k, l = 1,..., r, 27 / 42
28 The optimization problem inf η l 1 subject to P Ω Uη y δ. (3) η l 1 (N) 28 / 42
29 Theoretical Results Let U C N N be an isometry and β C N. Suppose that Ω = Ω N,m is a multilevel sampling scheme, where N = (N 1,..., N r ) N r and m = (m 1,..., m r ) N r. Let (s, M), where M = (M 1,..., M r ) N r, M 1 <... < M r, and s = (s 1,..., s r ) N r, be any pair such that the following holds: for ɛ > 0 and 1 k r, ( r ) 1 N k N k 1 log(ɛ 1 ) µ N,M (k, l) s l log (N). (4) m k Suppose that ξ C N is a minimizer of (3) with δ = δ K 1 and K = max 1 k r {(N k N k 1 )/m k }. Then, with probability exceeding 1 sɛ, where s = s s r, we have that ξ β C ( δ (1 + L s ) ) + σ s,m (f ), for some constant C, where σ s,m (f ) is as in (2), L = 1 + K = max k=1,...,r { Nk N k 1 m k }. l=1 log 2(6ɛ 1 ) log 2 (4KM and s) 29 / 42
30 Fourier to wavelets Number of samples in each level: m k log(ɛ 1 ) log(n) Nk N k 1 N k 1 ( k 2 ŝ k + s l 2 α(k l) + l=1 where ŝ k = max{s k 1, s k, s k+1 }. Note that r l=k+2 s l 2 v(l k) ) m s log(n), m = m m r, s = s s r. Note that this shows how the success of CS will be resolution dependent., 30 / 42
31 r-level Sampling Scheme Figure : The typical sampling pattern that will be used. 31 / 42
32 Resolution Dependence, 5% subsampling Size: , Error = 10.8% Original CS reconstruction Subsamp. map 32 / 42
33 full sampling and 5% subsampling (DB4) MRI Data courtesy of Andy Ellison, Boston University. Numerics taken from: On asymptotic structure in compressed sensing, B. Roman, B. Adcock, A. C. Hansen, arxiv: / 42
34 The GLPU-Phantom The Guerquin-Kern,Lejeune, Pruessmann, Unser-Phantom (ETH and EPFL) 34 / 42
35 Seeing further with compressed sensing Figure : The figure shows full sampling (= samples) with zero padding. Numerics taken from: On asymptotic structure in compressed sensing, B. Roman, B. Adcock, A. C. Hansen, arxiv: / 42
36 Seeing further with compressed sensing Figure : The figure shows 6.25% subsampling from (= samples) and DB4. Numerics taken from: On asymptotic structure in compressed sensing, B. Roman, B. Adcock, A. C. Hansen, arxiv: / 42
37 Key Question Question: Will compressed sensing ever be used in commercial MRI scanners? 37 / 42
38 Siemens Siemens has implemented our experiments and verified our theory experimentally on their scanners. See the Siemens report: Novel Sampling Strategies for Sparse MR Image Reconstruction, in Proceedings of the International Society for Magnetic Resonance in Medicine. 38 / 42
39 Siemens Conclusion: Significant differences in the spatial resolution can be observed. The image resolution has been greatly improved. Current results practically demonstrated that it is possible to break the coherence barrier by increasing the spatial resolution in MR acquisitions. This likewise implies that the full potential of the compressed sensing is unleashed only if asymptotic sparsity and asymptotic incoherence is achieved. Therefore, compressed sensing might better be used to increase the spatial resolution rather than accelerating the data acquisition in the context of non-dynamic 3D MR imaging. 39 / 42
40 Key Question Question: Will Compressed sensing ever be used in commercial MRI scanners? Answer: We are closer than ever. 40 / 42
41 Other work Structured sampling: Calderbank, Carin, Carson, Chen, Rodrigues (2012) Structured sparsity: Tsaig & Donoho (2006), Eldar (2009), He & Carin (2009), Baraniuk et al. (2010), Krzakala et al. (2011), Duarte & Eldar (2011), Som & Schniter (2012), Renna et al. (2013), Chen et al. (2013), Donoho & Montanari (2013) + others Existing algorithms: Model-based CS: Baranuik, Cevher, Duarte, Hegde (2010) Bayesian CS: Ji, Xue & Carin (2008), He & Carin (2009) Turbo AMP: Som & Schniter (2012) Variable density sampling: Krahmer & Ward (2013), Kutyniok & Lim (2014) 41 / 42
42 Related Papers For details see [4, 5, 2, 1, 3] B. Adcock and A. C. Hansen. Generalized sampling and infinite-dimensional compressed sensing. Found. Comp. Math., (in revision). B. Adcock, A. C. Hansen, C. Poon, and B. Roman. Breaking the coherence barrier: A new theory for compressed sensing. arxiv: , A. C. Hansen. On the solvability complexity index, the n-pseudospectrum and approximations of spectra of operators. J. Amer. Math. Soc., 24(1):81 124, B. Roman, B. Adcock, and A. C. Hansen. On asymptotic structure in compressed sensing. arxiv: , Q. Wang, M. Zenge, H. E. Cetingul, E. Mueller, and M. S. Nadar. Novel sampling strategies for sparse mr image reconstruction. Proc. Int. Soc. Mag. Res. in Med., (22), / 42
On the coherence barrier and analogue problems in compressed sensing
On the coherence barrier and analogue problems in compressed sensing Clarice Poon University of Cambridge June 1, 2017 Joint work with: Ben Adcock Anders Hansen Bogdan Roman (Simon Fraser) (Cambridge)
More informationCompressed sensing and imaging
Compressed sensing and imaging The effect and benefits of local structure Ben Adcock Department of Mathematics Simon Fraser University 1 / 45 Overview Previously: An introduction to compressed sensing.
More informationBeyond incoherence and beyond sparsity: compressed sensing in the real world
Beyond incoherence and beyond sparsity: compressed sensing in the real world Clarice Poon 1st November 2013 University of Cambridge, UK Applied Functional and Harmonic Analysis Group Head of Group Anders
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationCompressed sensing (CS), introduced by Candès, Romberg &
On asymptotic structure in compressed sensing Bogdan Roman, Ben Adcock and Anders Hansen Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK, and Department
More informationCompressed sensing (CS), introduced by Candès, Romberg &
On asymptotic structure in compressed sensing Bogdan Roman, Ben Adcock and Anders Hansen Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK, and Department
More informationLecture 22: More On Compressed Sensing
Lecture 22: More On Compressed Sensing Scribed by Eric Lee, Chengrun Yang, and Sebastian Ament Nov. 2, 207 Recap and Introduction Basis pursuit was the method of recovering the sparsest solution to an
More informationLEARNING DATA TRIAGE: LINEAR DECODING WORKS FOR COMPRESSIVE MRI. Yen-Huan Li and Volkan Cevher
LARNING DATA TRIAG: LINAR DCODING WORKS FOR COMPRSSIV MRI Yen-Huan Li and Volkan Cevher Laboratory for Information Inference Systems École Polytechnique Fédérale de Lausanne ABSTRACT The standard approach
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationIntroduction How it works Theory behind Compressed Sensing. Compressed Sensing. Huichao Xue. CS3750 Fall 2011
Compressed Sensing Huichao Xue CS3750 Fall 2011 Table of Contents Introduction From News Reports Abstract Definition How it works A review of L 1 norm The Algorithm Backgrounds for underdetermined linear
More informationInterpolation via weighted l 1 minimization
Interpolation via weighted l 1 minimization Rachel Ward University of Texas at Austin December 12, 2014 Joint work with Holger Rauhut (Aachen University) Function interpolation Given a function f : D C
More informationNear Ideal Behavior of a Modified Elastic Net Algorithm in Compressed Sensing
Near Ideal Behavior of a Modified Elastic Net Algorithm in Compressed Sensing M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas M.Vidyasagar@utdallas.edu www.utdallas.edu/ m.vidyasagar
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More informationAN INTRODUCTION TO COMPRESSIVE SENSING
AN INTRODUCTION TO COMPRESSIVE SENSING Rodrigo B. Platte School of Mathematical and Statistical Sciences APM/EEE598 Reverse Engineering of Complex Dynamical Networks OUTLINE 1 INTRODUCTION 2 INCOHERENCE
More informationCompressed Sensing and Neural Networks
and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications
More informationSensing systems limited by constraints: physical size, time, cost, energy
Rebecca Willett Sensing systems limited by constraints: physical size, time, cost, energy Reduce the number of measurements needed for reconstruction Higher accuracy data subject to constraints Original
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationA practical guide to the recovery of wavelet coefficients from Fourier measurements
A practical guide to the recovery of wavelet coefficients from Fourier measurements Milana Gataric Clarice Poon April 05; Revised March 06 Abstract In a series of recent papers Adcock, Hansen and Poon,
More informationLinear reconstructions and the analysis of the stable sampling rate
Preprint submitted to SAMPLING THEORY IN SIGNAL AND IMAGE PROCESSING DVI file produced on Linear reconstructions and the analysis of the stable sampling rate Laura Terhaar Department of Applied Mathematics
More informationDeterministic Compressed Sensing for Efficient Image Reconstruction
Chirp for Efficient Image University of Idaho Pacific Northwest Numerical Analysis Seminar 27 October 2012 Chirp Collaborators: Connie Ni (HRL Lab, Malibu, CA) Prasun Mahanti (Arizona State University)
More informationDeterministic Compressed Sensing for Efficient Image Reconstruction
Chirp for Efficient Image Joint with: Kangyu Ni, Prasun Mahanti, Svetlana Roudenko, Doug Cochran SCIMM 6 February 2010 Chirp - Overview R n { }} { y }{{} measured = n N { }} { A }{{} s sensing matrix }{{}
More informationCompressive Sensing under Matrix Uncertainties: An Approximate Message Passing Approach
Compressive Sensing under Matrix Uncertainties: An Approximate Message Passing Approach Asilomar 2011 Jason T. Parker (AFRL/RYAP) Philip Schniter (OSU) Volkan Cevher (EPFL) Problem Statement Traditional
More informationSparsity Regularization
Sparsity Regularization Bangti Jin Course Inverse Problems & Imaging 1 / 41 Outline 1 Motivation: sparsity? 2 Mathematical preliminaries 3 l 1 solvers 2 / 41 problem setup finite-dimensional formulation
More informationRandomness-in-Structured Ensembles for Compressed Sensing of Images
Randomness-in-Structured Ensembles for Compressed Sensing of Images Abdolreza Abdolhosseini Moghadam Dep. of Electrical and Computer Engineering Michigan State University Email: abdolhos@msu.edu Hayder
More informationCompressed Sensing: Extending CLEAN and NNLS
Compressed Sensing: Extending CLEAN and NNLS Ludwig Schwardt SKA South Africa (KAT Project) Calibration & Imaging Workshop Socorro, NM, USA 31 March 2009 Outline 1 Compressed Sensing (CS) Introduction
More informationCS 229r: Algorithms for Big Data Fall Lecture 19 Nov 5
CS 229r: Algorithms for Big Data Fall 215 Prof. Jelani Nelson Lecture 19 Nov 5 Scribe: Abdul Wasay 1 Overview In the last lecture, we started discussing the problem of compressed sensing where we are given
More informationUniqueness Conditions for A Class of l 0 -Minimization Problems
Uniqueness Conditions for A Class of l 0 -Minimization Problems Chunlei Xu and Yun-Bin Zhao October, 03, Revised January 04 Abstract. We consider a class of l 0 -minimization problems, which is to search
More informationCombining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation
UIUC CSL Mar. 24 Combining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation Yuejie Chi Department of ECE and BMI Ohio State University Joint work with Yuxin Chen (Stanford).
More informationMultiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing
Multiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing ibeltal F. Alem, Daniel H. Chae, and Rodney A. Kennedy The Australian National
More informationCompressed Sensing. 1 Introduction. 2 Design of Measurement Matrices
Compressed Sensing Yonina C. Eldar Electrical Engineering Department, Technion-Israel Institute of Technology, Haifa, Israel, 32000 1 Introduction Compressed sensing (CS) is an exciting, rapidly growing
More informationAn Overview of Compressed Sensing
An Overview of Compressed Sensing Nathan Schneider November 18, 2009 Abstract In a large number of applications, the system will be designed to sample at a rate equal to at least the frequency bandwidth
More informationLecture Notes 9: Constrained Optimization
Optimization-based data analysis Fall 017 Lecture Notes 9: Constrained Optimization 1 Compressed sensing 1.1 Underdetermined linear inverse problems Linear inverse problems model measurements of the form
More informationPart IV Compressed Sensing
Aisenstadt Chair Course CRM September 2009 Part IV Compressed Sensing Stéphane Mallat Centre de Mathématiques Appliquées Ecole Polytechnique Conclusion to Super-Resolution Sparse super-resolution is sometime
More informationKomprimované snímání a LASSO jako metody zpracování vysocedimenzionálních dat
Komprimované snímání a jako metody zpracování vysocedimenzionálních dat Jan Vybíral (Charles University Prague, Czech Republic) November 2014 VUT Brno 1 / 49 Definition and motivation Its use in bioinformatics
More informationTurbo-AMP: A Graphical-Models Approach to Compressive Inference
Turbo-AMP: A Graphical-Models Approach to Compressive Inference Phil Schniter (With support from NSF CCF-1018368 and DARPA/ONR N66001-10-1-4090.) June 27, 2012 1 Outline: 1. Motivation. (a) the need for
More informationCompressed Sensing and Sparse Recovery
ELE 538B: Sparsity, Structure and Inference Compressed Sensing and Sparse Recovery Yuxin Chen Princeton University, Spring 217 Outline Restricted isometry property (RIP) A RIPless theory Compressed sensing
More informationRecovery of Compressible Signals in Unions of Subspaces
1 Recovery of Compressible Signals in Unions of Subspaces Marco F. Duarte, Chinmay Hegde, Volkan Cevher, and Richard G. Baraniuk Department of Electrical and Computer Engineering Rice University Abstract
More informationA Survey of Compressive Sensing and Applications
A Survey of Compressive Sensing and Applications Justin Romberg Georgia Tech, School of ECE ENS Winter School January 10, 2012 Lyon, France Signal processing trends DSP: sample first, ask questions later
More informationCompressive Sensing of Sparse Tensor
Shmuel Friedland Univ. Illinois at Chicago Matheon Workshop CSA2013 December 12, 2013 joint work with Q. Li and D. Schonfeld, UIC Abstract Conventional Compressive sensing (CS) theory relies on data representation
More informationSparse recovery for spherical harmonic expansions
Rachel Ward 1 1 Courant Institute, New York University Workshop Sparsity and Cosmology, Nice May 31, 2011 Cosmic Microwave Background Radiation (CMB) map Temperature is measured as T (θ, ϕ) = k k=0 l=
More informationImproved Radon Based Imaging using the Shearlet Transform
Improved Radon Based Imaging using the Shearlet Transform Glenn R. Easley a, Flavia Colonna b, Demetrio Labate c a System Planning Corporation, Arlington, Virginia b George Mason University, Fairfax, Virginia
More informationarxiv: v1 [cs.it] 21 Feb 2013
q-ary Compressive Sensing arxiv:30.568v [cs.it] Feb 03 Youssef Mroueh,, Lorenzo Rosasco, CBCL, CSAIL, Massachusetts Institute of Technology LCSL, Istituto Italiano di Tecnologia and IIT@MIT lab, Istituto
More informationAn iterative hard thresholding estimator for low rank matrix recovery
An iterative hard thresholding estimator for low rank matrix recovery Alexandra Carpentier - based on a joint work with Arlene K.Y. Kim Statistical Laboratory, Department of Pure Mathematics and Mathematical
More informationExponential decay of reconstruction error from binary measurements of sparse signals
Exponential decay of reconstruction error from binary measurements of sparse signals Deanna Needell Joint work with R. Baraniuk, S. Foucart, Y. Plan, and M. Wootters Outline Introduction Mathematical Formulation
More informationIntroduction to Compressed Sensing
Introduction to Compressed Sensing Alejandro Parada, Gonzalo Arce University of Delaware August 25, 2016 Motivation: Classical Sampling 1 Motivation: Classical Sampling Issues Some applications Radar Spectral
More informationPhase Transition Phenomenon in Sparse Approximation
Phase Transition Phenomenon in Sparse Approximation University of Utah/Edinburgh L1 Approximation: May 17 st 2008 Convex polytopes Counting faces Sparse Representations via l 1 Regularization Underdetermined
More information2D X-Ray Tomographic Reconstruction From Few Projections
2D X-Ray Tomographic Reconstruction From Few Projections Application of Compressed Sensing Theory CEA-LID, Thalès, UJF 6 octobre 2009 Introduction Plan 1 Introduction 2 Overview of Compressed Sensing Theory
More informationSignal Recovery, Uncertainty Relations, and Minkowski Dimension
Signal Recovery, Uncertainty Relations, and Minkowski Dimension Helmut Bőlcskei ETH Zurich December 2013 Joint work with C. Aubel, P. Kuppinger, G. Pope, E. Riegler, D. Stotz, and C. Studer Aim of this
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationLecture 13 October 6, Covering Numbers and Maurey s Empirical Method
CS 395T: Sublinear Algorithms Fall 2016 Prof. Eric Price Lecture 13 October 6, 2016 Scribe: Kiyeon Jeon and Loc Hoang 1 Overview In the last lecture we covered the lower bound for p th moment (p > 2) and
More informationAnders Christian Hansen
University of Cambridge DAMTP, Cambridge, UK a.hansen@damtp.cam.ac.uk Education University of Cambridge, King s College, UK. Ph.D, Mathematics, June 2008. University of California, Berkeley, USA. M.A.,
More informationPerformance Regions in Compressed Sensing from Noisy Measurements
Performance egions in Compressed Sensing from Noisy Measurements Junan Zhu and Dror Baron Department of Electrical and Computer Engineering North Carolina State University; aleigh, NC 27695, USA Email:
More informationStochastic geometry and random matrix theory in CS
Stochastic geometry and random matrix theory in CS IPAM: numerical methods for continuous optimization University of Edinburgh Joint with Bah, Blanchard, Cartis, and Donoho Encoder Decoder pair - Encoder/Decoder
More informationCompressive Sensing (CS)
Compressive Sensing (CS) Luminita Vese & Ming Yan lvese@math.ucla.edu yanm@math.ucla.edu Department of Mathematics University of California, Los Angeles The UCLA Advanced Neuroimaging Summer Program (2014)
More informationIEEE SIGNAL PROCESSING LETTERS, VOL. 22, NO. 9, SEPTEMBER
IEEE SIGNAL PROCESSING LETTERS, VOL. 22, NO. 9, SEPTEMBER 2015 1239 Preconditioning for Underdetermined Linear Systems with Sparse Solutions Evaggelia Tsiligianni, StudentMember,IEEE, Lisimachos P. Kondi,
More informationDeterminis)c Compressed Sensing for Images using Chirps and Reed- Muller Sequences
Determinis)c Compressed Sensing for Images using Chirps and Reed- Muller Sequences Kangyu Ni Mathematics Arizona State University Joint with Somantika Datta*, Prasun Mahanti, Svetlana Roudenko, Douglas
More informationFast Hard Thresholding with Nesterov s Gradient Method
Fast Hard Thresholding with Nesterov s Gradient Method Volkan Cevher Idiap Research Institute Ecole Polytechnique Federale de ausanne volkan.cevher@epfl.ch Sina Jafarpour Department of Computer Science
More informationAn Adaptive Sublinear Time Block Sparse Fourier Transform
An Adaptive Sublinear Time Block Sparse Fourier Transform Volkan Cevher Michael Kapralov Jonathan Scarlett Amir Zandieh EPFL February 8th 217 Given x C N, compute the Discrete Fourier Transform (DFT) of
More informationGeneralized Orthogonal Matching Pursuit- A Review and Some
Generalized Orthogonal Matching Pursuit- A Review and Some New Results Department of Electronics and Electrical Communication Engineering Indian Institute of Technology, Kharagpur, INDIA Table of Contents
More informationCOMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY
COMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY Adriana González ICTEAM/UCL March 26th, 214 1 ISPGroup - ICTEAM - UCL http://sites.uclouvain.be/ispgroup Université catholique de Louvain, Louvain-la-Neuve,
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Mathematical Analysis and Applications Workshop in honor of Rupert Lasser Helmholtz
More informationCompressed Sensing: Lecture I. Ronald DeVore
Compressed Sensing: Lecture I Ronald DeVore Motivation Compressed Sensing is a new paradigm for signal/image/function acquisition Motivation Compressed Sensing is a new paradigm for signal/image/function
More informationApproximate Message Passing for Bilinear Models
Approximate Message Passing for Bilinear Models Volkan Cevher Laboratory for Informa4on and Inference Systems LIONS / EPFL h,p://lions.epfl.ch & Idiap Research Ins=tute joint work with Mitra Fatemi @Idiap
More informationApproximate Message Passing Algorithm for Complex Separable Compressed Imaging
Approximate Message Passing Algorithm for Complex Separle Compressed Imaging Akira Hirayashi, Jumpei Sugimoto, and Kazushi Mimura College of Information Science and Engineering, Ritsumeikan University,
More informationFast Reconstruction Algorithms for Deterministic Sensing Matrices and Applications
Fast Reconstruction Algorithms for Deterministic Sensing Matrices and Applications Program in Applied and Computational Mathematics Princeton University NJ 08544, USA. Introduction What is Compressive
More informationRapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization
Rapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization Shuyang Ling Department of Mathematics, UC Davis Oct.18th, 2016 Shuyang Ling (UC Davis) 16w5136, Oaxaca, Mexico Oct.18th, 2016
More informationCompressed sensing for radio interferometry: spread spectrum imaging techniques
Compressed sensing for radio interferometry: spread spectrum imaging techniques Y. Wiaux a,b, G. Puy a, Y. Boursier a and P. Vandergheynst a a Institute of Electrical Engineering, Ecole Polytechnique Fédérale
More informationCS on CS: Computer Science insights into Compresive Sensing (and vice versa) Piotr Indyk MIT
CS on CS: Computer Science insights into Compresive Sensing (and vice versa) Piotr Indyk MIT Sparse Approximations Goal: approximate a highdimensional vector x by x that is sparse, i.e., has few nonzero
More informationModel-Based Compressive Sensing for Signal Ensembles. Marco F. Duarte Volkan Cevher Richard G. Baraniuk
Model-Based Compressive Sensing for Signal Ensembles Marco F. Duarte Volkan Cevher Richard G. Baraniuk Concise Signal Structure Sparse signal: only K out of N coordinates nonzero model: union of K-dimensional
More informationOverview. Optimization-Based Data Analysis. Carlos Fernandez-Granda
Overview Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 1/25/2016 Sparsity Denoising Regression Inverse problems Low-rank models Matrix completion
More informationSparse Solutions of an Undetermined Linear System
1 Sparse Solutions of an Undetermined Linear System Maddullah Almerdasy New York University Tandon School of Engineering arxiv:1702.07096v1 [math.oc] 23 Feb 2017 Abstract This work proposes a research
More informationConstruction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling
Construction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling Rujie Yin Department of Mathematics Duke University USA Email: rujie.yin@duke.edu arxiv:1602.04882v1 [math.fa] 16 Feb
More informationBhaskar Rao Department of Electrical and Computer Engineering University of California, San Diego
Bhaskar Rao Department of Electrical and Computer Engineering University of California, San Diego 1 Outline Course Outline Motivation for Course Sparse Signal Recovery Problem Applications Computational
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationJoint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces
Aalborg Universitet Joint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces Christensen, Mads Græsbøll; Nielsen, Jesper Kjær Published in: I E E E / S P Workshop
More informationApproximation theory in neural networks
Approximation theory in neural networks Yanhui Su yanhui su@brown.edu March 30, 2018 Outline 1 Approximation of functions by a sigmoidal function 2 Approximations of continuous functionals by a sigmoidal
More informationCompressive Sensing Applied to Full-wave Form Inversion
Compressive Sensing Applied to Full-wave Form Inversion Felix J. Herrmann* fherrmann@eos.ubc.ca Joint work with Yogi Erlangga, and Tim Lin *Seismic Laboratory for Imaging & Modeling Department of Earth
More informationCOMPRESSED Sensing (CS) is a method to recover a
1 Sample Complexity of Total Variation Minimization Sajad Daei, Farzan Haddadi, Arash Amini Abstract This work considers the use of Total Variation (TV) minimization in the recovery of a given gradient
More informationSparsifying Transform Learning for Compressed Sensing MRI
Sparsifying Transform Learning for Compressed Sensing MRI Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and Coordinated Science Laborarory University of Illinois
More informationA Tutorial on Compressive Sensing. Simon Foucart Drexel University / University of Georgia
A Tutorial on Compressive Sensing Simon Foucart Drexel University / University of Georgia CIMPA13 New Trends in Applied Harmonic Analysis Mar del Plata, Argentina, 5-16 August 2013 This minicourse acts
More informationCompressed Sensing Using Reed- Solomon and Q-Ary LDPC Codes
Compressed Sensing Using Reed- Solomon and Q-Ary LDPC Codes Item Type text; Proceedings Authors Jagiello, Kristin M. Publisher International Foundation for Telemetering Journal International Telemetering
More informationAn Introduction to Sparse Representations and Compressive Sensing. Part I
An Introduction to Sparse Representations and Compressive Sensing Part I Paulo Gonçalves CPE Lyon - 4ETI - Cours Semi-Optionnel Méthodes Avancées pour le Traitement des Signaux 2014 Objectifs Part I The
More informationMessage Passing Algorithms for Compressed Sensing: II. Analysis and Validation
Message Passing Algorithms for Compressed Sensing: II. Analysis and Validation David L. Donoho Department of Statistics Arian Maleki Department of Electrical Engineering Andrea Montanari Department of
More informationCoSaMP. Iterative signal recovery from incomplete and inaccurate samples. Joel A. Tropp
CoSaMP Iterative signal recovery from incomplete and inaccurate samples Joel A. Tropp Applied & Computational Mathematics California Institute of Technology jtropp@acm.caltech.edu Joint with D. Needell
More informationCompressed sensing. Or: the equation Ax = b, revisited. Terence Tao. Mahler Lecture Series. University of California, Los Angeles
Or: the equation Ax = b, revisited University of California, Los Angeles Mahler Lecture Series Acquiring signals Many types of real-world signals (e.g. sound, images, video) can be viewed as an n-dimensional
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Instructor: Moritz Hardt Email: hardt+ee227c@berkeley.edu Graduate Instructor: Max Simchowitz Email: msimchow+ee227c@berkeley.edu
More informationConstrained optimization
Constrained optimization DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Compressed sensing Convex constrained
More informationCompressive sensing of low-complexity signals: theory, algorithms and extensions
Compressive sensing of low-complexity signals: theory, algorithms and extensions Laurent Jacques March 7, 9, 1, 14, 16 and 18, 216 9h3-12h3 (incl. 3 ) Graduate School in Systems, Optimization, Control
More informationInvertible Nonlinear Dimensionality Reduction via Joint Dictionary Learning
Invertible Nonlinear Dimensionality Reduction via Joint Dictionary Learning Xian Wei, Martin Kleinsteuber, and Hao Shen Department of Electrical and Computer Engineering Technische Universität München,
More informationA new method on deterministic construction of the measurement matrix in compressed sensing
A new method on deterministic construction of the measurement matrix in compressed sensing Qun Mo 1 arxiv:1503.01250v1 [cs.it] 4 Mar 2015 Abstract Construction on the measurement matrix A is a central
More informationSuper-resolution by means of Beurling minimal extrapolation
Super-resolution by means of Beurling minimal extrapolation Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Acknowledgements ARO W911NF-16-1-0008
More informationMathematical Methods in Machine Learning
UMD, Spring 2016 Outline Lecture 2: Role of Directionality 1 Lecture 2: Role of Directionality Anisotropic Harmonic Analysis Harmonic analysis decomposes signals into simpler elements called analyzing
More informationLEARNING OVERCOMPLETE SPARSIFYING TRANSFORMS FOR SIGNAL PROCESSING. Saiprasad Ravishankar and Yoram Bresler
LEARNING OVERCOMPLETE SPARSIFYING TRANSFORMS FOR SIGNAL PROCESSING Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and the Coordinated Science Laboratory, University
More informationCOMPRESSED SENSING IN PYTHON
COMPRESSED SENSING IN PYTHON Sercan Yıldız syildiz@samsi.info February 27, 2017 OUTLINE A BRIEF INTRODUCTION TO COMPRESSED SENSING A BRIEF INTRODUCTION TO CVXOPT EXAMPLES A Brief Introduction to Compressed
More informationCompressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery
Compressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery Jorge F. Silva and Eduardo Pavez Department of Electrical Engineering Information and Decision Systems Group Universidad
More informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationSparse solutions of underdetermined systems
Sparse solutions of underdetermined systems I-Liang Chern September 22, 2016 1 / 16 Outline Sparsity and Compressibility: the concept for measuring sparsity and compressibility of data Minimum measurements
More informationTree-Structured Compressive Sensing with Variational Bayesian Analysis
1 Tree-Structured Compressive Sensing with Variational Bayesian Analysis Lihan He, Haojun Chen and Lawrence Carin Department of Electrical and Computer Engineering Duke University Durham, NC 27708-0291
More informationSolution Recovery via L1 minimization: What are possible and Why?
Solution Recovery via L1 minimization: What are possible and Why? Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Eighth US-Mexico Workshop on Optimization
More informationIntroduction to Sparsity. Xudong Cao, Jake Dreamtree & Jerry 04/05/2012
Introduction to Sparsity Xudong Cao, Jake Dreamtree & Jerry 04/05/2012 Outline Understanding Sparsity Total variation Compressed sensing(definition) Exact recovery with sparse prior(l 0 ) l 1 relaxation
More information