python – 应用傅里叶变换后scipy.io波形文件处理的问题
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了python – 应用傅里叶变换后scipy.io波形文件处理的问题,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含1702字,纯文字阅读大概需要3分钟。
内容图文
![python – 应用傅里叶变换后scipy.io波形文件处理的问题](/upload/InfoBanner/zyjiaocheng/697/0d68be8daefa49c9b0a1a95d5c28ecf5.jpg)
我正在使用scipy.fftpack和scipy.io包在Python中使用一些音频文件.这意味着我有我正在导入的波形文件,使用傅里叶变换玩它们,然后输出到新的波形文件.然而,我遇到了一些问题,在运行这些转换之后,波形文件将无法播放,大约是原始文件大小的4倍.
目前我只是导入歌曲,获取速率,导入数据,执行ifft(fft(数据))然后输出.我已经尝试获取这些浮点值并将其转换为整数并计算初始导入数据和ifft(fft(数据))之间的差异,并看到它正好为零.
以下是我的问题:
>有人知道为什么我在试图进行傅立叶变换后试图播放的波形文件不能播放吗?
>对波浪文件是否有任何限制可能导致一般问题.或者数据必须是整数?
现在已经很晚了,我觉得我的帖子可能会分散,所以如果您需要更多信息,请询问,我会尽力给出更好的描述.
解决方法:
您需要在处理后将数据转换为适当位深度的整数类型.以this file为例:
>>> import scipy.io.wavfile
>>> rate, data = scipy.io.wavfile.read('Happy Tree Friends.wav')
>>> rate
8000
>>> data
array([ 5, -5, 5, ..., 0, -1, 0], dtype=int16)
>>> data_bis = np.fft.ifft(np.fft.fft(data))
>>> data_bis
array([ 5.00000000e+00 -1.55406753e-11j,
-5.00000000e+00 +1.95349676e-11j,
5.00000000e+00 +1.41131140e-11j, ...,
8.06674092e-12 -7.58643463e-13j,
-1.00000000e+00 -2.21611283e-12j, -2.04999489e-11 +4.55890751e-12j])
>>> data_bis.dtype
dtype('complex128')
尽管数据中的值非常接近data_bis中的值,但它们是非常不同的动物,如下所示:
>>> scipy.io.wavfile.write('test.wav', rate, data_bis)
>>> scipy.io.wavfile.read('test.wav')
TypeError: data type not understood
但是如果你只是将处理后的结果转换回原来的dtype,那么一切都很好用:
>>> scipy.io.wavfile.write('test.wav', rate, data_bis.astype(data.dtype))
__main__:1: ComplexWarning: Casting complex values to real discards the imaginary part
>>> scipy.io.wavfile.read('test.wav')
(8000, array([ 4, -5, 4, ..., 0, -1, 0], dtype=int16))
内容总结
以上是互联网集市为您收集整理的python – 应用傅里叶变换后scipy.io波形文件处理的问题全部内容,希望文章能够帮你解决python – 应用傅里叶变换后scipy.io波形文件处理的问题所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。