问如何使用healpy中的anafast函数画wmap数据的角功率谱？

提问于 2025-10-24 23:11:03

回答 0关注 0查看 6

我打算使用wmap的数据计算角功率谱，直接用下下来的数据对他使用anafast，结果输出的情况并不理想，请问是不是要对数据进行什么操作？

import healpy as hp
import matplotlib.pyplot as plt
import numpy as np

map_data = hp.read_map('/home/mzhang/Desktop/data/observation_data/wmap/wmap_band_forered_iqumap_r9_9yr_Q_v5/q_polarization_with_nobs.fits')

cl = hp.anafast(map_data)

plt.figure(figsize=(10, 6))
l = np.arange(len(cl))
plt.plot(l, cl * l * (l+1) / (2 * np.pi)) 
plt.xlabel('Multipole moment, l')
plt.ylabel('$D_l^{TT}$',fontsize=13)
plt.title('Angular Power Spectrum')
plt.show()

python