Joint Inverse Problems for Signal Reconstruction via Dictionary Splitting

Abstract

Sparse signal recovery from limited and/or degraded samples is fundamental to many applications, such as medical imaging, remote sensing, astronomical and seismic imaging. Discrete wavelet transform (DWT) has been commonly used for sparse representation of signals; nevertheless, due to its shift-variant nature, pseudo-Gibbs artifacts are present in the recovered signals. Using the redundant shift-invariant wavelet transform (SWT) is the ideal solution to obtain shift invariance; however, high redundancy factor of SWT limits its application in practical settings. We propose a dictionary splitting approach for sparse recovery from incomplete data, which leverages the ideas of cycle spinning in combination with Bregman splitting. The proposed method significantly improves the conventional signal reconstruction with DWT, offers the advantages of SWT, and overcomes high redundancy factor of SWT. We solve parallel sparse recovery problems with orthogonal dictionaries (DWT and its permuted versions), while we impose consistency between the results by updating the recovered image at each iteration. Our experiments demonstrate that few shifts are sufficient to achieve reconstruction accuracy as high as recovery with SWT, and significantly reduces its computational cost and redundancy factor.