2,137 research outputs found
Signal Recovery in Perturbed Fourier Compressed Sensing
In many applications in compressed sensing, the measurement matrix is a
Fourier matrix, i.e., it measures the Fourier transform of the underlying
signal at some specified `base' frequencies , where is the
number of measurements. However due to system calibration errors, the system
may measure the Fourier transform at frequencies
that are different from the base frequencies and where
are unknown. Ignoring perturbations of this nature can lead to major errors in
signal recovery. In this paper, we present a simple but effective alternating
minimization algorithm to recover the perturbations in the frequencies \emph{in
situ} with the signal, which we assume is sparse or compressible in some known
basis. In many cases, the perturbations can be expressed
in terms of a small number of unique parameters . We demonstrate that
in such cases, the method leads to excellent quality results that are several
times better than baseline algorithms (which are based on existing off-grid
methods in the recent literature on direction of arrival (DOA) estimation,
modified to suit the computational problem in this paper). Our results are also
robust to noise in the measurement values. We also provide theoretical results
for (1) the convergence of our algorithm, and (2) the uniqueness of its
solution under some restrictions.Comment: New theortical results about uniqueness and convergence now included.
More challenging experiments now include
PEAR: PEriodic And fixed Rank separation for fast fMRI
In functional MRI (fMRI), faster acquisition via undersampling of data can
improve the spatial-temporal resolution trade-off and increase statistical
robustness through increased degrees-of-freedom. High quality reconstruction of
fMRI data from undersampled measurements requires proper modeling of the data.
We present an fMRI reconstruction approach based on modeling the fMRI signal as
a sum of periodic and fixed rank components, for improved reconstruction from
undersampled measurements. We decompose the fMRI signal into a component which
a has fixed rank and a component consisting of a sum of periodic signals which
is sparse in the temporal Fourier domain. Data reconstruction is performed by
solving a constrained problem that enforces a fixed, moderate rank on one of
the components, and a limited number of temporal frequencies on the other. Our
approach is coined PEAR - PEriodic And fixed Rank separation for fast fMRI.
Experimental results include purely synthetic simulation, a simulation with
real timecourses and retrospective undersampling of a real fMRI dataset.
Evaluation was performed both quantitatively and visually versus ground truth,
comparing PEAR to two additional recent methods for fMRI reconstruction from
undersampled measurements. Results demonstrate PEAR's improvement in estimating
the timecourses and activation maps versus the methods compared against at
acceleration ratios of R=8,16 (for simulated data) and R=6.66,10 (for real
data). PEAR results in reconstruction with higher fidelity than when using a
fixed-rank based model or a conventional Low-rank+Sparse algorithm. We have
shown that splitting the functional information between the components leads to
better modeling of fMRI, over state-of-the-art methods
- …