Skip to content Skip to navigation Skip to collection information

Connexions

You are here: Home » Content » An Introduction to Compressive Sensing » Introduction to compressive sensing

Navigation

Recently Viewed

This feature requires Javascript to be enabled.
Reuse / Edit
x

Collection:

Module:

Add to a lens
x

Add collection to:

Add module to:

Add to Favorites
x

Add collection to:

Add module to:

 

Introduction to compressive sensing

Module by: Mark A. Davenport, Marco F. Duarte, Chinmay Hegde, Richard Baraniuk. E-mail the authors

Summary: Introduction to compressive sensing. This course introduces the basic concepts in compressive sensing. We overview the concepts of sparsity, compressibility, and transform coding. We then review applications of sparsity in several signal processing problems such as sparse recovery, model selection, data coding, and error correction. We overview the key results in these fields, focusing primarily on both theory and algorithms for sparse recovery. We also discuss applications of compressive sensing in communications, biosensing, medical imaging, and sensor networks.

We are in the midst of a digital revolution that is driving the development and deployment of new kinds of sensing systems with ever-increasing fidelity and resolution. The theoretical foundation of this revolution is the pioneering work of Kotelnikov, Nyquist, Shannon, and Whittaker on sampling continuous-time band-limited signals [14], [15], [18], [21]. Their results demonstrate that signals, images, videos, and other data can be exactly recovered from a set of uniformly spaced samples taken at the so-called Nyquist rate of twice the highest frequency present in the signal of interest. Capitalizing on this discovery, much of signal processing has moved from the analog to the digital domain and ridden the wave of Moore's law. Digitization has enabled the creation of sensing and processing systems that are more robust, flexible, cheaper and, consequently, more widely-used than their analog counterparts.

As a result of this success, the amount of data generated by sensing systems has grown from a trickle to a torrent. Unfortunately, in many important and emerging applications, the resulting Nyquist rate is so high that we end up with far too many samples. Alternatively, it may simply be too costly, or even physically impossible, to build devices capable of acquiring samples at the necessary rate. Thus, despite extraordinary advances in computational power, the acquisition and processing of signals in application areas such as imaging, video, medical imaging, remote surveillance, spectroscopy, and genomic data analysis continues to pose a tremendous challenge.

To address the logistical and computational challenges involved in dealing with such high-dimensional data, we often depend on compression, which aims at finding the most concise representation of a signal that is able to achieve a target level of acceptable distortion. One of the most popular techniques for signal compression is known as transform coding, and typically relies on finding a basis or frame that provides sparse or compressible representations for signals in a class of interest. By a sparse representation, we mean that for a signal of length NN, we can represent it with KNKN nonzero coefficients; by a compressible representation, we mean that the signal is well-approximated by a signal with only KK nonzero coefficients. Both sparse and compressible signals can be represented with high fidelity by preserving only the values and locations of the largest coefficients of the signal. This process is called sparse approximation, and forms the foundation of transform coding schemes that exploit signal sparsity and compressibility, including the JPEG, JPEG2000, MPEG, and MP3 standards.

Leveraging the concept of transform coding, compressive sensing (CS) has emerged as a new framework for signal acquisition and sensor design. CS enables a potentially large reduction in the sampling and computation costs for sensing signals that have a sparse or compressible representation. The Nyquist-Shannon sampling theorem states that a certain minimum number of samples is required in order to perfectly capture an arbitrary bandlimited signal, but when the signal is sparse in a known basis we can vastly reduce the number of measurements that need to be stored. Consequently, when sensing sparse signals we might be able to do better than suggested by classical results. This is the fundamental idea behind CS: rather than first sampling at a high rate and then compressing the sampled data, we would like to find ways to directly sense the data in a compressed form — i.e., at a lower sampling rate. The field of CS grew out of the work of Emmanuel Candès, Justin Romberg, and Terence Tao and of David Donoho, who showed that a finite-dimensional signal having a sparse or compressible representation can be recovered from a small set of linear, nonadaptive measurements [1], [5], [10]. The design of these measurement schemes and their extensions to practical data models and acquisition schemes are one of the most central challenges in the field of CS.

Although this idea has only recently gained significant attraction in the signal processing community, there have been hints in this direction dating back as far as the eighteenth century. In 1795, Prony proposed an algorithm for the estimation of the parameters associated with a small number of complex exponentials sampled in the presence of noise [16]. The next theoretical leap came in the early 1900's, when Carathéodory showed that a positive linear combination of any KK sinusoids is uniquely determined by its value at t=0t=0 and at any other 2K2K points in time [3], [4]. This represents far fewer samples than the number of Nyquist-rate samples when KK is small and the range of possible frequencies is large. In the 1990's, this work was generalized by George, Gorodnitsky, and Rao, who studied sparsity in the context of biomagnetic imaging and other contexts [13], [17], and by Bressler and Feng, who proposed a sampling scheme for acquiring certain classes of signals consisting of KK components with nonzero bandwidth (as opposed to pure sinusoids) [11], [12]. In the early 2000's Vetterli, Marziliano, and Blu proposed a sampling scheme for non-bandlimited signals that are governed by only KK parameters, showing that these signals can be sampled and recovered from just 2K2K samples [20].

A related problem focuses on recovery of a signal from partial observation of its Fourier transform. Beurling proposed a method for extrapolating these observations to determine the entire Fourier transform [2]. One can show that if the signal consists of a finite number of impulses, then Beurling's approach will correctly recover the entire Fourier transform (of this non-bandlimited signal) from any sufficiently large piece of its Fourier transform. His approach — to find the signal with smallest 11 norm among all signals agreeing with the acquired Fourier measurements — bears a remarkable resemblance to some of the algorithms used in CS.

More recently, Candès, Romberg, Tao [5], [6], [7], [8], [9], and Donoho [10] showed that a signal having a sparse representation can be recovered exactly from a small set of linear, nonadaptive measurements. This result suggests that it may be possible to sense sparse signals by taking far fewer measurements, hence the name compressive sensing. Note, however, that CS differs from classical sampling in two important respects. First, rather than sampling the signal at specific points in time, CS systems typically acquire measurements in the form of inner products between the signal and more general test functions. We will see throughout this course that randomness often plays a key role in the design of these test functions. Second, the two frameworks differ in the manner in which they deal with signal recovery, i.e., the problem of recovering the original signal from the compressive measurements. In the Nyquist-Shannon framework, signal recovery is achieved through cardinal sine (sinc) interpolation — a linear process that requires little computation and has a simple interpretation.

CS has already had notable impact on several applications. One example is medical imaging, where it has enabled speedups by a factor of seven in pediatric MRI while preserving diagnostic quality [19]. Moreover, the broad applicability of this framework has inspired research that extends the CS framework by proposing practical implementations for numerous applications, including sub-Nyquist analog-to-digital converters (ADCs), compressive imaging architectures, and compressive sensor networks.

This course introduces the basic concepts in compressive sensing. We overview the concepts of sparsity, compressibility, and transform coding. We overview the key results in the field, beginning by focusing primarily on the theory of sensing matrix design, 11-minimization, and alternative algorithms for sparse recovery. We then review applications of sparsity in several signal processing problems such as sparse regression and model selection, error correction, group testing, and compressive inference. We also discuss applications of compressive sensing in analog-to-digital conversion, biosensing, conventional and hyperspectral imaging, medical imaging, and sensor networks.

Acknowledgments

The authors would like to thank Ewout van den Berg, Yonina Eldar, Piotr Indyk, Gitta Kutyniok, and Yaniv Plan for their feedback regarding some portions of this course which now also appear in the introductory chapter of Compressed Sensing: Theory and Applications, Cambridge University Press, 2011.

References

  1. Baraniuk, R. (2007). Compressive sensing. IEEE Signal Processing Mag., 24(4), 118–120, 124.
  2. Beurling, A. (1938). Sur les intégrales de Fourier absolument convergentes et leur application à une transformation fonctionelle. In Proc. Scandinavian Math. Congress. Helsinki, Finland
  3. Carathéodory, C. (1907). Über den Variabilitätsbereich der Koeffizienten von Potenzreihen, die gegebene Werte nicht annehmen. Math. Ann., 64, 95–115.
  4. Carathéodory, C. (1911). Über den Variabilitätsbereich der Fourierschen Konstanten von positiven harmonischen Funktionen. Rend. Circ. Mat. Palermo, 32, 193–217.
  5. Candès, E. (2006, Aug.). Compressive sampling. In Proc. Int. Congress of Math. Madrid, Spain
  6. Candès, E. and Romberg, J. (2006). Quantitative robust uncertainty principles and optimally sparse decompositions. Found. Comput. Math., 6(2), 227–254.
  7. Candès, E. and Romberg, J. and Tao, T. (2006). Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inform. Theory, 52(2), 489–509.
  8. Candès, E. and Romberg, J. and Tao, T. (2006). Stable signal recovery from incomplete and inaccurate measurements. Comm. Pure Appl. Math., 59(8), 1207–1223.
  9. Candès, E. and Tao, T. (2006). Near Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? IEEE Trans. Inform. Theory, 52(12), 5406–5425.
  10. Donoho, D. (2006). Compressed sensing. IEEE Trans. Inform. Theory, 52(4), 1289–1306.
  11. Feng, P. and Bresler, Y. (1996, May). Spectrum-blind minimum-rate sampling and reconstruction of multiband signals. In Proc. IEEE Int. Conf. Acoust., Speech, and Signal Processing (ICASSP). Atlanta, GA
  12. Feng, P. (1997, Mar.). Universal spectrum blind minimum rate sampling and reconstruction of multiband signals. Ph. D. Thesis. University of Illinois at Urbana-Champaign.
  13. Gorodnitsky, I. and Rao, B. and George, J. (1992, Oct.). Source Localization in Magnetoencephalagraphy using an Iterative Weighted Minimum Norm Algorithm. In Proc. Asilomar Conf. Signals, Systems, and Computers. Pacific Grove, CA
  14. Kotelnikov, V. (1933). On the carrying capacity of the ether and wire in telecommunications. In Izd. Red. Upr. Svyazi RKKA. Moscow, Russia
  15. Nyquist, H. (1928). Certain topics in telegraph transmission theory. Trans. AIEE, 47, 617–644.
  16. Prony, R. (1795). Essai expérimental et analytique sur les lois de la Dilatabilité des fluides élastiques et sur celles de la Force expansive de la vapeur de l'eau et de la vapeur de l'alkool, à différentes températures. [R. Prony is Gaspard Riche, baron de Prony]. J. de l'École Polytechnique, Floréal et Prairial III, 1(2), 24–76.
  17. Rao, B. (1998, May). Signal Processing with the Sparseness Constraint. In Proc. IEEE Int. Conf. Acoust., Speech, and Signal Processing (ICASSP). Seattle, WA
  18. Shannon, C. (1949). Communication in the presence of noise. Proc. Institute of Radio Engineers, 37(1), 10–21.
  19. Vasanawala, S. and Alley, M. and Barth, R. and Hargreaves, B. and Pauly, J. and Lustig, M. (2009, Apr.). Faster Pediatric MRI Via Compressed Sensing. In Proc. Annual Meeting Soc. Pediatric Radiology (SPR). Carlsbad, CA
  20. Vetterli, M. and Marziliano, P. and Blu, T. (2002). Sampling signals with finite rate of innovation. IEEE Trans. Signal Processing, 50(6), 1417–1428.
  21. Whittaker, E. (1915). On the Functions Which are Represented by the Expansions of the Interpolation Theory. Proc. Royal Soc. Edinburgh, Sec. A, 35, 181–194.

Collection Navigation

Content actions

Download module as:

Add:

Collection to:

My Favorites (?)

'My Favorites' is a special kind of lens which you can use to bookmark modules and collections. 'My Favorites' can only be seen by you, and collections saved in 'My Favorites' can remember the last module you were on. You need an account to use 'My Favorites'.

| A lens I own (?)

Definition of a lens

Lenses

A lens is a custom view of the content in the repository. You can think of it as a fancy kind of list that will let you see content through the eyes of organizations and people you trust.

What is in a lens?

Lens makers point to materials (modules and collections), creating a guide that includes their own comments and descriptive tags about the content.

Who can create a lens?

Any individual member, a community, or a respected organization.

What are tags? tag icon

Tags are descriptors added by lens makers to help label content, attaching a vocabulary that is meaningful in the context of the lens.

| External bookmarks

Module to:

My Favorites (?)

'My Favorites' is a special kind of lens which you can use to bookmark modules and collections. 'My Favorites' can only be seen by you, and collections saved in 'My Favorites' can remember the last module you were on. You need an account to use 'My Favorites'.

| A lens I own (?)

Definition of a lens

Lenses

A lens is a custom view of the content in the repository. You can think of it as a fancy kind of list that will let you see content through the eyes of organizations and people you trust.

What is in a lens?

Lens makers point to materials (modules and collections), creating a guide that includes their own comments and descriptive tags about the content.

Who can create a lens?

Any individual member, a community, or a respected organization.

What are tags? tag icon

Tags are descriptors added by lens makers to help label content, attaching a vocabulary that is meaningful in the context of the lens.

| External bookmarks

Reuse / Edit:

Reuse or edit collection (?)

Check out and edit

If you have permission to edit this content, using the "Reuse / Edit" action will allow you to check the content out into your Personal Workspace or a shared Workgroup and then make your edits.

Derive a copy

If you don't have permission to edit the content, you can still use "Reuse / Edit" to adapt the content by creating a derived copy of it and then editing and publishing the copy.

| Reuse or edit module (?)

Check out and edit

If you have permission to edit this content, using the "Reuse / Edit" action will allow you to check the content out into your Personal Workspace or a shared Workgroup and then make your edits.

Derive a copy

If you don't have permission to edit the content, you can still use "Reuse / Edit" to adapt the content by creating a derived copy of it and then editing and publishing the copy.