No project description provided
Project description
hpinterp
A faster (~100x) and more accurate interpolation algorithm for HEALPix maps by Prof. Jonathan Sievers. The code pre-interpolates the map data in advance, saving time at the evaluation stage.
Example:
import numpy as np
import healpy as hp
from hpinterp import InterpMap
map_ = hp.read_map("your_map.fits")
# Set nest=True if your map has nested ordering
interp_map = InterpMap(map_)
# Generating example coordinate grid of co-latitude and longitude in radians
npoints = int(1e5)
theta = np.random.rand(npoints) * np.pi
phi = np.random.rand(2 * npoints) * 2 * np.pi
# Get interpolated values. Set lonlat=True if using longitude and latitude in degrees
interp_result = interp_map(theta, phi)
# or
interp_result = interp_map.get_interp_val(theta, phi)
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
hpinterp-0.1.0.tar.gz
(3.8 kB
view hashes)