[ ] convenient constructors, factories an/or helpers (e.g., construct a grid given a approx. fixed cell resolution)
[x] Python bindings
[ ] ghost node status: skip in flow routers, sink resolvers and eroders
[ ] node status utility methods (e.g., set continental mask from arbitrary polygons).
[ ] interpolation utility methods (some requires little work while others may require much more work, not all may be needed)
[x] return geodetic or geocentric coordinates of a given cell index, a vector of cell indices or all cells -> useful for mapping field values from an external mesh onto the healpix grid using 3rd-party interpolation
[ ] return an interpolated value at an arbitrary point (and/or vectorized version), using Healpix's (local) linear interpolation methods -> useful for mapping field values from the Healpix grid onto another mesh of similar resolution
[ ] return an aggregated field value (mean) for a given bounding box or arbitrary polygon, using Healpix's spatial query methods -> useful for mapping fields values from the Healpix grid onto a coarse grid or mesh.
[ ] return interpolated values for a given set of arbitrary points using global, point-wise interpolation (spherical harmonics, polynomial fitting) -> useful for mapping fields between different Healpix grids of different resolutions with fine control on the surface smoothing. Might require a lot of work (no high-level method for this available in Healpix, only command-line programs built on top basic functionality, e.g., see smoothing_cxx)
[x] configuration
[x] FindHealpix CMake module
[x] Cmake option to enable / disable Healpix support (opt-in)
smoothing_cxx)ducc0instead of healpix original library (the latter is not supported on Windows), see https://github.com/healpy/healpy/issues/717.