# Detection/estimation of the modulus of a vector. Application to point source detection in polarization data

## Abstract

Given a set of images, whose pixel values can be considered as the components of a vector, it is interesting to estimate the modulus of such a vector in some localised areas corresponding to a compact signal. For instance, the detection/estimation of a polarized signal in compact sources immersed in a background is relevant in some fields like astrophysics. We develop two different techniques, one based on the Neyman-Pearson lemma, the Neyman-Pearson filter (NPF), and another based on prefiltering-before-fusion, the filtered fusion (FF), to deal with the problem of detection of the source and estimation of the polarization given two or three images corresponding to the different components of polarization (two for linear polarization, three including circular polarization). For the case of linear polarization, we have performed numerical simulations on two-dimensional patches to test these filters following two different approaches (a blind and a non-blind detection), considering extragalactic point sources immersed in cosmic microwave background (CMB) and non-stationary noise with the conditions of the 70 GHz *Planck* channel. The FF outperforms the NPF, especially for low fluxes. We can detect with the FF extragalactic sources in a high noise zone with fluxes Jy for (blind/non-blind) detection and in a low noise zone with fluxes Jy for (blind/non-blind) detection with low errors in the estimated flux and position.

## 1Introduction

The detection and estimation of the intensity of compact objects (or small regions) embedded in a background plus instrumental noise is relevant in different contexts, e.g. astrophysics, cosmology, medicine, teledetection, radar, etc. Different techniques have proven useful in the literature. Some of the proposed techniques are frequentist, such as the standard matched filter [26], the matched multifilter [12] or the recently developed matched matrix filters [14] that correspond to scalar, vector or matrix filters, respectively. Other frequentist techniques include continuous wavelets like the standard Mexican Hat [39] and other members of its family [9] and, more generally, filters based on the Neyman-Pearson approach using the distribution of maxima [21]. All of these filters have been used in the literature, in particular, for the detection and estimation of the intensity of point-like sources (i.e. extragalactic objects) in Cosmic Microwave Background (CMB) maps. In addition, some of them have been applied to real data like those obtained by the WMAP satellite [24] and simulated data [23] for the upcoming experiment on board the *Planck* satellite [33]. Besides, Bayesian methods have also been recently developed [16].

In some applications it is important to measure not only the intensity of the light (signal) but also its polarization. An example is the study of CMB radiation. The polarization is given by the Stokes parameters , and the total intensity of polarization is [18]. For linear polarization, the previous expression reduces to . In such cases, three or two images are added quadratically followed by a square root.

Let us consider the case of linear polarization. In this case we have two images and different approaches can be used to deal with detection/estimation of point-like sources on these maps. On the one hand, one can try to get the polarization directly on the -map. In this approach, we will consider one filter, obtained through the Neyman-Pearson technique (NPF). On the other hand, we can operate with two matched filters, each one on and followed by a quadratic fusion and square root. We will call this procedure filtered fusion (FF). It is clear that from a formal point of view we are trying to ask which is the optimal way to find the modulus of a vector given the components. In the case we have only the map of the modulus of a vector and the components are unknown, the FF cannot be applied and the only possibility is the NPF.

We will develop the methodology for the cases of a 2-vector and a 3-vector because of the possible interesting applications to the 2-plane and 3-space . We will show the results when using numerical simulations on flat patches that are relevant for the component separation of sources (linear-polarized extragalactic sources) in CMB maps.

In Section 2, we develop the methodology. In Section 3, we describe the numerical simulations done in order to test the previous techniques. In Section 4 we present the main results and in Section 5 we give the main conclusions.

## 2Methodology

### 2.12-vector

To develop our methodology, we shall assume that we have a compact source, located for simplicity at the centre of two images , and characterised by amplitudes , and a profile , immersed in noises and that are Gaussian and independently distributed with zero mean and dispersions and . Again, for simplicity, we will consider that , a condition that is verified in most polarization detectors. In general, we will consider that the noise is non-stationary. We remark that the previous assumptions can be easily generalised to different profiles and different types of noise in the two images but we will not consider it in this paper. We will assume a linear model for the two images

The -map, , is characterised by a source at the centre of the image with amplitude immersed in non-additive noise which is correlated with the signal.

#### Neyman-Pearson filter (NPF) on the -map

If the noise is distributed normally and independently, then at any point the integration over the polar angle leads to the 2D-Rayleigh distribution of in absence of the source [27]

whereas if the source is present, with amplitude , one obtains the Rice distribution [29]

where is the modified Bessel function of zero order. If our image is pixelised, the different data , , with the number of pixels, will follow the two distributions

being the value of in the pixel, , the null hypothesis (absence of source) and the alternative (presence of source), respectively, and the profile at the pixel. The log-likelihood is defined by

The maximum likelihood estimator of the amplitude, , is given by the solution of the equation

This equation can be interpreted as a non-linear filter operating on the data that we shall call the Neyman-Pearson filter (NPF).

#### Filtered fusion (FF)

In this case we use the same matched filter (MF) operating over the two images , , respectively, as given by [3]:

Then, with the two filtered images we make the non-linear fusion pixel by pixel.

### 2.23-vector

Now, we shall assume that we have a compact source at the centre of three images , , characterised by amplitudes , , and a profile immersed in noise that is Gaussian and independently distributed with zero mean and dispersion . In general, we will consider that the noise is non-stationary. We will assume a linear model for the three images

The -map, , is characterised by a source at the centre of the image with amplitude immersed in non-additive noise correlated with the signal.

#### Neyman-Pearson filter (NPF) on the -map

If the noise is distributed normally and independently, then at any point the integration over the angles leads to the 3D-Rayleigh distribution in absence of the source, also called the Maxwell-Boltzmann distribution in Physics

whereas if the source is present, with amplitude , one obtains the distribution

If our image is pixelised, the different data , follow the two distributions

being , the null hypothesis (absence of source) and the alternative (presence of source), respectively and the profile at the pixel. The log-likelihood is defined by

The maximum likelihood estimator of the amplitude, , is given by the solution of the equation

This equation can be interpreted as a non-linear filter operating on the data that we shall also call the Neyman-Pearson filter (NPF).

#### Filtered fusion (FF)

In this case we use the same MF operating over the three images , , , respectively, as given by equation (Equation 6). Then, with the three filtered images , , we make the non-linear fusion pixel by pixel.

## 3Simulations

The European Space Agency (ESA) *Planck* satellite [33], to be launched in 2009, will measure the anisotropies of the CMB with unprecedented accuracy and angular resolution. It will also analyse the polarization of this primordial light. It is of great interest the detection and estimation of polarized sources in CMB maps [35]; since this radiation is linearly polarized [18], , we will apply the methods for the detection/estimation of the modulus of a 2-vector presented in Section 2.1. However, some cosmological models predict a possible circular polarization of CMB radiation [7]. Even if CMB is not circularly polarized, the extragalactic radio sources can indeed show circular polarization [2]. Besides, circular polarization occurs in many other astrophysical areas, from Solar Physics [36] to interestellar medium [8], just to put a few examples. In all these cases the results for the modulus of a 3-vector presented in Section 2.2 could be useful.

In order to compare and evaluate the performance of the filters presented in Section 2.1, we have carried out two-dimensional simulations with the characteristics of the 70 GHz *Planck* channel [35]. The simulated images have pixels with a pixel size of arcmin. We simulate the and components of the linear polarization as follows: each component consists of Gaussian uncorrelated detector noise plus the contribution of the CMB and a polarized point source filtered with a Gaussian-shaped beam whose full width at half maximum (FWHM) is , (the FWHM of the 70 GHz *Planck* channel beam). So the source polarization components can be written as

where is the beam dispersion (size) and we assume that the source is centred at the origin. The CMB simulation is based on the observed WMAP low multipoles, and on a Gaussian realisation assuming the WMAP best-fit at high multipoles. We do not include other possible foregrounds, since we are doing a first approximation to the detection/estimation problem and we assume that we apply our filters to relatively clean areas far away enough from the Galactic plane. Alternatively, we could consider a case in which the majority of foregrounds have been previously removed by means of some component separation method.

We consider realistic non-stationary detector noise. We have simulated the noise sky pattern for a *Planck* flight duration of 14 months, assuming a simple cycloidal scanning strategy with a 7 degree slow variation in the Ecliptic collatitude of the -axis. This scanning strategy implies that the sky will be covered inhomogeneously. The simulations have the same characteristics as the ones used in [4]. In order to illustrate the effects of non-stationary noise we have chosen two representative zones of the sky: one zone of high noise but quite isotropic and another zone of low noise but more anisotropic. The first zone is located in a region far from the Ecliptic poles, where the noise pattern is quite uniform and the number of hits per pixel of the detector is small. The average r.m.s deviation of the first zone (high noise) in units of (thermodynamic) is (for each component and ) and its standard deviation is . For this particular scanning strategy approximately of the sky has this kind of noise pattern. The second zone is close to one of the Ecliptic poles, where the sky is scanned more times (low noise level) but the hit pattern is very inhomogeneous. It has an average r.m.s. deviation and its standard deviation is . Then, in the second zone the noise is lower but proportionally more anisotropic. For this scanning strategy, of the sky has a noise pattern with these characteristics. Other zones of the sky would be intermediate cases between those considered here.

We take values of the source fluxes (before filtering with the Gaussian beam) and ranging from Jy to Jy with a step of Jy. A flux of 0.1 Jy corresponds for the 70 GHz channel to , so that it is of the order of the detector noise r.m.s deviation. The fluxes also correspond to typical polarization fluxes [37]. As we will see in Section 4, for some low-noise cases it is necessary to simulate even lower flux sources in order to study the behaviour of the filters in the low signal to noise ratio regime. We will explicitly report on this in the appropriate cases in Section 4. The number of simulations is for each combination of pairs of values of and . After carrying out the corresponding simulations for and , we add them quadratically and take the square root to calculate , the polarization modulus.

We study two different detection types: blind detection and non-blind detection. In the former case, we assume that we do not know the position of the source and then we place it at random in the image, in the latter case we know the source position and then we place it at the centre of the patch. We have considered images of pixels in order to do fast calculations. In order to avoid border effects, we simulate and filter pixel patches and, after the filtering step, we retain only the pixel central square. We have tested the case of larger patches but the obtained results do not change significantly.

We use two different filters: the filtered fusion (FF), which consists in the application of the matched filter to the images in and separately in a first moment, and then the calculation/construction of from the matched-filtered images and and the Neyman-Pearson filter (NPF), applied directly on , derived from the Neyman-Pearson lemma and presented with detail in Section 2. These filters are suitable for the case of uncorrelated noise, but we apply them to simulations including the CMB, which is correlated. However, we have checked that our results are similar when we consider simulations with and without the CMB: the relative differences of the errors in the estimated fluxes and positions are at most of a few percent. This is due to the low value of the CMB r.m.s. deviation as compared to that of the detector noise. Hence, the methods derived in Section 2.1 are also suitable for simulations including the CMB in the 70 GHz *Planck* case we consider.

In the blind case we apply these filters to each simulation, centering the filters successively at each pixel, since we do not know the source position. We estimate the source amplitude for the NPF, in this case we calculate the value of which maximises the log-likelihood, eqs. (Equation 4) and (Equation 5). For the FF we estimate separately and and obtain .

We compute the absolute maximum of in each filtered map and keep this value as the estimated value of the polarization of the source and the position of the maximum as the position of the source. Note that in the more realistic case where more than one source can be present in the images, it is still possible to proceed as described by looking for local peaks in the image. In the non-blind case we only centre the filters in the pixels included in one FWHM of the source centre (approximately of the total). In this way, we use the knowledge of the source position, then we calculate the absolute maximum of the estimated in these pixels. We also calculate the significance level of each detection. In order to do this, we carry out 1000 simulations with and , and we calculate the estimated value of the source polarization in this case for each filter. We consider the null hypothesis, , there is no polarized source, against the alternative hypothesis, , there is a polarized source. We set a significance level , this means that we reject the null hypothesis when a simulation has an estimated source amplitude higher than that of of the simulations without polarized source. We define the power of the test as , with the probability of accepting the null hypothesis when it is false, i.e. the power is the proportion of simulations with polarized source with an estimated amplitude higher than that of of the simulations without source. The higher the power the more efficient the filter is for detection. Note that the test can be performed in the same way in the blind and non-blind cases: in the second case we know the position of the source, but we do not know whether it is polarized or not.

## 4Results

### 4.1Blind case

We carry out simulations in the blind case for the high noise and low noise zones as explained above. We apply the filters to the images and calculate the absolute maximum of for each filtered image. If the detection has a significance higher than we consider it as a real detection, otherwise it would be a spurious one. For this significance level, we calculate the power of the detection test for the different pairs of and values. We also calculate the estimated value of the polarization amplitude , convert it to the estimated flux in Jy, , and compare it with the real value . We plot against for the NPF and the FF with high and low noise in Figure 1. We also compute the relative error of this estimation and its absolute value. For each pair of and values, we obtain the average and confidence intervals of the absolute value of the relative error, taking into account the performed simulations. These values are plotted against in Figure 2 for the same cases shown in Figure 1. We also calculate the estimated position of the detected source and obtain the position error expressed in terms of the number of pixels. The power and the average of the relative error in , of its absolute value and of the position error are presented in Table 1 for the two different types of noise. The rows in the Table are sorted in ascending order of . For the high noise case, we see in the Table that the power for the FF is higher than for the NPF. The position and polarization errors are also lower for the FF, this can also be seen in Figs. Figure 1 and Figure 2. For Jy, the FF and NPF perform similarly. For fluxes higher than Jy, the power with the FF is , the average relative error (bias) is , the average of the absolute value of the relative error is and the average position error is . From now on, we will use the flux limit for which the power is as a measure of the filter performance.

The results for the low noise case are logically much better, the FF also performs better than the NPF. The power of the two filters quickly reaches 100 for fluxes Jy and the errors in both flux and position remain stable from fluxes Jy. We therefore have cut the Table at Jy. In order to have a better sampling of the interesting signal to noise regime, we have simulated (using the same number of simulations as in the other cases) in addition the flux pairs , , , and Jy. This way the Table is much more informative. For fluxes higher than Jy, the power with the FF is , the average bias is , the average of the absolute value of the relative error is and the average position error is .

The FF is also much faster than the NPF. For instance the analysis of an image of 64x64 pixels in a personal computer takes 7 seconds with the FF and about 6 minutes with the NPF. This is due to the maximization involved. The computation time grows proportionally to the number of pixels, so that the NPF could be too slow if we want to analyze large images. This is the reason why we have considered small patches.

Finally, we show in Figure 3 four images corresponding to a polarized source with embedded in high noise. For the sake of a better visualisation, we show pixel images instead of the sized images used in the simulations. We show the original image in including noise, CMB and source, the image of the source, the image filtered with the NPF and the image treated with the FF method. By simple visual inspection, it is possible to see that the the FF performs better than the NPF.

### 4.2Non-blind case

We carry out simulations with a polarized source placed in the image centre and filter the images with the two different filters, centering them in the pixels at a distance from the source less than one FWHM. We simulate the same range of pairs as in the blind case (section Section 4.1). Then, as in the blind case, we calculate the maximum in , using only the selected pixels. This maximum is the estimated polarization of the source and the pixel where we find the maximum is the source position, we convert this amplitude to the estimated flux in Jy, . We have also performed a detection test, accepting as real sources only those detected with a significance lower than . In Table 2, we write the power of the test, the average of the relative polarization error, of its absolute value and of the position error. Our results are obtained from 500 simulations for each combination of pairs with high and low noise. The detection power is higher for the FF than for the other filter. For the high noise case, the improvement is very clear for Jy; we also obtain higher powers for the non-blind case than for the blind one. This could be expected, since we know the source position in the former case.

The position and polarization errors are also lower for the FF, this can also be seen in Figs. Figure 4 and Figure 5, which show the same quantities as Figs. Figure 1 and Figure 2 for non-blind detection. For fluxes higher than Jy, the power with the FF is , the bias is ,the average of the absolute value of the relative error is and the average position error is . The flux and position errors are lower in the non-blind case than in the blind one, especially for low fluxes.

The results for the low noise case are much better, the FF outperforms the NPF. From Jy on the results are quite similar. For fluxes higher than Jy, the power with the FF is , the average bias is , the average of the absolute value of the relative error is and the average position error is . Similarly to what happened with the low noise case of Table 1, we have cut the Table at flux Jy and we have added new cases in order to have a better sampling of the interesting signal to noise regime. In this case, we have simulated the pairs , , , , and Jy. The non-blind case gives better results than the blind one, especially for low fluxes.

## 5Conclusions

In this paper, we deal with the detection and estimation of the modulus of a vector, a problem of great interest in general and in particular in astrophysics when we consider the polarization of the cosmic microwave background (CMB). The polarization intensity is defined as , where , and are the Stokes parameters. We consider the case of images in , and consisting of a compact source with a profile immersed in Gaussian uncorrelated noise. We intend to detect the source and estimate its amplitude in by using two different methods: a Neyman-Pearson filter (NPF) operating in and based on the maximisation of the corresponding log-likelihood and a filtered fusion (FF) procedure, i.e. the application of the MF on the images of , and and the combination of the corresponding estimates by making the non-linear fusion . We present the two filters in Section 2 for two-dimensional, , and three-dimensional vectors, deriving the corresponding expressions for the estimation of the polarized source amplitude.

Since we are interested in applying these different filters to the CMB and this radiation is linearly polarized, we will only consider the two-dimensional vector case in our simulations.

Our goal is to compare the performance of the filters when applied to simulated images. Then, we carry out two-dimensional simulations with the characteristics of the 70 GHz *Planck* channel. The images have pixels with a pixel size of . We simulate the and components consisting of a compact source filtered with a Gaussian-shaped beam (FWHM of 14’) plus CMB and non-stationary detector noise. We consider two typical zones of the *Planck* survey: one with high noise and quite isotropic and another one with low noise but proportionally more anisotropic. These zones are extreme cases for the assumed scanning strategy we have chosen, and any other zone of the sky is an intermediate case between these two.

We study two types of detection: blind detection, in which we do not know the source position and non-blind, in which the position is known; we place the source at the centre of the image. We take values of the source fluxes in and , , and , ranging from 0.1 Jy to 0.5 Jy with a step of 0.1 Jy. Note that for extragalactic objects both and can take negative values, but since the sign of both components is irrelevant for the calculation of here we only simulate the positive case. We carry out 500 simulations for each pair and for all the cases of blind and non-blind detection with high and low noise.

We apply the filters to the simulated images, estimating the source amplitude as the maximum value of in the filtered images and the source position as the position of this maximum, then we convert the source amplitude to the source flux in Jy . We fix a significance for the detection and calculate a detection power for this significance, see Section 3 for more details. We also calculate the relative error of the estimated flux, its absolute value and the position error (in number of pixels), these errors together with the detection power are written in Tables Table 1 and Table 2 for the blind and non-blind case. We also show the estimated fluxes and the absolute value of the relative errors in Figs. Figure 1, Figure 2, Figure 4 and Figure 5 for the different cases.

The FF performs better than the NPF (as can be seen in the Tables and Figures), especially for low fluxes and it is also much faster than the NPF. However NPF is still interesting in a case where only the modulus of a vector is known and not its components (for example, if we had a map of polarization but not the and maps separately). The powers are much higher and the errors much lower for the low noise case than for the high noise one. The filters also perform better in the non-blind case than in the blind one, especially for low fluxes.

We can detect extragalactic point sources in polarization images (at 100% power) with the FF in the high noise zone with fluxes Jy for (blind/non-blind) detection and in the low noise zone with fluxes Jy for (blind/non-blind) detection. The bias and the position error are very low and pixel, respectively for all these fluxes.

High noise (,;) |
pow |
pow |
err |
err |
err |
err |
pos |
pos |
---|---|---|---|---|---|---|---|---|

(0.00,0.10;0.10) |
5 | 9 | 3.00 | 1.65 | 3.00 | 1.65 | 7.84 | 4.31 |

(0.10,0.10;0.14) |
8 | 17 | 1.87 | 0.97 | 1.87 | 0.97 | 7.97 | 2.74 |

(0.00,0.20;0.20) |
9 | 43 | 1.06 | 0.45 | 1.06 | 0.45 | 4.57 | 1.47 |

(0.10,0.20;0.22) |
13 | 54 | 0.86 | 0.32 | 0.86 | 0.32 | 4.40 | 1.13 |

(0.20,0.20;0.28) |
22 | 81 | 0.48 | 0.13 | 0.48 | 0.16 | 2.32 | 0.88 |

(0.00,0.30;0.30) |
29 | 91 | 0.43 | 0.12 | 0.43 | 0.16 | 1.73 | 0.73 |

(0.10,0.30;0.32) |
37 | 94 | 0.35 | 0.10 | 0.35 | 0.15 | 2.23 | 0.63 |

(0.20,0.30;0.36) |
52 | 99 | 0.22 | 0.05 | 0.22 | 0.13 | 1.26 | 0.48 |

(0.00,0.40;0.40) |
64 | 99 | 0.14 | 0.02 | 0.15 | 0.13 | 1.11 | 0.46 |

(0.10,0.40;0.41) |
73 | 99 | 0.11 | 0.01 | 0.12 | 0.11 | 1.18 | 0.40 |

(0.30,0.30;0.42) |
78 | 100 | 0.10 | 0.03 | 0.13 | 0.12 | 0.92 | 0.36 |

(0.20,0.40;0.45) |
87 | 100 | 0.07 | 0.03 | 0.11 | 0.11 | 0.97 | 0.28 |

(0.00,0.50;0.50) |
93 | 100 | 0.02 | 0.01 | 0.10 | 0.10 | 0.59 | 0.28 |

(0.30,0.40;0.50) |
97 | 100 | 0.03 | 0.02 | 0.11 | 0.10 | 0.67 | 0.29 |

(0.10,0.50;0.51) |
95 | 100 | 0.02 | 0.01 | 0.11 | 0.10 | 0.58 | 0.27 |

(0.20,0.50;0.54) |
98 | 100 | 0.02 | 0.01 | 0.10 | 0.09 | 0.51 | 0.17 |

(0.40,0.40;0.57) |
99 | 100 | 0.01 | 0.00 | 0.10 | 0.09 | 0.39 | 0.17 |

(0.30,0.50;0.58) |
100 | 100 | 0.01 | 0.01 | 0.10 | 0.09 | 0.41 | 0.19 |

(0.40,0.50;0.64) |
100 | 100 | 0.01 | 0.01 | 0.09 | 0.08 | 0.27 | 0.10 |

(0.50,0.50;0.71) |
100 | 100 | 0.00 | 0.00 | 0.08 | 0.07 | 0.20 | 0.10 |

Low noise (,;) |
pow |
pow |
err |
err |
err |
err |
pos |
pos |

(0.00,0.10;0.10) |
5 | 39 | 1.01 | 0.37 | 1.01 | 0.37 | 4.53 | 1.57 |

(0.05,0.10;0.11) |
9 | 53 | 0.81 | 0.23 | 0.81 | 0.23 | 5.04 | 1.16 |

(0.10,0.10;0.14) |
17 | 91 | 0.45 | 0.08 | 0.45 | 0.13 | 2.54 | 0.55 |

(0.00,0.15;0.15) |
19 | 93 | 0.39 | 0.07 | 0.39 | 0.13 | 1.94 | 0.48 |

(0.05,0.15;0.16) |
27 | 97 | 0.32 | 0.05 | 0.32 | 0.13 | 2.13 | 0.45 |

(0.10,0.15;0.18) |
55 | 99 | 0.16 | 0.03 | 0.16 | 0.12 | 1.30 | 0.34 |

(0.00,0.20;0.20) |
79 | 99 | 0.08 | 0.03 | 0.10 | 0.10 | 0.71 | 0.24 |

(0.10,0.20;0.22) |
92 | 100 | 0.04 | 0.01 | 0.09 | 0.09 | 0.45 | 0.21 |

(0.00,0.25;0.25) |
99 | 100 | 0.00 | 0.00 | 0.08 | 0.08 | 0.43 | 0.16 |

(0.20,0.20;0.28) |
100 | 100 | 0.01 | 0.01 | 0.08 | 0.07 | 0.18 | 0.09 |

(0.00,0.30;0.30) |
100 | 100 | 0.00 | 0.00 | 0.08 | 0.07 | 0.20 | 0.10 |

(0.10,0.30;0.32) |
100 | 100 | 0.01 | 0.01 | 0.07 | 0.06 | 0.15 | 0.08 |

(0.20,0.30;0.36) |
100 | 100 | 0.00 | 0.00 | 0.06 | 0.06 | 0.10 | 0.05 |

(0.00,0.40;0.40) |
100 | 100 | 0.00 | 0.00 | 0.06 | 0.05 | 0.05 | 0.03 |

High noise (,;) |
pow |
pow |
err |
err |
err |
err |
pos |
pos |
---|---|---|---|---|---|---|---|---|

(0.00,0.10;0.10) |
9 | 19 | 2.53 | 1.33 | 2.53 | 1.33 | 1.90 | 1.52 |

(0.10,0.10;0.14) |
13 | 35 | 1.52 | 0.74 | 1.52 | 0.74 | 1.70 | 1.37 |

(0.00,0.20;0.20) |
22 | 68 | 0.84 | 0.32 | 0.84 | 0.32 | 1.53 | 1.09 |

(0.10,0.20;0.22) |
32 | 78 | 0.66 | 0.23 | 0.66 | 0.24 | 1.44 | 0.90 |

(0.20,0.20;0.28) |
52 | 95 | 0.35 | 0.10 | 0.35 | 0.17 | 1.26 | 0.81 |

(0.00,0.30;0.30) |
53 | 99 | 0.28 | 0.06 | 0.28 | 0.16 | 1.18 | 0.67 |

(0.10,0.30;0.32) |
63 | 98 | 0.26 | 0.07 | 0.26 | 0.15 | 1.13 | 0.66 |

(0.20,0.30;0.36) |
77 | 100 | 0.14 | 0.03 | 0.16 | 0.13 | 0.92 | 0.55 |

(0.00,0.40;0.40) |
90 | 100 | 0.08 | 0.03 | 0.14 | 0.13 | 0.79 | 0.38 |

(0.10,0.40;0.41) |
92 | 100 | 0.06 | 0.02 | 0.13 | 0.12 | 0.77 | 0.39 |

(0.30,0.30;0.42) |
94 | 100 | 0.06 | 0.03 | 0.13 | 0.12 | 0.77 | 0.39 |

(0.20,0.40;0.45) |
97 | 100 | 0.04 | 0.02 | 0.12 | 0.11 | 0.72 | 0.40 |

(0.00,0.50;0.50) |
99 | 100 | 0.02 | 0.02 | 0.11 | 0.10 | 0.58 | 0.25 |

(0.30,0.40;0.50) |
100 | 100 | 0.02 | 0.02 | 0.11 | 0.10 | 0.59 | 0.29 |

(0.10,0.50;0.51) |
99 | 100 | 0.02 | 0.01 | 0.11 | 0.09 | 0.59 | 0.30 |

(0.20,0.50;0.54) |
100 | 100 | 0.02 | 0.02 | 0.11 | 0.09 | 0.49 | 0.22 |

(0.40,0.40;0.57) |
100 | 100 | 0.01 | 0.00 | 0.10 | 0.09 | 0.40 | 0.18 |

(0.30,0.50;0.58) |
100 | 100 | 0.01 | 0.01 | 0.10 | 0.09 | 0.43 | 0.16 |

(0.40,0.50;0.64) |
100 | 100 | 0.01 | 0.01 | 0.09 | 0.08 | 0.30 | 0.14 |

(0.50,0.50;0.71) |
100 | 100 | 0.01 | 0.01 | 0.08 | 0.08 | 0.23 | 0.11 |

Low noise (,;) |
pow |
pow |
err |
err |
err |
err |
pos |
pos |

(0.00,0.10;0.10) |
12 | 69 | 0.77 | 0.20 | 0.77 | 0.21 | 1.44 | 0.78 |

(0.05,0.10;0.11) |
22 | 88 | 0.58 | 0.15 | 0.58 | 0.18 | 1.52 | 0.72 |

(0.00,0.12;0.12) |
24 | 93 | 0.47 | 0.10 | 0.47 | 0.15 | 1.18 | 0.65 |

(0.00,0.13;0.13) |
34 | 95 | 0.38 | 0.09 | 0.38 | 0.16 | 1.14 | 0.57 |

(0.10,0.10;0.14) |
48 | 98 | 0.26 | 0.06 | 0.26 | 0.14 | 0.92 | 0.51 |

(0.00,0.15;0.15) |
57 | 99 | 0.21 | 0.05 | 0.21 | 0.14 | 0.82 | 0.52 |

(0.10,0.15;0.18) |
87 | 100 | 0.07 | 0.01 | 0.11 | 0.11 | 0.61 | 0.39 |

(0.00,0.20;0.20) |
94 | 100 | 0.03 | 0.02 | 0.12 | 0.11 | 0.56 | 0.28 |

(0.10,0.20;0.22) |
99 | 100 | 0.03 | 0.02 | 0.10 | 0.09 | 0.39 | 0.23 |

(0.00,0.25;0.25) |
100 | 100 | 0.02 | 0.01 | 0.09 | 0.08 | 0.32 | 0.15 |

(0.20,0.20;0.28) |
100 | 100 | 0.01 | 0.01 | 0.08 | 0.07 | 0.22 | 0.11 |

(0.00,0.30;0.30) |
100 | 100 | 0.01 | 0.01 | 0.08 | 0.07 | 0.21 | 0.09 |

(0.10,0.30;0.32) |
100 | 100 | 0.00 | 0.00 | 0.08 | 0.07 | 0.16 | 0.08 |

(0.20,0.30;0.36) |
100 | 100 | 0.00 | 0.01 | 0.06 | 0.06 | 0.09 | 0.04 |

(0.00,0.40;0.40) |
100 | 100 | 0.00 | 0.01 | 0.05 | 0.05 | 0.08 | 0.04 |

## Acknowledgements

The authors acknowledge partial financial support from the Spanish Ministry of Education (MEC) under project ESP2004-07067-C03-01 and from the joint CNR-CSIC research project 2006-IT-0037. JLS acknowledges partial financial support by the Spanish MEC and thanks the CNR ISTI in Pisa for their hospitality during his sabbatical leave. FA also wishes to thank the CNR ISTI in Pisa for their hospitality during two short stays when part of this work was done. MLC acknowledges the Spanish MEC for a postdoctoral fellowship. Some of the results in this paper have been derived using the HEALPix [10] package.

### References

- Agarwal N., Jain P., McKay D. W., Ralston J. P., 2008, Physical Review D, 78, 085028
- Aller M. F., Aller H. D., Plotkin R. M., 2005, in Bulletin of the American Astronomical Society Vol. 37 of Bulletin of the American Astronomical Society, Circular Polarization in AGNs: Polarity and Spectra. pp 1371–+
- Argüeso F., Sanz J. L., 2008, in Proceedings of the 16th European Signal Processing Conference (2008) EUSIPCO 2008 Conference, Filter design for the detection of compact sources embedded in non-stationary noise plus a deterministic background. Lausanne, Switzerland, pp 1–5
- Argüeso F., Sanz J. L., Barreiro R. B., Herranz D., González-Nuevo J., 2006, MNRAS, 373, 311
- Barnard V. E., Vielva P., Pierce-Price D. P. I., Blain A. W., Barreiro R. B., Richer J. S., Qualtrough C., 2004, MNRAS, 352, 961
- Carvalho P., Rocha G., Hobson M. P., 2009, MNRAS, 393, 681
- Cooray A., Melchiorri A., Silk J., 2003, Physics Letters B, 554, 1
- Cox N. L. J., Boudin N., Foing B. H., Schnerr R. S., Kaper L., Neiner C., Henrichs H., Donati J.-F., Ehrenfreund P., 2007, A&A, 465, 899
- González-Nuevo J., Argüeso F., López-Caniego M., Toffolatti L., Sanz J. L., Vielva P., Herranz D., 2006, MNRAS, 369, 1603
- Górski K. M., Hivon E., Banday A. J., Wandelt B. D., Hansen F. K., Reinecke M., Bartelmann M., 2005, ApJ, 622, 759
- Herranz D., Gallegos J., Sanz J. L., Martínez-González E., 2002a, MNRAS, 334, 533
- Herranz D., Sanz J. L., Hobson M. P., Barreiro R. B., Diego J. M., Martínez-González E., Lasenby A. N., 2002b, MNRAS, 336, 1057
- Herranz D., Sanz J. L., Barreiro R. B., López-Caniego M., 2005, MNRAS, 356, 944
- Herranz D., Sanz J. L., 2008, IEEE J. Sel. Topics in Signal Processing, 2(5), 727
- Herranz D., López-Caniego M., Sanz J. L., González-Nuevo J., 2009, MNRAS, 394, 510
- Hobson M. P., McLachlan C., 2003, MNRAS, 338, 765
- Homan D. C., Lister M. L., Aller H. D., Aller M. F., Wardle J. F., 2006, in Bulletin of the American Astronomical Society Vol. 38 of Bulletin of the American Astronomical Society, Multi-frequency VLBA Observations of Circular Polarization from Extragalactic Radio Jets. pp 904–+
- Kamionkowski M., Kosowsky A., Stebbins A., 1997, Phys. Rev. D, 55, 7368
- Kirk J. G., Tsang O., 2006, A&A, 447, L13
- Leach S. M., Cardoso J. F., Baccigalupi C., Barreiro R. B., Betoule M., Bobin J., Bonaldi A., de Zotti G., Delabrouille J., Dickinson C., 20 coauthors more, 2008, A&A, 481, 597
- López-Caniego M., Herranz D., Barreiro R. B., Sanz J. L., 2005a, MNRAS, 359, 993
- López-Caniego M., Herranz D., Sanz J. L., Barreiro R. B., 2005b, EURASIP Journal on Applied Signal Processing, 15, 2426
- López-Caniego M., Herranz D., González-Nuevo J., Sanz J. L., Barreiro R. B., Vielva P., Argüeso F., Toffolatti L., 2006, MNRAS, 370, 2047
- López-Caniego M., González-Nuevo J., Herranz D., Massardi M., Sanz J. L., De Zotti G., Toffolatti L., Argüeso F., 2007, ApJS, 170, 108
- Malik R. K., Subramanian K., 1997, A&A, 317, 318
- Nailong W., 1992, in Worrall D. M., Biemesderfer C., Barnes J., eds, Astronomical Data Analysis Software and Systems I Vol. 25 of Astronomical Society of the Pacific Conference Series, Using a matched filter to improve SNR of radio maps.. pp 291–293
- Papoulis A., 1984, Probability, Random Variables, and Stochastic Processes. Mc-Graw Hill
- Reiner M. J., Fainberg J., Kaiser M. L., Bougeret J.-L., 2007, Solar Physics, 241, 351
- Rice S. O., 1954, Selected Papers on Noise and Stochastic Processes. Dover Pubns, pp 133–294
- Sanz J. L., Herranz D., López-Caniego M., Argüeso F., 2006, in Proceedings of the 14th European Signal Processing Conference (2006). EUSIPCO 2006 Conference, Wavelets on the sphere. Application to the detection problem. pp 1–5
- Sanz J. L., Herranz D., Martínez-Gónzalez E., 2001, ApJ, 552, 484
- Stewart I. M., 2006, A&A, 454, 997
- Tauber J. A., 2005, in Lasenby A. N., Wilkinson A., eds, New Cosmological Data and the Values of the Fundamental Parameters Vol. 201 of IAU Symposium, The Planck Mission. pp 86–+
- Tegmark M., de Oliveira-Costa A., 1998, ApJL, 500, L83+
- The
*Planck*Collaboration 2006, The Scientific Programme of Planck - Tritschler A., Müller D. A. N., Schlichenmaier R., Hagenaar H. J., 2007, ApJL, 671, L85
- Tucci M., Martínez-González E., Toffolatti L., González-Nuevo J., De Zotti G., 2004, MNRAS, 349, 1267
- Tucci M., Martínez-González E., Vielva P., Delabrouille J., 2005, MNRAS, 360, 935
- Vielva P., Barreiro R. B., Hobson M. P., Martínez-González E., Lasenby A. N., Sanz J. L., Toffolatti L., 2001, MNRAS, 328, 1
- Vielva P., Martínez-González E., Gallegos J. E., Toffolatti L., Sanz J. L., 2003, MNRAS, 344, 89
- Vikhlinin A., Forman W., Jones C., Murray S., 1995, ApJ, 451, 542