python提取图片频谱,如何从频谱图中提取频率?

import cv2

import matplotlib.pyplot as plt

import sys

sys.path.append('/usr/local/lib/python2.7/site-packages')

import glob, os

import numpy as np

import wave from scipy.io

import wavfile

import time

results_dir = "/home/undead/Documents/MetricTestFolder/"

for root, dirs, files in os.walk("/home/undead/ShareWindowsDesktop/SampleWavs"): #LOADS SONGS FROM DIRECTORY

for file in files:

if file.endswith(".wav"):

filepath = os.path.join(root, file)

fs, frames = wavfile.read(filepath)

channels = [

np.array(frames[:, 0]),

np.array(frames[:, 1])

]

fig, (ax2) = plt.subplots(nrows=1)

data, freqs, bins, im = ax2.specgram(channels[0], NFFT=8192, Fs=44100, noverlap=int(8192 * 0.75))

start = time.time()

ax2.pcolormesh(bins, freqs, 10. * np.log10(data))

ax2.axis('tight')

plt.axis('off')

plt.subplots_adjust(left=0.0, right=1.0, bottom=0.0, top=1.0)

plt.xlim([0, 250])

plt.ylim([0, 15000])

plt.savefig(results_dir + "%s.png" % os.path.splitext(file)[0], dpi=450)

plt.close("all")

endt = time.time()

print "%f" % (endt - start)

#end

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值