Brevet FR3060753A1 — publié le 20180622

La présente invention concerne un procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore pour réaliser de l'imagerie ou de la focalisation adaptative et sélective. Elle concerne également un programme d'ordinateur et un dispositif de sondage à ultrasons correspondants. L'invention s'applique en particulier au domaine du contrôle non destructif par ultrasons, dans lequel l'acquisition de signaux ultrasonores permet de détecter et/ou visualiser des défauts dans des structures, mais elle peut aussi s'appliquer à tout type de détection ou imagerie échographique ultrasonore, notamment au domaine médical pour l'inspection de zones d'intérêt dans le corps humain ou animal. Elle porte plus particulièrement sur un procédé de traitement acquérant les signaux ultrasonores de la façon suivante : - commande de L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission, - commande de N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, - obtention d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR, J(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission. Une telle acquisition se fait généralement à l'aide d'un dispositif de sondage à capteur multiéléments, dans lequel chaque transducteur est à la fois émetteur et récepteur, une commutation entre ces deux modes pouvant être commandée électroniquement. Le capteur peut être mis au contact de l'objet à sonder ou à distance, mais dans ce dernier cas il doit être immergé pour assurer la transmission des ondes ultrasonores dans l'objet à sonder. Ce capteur peut être linéaire (1D) ou matriciel (2D), à éléments rigides ou flexibles. La matrice [MR(t)] de signaux temporels obtenue par ce type d'acquisition, généralement qualifiée de matrice des ondes planes, peut alors faire l'objet d'un traitement, notamment pour la fourniture d'une image de la zone d'intérêt inspectée ou pour l'extraction de paramètres significatifs de défauts de structure dans la zone d'intérêt inspectée. Compte tenu des capacités de calcul actuelles des processeurs, 3060753 2 ce traitement peut être embarqué dans les instruments de contrôle pour des traitements en temps réel. Ce type d'acquisition, généralement qualifié de « plane-wave compounding » ou « plane-wave imaging », est par exemple décrit dans l'article de Montaldo et al, 5 intitulé « Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography », publié dans IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 56, n° 3, pages 489-506, mars 2009. Il est essentiellement envisagé dans le domaine médical et mis en oeuvre dans certains échographes ultra-rapides permettant d'imager des cartographies de l'élasticité du 10 corps humain. L'image résultante est de haute qualité et ne nécessite que quelques dizaines de tirs ultrasonores (généralement M est compris entre 10 et 30) pour un capteur de N = 128 transducteurs d'émission/réception, contre 128 tirs ou presque pour d'autres techniques d'imagerie telles que la focalisation synthétique en tous points basée sur une acquisition de type FMC (de l'anglais « Full Matrix Capture ») 15 ou équivalent. De plus, les algorithmes résultants et les applications associées se prêtent particulièrement bien à une parallélisation des calculs sur des processeurs de type GPU (de l'anglais « Graphic Processing Unit ») implantés dans les cartes graphiques des ordinateurs. Les performances d'un échographe implémentant une technique de « plane-wave compounding » peuvent ainsi atteindre en pratique 20 10 000 images/s. Un autre atout de cette technique d'acquisition réside dans le fait que chaque tir est réalisé en sollicitant tous les transducteurs d'émission de sorte que l'énergie émise est élevée, rendant cette méthode moins sensible aux phénomènes d'atténuation, de bruit électronique ou de structure. Ce type d'acquisition est également exploité dans la demande de brevet WO 25 2015/092250 Al en lui adaptant astucieusement le principe de focalisation synthétique en tous points, de manière à tirer profit de la simplicité de la technique de « plane-wave compounding » en vue d'atteindre une cadence d'acquisition élevée et de la qualité d'image, en termes de résolution spatiale et de contraste, liée à une focalisation synthétisée en tous points de l'image souhaitée. Cela a permis 30 d'envisager la technique de « plane-wave compounding » pour des applications de contrôle non destructif. Mais quelle que soit la technique d'imagerie choisie, les images obtenues en contrôle non destructif peuvent présenter un niveau de bruit significatif en fonction des propriétés du milieu sondé. Par exemple, il s'agit d'un bruit électronique lorsque le matériau est homogène et viscoélastique, ou d'un bruit de structure lorsque les 3060753 3 ondes sont diffusées par des hétérogénéités du matériau. Les images calculées avec les techniques de focalisation en tous points sont aussi affectées par des artefacts, ou échos fantômes, liés à des échos de géométrie, par exemple l'écho d'un fond de pièce à proximité de celui d'un défaut. 5 Il peut ainsi être souhaité de concevoir un procédé de traitement de signaux ultrasonores qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités tout en exploitant l'acquisition avantageuse des signaux par émissions successives d'ondes planes. Il est donc proposé un procédé de traitement de signaux issus d'une 10 acquisition par sondage ultrasonore comportant les étapes suivantes : - commande de L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission, - commande de N transducteurs de réception de manière à recevoir 15 simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et - obtention d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR, J(t) de cette matrice représentant le signal de 20 mesure reçu par le i-ième transducteur de réception dû à la j-ième émission, puis - transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels, 25 - élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et - reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés. Bien que la matrice des ondes planes [MR(t)] ne soit pas de même nature que 30 les matrices de réponses impulsionnelles inter-éléments habituellement obtenues par les techniques classiques d'acquisition exploitant la focalisation synthétique en tous points, il a été constaté de façon inattendue qu'une méthode de filtrage de bruit basée sur une décomposition en valeurs singulières d'une transformée de la matrice des ondes planes dans le domaine fréquentiel fournit des résultats surprenants en 35 termes d'atténuation du bruit. La matrice débruitée [MRu(t)] ainsi obtenue permet 3060753 4 notamment une reconstitution d'image de qualité nettement améliorée par rapport à ce que produit la matrice des ondes planes [MR(t)] sans ce traitement. De façon optionnelle : - l'élimination d'une partie des valeurs singulières et vecteurs singuliers 5 associés se fait selon un critère prédéterminé de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, et - la reconstitution de la matrice débruitée [MRu(t)] de signaux temporels comporte une reconstruction d'une matrice débruitée [FTMRu(f)] de signaux fréquentiels à partir des valeurs singulières et vecteurs singuliers 10 non éliminés, puis une transformation inverse de cette matrice débruitée [FTMRu(f)] de signaux fréquentiels en la matrice débruitée [MRu(t)] de signaux temporels. De façon optionnelle également, la transformation et la transformation inverse sont des transformations de Fourier discrètes. 15 De façon optionnelle également, l'élimination d'une partie des valeurs singulières et vecteurs singuliers associés se fait par : - comparaison d'une décroissance d'amplitudes des valeurs singulières avec une courbe théorique de décroissance de valeurs singulières issues d'une matrice de bruit théorique dont les composantes sont des variables 20 aléatoires gaussiennes indépendantes, puis - suppression des valeurs singulières appartenant, à un coefficient de proportionnalité commun près, à la courbe théorique de décroissance. De façon optionnelle également, la courbe théorique de décroissance est définie par une fonction réciproque F-'(1-a) elle-même définie par une fonction F(a), 25 dite de répartition de valeurs singulières aléatoires, telle que : F(o-) = - [-a \ 14- a2 + 2 arcsin (I)] n 13 --[ - 2 2 2 où 11[0;2] désigne la fonction porte sur l'intervalle [0;2]. De façon optionnelle également, un procédé de traitement de signaux ultrasonores selon l'invention peut en outre comporter une reconstitution d'une zone imagée par calcul, en chaque point d'une pluralité de points prédéterminés de cette 30 zone imagée, d'une valeur résultant d'une sommation cohérente de valeurs instantanées prises respectivement par au moins une partie des NxM signaux temporels de la matrice [MRu(t)] à des temps de vol correspondant respectivement à un passage par le point considéré selon un mode de propagation prédéterminé. 3060753 5 De façon optionnelle également, le calcul se fait sur une partie des NxM signaux temporels de la matrice [MRu(t)] dans un secteur angulaire restreint dans l'ensemble des émissions successives, ce secteur angulaire restreint étant sélectionné de telle sorte que les ondes planes qui en sont exclues n'interagissent 5 pas avec au moins un défaut lié aux valeurs singulières et vecteurs singuliers non éliminés. De façon optionnelle également, le secteur angulaire restreint est sélectionné sur la base d'une comparaison, pour au moins l'une des valeurs singulières non éliminées, d'une valeur de phase expérimentale du vecteur singulier qui lui est 10 associé en émission avec une valeur de phase théorique en présence dudit au moins un défaut lié à ce vecteur singulier. Il est également proposé un programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions pour l'exécution des 15 étapes d'un procédé de traitement de signaux ultrasonores selon l'invention, lorsque ledit programme est exécuté sur un ordinateur. Il est également proposé un dispositif de sondage à ultrasons comportant : - une sonde comprenant L transducteurs d'émission à ultrasons et N transducteurs de réception à ultrasons, 20 - des moyens de commande des L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission, - des moyens de commande des N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque 25 émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et - un processeur de reconstitution d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR, J(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de 30 réception dû à la j-ième émission, le processeur étant en outre configuré pour effectuer les traitements suivants : - transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels, 3060753 6 - élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et - reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés. 5 L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : - la figure 1 représente schématiquement la structure générale d'un dispositif de sondage à ultrasons selon un mode de réalisation de l'invention, 10 - la figure 2 illustre un principe d'émissions successives d'ondes ultrasonores planes mis en oeuvre par le dispositif de la figure 1, - la figure 3 illustre graphiquement le résultat d'une décomposition en valeurs singulières d'une matrice de signaux fréquentiels obtenue par transformation d'une matrice d'ondes planes elle-même obtenue à l'aide 15 d'émissions réalisées selon le principe de la figure 2, - les figures 4 et 5 illustrent schématiquement et graphiquement le principe d'une sélection de secteur angulaire restreint pour améliorer la détection d'un défaut, selon un mode de réalisation de l'invention, - la figure 6 illustre trois exemples d'images reconstituées sans ou avec 20 application des principes de l'invention, et - la figure 7 illustre les étapes successives d'un procédé d'acquisition et de traitement de signaux ultrasonores mis en oeuvre par le dispositif de la figure 1, selon un mode de réalisation de l'invention. En référence à la figure 1, un dispositif de sondage 100 d'un objet 102 selon 25 un mode de réalisation de l'invention comporte une sonde à ultrasons 104 présentant un boîtier 106, c'est-à-dire un élément de structure indéformable qui sert de référentiel attaché à la sonde 104, dans lequel sont disposés, par exemple linéairement ou matriciellement, N transducteurs 1081, ... ,108N fixes ou mobiles disposés en réseau. 30 L'objet 102 est par exemple une pièce mécanique que l'on souhaite examiner par contrôle non destructif ou bien, dans un contexte médical, une partie de corps humain ou animal que l'on souhaite contrôler de manière non invasive. Dans le mode de réalisation de la figure 1, l'objet 102 est immergé dans un liquide, tel que de l'eau 110, et la sonde 104 est maintenue à distance de l'objet 102 afin que l'eau 110 les 3060753 7 sépare. Mais dans un autre mode de réalisation équivalent, la sonde 104 pourrait être au contact direct de l'objet 102. Les transducteurs 1081, ,108N sont conçus pour émettre individuellement des ondes ultrasonores en direction de l'objet 102 en réponse à des signaux de 5 commande identifiés sous la référence générale C, selon des directions principales parallèles les unes aux autres, indiquées par des flèches en pointillés sur la figure 1, et dans un plan principal qui est celui de la figure. Les transducteurs 1081, ... ,108N sont en outre conçus pour détecter des échos des ondes ultrasonores se réfléchissant sur et dans l'objet 102 et pour fournir des 10 signaux de mesure identifiés sous la référence générale S et correspondant à ces échos. Ainsi, dans l'exemple non limitatif de la figure 1, les transducteurs 1081, ,108N remplissent à la fois les fonctions d'émission et de réception, mais des récepteurs différents des émetteurs pourraient également être prévus dans des boîtiers différents et indépendants tout en restant conformes avec les principes de 15 l'invention. De plus, le nombre L d'émetteurs pourrait tout à fait être différent du nombre N de récepteurs. Le dispositif de sondage 100 comporte en outre un circuit électronique 112 de commande des transducteurs 1081, ,108N de la sonde 104 et de traitement des signaux de mesure S. Ce circuit électronique 112 est connecté à la sonde 104 afin de 20 lui transmettre les signaux de commande C et afin de recevoir les signaux de mesure S. Le circuit électronique 112 est par exemple celui d'un ordinateur. Il présente une unité centrale de traitement 114, telle qu'un microprocesseur conçu pour émettre vers la sonde 104 les signaux de commande C et pour recevoir de la sonde 104 les signaux de mesure S, et une mémoire 116 dans laquelle est notamment enregistré 25 un programme d'ordinateur 118. Le programme d'ordinateur 118 comporte tout d'abord des instructions 120 pour générer les signaux C de commande des transducteurs 1081, ,108N et recevoir leurs échos. Ces instructions sont plus précisément programmées de manière à : 30 - activer les L = N transducteurs 1081, ,108N en tant qu'émetteurs pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission de l'objet 102, - activer les transducteurs 1081, ,108N en tant que récepteurs pour, suite à chaque émission, recevoir simultanément, par ces N récepteurs et pendant 35 une durée prédéterminée de la profondeur d'inspection souhaitée, N 3060753 8 signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de chaque émission considérée. Les ondes ultrasonores planes sont obtenues à l'émission en appliquant aux transducteurs 1081, ...,108N des lois de retards enregistrées en mémoire 116 dans 5 une base 122 de lois de retards. Chaque loi de retards définit des retards à appliquer aux transducteurs 1081, ... ,108N en émission, de manière à engendrer une onde ultrasonore plane à un angle d'émission souhaité parmi les M angles d'émission successifs différents. Il est donc prévu autant de lois de retards que d'émissions successives souhaitées. 10 Comme illustré sur la figure 2 dans un cas où le nombre M d'émissions successives est impair et où les angles d'émissions se succèdent avec un pas constant dans un secteur angulaire symétrique par rapport à la direction z orthogonale au réseau de transducteurs 1081, ... ,108N, la première émission d'onde plane est associée à une loi de retards T1 portant sur des signaux émis par les 15 transducteurs 1081, ... ,108N, permettant l'émission d'une onde plane d'angle d'émission 01 par rapport à la direction z dans une première zone d'émission ZE1 partiellement située en dehors de l'ouverture de la sonde 104. La (M+1)/2-ième émission d'onde plane est associée à une loi de retards T(m+1)/2 uniforme pour l'émission d'une onde plane d'angle d'émission nul par rapport à la direction z dans 20 une (M+1)/2-ième zone d'émission ZE(M+l)/2 couvrant l'ouverture de la sonde 104. Enfin, la dernière émission d'onde plane est associée à une loi de retards Tm permettant l'émission d'une onde plane d'angle d'émission eM = - e1 par rapport à la direction z dans une dernière zone d'émission ZEM partiellement située en dehors de l'ouverture de la sonde 104. D'une façon générale, la j-ième émission d'onde plane 25 est associée à une loi de retards T, permettant l'émission d'une onde plane d'angle d'émission e, = 01 + (j - 1).( eM - ei)/(M - 1) par rapport à la direction z. En réalité, le plus souvent M est pair : il n'y a alors pas d'émission à 0° et le pas angulaire n'est pas constant. Compte tenu de la technique d'acquisition employée, la zone à imager doit 30 être contenue dans l'union des M zones d'émissions successives. Il en résulte que cette zone peut s'étendre au-delà de l'ouverture de la sonde 104, comme cela est visible sur la figure 2. En particulier, la zone imagée peut prendre la forme d'une zone sectorielle délimitée par les extrémités des zones d'émission d'angles maximal et minimal. On peut ainsi obtenir une image de type S-scan. 3060753 9 En variante, et comme avantageusement rendu possible par le principe d'acquisition des signaux par émissions successives d'ondes planes, les M angles d'émission successifs différents 01 à eM peuvent être définis autour d'une direction moyenne e(m+1)/2 non perpendiculaire au réseau de transducteurs 1081, ... ,108N. En 5 particulier, lorsqu'il s'agit de détecter des défauts tels qu'une fissure disposée au fond d'un objet à inspecter en contrôle non destructif, cette fissure étant en outre perpendiculaire au réseau de transducteurs, il est préférable de décaler latéralement la zone à inspecter par rapport à la sonde 104 et d'émettre autour d'une moyenne de 45° par exemple. La zone à inspecter peut même être décalée au point de sortir 10 complètement de l'ouverture de la sonde 104. De manière à améliorer la qualité des signaux de mesure exploités pour reconstituer la zone imagée, il est en outre possible d'appliquer une apodisation des signaux ultrasonores émis par les transducteurs 1081, ... ,108N pour former une onde ultrasonore plane de meilleure qualité, sans distorsion subie à cause des effets de 15 bords. Une telle apodisation est réalisée à l'occasion de chaque émission spatialement sur l'ensemble des transducteurs à l'aide d'une fenêtre d'apodisation telle qu'une loi d'amplitude trapézoïdale, de Hamming ou de Blackman-Harris. Elle a pour résultat de fournir une meilleure définition des zones d'émission successives. En référence de nouveau à la figure 1, à la réception des signaux résultant de 20 chacune des M émissions successives, l'ensemble S des NxM signaux temporels de mesure reçus par les N transducteurs 1081, ... ,108N est renvoyé par la sonde 104 à l'unité centrale de traitement 114. Le programme d'ordinateur 118 comporte alors en outre des instructions 124 pour construire une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, 25 qualifiée de matrice des ondes planes. Chaque coefficient MR, J(t) de cette matrice représente le signal de mesure reçu par le transducteur 108, en réponse à la j-ième émission. De façon optionnelle, le programme d'ordinateur 118 comporte en outre des instructions 126 pour effectuer un filtrage temporel de la matrice [MR(t)], ce filtrage 30 visant à supprimer toute information se trouvant à des temps de vol exclus de la zone d'intérêt dans l'objet 102. Le programme d'ordinateur 118 comporte en outre des instructions 128 pour transformer la matrice [MR(t)] en une matrice [FTMR(f)] de signaux fréquentiels par transformée de Fourier, avantageusement par transformée de Fourier discrète après 35 échantillonnage temporel des signaux ultrasonores formant les coefficients de la 3060753 10 matrice [MR(t)], ou, encore plus avantageusement encore, par calcul de FFT (de l'anglais « Fast Fourier Transform ») si le nombre d'échantillons de chaque coefficient de la matrice [MR(t)] le permet. Le programme d'ordinateur 118 comporte en outre des instructions 130 pour 5 décomposer la matrice [FTMR(f)] de signaux fréquentiels en valeurs singulières sur une bande de fréquences. Bien qu'il soit connu de décomposer en valeurs singulières une matrice de réponses impulsionnelles inter-éléments habituellement obtenues par les techniques classiques d'acquisition exploitant la focalisation synthétique en tous points, cette opération n'est pas équivalente lorsqu'elle est appliquée sur une matrice 10 telle que la matrice [FTMR(f)]. En effet, la matrice notée [PFTMR(f)] de dimensions MxM et définie par le produit [FTMR(f)]t.[FTMR(f)], où « t » est le symbole de la transposée conjuguée d'une matrice, ne représente pas l'opérateur de retournement temporel en émission comme pour la matrice des réponses impulsionnelles inter-éléments. L'interprétation physique des valeurs singulières et vecteurs singuliers 15 n'est alors pas la même que dans la méthode de Décomposition de l'Opérateur de Retournement Temporel (méthode DORT) telle qu'enseignée dans l'article de Prada et al, intitulé « Eigenmodes of the time reversai operator: a solution to selective multiple-target media », publié dans Wave Motion 20, pages 151-163 (1994). De plus, la matrice [FTMR(f)] n'est généralement ni carrée ni symétrique, M étant souvent 20 même nettement inférieur à N en acquisition par émission d'ondes planes. De la sorte, une décomposition en valeurs singulières ne produit pas de diagonalisation comme dans la plupart des travaux portant sur la matrice des réponses impulsionnelles inter-éléments. Plus précisément, l'opération de décomposition en valeurs singulières permet 25 d'estimer les matrices U(f), S(f) et V(f) telles que : [FTMR(0]=U(f)S(f)Vt(f)= 61(f). (f)v:(f), où U(f) = [u1(f), uNW] et -1 V=[v1(f), vM(f)] sont des matrices unitaires orthogonales de tailles respectives NxN et MxM, qui contiennent respectivement les vecteurs singuliers en réception et émission, où S est une matrice réelle de taille NxM contenant les K = min(N;M) 30 valeurs singulières ce de la matrice [FTMR(f)], ordonnées de façon décroissante à une fréquence f donnée ai(f) 6K(f) k 0, et où min(N;M) est la fonction qui retourne la valeur minimale entre N et M (généralement c'est M). Le programme d'ordinateur 118 comporte en outre des instructions 132 pour réduire le rang de la matrice [FTMR(f)], en éliminant une partie des valeurs 3060753 11 singulières a,(f). Cette élimination se fait selon un critère de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, les premières étant d'amplitudes supérieures aux secondes. Compte tenu du fait que ai(f) GK(f) k 0, il s'agit de trouver la valeur entière P(f) comprise entre 1 et K telle que ai (f), 5 ..., ap(o(f) peuvent être considérées comme liées à des défauts à détecter dans l'objet 102 et ap(0.1(f) ---, aK(t) peuvent être éliminées car considérées comme liées à du bruit. Dans le cas de petits défauts idéalement espacés les uns des autres, P(f) est égal au nombre de défauts présents dans la zone d'intérêt inspectée. On peut alors écrire la matrice [FTMR(f)] sous la forme d'une addition de deux matrices [FTMRu(f)] 10 et [FTMRN(f)] [FTMR( f )]= [FTMRu ( f )]+ [FTMRN ( f )], avec P(f) [FTMRu = 6 f f (f ), appelée matrice de signal utile, et [FTMRN f )] = Io-,(f).11,(f).v:(f), appelée matrice de bruit. =P( f)+1 Réduire le rang de la matrice [FTMR(f)] revient ainsi à ne conserver que 15 [FTMRuM]. En pratique, la détermination de P(f) se fait par une étude de la courbe de décroissance des amplitudes des valeurs singulières. Cette étude peut se faire sur la base de la théorie des matrices aléatoires telle qu'enseignée dans : - l'article de Martenko et al, intitulé « Distribution of eigenvalues for some 20 sets of random matrices », publié dans Mathematics of the USSR-Sbornik, vol. 1, n°4, pages 457-483 (1967), et - l'article de Aubry et al, intitulé « Detection and imaging in a random medium: a matrix method to overcome multiple scattering and aberration », publié dans Journal of Applied Physics, 106(4), 044903 (2009). 25 En supposant que toutes les composantes de la matrice de bruit [FTMRN(f)] sont des variables aléatoires gaussiennes indépendantes, il peut en effet être montré que les valeurs singulières liées au bruit Opm-re aK(f) appartiennent, à un coefficient de proportionnalité commun près, à une courbe théorique définie par une fonction réciproque F-1(1-a) elle-même définie par une fonction F(a), dite de 30 répartition de valeurs singulières aléatoires, avec : F(6) = - [-a,\/4- a2 + 2 arcsin (I)]132 n -[ - ], 2 2 OÙ 11[0;2] désigne la fonction porte sur l'intervalle [0;2]. 3060753 12 La fonction F(a) donne des valeurs comprises entre 0 et 1 sur l'intervalle support G E [0;2], où G est la variable de valeur singulière. Elle est par ailleurs strictement croissante. Sa fonction réciproque F-1(a) est donc également strictement croissante, de sorte que la fonction F-1(1-a) donne la courbe de décroissance des 5 valeurs singulières recherchée, à une constante de proportionnalité près pour l'ajuster à la courbe expérimentale. Ainsi, le nombre de valeurs singulières calculées par exécution des instructions 130 n'appartenant pas à cette courbe théorique définit le rang P(f) de la matrice [FTMRu(f)] de signal utile. A titre d'exemple, la figure 3 illustre en unité 10 arbitraire a.u. (de l'anglais « arbitrary unit ») une distribution de valeurs singulières calculées à une fréquence f = 5 MHz, pour une acquisition réalisée avec un échantillon de polyéthylène en émettant M = 64 ondes planes entre -31,5° et +31,5°. Deux valeurs singulières se détachent de la courbe théorique, dont une assez nettement, indiquant dans cet exemple que P(f) = 2 à 5 MHz. 15 Pour calculer la matrice [FTMRu(f)] de signal utile sur plusieurs fréquences f, il suffit d'identifier à chaque fréquence souhaitée le nombre de valeurs singulières n'appartenant pas à la courbe théorique. La matrice [FTMRu(f)] ainsi reconstituée est une matrice débruitée de signaux fréquentiels, le sous-espace de bruit représenté par la matrice [FTMRN(f)] ayant été 20 éliminé. Le programme d'ordinateur 118 comporte en outre des instructions 134 pour transformer la matrice [FTMRu(f)] en une matrice débruitée [MRu(t)] de signaux temporels par transformée de Fourier inverse, avantageusement par transformée de Fourier discrète inverse, ou, encore plus avantageusement encore, par calcul de 25 IFFT (de l'anglais « Inverse Fast Fourier Transform ») si le nombre d'échantillons de chaque coefficient de la matrice [FTMRu(f)] le permet. Enfin, le programme d'ordinateur 118 comporte des instructions, désignées par la référence générale 136, de traitement de la matrice [MRu(t)]. Le traitement réalisé par les instructions 136 peut inclure une reconstitution d'image numérique de 30 la zone d'intérêt dans l'objet 102 par adaptation du principe de focalisation synthétique en tous points, comme enseigné par exemple dans le document WO 2015/092250 Al. On reconstitue ainsi une image numérique de la zone d'intérêt dont la qualité est meilleure que si la reconstitution avait été réalisée sur la matrice non débruitée [MR(t)]. Notamment, le Rapport Signal sur Bruit (RSB) est amélioré. En 3060753 13 variante ou en complément, le traitement réalisé par les instructions 136 pourrait inclure une focalisation adaptative et sélective. Selon une adaptation relativement simple mais pas toujours optimale de la focalisation synthétique en tous points, en chaque pixel P de l'image numérique 5 reconstituée, on calcule le module A(P) d'une sommation cohérente impliquant les NxM signaux temporels de la matrice [MRu(t)] à NxM temps de vols calculés selon un mode de propagation prédéterminé, chaque temps de vol tij étant le temps mis par la j-ième onde plane pour être reçue par le i-ième transducteur de réception en passant par le pixel considéré selon le mode de propagation prédéterminé : A(P) =Iii(13)Pi(P)MRu-,].[ti'l(131, j =1 i=1 10 où µj(P) et pi(P) sont des coefficients de pondération respectivement en émission et en réception dont les expressions dépendent de l'application considérée pour prendre en compte des phénomènes ou traitements tels qu'un filtrage d'échos de géométrie, une compensation d'atténuation due à un étalement spatial des ondes, etc. Ce mode de calcul n'est souvent pas optimal parce que, comme l'illustre la 15 figure 4, toutes les ondes planes émises par la sonde 104 n'interagissent pas avec le ou les défauts D liés aux valeurs singulières de la matrice de signal utile [FTMRu(f)], c'est-à-dire les valeurs singulières Ge, --- , Gp(0(f). Il peut alors être avantageux de sélectionner un secteur angulaire restreint [0,1;0,2] dans l'ensemble des émissions successives à l'extérieur duquel les ondes planes n'interagissent pas avec le ou les 20 défauts précités. Ce secteur angulaire retreint est défini par ses indices d'émissions minimal (m1) et maximal (m2) où 1erl1m2M. Le module A(P) est alors de préférence calculé par l'addition cohérente suivante : A(P) = /1712 /N iii(P)Pi(P)MRu J=mi i=1 Dans le cas d'un défaut assimilable à un réflecteur ponctuel, une méthode astucieuse pour déterminer les valeurs de m1 et m2 va maintenant être détaillée en 25 référence à la figure 5. Elle consiste à comparer, pour au moins l'une des valeurs singulières ai (f), Gp(0(f), et notamment au moins pour la première valeur singulière ai(f) qui est la plus élevée, une valeur de phase expérimentale du vecteur singulier qui lui est associé en émission avec une valeur de phase théorique en présence du ou des défauts D. Cette comparaison est réalisée à une fréquence fc choisie qui peut 30 être la fréquence centrale de fonctionnement de la sonde 104, une fréquence pour laquelle la valeur singulière a1(f) prend sa valeur la plus élevée, ou toute autre 3060753 14 fréquence prédéterminée. Sur la figure 5, la comparaison est ainsi réalisée à 5 MHz pour le vecteur singulier en émission v1(f) correspondant à un emplacement de défaut(s) D de coordonnées (XD,Zo). Dans les conditions expérimentales indiquées ci-dessus, la valeur de phase 5 théorique est une phase corrigée calculée de la façon suivante pour être toujours négative : 2n-fc.(XD sin Oi + ZD cos e) rn-fc.(XD sin Oi + ZD cos el t j,1 j M,(pjh = max avec c la célérité de l'onde plane dans le milieu considéré et son angle d'émission. Les coordonnées (XD,ZD) peuvent être déterminées avec la matrice de signal utile [FTMRu(f)] en calculant la rétro-propagation du vecteur singulier en réception 10 u1 (f) à la fréquence fc, comme par exemple enseigné dans l'article de Lopez Villaverde et al, intitulé « Ultrasonic imaging of defects in coarse-grained steels with the decomposition of the time reversai operator », publié dans Journal of the Acoustical Society of America, volume 140, n°1, pages 541-550 (2016). La définition ci-dessus de la phase théorique corrigée donne la courbe en trait 15 continu de la figure 5. La valeur de phase expérimentale du vecteur singulier en émission vi(fc) est également une phase corrigée calculée de la façon suivante pour être toujours négative : ],1 REVENDICATIONS1. Procédé (700) de traitement de signaux (S) issus d'une acquisition par sondage ultrasonore comportant les étapes suivantes : - commande (702) de L transducteurs d'émission (1081, ..., 108N) pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents (01, eM) dans M zones d'émission (ZE1, ZEM), - commande (702) de N transducteurs de réception (1081, ..., 108N) de manière à recevoir (704) simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, - obtention (704) d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR, J(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission, caractérisé en ce qu'il comporte en outre les étapes suivantes : - transformation (708) de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières (710) de la matrice [FTMR(f)] de signaux fréquentiels, - élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières (710), et - reconstitution (712, 714) d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés. 2. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 1, dans lequel : - l'élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés se fait selon un critère prédéterminé de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, et 3060753 19 - la reconstitution (712, 714) de la matrice débruitée [MRu(t)] de signaux temporels comporte une reconstruction (712) d'une matrice débruitée [FTMRu(f)] de signaux fréquentiels à partir des valeurs singulières et vecteurs singuliers non éliminés, puis une 5 transformation inverse (714) de cette matrice débruitée [FTMRu(f)] de signaux fréquentiels en la matrice débruitée [MRu(t)] de signaux temporels. 3. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 2, dans lequel la transformation (708) et la transformation inverse (714) 10 sont des transformations de Fourier discrètes. 4. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 2 ou 3, dans lequel l'élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés se fait par : - comparaison d'une décroissance d'amplitudes des valeurs 15 singulières avec une courbe théorique de décroissance de valeurs singulières issues d'une matrice de bruit théorique dont les composantes sont des variables aléatoires gaussiennes indépendantes, puis - suppression des valeurs singulières appartenant, à un coefficient de 20 proportionnalité commun près, à la courbe théorique de décroissance. 5. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 4, dans lequel la courbe théorique de décroissance est définie par une fonction réciproque F-'(1-a) elle-même définie par une fonction F(a), dite de 25 répartition de valeurs singulières aléatoires, telle que : F(o-) = - [-a \ 14- a2 + 2 arcsin (I)] n 13 --[ - 2 2 2 où 11[0;2] désigne la fonction porte sur l'intervalle [0;2]. 6. Procédé (700) de traitement de signaux ultrasonores (S) selon l'une quelconque des revendications 1 à 5, comportant en outre une reconstitution (716, 718) d'une zone imagée par calcul, en chaque point d'une pluralité de points 30 prédéterminés de cette zone imagée, d'une valeur résultant d'une sommation cohérente de valeurs instantanées prises respectivement par au moins une partie des NxM signaux temporels de la matrice [MRu(t)] à des temps de vol correspondant respectivement à un passage par le point considéré selon un mode de propagation prédéterminé. 3060753 20 7. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 6, dans lequel le calcul se fait sur une partie des NxM signaux temporels de la matrice [MRu(t)] dans un secteur angulaire restreint (m1, m2) dans l'ensemble des émissions successives, ce secteur angulaire restreint (m1, m2) étant 5 sélectionné (716) de telle sorte que les ondes planes qui en sont exclues n'interagissent pas avec au moins un défaut lié aux valeurs singulières et vecteurs singuliers non éliminés. 8. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 7, dans lequel le secteur angulaire restreint (m1, m2) est sélectionné 10 (716) sur la base d'une comparaison, pour au moins l'une des valeurs singulières non éliminées, d'une valeur de phase expérimentale du vecteur singulier qui lui est associé en émission avec une valeur de phase théorique en présence dudit au moins un défaut lié à ce vecteur singulier. 9. Programme d'ordinateur (118) téléchargeable depuis un réseau de 15 communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur (114), caractérisé en ce qu'il comprend des instructions (120, 124, 126, 128, 130, 132, 134, 136) pour l'exécution des étapes d'un procédé (700) de traitement de signaux ultrasonores (S) selon l'une quelconque des revendications 1 à 8, lorsque ledit programme est exécuté sur un ordinateur (112). 20 10. Dispositif (100) de sondage à ultrasons comportant : - une sonde (104) comprenant L transducteurs d'émission à ultrasons (1081, ..., 108N) et N transducteurs de réception à ultrasons (1081, , 108N), - des moyens (112) de commande des L transducteurs d'émission 25 pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents (01, eM) dans M zones d'émission (ZE1, ZEM), - des moyens (112) de commande des N transducteurs de réception de manière à recevoir simultanément et pendant une durée 30 prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et - un processeur (114) de reconstitution d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient 3060753 21 MR, J(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission, caractérisé en ce que le processeur (114) est en outre configuré pour effectuer les traitements suivants : 5 - transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels - élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et 10 - reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.