使用Healpy从fit文件绘制极化

时间:2014-03-07 01:30:28

标签: python matplotlib healpy

我想从拟合文件中获取极化并使用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可以绘制极化。

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

1 个答案:

答案 0 :(得分:0)

您要查找的关键字field位于read_map,如果您将其设置为(0,1,2),则会同时读取温度和极化。 然后你可以用[1]绘制Q切片地图,用[2]绘制U切片。 参见:

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