Sparse linear models
|
|
- Harriet Houston
- 5 years ago
- Views:
Transcription
1 Sparse linear models Optimization-Based Data Analysis Carlos Fernandez-Granda 2/22/2016
2 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
3 Linear model A linear model for a signal x R n is a representation of the form x = m c i φ i i=1 {φ 1,..., φ m } is a family of atoms in R n c R m is an alternative representation or transform of x
4 Sparse linear model A sparse linear model contains a small number of coefficients x = i S c i φ i S m How do we choose a transform that sparsifies a class of signals? Intuition / Domain knowledge (this lecture) Learning it from the data (later on)
5 Why? Signals of interest (speech, natural images, biomedical activity, etc.) are often highly structured Sparse linear models are able to exploit this structure to enhance data analysis and processing Applications: Compression Denoising Inverse problems
6 Sparse representation in an orthonormal basis If the atoms {φ 1,..., φ n } form an orthonormal basis, then U := [ ] φ 1 φ 2 φ n UU T = I Coefficients are obtained by computing inner products with the atoms x = UU T x = m φ i, x φ i i=1
7 Sparse representation in a basis If the atoms {φ 1,..., φ n } form a basis, then B := [ φ 1 φ 2 φ n ] BB 1 = I Coefficients are obtained by computing inner products with dual atoms x = BB 1 x = m θ i, x φ i i=1 where θ T 1 B 1 θ = T 2 θ T n
8 Overcomplete dictionaries If the atoms {φ 1,..., φ m } are linearly independent and m > n D := [ φ 1 φ 2 φ m ] Two alternative sparse models 1. Synthesis sparse model x = Dc Problem: Given x find a sparse c 2. Analysis sparse model: where c is sparse D T x is sparse Both are equivalent if D is an orthonormal basis
9 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
10 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
11 Fourier series Fourier series coefficients of f : [0, 1] C Fourier series S n (t) := c k := n k= n 1 0 f (t) e i2πkt dt c k e i2πkt = n k= n φ k, f φ k Sinusoidal atoms φ k (t) = e i2πkt = cos (2πkt) + i sin (2πkt) Orthonormal basis of L 2 under the usual inner product lim n f (t) S n (t) = 0 for all f L 2
12 Discrete Fourier transform (DFT) Discretized frequency representation for vectors in C n Sinusoidal atoms are a basis for C n φ k = 1 n 1 e i2πk n e i2πk2 n e i2πk(n 1) n DFT {x} k := (Fx) k = φ k, x F := [ φ 0 φ 1 φ n 1 ] The fast Fourier transform (FFT) computes the DFT in O (n log n) The discrete cosine transform (DCT) is a related transformation designed for real vectors
13 Discrete cosine transform Signal DCT coefficients
14 Electrocardiogram
15 Electrocardiogram (spectrum)
16 Electrocardiogram (spectrum)
17 2D DFT Discretized frequency representation for 2D arrays in C n n Sinusoidal atoms are a basis for C n n φ 2D k 1,k 2 = 1 n 1 e i2πk 2 n e i2πk 2 (n 1) n e i2πk 1 n e i2π(k 1 +k 2) n e i2π(k 1 +k 2 (n 1)) n e i2πk 1 (n 1) n e i2π(k 1 (n 1)+k 2) n e i2π(k 1 (n 1)+k 2 (n 1)) n ( ) T = φ 1D k 1 φ 1D k 2 DFT 2D {X } := FXF = n 1 n 1 k 1 =0 k 2 =0 φ 2D k 1,k 2, X φ 2D k 1,k 2 Generalizes to C m n, m n, and to higher dimensions
18 Compression via frequency representation The 2D frequency representation of images tends to be sparse Thresholding the coefficients yields a compressed representation The JPEG compression standard is based on the 2D DCT High-frequency coefficients are discarded according to a perceptual model
19 Compression via frequency representation Original
20 Compression via frequency representation 10 % largest DCT coeffs
21 Compression via frequency representation 2 % largest DCT coeffs
22 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
23 Motivation Spectrum of speech, music, etc. varies over time Idea: Compute frequency representation of time segments of the signal We must use a window to avoid introducing spurious high frequencies
24 The need for windowing Signal Window = Spectrum =
25 The need for windowing Signal Window = Spectrum =
26 Short-time Fourier transform Let w : [0, 1] C be a window function localized in time and frequency STFT {f } (k, τ) := 1 0 f (t) w (t τ)e i2πkt dt = φ k,τ, f Each atom φ k,τ (t) := w (t τ) e i2πkt corresponds to w shifted by τ in time and by k in frequency In discrete time, pointwise multiplication by a shifted window followed by a DFT, equivalent to D T x where D C n m, m > n The STFT of speech tends to be sparse (analysis sparse model) Including dilations of w (in addition to time and frequency translations) yields a dictionary of Gabor atoms
27 Atom τ = 0, k = 0 Real part Imaginary part Spectrum
28 Atom τ = 1/32, k = 0 Real part Imaginary part Spectrum
29 Atom τ = 0, k = 64 Real part Imaginary part Spectrum
30 Atom τ = 1/32, k = 64 Real part Imaginary part Spectrum
31 Speech signal
32 Spectrum
33 Spectrogram (log magnitude of STFT coefficients) Frequency Time
34 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
35 Wavelets Aim: Approximate signals at different scales A wavelet ψ is a unit-norm, zero-mean function in L 2 Wavelet transform W {f } (s, τ) := 1 1 ( t τ f (t) ψ s s 0 ) dt = φ s,τ, f Atoms are dilations and translations of the mother wavelet φ s,τ (t) = 1 ( ) t τ ψ s s We can build orthonormal basis for L 2 using wavelets
36 Multiresolution approximation Sequence {V j, j Z} of closed subspaces of L 2 (R) such that P Vj (f ) is an approximation of f at scale 2 j Conditions: Dilating functions in V j by 2 yields functions in V j+1 f (t) V j f ( t V j+1 2) Approximations at a scale 2 j are always better than at 2 j+1 V j+1 V j
37 Multiresolution approximation V j is invariant to translations at the scale 2 j f (t) V j f ( t 2 j k ) V j for all k Z As j the approximation loses all information lim V j = {0} j As j the approximation is perfect lim j V j = L 2 There exists a scaling function ζ V 0 such that { ζ0,k (t) := ζ (t k), k Z } is an orthonormal basis for V 0
38 { ψ2 j,k, k Z } is an orthonormal basis for V j V j { ζ0,k (t), ψ 2 1,k, ψ 2 2,k,..., ψ 2 j,k, k Z } is an orthonormal basis for V j Wavelet basis Mallat and Meyer prove that there exists a wavelet ψ such that P Vj (f ) = P Vj+1 (f ) + k Z ψ2 j,k, f ψ 2 j,k. Many different wavelet bases: Meyer, Daubechies, Battle-Lemarie,... Discrete wavelet transform can be computed in O (n) Signal processing interpretation: Wavelets act as band-pass filters, scaling functions act as low-pass filters
39 Haar wavelet Scaling function Mother wavelet
40 Electrocardiogram Signal Haar transform
41 Scale 2 9 Contribution Approximation
42 Scale 2 8 Contribution Approximation
43 Scale 2 7 Contribution Approximation
44 Scale 2 6 Contribution Approximation
45 Scale 2 5 Contribution Approximation
46 Scale 2 4 Contribution Approximation
47 Scale 2 3 Contribution Approximation
48 Scale 2 2 Contribution Approximation
49 Scale 2 1 Contribution Approximation
50 Scale 2 0 Contribution Approximation
51 2D Wavelets Extension to 2D by using outer products of 1D atoms φ 2D s 1,s 2,k 1,k 2 := φ 1D s 1,k 1 ( φ 1D s 2,k 2 ) T Yields sparse representation of natural images The JPEG 2000 compression standard is based on 2D wavelets Many extensions: Steerable pyramid, ridgelets, curvelets, bandlets,...
52 2D wavelet transform
53 2D wavelet transform
54 Sorted coefficients
55 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
56 Overcomplete dictionaries Atoms {φ 1,..., φ m } are linearly independent and m > n D := [ φ 1 φ 2 φ m ] Synthesis sparse model: x = Dc where c is sparse Problem: There are infinite choices of c such that x = Dc Some may not be sparse at all! Example: Dictionary of sinusoids
57 Dictionary of sinusoids x = Dc c
58 First idea Apply pseudoinverse ĉ := D x = D T ( DD T ) 1 x Interpretations: Projection of c onto row space of D Solution to minimize c 2 subject to x = D c,
59 Dictionary of sinusoids: Minimum l 2 -norm coefficients
60 Computing sparse representations We would like to solve Computationally intractable minimize c 0 subject to x = D c Two possibilities: Greedy methods: Select atoms one by one l 1 -norm minimization: min c R m c 1 such that x = Dc
61 Matching pursuit (MP) Iteratively choose atoms that are most correlated with the signal Initialization: Iterations: k = 1, 2,... φ (k) = arg max j ˆx (k) = ˆx (k 1) + r (0) = x ˆx (0) = 0 r (k 1), φ k r (k 1), φ (k) φ (k) r (k) = r (k 1) r (k 1), φ (k) φ (k)
62 Dictionary of sinusoids: Coefficients Original MP
63 Dictionary of sinusoids: Approximation Original MP
64 Orthogonal matching pursuit (OMP) Makes sure approximation is orthogonal to residual at Initialization: Iterations: k = 1, 2,... φ (k) = arg max j r (0) = x r (k 1), φ k A (k) = [ φ (1) φ (2) φ (k)] ĉ (k) = A (k) x = ˆx (k) = A (k) ĉ (k) r (k) = x ˆx (k) ( A (k) A (k)t ) 1 A (k)t x
65 Dictionary of sinusoids: Coefficients Original OMP
66 Dictionary of sinusoids: Approximation Original OMP
67 l 1 -norm minimization Estimate coefficients by solving minimize c 1 subject to x = D c Computationally tractable (convex program) Known as basis pursuit in the literature
68 Geometric intuition c 2 Min. l 1 -norm solution Min. l 2 -norm solution c 1 Dc = x
69 Dictionary of sinusoids: Coefficients Original l 1 -norm min.
70 Dictionary of sinusoids: Approximation Original l 1 -norm min.
71 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
72 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
73 Denoising Aim: Extracting information (signal) from data in the presence of uninformative perturbations (noise) Additive noise model data = signal + noise y = x + z Prior knowledge about structure of signal vs structure of noise is required
74 Electrocardiogram Spectrum
75 Electrocardiogram: High-frequency noise (power line hum) Original spectrum Low-pass filtered spectrum
76 Electrocardiogram: High-frequency noise (power line hum) Original spectrum Low-pass filtered spectrum
77 Electrocardiogram: High-frequency noise (power line hum) Original Low-pass filtered
78 Electrocardiogram: High-frequency noise (power line hum) Original Low-pass filtered
79 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
80 Thresholding Prior knowledge: Signal is a sparse superposition of dictionary atoms Noise is not (incoherence between atoms and noise) Hard-thresholding operator H η (x) i := { x i if x i > η, 0 otherwise
81 Denoising via thresholding Data Signal
82 Denoising via thresholding Estimate Signal
83 Sparsity in a basis Assumption: x = Bc, where c is sparse Threshold B 1 y ĉ = H η ( B 1 y ) = H η ( c + B 1 z ) ŷ = Bĉ Noise and sparsifying atoms should be incoherent, i.e. B 1 z is not sparse Example: Orthogonal sparsifying basis and Gaussian noise
84 Denoising via thresholding in DCT basis DCT coefficients Data Signal Data
85 Denoising via thresholding in DCT basis DCT coefficients Estimate Signal Estimate
86 Denoising via thresholding in a wavelet basis
87 Denoising via thresholding in a wavelet basis
88 2D wavelet coefficients
89 Original coefficients
90 Thresholded coefficients
91 Estimate
92 Estimate
93 Denoising via thresholding in a wavelet basis Original Noisy Estimate
94 Analysis model Assumption: D T x is sparse Threshold, then use left inverse of D T L ) ĉ = H η (D T y ) = H η (D T x + D T z ŷ = Lĉ Example: Thresholding STFT coefficients for speech denoising
95 Block thresholding Assumption: Coefficients are group sparse, nonzero coefficients cluster together Block thresholding: Partition coefficients into blocks I 1, I 2,..., I k and threshold whole blocks { x i if i I j such that x Ij 2 > η, B η (x) i := 0 otherwise
96 Haar transform
97 2D wavelet transform
98 Speech denoising
99 Time thresholding
100 Spectrum
101 Frequency thresholding
102 Frequency thresholding Data DFT thresholding
103 Spectrogram (STFT) Frequency Time
104 STFT thresholding Frequency Time
105 STFT thresholding Data STFT thresholding
106 STFT block thresholding Frequency Time
107 STFT block thresholding Data STFT block thresh.
108 Wavelets
109 Sorted wavelet coefficients
110 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
111 Denoising via l 1 -norm regularized least squares Synthesis sparse model x = Dc where c is sparse We would like to solve minimize c 0 subject to y D c Computationally intractable if D is overcomplete Basis-pursuit denoising ĉ = arg min x R m y D c λ c 1 ˆx = Dĉ
112 Sines and spikes x = Dc c DCT subdictionary Spike subdictionary
113 Denoising via l 1 -norm regularized least squares
114 Denoising via l 1 -norm regularized least squares Signal Estimate
115 Denoising via l 1 -norm regularized least squares Signal Estimate
116 Introduction Linear transforms Frequency representation Short-time Fourier transform (STFT) Wavelets Overcomplete sparse models Denoising The denoising problem Thresholding Synthesis model Analysis model
117 Denoising via l 1 -norm regularized least squares Analysis sparse model We would like to solve D T x minimize subject to is sparse D T x y x 0 Computationally intractable if D is overcomplete Instead, we solve ˆx = arg min y x R x 2 m 2 + λ AT x Significantly more challenging to solve than synthesis formulation 1
118 Total variation Images and some time series tend to be piecewise constant Equivalently: Sparse gradient (or derivative) The total variation of an image I is the l 1 -norm of the horizontal and vertical components of the gradient TV (I ) := x I 1 + y I 1 Equivalent to l 1 -norm regularization with an overcomplete analysis operator Denoising via TV regularization Y Î = arg min Ĩ 2 ) (Ĩ + λ TV Ĩ R n n F
119 Denoising via TV regularization Signal Data
120 Denoising via TV regularization Signal TV reg. (small λ)
121 Denoising via TV regularization Signal TV reg. (medium λ)
122 Denoising via TV regularization Signal TV reg. (large λ)
123 Denoising via TV regularization
124 Denoising via TV regularization
125 Small λ
126 Small λ
127 Medium λ
128 Medium λ
129 Large λ
130 Large λ
131 Denoising via TV regularization Original Noisy Estimate
Sparse linear models and denoising
Lecture notes 4 February 22, 2016 Sparse linear models and denoising 1 Introduction 1.1 Definition and motivation Finding representations of signals that allow to process them more effectively is a central
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
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 informationApplied Machine Learning for Biomedical Engineering. Enrico Grisan
Applied Machine Learning for Biomedical Engineering Enrico Grisan enrico.grisan@dei.unipd.it Data representation To find a representation that approximates elements of a signal class with a linear combination
More informationMultiresolution image processing
Multiresolution image processing Laplacian pyramids Some applications of Laplacian pyramids Discrete Wavelet Transform (DWT) Wavelet theory Wavelet image compression Bernd Girod: EE368 Digital Image Processing
More informationSparse regression. Optimization-Based Data Analysis. Carlos Fernandez-Granda
Sparse regression Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 3/28/2016 Regression Least-squares regression Example: Global warming Logistic
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 informationIntroduction to Wavelet. Based on A. Mukherjee s lecture notes
Introduction to Wavelet Based on A. Mukherjee s lecture notes Contents History of Wavelet Problems of Fourier Transform Uncertainty Principle The Short-time Fourier Transform Continuous Wavelet Transform
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationSparse analysis Lecture III: Dictionary geometry and greedy algorithms
Sparse analysis Lecture III: Dictionary geometry and greedy algorithms Anna C. Gilbert Department of Mathematics University of Michigan Intuition from ONB Key step in algorithm: r, ϕ j = x c i ϕ i, ϕ j
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationTowards a Mathematical Theory of Super-resolution
Towards a Mathematical Theory of Super-resolution Carlos Fernandez-Granda www.stanford.edu/~cfgranda/ Information Theory Forum, Information Systems Laboratory, Stanford 10/18/2013 Acknowledgements This
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
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 informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 12 Introduction to Wavelets Last Time Started with STFT Heisenberg Boxes Continue and move to wavelets Ham exam -- see Piazza post Please register at www.eastbayarc.org/form605.htm
More informationAn Introduction to Sparse Approximation
An Introduction to Sparse Approximation Anna C. Gilbert Department of Mathematics University of Michigan Basic image/signal/data compression: transform coding Approximate signals sparsely Compress images,
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationDigital Image Processing Lectures 15 & 16
Lectures 15 & 16, Professor Department of Electrical and Computer Engineering Colorado State University CWT and Multi-Resolution Signal Analysis Wavelet transform offers multi-resolution by allowing for
More informationSparsity in Underdetermined Systems
Sparsity in Underdetermined Systems Department of Statistics Stanford University August 19, 2005 Classical Linear Regression Problem X n y p n 1 > Given predictors and response, y Xβ ε = + ε N( 0, σ 2
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 informationMultiresolution analysis
Multiresolution analysis Analisi multirisoluzione G. Menegaz gloria.menegaz@univr.it The Fourier kingdom CTFT Continuous time signals + jωt F( ω) = f( t) e dt + f() t = F( ω) e jωt dt The amplitude F(ω),
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing () Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids Subband coding
More informationSparse & Redundant Signal Representation, and its Role in Image Processing
Sparse & Redundant Signal Representation, and its Role in Michael Elad The CS Department The Technion Israel Institute of technology Haifa 3000, Israel Wave 006 Wavelet and Applications Ecole Polytechnique
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationMachine Learning for Signal Processing Sparse and Overcomplete Representations
Machine Learning for Signal Processing Sparse and Overcomplete Representations Abelino Jimenez (slides from Bhiksha Raj and Sourish Chaudhuri) Oct 1, 217 1 So far Weights Data Basis Data Independent ICA
More informationWavelets in Pattern Recognition
Wavelets in Pattern Recognition Lecture Notes in Pattern Recognition by W.Dzwinel Uncertainty principle 1 Uncertainty principle Tiling 2 Windowed FT vs. WT Idea of mother wavelet 3 Scale and resolution
More informationRecent developments on sparse representation
Recent developments on sparse representation Zeng Tieyong Department of Mathematics, Hong Kong Baptist University Email: zeng@hkbu.edu.hk Hong Kong Baptist University Dec. 8, 2008 First Previous Next Last
More information2.3. Clustering or vector quantization 57
Multivariate Statistics non-negative matrix factorisation and sparse dictionary learning The PCA decomposition is by construction optimal solution to argmin A R n q,h R q p X AH 2 2 under constraint :
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationEE 381V: Large Scale Optimization Fall Lecture 24 April 11
EE 381V: Large Scale Optimization Fall 2012 Lecture 24 April 11 Lecturer: Caramanis & Sanghavi Scribe: Tao Huang 24.1 Review In past classes, we studied the problem of sparsity. Sparsity problem is that
More informationSGN Advanced Signal Processing Project bonus: Sparse model estimation
SGN 21006 Advanced Signal Processing Project bonus: Sparse model estimation Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 12 Sparse models Initial problem: solve
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationNiklas Grip, Department of Mathematics, Luleå University of Technology. Last update:
Some Essentials of Data Analysis with Wavelets Slides for the wavelet lectures of the course in data analysis at The Swedish National Graduate School of Space Technology Niklas Grip, Department of Mathematics,
More informationWavelets, Filter Banks and Multiresolution Signal Processing
Wavelets, Filter Banks and Multiresolution Signal Processing It is with logic that one proves; it is with intuition that one invents. Henri Poincaré Introduction - 1 A bit of history: from Fourier to Haar
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
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 informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationSparse Time-Frequency Transforms and Applications.
Sparse Time-Frequency Transforms and Applications. Bruno Torrésani http://www.cmi.univ-mrs.fr/~torresan LATP, Université de Provence, Marseille DAFx, Montreal, September 2006 B. Torrésani (LATP Marseille)
More informationMachine Learning for Signal Processing Sparse and Overcomplete Representations. Bhiksha Raj (slides from Sourish Chaudhuri) Oct 22, 2013
Machine Learning for Signal Processing Sparse and Overcomplete Representations Bhiksha Raj (slides from Sourish Chaudhuri) Oct 22, 2013 1 Key Topics in this Lecture Basics Component-based representations
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationStructured matrix factorizations. Example: Eigenfaces
Structured matrix factorizations Example: Eigenfaces An extremely large variety of interesting and important problems in machine learning can be formulated as: Given a matrix, find a matrix and a matrix
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 informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing (Wavelet Transforms) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids
More informationMIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications. Class 19: Data Representation by Design
MIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications Class 19: Data Representation by Design What is data representation? Let X be a data-space X M (M) F (M) X A data representation
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e
More informationBayesian Paradigm. Maximum A Posteriori Estimation
Bayesian Paradigm Maximum A Posteriori Estimation Simple acquisition model noise + degradation Constraint minimization or Equivalent formulation Constraint minimization Lagrangian (unconstraint minimization)
More informationIntroduction to Signal Processing
to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence
More informationSupport Detection in Super-resolution
Support Detection in Super-resolution Carlos Fernandez-Granda (Stanford University) 7/2/2013 SampTA 2013 Support Detection in Super-resolution C. Fernandez-Granda 1 / 25 Acknowledgements This work was
More informationReview: Learning Bimodal Structures in Audio-Visual Data
Review: Learning Bimodal Structures in Audio-Visual Data CSE 704 : Readings in Joint Visual, Lingual and Physical Models and Inference Algorithms Suren Kumar Vision and Perceptual Machines Lab 106 Davis
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 informationWavelets and Signal Processing
Wavelets and Signal Processing John E. Gilbert Mathematics in Science Lecture April 30, 2002. Publicity Mathematics In Science* A LECTURE SERIES FOR UNDERGRADUATES Wavelets Professor John Gilbert Mathematics
More informationECE472/572 - Lecture 13. Roadmap. Questions. Wavelets and Multiresolution Processing 11/15/11
ECE472/572 - Lecture 13 Wavelets and Multiresolution Processing 11/15/11 Reference: Wavelet Tutorial http://users.rowan.edu/~polikar/wavelets/wtpart1.html Roadmap Preprocessing low level Enhancement Restoration
More informationOptimization-based sparse recovery: Compressed sensing vs. super-resolution
Optimization-based sparse recovery: Compressed sensing vs. super-resolution Carlos Fernandez-Granda, Google Computational Photography and Intelligent Cameras, IPAM 2/5/2014 This work was supported by a
More informationDesigning Information Devices and Systems I Discussion 13B
EECS 6A Fall 7 Designing Information Devices and Systems I Discussion 3B. Orthogonal Matching Pursuit Lecture Orthogonal Matching Pursuit (OMP) algorithm: Inputs: A set of m songs, each of length n: S
More informationTutorial: Sparse Signal Processing Part 1: Sparse Signal Representation. Pier Luigi Dragotti Imperial College London
Tutorial: Sparse Signal Processing Part 1: Sparse Signal Representation Pier Luigi Dragotti Imperial College London Outline Part 1: Sparse Signal Representation ~90min Part 2: Sparse Sampling ~90min 2
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationMLISP: Machine Learning in Signal Processing Spring Lecture 10 May 11
MLISP: Machine Learning in Signal Processing Spring 2018 Lecture 10 May 11 Prof. Venia Morgenshtern Scribe: Mohamed Elshawi Illustrations: The elements of statistical learning, Hastie, Tibshirani, Friedman
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 informationTransforms and Orthogonal Bases
Orthogonal Bases Transforms and Orthogonal Bases We now turn back to linear algebra to understand transforms, which map signals between different domains Recall that signals can be interpreted as vectors
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Multi-dimensional signal processing Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione
More informationSparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images
Sparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images Alfredo Nava-Tudela ant@umd.edu John J. Benedetto Department of Mathematics jjb@umd.edu Abstract In this project we are
More informationTranslation-Invariant Shrinkage/Thresholding of Group Sparse. Signals
Translation-Invariant Shrinkage/Thresholding of Group Sparse Signals Po-Yu Chen and Ivan W. Selesnick Polytechnic Institute of New York University, 6 Metrotech Center, Brooklyn, NY 11201, USA Email: poyupaulchen@gmail.com,
More informationSparse Optimization Lecture: Basic Sparse Optimization Models
Sparse Optimization Lecture: Basic Sparse Optimization Models Instructor: Wotao Yin July 2013 online discussions on piazza.com Those who complete this lecture will know basic l 1, l 2,1, and nuclear-norm
More informationLecture 16: Multiresolution Image Analysis
Lecture 16: Multiresolution Image Analysis Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu November 9, 2004 Abstract Multiresolution analysis
More informationLinear Algebra and Dirac Notation, Pt. 1
Linear Algebra and Dirac Notation, Pt. 1 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 1 February 1, 2017 1 / 13
More informationMULTIRATE DIGITAL SIGNAL PROCESSING
MULTIRATE DIGITAL SIGNAL PROCESSING Signal processing can be enhanced by changing sampling rate: Up-sampling before D/A conversion in order to relax requirements of analog antialiasing filter. Cf. audio
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 informationRandom projections. 1 Introduction. 2 Dimensionality reduction. Lecture notes 5 February 29, 2016
Lecture notes 5 February 9, 016 1 Introduction Random projections Random projections are a useful tool in the analysis and processing of high-dimensional data. We will analyze two applications that use
More informationMorphological Diversity and Source Separation
Morphological Diversity and Source Separation J. Bobin, Y. Moudden, J.-L. Starck, and M. Elad Abstract This paper describes a new method for blind source separation, adapted to the case of sources having
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 27th, 21 From separation of variables, we move to linear algebra Roughly speaking, this is the study
More informationSparse Approximation and Variable Selection
Sparse Approximation and Variable Selection Lorenzo Rosasco 9.520 Class 07 February 26, 2007 About this class Goal To introduce the problem of variable selection, discuss its connection to sparse approximation
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationIntroduction to Sparsity in Signal Processing
1 Introduction to Sparsity in Signal Processing Ivan Selesnick Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 212 2 Under-determined linear equations Consider a system
More informationModule 4 MULTI- RESOLUTION ANALYSIS. Version 2 ECE IIT, Kharagpur
Module 4 MULTI- RESOLUTION ANALYSIS Lesson Theory of Wavelets Instructional Objectives At the end of this lesson, the students should be able to:. Explain the space-frequency localization problem in sinusoidal
More informationInvariant Scattering Convolution Networks
Invariant Scattering Convolution Networks Joan Bruna and Stephane Mallat Submitted to PAMI, Feb. 2012 Presented by Bo Chen Other important related papers: [1] S. Mallat, A Theory for Multiresolution Signal
More informationInverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France.
Inverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Overview of the course Introduction sparsity & data compression inverse problems
More informationCh. 15 Wavelet-Based Compression
Ch. 15 Wavelet-Based Compression 1 Origins and Applications The Wavelet Transform (WT) is a signal processing tool that is replacing the Fourier Transform (FT) in many (but not all!) applications. WT theory
More informationSparse Recovery Beyond Compressed Sensing
Sparse Recovery Beyond Compressed Sensing Carlos Fernandez-Granda www.cims.nyu.edu/~cfgranda Applied Math Colloquium, MIT 4/30/2018 Acknowledgements Project funded by NSF award DMS-1616340 Separable Nonlinear
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
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 informationImage Processing by the Curvelet Transform
Image Processing by the Curvelet Transform Jean Luc Starck Dapnia/SEDI SAP, CEA Saclay, France. jstarck@cea.fr http://jstarck.free.fr Collaborators: D.L. Donoho, Department of Statistics, Stanford E. Candès,
More information1. Fourier Transform (Continuous time) A finite energy signal is a signal f(t) for which. f(t) 2 dt < Scalar product: f(t)g(t)dt
1. Fourier Transform (Continuous time) 1.1. Signals with finite energy A finite energy signal is a signal f(t) for which Scalar product: f(t) 2 dt < f(t), g(t) = 1 2π f(t)g(t)dt The Hilbert space of all
More informationc 1998 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 2, No., pp. 33 6 c 998 Society for Industrial and Applied Mathematics ATOMIC DECOMPOSITION BY BASIS PURSUIT SCOTT SHAOBING CHEN, DAVID L. DONOHO, AND MICHAEL A. SAUNDERS Abstract.
More informationEE67I Multimedia Communication Systems
EE67I Multimedia Communication Systems Lecture 5: LOSSY COMPRESSION In these schemes, we tradeoff error for bitrate leading to distortion. Lossy compression represents a close approximation of an original
More informationWavelets. Introduction and Applications for Economic Time Series. Dag Björnberg. U.U.D.M. Project Report 2017:20
U.U.D.M. Project Report 2017:20 Wavelets Introduction and Applications for Economic Time Series Dag Björnberg Examensarbete i matematik, 15 hp Handledare: Rolf Larsson Examinator: Jörgen Östensson Juni
More informationMultiple Change Point Detection by Sparse Parameter Estimation
Multiple Change Point Detection by Sparse Parameter Estimation Department of Econometrics Fac. of Economics and Management University of Defence Brno, Czech Republic Dept. of Appl. Math. and Comp. Sci.
More informationWavelets For Computer Graphics
{f g} := f(x) g(x) dx A collection of linearly independent functions Ψ j spanning W j are called wavelets. i J(x) := 6 x +2 x + x + x Ψ j (x) := Ψ j (2 j x i) i =,..., 2 j Res. Avge. Detail Coef 4 [9 7
More informationA powerful method for feature extraction and compression of electronic nose responses
Sensors and Actuators B 105 (2005) 378 392 A powerful method for feature extraction and compression of electronic nose responses A. Leone a, C. Distante a,, N. Ancona b, K.C. Persaud c, E. Stella b, P.
More informationSparsity Models. Tong Zhang. Rutgers University. T. Zhang (Rutgers) Sparsity Models 1 / 28
Sparsity Models Tong Zhang Rutgers University T. Zhang (Rutgers) Sparsity Models 1 / 28 Topics Standard sparse regression model algorithms: convex relaxation and greedy algorithm sparse recovery analysis:
More informationDigital Image Processing Lectures 13 & 14
Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2013 Properties of KL Transform The KL transform has many desirable properties which makes
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 informationMultiscale Geometric Analysis: Thoughts and Applications (a summary)
Multiscale Geometric Analysis: Thoughts and Applications (a summary) Anestis Antoniadis, University Joseph Fourier Assimage 2005,Chamrousse, February 2005 Classical Multiscale Analysis Wavelets: Enormous
More informationA Review of Sparsity-Based Methods for Analysing Radar Returns from Helicopter Rotor Blades
A Review of Sparsity-Based Methods for Analysing Radar Returns from Helicopter Rotor Blades Ngoc Hung Nguyen 1, Hai-Tan Tran 2, Kutluyıl Doğançay 1 and Rocco Melino 2 1 University of South Australia 2
More information