1,910 research outputs found
On sparsity averaging
Recent developments in Carrillo et al. (2012) and Carrillo et al. (2013)
introduced a novel regularization method for compressive imaging in the context
of compressed sensing with coherent redundant dictionaries. The approach relies
on the observation that natural images exhibit strong average sparsity over
multiple coherent frames. The associated reconstruction algorithm, based on an
analysis prior and a reweighted scheme, is dubbed Sparsity Averaging
Reweighted Analysis (SARA). We review these advances and extend associated
simulations establishing the superiority of SARA to regularization methods
based on sparsity in a single frame, for a generic spread spectrum acquisition
and for a Fourier acquisition of particular interest in radio astronomy.Comment: 4 pages, 3 figures, Proceedings of 10th International Conference on
Sampling Theory and Applications (SampTA), Code available at
https://github.com/basp-group/sopt, Full journal letter available at
http://arxiv.org/abs/arXiv:1208.233
PURIFY: a new algorithmic framework for next-generation radio-interferometric imaging
In recent works, compressed sensing (CS) and convex opti- mization techniques have been applied to radio-interferometric imaging showing the potential to outperform state-of-the-art imaging algorithms in the field. We review our latest contributions [1, 2, 3], which leverage the versatility of convex optimization to both handle realistic continuous visibilities and offer a highly parallelizable structure paving the way to significant acceleration of the reconstruction and high-dimensional data scalability. The new algorithmic structure promoted in a new software PURIFY (beta version) relies on the simultaneous-direction method of multipliers (SDMM). The performance of various sparsity priors is evaluated through simulations in the continuous visibility setting, confirming the superiority of our recent average sparsity approach SARA
Cygnus A super-resolved via convex optimisation from VLA data
We leverage the Sparsity Averaging Reweighted Analysis (SARA) approach for
interferometric imaging, that is based on convex optimisation, for the
super-resolution of Cyg A from observations at the frequencies 8.422GHz and
6.678GHz with the Karl G. Jansky Very Large Array (VLA). The associated average
sparsity and positivity priors enable image reconstruction beyond instrumental
resolution. An adaptive Preconditioned Primal-Dual algorithmic structure is
developed for imaging in the presence of unknown noise levels and calibration
errors. We demonstrate the superior performance of the algorithm with respect
to the conventional CLEAN-based methods, reflected in super-resolved images
with high fidelity. The high resolution features of the recovered images are
validated by referring to maps of Cyg A at higher frequencies, more precisely
17.324GHz and 14.252GHz. We also confirm the recent discovery of a radio
transient in Cyg A, revealed in the recovered images of the investigated data
sets. Our matlab code is available online on GitHub.Comment: 14 pages, 7 figures (3/7 animated figures), accepted for publication
in MNRA
Compressed Sensing Applied to Weather Radar
We propose an innovative meteorological radar, which uses reduced number of
spatiotemporal samples without compromising the accuracy of target information.
Our approach extends recent research on compressed sensing (CS) for radar
remote sensing of hard point scatterers to volumetric targets. The previously
published CS-based radar techniques are not applicable for sampling weather
since the precipitation echoes lack sparsity in both range-time and Doppler
domains. We propose an alternative approach by adopting the latest advances in
matrix completion algorithms to demonstrate the sparse sensing of weather
echoes. We use Iowa X-band Polarimetric (XPOL) radar data to test and
illustrate our algorithms.Comment: 4 pages, 5 figrue
Ab initio compressive phase retrieval
Any object on earth has two fundamental properties: it is finite, and it is
made of atoms. Structural information about an object can be obtained from
diffraction amplitude measurements that account for either one of these traits.
Nyquist-sampling of the Fourier amplitudes is sufficient to image single
particles of finite size at any resolution. Atomic resolution data is routinely
used to image molecules replicated in a crystal structure. Here we report an
algorithm that requires neither information, but uses the fact that an image of
a natural object is compressible. Intended applications include tomographic
diffractive imaging, crystallography, powder diffraction, small angle x-ray
scattering and random Fourier amplitude measurements.Comment: 7 pages, 4 figures, presented at the XXI IUCr Congress, Aug. 2008,
Osaka Japa
- …