-
The Mexican Hat Wavelet Family. Application to point source detection in CMB maps
Authors:
J. Gonzalez-Nuevo,
F. Argueso,
M. Lopez-Caniego,
L. Toffolatti,
J. L. Sanz,
P. Vielva,
D. Herranz
Abstract:
We propose a new detection technique in the plane based on an isotropic wavelet family. This family is naturally constructed as an extension of the Gaussian-Mexican Hat Wavelet pair and for that reason we call it the Mexican Hat Wavelet Family (MHWF). We show the performance of these wavelets when dealing with the detection of extragalactic point sources in cosmic microwave background (CMB) maps…
▽ More
We propose a new detection technique in the plane based on an isotropic wavelet family. This family is naturally constructed as an extension of the Gaussian-Mexican Hat Wavelet pair and for that reason we call it the Mexican Hat Wavelet Family (MHWF). We show the performance of these wavelets when dealing with the detection of extragalactic point sources in cosmic microwave background (CMB) maps: a very important issue within the most general problem of the component separation of the microwave sky. Specifically, flat two-dimensional simulations of the microwave sky comprising all astrophysical components plus instrumental noise have been analyzed for the channels at 30, 44 and 70 GHz of the forthcoming ESA's Planck mission Low Frequency Instrument (LFI). We adopt up-to-date cosmological evolution models of extragalactic sources able to fit well the new data on high-frequency radio surveys and we discuss our current results on point source detection by comparing them with those obtained using the Mexican Hat Wavelet (MHW) technique, which has been already proven a suitable tool for detecting point sources. By assuming a 5% reliability level, the first new members of the MHWF, at their ``optimal scale'', provide three point source catalogues on half of the sky (at galactic latitude |b|> 30) at 30, 44 and 70 GHz of 639, 387 and 340 extragalactic sources, respectively. The corresponding flux detection limits are 0.38, 0.45 and 0.47 Jy . By using the same simulated sky patches and at the same frequencies as before, the MHW at its optimal scale provides 543, 322 and 311 sources with flux detection limits of 0.44, 0.51 and 0.50 Jy, respectively (5% reliability level). These results show a clear improvement when we use the new members of the MHWF and, in particular, the MHW2 with respect to the MHW.
△ Less
Submitted 21 April, 2006; v1 submitted 18 April, 2006;
originally announced April 2006.
-
Some comments on the note "Some comments on the paper 'Filter design for the detection of compact sources based on the Neyman-Pearson detector' by M. Lopez-Caniego et. al (2005, MNRAS 359, 993)" by R. Vio and P. Andreani (astroph/0509394)
Authors:
M. Lopez-Caniego,
D. Herranz,
R. B. Barreiro,
J. L. Sanz
Abstract:
In this note we stress the necessity of a careful check of the arguments used by Vio & Andreani (2005) (VA hereinafter) to criticise the superior performance of the biparametric scale adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random Gaussian background. In particular, we point out that a defective reading and understanding of previous…
▽ More
In this note we stress the necessity of a careful check of the arguments used by Vio & Andreani (2005) (VA hereinafter) to criticise the superior performance of the biparametric scale adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random Gaussian background. In particular, we point out that a defective reading and understanding of previous works in the literature (Rice 1954; Barreiro et al. 2003; Lopez-Caniego et al. 2005) leads the authors of VA to the derivation of an incorrect formula and to some misleading conclusions.
△ Less
Submitted 23 September, 2005;
originally announced September 2005.
-
Detection of point sources on two-dimensional images based on peaks
Authors:
M. Lopez-Caniego,
D. Herranz,
J. L. Sanz,
R. B. Barreiro
Abstract:
This article considers the detection of point sources in two dimensional astronomical images. The detection scheme we propose is based on peak statistics. We discuss the example of the detection of far galaxies in Cosmic Microwave Background experiments throughout the paper, although the method we present is totally general and can be used in many other fields of data analysis. We assume sources…
▽ More
This article considers the detection of point sources in two dimensional astronomical images. The detection scheme we propose is based on peak statistics. We discuss the example of the detection of far galaxies in Cosmic Microwave Background experiments throughout the paper, although the method we present is totally general and can be used in many other fields of data analysis. We assume sources with a Gaussian profile --that is a fair approximation of the profile of a point source convolved with the detector beam in microwave experiments-- on a background modeled by a homogeneous and isotropic Gaussian random field characterized by a scale-free power spectrum. Point sources are enhanced with respect to the background by means of linear filters. After filtering, we identify local maxima and apply our detection scheme, a Neyman-Pearson detector that defines our region of acceptance based on the a priori pdf of the sources and the ratio of number densities. We study the different performances of some linear filters that have been used in this context in the literature: the Mexican Hat wavelet, the matched filter and the scale-adaptive filter. We consider as well an extension to two dimensions of the biparametric scale adaptive filter (BSAF). The BSAF depends on two parameters which are determined by maximizing the number density of real detections while fixing the number density of spurious detections. For our detection criterion the BSAF outperforms the other filters in the interesting case of white noise.
△ Less
Submitted 7 March, 2005;
originally announced March 2005.
-
Filter design for the detection of compact sources based on the Neyman-Pearson detector
Authors:
M. Lopez-Caniego,
D. Herranz,
R. B. Barreiro,
J. L. Sanz
Abstract:
This paper considers the problem of compact source detection on a Gaussian background in 1D. Two aspects of this problem are considered: the design of the detector and the filtering of the data. Our detection scheme is based on local maxima and it takes into account not only the amplitude but also the curvature of the maxima. A Neyman-Pearson test is used to define the region of acceptance, that…
▽ More
This paper considers the problem of compact source detection on a Gaussian background in 1D. Two aspects of this problem are considered: the design of the detector and the filtering of the data. Our detection scheme is based on local maxima and it takes into account not only the amplitude but also the curvature of the maxima. A Neyman-Pearson test is used to define the region of acceptance, that is given by a sufficient linear detector that is independent on the amplitude distribution of the sources. We study how detection can be enhanced by means of linear filters with a scaling parameter and compare some of them (the Mexican Hat wavelet, the matched and the scale-adaptive filters). We introduce a new filter, that depends on two free parameters (biparametric scale-adaptive filter). The value of these two parameters can be determined, given the a priori pdf of the amplitudes of the sources, such that the filter optimizes the performance of the detector in the sense that it gives the maximum number of real detections once fixed the number density of spurious sources. The combination of a detection scheme that includes information on the curvature and a flexible filter that incorporates two free parameters (one of them a scaling) improves significantly the number of detections in some interesting cases. In particular, for the case of weak sources embedded in white noise the improvement with respect to the standard matched filter is of the order of 40%. Finally, an estimation of the amplitude of the source is introduced and it is proven that such an estimator is unbiased and it has maximum efficiency. We perform numerical simulations to test these theoretical ideas and conclude that the results of the simulations agree with the analytical ones.
△ Less
Submitted 11 April, 2005; v1 submitted 7 March, 2005;
originally announced March 2005.
-
A Bayesian approach to filter design: detection of compact sources
Authors:
M. Lopez-Caniego,
D. Herranz,
R. B. Barreiro,
J. L. Sanz
Abstract:
We consider filters for the detection and extraction of compact sources on a background. We make a one-dimensional treatment (though a generalization to two or more dimensions is possible) assuming that the sources have a Gaussian profile whereas the background is modeled by an homogeneous and isotropic Gaussian random field, characterized by a scale-free power spectrum. Local peak detection is…
▽ More
We consider filters for the detection and extraction of compact sources on a background. We make a one-dimensional treatment (though a generalization to two or more dimensions is possible) assuming that the sources have a Gaussian profile whereas the background is modeled by an homogeneous and isotropic Gaussian random field, characterized by a scale-free power spectrum. Local peak detection is used after filtering. Then, a Bayesian Generalized Neyman-Pearson test is used to define the region of acceptance that includes not only the amplification but also the curvature of the sources and the a priori probability distribution function of the sources. We search for an optimal filter between a family of Matched-type filters (MTF) modifying the filtering scale such that it gives the maximum number of real detections once fixed the number density of spurious sources. We have performed numerical simulations to test theoretical ideas.
△ Less
Submitted 7 March, 2005;
originally announced March 2005.
-
The estimation of the SZ effects with unbiased multifilters
Authors:
D. Herranz,
J. L. Sanz,
R. B. Barreiro,
M. Lopez-Caniego
Abstract:
In this work we study the performance of linear multifilters for the estimation of the amplitudes of the thermal and kinematic Sunyaev-Zel'dovich effects. We show that when both effects are present, estimation of these effects with standard matched multifilters is intrinsically biased. This bias is due to the fact that both signals have basically the same spatial profile. We find a new family of…
▽ More
In this work we study the performance of linear multifilters for the estimation of the amplitudes of the thermal and kinematic Sunyaev-Zel'dovich effects. We show that when both effects are present, estimation of these effects with standard matched multifilters is intrinsically biased. This bias is due to the fact that both signals have basically the same spatial profile. We find a new family of multifilters related to the matched multifilters that cancel this systematic bias, hence we call them Unbiased Matched Multifilters. We test the unbiased matched multifilters and compare them with the standard matched multifilters using simulations that reproduce the future Planck mission's observations. We find that in the case of the standard matched multifilters the systematic bias in the estimation of the kinematic Sunyaev-Zel'dovich effect can be very large, even greater than the statistical error bars. Unbiased matched multifilters cancel effectively this kind of bias. In concordance with other works in the literature, our results indicate that the sensitivity and resolution of Planck will not be enough to give reliable estimations of the kinematic Sunyaev-Zel'dovich of individual clusters. However, since the estimation with the unbiased matched multifilters is not intrinsically biased, it can be possible to use them to statistically study peculiar velocities in large scales using large sets of clusters.
△ Less
Submitted 9 June, 2004;
originally announced June 2004.