Discussion :

[Sebsheep]

Bonjour,

Pour tracer de jolies fft avec scipy, j'ai besoin de récupérer le signal sonore d'un fichier audio (mp3,wav...). Après quelques recherches, j'ai vu qu'il était possible d'utiliser pymedia pour acquérir le son puis le jouer, mais je n'ai pas trouvé comment obtenir le signal sous forme d'un tableau ou d'une liste afin de passer le boulot ensuite à scipy.

Quelqu'un peut-il m'aider ?

Merci,
Sébastien

[Desintegr]

Tu peux peut-être simplement utiliser le module wave [1] et la méthode Wave_read.readframes(n)

Wave_read.readframes(n)
Reads and returns at most n frames of audio, as a string of bytes.

[1] http://docs.python.org/library/wave.html

[Sebsheep]

merci, ça à l'air de marcher ...


mais (sinon ce serait trop simple), j'ai des petits soucis:
j'ouvre mon fichier, pas de problèmes : son=wave.open("test.wav",'r')
Je récupère la string de bit : chaine = son.readframes(100)

Et là, je sais pas si je fais ce qu'il faut, pour avoir une liste de valeurs, je fais:
liste = [ord(c)-128 for c in chaine]

J'ai mis le "-128", car j'ai remarqué que sinon les valeurs allaient de 0 à 255.

Après, j'affiche la liste sur un graphe grâce à un plot (from pylab, mais ca on s'en tape), et le graphe que j'obtiens ne correspond pas du tout au graphe que me donne audacity.

Autre curiosité, je lui demande 100 frames, et la liste que j'ai est de 200 éléments.

Merci de ton aide

[pacificator]

Autre curiosité, je lui demande 100 frames, et la liste que j'ai est de 200 éléments.

Au hasard, ton fichier son est en stéreo?

[Sebsheep]

Non il est bien en mono.

Mais j'ai un peu avancé. Je me suis dit "Et si j'arretai de faire le kéké à vouloir prendre n'importe quel son, si je prenais un son sinusoidal" ... et bizarement, on voit plus de choses.

Je me suis rendu compte que le signal qui m'interessait était en fait sur les cases 'impaires', donc pour avoir la courbe qui m'interesse, je dois faire :
liste[1::2]
Mais, là encore un petit détail, la conversion caractère -> int ne se fait pas comme il faut, j'ai du appliquer : if x >128 : return x-255 ; else return x à ma liste. (en gros le bit de poid fort n'est pas interprété comme un bit de signe). Y a t il une solution à ça ?

Mais dans ce cas là, à quoi correspond le cases paires ? Je l'ai ploté, et ca donne un truc comme du bruit (aucune période, truc qui parait relativement aléatoire), bruit que je n'ai absolument pas dans mon fichier : ce fichier contient un son formé d'une sinusoïde pure

[Desintegr]

À mon avis, c'est un peu plus compliqué que ça.

Quand tu utilises readframes(n) ça te renvoie une chaîne d'octets pour n frames.
Cette chaîne est, je pense, de taille nombre de frames * nombre de canaux * format d'échantillonnage (en octets).
Par exemple, un fichier avec 10 frames en stéréo échantillonné sur 16 bits renverra une chaîne de taille 10 * 2 * 2 = 40.

D'après [1] :

8-bit samples are stored as unsigned bytes, ranging from 0 to 255. 16-bit samples are stored as 2's-complement signed integers, ranging from -32768 to 32767.

Il faut donc extraire les données de cette chaîne. Pour les fichiers wav échantillonnés sur 8 bits, il faut utiliser un unsigned char (0 à 255), et pour les wav échantillonnés sur 16 bits, il faut utiliser un signed short (-32768 à 32767). Tu peux utiliser struct pour faire ça :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
import wave, struct
 
def readwave(filename):
	wav = wave.open(filename,'rb')
	(nchannels, sampwidth, framerate, nframes, comptype, compname) = wav.getparams()
	frames = wav.readframes(nframes)
        # http://docs.python.org/library/struct.html
        if sampwidth == 1:
          # 8 bit : unsigned char
          data = struct.unpack('%sB' % (nframes * nchannels), frames)
        elif sampwidth == 2:
          # 16 bits : signed short
	  data = struct.unpack('%sh' % (nframes * nchannels), frames)
        return data

Ensuite, il reste à découper les données restantes selon le nombre de canaux. Si c'est du mono, il n'y a rien à faire. Si c'est du stéréo, il faut lire un entier sur deux pour chaque piste :
[gauche, droit, gauche, droit, gauche, droit...]

Enfin, il ya peut-être également une histoire de little-endian et big-endian.

[1] http://ccrma.stanford.edu/courses/42...ts/WaveFormat/

[Sebsheep]

Merci beaucoup, ça marche très bien ! Je connaissais pas ce module struct de python, très utile pourtant !

Par contre, je ne comprend pas pourquoi tu prends des unsigned char pour les 8 bits (pas de fichier 8 bits sous la main pour tester). Si on prend unsigned, ne va t on pas avoir un signal mal interprété ? (un peu comme j'avais en fait?)

[Desintegr]

J'ai simplement suivi ce que dit la doc :

8-bit samples are stored as unsigned bytes, ranging from 0 to 255

À tester pour voir si ça fonctionne bien avec des fichiers échantillonnés sur 8 bits.

[Desintegr]

Pour info dans Scipy 0.7 (qui n'est pas encore sorti), il y aura une méthode io.wavfile.read [1] peut-être plus simple à utiliser.

[1] http://docs.scipy.org/doc/scipy/refe...file.read.html

[Sebsheep]

oui mais j'en ai besoin maintenant ... et cette méthode marche, donc je ne vais pas attendre la sortie de scipy 0.7 :p

Merci du tuyau !

Seb