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 is microwave radiation that fills all space in the observable universe. It is a remnant that provides an important source of data on the primordial 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">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.

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.

In modern physical cosmology, the cosmological principle is the notion that the spatial distribution of matter in the universe is equally distributed and isotropic when viewed on a large enough scale, since the forces are expected to act equally throughout the universes 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">BOOMERanG experiment</span> High-altitude balloon package measuring the universes geometry

In astronomy and observational cosmology, the BOOMERanG experiment was an experiment which measured the cosmic microwave background radiation of a part of the sky during three sub-orbital (high-altitude) balloon flights. It was the first experiment to make large, high-fidelity images of the CMB temperature anisotropies, and is best known for the discovery in 2000 that the geometry of the universe is close to flat, with similar results from the competing MAXIMA experiment.

<span class="mw-page-title-main">Dark matter halo</span> Theoretical cosmological structure

According to modern models of physical cosmology, a dark matter halo is a basic unit of cosmological structure. It is a hypothetical region that has decoupled from cosmic expansion and contains gravitationally bound matter. A single dark matter halo may contain multiple virialized clumps of dark matter bound together by gravity, known as subhalos. Modern cosmological models, such as ΛCDM, propose that dark matter halos and subhalos may contain galaxies. The dark matter halo of a galaxy envelops the galactic disc and extends well beyond the edge of the visible galaxy. Thought to consist of dark matter, halos have not been observed directly. Their existence is inferred through observations of their effects on the motions of stars and gas in galaxies and gravitational lensing. Dark matter halos play a key role in current models of galaxy formation and evolution. Theories that attempt to explain the nature of dark matter halos with varying degrees of success include cold dark matter (CDM), warm dark matter, and massive compact halo objects (MACHOs).

<span class="mw-page-title-main">Reflection mapping</span>

In computer graphics, environment mapping, or reflection 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">South Pole Telescope</span> Telescope at the South Pole

The South Pole Telescope (SPT) is a 10-metre (390 in) diameter telescope located at the Amundsen–Scott South Pole Station, Antarctica. The telescope is designed for observations in the microwave, millimeter-wave, and submillimeter-wave regions of the electromagnetic spectrum, with the particular design goal of measuring the faint, diffuse emission from the cosmic microwave background (CMB). The first major survey with the SPT—designed to find distant, massive, clusters of galaxies through their interaction with the CMB, with the goal of constraining the dark energy equation of state—was completed in October 2011. In early 2012, a new camera (SPTpol) was installed on the SPT with even greater sensitivity and the capability to measure the polarization of incoming light. This camera operated from 2012–2016 and was used to make unprecedentedly deep high-resolution maps of hundreds of square degrees of the Southern sky. In 2017, the third-generation camera SPT-3G was installed on the telescope, providing nearly an order-of-magnitude increase in mapping speed over SPTpol.

<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. 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. The merger rate also provides astronomers with clues about how galaxies bulked up over time.

An inhomogeneous cosmology is a physical cosmological theory which, unlike the currently widely accepted cosmological concordance model, assumes that inhomogeneities in the distribution of matter across the universe affect local gravitational forces enough to skew our view of the Universe. When the universe began, matter was distributed homogeneously, but over billions of years, galaxies, clusters of galaxies, and superclusters have coalesced, and must, according to Einstein's theory of general relativity, warp the space-time around them. While the concordance model acknowledges this fact, it assumes that such inhomogeneities are not sufficient to affect large-scale averages of gravity in our observations. When two separate studies claimed in 1998-1999 that high redshift supernovae were further away than our calculations showed they should be, it was suggested that the expansion of the universe is accelerating, and dark energy, a repulsive energy inherent in space, was proposed to explain the acceleration. Dark energy has since become widely accepted, but it remains unexplained. Accordingly, some scientists continue to work on models that might not require dark energy. Inhomogeneous cosmology falls into this class.

<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.

<span class="mw-page-title-main">Astropy</span> Python language software

Astropy is a collection of software packages written in the Python programming language and designed for use in astronomy. The software is a single, free, core package for astronomical utilities due to the increasingly widespread usage of Python by astronomers, and to foster interoperability between various extant Python astronomy packages. Astropy is included in several large Python distributions; it is part of package managers for Linux and macOS, the Anaconda Python Distribution, Enthought Canopy and Ureka.

<span class="mw-page-title-main">Void (astronomy)</span> Vast empty spaces between filaments with few or no galaxies

Cosmic voids are vast spaces between filaments, which contain very few or no galaxies. Most galaxies are not located in voids, despite their size, due to most galaxies being gravitationally bound together, creating huge cosmic structures known as galaxy filaments. The cosmological evolution of the void regions differs drastically from the evolution of the Universe as a whole: there is a long stage when the curvature term dominates, which prevents the formation of galaxy clusters and massive galaxies. Hence, although even the emptiest regions of voids contain more than ~15% of the average matter density of the Universe, the voids look almost empty to an observer.

<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 . Oxford University Press. 381 (2): 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.