Brevet FR3022391A1 — publié le 20151218

PROCEDE DE TRAITEMENT DE DONNEES ACOUSTIQUES CORRESPONDANT A UN SIGNAL ENREGISTRE DOMAINE TECHNIQUE Il est question du domaine technique du traitement des signaux acoustiques. Il est en particulier question d'un procédé et d'un système 400 permettant un tel traitement. ETAT DE LA TECHNIQUE Les signaux acoustiques, en particulier les signaux sonores font l'objet d'enregistrements afin d'être analysés. Afin de tirer des informations de ces enregistrements, il est connu d'appliquer une transformée de Fourier convertissant le signal du domaine temporel au domaine spectral. De nombreux types de transformée ont été élaborées afin de rendre plus accessibles certaines informations recherchées. Cependant, les techniques d'analyse fréquentielles connues offrent une résolution limitée, notamment en raison du principe d'incertitude temps-fréquence. Ceci limite d'autant les analyses possibles des enregistrements, et donc les applications envisageables. Par exemple, dans le domaine de l'analyse sonore, il est difficile d'estimer précisément des hauteurs de notes à partir de sons réels. De plus, la hauteur d'une note est un critère qui devient difficilement définissable dès que le signal correspondant n'est pas parfaitement périodique, par exemple en cas d'inharnnonicité. Il en résulte qu'il est difficile de fournir des informations précises, en particulier de représenter de telles informations. Il en résulte notamment qu'il est difficile de réaliser des accordeurs électroniques fournissant des informations fiables et aisées à lire. PRESENTATION GENERALE Un but de l'invention est de pallier au moins un des inconvénients présentés ci-avant. Un but de l'invention est de permettre une analyse fréquentielle plus fine. A cet effet, il est prévu un procédé de traitement de données acoustiques, comprenant des étapes de: réception de données acoustiques correspondant à un signal acoustique transformé dans le domaine fréquentiel à partir d'un signal acoustique dans le domaine temporel enregistré par au moins un capteur acoustique, pour au moins une fréquence du domaine fréquentiel, démodulation de la phase du signal acoustique transformé, et calcul d'un indice de constance de la phase démodulée. L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leurs combinaisons techniquement possibles : - la phase démodulée est obtenue par une formule du type : (1)demod(f, , t) = I`I 1) spectr (f , t) - 271- f t OÙ Odemod représente la phase démodulée pour la fréquence f et au temps t, et où (Dspectr représente la phase du signal transformé pour la pour la fréquence f et au temps t; le calcul d'un indice de constance de la phase démodulée comprend l'application d'un filtre passe-bas, le filtre passe-bas étant adapté pour réaliser une atténuation en fonction de la vitesse de rotation de la phase démodulée ; le calcul d'un indice de constance de la phase démodulée comprend, préalablement à l'application du filtre passe-bas la conversion de la phase démodulée en une exponentielle complexe de module constant, le filtre passe-bas étant appliqué à l'exponentielle complexe, l'indice de constance de la phase démodulée résultant du module d'une sortie complexe du filtre passe-bas ; des étapes de : o réception d'une valeur de variation de fréquence à détecter ; o réception de données acoustiques correspondant au signal acoustique dans le domaine temporel enregistré par l'au moins un capteur acoustique ; o transformation du signal acoustique dans le domaine temporel en le signal dans le domaine fréquentiel par transformée de Fourier discrète pondérée par une fenêtre d'analyse ; la transformée étant réalisée avec un pas d'avancement déterminé de sorte que l'inverse du double du pas d'avancement est supérieur a: o la largeur du lobe principal de la fenêtre d'analyse dans le domaine de Fourier, et à o la valeur de variation à détecter ; le signal acoustique enregistré est un signal sonore et dans lequel le signal acoustique dans le domaine temporel enregistré est échantillonné ; une étape de construction d'une image adaptée pour être affichée sur un terminal, l'image étant construite à partir des données acoustiques et de l'indice de constance de phase démodulée. la construction de l'image comprend un élément visuel correspondant à la fréquence du domaine fréquentiel, la représentation de l'élément visuel étant fonction de l'amplitude de la fréquence et de l'indice de constance de phase démodulée. L'invention concerne également un système de traitement de données acoustiques, adapté pour mettre en oeuvre un tel procédé, comprenant : un module de réception de données acoustiques correspondant à un signal acoustique transformé dans le domaine fréquentiel à partir d'un signal acoustique dans le domaine temporel enregistré par au moins un capteur acoustique ; un module de démodulation de la phase du signal acoustique transformé pour au moins une fréquence du domaine fréquentiel ; et un module de calcul d'un indice de constance de la phase démodulée. FIGURES D'autres objectifs, caractéristiques et avantages apparaitront à la lecture de la description qui suit donnée à titre illustratif et non limitatif en référence aux figures, parmi lesquels : la figure 1 représente un procédé selon un autre exemple de mode de réalisation de l'invention, la figure 2 représente une étape de réalisation d'une image selon un autre exemple de mode de réalisation de l'invention, la figure 3 représente une représentation visuelle réalisée par l'étape de la figure 2, la figure 4 représente un système selon un autre exemple de mode de réalisation de l'invention, les figures 5a et 5b représentent un traitement des harmoniques selon un autre exemple de mode de réalisation de l'invention. DESCRIPTION DETAILLEE Procédé Réception de données acoustiques En référence à la figure 1, il est décrit un procédé de traitement de données acoustiques. Le procédé est par exemple mis en oeuvre par un système 400 de traitement de données acoustiques, par exemple un système 400 tel que représenté figure 4 et tel que décrit plus en détail plus loin. Le procédé comprend une étape 110 de réception de données acoustiques transformées. Les données acoustiques correspondent à un signal acoustique transformé dans le domaine fréquentiel à partir d'un signal acoustique dans le domaine temporel. Le signal acoustique dans le domaine temporel est par exemple un signal acoustique enregistré dans le domaine temporel par au moins un capteur acoustique. Le signal acoustique enregistré peut être un signal infrasonore, sonore ou ultrasonore. Le signal acoustique enregistré est ainsi par exemple un signal sonore. Le signal acoustique enregistré peut être un signal monophonique. Le signal acoustique dans le domaine temporel enregistré peut être échantillonné. L'échantillonnage peut être réalisé à une fréquence déterminée, par exemple 44100 Hz. Le signal acoustique dans le domaine temporel comprend par exemple une pluralité de trames, par exemple de trames successives ayant chacune une durée déterminée, par exemple une même durée de l'ordre de 50 ms. Etapes préalables à la réception Réception de données acoustiques correspondant au signal acoustique dans le domaine temporel Les données reçues à l'étape 110 peuvent être des données stockées au préalable (non représenté). Alternativement, comme illustré figure 1, le procédé peut comprendre une étape 104 de réception de données acoustiques correspondant au signal acoustique dans le domaine temporel enregistré par l'au moins un capteur acoustique. Le procédé peut alors comprendre une étape 106 d'application d'un gain au signal acoustique dans le domaine temporel. Le gain est par exemple un gain multiplicatif. Le gain peut être appliqué de sorte à constituer un réglage de niveau pour le confort d'un utilisateur utilisant ces données ou destiné à utiliser ces données. Le procédé peut comprendre une étape 108 de transformation du signal acoustique dans le domaine temporel pour obtenir le signal dans le domaine fréquentiel. L'étape 108 de transformation peut être réalisée suite à l'étape 106 d'application de gain ou en l'absence d'étape d'application de gain. Une telle transformation 108 est par exemple réalisée par transformée de Fourier, par exemple par transformée de Fourier discrète. La transformée de Fourier discrète est par exemple appliquée à chaque trame de signal acoustique. La transformée de Fourier est par exemple une transformée de Fourier rapide. La transformée de Fourier peut être pondérée par une fenêtre d'analyse. La fenêtre d'analyse est par exemple une fenêtre de Nanning. La fenêtre d'analyse peut être une fenêtre dépourvue de lobe secondaire, par exemple sur un nombre de points étant une puissance de deux, par exemple 4096, par exemple une fenêtre telle que décrite dans T. Hélie - Amélioration de l'extraction de partiels dans les signaux sonores », Rapport de Stage Ingénieur, 1996. Une telle fenêtre permet de ne plus perturber l'interprétation visuelle des lobes puisque ceux-ci ne peuvent plus correspondre à des lobes secondaires. La transformée de Fourier peut être réalisée avec un pas d'avancement p déterminé. Le pas d'avancement est par exemple déterminé de sorte à vérifier une condition du type : 1 1 Le procédé peut comprendre une étape de réception d'une variation de fréquence (vibrato, glissando, etc.) de valeur de variation de fréquence y à détecter. Le pas d'avancement peut alors en complément, être déterminé de sorte à vérifier une condition du type : 1 v Une telle analyse temps/fréquence, en particulier une telle transformation, peut être réalisée conformément à L. Cohen - Tinne-Frequency Analysis », Prentice-Hall, New York, 1995. Enregistrement d'un signal acoustique dans le domaine temporel Les données acoustiques correspondant au signal acoustique dans le domaine temporel enregistré par l'au moins un capteur acoustique reçues à l'étape 104 peuvent être des données stockées au préalable (non représenté). Alternativement, comme illustré figure 1, le procédé peut comprendre une étape 100 d'enregistrement d'un signal acoustique dans le domaine temporel par l'au moins un capteur acoustique. Il peut s'agir par exemple d'un enregistrement à la volée. Le procédé peut également comprendre une étape 102 d'échantillonnage du signal acoustique enregistré. L'échantillonnage peut être réalisé à une fréquence déterminée, par exemple 44100 Hz. Démodulation de phase Le procédé comprend une étape 118 de démodulation de la phase du signal acoustique transformé. La démodulation 118 est réalisée pour une ou plusieurs fréquences du domaine fréquentiel. L'étape 118 de démodulation peut comprendre une étape 1180 d'encodage de la phase, par exemple d'encodage par une valeur, par exemple d'encodage par une valeur rationnelle, par exemple par un entier, par exemple un entier positif. La phase du signal acoustique transformé est par exemple encodée avec un nombre de bits kb,t' par exemple 32. La formule de passage entre la valeur en radians (D rad rad de la phase, comprise par exemple entre 0 et 2u, et l'entier encodant la phase (1),'t est par exemple du type : (I) int = FlOOr ((1) Int X 21 blts 71") où Floor(.) est à la partie entière arrondie par défaut. Il est ainsi possible de conserver la classe d'équivalence « modulo 2u » par cet encodage. La phase démodulée peut être obtenue par une formule du type : (1)demod(f, t) (I)spectr(f t) 271" f t mod 2n- Ù Odemod représente la phase démodulée pour la fréquence f et au temps t, et oÙ Ospectr représente la phase du signal transformé pour la fréquence f et au temps t, par exemple le temps t associé à une trame en cours, par exemple le temps de début de trame, le temps t avançant par exemple au temps d'avancement. La phase démodulée peut être obtenue par une telle formule appliqué de manière discrète à un temps t' obtenu par incrénnentation d'un pas d'avancement. Le temps t' pour la nèrne incrénnentation est par exemple obtenu par une formule du type : tn = n X t' où tav est le pas d'avancement, et est par exemple constant à chaque incrénnentation. L'étape de démodulation peut ainsi comprendre une étape d'application 30 d'une formule de démodulation pour cet encodage, du type : (1)demod(f, en) = (1)spectr(f, tn) ,Ademod où le vecteur Ademod est rafraichi à chaque temps tr, par un formule du type : ,Ademod= ,Ademod où à est l'encodage entier positif de -2nftav, pour chaque fréquence analysée f, le pas d'avancement ta, étant par exemple fourni par un paramètre stocké éventuellement converti en une unité de temps utilisée, par exemple en secondes. Le terme -2nfta, peut n'être calculé qu'une seule fois, par exemple à chaque nouveau réglage d'une fréquence de référence par l'utilisateur. De plus, lorsque l'encodage des phases est réalisé sur au moins 32 bits, l'erreur cumulée dans IDdemod via sa formule de mise à jour n'entraîne une déviation de 2n qu'après une durée supérieure à 2khns ±1p, soit au moins 11930 heures pour un pas d'avancement p de 5nns et un ka,t, de 32 bits. Une telle démodulation a donc une précision qui dépasse très largement l'ordre de grandeur requis qui est par exemple de l'ordre de la minute. Calcul d'un indice de constance de la phase démodulée Le procédé comprend une étape 120 de calcul d'un indice de constance de la phase démodulée. L'indice de constance de la phase démodulée peut ainsi être fourni en sortie d'une telle étape 120. Le procédé permet ainsi de construire un réhausseur de stationnarité fréquentielle des composantes spectrales de données acoustiques, en particulier d'un son. En effet, le calcul d'un indice de constance de phase démodulée permet de mettre en évidence de manière efficace qu'une composante spectrale correspond parfaitement à une fréquence d'analyse. La déviation dans le temps de la phase démodulée est par exemple proportionnelle à la déviation de fréquence entre la fréquence de la composante analysée et la fréquence de démodulation. En effet, pour une composante xi(t) = Aie2i1l1t+`P1, d'amplitude A1, de fréquence fixe fl, et de phase (pi à l'instant t = 0 présente dans le signal complet, la phase instantanée (12.1(t) = 2inf1t + (pi démodulée à une fréquence d'analyse f est (12. 1 demod(t, f) = Arexi(t) x e'in-ft) = 2n-t(f1 - n + (pi . Un tel résultat met en évidence que 'PI 1 demod est cyclique, par exemple modulo 2n, de fréquence fi-f, c'est-à-dire la déviation de fréquence entre la composante analysée fi et la fréquence de démodulation f. Le calcul 120 de l'indice de constance de la phase démodulée peut comprendre l'application 1202 d'un filtre passe-bas. Le filtre passe-bas est par exemple adapté pour réaliser une atténuation en fonction de la vitesse de variation de la phase démodulée, par exemple la vitesse de rotation de la phase démodulée. Le calcul 120 de l'indice de constance de la phase démodulée peut comprendre, préalablement à l'application 1202 du filtre passe-bas, la conversion 1200 de la phase démodulée en une exponentielle complexe de module constant. Le filtre passe-bas est par exemple appliqué à l'exponentielle complexe. Le signal complexe ainsi construit reste constant si la phase démodulée est constante et tourne sur le cercle unité si la phase démodulée évolue. La conversion 1200 de la phase démodulée en une exponentielle complexe peut être réalisée via une table de conversion. La table de conversion est par exemple une table de valeurs pré-calculées avec une résolution d'encodage des phases déterminée, par exemple de l'ordre de 10 bits. A cet effet, des bits de poids fort, par exemple 10 bits, sont prélevés de l'entier, par exemple encodé sur 32 bits. Pour les 210 valeurs N possibles entre 0 et 210-1, on a par exemple pré-calculé et stocké dans une telle table de valeurs pré-calculées des valeurs complexes e2urN / 121° encodées par exemple par deux flottants ou deux doubles. Il est ainsi possible d'éviter le calcul de l'exponentielle complexe et de limiter à 210 le nombre de paires de nombres, par exemple de flottants, à stocker en mémoire. Le filtre passe-bas est par exemple un filtre de type Butterworth, par exemple tel que décrit dans : A. V. Oppenheinn, R. W. Schafer, Digital Signal Processing », 3rd edition, 2011. Prentice-Hall. Un tel filtre est par exemple de gain le plus plat possible en dessous de la fréquence de coupure, par exemple un gain proche de 1. La fréquence de coupure est par exemple de l'ordre de 6 Hz. La fréquence de coupure peut ainsi permettre de régler la qualité de constance de phase souhaitée. Le filtre passe-bas peut être appliqué de sorte que l'indice de constance de la phase démodulée résulte du module d'une valeur complexe formant une sortie complexe fournie en sortie du filtre passe-bas. L'indice de constance de la phase démodulée peut ainsi par exemple être le module de la sortie du filtre. La valeur de l'indice peut par exemple être de l'ordre du module de l'entrée (valeur 1) si le signal d'entrée tourne à une fréquence inférieure à la fréquence de coupure, et être atténuée au-delà, d'autant plus que l'argument du signal complexe d'entrée tourne vite et que l'ordre du filtre est élevé. La sortie du filtre passe-bas peut être saturée si son module dépasse une valeur déterminée, par exemple le module d'entrée. La phase étant définie modulo 2u, une plus grande déviation de fréquence détectable sans repliement spectral par le calcul de la déviation de phase démodulée est par exemple de l'ordre de la moitié de l'inverse du pas d'avancement p, par exemple une déviation de l'ordre de 100Hz pour un pas d'avancement de 5 ms. La plus grande déviation détectable est par exemple déterminée de sorte à être supérieure aux variations musicales que l'utilisateur souhaite considérer. Ces variations sont par exemple les vibratos et/ou les battements. La plus grande déviation détectable est par exemple déterminée de sorte à être supérieure à la largeur du lobe principal de la fenêtre d'analyse dans le domaine de Fourier. Ainsi, le pas d'avancement peut être déterminé de sorte à vérifier une condition du type : 1 1 Par exemple, en permettant de construire un indice descripteur de constance de phase proche de 1 lorsque les phases démodulées ont des cycles, par exemple de type modulo 2u, par exemple de fréquence fi-f, inférieure à la résolution fréquentielle finale souhaitée, et s'approchant rapidement de 0 lorsque la fréquence est supérieure à la résolution fréquentielle finale souhaitée, il est possible de caractériser la plage prédéterminée et d'atteindre la résolution fréquentielle cible . En particulier, pour chaque fréquence f calculée et à analyser, le procédé permet de construire la phase démodulée du spectre estimé à la fréquence f pour chaque temps de trame t'. La fréquence f appartenant par exemple à une grille des fréquences, par exemple une grille de fréquences accordées sur un diapason. Lorsque la composante spectrale pour f est parfaitement stable et accordée sur la fréquence f, la vitesse de rotation de la phase démodulée correspondante calculée est nulle. Lorsque la composante spectrale de fréquence stable est accordée sur f-i-5f, c'est-à-dire qu'elle est à peine décalée, avec par exemple df de l'ordre de 1 Hz, l'amplitude spectrale correspondante obtenue par transformée de Fourier reste significativement grande et peu modifiée, notamment à cause de la taille du lobe principal et du principe d'incertitude. Cependant, la vitesse de rotation de la phase démodulée, a une valeur fonction du décalage, par exemple fonction linéaire du décalage, par exemple égale au décalage. Le procédé permet ainsi de construire un indice caractéristique de la résolution fréquentielle souhaitée. Selon un exemple, un signal complexe de module 1 avec une phase démodulée tournant sur le cercle unité à la fréquence df peut être filtré par un filtre passe-bas. Si la fréquence de coupure est plus haute que df, le signal de sortie est par exemple similaire à celui de l'entrée. Dans le cas contraire, le signal de sortie a par exemple un module significativement atténué. Le module, forme par exemple l'indice de constance de phase démodulée, et forme ainsi un indicateur permettant de déterminer si une composante spectrale est ou non dans une plage de résolution fréquentielle déterminée. Il est ainsi possible d'obtenir une résolution fréquentielle affinée qui permet de dépasser le problème d'incertitude temps/fréquence. L'indice de constance de phase permet d'être plus sélectif que la largeur du lobe principal de la fenêtre d'analyse de Fourier tout en indiquant quelles sont les composantes spectrales pour lesquelles la phase démodulée est cohérente avec la fréquence d'analyse. Interpolation par un vecteur de fréquences Le procédé peut en outre comprendre une étape 112 de réception d'un vecteur de fréquences. Les fréquences du vecteur de fréquence sont par exemple des fréquences accordées sur un diapason, c'est-à-dire que les fréquences sont accordées sur une note de référence. La fréquence de la note de référence est par exemple 440 Hz. La fréquence de la note référence peut par exemple être modifiée par l'utilisateur. Les fréquences du vecteur de fréquence peuvent par exemple correspondre aux notes de musique de l'échelle tempérée (musique européenne classique). Chaque note de musique peut par exemple être représentée par son code midi, par exemple selon la spécification internationale, par exemple telle que définie par - Complete MIDI 1.0 Detailed Specification La correspondance entre une fréquence f et le code midi m d'une note est par exemple définie par la formule suivante : (m-6" f = f69 X 2 /12 où 69 est le code midi de la note 1a3 et f69 vaut est une valeur de référence valant par exemple 440 Hz Les fréquences du vecteur de fréquence peuvent alternativement correspondre aux notes de musique d'autres échelles musicales, issues par exemple de la musique orientale, ou d'échelles micro-tonales. Le vecteur de fréquence est par exemple calculé une seule fois à l'initialisation du système 400. En cours de fonctionnement du système 400 le vecteur de fréquences peut n'être réinitialisé qu'a chaque nouveau réglage de la fréquence de référence par l'utilisateur. Le procédé peut comprendre une étape 114 de calcul par interpolation des valeurs du signal acoustique transformé aux fréquences du vecteur de fréquences. Une telle étape 114 peut être réalisée préalablement à la démodulation 118 de phase. Le procédé permet ainsi de construire un réhausseur de stationnarité 15 fréquentielle des composantes spectrales de données acoustiques, en particulier d'un son, pour une grille de fréquences, par exemple une grille de fréquences accordées sur un diapason. Un tel calcul 114 peut être réalisé de sorte que la démodulation 118 et le calcul 120 d'indice de constance de phase sont réalisés pour chaque 20 fréquence du vecteur de fréquences. Ainsi, l'indice de constance de phase permet de marquer, pour chaque fréquence f, la pertinence de ladite fréquence, par exemple pour un indice proche de 1, ou son absence de pertinence, par exemple pour un indice proche de 0, en tant que composante sinusoïdale de fréquence f, et ceci, quelle que 25 soit l'amplitude associée à la fréquence f. L'interpolation 114 est par exemple une interpolation affine par morceaux. Les valeurs peuvent être des valeurs complexes obtenues par interpolation dans un plan complexe. 30 Calcul de coordonnées polaires Le procédé peut comprendre une étape 116 de calcul des coordonnées polaires du signal acoustique transformé. L'étape 116 de calcul des coordonnées polaires peut être réalisée préalablement à l'étape 118 de démodulation de phase. L'étape 116 de calcul des coordonnées polaires peut être réalisée suite à l'étape 114 de calcul par interpolation. Les coordonnées polaires comprennent par exemple le module et la phase modulo 2u des valeurs complexes du signal acoustique transformé associées à chaque fréquence du vecteur de fréquences. Réalisation d'une image Le procédé décrit peut ainsi par exemple permettre de fournir au moins un vecteur de données, par exemple une pluralité de vecteurs de données, par exemple des vecteurs de données correspondant au vecteur de fréquences et/ou au module des valeurs complexes du signal acoustique transformé associées à chaque fréquence du vecteur de fréquences et/ou à la phase démodulée, et/ou à l'indice de constance de phase. Au moins un des vecteurs de données, par exemple plusieurs ou tous, peut par exemple correspondre à un tableau de valeurs d'analyse de trames de signal acoustique ayant une durée déterminée, par exemple de l'ordre de 50 20 ms. Un ou plusieurs des vecteurs de données peut par exemple être rafraîchi à une cadence déterminée, par exemple par le pas d'avancement p du procédé, qui est par exemple de l'ordre de 5 ms. Le procédé peut comprendre en outre une étape 122 de réalisation d'au 25 moins une image adaptée pour être affichée sur un terminal, l'image étant réalisée au moins à partir des données acoustiques et de l'indice de constance de phase démodulée. Il est ainsi possible de combiner l'analyse fréquentielle des données acoustiques et un affichage, par exemple un affichage en temps réel. Il est 30 ainsi possible de visualiser l'information spectrale obtenue en sortie des étapes précédemment décrites. Le procédé peut ainsi comprendre une étape 122 de réalisation de plusieurs telles images. Chaque image peut par exemple correspondre à des données calculée pour un temps t, par exemple un avancement t'. Le système 400 comprend par exemple le terminal. Le terminal comprend par exemple un afficheur. L'afficheur peut par exemple fonctionner avec une fréquence de rafraîchissement déterminée propre à l'afficheur, la fréquence de rafraîchissement étant par exemple comprise entre 50 Hz et 122 Hz. Réception de données En référence à la figure 2, l'étape 122 de réalisation de l'image comprend par exemple une étape 200 de réception de données. Les données comprennent par exemple le vecteur de fréquences reçu, et/ou le module des valeurs complexes du signal acoustique transformé associées à chaque fréquence du vecteur de fréquences, et/ou la phase démodulée, et/ou l'indice de constance de phase démodulée. Les données sont par exemple reçues par paquets, par exemple par paquets ordonnés. Chaque paquet correspond par exemple à une trame 20 analysée au cours des étapes décrites ci-avant. Chaque paquet peut comprendre les quatre vecteurs de sortie de l'analyseur. Ces données peuvent alors être traitées pour construire l'image à afficher. Les données de plusieurs paquet peuvent être simultanément prises en compte pour construire une même image, par exemple pour un afficheur à 25 mémoire, par exemple un tel afficheur de type sonagramme. Les données d'un unique paquet peuvent être prises en compte pour construire une image, par exemple pour un afficheur instantané de type spirale tel que décrit plus loin. Normalisation de l'amplitude 30 L'étape 122 de réalisation de l'image comprend par exemple une étape 202 de conversion du module des valeurs complexes du signal acoustique transformé associées à chaque fréquence du vecteur de fréquences sur une échelle et/ou une étape 204 de normalisation du module des valeurs complexes du signal acoustique transformé associées à chaque fréquence du vecteur de fréquences. La conversion 202 est par exemple une conversion du module sur une échelle adaptée à la perception pour chaque fréquence du vecteur de fréquences. L'échelle est par exemple une échelle en sonies comprenant par exemple des courbes isosoniques normales de la norme ISO 226:2003. Alternativement, l'échelle est par exemple une échelle en décibels. Le système 400 peut permettre de sélectionner une échelle parmi l'échelle en sonie et l'échelle en décibels. La normalisation 204 comprend par exemple une normalisation entre une valeur minimale, par exemple 0, et une valeur maximale, par exemple 1. La normalisation comprend par exemple en une saturation à la valeur minimale de normalisation au-dessous d'un seuil bas, par exemple 0 phones ou 20dB-spl, et/ou une saturation à la valeur maximale de normalisation au-dessus d'un seuil haut, par exemple 100 phones ou 122dB-spl. Elément visuel En référence à la figure 3, l'image 300 peut comprendre un élément visuel 302 correspondant à la fréquence du domaine fréquentiel. L'élément visuel 302 peut correspondre à la fréquence du domaine fréquentiel pour un temps t. La représentation de l'élément visuel 302 peut être fonction de l'amplitude de la fréquence et de l'indice de constance de phase démodulée. Un tel élément visuel 302 permet une visualisation des données traitées, en particulier de l'indice de constance de phase démodulée. Une telle visualisation peut ainsi être reliée à un percept auditif. ill est ainsi possible de stimuler des mises en relations sonores/visuelles, et donc de permettre des applications du procédé à la pédagogie, par exemple la pédagogie musicale en vue d'affiner la justesse de l'oreille de façon interactive. Un tel élément visuel 302 permet l'obtention d'une représentation aisément intelligible pour la personne visualisant l'image, qui peut par exemple être appréhendée par des enfants et/ou sans mode d'emploi. En outre, l'utilisation d'une résolution fréquentielle plus précise qu'une analyse spectrale standard pour des notes établies stables permet une utilisation efficace du procédé en fonction accordeur. Il est ainsi possible de réaliser au même moment un stimulus auditif et un stimulus visuel. La mise en oeuvre d'un tel procédé peut ainsi être en cohérence avec l'échelle perceptive selon laquelle on ne voit que ce qu'on entend. Il est ainsi possible d'obtenir un procédé aisé à appréhender pour l'utilisateur et permettant une compréhension instinctive des données calculées précédemment, en particulier par l'utilisation de l'indice de constance de phase démodulée. Il est ainsi possible d'obtenir une visualisation spectrale du son, qui reste discriminante si elle organisée par hauteurs, par exemple par notes, alors qu'une analyse de Fourier standard serait trop grossière. Il est ainsi possible d'accéder à des résolutions suffisamment fines pour permettre d'accorder un instrument de musique. L'image peut comprendre une grille 304. L'élément visuel 302 est par exemple disposé devant la grille 304. La représentation de l'élément visuel 302 peut comprendre différentes caractéristiques visuelles. Les caractéristiques peuvent comprendre une épaisseur et/ou une couleur et/ou une transparence. L'élément visuel 302 comprend par exemple une forme 302, par exemple un trait, par exemple un trait 302 d'une spirale 306 comprenant une pluralité de traits 302, chaque trait 302 étant associé à une fréquence du vecteur de fréquences. Dans le cas d'une spirale 306, la grille 304 est par exemple une grille 304 circulaire ou en forme de spirale 306. L'étape 122 de réalisation de l'image 300 comprend par exemple une étape 206 de calcul de l'épaisseur de la forme 302, par exemple du trait. L'épaisseur est par exemple normalisée entre une valeur minimale ou nulle, par exemple 0, et une valeur maximale, par exemple 1. L'épaisseur peut ainsi permettre de faire ressortir les composantes sinusoïdales présentes, avec une résolution fréquentielle suffisamment fine pour être au moins proche de celle de la justesse détectable par l'oreille. L'épaisseur de la forme 302, par exemple du trait est par exemple une fonction croissante du module normalisé et/ou de l'indice de constance de la phase démodulée. Pour chaque fréquence f du vecteur de fréquence, l'épaisseur T peut par exemple être donnée par une formule du type : Stress T = florin X (,-f Stress où M est le module normalisé,CD - const l'indice de constance de phase démodulée, et Stress un paramètre déterminé compris par exemple entre 0 et 2, par exemple entre 0 exclus et 2, par exemple réglable entre 0 et 2, par exemple de valeur 1. Lorsque Stress vaut 0, l'épaisseur correspond exactement au module normalisé. Lorsque le paramètre Stress est strictement positif, l'indice de constance de phase démodulée vient réhausser les composantes spectrales à fréquence stationnaires dans la plage tolérée par le filtre passe-bas. Par exemple, pour un paramètre Stress de valeur 2, l'épaisseur des composantes non stationnaires est significativement réduite. L'étape 122 de réalisation de l'image 300 comprend par exemple une étape 208 de calcul de couleur de forme 302, par exemple du trait, par exemple pour chaque fréquence du vecteur de fréquence. L'étape de calcul de couleur peut comprendre le calcul 2080 d'une brillance et/ou le calcul 2084 d'une teinte et/ou le calcul 2082 d'une saturation. En particulier, la brillance est par exemple calculée fonction du module normalisé. La forme 302, par exemple le trait, est par exemple d'autant plus visible que l'amplitude pour l'échelle associée, par exemple en décibel ou en sonie, est grande. La brillance est par exemple calculée sans prendre en compte l'effet de la stationnarité de la fréquence, en particulier indépendamment de l'indice de constance de phase démodulée. Ainsi, quel que soit le type de composante (sinusoïdal, bruit, impact, etc), un évènement deux fois plus fort (échelle de la sonie) sera deux fois plus brillant (échelle perceptive de la vision). A sonie donnée, la brillance est la même quel que soit le type de composante. Par ailleurs, une composante sinusoïdale sera plus marquée géométriquement (- épaisseur du trait plus ample ») qu'une composante d'un autre type. En particulier, la teinte est par exemple calculée en fonction de la phase démodulée. En outre, la teinte est par exemple calculée en respectant une condition de circularité, de sorte qu'une phase nulle et une phase de 3600 sont représentées par la même teinte. En particulier, la saturation est par exemple calculée en fonction de l'indice de constance de phase démodulée. Ainsi, une composante stationnaire peut être très colorée et une composante peu stationnaire peut être d'autant plus décolorée et/ou blanchie que l'indice de constance de phase démodulée est faible. Les zones colorées peuvent ainsi isoler les composantes sinusoïdales des autres types de composantes, correspondant par exemple à des bruits ou des impacts. La constance de la couleur peut ainsi par exemple correspondre à une phase démodulée constante. La constance d'une couleur associée à une fréquence cible d'accordage peut ainsi permettre d'identifier la perfection d'un accord sur cette fréquence, par exemple au-delà du seuil de détectabilité de la justesse de l'oreille humaine, et de garantir ainsi, par exemple, que deux instruments accordés sur cette base ne produiront pas de battement s'ils jouent à l'unisson. Par exemple, une variation de couleur cyclique sur cinq secondes peut indiquer que des battements cycliques de cinq secondes sont susceptibles se produire. Alternativement, la teinte peut être exclusivement contrôlée par le module normalisé. Une telle coloration ne représente plus d'information de phase mais surcharge l'information de représentation de l'amplitude. Une telle surcharge permet de simplifier la lecture sur l'image 300 de situations acoustiques, par exemple sonores, complexes et d'alléger l'effort d'analyse visuelle. La qualité de la stationnarité qui raffine la précision fréquentielle visualisée reste encore prise en charge par l'épaisseur, par exemple l'épaisseur et/ou par la saturation de la couleur de la forme 302, par exemple du trait de la spirale 306. L'étape 122 de réalisation de l'image 300 comprend par exemple une étape 210 de calcul de transparence de la forme 302, par exemple du trait, par exemple pour chaque fréquence du vecteur de fréquence. La transparence peut être calculée de sorte que si la forme 302, par exemple le trait, devient très fine, il n'en résulte pas un écrasement de l'abaque en spirale 306. En effet, il est ainsi possible d'éviter que des évènements quasi-inaudibles, tels qu'un faible bruit, construisent un trait par exemple peu épais et par exemple peu brillant, c'est-à-dire presque noir, qui viendrait masquer le trait de l'abaque, et ainsi d'éviter l'effet de scintillement du trait de l'abaque qui en résulte et a un effet visuellement perturbant. Il est ainsi possible que la grille 304 reste visible lors de l'affichage de l'image 300, sans effet de scintillement pour les éléments visuels 302 correspondant à des fréquences à faible pertinence tels que des fréquences correspondant à des modules faibles ou très instationnaires. La transparence est par exemple une fonction croissante de l'épaisseur de la forme 302, par exemple du trait, par exemple linéairement croissante, par exemple entre une valeur minimale, par exemple 0, et une valeur maximale, par exemple 1. La transparence peut être saturée à la valeur minimale en-dessous d'une certaine valeur d'épaisseur, et/ou saturée à la valeur maximale au-dessus d'une certaine valeur d'épaisseur. Représentations graphiques L'étape 122 de réalisation d'image comprend par exemple la construction 212 et/ou 214 de représentations graphiques à afficher. En particulier, l'étape 122 de réalisation d'image peut comprendre par exemple le tracé 2120 de la grille 304, par exemple de la grille 304 de la spirale 306. L'étape 122 de réalisation d'image peut comprendre le tracé 2122 de la spirale colorée formée des éléments visuels et l'affichage 2124 de la spirale 306 et de la grille 304. Par exemple, pour des données acoustiques comprenant une plage de notes à afficher, par exemple encodées via des codes midi compris entre une valeur minimale et une valeur maximale, on construit une grille 304 dont le rayon est par exemple calculé en coordonnées polaires par une formule du type : p(0) = 1 + (0 - Omin) où 0 est l'angle en coordonnées polaires, p le rayon en coordonnées polaires, et la valeur minimale d'angle, par exemple en prenant la valeur minimale en code MIDI, et dont l'angle est par exemple calculé en degré par une formule du type : 0 = 360 (Midicode - 69) 12 _L uref O (nnidi_code)=360. (nnidi_code-69)/12+A4ReferenceAngle où nnidicode représente le code midi de la note considérée, et où Oref représente l'angle de référence. Ainsi, la construction d'une telle grille 304 revient par exemple à construire une spirale 306 en effectuant un tour par octave de sorte que toutes les notes en relation d'octave sont alignées et sont associées au même angle sur la spirale. Dans le cas d'un affichage sur un écran tactile, l'angle peut par exemple être modifié en posant un ou plusieurs doigts sur la représentation de la spirale 306 sur l'écran, et en effectuant un mouvement de rotation avec le ou les doigts. La grille 304 peut par exemple représenter la position des notes tempérées par des lignes appuyées et/ou des subdivisions par des lignes moins appuyées. L'étape 122 de réalisation d'image peut comprendre par exemple une étape 2124 successive de superposition sur la grille 304 de la spirale 306. La superposition 2124 est par exemple réalisée avec le niveau de transparence local réglé par la transparence calculée au préalable. Une telle superposition 2124 est par exemple réalisée avec une interpolation adaptée aux nombres de pixels disponibles, par construction vectorielle de l'image 300 colorée. Une telle représentation permet de fournir une représentation limitant la perte des informations obtenues aux étapes précédentes. Une telle représentation peut ainsi être compatible avec un usage professionnel en production, par exemple dans le cadre d'activités d'ingénieur du son, de nnastering, ou de restauration d'enregistrements anciens. Ainsi, plutôt que sur un axe, le spectre peut être représenté sur un squelette en spirale qui assure l'alignement des fréquences si un alignement est souhaitable. Par exemple, dans le cas de données sonores, les notes peuvent être alignées en relation d'octave. Il est ainsi possible de réaliser une topologie cohérente avec la notion de note musicale sur une échelle tempérée. Dans le cadre d'une utilisation en tant qu'accordeur un tel procédé, il est possible de ne recourir à aucune estimation de hauteur de note, estimation qui est souvent mal définie pour les sons réels. Le procédé permet de se passer d'estimation en permettant à l'utilisateur, par sa perception visuelle de l'image 300 et en particulier de l'élément visuel 302, d'être lui-même l'accordeur. Si un son est bien harmonique, accorder se résumera à aligner correctement une composante fondamentale sur une grille 304 de la note visée, et éventuellement les harmoniques en relation d'octaves. Si l'harnnonicité n'est pas idéale, la notion de hauteur ne coïncide plus avec la fréquence fondamentale et devient mal posée de sorte que les accordeurs standards donnent des estimations dépendantes des méthodes choisies par les constructeurs. En effet, dans un tel cas, le signal n'est pas un signal idéalement périodique. Dans un tel cas, les composantes fréquentielles du signal acoustique, par exemple du son, ne suivent par exemple pas exactement une progression harmonique avec une le" harmonique ayant une fréquence fk multiple d'une fréquence fondamentale fl, par exemple fk=k*fi. Par exemple des cordes (spécialement pour les notes graves de piano), les barres (vibraphone, nnarinnba, etc) ou des cloches ne suivent pas de progression harmonique. Le procédé et le système 400 décrits peuvent permettre à l'utilisateur de visualiser le résultat spectral pour chaque composante, et d'apprécier lui- même le compromis à ajuster. Un tel procédé et un tel système 400 peuvent donc s'avérer particulièrement utiles dans le cas de relations de quasi-octaves, où les composantes ne peuvent plus être alignées exactement avec le centre de la spirale. Le procédé décrit peut donc permettre de comprendre la source de la difficulté de l'accordage dans un tel cas, source qui est complètement masquée pour un accordeur standard, et d'y remédier par les retours auditifs et visuels à partir desquels le choix par l'utilisateur d'un compromis sera guidé par sa propre expérience psycho-acoustique plutôt qu'une estimation au sens d'un modèle idéalisant. En complément, le procédé peut comprendre en réponse à la sélection d'une note par un utilisateur, une étape de représentation agrandie et redressée de la portion de la spirale 306 pour cette fréquence, par exemple pour cette note, et par exemple également des zones voisines, par exemple des zones voisines sur une octave. Une telle fonctionnalité est particulièrement adaptée à la mise en oeuvre d'un système 400 accordeur. Dans le cas d'un affichage sur un écran tactile, en complément, le procédé peut comprendre, en réponse à la sélection d'une note par un utilisateur, une étape de génération d'un son qui correspond à cette note accordée sur le diapason réglé, par exemple sinusoïde pure de la fréquence de la note pointée. Alternativement ou en complément, la construction de représentations graphiques à afficher comprend par exemple le tracé 214 d'une image colorée de type - sonagramme », et par exemple son affichage. Un tel tracé 214, utilisé en complément de la représentation en spirale 306, permet une préservation de la mémoire des images déjà affichées, en représentant un historique de ces images 306. L'abscisse du sonagramme correspond par exemple au temps t' de chaque trame. L'ordonnée du sonagramme correspond par exemple aux fréquences du vecteur de fréquence, par exemple aux notes des fréquences analysées. La couleur de chaque point est par exemple calculée comme indiqué précédemment. De même, un tel affichage peut être réalisé avec une interpolation adaptée aux nombres de pixels disponibles, par construction vectorielle de l'image colorée. Un tel procédé peut être mis en oeuvre dans le cadre de différents modes de fonctionnement possibles. Par exemple, lorsque les données sont des données sonores, des modes de fonctionnement possibles sont par exemple chants d'enfants, chorale, jeu d'instruments, accordeur, post-production et analyse du signal musical et sonore. En mode accordeur, l'accord est par exemple effectué par le musicien en modifiant son instrument pour centrer sur la grille 304 accordée le percept visuel associé à la composante spectrale détectée. Si besoin pour des configurations particulières, le procédé peut permettre d'affiner encore plus la qualité d'accordage, par exemple au-delà des capacités de l'oreille sur un son nnonodique, de sorte à assurer l'absence de - battements » avec d'autres instruments mais sans nécessiter la présence d'un second instrument. Pour cela, l'instrumentiste doit par exemple assurer que la couleur, si elle est associée à la phase démodulée, ou bien ne varie pas (accordage idéal), ou bien dérive très lentement (la période d'un - tour de couleur » correspondant par exemple à la période du battement avec un autre instrument parfaitement accordé avec une phase démodulée idéalement constante. Traitement des harmoniques En référence à la figure 5a, alternativement ou en complément des étapes décrites ci-dessus, le procédé peut comprendre un traitement d'harmoniques, un tel traitement pouvant être suivi par la construction d'une représentation graphique à partir d'un tel traitement. Etape de répartition Ainsi, en référence à la figure 5b, le procédé peut comprendre une étape 502 de répartition d'une partie d'une énergie d'une fréquence f traitée 30 vers des énergies de fréquences sous-harmoniques de la fréquence traitée, par exemple selon des proportions déterminées. Une telle étape 502 de répartition peut être réalisée pour plusieurs fréquences f d'une plage fréquentielle d'analyse donnée, par exemple fixée par 25 les bornes du vecteur de fréquence. Une telle étape de répartition peut être réalisée pour chaque trame donnée. Une telle étape 502 de répartition peut être réalisée itérativement, par exemple itérativement par octave, par exemple itérativement par octave décroissante. L'étape 502 de répartition peut par exemple être mise en oeuvre par le système 400 de traitement de données acoustiques. La répartition peut être déterminée au prorata d'une fonction d'énergie de sous-harmoniques ramenée à l'énergie de la fréquence f testée. Une telle 10 fonction est par exemple une fonction monotone, par exemple une fonction croissante. Alternativement ou en complément, la répartition peut être déterminée en fonction des indices de constance de phase de sous-harmoniques et de la fréquence f testée, par exemple avec une activation qui est fonction des 15 indices de constance de phase de sous-harmoniques et de la fréquence f testée, la fonction étant par exemple une fonction monotone des indices de constance de phase, par exemple une fonction croissante. Il est ainsi possible de sélectionner les harmoniques susceptibles de participer à un timbre perçu, par rapport à d'autres signaux tels que des bruits. 20 L'étape de répartition 502 peut ainsi comprendre, par exemple pour chaque trame, une étape 5020 d'obtention de contributions harmoniques, les contributions harmoniques permettant par exemple de quantifier une marque de présence de contribution harmonique associée à une hauteur de structure quasi-périodique telle qu'une note. 25 Une telle répartition permet de quantifier la présence de notes potentielles parmi les harmoniques, les notes potentielles étant celles qui ont reçu le plus d'énergie à l'issue de l'étape de répartition. Représentation graphique de contributions harmoniques 30 Les contributions harmoniques peuvent être mises en oeuvre dans la construction d'une représentation graphique. Les contributions harmoniques sont par exemple obtenues sous la forme d'amplitudes de contribution harmonique. De telles amplitudes sont ainsi aisées à intégrer dans une représentation graphique. Le procédé peut ainsi comprendre une étape 504 subséquente de construction d'une image dont un élément visuel au moins est construit en fonction des contributions harmoniques, par exemple des amplitudes de contribution harmonique. Le procédé peut également comprendre une étape 506 subséquente d'affichage d'une telle image, par exemple par le système 400, par exemple par un afficheur. En particulier, le procédé peut ainsi permettre de mettre en oeuvre une extension en trois dimensions des images décrites ci-avant réalisées par exemple en deux dimensions. L'extension en trois dimensions peut par exemple être mise en oeuvre de sorte qu'a harmonique h fixée, au moins une partie des caractéristiques visuelles déjà mises en oeuvre sont réutilisées, par exemple la totalité de ces caractéristiques visuelles, les caractéristiques visuelles réutilisées comprenant par exemple attributs de couleur et d'épaisseur de forme, par exemple de trait, telles que définies précédemment. L'extension en trois dimensions peut par exemple représenter plusieurs images en deux dimensions de données associées à différentes valeurs d'harmoniques. L'extension en trois dimensions peut ainsi par exemple prendre en compte une troisième coordonnée spatiale z. La troisième coordonnée spatiale est par exemple utilisée pour représenter les images en deux dimensions des données associées à différentes valeurs d'harmoniques. Les images en deux dimensions sont par exemple aux altitudes z=log2(h), en particulier pour des valeurs telles que0, 1, log2(3), 2, log2(5), etc. . Le procédé peut ainsi comprendre la construction d'un vecteur qui quantifie la présence de notes à des fréquences déterminées. Un tel vecteur peut être utilisé en tant qu'amplitude dans l'afficheur en deux dimensions décrit précédemment, de sorte que l'image indique la présence de notes. L'étape 504 de construction d'image peut comprendre l'application 5040 d'une interpolation peut être appliquée, par exemple entre les altitudes décrites, par exemple pour au moins une caractéristique visuelle telle que l'épaisseur de la forme, par exemple l'épaisseur du trait, la couleur et/ou la transparence. L'interpolation 5040 est par exemple constante ou affine par morceaux. Il est ainsi possible d'obtenir visuellement des structures en "chaînes de montagnes" représentant le nombre et la présence des harmoniques. En particulier, la hauteur totale et la largeur en fonction de l'altitude peuvent être fonction respectivement du nombre et de la présence des harmoniques. De telles "chaînes de montagnes" peuvent ainsi par exemple évoluer avec le timbre. Il est ainsi possible, sans établir d'estimations de notes, de fournir par une image en trois dimensions une aide à la localisation visuelle de notes avec une caractérisation du timbre, tel qu'un rehaussement visuel pour l'indication de la présence de notes. Une telle image en trois dimensions peut ainsi constituer un indicateur élaboré, plus adapté à la mise en relation avec la perception de la musique. De plus, il est ainsi possible de répartir d'une façon équilibrée l'énergie de chaque fréquence sur un ensemble de contributions harmoniques, par exemple vers des fréquences ciblées et déterminées par leur activation énergétique et leur qualité de stationnarité. De même, l'inharnnonicité n'est pas traitée au sens d'un gabarit ou d'un modèle paramétrique de déviation. L'inharnnonicité est toutefois prise en charge par le procédé dans la mesure où la fréquence quasi-harmonique traitée divisée par un index harmonique h est localisé dans un lobe principal d'une fréquence d'une note testée par le procédé. Un exemple de mode de réalisation détaillé d'un tel traitement des harmoniques est détaillé ci-après. Tableau de phases démodulées Le procédé peut comprendre, en amont de l'étape 502 de répartition, une étape 500 de réception de données, par exemple comprenant la réception 5000 de la phase démodulée, par exemple d'une pluralité de phases démodulées, par exemple d'un vecteur de données correspondant à la phase démodulée pour une ou plusieurs trames. Le procédé peut comprendre une étape de réception ou de calcul 5002 d'un tableau de phases démodulées comprenant la phase démodulée d'au moins une fréquence du vecteur de fréquences et d'au moins une harmonique de ladite fréquence du vecteur de fréquence. La hème harmonique fk,h, h étant par exemple un entier naturel, de la le" fréquence fk, k étant par exemple un entier naturel, du vecteur de fréquences est par exemple déterminée par une formule du type : fk,h = fk x h En particulier, le procédé peut comprendre une telle étape 5002, de sorte à ce que le tableau de phases démodulées comprennent de telles phases 10 démodulées pour toutes les fréquences du vecteur de fréquence et une pluralité d'harmoniques pour chaque fréquence du vecteur de fréquence. Une telle étape 500 peut par exemple comprendre le calcul de la phase démodulée pour la ou les harmoniques de fréquence du vecteur de fréquence par interpolation. L'interpolation est par exemple une interpolation - au plus 15 proche voisin » ou affine par morceau. L'encodage peut être un encodage entier, par exemple sur 32 bits, pour garantir la précision et le modulo décrits précédemment. Les fréquences fk,h disposées au-delà des valeurs extrémales du vecteur de fréquences peuvent ne pas être prises en compte. 20 Tableau d'indices de phase démodulée De même, le procédé peut comprendre une étape de réception 5004 de l'indice de phase démodulée, par exemple d'une pluralité d'indices de phase 25 démodulée, par exemple d'un vecteur de données correspondant à l'indice de phase démodulée pour une ou plusieurs trames. Le procédé peut comprendre une étape de réception ou de calcul 5006 d'un tableau d'indices de phases démodulées comprenant la phase démodulée d'au moins une fréquence du vecteur de fréquences et d'au moins une 30 harmonique de ladite fréquence du vecteur de fréquence. En particulier, le procédé peut comprendre une telle étape 5006, de sorte à ce que le tableau d'indices de phase démodulée comprenne de tels indices pour toutes les fréquences du vecteur de fréquence et une pluralité d'harmoniques pour chaque fréquence du vecteur de fréquence. Une telle étape 500 peut par exemple comprendre le calcul de l'indice de phase démodulée pour la ou les harmoniques de fréquence du vecteur de fréquence par interpolation. L'interpolation est par exemple une interpolation - au plus proche voisin » ou affine par morceau. L'encodage peut être un encodage entier, par exemple sur 32 bits, pour garantir la précision et le modulo décrits précédemment. Les fréquences fk,h disposées au-delà des valeurs extrémales du vecteur 10 de fréquences peuvent ne pas être prises en compte. Tableau de contributions harmoniques La contribution harmonique est par exemple une amplitude de 15 contribution harmonique. De même, le procédé peut comprendre une étape de réception 5018 de contributions harmoniques, par exemple d'une pluralité de contributions harmoniques, par exemple d'un vecteur de données correspondant à des contributions harmoniques pour une ou plusieurs trannes.Le procédé peut 20 comprendre une étape de réception ou de calcul 5020 d'un tableau de contributions harmoniques comprenant la contribution harmonique associée à au moins une fréquence du vecteur de fréquences et d'au moins une harmonique de ladite fréquence du vecteur de fréquence. Il est ainsi possible d'identifier les fréquences d'origine de l'énergie 25 reçue par une harmonique. Une telle étape permet ainsi de construire une représentation visuelle caractéristique du timbre. En particulier, le procédé peut comprendre une telle étape 5020, de sorte à ce que le tableau de contributions harmoniques comprenne de telles 30 amplitudes pour toutes les fréquences du vecteur de fréquence et une pluralité d'harmoniques pour chaque fréquence du vecteur de fréquence. L'étape de calcul 5020 d'un tableau d'amplitudes est par exemple une étape de calcul par itération. L'étape de calcul 5020 d'un tableau de contributions harmoniques comprend par exemple une étape d'initialisation 50200. L'étape d'initialisation 50200 comprend par exemple une étape 502000 de réception d'un vecteur d'énergies des fréquences du vecteur de fréquence. Les énergies sont par exemples des énergies normalisées. Ainsi, l'étape d'initialisation 50200 peut par exemple comprendre une étape 502002 de normalisation des énergies. La normalisation 501002 consiste par exemple à rendre la somme des énergies égale à une valeur donnée, par exemple à une valeur d'unité. La normalisation 502002 est par exemple obtenue par une formule du type : E(f) Enorm(fk) = '-tot OÙ Enorm(fk) est l'énergie normalisée pour la kème fréquence du vecteur de fréquence, où E(fk) est l'énergie pour la la kème fréquence du vecteur de fréquence, et où Et, est l'énergie totale. L'énergie d'une fréquence, est par exemple l'énergie stockée dans une bande fréquentielle associée à cette fréquence, par exemple une bande fréquentielle stockée sur cette fréquence. Les bandes fréquentielles sont par exemple des bandes connexes, par exemple formant un recouvrement d'une plage de fréquences donnée. Les fréquences du vecteur de fréquences suivent par exemple une progression en suite, par exemple en suite géométrique. Chaque fréquence peut être associée à une bande de fréquence, de sorte que les bandes de fréquences correspondant à des fréquences qui se suivent sont connexes. Chaque bande de fréquences est par exemple centrée sur la fréquence 25 correspondante, par exemple centrée au sens d'une moyenne géométrique. La largeur Lf d'une bande associée à une fréquence f est par exemple fournie par une formule du type : L =Kxf avec K un coefficient indépendant de la fréquence, dont la valeur est par exemple fournie par une formule du type : K = 2x sinh(ln(2) x RmIDI/24) où RmiDi est par exemple une valeur de pas fractionnaire entre deux valeurs de code MIDI correspondant à deux fréquences consécutives du vecteur de fréquence. L'énergie E(fk) dans la bande de largeur Lfk associée à la fréquence donnée fk peut ainsi par exemple être fournie par une formule du type : E(fk) =Kxfkx AÏ, OÙ Ak est l'amplitude du signal pour la fréquence fk. L'énergie totale Et, peut alors être définie comme la somme des énergies pour toutes les fréquences fk. L'étape d'initialisation 50200 comprend par exemple en outre une étape 502004 de construction d'un tableau initial des contributions harmoniques. L'étape 502004 correspond par exemple à la construction d'un tableau initial des amplitudes de contribution harmonique AHarmo,k(h), par exemple de sorte que AHarmo,k(1) = Ak, et AHarmo,k(h)=0 pour des harmoniques d'index h compris entre 2 et Hmax, où Hmax est par exemple un vecteur des indices harmoniques maximaux associés aux fréquences du vecteur de fréquence. Le vecteur des indices harmoniques maximaux est par exemple défini par : Hmax(k) = Floor(' 7X) x) lk où Floor est la partie entière, par exemple arrondie par défaut, Fma, la borne supérieure de la plage de fréquences traitées. La plage de fréquences traitées est par exemple du type [F,,,,,Fmax] avec Fm,,, la valeur minimale parmi les fréquences du vecteur de fréquences et Fmax la valeur maximale parmi les fréquences du vecteur de fréquences. La plage de fréquences traitées comprend par exemple une pluralité d'octaves pleines, le nombre d'octaves pleines Noctaves étant par exemple déterminé par une formule du type : Noctaves = F100r(10g2 (Fmax/Fmin)) où Floor définit par exemple la partie entière, par exemple la partie entière par défaut. Le procédé peut alors comprendre la réception ou la construction d'une séquence décroissante de fréquences de plage Fimage, par exemple pour n entier naturel allant de 1 à Noctave' par exemple par une formule du type Fplages(n) 21-11 X Fniax par exemple de sorte que les octaves traitées par l'étape de ré-assignation 50204 des énergies de bandes fréquentielles décrite ci-après correspondent aux plages fréquentielles Plage(n), par exemple ouvertes à gauche et fermées à droite, du type : Plage(n) = 1F plages(n + 1) ;Fplages(n)1 L'étape de calcul 5020 par itération comprend par exemple une étape 50102 de calcul de coefficients de transfert pour chaque fréquence du vecteur de fréquence appartenant la plage de fréquences traitées, par exemple de calcul d'un vecteur de tels coefficients. L'étape 50202 de calcul de coefficients de transfert est par exemple réalisée au prorata d'une fonction croissante de l'énergie E(fk/sh), où sh est un indice sous-harmonique' ladite énergie étant ramenée à l'énergie E(fk) de la fréquence fk testée. sh est ainsi un indice sous harmonique correspondant à une fréquence réceptrice d'énergie pour l'étape de ré-assignation, l'énergie étant pompée depuis la fréquence f testée, par exemple localisée par construction hors de l'octave analysée. Une telle étape est par exemple réalisée pour des indices sous-harmoniques 5h=2,3,..,SHmax(k), où SHmax(k) est le le" coefficient d'un vecteur des indices sous-harmoniques maximaux associés aux fréquences du vecteur de fréquences, ledit coefficient étant par exemple défini par une formule du type : Sliniax(k) = Floor(fk/Fmin) L'étape 50202 de calcul de coefficients de transfert est par exemple réalisée avec une activation qui est fonction croissante des indices de constance de phase de l' émetteur d'énergie, c'est-à-dire de la fréquence testée, et du récepteur d'énergie c'est-à-dire des indices sous-harmoniques correspondants. Une telle étape 50202 de calcul de coefficients de transfert est par exemple réalisée de sorte à pomper l'énergie Enorm(fk) et n'en laisser qu'une proportion Coeff(1). En outre, une telle étape 50202 de calcul de coefficients de transfert est par exemple réalisée de sorte à répartir le reste de l'énergie Enorm(fk) vers les sous-harmoniques d'indices sh, compris entre 2 et SHmax(k), par exemple en transférant les proportions Coeff(sh) vers les fréquences d'indice KSH(k,sh). Par construction de Plage(n), toutes les fréqences réceptrices d'énergie sont inférieures à Plage(n). L'entier KSH(k,sh) est par exemple défini de sorte que fKsH(k,sh) est la valeur du vecteur de fréquences la plus proche de la fréquence sous- harmonique fk/sh. L'étape 500 comprend par exemple une étape 5008 de réception ou de construction d'un tableau des entiers KSH(k,sh). Puisque la valeur maximale de sh dépend de k, un tel tableau n'est pas nécessairement rectangulaire et seule la partie complétée est utile. De ce fait, un encodage par liste peut être choisi pour un stockage en mémoire plus compact. Le vecteur de coefficients de transfert positifs, et par exemple de somme 1, est donné par une formule du type : Coe f f (sh) = C(sh)/C0t où Ctot est la somme des composantes C(sh) et où, C(sh) est par exemple défini comme suit : C(1) = 1 et pour sh compris entre 2 et SHmax(k) : C(sh) = (cl)const(k) x cl)const(k, sh))a/2 X Enorm(K S H (k, sh) / ( / Enorm(k))b 1 2 x((econst (k)xeconst(k,sh))72où Oconst(k) est l'indice de constance de phase démodulée, (D - const(k,sh) est par exemple défini comme (Dconst(fKSH(k,sh)), a, b c sont des coefficients, par exemple compris entre 0 et 100, par exemple tels que a=0,1, b=10 et c=0,1. Les coefficients a, b et c sont par exemple modifiables. Le choix a=c=0 rend par exemple l'attribution insensible à la constance de phase. Ré-assignation des énergies et amplitudes Le procédé peut par exemple comprendre une étape 50204 de réassignation des énergies et amplitudes. Le procédé peut ainsi comprendre une étape 502040 de mise à jour du vecteur d'énergie, par exemple du vecteur d'énergie normalisée et du tableau d'amplitudes de contributions harmoniques. Une telle étape 502040 permet d'effectuer le transfert des amplitudes en garantissant la préservation de l'énergie. L'étape 502040 de mise à jour comprend par exemple une étape de transfert des amplitudes de contributions harmoniques vers les sous-harmoniques ("récepteur" d'énergie). Par exemple, l'étape de transfert peut comprendre, pour chaque indice sous-harmonique sh allant de 2 à SHmax (k), l'application de formules du type: ksh = KSH(k,sh) où ksh est l'indice de la fréquence sous-harmonique traitée, Enorm (ksh) = Enorm(kSh) + Coeff (sh) x Enorm(k) h=sh où h est l'indice harmonique de base traité pour la fréquence sous-harmonique traitée. L'étape de transfert peut en outre comprendre, pour hbis allant de 1 à min( Hmax (k) , Floor(Fmax/(F(ksh) x h)) ), l'application d'une formule du type : AHarmo,ksh(h X hbis) = (AHarmo,ksh(h x hb15)2 + Coeff(sh) x AHarmo,k(hms)2)0'5 où hbis correspond aux indices harmoniques de la fréquence fk déjà traités et que l'on doit aussi tous transférer, et où la puissance 0.5 dans la formule provient du fait que l'on traite des amplitudes et non des énergies, une telle puissance étant est nécessaire pour garantir la préservation d'énergie totale.] L'étape 502040 de mise à jour peut comprendre en outre une étape dite de pompage, par exemple de pompage sur l' "émetteur" d'énergie. L'étape de pompage comprend par exemple l'application de formules du type : Enorm(k) = Coeff(1) x Enorm(k) et pour hm, allant de 1 à HmaAk) : AHarmo,k(hbis) = Coeff (1)0'5 x AHarmo,k(hbiS) Le procédé comprend par exemple l'itération de l'étape 502040 de mise à jour pour chaque indice k de sorte que fk appartient à Plage(n). Construction de vecteur d'amplitude Le procédé peut en outre comprendre une étape 508 de construction du vecteur d'amplitudes caractérisant la présence de notes sans information de timbre. Le vecteur AnnpNote associé aux fréquences f est par exemple donné par: AnnpNote(k) = (Enorm(k) x Et0til-fk)°.5 Système En référence à la figure 4, il est décrit un système 400 de traitement de données acoustiques. Le système 400 est par exemple un système électronique de traitement et de stockage de données, par exemple de type serveur ou poste de travail 15 informatique. Le système 400 est adapté pour mettre en oeuvre un procédé tel que décrit précédemment. Le système 400 comprend un module 402 de réception de données acoustiques correspondant à un signal acoustique transformé dans le domaine 20 fréquentiel à partir d'un signal acoustique dans le domaine temporel enregistré par au moins un capteur acoustique, Le système 400 comprend un module 404 de démodulation de la phase du signal acoustique transformé pour au moins une fréquence du domaine fréquentiel, et 25 Le système 400 comprend un module 406 de calcul d'un indice de constance de la phase démodulée. Le système peut comprendre un module (non représenté) de rehaussement visuel de présence de notes. Le module de rehaussement visuel peut comprendre un indicateur de timbre. 30 On comprendra que selon un mode de réalisation, le système comprend un moyen de traitement de données de type processeur adapté pour mettre en oeuvre les modules 402, 404, 406 précités qui peuvent donc être des modules logiciels. Alternativement, on comprendra que les modules 402, 404 et 406 peuvent être des modules physiques. Le système 400 peut en outre comprendre des moyens de stockage de données, tels qu'un disque dur, pour le stockage des données acoustiques et/ou des résultats. Le système 400 comprend en outre une interface, par exemple un afficheur, par exemple un écran permettant d'afficher des graphiques et/ou des images résultant du procédé décrit ci-dessus, de sorte à restituer les données associées à destination d'un utilisateur. Produit programme d'ordinateur Selon un autre aspect, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé de traitement tel que décrit précédemment. Moyens de stockage Selon encore un autre aspect, l'invention concerne des moyens de stockage lisibles par un équipement informatique sur lesquels un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé de traitement tel que décrit précédemment. REVENDICATIONS1. Procédé de traitement de données acoustiques, comprenant des étapes de: réception (110) de données acoustiques correspondant à un signal acoustique transformé dans le domaine fréquentiel à partir d'un signal acoustique dans le domaine temporel enregistré par au moins un capteur acoustique, pour au moins une fréquence du domaine fréquentiel, démodulation (118) de la phase du signal acoustique transformé, et calcul (120) d'un indice de constance de la phase démodulée. 2. Procédé selon la revendication 1, comprenant en outre des étapes de : réception (112) d'un vecteur de fréquences, calcul (114) par interpolation des valeurs du signal acoustique transformé aux fréquences du vecteur de fréquences, de sorte que la démodulation (118) est réalisée pour chaque fréquence du vecteur de fréquences. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel la phase démodulée est obtenue par une formule du type : (1)demod(f, , t) - I`I 1) spectr (f , t) - 271- f t OÙ Odemod représente la phase démodulée pour la fréquence f et au temps t, et OÙ Ospectr représente la phase du signal transformé pour la pour la fréquence f et au temps t. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le calcul (120) d'un indice de constance de la phase démodulée comprend l'application (1202) d'un filtre passe-bas, le filtre passe-bas étant adapté pour réaliser une atténuation en fonction de la vitesse de rotation de la phase démodulée. 5. Procédé selon la revendication 4, dans lequel le calcul d'un indice de constance de la phase démodulée comprend, préalablement à l'application (1202) du filtre passe-bas la conversion (1200) de la phase démodulée en une exponentielle complexe de module constant, le filtre passe-bas étant appliqué à l'exponentielle complexe, l'indice de constance de la phase démodulée résultant du module d'une sortie complexe du filtre passe-bas. 6. Procédé selon l'une quelconque des revendications 1 à 5, comprenant en outre des étapes de : réception d'une valeur de variation de fréquence à détecter, réception (104) de données acoustiques correspondant au signal acoustique dans le domaine temporel enregistré par l'au moins un capteur acoustique, transformation (108) du signal acoustique dans le domaine temporel en le signal dans le domaine fréquentiel par transformée de Fourier discrète pondérée par une fenêtre d'analyse, la transformée étant réalisée avec un pas d'avancement déterminé de sorte que l'inverse du double du pas d'avancement est supérieur à : la largeur du lobe principal de la fenêtre d'analyse dans le domaine de Fourier, et à la valeur de variation de fréquence à détecter. 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel le signal acoustique enregistré est un signal sonore et dans lequel le signal acoustique dans le domaine temporel enregistré est échantillonné. 8. Procédé selon l'une quelconque des revendications 1 à 7, comprenant en outre une étape (122) de construction d'une image (300) représentative de la constance de phase démodulée issue des données acoustiques, l'image étant adaptée pour être affichée sur un terminal, l'image (300) étant construite à partir des données acoustiques et de l'indice de constance de phase démodulée. 9. Procédé selon la revendication 8, dans lequel la construction (122) de l'image comprend un élément visuel (302) correspondant à la fréquence du domaine fréquentiel, la représentation de l'élément visuel (302) étant fonction de l'amplitude de la fréquence et de l'indice de constance de phase démodulée. 10. Système de traitement de données acoustiques, adapté pour mettre en oeuvre un procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce qu'il comprend : un module (402) de réception de données acoustiques correspondant à un signal acoustique transformé dans le domaine fréquentiel à partir d'un signal acoustique dans le domaine temporel enregistré par au moins un capteur acoustique, un module (404) de démodulation de la phase du signal acoustique transformé pour au moins une fréquence du domaine fréquentiel, et un module (406) de calcul d'un indice de constance de la phase démodulée.