HEALPix

Last updated
HEALPix H=4, K=3 projection of the world. The lines on the map are a graticule of latitudes and longitudes. HEALPix projection SW.svg
HEALPix H=4, K=3 projection of the world. The lines on the map are a graticule of latitudes and longitudes.
The grid used by HEALPix and its subdivision of the sphere, in four different grid refinements. Notice the similarity between the coarsest grid and the rhombic dodecahedron. HealpixGridRefinement.jpg
The grid used by HEALPix and its subdivision of the sphere, in four different grid refinements. Notice the similarity between the coarsest grid and the rhombic dodecahedron.

HEALPix (sometimes written as Healpix), an acronym for Hierarchical Equal Area isoLatitude Pixelisation of a 2-sphere, is an algorithm for pixelisation of the 2-sphere based on subdivision of a distorted rhombic dodecahedron, and the associated class of map projections. [1] The pixelisation algorithm was devised in 1997 by Krzysztof M. Górski at the Theoretical Astrophysics Center in Copenhagen, Denmark, [2] and first published as a preprint in 1998. [3] [4]

Contents

Projection and pixelisation

The HEALPix projection is a general class of spherical projections, sharing several key properties, which map the 2-sphere to the Euclidean plane. [1] Any of these can be followed by partitioning (pixelising) the resulting region of the 2-plane. In particular, when one of these projections (the H=4, K=3 HEALPix projection) is followed by a pixelisation of the 2-plane, the result is generally known as the HEALPix pixelisation, [3] [4] which is widely used in physical cosmology for maps of the cosmic microwave background. This pixelisation can be thought of as mapping the sphere to twelve square facets (diamonds) on the plane followed by the binary division of these facets into pixels, [5] [6] [1] though it can be derived without using the projection. [3] [4] [7] The associated software package HEALPix implements the algorithm. [3] [7] The HEALPix projection (as a general class of spherical projections) is represented by the keyword HPX in the FITS standard for writing astronomical data files. It was approved as part of the official FITS World Coordinate System (WCS) by the International Astronomical Union FITS Working Group on April 26, 2006. [8]

The spherical projection combines a cylindrical equal area projection, the Lambert cylindrical equal-area projection, for the equatorial regions of the sphere and a pseudocylindrical equal area projection, an interrupted Collignon projection, for the polar regions. [1]

At a given level in the hierarchy the pixels are of equal area (which is done by bisecting the square in the case of the H=4, K=3 projection) and their centers lie on a discrete number of circles of latitude, with equal spacing on each circle. The scheme has a number of mathematical properties which make it efficient for certain computations, e.g. spherical harmonic transforms. In the case of the H=4, K=3 projection, the pixels are squares in the plane (which can be inversely projected back to quadrilaterals with non-geodesic sides on the 2-sphere) and every vertex joins four pixels, with the exception of eight vertices which each join only three pixels.

The latitude of transition between equatorial-orthogonal and polar-convergent longitude lines has been selected to allow the folding of the projection into a perfect cube — "cubing the sphere"; indeed in this way the Arctic Circle becomes a square.

Usage and alternatives

The pixelisation related to the H=4, K=3 projection has become widely used in cosmology for storing and manipulating maps of the cosmic microwave background.

Gaia mission uses HEALPix as the basis for source identification. [9]

An alternative hierarchical grid is the Hierarchical Triangular Mesh (HTM). [10] [11] The pixels at a given level in the hierarchy are of similar but not identical size. The scheme is good at representing complex shapes because the boundaries are all segments of circles of the sphere. Another alternative hierarchical grid is the Quadrilateralized Spherical Cube.

The 12 "base resolution pixels" of H=4, K=3 HEALPix projection may be thought of as the facets of a rhombic dodecahedron.

The H=6 HEALPix has similarities to another alternative grid based on the icosahedron. [12]

See also

Related Research Articles

<span class="mw-page-title-main">Cosmic microwave background</span> Trace radiation from the early universe

The cosmic microwave background, or relic radiation, is microwave radiation that fills all space in the observable universe. With a standard optical telescope, the background space between stars and galaxies is almost completely dark. However, a sufficiently sensitive radio telescope detects a faint background glow that is almost uniform and is not associated with any star, galaxy, or other object. This glow is strongest in the microwave region of the radio spectrum. The accidental discovery of the CMB in 1965 by American radio astronomers Arno Penzias and Robert Wilson was the culmination of work initiated in the 1940s.

<span class="mw-page-title-main">Latitude</span> Geographic coordinate specifying north–south position

In geography, latitude is a coordinate that specifies the north–south position of a point on the surface of the Earth or another celestial body. Latitude is given as an angle that ranges from −90° at the south pole to 90° at the north pole, with 0° at the Equator. Lines of constant latitude, or parallels, run east–west as circles parallel to the equator. Latitude and longitude are used together as a coordinate pair to specify a location on the surface of the Earth.

<span class="mw-page-title-main">Universe</span> Everything in space and time

The universe is all of space and time and their contents. It comprises all of existence, any fundamental interaction, physical process and physical constant, and therefore all forms of matter and energy, and the structures they form, from sub-atomic particles to entire galactic filaments. Since the early 20th century, the field of cosmology establishes that space and time emerged together at the Big Bang 13.787±0.020 billion years ago and that the universe subsequently expanded. Today, the universe has expanded into an age and size that is only partially observable from Earth; while the spatial size of the entire universe is unknown, the smaller observable universe is approximately 93 billion light-years in diameter at present.

<span class="mw-page-title-main">Photosphere</span> Stars outer shell from which light is radiated

The photosphere is a star's outer shell from which light is radiated. It extends into a star's surface until the plasma becomes opaque, equivalent to an optical depth of approximately 23, or equivalently, a depth from which 50% of light will escape without being scattered.

<span class="mw-page-title-main">Map projection</span> Systematic representation of the surface of a sphere or ellipsoid onto a plane

In cartography, a map projection is any of a broad set of transformations employed to represent the curved two-dimensional surface of a globe on a plane. In a map projection, coordinates, often expressed as latitude and longitude, of locations from the surface of the globe are transformed to coordinates on a plane. Projection is a necessary step in creating a two-dimensional map and is one of the essential elements of cartography.

<span class="mw-page-title-main">Shape of the universe</span> Local and global geometry of the universe

In physical cosmology, the shape of the universe refers to both its local and global geometry. Local geometry is defined primarily by its curvature, while the global geometry is characterised by its topology. General relativity explains how spatial curvature is constrained by gravity. The global topology of the universe cannot be deduced from measurements of curvature inferred from observations within the family of homogeneous general relativistic models alone, due to the existence of locally indistinguishable spaces with varying global topological characteristics. For example; a multiply connected space like a 3 torus has everywhere zero curvature but is finite in extent, whereas a flat simply connected space is infinite in extent.

<span class="mw-page-title-main">Stereographic projection</span> Particular mapping that projects a sphere onto a plane

In mathematics, a stereographic projection is a perspective projection of the sphere, through a specific point on the sphere, onto a plane perpendicular to the diameter through the point. It is a smooth, bijective function from the entire sphere except the center of projection to the entire plane. It maps circles on the sphere to circles or lines on the plane, and is conformal, meaning that it preserves angles at which curves meet and thus locally approximately preserves shapes. It is neither isometric nor equiareal.

<span class="mw-page-title-main">Cosmological principle</span> Theory that the universe is the same in all directions

In modern physical cosmology, the cosmological principle is the notion that the spatial distribution of matter in the universe is uniformly isotropic and homogeneous when viewed on a large enough scale, since the forces are expected to act equally throughout the universe on a large scale, and should, therefore, produce no observable inequalities in the large-scale structuring over the course of evolution of the matter field that was initially laid down by the Big Bang.

Flexible Image Transport System (FITS) is an open standard defining a digital file format useful for storage, transmission and processing of data: formatted as multi-dimensional arrays, or tables. FITS is the most commonly used digital file format in astronomy. The FITS standard was designed specifically for astronomical data, and includes provisions such as describing photometric and spatial calibration information, together with image origin metadata.

<span class="mw-page-title-main">Lambda-CDM model</span> An anomaly in astronomical observations of the Cosmic Microwave Background

The Lambda-CDM, Lambda cold dark matter, or ΛCDM model is a mathematical model of the Big Bang theory with three major components:

  1. a cosmological constant, denoted by lambda (Λ), associated with dark energy
  2. the postulated cold dark matter, denoted by CDM
  3. ordinary matter
<span class="mw-page-title-main">Very Small Array</span> Radio telescope in the Canary Islands

The Very Small Array (VSA) was a 14-element interferometric radio telescope operating between 26 and 36 GHz that is used to study the cosmic microwave background radiation. It was a collaboration between the University of Cambridge, University of Manchester and the Instituto de Astrofisica de Canarias (Tenerife), and was located at the Observatorio del Teide on Tenerife. The array was built at the Mullard Radio Astronomy Observatory by the Cavendish Astrophysics Group and Jodrell Bank Observatory, and was funded by PPARC. The design was strongly based on the Cosmic Anisotropy Telescope.

<span class="mw-page-title-main">Reflection mapping</span> Technique in computer graphics to represent reflective surfaces

In computer graphics, reflection mapping or environment mapping is an efficient image-based lighting technique for approximating the appearance of a reflective surface by means of a precomputed texture. The texture is used to store the image of the distant environment surrounding the rendered object.

<span class="mw-page-title-main">Geodesic grid</span> Spatial grid based on a geodesic polyhedron

A geodesic grid is a spatial grid based on a geodesic polyhedron or Goldberg polyhedron.

<span class="mw-page-title-main">Hubble volume</span> Region of the observable universe

In cosmology, a Hubble volume (named for the astronomer Edwin Hubble) or Hubble sphere, Hubble bubble, subluminal sphere, causal sphere and sphere of causality is a spherical region of the observable universe surrounding an observer beyond which objects recede from that observer at a rate greater than the speed of light due to the expansion of the universe. The Hubble volume is approximately equal to 1031 cubic light years (or about 1079 cubic meters).

<span class="mw-page-title-main">Galaxy merger</span> Merger whereby at least two galaxies collide

Galaxy mergers can occur when two galaxies collide. They are the most violent type of galaxy interaction. The gravitational interactions between galaxies and the friction between the gas and dust have major effects on the galaxies involved, but the exact effects of such mergers depend on a wide variety of parameters such as collision angles, speeds, and relative size/composition, and are currently an extremely active area of research. Galaxy mergers are important because the merger rate is a fundamental measurement of galaxy evolution and also provides astronomers with clues about how galaxies grew into their current forms over long stretches of time.

<span class="mw-page-title-main">Quadrilateralized spherical cube</span> Polyhedral equal-area map projection

In mapmaking, a quadrilateralized spherical cube, or quad sphere for short, is an equal-area polyhedral map projection and discrete global grid scheme for data collected on a spherical surface. It was first proposed in 1975 by Chan and O'Neill for the Naval Environmental Prediction Research Facility. This scheme is also often called the COBE sky cube, because it was designed to hold data from the Cosmic Background Explorer (COBE) project.

In cosmology, intensity mapping is an observational technique for surveying the large-scale structure of the universe by using the integrated radio emission from unresolved gas clouds.

<span class="mw-page-title-main">Discrete global grid</span> Partition of Earths surface into subdivided cells

A discrete global grid (DGG) is a mosaic that covers the entire Earth's surface. Mathematically it is a space partitioning: it consists of a set of non-empty regions that form a partition of the Earth's surface. In a usual grid-modeling strategy, to simplify position calculations, each region is represented by a point, abstracting the grid as a set of region-points. Each region or region-point in the grid is called a cell.

<span class="mw-page-title-main">Hierarchical triangular mesh</span> Data structure for storing geometric information

Hierarchical Triangular Mesh (HTM) is a kind of quad tree based on subdivision of a distorted octahedron, used for mesh generation in 3-D computer graphics and geometric data structures.

References

  1. 1 2 3 4 Calabretta, Mark R.; Roukema, Boudewijn F. (2007). "Mapping on the HEALPix grid". MNRAS . 381 (2). Oxford University Press: 865–872. Bibcode:2007MNRAS.381..865C. doi: 10.1111/j.1365-2966.2007.12297.x .
  2. "HEALPix Background - History". healpix.jpl.nasa.gov. Retrieved 2019-06-08.
  3. 1 2 3 4 Górski, Krzysztof M.; Hivon, Éric; Wandelt, Benjamin D. (1999). "Analysis Issues for Large CMB Data Sets". Evolution of Large Scale Structure: From Recombination to Garching. Proceedings of the MAP/ESO Cosmology Conference 'Evolution of Large-Scale Structure'. Netherlands: PrintPartners Ipskamp. p. 37. arXiv: astro-ph/9812350 . Bibcode:1999elss.conf...37G.
  4. 1 2 3 Górski, Krzysztof M.; Wandelt, Benjamin D.; Hansen, Frode K.; Hivon, Éric; Banday, Anthony J. (1999-05-21). "The HEALPix Primer". arXiv: astro-ph/9905275 .
  5. Roukema, Boudewijn F.; Lew, Bartosz (2004-09-08). "A Solution to the Isolatitude, Equi-area, Hierarchical Pixel-Coordinate System". Public Draft. arXiv: astro-ph/0409533 . Bibcode:2004astro.ph..9533R. Archived from the original on 2019-08-04. Retrieved 2004-09-08.
  6. Roukema, Boudewijn F.; Lew, Bartosz (2004-09-22). "A Solution to the Isolatitude, Equi-area, Hierarchical Pixel-Coordinate System". arXiv: astro-ph/0409533 .
  7. 1 2 Górski, Krzysztof M.; Hivon, Éric; Banday, Anthony J.; Hansen, Frode K.; Wandelt, Benjamin D.; Reinecke, M.; Bartelmann, M. (2005). "HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere". Astrophysical Journal . 622 (2): 759–771. arXiv: astro-ph/0409513 . Bibcode:2005ApJ...622..759G. doi:10.1086/427976. S2CID   18743679.
  8. Pence, William D. "FITS World Coordinate System (WCS)". High Energy Astrophysics Science Archive Research Center (HEASARC). Archived from the original on 2019-08-04. Retrieved 2007-01-09.
  9. "Gaia Data Release 1: Datamodel description Documentation release 1.2". gea.esac.esa.int. Retrieved 2021-05-31.
  10. "SkyServer.org - HTM: Hierarchical Triangular Mesh". SkyServer. June 6, 2006. Retrieved 2007-02-05.
  11. Szalay, Alex; Jim Gray; Gyorgy Fekete; Peter Kunszt; Peter Kukol; Ani Thakar (September 2005). "Indexing the Sphere with the Hierarchical Triangular Mesh". Microsoft Research. arXiv: cs/0701164 . Bibcode:2007cs........1164S . Retrieved 2007-02-05.
  12. Tegmark, Max. "Welcome to the icosahedron home page". space.mit.edu.