pymf.filter_map_mmf

pymf.filter_map_mmf(images, source, spec, noise_maps=None, sigma_noise=None, no_cross=False, global_cross=False)

Computes and applies a matched multifilter to a set of multi-frequency map. The filter is build from provided spatial and spectral source templates and the cross spectrum of a set of maps at different frequencies.

Parameters:

images: 3D float array

Data cube containing the input images that are provided at n_f frequencies. The dimensions are supposed to be n_f * n_x * n_y, where n_x and n_y are the number of pixels along the * and y axis of the images.

source: 2D or 3D float array

Source template. Can either be a single image or alternatively a data cube if each image comes at a different spatial resolution. Needs to have identical dimensions as images if a cube is provided. Otherwise the single source image needs to match the dimensions of a single image.

spec: float array

SED of the source. Needs to be provided at the same n_f frequencies as the input images.

noise_maps: 3D float array, optional

Data cube containing provided noise maps that will be used to compute noise power spectrum matrix. Otherwise the noise power spectrum matrix is directly computed from the images. noise_maps need to have the same dimensions as images. Default: None

sigma_noise: float array, optional

Array containing the standard deviation of the noise for the n_f maps. Does allow the analytic computation of the noise power spectrum. Useful for tests with white noise. Default: None

no_cross: bool, optional

If True the frequency-to-frequency cross terms of the n_f maps will be ignored, i.e. non-diagonal elements of the noise power spectrum matrix are set to zero. Default: False

global_cross: bool, optional

If set to True, the cross power of the maps is computed globally (i.e. averaged over the full maps / all scales). If False, the cross power is computed at each spatial frequency. Default: False

Returns:

filtered_image: 2D float array

Linear combination of the images convolved with the computed matched multifilter.

filter_ft: 3D float array

Data cube containing the 2D window functions of the matched multifilter.

noise: float

Noise level of the filtered map.