+ View the NASA Portal

Search JPL

Pixel window functions

A pixelised signal f(p) is the average within each pixel p (with surface area ${\Omega_{\rm pix}}$ ) of the underlying signal

f(p)

$\textstyle \myequal$

$\displaystyle \int d\u w_p(\u )f(\u )$

(26)

where w_p is equal to $1/{\Omega_{\rm pix}}$ within the pixel, and equal to 0 outside, so that $\int d\u w_p(\u ) = 1$ . Eq. (1) then becomes

f(p)

$\textstyle \myequal$

$\displaystyle \sum_{l=0}^{l_{max}}\sum_{m}a_{lm}w_{lm}(p),$

(27)

where

w_lm(p)

$\textstyle \myequal$

$\displaystyle \int d\u w_p(\u ) Y_{lm}(\u ),$

(28)

is the Spherical Harmonic Transform of the pixel p.

However, complete analysis of a pixelised map with the exact w_lm(p) defined above would be computationally intractable (because of azimutal variation of pixel shape over the polar caps of the HEALPix grid), and some simplifying asumptions have to be made. If the pixel is small compared to the signal correlation length (determined by the beam size), the exact structure of the pixel can be ignored in the subsequent analysis and we can assume

w_lm(p) = w_l(p) Y_lm(p)

(29)

where we introduced the m-averaged window function

w_l(p)

$\textstyle \myequal$

$\displaystyle \left(\frac{4 \pi}{2l+1}\sum_{m=-l}^{l} \left\vert w_{lm}(p)\right\vert^2\right)^{1/2},$

(30)

which is independent of the pixel location on the sky.

If we assume all the pixels to be identical, the power spectrum of the pixelized map, $C_{l}^{\rm pix}$ , is related to the hypothetical unpixelized one, $C_{l}^{\rm unpix}$ , by

$\displaystyle C_{l}^{\rm pix}$

$\textstyle \myequal$

$\displaystyle w^2_{l} C_{l}^{\rm unpix}$

(31)

where the effective pixel window function w_l is defined as

w_l

$\textstyle \myequal$

$\displaystyle \left(\frac{1}{{N_{\rm pix}}}\sum_{p=0}^{{N_{\rm pix}}-1} w^2_{l}(p)\right)^{1/2}.$

(32)

This function is provided with the HEALPix package for $l\le 4{N_{\rm side}}$ for each resolution parameter ${N_{\rm side}}$ .

The pixel window functions are now available for both temperature and polarization.

For ${N_{\rm side}}\le 128$ , those window functions are computed exactly using Eqs. (30) and (32). For ${N_{\rm side}}> 128$ the calculations are too costly to be done exactly at all l. The temperature windows are extrapolated from the case ${N_{\rm side}}= 128$ assuming a scaling in l similar to the one exhibited by the window of a tophat pixel. The polarization windows are assumed to be proportional to those for temperature, with a proportionality factor given by the exact calculation of w_l at low l.

Because of a change of the extrapolation scheme used, the temperature window functions provided with HEALPix 1.2 and higher for ${N_{\rm side}}> 128$ are slighty different from those provided with HEALPix 1.1. For a given ${N_{\rm side}}$ , the relative difference increases almost linearly with l, and is of the order of $\Delta w/w < 7\ 10^{-4}$ at $l=2{N_{\rm side}}$ and $\Delta w/w < 1.7\ 10^{-3}$ at $l=4{N_{\rm side}}$ .

Eric Hivon 2010-06-18

Contact:

Site Manager:
Webmaster:
CL 03-2650