Scipy文件输入/输出mat,wav,mp3

moon (df)

降噪后的moon_cleaned &＃xff08;df&＃xff09;

导入&＃xff1a; import scipy.io as spio

写入 .savemat()

将上面两个df数据&＃xff0c;存储到“moon_cleaned.mat”文件中&＃xff0c;分布命名为{‘moon_cleaned’和’moon_origin’

import scipy.io as spiospio.savemat(&＃39;moon_cleaned.mat&＃39;, mdict&＃61;{&＃39;moon_cleaned&＃39;: moon_cleaned, &＃39;moon_origin&＃39;: moon})


读取 .loadmat()

读取上面的“moon_cleaned.mat”文件

data &＃61; spio.loadmat(&＃39;./moon_cleaned.mat&＃39;)
data


可分别取df,这里取‘moon_origin’为例&＃xff1a;

data[&＃39;moon_origin&＃39;]


使用scipy处理wav文件

• WAV为微软公司&＃xff08;Microsoft)开发的一种声音文件格式&＃xff0c;音质好
• 单声道是一维矩阵&＃xff0c;立体声&＃xff08;双声道&＃xff09;是二维
• 标准的采样频率为&＃xff1a;11.025kHz&＃xff08;语音效果&＃xff09;、22.05kHz&＃xff08;音乐效果&＃xff09;、44.1kHz&＃xff08;高保真效果&＃xff09;
• 导包&＃xff1a;from scipy.io import wavfile

读 wavfile.read(&＃xff09;

from scipy.io import wavfile# 返回值为元组&＃xff0c;分别为&＃xff1a;“采样频率&＃xff1a; 和 ”数据“_, deng &＃61; wavfile.read(&＃39;../data/邓紫棋-喜欢你.wav&＃39;)


wav的音乐文件, 数据类型必须是int16&＃xff0c;取值范围&＃xff1a; -32768 ~ &＃43;32767

查看矩阵形状

deng.shape


10776444 / 44100 剩下的是秒数

#算分钟数
10776444 / 44100 / 60


4.072分钟 相当于 4分42秒

截取一半的音频只需要 len(deng)//2 (整除)

再读取另一首歌

_, lin &＃61; wavfile.read(&＃39;../data/林俊杰-爱不会绝迹.wav&＃39;)
lin


拼接 np.vstack()

将两首歌各自截取一半&＃xff0c;再拼接成一首歌

# 拼接
mix_music &＃61; np.vstack((deng[:len(deng)// 2], lin[:len(lin)//2]))


写 wavfile.write(&＃xff09;

# 导出
wavfile.write(&＃39;mix_music.wav&＃39;, 44100, mix_music)

MP3

使用pydub&＃43; ffmpeg处理mp3文件

• 需要安装pydub
• ffmpeg需要下载,并解压将3个文件放至Python安装目录中的bin目录下,然后还要配置环境变量
• pydup中的AudioSegment可以处理mp3文件

导包&＃xff1a;

from pydub import AudioSegment


读&＃xff1a; AudioSegment.from_mp3()

jing &＃61; AudioSegment.from_mp3(&＃39;../data/难念的经.mp3&＃39;)
type(jing)


返回一个AudioSegment对象

可直接运行播放音频&＃xff08;需要等一阵子&＃xff09;

jing


获取原始数据

jing_data &＃61; jing.get_array_of_samples()


注意&＃xff1a;这里不要直接打印它&＃xff08;因为数据特别长&＃xff0c;可能会卡死&＃xff09;

查看一下类型

type(jing_data)


转换成 ndarray

# 转化为ndarray
jing_ndarray &＃61; np.array(jing_data.tolist())


查看形状

jing_ndarray.shape


查看类型

jing_ndarray.dtype


转化为wav&＃xff1a; wavfile.write&＃xff08;&＃xff09;

wavfile.write(&＃39;经.wav&＃39;, 44100, jing_ndarray.reshape(-1, 2).astype(np.int16))


转换后 4M --> 40&＃43;M


导入另一首歌 并 查看时长

love &＃61; AudioSegment.from_mp3(&＃39;../data/Love Story.mp3&＃39;)
love.duration_seconds


返回的是秒

切片的时候切的是音乐的时长, 单位是毫秒.
先算出两首歌的时长毫秒数

jing.duration_seconds * 1000


love.duration_seconds * 1000


裁剪&＃xff08;就是按毫秒切片&＃xff09;

# 切片的时候切的是音乐的时长, 单位是毫秒.
jing[: 20000]


剪辑合并各一半的两首歌

mix &＃61; jing[: 269609 // 2] &＃43; love[: 234527 // 2]
mix


导出 export&＃xff08;&＃xff09;

mix.export(&＃39;./mix.mp3&＃39;)