1
我试图在healpix地图上使用healpy来产生一个波束。对于初学者,我希望能够在mollweide投影中产生2D高斯,但我真的不知道从哪里开始。在healpy中绘制一个numpy数组
我可以定义一个二维高斯:
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))
,这样我可以建立类似的三维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像素化,所以a * map *是一维数组,其中索引对应于像素。如果你只需要一个Mollweide投影,你可以使用'matplotlib',参见http://matplotlib.org/examples/pylab_examples/geo_demo.html – 2014-08-28 22:28:56