音频去除前后端的静音

Audio Silence Trimmer Pro 这个软件是可以的，网上有破解版，但在我电脑上不好用
去除前后端的静音部分，顺便再看看波形图对比这样

import librosa
import wave
import pylab as pl
import numpy as np
import os
from pydub import AudioSegment
from pydub.silence import detect_silence

from tqdm import tqdm
l = os.listdir('F:\我的程序练习\BZNSYP\sox')
for i in tqdm(l):
    f = wave.open('F:\我的程序练习\BZNSYP\sox\\'+i, 'rb')
    sound = AudioSegment.from_file('F:\我的程序练习\BZNSYP\sox\\'+i, format="wav")
    start_end = detect_silence(sound, 250, -40, 1)
    soundstat = 0
    soundend = len(sound)
    for o in start_end:
        if o[1] == len(sound):
            soundend = o[0]
        if o[0] == 0:
            soundstat = o[1]
    final_name = 'F:\我的程序练习\BZNSYP\data\\'+i
    sound[soundstat:soundend].export(final_name, format='wav')
    f1 = wave.open('F:\我的程序练习\BZNSYP\Wave\\' + i, 'rb')


    # f = wave.open('E:\\tacotron2\LJ050-0258.wav','rb')
    params = f.getparams()
    nchannels, sampwidth, framerate, nframes = params[:4]
    # 读取波形数据
    str_data = f.readframes(nframes)
    f.close()
    # 将波形数据转换为数组
    wave_data = np.fromstring(str_data, dtype=np.short)

    # f = wave.open('E:\\tacotron2\LJ050-0258.wav','rb')
    params1 = f1.getparams()
    nchannels1, sampwidth1, framerate1, nframes1 = params1[:4]
    # 读取波形数据
    str_data1 = f1.readframes(nframes1)
    f1.close()
    # 将波形数据转换为数组
    wave_data1 = np.fromstring(str_data1, dtype=np.short)

    if nchannels == 2:
        wave_data.shape = -1, 2
        wave_data = wave_data.T
        time = np.arange(0, nframes) * (1.0 / framerate)
        # 绘制波形
        pl.subplot(211)
        pl.plot(time, wave_data[0])
        pl.subplot(212)
        pl.plot(time, wave_data[1], c="g")
        pl.xlabel("time (seconds)")
        pl.show()
    elif nchannels == 1:
        wave_data.shape = -1, 1
        wave_data = wave_data.T
        time = np.arange(0, nframes) * (1.0 / framerate)
        # 绘制波形

        wave_data1.shape = -1, 1
        wave_data1 = wave_data1.T
        time1 = np.arange(0, nframes1) * (1.0 / framerate1)

        pl.subplot(211)
        pl.plot(time, wave_data[0])
        pl.subplot(212)
        pl.plot(time1, wave_data1[0])
        pl.xlabel("time (seconds)")
        pl.show()