我想从拟合文件中获取极化并使用Healpy绘制它。我看到Healpix有一个名为POLARIZATION的内置选项,但我没有看到Healpy的相同选项。
我如何绘制这些数据?
谢谢!
编辑: @AndreaZonca 以下是我目前用于绘制温度的代码:
t = "%s Field - WMAP at W Band (90 GHz) Foreground Reduced iquma"%(name)
map = hp.read_map(band_iquma)
hp.cartview(map,title = t, coord = ['G','C'],max = .25, min = -.25, unit = 'mk',\
lonra = range[0], latra = range[1])
hp.graticule()
plt.show()
我收到此文件的网站说,极化数据包含在拟合文件中。
这个网站说如果给出关键字POLARIZATION和数字1,2或3,HealPix可以绘制极化。
答案 0 :(得分:0)
您要查找的关键字field
位于read_map
,如果您将其设置为(0,1,2),则会同时读取温度和极化。
然后你可以用[1]绘制Q切片地图,用[2]绘制U切片。
参见:
http://healpy.readthedocs.org/en/latest/generated/healpy.fitsfunc.read_map.html