使用Healpy从fits文件绘制极化图

Question

我想从一个fits文件中提取极化数据，并使用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()

我从这个网站上获取的文件说极化数据包含在fits文件中。

而这个网站则表示，如果给HealPix提供关键词POLARIZATION和数字1、2或3，它可以绘制极化数据。

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

Answer 1

你要找的关键词是 field，在 read_map 里。如果你把它设置为 (0,1,2)，那么你就可以读取温度和极化数据。接着，你可以用 [1] 来绘制 Q 图，用 [2] 来绘制 U 图。

详细信息可以查看这里：

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