cdshealpix-rust

CDS implementation of the HEALPix tesselation in Rust and modules to generate libraries in WebAssembly, Python, ...

About

This library is an implementation in Rust of the HEALPix tesselation. This implementation has been made by the Strasbourg astronomical Data Centre (Centre de Données astronomique de Strasbourg, CDS).

Initially, it is a port of a part of the CDS Java library available here, but improvement have been added while porting the code.

For informations on HEALPix in general, see: * The official web site * The Wikipedia page * The two main reference papers: Gorsky (2005) and Calabretta (2007)

Official implementations, are available here. It contains GPL v2 codes in Fortran, C++, Java, IDL, Python, ...

Other independant HEALPix implementations: * Astropy-healpix python wrapper using a C code (C code by Dustin Lang) * Javascript/Typescript implementation by Koike Michitaro * ... (Help me to add links to other HEALPix resources and codes).

Features

• Supports the HEALix Nested scheme
• Supports approximated cone and elliptical cone and exact polygon queries
• Supports BMOC (MOC with a flag telling if a cell is fully or partially covered by a surface) as a result of cone and polygon queries
• Supports logical operations on BMOCs and BMOC creation from a list of cell number at a given depth

Missing Features

• Not supported
  • RING scheme
  • Spherical Harmonics computations
  • (Help me fill this)
• Not yet implemented
  • Exact cone and ellipse solution (but using the custom approx methods, one can handle the rate of false positives)

Examples

Compute the cell number of a given position on the unit-sphere at a given HEALPix depth.

```rust use cdshealpix::{nside}; use cdshealpix::nested::{getorcreate, Layer};

let depth = 12u8; let lon = 12.5f64.toradians(); let lat = 89.99999f64.to_radians();

let nestedd12 = getorcreate(depth); let nside = nside(depth) as u64; let expectedcellnumber = nside * nside - 1

asserteq!(expectedcellnumber, nestedd12.hash(lon, lat)); ```

Get the spherical coorinates of the 4 vertices of a given cell at a given depth:

```rust use cdshealpix::nested::{getorcreate, Layer};

let depth = 12u8; let cellnumber= 10_u64;

let nestedd12 = getor_create(depth);

let [ (lonsouth, latsouth), (loneast, lateast), (lonnorth, latnorth), (lonwest, latwest) ] = nestedd12.vertices(cellnumber);

```

Get a hierarchical view (a MOC) on the cells overlapped by a given cone:

```rust use cdshealpix::nested::{getorcreate, Layer};

let depth = 6u8; let nestedd6 = getorcreate(depth);

let lon = 13.158329f64.toradians(); let lat = -72.80028f64.toradians(); let radius = 5.64323f64.toradians();

let moc = nestedd6.coneoverlapapprox(lon, lat, radius); for cell in moc.intoiter() { println!("cell: {:?}", cell); } ```

Standalone

(Not on crates.io, but on github) The code source of the very beginning of a standalone exec can be found in cli/src/bin.rs.

WebAssembly

(Not on crates.io, but on github) To build and use the WebAssembly (and Javascript) files, the libwasmbingen directory. We rely on wasm-bingen.

Python

(Not on crates.io, but on github) See the libpython directory containing a very first integration in python using CFFI.

For a clean Python wrapper and associated Wheels, see Matthieu Baumann's project cds-healpix-python. To use the library in python, install it through pip (examples are provided on github cds-healpix-python): bash pip install cdshealpix

License

Like most projects in Rust, this project is licensed under either of

• Apache License, Version 2.0, (LICENSE-APACHE or http://www.apache.org/licenses/LICENSE-2.0)
• MIT license (LICENSE-MIT or http://opensource.org/licenses/MIT)

at your option.

Contribution

Unless you explicitly state otherwise, any contribution intentionally submitted for inclusion in this project by you, as defined in the Apache-2.0 license, shall be dual licensed as above, without any additional terms or conditions.

Disclaimer

It a first code in Rust, feel free to give some advice/feedback.