sphere_healpy.md

Plot a sphere and pixel centers on it using Healpix package

tags: python, healpix, healpy

Snippet

from mayavi import mlab
import numpy as np
import healpy as hp

# Create a sphere
r = 1.0
pi = np.pi
cos = np.cos
sin = np.sin
phi, theta = np.mgrid[0:pi:101j, 0:2 * pi:101j]

x = r*sin(phi)*cos(theta)
y = r*sin(phi)*sin(theta)
z = r*cos(phi)

mlab.figure(1, bgcolor=(1, 1, 1), fgcolor=(0, 0, 0), size=(400, 300))
mlab.clf()

pix_centers = hp.pix2vec(1, np.arange(12))
xx, yy, zz = pix_centers


mlab.mesh(x , y , z, color=(0.0,0.5,0.5))
mlab.points3d(xx, yy, zz, scale_factor=0.2)


mlab.show()