Researcher profile

André Ferrari

André Ferrari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2016arXiv

Distributed multi-frequency image reconstruction for radio-interferometry

The advent of enhanced technologies in radio interferometry and the perspective of the SKA telescope bring new challenges in image reconstruction. One of these challenges is the spatio-spectral reconstruction of large (Terabytes) data cubes with high fidelity. This contribution proposes an alternative implementation of one such 3D prototype algorithm, MUFFIN (MUlti-Frequency image reconstruction For radio INterferometry), which combines spatial and spectral analysis priors. Using a recently proposed primal dual algorithm, this new version of MUFFIN allows a parallel implementation where computationally intensive steps are split by spectral channels. This parallelization allows to implement computationally demanding translation invariant wavelet transforms (IUWT), as opposed to the union of bases used previously. This alternative implementation is important as it opens the possibility of comparing these efficient dictionaries, and others, in spatio-spectral reconstruction. Numerical results show that the IUWT-based version can be successfully implemented at large scale with performances comparable to union of bases.

preprint2016arXiv

Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach

We describe a "spatio-spectral" deconvolution algorithm for wide-band imaging in radio interferometry. In contrast with the existing multi-frequency reconstruction algorithms, the proposed method does not rely on a model of the sky-brightness spectral distribution. This non-parametric approach can be of particular interest for the new generation of low frequency radiotelescopes. The proposed solution formalizes the reconstruction problem as a convex optimization problem with spatial and spectral regularizations. The efficiency of this approach has been already proven for narrow-band image reconstruction and the present contribution can be considered as its extension to the multi-frequency case. Because the number of frequency bands multiplies the size of the inverse problem, particular attention is devoted to the derivation of an iterative large scale optimization algorithm. It is shown that the main computational bottleneck of the approach, which lies in the resolution of a linear system, can be efficiently overcome by a fully parallel implementation w.r.t. the frequencies, where each processor reconstructs a narrow-band image. All the other optimization steps are extremely fast. A parallel implementation of the algorithm in Julia is publicly available at https://github.com/andferrari. Preliminary simulations illustrate the performances of the method and its ability to reconstruct complex spatio-spectral structures.

preprint2016arXiv

Multitask diffusion adaptation over asynchronous networks

The multitask diffusion LMS is an efficient strategy to simultaneously infer, in a collaborative manner, multiple parameter vectors. Existing works on multitask problems assume that all agents respond to data synchronously. In several applications, agents may not be able to act synchronously because networks can be subject to several sources of uncertainties such as changing topology, random link failures, or agents turning on and off for energy conservation. In this work, we describe a model for the solution of multitask problems over asynchronous networks and carry out a detailed mean and mean-square error analysis. Results show that sufficiently small step-sizes can still ensure both stability and performance. Simulations and illustrative examples are provided to verify the theoretical findings. The framework is applied to a particular application involving spectral sensing.

preprint2015arXiv

Distributed image reconstruction for very large arrays in radio astronomy

Current and future radio interferometric arrays such as LOFAR and SKA are characterized by a paradox. Their large number of receptors (up to millions) allow theoretically unprecedented high imaging resolution. In the same time, the ultra massive amounts of samples makes the data transfer and computational loads (correlation and calibration) order of magnitudes too high to allow any currently existing image reconstruction algorithm to achieve, or even approach, the theoretical resolution. We investigate here decentralized and distributed image reconstruction strategies which select, transfer and process only a fraction of the total data. The loss in MSE incurred by the proposed approach is evaluated theoretically and numerically on simple test cases.

preprint2015arXiv

Large Scale 3D Image Reconstruction in Optical Interferometry

Astronomical optical interferometers (OI) sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid atmospheric perturbations, the phases of the complex Fourier samples (visibilities) cannot be directly exploited , and instead linear relationships between the phases are used (phase closures and differential phases). Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic OI instruments are now paving the way to multiwavelength imaging. This paper presents the derivation of a spatio-spectral ("3D") image reconstruction algorithm called PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm is able to solve large scale problems. It relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also from differential phases, which help to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm.

preprint2014arXiv

A graph Laplacian regularization for hyperspectral data unmixing

This paper introduces a graph Laplacian regularization in the hyperspectral unmixing formulation. The proposed regularization relies upon the construction of a graph representation of the hyperspectral image. Each node in the graph represents a pixel's spectrum, and edges connect spectrally and spatially similar pixels. The proposed graph framework promotes smoothness in the estimated abundance maps and collaborative estimation between homogeneous areas of the image. The resulting convex optimization problem is solved using the Alternating Direction Method of Multipliers (ADMM). A special attention is given to the computational complexity of the algorithm, and Graph-cut methods are proposed in order to reduce the computational burden. Finally, simulations conducted on synthetic data illustrate the effectiveness of the graph Laplacian regularization with respect to other classical regularizations for hyperspectral unmixing.

preprint2014arXiv

Blind and fully constrained unmixing of hyperspectral images

This paper addresses the problem of blind and fully constrained unmixing of hyperspectral images. Unmixing is performed without the use of any dictionary, and assumes that the number of constituent materials in the scene and their spectral signatures are unknown. The estimated abundances satisfy the desired sum-to-one and nonnegativity constraints. Two models with increasing complexity are developed to achieve this challenging task, depending on how noise interacts with hyperspectral data. The first one leads to a convex optimization problem, and is solved with the Alternating Direction Method of Multipliers. The second one accounts for signal-dependent noise, and is addressed with a Reweighted Least Squares algorithm. Experiments on synthetic and real data demonstrate the effectiveness of our approach.

preprint2014arXiv

PAINTER: a spatio-spectral image reconstruction algorithm for optical interferometry

Astronomical optical interferometers sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid perturbations caused by atmospheric turbulence, the phases of the complex Fourier samples (visibilities) cannot be directly exploited. Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic optical interferometric instruments are now paving the way to multiwavelength imaging. This paper is devoted to the derivation of a spatio-spectral (3D) image reconstruction algorithm, coined PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also differential phases, which helps to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm and in particular the relevance of injecting a differential phases model in the reconstruction.

preprint2013arXiv

A normalized scaled gradient method to solve non-negativity and equality constrained linear inverse problem - Application to spectral mixture analysis

This paper addresses the problem of minimizing a convex cost function under non-negativity and equality constraints, with the aim of solving the linear unmixing problem encountered in hyperspectral imagery. This problem can be formulated as a linear regression problem whose regression coefficients (abundances) satisfy sum-to-one and positivity constraints. A normalized scaled gradient iterative method (NSGM) is proposed for estimating the abundances of the linear mixing model. The positivity constraint is ensured by the Karush Kuhn Tucker conditions whereas the sum-to-one constraint is fulfilled by introducing normalized variables in the algorithm. The convergence is ensured by a one-dimensional search of the step size. Note that NSGM can be applied to any convex cost function with non negativity and flux constraints. In order to compare the NSGM with the well-known fully constraint least squares (FCLS) algorithm, this latter is reformulated in term of a penalized function, which reveals its suboptimality. Simulations on synthetic data illustrate the performances of the proposed algorithm in comparison with other unmixing algorithms and, more particulary, demonstrate its efficiency when compared to the popular FCLS. Finally, results on real data are given.