CosmoCoffee

[Laurent Jacques]

Hi,

If there are some persons specialized in HEALpix on this forum, can somebody explain me why this "Hierarchical Equal Area isoLatitude Pixelisation" contains the two words "Equal Area" ?

IMHO, there is no truly equal area in this grid, excepted if you define a new definition of area on the sphere, as it is defined in the HEALpix original paper, i.e.
Dz * Dphi !

For me, this equal area concept is just a rough approximation. What do you think ?

Best,
Laurent.

[Anthony Banday]

I think you are wrong.

Please take a look at the recent paper available on the astro-ph
archive. You will find a prescription for the pixel boundaries there.
Formal integration will demonstrate that the pixels are exactly equal area.
This was a driving requirement in the development of HEALPix, and
great care was taken to ensure this condition was met.

Tony Banday

[Laurent Jacques]

If you speak about the paper astro-ph/0409513

Title: HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere. K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, M. Bartelman.

I didn't find any clues about the integration realized: Is-it done in (z=cos(theta),phi) coordinate or truly in (theta,phi) ? The former is not an area computing on the sphere by on the cylinder even if we take into account the non-geodesic pixel boundaries equation.

Do you know a paper where I can find a true proof of this "equal area" theorem ?
astro-ph/0409513 is too much summarized for that.

Best,
Laurent

[Antony Lewis]

HealPix is constructed by division of the sphere into 12 equal large pixels, then recursiverly diving each pixel into 4 equal size parts. So you get 12, 12 x 4, 12 x 4 x 4, 12 x 4x 4x 4, ... pixels at each higher devision, and by construction they are equal area. Just look at Fig 4 in the paper of how they are constructing higher resolutions.

[Eric Hivon]

Hi Laurent,

what is your problem with the dz * dphi integration ?
You know that the integration element on the sphere is
d omega = sin (theta) * d theta * d phi,
which is equal to d omega = d z * d phi if you define z = cos (theta).
Since -1 < z < 1, the total area available in (z,phi) is 4 Pi, llike on the sphere.

In other terms, the mapping from the sphere (theta, phi) to cylindrical coordinates (z, phi) preserves the area.
Therefore, any area measured in (z,phi) will be the same on the sphere.

As for the exact calculation, in the equatorial region,
the pixel boundaries have the form
zp = 2/3 - a*k + b*phi,
zm = 2/3 - a*k' - b*phi,
where a=4/3/Nside and b=8/3/Pi and k and k' are integers
(see http://arxiv.org/abs/astro-ph/0409513 Eq 22) .
The pixel boundaries are therefore straight lines (in (z,phi)) which are either increasing (zp) or decreasing (zm) with phi, and the pixel is a diamond, with its diagonals along
the North-South axis and East-West axis respectively (ie, crossing at a right angle).

A pixel enclosed in the lines zp and zm with the parameters k, k+1 and k', k'+1 respectively,
will have a distance Dz = a between its North and South corners,
and a distance Dphi = a/b between its East and West corners.

As you know, the surface of a diamond with orthogonal diagonals
if a half of the products of its diagonal. Therefore
S = Dz * Dphi / 2 = a * a / b / 2 = 2 Pi / (3 * Nside * Nside) = 4 Pi / Npix

In the polar regions, the pixel boundary equations are slightly more complicated
(too much for an e-mail),
requiring the integration to be done in details (it is easier to divide the pixel in 4 quadrants
along its diagonals, compute the 4 integrals and add them) to find a surface of
4 Pi / Npix per pixel.

Hope this helps,

Eric

[Laurent Jacques]

[Laurent Jacques]

[Boud Roukema]

[Boud Roukema]

i guess i forgot to say: there's still the question of whether or not the spherical projection is an equal area projection - e.g. the Mercator projection is clearly not an equal area projection.

In (1) the spherical projection for healpix:

- In the equatorial area, this is done by using sin(theta) instead of theta e.g. see eq. (7) of astro-ph/0409533 .

- In the polar area, you need to check that the partial derivatives give the right values - eqs. (29)-(31) give this in astro-ph/0409533 .

hope this helps.