Python如何可视化包裹路径

Question

我在用DTW算法对齐两份文件的语音时，有一步任务是可视化包裹路径，但是效果图应该是这样的，我试着操作了几遍没有成功，不知道该怎么做？这样的图片要用什么方式获得？

我目前的代码如下：

import librosa
import numpy as np
from fastdtw import fastdtw
import matplotlib.pyplot as plt
#First, upload the file
example_path="./LJ001-0048.wav"
myself_path="./Re_MyVoice.wav"
​
exam_audio,exam_rate=librosa.load(example_path, sr=16000)
my_audio,my_rate=librosa.load(myself_path, sr=16000)
​
#calculate the spec
exam_spect=np.abs(librosa.stft(y=exam_audio,n_fft=exam_rate))
my_spect=np.abs(librosa.stft(y=my_audio,n_fft=my_rate))
​
#use DTW algorithm
distance, path=fastdtw(exam_spect.T, my_spect.T)
​
#Visualize the aligned power spectrogram of the reference speech and my own recording
plt.figure(figsize=(6, 6))
plt.subplot(2, 1, 1)
plt.title('Spectrogram of LJ001-0048.wav')
librosa.display.specshow(librosa.amplitude_to_db(exam_spect), y_axis='log')
plt.colorbar(format='%+2.0f dB')
​
plt.subplot(2, 1, 2)
plt.title('Spectrogram of Re_MyVoice.wav')
librosa.display.specshow(librosa.amplitude_to_db(my_spect), y_axis='log')
plt.colorbar(format='%+2.0f dB')
plt.tight_layout()
plt.show()
​
​
#Visualize the wraped path/map from the DTW algorithm. Here is how the wraped path may look like.
plt.figure(figsize=(6, 6))
b=[]
for i in range(len(path)):
    b.append(np.array( path[i]))
c=np.array(b)
librosa.display.specshow(c, x_axis='time', y_axis='time')
plt.plot(c[:, 1],c[:, 0], color='r')
plt.xlim([0, c.shape[1]])
plt.ylim([0, c.shape[0]])
plt.title('DTW')
plt.colorbar(format='%+2.0f')
plt.tight_layout()
plt.show()

Answer 1

【以下回答由 GPT 生成】

**优化点：**
当前代码已经可以实现所需的功能，但还可以进一步优化：
1. 添加必要的注释说明代码的功能和实现步骤。
2. 使用更具描述性的变量名，减少代码中的硬编码。
3. 将代码封装为可以复用的函数，以提高代码的可读性和维护性。
4. 考虑处理异常情况，例如文件不存在或加载失败。
5. 提供更详细的文档说明，包括函数的输入和输出，以及对所用库和算法的简要介绍。

【相关推荐】

• 这篇博客: 初识DTW(动态时间规整)算法及Python实现例中的 3. DTW基本算法 部分也许能够解决你的问题。

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 写成博客, 将相关链接放在评论区, 以帮助更多的人 ^-^