This work presents a method for determining the accuracy of a source finder algorithm for spectral line radio astronomy data and the Source Finder Accuracy Evaluator (SFAE), a program that implements this method. The accuracy of a source finder is defined in terms of its completeness, reliability, and accuracy of the parameterisation of the sources that were found. These values are calculated by executing the source finder on an image with a known source catalogue, then comparing the output of the source finder to the known catalogue. The intended uses of SFAE include determining the most accurate source finders for use in a survey, determining the types of radio sources a particular source finder is capable of accurately locating, and identifying optimum parameters and areas of improvement for these algorithms. This paper demonstrates a sample of accuracy information that can be obtained through this method, using a simulated ASKAP data cube and the duchamp source finder.

This project uses the 2MASS all-sky image database to study the structure of galaxies over a range of luminosities, sizes and morphological types. This first paper in this series will outline the techniques, reliability and data products to our surface photometry program. Our program will analyze all acceptable galaxies (meeting our criteria for isolation from companions and bright stars) from the Revised Shapley–Ames and Uppsala galaxy catalogs. Resulting photometry and surface brightness profiles are released using a transparent scheme of data storage which includes not only all the processed data but knowledge of the processing steps and calibrating parameters.

While automatic detection of point sources in astronomical images has experienced a great degree of success, less effort has been directed towards the detection of extended and low-surface-brightness features. At present, existing telescopes still rely on human expertise to reduce the raw data to usable images and then to analyse the images for non-pointlike objects. However, the next generation of radio telescopes will generate unprecedented volumes of data making manual data reduction and object extraction infeasible. Without developing new methods of automatic detection for extended and diffuse objects such as supernova remnants, bent-tailed galaxies, radio relics and halos, a wealth of scientifically important results will not be uncovered. In this paper we explore the response of the Circle Hough Transform to a representative sample of different extended circular or arc-like astronomical objects. We also examine the response of the Circle Hough Transform to input images containing noise alone and inputs including point sources.

Chromospheric brightening and Hα surges are evident and common phenomena along sunspot light bridges. In this paper, a coronal jet ejection from a sunspot light bridge is presented. Using data from the Solar Dynamics Observatory and Hinode satellites, it is confirmed that the jet has its root near the light bridge. This suggests that the jet may be a result of reconnection between the main sunspot and the light bridge. Due to the processing of jet ejecta, the intensity and width of the light bridge show changes to some extent. This also suggests that the jet is related to the interaction between light bridge and umbra, possibly owing to magnetic reconnection or heated plasma trapped in the light bridge escaping and moving along field lines.

This is a design study into the capabilities of the Australian Square Kilometre Array Pathfinder in performing a full-sky low redshift neutral hydrogen survey, termed WALLABY, and the potential cosmological constraints one can attain from measurement of the galaxy power spectrum. We find that the full sky survey will likely attain 6 × 105 redshifts which, when combined with expected Planck CMB data, will constrain the Dark Energy equation of state to 20%, representing a coming of age for radio observations in creating cosmological constraints.

The large size of the ASKAP Hi surveys DINGO and WALLABY necessitates automated 3D source finding. A performance difference of a few percent corresponds to a significant number of galaxies being detected or undetected. As such, the performance of the automated source finding is of paramount importance to both of these surveys. We have analysed the performance of various source finders to determine which will allow us to meet our survey goals during the DINGO and WALLABY design studies. Here we present a comparison of the performance of five different methods of automated source finding. These source finders are duchamp, gamma-finder, a cnhi finder, a 2d–1d wavelet reconstruction finder and a sigma clipping method (s + c finder). Each source finder was applied to the same three-dimensional data cubes containing (a) point sources with a Gaussian velocity profile and (b) spatially extended model-galaxies with inclinations and rotation profiles. We focus on the completeness and reliability of each algorithm when comparing the performance of the different source finders.

We present a framework to volume-render three-dimensional data cubes interactively using distributed ray-casting and volume-bricking over a cluster of workstations powered by one or more graphics processing units (GPUs) and a multi-core central processing unit (CPU). The main design target for this framework is to provide an in-core visualization solution able to provide three-dimensional interactive views of terabyte-sized data cubes. We tested the presented framework using a computing cluster comprising 64 nodes with a total of 128 GPUs. The framework proved to be scalable to render a 204 GB data cube with an average of 30 frames per second. Our performance analyses also compare the use of NVIDIA Tesla 1060 and 2050 GPU architectures and the effect of increasing the visualization output resolution on the rendering performance. Although our initial focus, as shown in the examples presented in this work, is volume rendering of spectral data cubes from radio astronomy, we contend that our approach has applicability to other disciplines where close to real-time volume rendering of terabyte-order three-dimensional data sets is a requirement.

The Square Kilometre Array and its pathfinders ASKAP and MeerKAT will produce prodigious amounts of data that necessitate automated source finding. The performance of automated source finders can be improved by pre-processing a dataset. In preparation for the WALLABY and DINGO surveys, we have used a test Hi datacube constructed from actual Westerbork Telescope noise and WHISP Hi galaxies to test the real world improvement of linear smoothing, the Duchamp source finder's wavelet denoising, iterative median smoothing and mathematical morphology subtraction, on intensity threshold source finding of spectral line datasets. To compare these pre-processing methods we have generated completeness-reliability performance curves for each method and a range of input parameters. We find that iterative median smoothing produces the best source finding results for ASKAP Hi spectral line observations, but wavelet denoising is a safer pre-processing technique. In this paper we also present our implementations of iterative median smoothing and mathematical morphology subtraction.

Here I will outline successes and challenges for finding spectral line sources in large data cubes that are dominated by noise. This is a 3D challenge as the sources we wish to catalog are spread over several spatial pixels and spectral channels. While 2D searches can be applied, e.g. channel by channel, optimal searches take into account the 3-dimensional nature of the sources. In this overview I will focus on Hi 21-cm spectral line source detection in extragalactic surveys, in particular HIPASS, the Hi Parkes All-Sky Survey and WALLABY, the ASKAP Hi All-Sky Survey. I use the original HIPASS data to highlight the diversity of spectral signatures of galaxies and gaseous clouds, both in emission and absorption. Among others, I report the discovery of a 680 km s−1 wide Hi absorption trough in the megamaser galaxy NGC 5793. Issues such as source confusion and baseline ripples, typically encountered in single-dish Hi surveys, are much reduced in interferometric Hi surveys. Several large Hi emission and absorption surveys are planned for the Australian Square Kilometre Array Pathfinder (ASKAP): here we focus on WALLABY, the 21-cm survey of the sky (δ < + 30°; z < 0.26) which will take about one year of observing time with ASKAP. Novel phased array feeds (‘radio cameras’) will provide 30 deg2 instantaneous field-of-view. WALLABY is expected to detect more than 500 000 galaxies, unveil their large-scale structures and cosmological parameters, detect their extended, low-surface-brightness disks as well as gas streams and filaments between galaxies. It is a precursor for future Hi surveys with SKA Phase I and II, exploring galaxy formation and evolution. The compilation of highly reliable and complete source catalogs will require sophisticated source-finding algorithms as well as accurate source parametrisation.

The Australian Square Kilometre Array Pathfinder (ASKAP) presents a number of challenges in the area of source finding and cataloguing. The data rates and image sizes are very large, and require automated processing in a high-performance computing environment. This requires development of new tools, that are able to operate in such an environment and can reliably handle large datasets. These tools must also be able to accommodate the different types of observations ASKAP will make: continuum imaging, spectral-line imaging, transient imaging. The ASKAP project has developed a source-finder known as selavy, built upon the duchamp source-finder. selavy incorporates a number of new features, which we describe here.

Since distributed processing of large images and cubes will be essential, we describe the algorithms used to distribute the data, find an appropriate threshold and search to that threshold and form the final source catalogue. We describe the algorithm used to define a varying threshold that responds to the local, rather than global, noise conditions, and provide examples of its use. And we discuss the approach used to apply two-dimensional fits to detected sources, enabling more accurate parameterisation. These new features are compared for timing performance, where we show that their impact on the pipeline processing will be small, providing room for enhanced algorithms.

We also discuss the development process for ASKAP source finding software. By the time of ASKAP operations, the ASKAP science community, through the Survey Science Projects, will have contributed important elements of the source finding pipeline, and the mechanisms in which this will be done are presented.