Simulating full-sky interferometric observations

Aperture array interferometers, such as that proposed for the Square Kilometre Array (SKA), will see the entire sky, hence the standard approach to simulating visibilities will not be applicable since it relies on a tangent plane approximation that is valid only for small fields of view. We derive interferometric formulations in real, spherical harmonic and wavelet space that include contributions over the entire sky and do not rely on any tangent plane approximations. A fast wavelet method is developed to simulate the visibilities observed by an interferometer in the full-sky setting. Computing visibilities using the fast wavelet method adapts to the sparse representation of the primary beam and sky intensity in the wavelet basis. Consequently, the fast wavelet method exhibits superior computational complexity to the real and spherical harmonic space methods and may be performed at substantially lower computational cost, while introducing only negligible error to simulated visibilities. Low-resolution interferometric observations are simulated using all of the methods to compare their performance, demonstrating that the fast wavelet method is approximately three times faster that the other methods for these low-resolution simulations. The computational burden of the real and spherical harmonic space methods renders these techniques computationally infeasible for higher resolution simulations. High-resolution interferometric observations are simulated using the fast wavelet method only, demonstrating and validating the application of this method to realistic simulations. The fast wavelet method is estimated to provide a greater than ten-fold reduction in execution time compared to the other methods for these high-resolution simulations.Comment: 16 pages, 9 figures, replaced to match version accepted by MNRAS (major additions to previous version including new fast wavelet method

On The Continuous Steering of the Scale of Tight Wavelet Frames

In analogy with steerable wavelets, we present a general construction of adaptable tight wavelet frames, with an emphasis on scaling operations. In particular, the derived wavelets can be "dilated" by a procedure comparable to the operation of steering steerable wavelets. The fundamental aspects of the construction are the same: an admissible collection of Fourier multipliers is used to extend a tight wavelet frame, and the "scale" of the wavelets is adapted by scaling the multipliers. As an application, the proposed wavelets can be used to improve the frequency localization. Importantly, the localized frequency bands specified by this construction can be scaled efficiently using matrix multiplication

Sparse image reconstruction on the sphere: analysis and synthesis

We develop techniques to solve ill-posed inverse problems on the sphere by sparse regularisation, exploiting sparsity in both axisymmetric and directional scale-discretised wavelet space. Denoising, inpainting, and deconvolution problems, and combinations thereof, are considered as examples. Inverse problems are solved in both the analysis and synthesis settings, with a number of different sampling schemes. The most effective approach is that with the most restricted solution-space, which depends on the interplay between the adopted sampling scheme, the selection of the analysis/synthesis problem, and any weighting of the l1 norm appearing in the regularisation problem. More efficient sampling schemes on the sphere improve reconstruction fidelity by restricting the solution-space and also by improving sparsity in wavelet space. We apply the technique to denoise Planck 353 GHz observations, improving the ability to extract the structure of Galactic dust emission, which is important for studying Galactic magnetism.Comment: 11 pages, 6 Figure

Algorithmic options for joint time-frequency analysis in structural dynamics applications

The purpose of this paper is to present recent research efforts by the authors supporting the superiority of joint time-frequency analysis over the traditional Fourier transform in the study of non-stationary signals commonly encountered in the fields of earthquake engineering, and structural dynamics. In this respect, three distinct signal processing techniques appropriate for the representation of signals in the time-frequency plane are considered. Namely, the harmonic wavelet transform, the adaptive chirplet decomposition, and the empirical mode decomposition, are utilized to analyze certain seismic accelerograms, and structural response records. Numerical examples associated with the inelastic dynamic response of a seismically-excited 3-story benchmark steel-frame building are included to show how the mean-instantaneous-frequency, as derived by the aforementioned techniques, can be used as an indicator of global structural damage

Automated analysis of quantitative image data using isomorphic functional mixed models, with application to proteomics data

Image data are increasingly encountered and are of growing importance in many areas of science. Much of these data are quantitative image data, which are characterized by intensities that represent some measurement of interest in the scanned images. The data typically consist of multiple images on the same domain and the goal of the research is to combine the quantitative information across images to make inference about populations or interventions. In this paper we present a unified analysis framework for the analysis of quantitative image data using a Bayesian functional mixed model approach. This framework is flexible enough to handle complex, irregular images with many local features, and can model the simultaneous effects of multiple factors on the image intensities and account for the correlation between images induced by the design. We introduce a general isomorphic modeling approach to fitting the functional mixed model, of which the wavelet-based functional mixed model is one special case. With suitable modeling choices, this approach leads to efficient calculations and can result in flexible modeling and adaptive smoothing of the salient features in the data. The proposed method has the following advantages: it can be run automatically, it produces inferential plots indicating which regions of the image are associated with each factor, it simultaneously considers the practical and statistical significance of findings, and it controls the false discovery rate.Comment: Published in at http://dx.doi.org/10.1214/10-AOAS407 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org