Today I write about sparse recovery of “multidimensional” signals. With “multidimensional” I mean something like this: A one-dimensional signal is a vector while a two-dimensional signal is a matrix . Similarly, -dimensional signal is . Of course, images as two-dimensional signals come time mind. Moreover, movies are three-dimensional, a hyperspectral 2D image (which has a whole spectrum attached to any pixel) is also three-dimensional, and time-dependent volume-data is four-dimensional.

Multidimensional data is often a challenge due the large amount of data. While the size of the signals is usually not the problem it is more the size of the *measurement matrices*. In the context of compressed sensing or sparse recovery the signal is measured with a linear operator, i.e. one applies a number of linear functionals to the signal. In the -dimensional case this can be encoded as a matrix and this is where the trouble with the data comes in: If you have megapixel image (which is still quite small) the matrix has a million of columns and if you have a dense matrix, storage becomes an issue.

One approach (which is indeed quite old) to tackle this problem is, to consider special measurement matrices: If the signal has a sparse structure in every slice, i.e. every vector is of the form where we fix all but the -th component, then the *Kronecker product* of measurement matrices for each dimension is the right thing.

The Kronecker product of two matrices and is the matrix

This has a lot to do with the tensor product and you should read the Wikipedia entry. Moreover, it is numerically advantageous not to build the Kronecker product of dense matrices if you only want to apply it to a given signal. To see this, we introduce the vectorization operator which takes a matrix and stacks its columns into a tall column vector. For matrices and (of fitting sizes) it holds that

So, multiplying from the left and from the right gives the application of the Kronecker product.

The use of Kronecker products in numerical linear algebra is fairly old (for example they are helpful for multidimensional finite difference schemes where you can build Kronecker products of sparse difference operators in respective dimensions). Recently, they have been discovered for compressed sensing and sparse recovery in these two papers: Sparse solutions to underdetermined Kronecker product systems by Sadegh Jokar and Volker Mehrmann and the more recent Kronecker Compressed Sensing by Marco Duarte and Rich Baraniuk.

From these papers you can extract some interestingly simple and nice theorems:

Theorem 1For matrices with restricted isometry constant of order it holds that the restricted isometry constant of the Kronecker product fulfills

Basically, the RIP constant of a Kronecker product is not better than the worst one but still not too large.

Theorem 2For matrices with columns normalized to one, it hold that the spark of their Kronecker product fulfills

Theorem 3For matrices with columns normalized to one, it hold that the mutual coherence of their Kronecker product fulfills

February 26, 2012 at 11:42 am

Earlier applications for compressed sensing

:

Y. Rivenson and A. Stern, ”Compressed imaging with separable sensing operator”, IEEE Signal Processing Letters, vol. 16( 6), 449-452, (June 2009).

Y. Rivenson and A. Stern, “An Efficient Method for Multi-Dimensional Compressive Imaging,” in Computational Optical Sensing and Imaging, OSA Technical Digest (CD) (Optical Society of America, 2009), paper CTuA4.

February 26, 2012 at 6:44 pm

Thanks Adrian for the additional references!