Détecter robinet avec pyaudio de vivre mic

Comment puis-je utiliser pyaudio pour détecter un soudain bruit de tapotement à partir d'un micro?

InformationsquelleAutor a sandwhich | 2010-11-11
  1. 95

    Un sens, je l'ai fait:

    • lire un bloc d'échantillons à la fois,
      dire 0.05 secondes
    • calculer la
      Amplitude RMS du bloc (carré
      la racine de la moyenne des carrés des
      les échantillons individuels)
    • si le bloc de l'amplitude RMS est supérieure à un seuil, c'est un "bruit de bloc" d'autre, c'est un "bloc de calme"
    • un coup du robinet serait un bloc de calme suivie par un petit nombre de noisy blocs suivi par un bloc de calme
    • si vous n'avez jamais obtenir un bloc de calme, votre seuil est trop bas
    • si vous n'avez jamais obtenir un bruit de bloc, votre seuil est trop élevé

    Ma demande a été d'enregistrement "intéressant" bruits sans surveillance, il serait dossier tant qu'il n'y étaient bruyants blocs. Il serait de multiplier le seuil de 1,1 si il y a 15 secondes bruyant période ("couvrant ses oreilles") et de multiplier le seuil de 0,9 si il y avait 15minutes période de calme ("écoute plus difficile"). Votre application aura des besoins différents.

    J'ai aussi remarqué que certains commentaires dans mon code concernant observé valeurs RMS. Sur le construit en mic sur un Macbook Pro, avec +/- 1.0 audio normalisé de la plage de données, en collaboration avec le volume réglé sur max, certains points de données:

    • De 0,003-0.006 (-50dB à -44dB) un odieusement fort chauffage central ventilateur dans ma maison
    • 0.010-0.40 (-40dB à -8dB) en tapant sur le même ordinateur portable
    • 0.10 (-20dB) de claquer des doigts doucement à 1' distance
    • 0.60 (-4.4 dB) à claquer des doigts très fort à 1'

    Mise à jour: voici un exemple pour vous aider à démarrer.

    #!/usr/bin/python
# open a microphone in pyAudio and listen for taps
import pyaudio
import struct
import math
INITIAL_TAP_THRESHOLD = 0.010
FORMAT = pyaudio.paInt16 
SHORT_NORMALIZE = (1.0/32768.0)
CHANNELS = 2
RATE = 44100  
INPUT_BLOCK_TIME = 0.05
INPUT_FRAMES_PER_BLOCK = int(RATE*INPUT_BLOCK_TIME)
# if we get this many noisy blocks in a row, increase the threshold
OVERSENSITIVE = 15.0/INPUT_BLOCK_TIME                    
# if we get this many quiet blocks in a row, decrease the threshold
UNDERSENSITIVE = 120.0/INPUT_BLOCK_TIME 
# if the noise was longer than this many blocks, it's not a 'tap'
MAX_TAP_BLOCKS = 0.15/INPUT_BLOCK_TIME
def get_rms( block ):
# RMS amplitude is defined as the square root of the 
# mean over time of the square of the amplitude.
# so we need to convert this string of bytes into 
# a string of 16-bit samples...
# we will get one short out for each 
# two chars in the string.
count = len(block)/2
format = "%dh"%(count)
shorts = struct.unpack( format, block )
# iterate over the block.
sum_squares = 0.0
for sample in shorts:
# sample is a signed short in +/- 32768. 
# normalize it to 1.0
n = sample * SHORT_NORMALIZE
sum_squares += n*n
return math.sqrt( sum_squares / count )
class TapTester(object):
def __init__(self):
self.pa = pyaudio.PyAudio()
self.stream = self.open_mic_stream()
self.tap_threshold = INITIAL_TAP_THRESHOLD
self.noisycount = MAX_TAP_BLOCKS+1 
self.quietcount = 0 
self.errorcount = 0
def stop(self):
self.stream.close()
def find_input_device(self):
device_index = None            
for i in range( self.pa.get_device_count() ):     
devinfo = self.pa.get_device_info_by_index(i)   
print( "Device %d: %s"%(i,devinfo["name"]) )
for keyword in ["mic","input"]:
if keyword in devinfo["name"].lower():
print( "Found an input: device %d - %s"%(i,devinfo["name"]) )
device_index = i
return device_index
if device_index == None:
print( "No preferred input found; using default input device." )
return device_index
def open_mic_stream( self ):
device_index = self.find_input_device()
stream = self.pa.open(   format = FORMAT,
channels = CHANNELS,
rate = RATE,
input = True,
input_device_index = device_index,
frames_per_buffer = INPUT_FRAMES_PER_BLOCK)
return stream
def tapDetected(self):
print("Tap!")
def listen(self):
try:
block = self.stream.read(INPUT_FRAMES_PER_BLOCK)
except IOError as e:
# dammit. 
self.errorcount += 1
print( "(%d) Error recording: %s"%(self.errorcount,e) )
self.noisycount = 1
return
amplitude = get_rms( block )
if amplitude > self.tap_threshold:
# noisy block
self.quietcount = 0
self.noisycount += 1
if self.noisycount > OVERSENSITIVE:
# turn down the sensitivity
self.tap_threshold *= 1.1
else:            
# quiet block.
if 1 <= self.noisycount <= MAX_TAP_BLOCKS:
self.tapDetected()
self.noisycount = 0
self.quietcount += 1
if self.quietcount > UNDERSENSITIVE:
# turn up the sensitivity
self.tap_threshold *= 0.9
if __name__ == "__main__":
tt = TapTester()
for i in range(1000):
tt.listen()
    • Pourriez-vous poster un simple exemple de code? Je n'ai jamais travaillé avec l'audio avant.
    • +1 - C'est fascinant.
    • Je vous remercie. Cela m'aide grandement.
    • merci lottt!! cela me permet de tonnes !! très instructif. Bien que, est-il possible de prendre le tout automatique notion de seuil et calibrer manuellement il? Comme par exemple, Si je enregistrer des robinets, le bruit, les ruptures, les claquements de mic et de les voir en logiciel, le son a clairement un niveau allant jusqu'à -12 dB, alors que les robinets sont beaucoup plus grandes que -12 dB ou plus de 0 dB ou même plus. Si je veux mettre mon seuil de tolérance à -12 dB. Comment puis-je le faire?
    • il suffit de retirer la logique que les changements de soi.tap_threshold. Selon ce que votre '-12dB' est relatif, il peut ou peut ne pas correspondre à un seuil de 0,25, alors essayez de l'initialisation tap_threshold à cette valeur à la place de l'0,01 dans mon exemple.
    • Python dispose d'un moyen standard de calcul de l'amplitude RMS, croyez-le ou pas: audioop. Vous pouvez remplacer le get_rms fonction ci-dessus avec ceci: def get_rms(block): return audioop.rms(block, 2).
    • Wow, je ne savais pas à propos de audioop. C'est une grave import antigravity action. (Pour une compatibilité totale, il vous reste à redimensionner, mais oui.)
    • Incroyable, je vous remercie!
    • Ce, aussi loin que je peux trouver, c'est le seul exemple montrant comment extraire les valeurs d'échantillon de pyaudio. Cheers!

    InformationsquelleAutor Russell Borogove
  2. 14

    une version simplifiée du code ci-dessus...

    import pyaudio
import struct
import math
INITIAL_TAP_THRESHOLD = 0.010
FORMAT = pyaudio.paInt16 
SHORT_NORMALIZE = (1.0/32768.0)
CHANNELS = 2
RATE = 44100  
INPUT_BLOCK_TIME = 0.05
INPUT_FRAMES_PER_BLOCK = int(RATE*INPUT_BLOCK_TIME)
OVERSENSITIVE = 15.0/INPUT_BLOCK_TIME                    
UNDERSENSITIVE = 120.0/INPUT_BLOCK_TIME # if we get this many quiet blocks in a row, decrease the threshold
MAX_TAP_BLOCKS = 0.15/INPUT_BLOCK_TIME # if the noise was longer than this many blocks, it's not a 'tap'
def get_rms(block):
# RMS amplitude is defined as the square root of the 
# mean over time of the square of the amplitude.
# so we need to convert this string of bytes into 
# a string of 16-bit samples...
# we will get one short out for each 
# two chars in the string.
count = len(block)/2
format = "%dh"%(count)
shorts = struct.unpack( format, block )
# iterate over the block.
sum_squares = 0.0
for sample in shorts:
# sample is a signed short in +/- 32768. 
# normalize it to 1.0
n = sample * SHORT_NORMALIZE
sum_squares += n*n
return math.sqrt( sum_squares / count )
pa = pyaudio.PyAudio()                                 #]
#|
stream = pa.open(format = FORMAT,                      #|
channels = CHANNELS,                          #|---- You always use this in pyaudio...
rate = RATE,                                  #|
input = True,                                 #|
frames_per_buffer = INPUT_FRAMES_PER_BLOCK)   #]
tap_threshold = INITIAL_TAP_THRESHOLD                  #]
noisycount = MAX_TAP_BLOCKS+1                          #|---- Variables for noise detector...
quietcount = 0                                         #|
errorcount = 0                                         #]         
for i in range(1000):
try:                                                    #]
block = stream.read(INPUT_FRAMES_PER_BLOCK)         #|
except IOError, e:                                      #|---- just in case there is an error!
errorcount += 1                                     #|
print( "(%d) Error recording: %s"%(errorcount,e) )  #|
noisycount = 1                                      #]
amplitude = get_rms(block)
if amplitude > tap_threshold: # if its to loud...
quietcount = 0
noisycount += 1
if noisycount > OVERSENSITIVE:
tap_threshold *= 1.1 # turn down the sensitivity
else: # if its to quiet...
if 1 <= noisycount <= MAX_TAP_BLOCKS:
print 'tap!'
noisycount = 0
quietcount += 1
if quietcount > UNDERSENSITIVE:
tap_threshold *= 0.9 # turn up the sensitivity
    • Sans input_device_index dans pyaudio.PyAudio().open(... ) voulez-vous obtenir le silence ou serait pyaudio en quelque sorte recherchez un travail mic?
    InformationsquelleAutor user1405612