1

我有一个频谱图：

在此处输入图像描述

我想清理频谱图，所以我只捕获特定范围内的频率（即在此示例中，在 2627 - 3939 之间）并删除低于此频率的所有块。我的总体目标是只留下这个频率范围内的 4 个片段，并且可以识别。

到目前为止，这是我的代码：

import wave, struct, numpy as np, matplotlib.mlab as mlab, pylab as pl
def wavToArr(wavefile):
    w = wave.open(wavefile,"rb")
    p = w.getparams()
    s = w.readframes(p[3])
    w.close()
    sd = np.fromstring(s, np.int16)
    return sd,p

def wavToSpec(wavefile,log=False,norm=False):
    wavArr,wavParams = wavToArr(wavefile)
    print wavParams
    return  mlab.specgram(wavArr, NFFT=256,Fs=wavParams[2],window=mlab.window_hanning,noverlap=128,sides='onesided',scale_by_freq=True)

wavArr,wavParams = wavToArr("4bats.wav")
Pxx, freqs, bins = wavToSpec("4bats.wav")
Pxx += 0.0001

freqs += (len(wavArr) / wavParams[2]) / 2.
hf=pl.figure(figsize=(12,12));
ax = hf.add_subplot(2,1,1);
#plot spectrogram as decibals
hm = ax.imshow(10*np.log10(Pxx),interpolation='nearest',origin='lower',aspect='auto')
hf.colorbar(hm)
ylcnt = len(ax.get_yticklabels())
ycnt = len(freqs)
ylstep = int(ycnt / ylcnt)
ax.set_yticklabels([ int(freqs[f]) for f in xrange(0,ycnt,ylstep) ])
pl.show()

问题是，我不知道如何使用 Python 做到这一点。我知道范围（2627 - 3939）但是，我会遍历整个二维数组并总结所有块，或者，对于频谱图中的每个块，计算频率，如果它高于阈值，保持它，否则值变为0.0？

如果我总结每个垃圾箱，我会得到以下信息：

在此处输入图像描述

我需要保留这些块，但是，除了这些块之外，我想删除所有其他块。

我希望有一个人可以帮助我！

4

1 回答 1

1

也许你想要类似的东西：

Pxx[np.greater(np.sum(Pxx, axis=1), 2955), :] = 0.

或者，我可能已经切换了你的轴，所以也试试：

Pxx[:, np.greater(np.sum(Pxx, axis=0), 2955)] = 0.

或者也许你想要别的东西......我觉得这个问题有点不清楚。

于 2014-02-03T05:28:38.907 回答