用Healpy绘制fits文件的极化曲线

Question

我想从一个fits文件访问极化，并用Healpy绘制它。我看到Healpix有一个内置的选项叫做POLARIZATION，但是Healpy没有相同的选项。在

如何绘制这些数据？在

谢谢！在

编辑： @安德烈亚森卡 以下是我当前用于绘制温度图的代码：

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()

我得到这个文件的网站说极化数据包含在fits文件中。在

这个网站说HealPix可以绘制极化，如果给定关键字polarization和数字1、2或3。在

http://healpix.jpl.nasa.gov/html/idlnode10.htm

Answer 1

你要找的关键字是field中的read_map，如果你把它设为（0,1,2），你将同时读取温度和极化。 然后你可以用[1]和[2]来绘制地图的Q切片。 参见：

http://healpy.readthedocs.org/en/latest/generated/healpy.fitsfunc.read_map.html