我想从一个fits文件中访问极化,并使用Healpy绘制它。我看到Healpix有一个内置的选项叫做极化,但我不认为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文件中。
这个网站说,如果给定关键字极化和数字1、2或3,HealPix可以绘制极化图。
http://healpix.jpl.nasa.gov/html/idlnode10.htm
发布于 2014-03-10 15:13:17
您要查找的关键字是field
in read_map
,如果将其设置为(0,1,2),您将同时读取温度和极化。然后,您可以用1绘制Q,用2绘制U。参见:
map.html
https://stackoverflow.com/questions/22239720
复制相似问题