introduction to compressive sensing
DESCRIPTION
Introduction to Compressive Sensing. Richard Baraniuk , Compressive sensing . IEEE Signal Processing Magazine, 24(4), pp. 118-121, July 2007) Emmanuel Candès and Michael Wakin , An introduction to compressive sampling . IEEE Signal Processing Magazine, 25(2), pp. 21 - 30, March 2008 - PowerPoint PPT PresentationTRANSCRIPT
Introduction to Compressive Sensing
Richard Baraniuk, Compressive sensing . IEEE Signal Processing Magazine, 24(4), pp. 118-121, July 2007)
Emmanuel Candès and Michael Wakin, An introduction to compressive sampling . IEEE Signal Processing Magazine, 25(2), pp. 21 - 30, March 2008
A course on compressive sensing, http://w3.impa.br/~aschulz/CS/course.html
Outline
• Introduction to compressive sensing (CS)– First CS theory– Concepts and applications– Theory• Compression• Reconstruction
Introduction• Compressive sensing
– Compressed sensing– Compressive sampling
• First CS theory– E. Cand`es, J. Romberg, and T. Tao, “Robust uncertainty principles: Exact signal
reconstruction from highly incomplete frequency information,” IEEE Trans. Inform. Theory, vol. 52, no. 2, pp. 489–509, Feb. 2006.
Cand`es Romberg Tao
Compressive Sensing: concept and applications
Compression/Reconstruction
CS samplingX RNx1
RMxN
Measurement matrix
yRMx1 Quantizationhuman coding
Transmit
OptimizationInverse
transform (e.g., IDCT)
s X’InverseQuantization
human coding
y’
syss ' subjet to ||||min 0 sX '
: transform basis (e.g., DCT basis)
CS Reconstruction
Theory and Core Technologycompression
• K-sparse– most of the energy is at low frequencies– K non-zero wavelet (DCT) coefficients
sX
Compression
y XMeasurement matrix
Compression
y transform basis
scoefficient
X
Compression
y
transform basis
scoefficient
Reconstruction
Reconstruction: optimization
syss
osubject t ||||min 0
syss
osubject t ||||min 1
NP-hard problem
Linear programming [1][2]Orthogonal matching pursuit (OMP)
syss
osubject t ||||min 2
Minimum energy ≠ k-sparse
(1)
(2)
(3)
(4) Greedy algorithm [3]
Compressive sensing: significant parameters
1. What measurement matrix should we use?2. How many measurements? (M=?)3. K-sparse?
Measurement Matrix Incoherence
(1) Correlation between and ],1[),( )2( n
Examples
= noiselet, = Haar wavelet (,)=2= noiselet, = Daubechies D4 (,)=2.2= noiselet, = Daubechies D8 (,)=2.9– Noiselets are also maximally incoherent with spikes and
incoherent with the Fourier basis
= White noise (random Gaussian)
Restricted Isometry Property (RIP)preserving length
• RIP: For each integer k = 1, 2, …, define the isometry constant k of a matrix A as the smallest number such that
(1) A approximately preserves the Euclidean length of k-sparse signals
(2) Imply that k-sparse vectors cannot be in the nullspace of A
(3) All subsets of s columns taken from A are in fact nearly orthogonal– To design a sensing matrix , so that any subset of columns of size k be
approximately orthogonal.
)1(||||||||
)1( 2
2
2
2k
l
lk x
Ax
How many measurements ?
)log(~log),(2
nkOMnkCM
Single-Pixel CS Camera[Baraniuk and Kelly, et al.]
On the Interplay Between Routing and SignalRepresentation for Compressive Sensing in
Wireless Sensor Networks
G. Quer, R. Masiero, D. Munaretto, M. Rossi, J. Widmer and M. Zorzi
University of Padova, Italy.DoCoMo Euro-Labs, Germany
Information Theory and Applications Workshop (ITA 2009)
Network Scenario Settingx11 x12 x13 x14
x21 x22 x23 x24
… … .. ..… .. … …
X
Example of the considered multi-hop topology.
Irregular network setting [4](1) Graph wavelet(2) Diffusion wavelet
Measurement matrix Built on routing path
Routing path mm xy 111
random }1,1{1
mmm xyy 2212
mmm xyy 3323
mmm xyy 4434
…………………………………………
……………… ……………………
mlm
y ,........ ml321
nx
xxx
.3
2
1
Measurement matrix
• R1: is built according to routing protocol, – randomly selected from {+1, -1}
• R2: is built according to routing protocol – randomly selected from (0, 1]
• R3: has all coefficients in randomly selected from {+1, -1}• R4: has all coefficients in randomly selected from(0, 1]
Transform basis
• T1: DCT• T2: Haar Wavelet• T3: Horizontal difference• T4: Vertical difference + Horizontal difference
Degree of sparsity
DCT
Haar
H-diff
VH-diff
Incoherence
DCT
Haar
H-diffVH-diff
Performance Comparison
• Random sampling (RS)– each node sends its data with probability P = M/N,
the data packets are not processed at internal nodes but simply forwarded.
• RS-CS– the data values are combined
with that of any other node encountered along the path.
Routing path mm xy 111
random }1,1{1
mmm xyy 2212
mmm xyy 3323
mmm xyy 4434
Reconstruction Error
Reconstruction Errorpre-distribution for T3 and T4 [5]
Research issues when applying CS in Sensor Networks
1. How to construct measurement matrix – Incoherent with transform basis – Distributed – M=?
2. How to choose transformation basis – Sparsity– Incoherent with measurement matrix
3. Irregular sensor deployment– Graph wavelet– Diffusion wavelet
References[1] Bloomfield, P., Steiger, W., Least Absolute Deviations: Theory, Applications, and Algorithms. Progr. Probab. Statist. 6, Birkhäuser, Boston, MA, 1983.[2] Chen, S. S., Donoho, D. L., Saunders, M. A, Atomic decomposition by basis pursuit. SIAM J. Sci. Comput. 20 (1999), 33–61.[3] J. Tropp and A. C. Gilbert, “Signal recovery from partial information via orthogonal matching pursuit,” Apr. 2005, Preprint.[4] J. Haupt, W.U. Bajwa, M. Rabbat, and R. Nowak, “Compressed
sensing for networked data,” IEEE Signal Processing Mag., vol. 25, no. 2, pp. 92-101, Mar. 2008.
[5] M. Rabbat, J. Haupt, A. Singh, and R. Novak, “Decentralized Compression and Predistribution via Randomized Gossiping,” in IPSN, 2006.