A Mexican Hat with holes: calculating low resolution power spectra from data with gaps

Abstract

A simple method for calculating a low-resolution power spectrum from data with gaps is described. The method is a modification of the $\Delta$-variance method previously described by Stutzki and Ossenkopf. A Mexican Hat filter is used to single out fluctuations at a given spatial scale and the variance of the convolved image is calculated. The gaps in the image, defined by the mask, are corrected for by representing the Mexican Hat filter as a difference between two Gaussian filters with slightly different widths, convolving the image and mask with these filters and dividing the results before calculating the final filtered image. This method cleanly compensates for data gaps even if these have complicated shapes and cover a significant fraction of the data. The method was developed to deal with problematic 2D images, where irregular detector edges and masking of contaminating sources compromise the power spectrum estimates, but it can also be straightforwardly applied to 1D timing analysis or 3D data cubes from numerical simulations.

accepted for publication in MNRAS

References in corpus (2)

• Structure analysis of interstellar clouds: I. Improving the Delta-variance method
• Long-term optical variability of PKS 2155-304

Cited by in corpus (6)

• Turbulent Heating in Galaxy Clusters Brightest in X-rays
• The 2-79 keV X-ray Spectrum of the Circinus Galaxy with NuSTAR, XMM-Newton and Chandra: a Fully Compton-Thick AGN
• The relation between gas density and velocity power spectra in galaxy clusters: high-resolution hydrodynamic simulations and the role of conduction
• The variable hard X-ray emission of NGC4945 as observed by NuSTAR
• Suppressed effective viscosity in the bulk intergalactic plasma
• Deep Chandra observations of the stripped galaxy group falling into Abell 2142