pip install moviepy -i http://pypi.douban.com/simple --trusted-host pypi.douban.com

import moviepy.editor as mpy


# 截取背景音樂(lè)
audio_background = mpy.AudioFileClip(r'D:\python\pycharm2020\my_program\video_process\test_01.mp4').subclip(1, 60)
audio_background.write_audiofile(r'D:\python\pycharm2020\my_program\video_process\vst01.mp3')

import numpy as np  # 科學(xué)計(jì)算所用的numpy庫(kù)
import matplotlib.pyplot as plt  # 繪圖所用的庫(kù)matplotlib
from scipy.io import wavfile  # 讀取波形文件所用的庫(kù)


rate_h, hstrain = wavfile.read(r"H1_Strain.wav", "rb")   # 讀取音頻文件
rate_l, lstrain = wavfile.read(r"L1_Strain.wav", "rb")
# reftime, ref_H1 = np.genfromtxt('GW150914_4_NR_waveform_template.txt').transpose()
reftime, ref_H1 = np.genfromtxt('wf_template.txt').transpose()  # 使用python123.io下載txt文件

htime_interval = 1 / rate_h
ltime_interval = 1 / rate_l
fig = plt.figure(figsize=(12, 6))  # 創(chuàng)建大小為12*6的繪圖空間

# 丟失信號(hào)起始點(diǎn)
htime_len = hstrain.shape[0] / rate_h  # 讀取數(shù)據(jù)第一維的長(zhǎng)度，得到函數(shù)在坐標(biāo)軸上總長(zhǎng)度
htime = np.arange(-htime_len / 2, htime_len / 2, htime_interval)  # （起點(diǎn)，終點(diǎn)，時(shí)間間隔）

plth = fig.add_subplot(221)  # 設(shè)置繪圖區(qū)域
plth.plot(htime, hstrain, 'r')  # 畫(huà)出以時(shí)間為x軸，應(yīng)變數(shù)據(jù)為y軸的圖像，‘y'為黃色
plth.set_xlabel('Time (seconds)')
plth.set_ylabel('H1 Strain')
plth.set_title('H1 Strain')

ltime_len = lstrain.shape[0] / rate_l
ltime = np.arange(-ltime_len / 2, ltime_len / 2, ltime_interval)
pltl = fig.add_subplot(222)
pltl.plot(ltime, lstrain, 'k')
pltl.set_xlabel('Time (seconds)')
pltl.set_ylabel('L1 Strain')
pltl.set_title('L1 Strain')

pltref = fig.add_subplot(212)
pltref.plot(reftime, ref_H1, 'purple')
pltref.set_xlabel('Time (seconds)')
pltref.set_ylabel('Template Strain')
pltref.set_title('Template')
fig.tight_layout()  # 自動(dòng)調(diào)整圖像外部邊緣

plt.savefig("Gravitational_Waves_Original.png")  # 保存圖像為png格式
plt.show()
plt.close(fig)