Abstract:We propose a probabilistic model for sparse signal reconstruction and develop several novel algorithms for computing the maximum likelihood (ML) parameter estimates under this model. The measurements follow an underdetermined linear model where the regression-coefficient vector is the sum of an unknown deterministic sparse signal component and a zero-mean white Gaussian component with an unknown variance. Our reconstruction schemes are based on an expectation-conditional maximization either (ECME) iteration that aims at maximizing the likelihood function with respect to the unknown parameters for a given signal sparsity level. Compared with the existing iterative hard thresholding (IHT) method, the ECME algorithm contains an additional multiplicative term and guarantees monotonic convergence for a wide range of sensing (regression) matrices. We propose a double overrelaxation (DORE) thresholding scheme for accelerating the ECME iteration. We prove that, under certain mild conditions, the ECME and DORE iterations converge to local maxima of the likelihood function. The ECME and DORE iterations can be implemented exactly in small-scale applications and for the important class of large-scale sensing operators with orthonormal rows used e. g., partial fast Fourier transform (FFT). If the signal sparsity level is unknown, we introduce an unconstrained sparsity selection (USS) criterion and a tuning-free automatic double overrelaxation (ADORE) thresholding method that employs USS to estimate the sparsity level. We compare the proposed and existing sparse signal reconstruction methods via one-dimensional simulation and two-dimensional image reconstruction experiments using simulated and real X-ray CT data.

Estimating the Signal Reconstruction Error from Threshold-Based Sampling Without Knowing the Original Signal

Estimation-free Spatial-Domain Image Reconstruction of Structured Illumination Microscopy

Sparse Bayesian Learning Approach for Discrete Signal Reconstruction

Reducing Identification Error Using Over-Sampling Technique

Random Sampling Using Shannon Interpolation and Poisson Summation Formulae

ECME Thresholding Methods for Sparse Signal Reconstruction

Approximate reconstruction of bandlimited functions for the integrate and fire sampler

Signal Number Estimation Based on Reconstructed Matrix

Signal reconstruction using determinantal sampling

Sparsity-based Reconstruction Method for Signals with Finite Rate of Innovation

Error Analysis of Reconstruction from Linear Canonical Transform-Based Sampling

Minimum Rate Sampling and Spectrum Blind Reconstruction in Random Equivalent Sampling

A Robust Reconstruction Algorithm for Sparse Signals

Signal Reconstruction from Nonideal Samples in Fractional Fourier Transform Domain

Sparse Signal Reconstruction Based on Multiparameter Approximation Function with Smoothedl0Norm

Randomized Iterative Hard Thresholding for Sparse Approximations

A Reconstruction Method of Generalized Sampling Based on Generalized Inverse

Xampling at the Rate of Innovation

Sparse Signal Reconstruction from Quantized Noisy Measurements Via GEM Hard Thresholding

Sparse Signal Reconstruction Via ECME Hard Thresholding.

Learning-Based Reconstruction of FRI Signals