Automated detection of extended sources in radio maps: progress from the SCORPIO survey

Automated source extraction and parameterization represents a crucial challenge for the next-generation radio interferometer surveys, such as those performed with the Square Kilometre Array (SKA) and its precursors. In this paper we present a new algorithm, dubbed CAESAR (Compact And Extended Source Automated Recognition), to detect and parametrize extended sources in radio interferometric maps. It is based on a pre-filtering stage, allowing image denoising, compact source suppression and enhancement of diffuse emission, followed by an adaptive superpixel clustering stage for final source segmentation. A parameterization stage provides source flux information and a wide range of morphology estimators for post-processing analysis. We developed CAESAR in a modular software library, including also different methods for local background estimation and image filtering, along with alternative algorithms for both compact and diffuse source extraction. The method was applied to real radio continuum data collected at the Australian Telescope Compact Array (ATCA) within the SCORPIO project, a pathfinder of the ASKAP-EMU survey. The source reconstruction capabilities were studied over different test fields in the presence of compact sources, imaging artefacts and diffuse emission from the Galactic plane and compared with existing algorithms. When compared to a human-driven analysis, the designed algorithm was found capable of detecting known target sources and regions of diffuse emission, outperforming alternative approaches over the considered fields.

[1]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[2]  T. Budavari,et al.  Radio Continuum Surveys with Square Kilometre Array Pathfinders , 2012, Publications of the Astronomical Society of Australia.

[3]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[4]  Chandan Singh,et al.  Algorithms for fast computation of Zernike moments and their numerical stability , 2011, Image Vis. Comput..

[5]  Vyacheslav V. Kitaeff,et al.  The impact of JPEG2000 lossy compression on the scientific quality of radio astronomy imagery , 2014 .

[6]  Bülent Sankur,et al.  Survey over image thresholding techniques and quantitative performance evaluation , 2004, J. Electronic Imaging.

[7]  G. Dubner,et al.  A new study of the supernova remnant G344.7-0.1 located in the vicinity of the unidentified TeV source HESS J1702-420 , 2011 .

[8]  David Taubman,et al.  Astronomical imagery: Considerations for a contemporary approach with JPEG2000 , 2014, Astron. Comput..

[9]  E. Bertin,et al.  SExtractor: Software for source extraction , 1996 .

[10]  Matthew T. Whiting,et al.  duchamp: a 3D source finder for spectral‐line data , 2012, 1201.2710.

[11]  Gene H. Golub,et al.  Matrix computations , 1983 .

[12]  Sergey Brin,et al.  The Anatomy of a Large-Scale Hypertextual Web Search Engine , 1998, Comput. Networks.

[13]  A. L. Busfield,et al.  The RMS survey: radio observations of candidate massive YSOs in the southern hemisphere , 2006, Proceedings of the International Astronomical Union.

[14]  F. Rademakers,et al.  ROOT — An object oriented data analysis framework , 1997 .

[15]  M. J. Wolff,et al.  The Bubbling Galactic Disk , 2006 .

[16]  Andrew Hopkins,et al.  The Completeness and Reliability of Threshold and False-discovery Rate Source Extraction Algorithms for Compact Continuum Sources , 2011, Publications of the Astronomical Society of Australia.

[17]  R. Norris,et al.  SCORPIO: a deep survey of radio emission from the stellar life-cycle , 2015, 1506.04531.

[18]  David Mary,et al.  MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators - A sparse deconvolution algorithm for radio interferometric imaging , 2014, 1412.5387.

[19]  Andrew Hopkins,et al.  Compact continuum source finding for next generation radio surveys , 2012, 1202.4500.

[20]  Jian Sun,et al.  Guided Image Filtering , 2010, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[21]  James R. Curran,et al.  blobcat: software to catalogue flood‐filled blobs in radio images of total intensity and linear polarization , 2012, 1205.5313.

[22]  Ben Humphreys,et al.  Source-Finding for the Australian Square Kilometre Array Pathfinder , 2012, Publications of the Astronomical Society of Australia.

[23]  A. Green,et al.  The MOST supernova remnant catalogue (MSC). , 1996 .

[24]  Pascal Fua,et al.  SLIC Superpixels Compared to State-of-the-Art Superpixel Methods , 2012, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[25]  Christopher Hollitt,et al.  Feature Detection in Radio Astronomy using the Circle Hough Transform , 2012, Publications of the Astronomical Society of Australia.

[26]  Jianping Shi,et al.  Hierarchical Saliency Detection on Extended CSSD , 2014, ArXiv.

[27]  Christopher R. Genovese,et al.  A New Source Detection Algorithm Using the False-Discovery Rate , 2001, astro-ph/0110570.

[28]  D. J. Saikia,et al.  EMU: Evolutionary Map of the Universe , 2011, Publications of the Astronomical Society of Australia.

[29]  J. R. Sanchez-Sutil,et al.  Exploring three faint source detections methods for aperture synthesis radio images , 2015 .

[30]  A. Shulevski,et al.  The ASKAP/EMU Source Finding Data Challenge , 2015, Publications of the Astronomical Society of Australia.

[31]  P. Rousseeuw,et al.  Unmasking Multivariate Outliers and Leverage Points , 1990 .

[32]  Tony F. Chan,et al.  A Multiphase Level Set Framework for Image Segmentation Using the Mumford and Shah Model , 2002, International Journal of Computer Vision.

[33]  Kidiyo Kpalma,et al.  Affine invariance contour descriptor based on iso-area normalisation , 2007 .