基于CWT系数的PCG信号尺度图生成问题排查求助

我尝试生成基于CWT系数的心音图（PCG）尺度图，但无法复现目标论文中的结果。已复现论文中相同记录的所有预处理步骤，处理后的PCG信号与论文中的图一致：

CWT计算代码 #

使用pywavelets库的Morlet小波，代码如下：

sampling_rate = fs = 2000
wavelet = pywt.ContinuousWavelet("morl")
scales = pywt.central_frequency(wavelet) * sampling_rate / np.arange(10, 150, 1)
pcg = sig[:10000]

pcg = preprocess_pcg(pcg) # 包含所有预处理步骤，包括带通滤波
data, f =  pywt.cwt(pcg, scales, wavelet, sampling_period= 1/fs)

#也曾尝试scales = np.arange(7,160,1)，但未解决问题

尺度图绘制代码 #

生成并导出尺度图的代码：

data=np.abs(data)
cmap = plt.get_cmap('jet', 256)
fig = plt.figure(figsize=(5,5))
ax = fig.add_subplot(111)
t = np.arange(data.shape[1]) / fs  

ax.pcolormesh(t, f, data, cmap=cmap, vmin=data.min(), vmax=data.max(), shading='auto')

plt.show()

当前生成结果 #

生成的尺度图：

对数Y轴版本：

目标论文结果 #

上述结果与论文中的图像差异显著：

请问我哪里操作有误？

内容的提问来源于stack exchange，提问作者Throwaway123