如何解决将等距频率轴更改为等距波长轴-python 3.8
是否有一种简便的方法将 等频率轴更改为等距波长轴 当频谱以spec = {ndarray:(101,Fiveormoredigits)}的形式给出时?
光谱范围
#spectral extend in THz
v_min = (fd.frequency_0 - fd.freq_intervall / 2.0) * 1e-12
v_max = (fd.frequency_0 + fd.freq_intervall / 2.0) * 1e-12
绘制光谱
fig2 = plt.figure()
bx2 = fig2.add_subplot(122)
bx2.imshow(spec,aspect='auto',interpolation='bicubic',cmap='jet',vmin=spec_max-db_range,vmax=spec_max,origin='lower',extent=[l_min,l_max,fd.distance])
bx2.set_xlabel('Frequenz (THz)',fontsize=20)
bx2.set_ylabel('Propagation (m)',fontsize=20)
fig2.suptitle('PM 980 XP/HP: 5cm,200nJ',fontsize=25)
matplotlib.pyplot.xticks(fontsize=15)
matplotlib.pyplot.yticks(fontsize=15)
plt.show()`
#spectral extend in THz
l_min = 3e8/(fd.frequency_0 - fd.freq_intervall / 2.0) * 1e9
l_max = 3e8/(fd.frequency_0 + fd.freq_intervall / 2.0) * 1e9
问题是,轴是等距的,但中心波长不正确。 使用这种解决方案,每个绘图轴都需要移动。
fig2 = plt.figure()
bx2 = fig2.add_subplot(122)
bx2.imshow(spec,fontsize=25)
x_label_list = ['800','900','1035','1200','1300','1400']
bx2.set_xticks([375,333,289,250,230,214])
bx2.set_xticklabels(x_label_list)
matplotlib.pyplot.xticks(fontsize=15)
matplotlib.pyplot.yticks(fontsize=15)
plt.show()`
问题是,一方面中心波长正确,但轴不等距。
其他。我试图进行插值:
变量:spec = {ndarray:(101,32768)}
X = np.linspace(0,100,101)
Y = np.linspace(3e8/v_max,3e8/v_min,np.size(spec,1))
x,y = np.meshgrid(X,Y)
spec_new = interpolate.interp2d(x,y,spec,kind='linear')
因此,基本上,规格由101行组成,分别代表框架和32768列。每列包含一个能量值。我试图在新的等距轴上插值能量,但发生错误:
raise OverflowError(msg)
OverflowError: Too many data points to interpolate
我尝试重新缩放能量并在新的等距轴上进行插值,但发生错误:
spec_new = interpolate.interp2d(x,spec/(X*X),kind='linear')
ValueError: operands could not be broadcast together with shapes (101,8192) (8192,101)
解决方法
我尝试逐步插值,但这不起作用:
x = np.linspace(3e8/v_max,3e8/v_min,np.size(spec,1))
scale = x*x
s = 0
n = 100
for i in range(1,n):
spec[i,:] = interpolate.interp1d(x,spec[i,:]*scale,'linear')
print(spec);
发生错误。为什么?
spec[i,:],'linear')
TypeError: float() argument must be a string or a number,not 'interp1d'
为什么这不可能呢?:
for i in range(1,n):
spec_neu[i,'linear')
print(spec_neu);
spec_neu[i,'linear')
NameError: name 'spec_neu' is not defined
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。