我一直在尝试用CAMB和healpy来制作E和B模式的映射,但是到目前为止,我所做的并不太像。
我使用CAMB在以下链接上生成模拟映射:form.cfm,除了选中张量cl和变量r的框外,我保留所有默认参数。
因此,我一直在使用以下代码来生成地图:
# read in the power spectra produced by camb, ordering is: [TT, EE, BB, TE]
cl = hp.read_cl('scalar_and_tensor_power_spectra_produced_by_CAMB.fits')
# set the random seed for phase information
numpy.random.seed(some number) #Im using 1 for now
#Produce E mode and B mode maps
EEmap = hp.synfast(cl[1], NSIDE)
BBmap = hp.synfast(cl[2], NSIDE)
一个快速的谷歌搜索显示了与我得到的非常不同的结果,你可以在下面的链接中看到:http://tinypic.com/view.php?pic=ifvrpt&s=9#.VmY7Jx8So8o http://tinypic.com/view.php?pic=a5ngx&s=9#.VmY7KR8So8o
我希望在E和B中看到某些模式,因为它们是渐变和卷曲成分,这在我在网上找到的图像中很明显,但在我的产品中却没有。
所以我有三个问题:
- something like at the following link: [http://cosmologist.info/polar/EBsupport.html](http://cosmologist.info/polar/EBsupport.html)
发布于 2015-12-08 07:06:29
synalm
http://healpy.readthedocs.org/en/1.x/generated/healpy.sphtfunc.synalm.html#healpy.sphtfunc.synalm的帮助,如果您使用hp.synfast(cl)
来创建温度图,则使用温度中的cl
。您应该为每个组件创建一个Cl列表,即TT、EE、BB、TE,然后所有的Cl都为零,而只有一个非零。healpy
中绘制向量https://stackoverflow.com/questions/34146794
复制相似问题