Я пытаюсь создать луч на карте healpix, используя healpy. Для начала я хотел бы иметь возможность создавать 2D-гауссово в проекции mollweide, но я действительно не знаю, с чего начать.Построение массива numpy в healpy
я могу определить 2D гауссово:
import numpy as np
def gaussian_2D(x,y,mu_x=0.,mu_y=0.,sig_x=1.,sig_y=1.):
return np.exp(-0.5*(((x-mu_x)/sig_x)**2 + ((y-mu_y)/sig_y)**2))
таким образом, что я могу построить 3D X, Y, Z пространство как:
delta = 0.025
x = np.arange(-4, 4, delta)
y = np.arange(-4, 4, delta)
X, Y = np.meshgrid(x,y)
Z = gaussian_2D(X,Y)
, но здесь я довольно потерял, и не может найти много полезной документации, касающейся того, как и как проектировать. Любые предложения относительно направления атаки будут высоко оценены!
'healpy' использует пикселизацию HEALPix, поэтому * map * представляет собой 1D-массив, где индексы соответствуют пикселям. Если вам просто нужна проекция Mollweide, вы можете просто использовать «matplotlib», см. Http://matplotlib.org/examples/pylab_examples/geo_demo.html –