Data compression on the sphere

Large data-sets defined on the sphere arise in many fields. In particular, recent and forthcoming observations of the anisotropies of the cosmic microwave background (CMB) made on the celestial sphere contain approximately three and fifty mega-pixels respectively. The compression of such data is therefore becoming increasingly important. We develop algorithms to compress data defined on the sphere. A Haar wavelet transform on the sphere is used as an energy compression stage to reduce the entropy of the data, followed by Huffman and run-length encoding stages. Lossless and lossy compression algorithms are developed. We evaluate compression performance on simulated CMB data, Earth topography data and environmental illumination maps used in computer graphics. The CMB data can be compressed to approximately 40% of its original size for essentially no loss to the cosmological information content of the data, and to approximately 20% if a small cosmological information loss is tolerated. For the topographic and illumination data compression ratios of approximately 40:1 can be achieved when a small degradation in quality is allowed. We make our SZIP program that implements these compression algorithms available publicly.

[1]  W. Sweldens.  The Lifting Scheme: A Custom - Design Construction of Biorthogonal Wavelets "Industrial Mathematics , 1996 .

[2]  Syed Twareque Ali,et al.  Two-Dimensional Wavelets and their Relatives , 2004 .

[3]  D. Healy,et al.  Computing Fourier Transforms and Convolutions on the 2-Sphere , 1994 .

[4]  Bruno Torrésani,et al.  Position-frequency analyis for signals defined on spheres , 1995, Signal Process..

[5]  Willi Freeden,et al.  Combined Spherical Harmonic and Wavelet Expansion—A Future Concept in Earth's Gravitational Determination , 1997 .

[6]  Laurent Jacques,et al.  Wavelets on the sphere: implementation and approximations , 2002 .

[7]  Graham J. L. Kemp,et al.  Fast computation, rotation, and comparison of low resolution spherical harmonic molecular surfaces , 1999, J. Comput. Chem..

[8]  K. Gorski,et al.  HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere , 2004, astro-ph/0409513.

[9]  R. J. Willemann,et al.  Role of membrane stresses in the support of planetary topography , 1981 .

[10]  Willi Freeden,et al.  Constructive Approximation on the Sphere: With Applications to Geomathematics , 1998 .

[11]  Mark A. Wieczorek,et al.  Spatiospectral Concentration on a Sphere , 2004, SIAM Rev..

[12]  Laurent Jacques,et al.  Stereographic wavelet frames on the sphere , 2005 .

[13]  Jean-Luc Starck,et al.  Wavelets, ridgelets and curvelets on the sphere , 2006 .

[14]  Edward J. Wollack,et al.  First-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results , 2003, astro-ph/0302207.

[15]  Matthias Holschneider,et al.  Continuous wavelet transforms on the sphere , 1996 .

[16]  S. Swenson,et al.  Methods for inferring regional surface‐mass anomalies from Gravity Recovery and Climate Experiment (GRACE) measurements of time‐variable gravity , 2002 .

[17]  Roger J. Phillips,et al.  Potential anomalies on a sphere: Applications to the thickness of the lunar crust , 1998 .

[18]  Pat Hanrahan,et al.  A signal-processing framework for reflection , 2004, ACM Trans. Graph..

[19]  George G. Lorentz,et al.  Constructive Approximation , 1993, Grundlehren der mathematischen Wissenschaften.

[20]  Mark S. Gordon,et al.  Rapid and stable determination of rotation matrices between spherical harmonics by direct recursion , 1999 .

[21]  P. Hanrahan,et al.  Triple product wavelet integrals for all-frequency relighting , 2004, SIGGRAPH 2004.

[22]  O. V. Verkhodanov,et al.  Gauss - Legendre sky pixelization (GLESP) for CMB maps , 2005 .

[23]  Kathy Whaler,et al.  Downward continuation of Magsat lithospheric anomalies to the Earth's surface , 1994 .

[24]  Belgium,et al.  Correspondence principle between spherical and euclidean wavelets , 2005, astro-ph/0502486.

[25]  Wim Sweldens,et al.  The lifting scheme: a construction of second generation wavelets , 1998 .

[26]  David Assaf,et al.  Data Compression on the Sphere Using Faber Decomposition , 2022 .

[27]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[28]  P. Vandergheynst,et al.  Wavelets on the 2-sphere: A group-theoretical approach , 1999 .

[29]  O. Blanc,et al.  Exact reconstruction with directional wavelets on the sphere , 2007, 0712.3519.

[30]  James A. Storer,et al.  DATA COMPRESSION CONFERENCE , 2001 .

[31]  Pierre Vandergheynst,et al.  Wavelets on the n-sphere and related manifolds , 1998 .

[32]  M. Halpern,et al.  Three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: temperature analysis , 2006 .

[33]  C. H. Lineweaver,et al.  Applications of Wavelets to the Analysis of Cosmic Microwave Background Maps , 1999, astro-ph/9903206.

[34]  S. D. M. White,et al.  Understanding the halo-mass and galaxy-mass cross-correlation functions , 2007, 0709.3933.

[35]  Laurent Demanet,et al.  Gabor wavelets on the sphere , 2003, SPIE Optics + Photonics.

[36]  Anthony N. Lasenby,et al.  Testing the Gaussianity of the COBE DMR data with spherical wavelets , 2000 .