Multi-frequency image reconstruction for radio-interferometry with self-tuned regularization parameters

As the world's largest radio telescope, the Square Kilometer Array (SKA) will provide radio interferometric data with unprecedented detail. Image reconstruction algorithms for radio interferometry are challenged to scale well with TeraByte image sizes never seen before. In this work, we investigate one such 3D image reconstruction algorithm known as MUFFIN (MUlti-Frequency image reconstruction For radio INterferometry). In particular, we focus on the challenging task of automatically finding the optimal regularization parameter values. In practice, finding the regularization parameters using classical grid search is computationally intensive and nontrivial due to the lack of ground-truth. We adopt a greedy strategy where, at each iteration, the optimal parameters are found by minimizing the predicted Stein unbiased risk estimate (PSURE). The proposed self-tuned version of MUFFIN involves parallel and computationally efficient steps, and scales well with large-scale data. Finally, numerical results on a 3D image are presented to showcase the performance of the proposed approach.

[1]  Laurent Condat,et al.  A Generic Proximal Algorithm for Convex Optimization—Application to Total Variation Minimization , 2014, IEEE Signal Processing Letters.

[2]  C. Stein Estimation of the Mean of a Multivariate Normal Distribution , 1981 .

[3]  David Mary,et al.  Distributed multi-frequency image reconstruction for radio-interferometry , 2016, 2016 24th European Signal Processing Conference (EUSIPCO).

[4]  Bang Công Vu,et al.  A splitting algorithm for dual monotone inclusions involving cocoercive operators , 2011, Advances in Computational Mathematics.

[5]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[6]  Jeffrey A. Fessler,et al.  Regularization Parameter Selection for Nonlinear Iterative Image Restoration and MRI Reconstruction Using GCV and SURE-Based Methods , 2012, IEEE Transactions on Image Processing.

[7]  Teresa Reginska,et al.  A Regularization Parameter in Discrete Ill-Posed Problems , 1996, SIAM J. Sci. Comput..

[8]  W. Clem Karl,et al.  3.6 – Regularization in Image Restoration and Reconstruction , 2005 .

[9]  V. Morozov On the solution of functional equations by the method of regularization , 1966 .

[10]  Henk Corporaal,et al.  An End-to-End Computing Model for the Square Kilometre Array , 2014, Computer.

[11]  Alexandru Onose,et al.  A low-rank and joint-sparsity model for hyper-spectral radio-interferometric imaging , 2016, 2016 24th European Signal Processing Conference (EUSIPCO).

[12]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[13]  Rainer Beck,et al.  Square kilometre array , 2010, Scholarpedia.

[14]  Alle-Jan van der Veen,et al.  Signal Processing Tools for Radio Astronomy , 2013, Handbook of Signal Processing Systems.

[15]  Jean-Philippe Thiran,et al.  Scalable splitting algorithms for big-data interferometric imaging in the SKA era , 2016, 1601.04026.

[16]  Yves Wiaux,et al.  PURIFY: a new approach to radio-interferometric imaging , 2013, 1307.4370.

[17]  Yonina C. Eldar Generalized SURE for Exponential Families: Applications to Regularization , 2008, IEEE Transactions on Signal Processing.

[18]  Robert C. Kennicutt,et al.  MULTI-SCALE CLEAN: A COMPARISON OF ITS PERFORMANCE AGAINST CLASSICAL CLEAN ON GALAXIES USING THINGS , 2008, 0810.2825.

[19]  T. Murphy,et al.  wsclean: an implementation of a fast, generic wide-field imager for radio astronomy , 2014, 1407.1943.

[20]  M. Hutchinson A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines , 1989 .

[21]  David Mary,et al.  Distributed image reconstruction for very large arrays in radio astronomy , 2014, 2014 IEEE 8th Sensor Array and Multichannel Signal Processing Workshop (SAM).

[22]  H. Junklewitz,et al.  A new approach to multifrequency synthesis in radio interferometry , 2014, 1401.4711.

[23]  Yonina C. Eldar,et al.  The Projected GSURE for Automatic Parameter Tuning in Iterative Shrinkage Methods , 2010, Applied and Computational Harmonic Analysis.

[24]  Tim J. Cornwell,et al.  The Noncoplanar Baselines Effect in Radio Interferometry: The W-Projection Algorithm , 2008, IEEE Journal of Selected Topics in Signal Processing.

[25]  Alexandru Onose,et al.  A randomised primal-dual algorithm for distributed radio-interferometric imaging , 2016, 2016 24th European Signal Processing Conference (EUSIPCO).

[26]  Pascal Bianchi,et al.  Two distributed algorithms for the deconvolution of large radio-interferometric multispectral images , 2016, 2016 24th European Signal Processing Conference (EUSIPCO).

[27]  Stephanie Thalberg,et al.  Interferometry And Synthesis In Radio Astronomy , 2016 .

[28]  J. Högbom,et al.  APERTURE SYNTHESIS WITH A NON-REGULAR DISTRIBUTION OF INTERFEROMETER BASELINES. Commentary , 1974 .

[29]  J. Anderson,et al.  LOFAR sparse image reconstruction , 2014, 1406.7242.

[30]  O. Smirnov,et al.  The MeqTrees software system and its use for third-generation calibration of radio interferometers , 2010, 1101.1745.