记录帖子,首先要确定每个通道对应的频率,然后他需要一个频率对应不同分辨率的文件,然后线性拟合得到对应通道的分辨率。
from scipy import interpolate
freq0 = 1000.00357627869
df = 0.00762939453125
chan=5000
freq = (freq0 + chan * df)/1.e3 # in GHz
beam_ave_file = 'beam_ave.dat'
tmp1 = np.fromfile(beam_ave_file, dtype=float,count=-1,sep='\n')
tmp1 = tmp1.reshape((tmp1.size//2, 2))
f1 = interpolate.interp1d(tmp1[:,0], tmp1[:,1], kind='linear')
hpbw_o = f1(freq)