Brevet FR3060768A1 — publié le 20180622

La présente invention concerne un procédé d'acquisition de signaux par sondage ultrasonore, par exemple 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é d'acquisition de signaux ultrasonores opérant de la façon suivante : - commande de L transducteurs d'émission pour M émissions successives d'ondes ultrasonores vers une zone d'intérêt, - commande de N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chacune des M émissions successives, N signaux de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée dans la zone d'intérêt, - obtention d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MRi,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 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, ce traitement peut être embarqué dans les instruments de contrôle pour des traitements en temps réel. 3060768 2 Dans la pratique, la commande des L transducteurs d'émission pour les M émissions successives d'ondes ultrasonores vers la zone d'intérêt peut être codée à l'aide d'une matrice de codage [MC], chaque coefficient MC,,, de cette matrice représentant un facteur multiplicatif appliqué à un signal temporel d'excitation 5 commun e(t) pour son émission par le i-ième transducteur d'émission à l'occasion de la j-ième émission. Des lois de retards peuvent en outre être appliquées aux émissions successives. Lorsque [MC] est la matrice identité et qu'aucune loi de retard n'est appliquée, l'acquisition ultrasonore précédemment définie est qualifiée d'acquisition FMC (de 10 l'anglais « Full Matrix Capture »). Elle consiste à émettre une onde ultrasonore en excitant le premier transducteur d'émission et à recevoir les échos de cette émission avec l'ensemble des N transducteurs de réception, puis à commuter électroniquement sur l'ensemble des transducteurs d'émission pour les exciter successivement. Quand ce sont les mêmes transducteurs qui remplissent les 15 fonctions d'émission et de réception, on obtient une matrice [MR(t)] notée [K(t)] de signaux temporels ultrasonores de taille NxN. En imagerie, les signaux temporels ultrasonores formant les coefficients de la matrice [K(t)] sont exploités pour réaliser une focalisation synthétique de type « focalisation en tous points » qui permet d'obtenir une image de résolution optimale 20 dans toute la zone d'intérêt. Cependant, en présence d'un fort bruit électronique ou de structure, la reconstruction par focalisation en tous points peut fournir des images de qualité moindre comparativement aux méthodes échographiques classiques. En effet, dans ces dernières, tous les transducteurs émettent de manière simultanée par application d'une loi de retards prédéterminée de manière à focaliser en un point 25 donné. Or, conformément au procédé d'acquisition FMC généralement mis en oeuvre pour réaliser ensuite la reconstruction par focalisation synthétique, chaque émission est effectuée par un seul transducteur ce qui limite l'énergie transmise et la profondeur de pénétration des ondes dans la pièce inspectée. Cela se traduit finalement par une dégradation du Rapport Signal sur Bruit (RSB) car les amplitudes 30 des signaux échos peuvent être comparables au niveau de bruit électronique ou de structure. Cette dégradation du RSB est d'autant plus forte que l'atténuation des ondes ultrasonores est élevée (atténuation viscoélastique, ou atténuation par diffusion due à des hétérogénéités du milieu), rendant difficile la détection et la caractérisation d'éventuels défauts. 3060768 3 Une solution partielle à ce problème de dégradation du RSB est apportée dans l'article de Karaman et al, intitulé « Synthetic aperture imaging for small scale systems », publié dans IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 42, n°3 (mai 1995), pages 429-442. 5 Elle consiste à utiliser, pour chaque émission, non plus un transducteur mais plusieurs transducteurs adjacents. Une loi de retards est appliquée aux transducteurs d'émission sollicités pour qu'ils transmettent dans le milieu une onde ultrasonore sphérique, proche de celle qui serait émise par une source virtuelle située à une certaine distance du capteur. L'onde ultrasonore ainsi émise par la source virtuelle 10 est plus intense puisque son énergie est proportionnelle à la racine carrée du nombre de transducteurs d'émission formant cette source. Le RSB est amélioré d'autant, en supposant que le bruit engendré est principalement du bruit électronique décorrélé. La matrice de codage [MC] correspondante comporte plusieurs coefficients non nuls dans chacune de ses colonnes, définissant précisément le nombre de transducteurs 15 adjacents sollicités en même temps à chaque émission et les lois de retards éventuellement appliquées. Cependant, dans le cas de pièces inspectées présentant un bruit très important, l'amélioration de la qualité des images finalement obtenues par focalisation synthétique en tous points est plus limitée, l'augmentation du RSB est moindre et 20 l'impact sur la détection n'est pas aussi positif qu'on pourrait l'espérer. Cette solution compense en partie le problème précité mais ne l'élimine pas. En effet, l'augmentation du RSB est proportionnelle à la racine carrée du nombre de transducteurs formant chaque source virtuelle, et ce nombre est beaucoup plus petit que le nombre total de transducteurs N du capteur. Par ailleurs, émettre au moyen de 25 sources virtuelles ne permet pas de s'affranchir du problème que peuvent poser les artefacts de reconstruction essentiellement dus aux échos parasites tels que les échos de géométrie ou les échos complexes incluant des réflexions multiples sur les frontières de l'objet et des conversions de modes. Une autre solution à ce problème de dégradation du RSB est apportée dans 30 l'article de Lopez Villaverde et al, intitulé « Ultrasonic imaging in highly attenuating materials with VVhalsh-Hadamard codes and the decomposition of the time reversai operator », publié à l'occasion de la conférence IEEE International Ultrasonics Symposium qui s'est tenue à Tours (FR) du 18 au 21 septembre 2016. Elle consiste à coder les émissions à l'aide d'une matrice [MC] de Hadamard notée HN ou à l'aide de matrices issues de cette matrice de Hadamard (Ie, et/ou HN-). Il a été démontré 3060768 4 donc possible d'accroître la taille des évents sans risquer de déformer le récipient. Diverses caractéristiques supplémentaires peuvent être prévues, seules ou en combinaison : 5 - en position déployée, la face terminale de la (ou de chaque) saillie s'étend dans le prolongement de la face de moulage ; en position déployée, la face terminale de la (ou de chaque) saillie s'étend en retrait, par rapport à l'ouverture interne, d'une distance comprise entre 0,5 mm et 5 mm ; 10 - la face terminale est concave, ou convexe ; l'insert comprend un piston monté en translation dans une chemise solidaire du bloc de fond ; le ou chaque évent comprend une portion interne, qui s'étend jusqu'à l'ouverture, et une portion externe de plus grande section, qui 15 s'étend dans le prolongement de la portion interne en étant séparée de celle-ci par un épaulement ; le fond de moule comprend un conduit primaire de mise à l'air, qui débouche dans l'évent au niveau de l'épaulement ; le fond de moule comprend un conduit secondaire de mise à l'air, 20 qui débouche dans l'évent au niveau de la portion externe ; le bloc est pourvu d'une pluralité d'évents espacés les uns des autres, et l'insert comprend une série périphérique de saillies espacées les unes des autres ; la face de moulage est à l'empreinte d'un fond pétaloïde de récipient 25 et comprend une alternance de nervures à l'empreinte de vallées du fond pétaloïde, qui rayonnent à partir d'une zone centrale, et de réserves en creux à l'empreinte de pieds du fond pétaloïde, qui s'étendent entre les nervures, les évents sont pratiqués dans les réserves en creux, et la face terminale de chaque saillie est à 30 l'empreinte des pieds du récipient. Il est proposé, en deuxième lieu, un moule pour la fabrication d'un récipient à partir d'une ébauche en matière plastique, qui comprend une paroi latérale à l'empreinte du corps du récipient, et un fond de moule tel que présenté ci-dessus, qui complète avec la paroi latérale l'empreinte 35 du récipient. 3060768 5 - la commande des L transducteurs d'émission pour les M émissions successives d'ondes ultrasonores vers la zone d'intérêt est codée à l'aide de la matrice réduite de codage [MC] appliquée au signal temporel d'excitation commun e(t). 5 Ainsi, en partant d'une matrice de codage [MC'] qui peut être l'une de celles indiquées précédemment ou autre, notamment une matrice de codage assurant un RSB satisfaisant selon l'application visée et les traitements ultérieurs envisagés, l'invention permet de réduire le nombre d'émissions en éliminant une ou plusieurs émissions sur la base d'un critère pertinent puisqu'il est lié à des calculs de champs 10 acoustiques. Cette manière de procéder permet notamment de réduire le nombre des émissions en limitant au maximum l'impact de cette réduction sur le RSB de la matrice [MR(t)]. De façon optionnelle, la matrice [MR(t)] de signaux temporels ultrasonores est décodée pour l'obtention d'une matrice décodée [MR'(t)] calculée par produit matriciel 15 de la façon suivante : [MR'(t)] = [MR(t)] - [MC]T - ([MC] - [MC]T) 1, où « T » est le symbole de transposition matricielle. On obtient ainsi une matrice de taille NxL à partir de laquelle peuvent être envisagés les mêmes traitements que sur une matrice de signaux temporels obtenue par une acquisition classique de type FMC. 20 De façon optionnelle également, le calcul de champ acoustique réalisé pour chacune des M' émissions initiales successives comporte le calcul d'un modèle simplifié de champ Ei',(f,0) défini pour chaque colonne d'indice m' de la matrice initiale de codage [MC'] de la façon suivante : fL =Im c, on, . si(f) . Diu 0) e)2-Jkla sin e 1=1 où : 25 - f est une fréquence temporelle, - 0 est un angle par rapport à une normale à un plan ou axe principal des L transducteurs d'émission lorsque ceux-ci sont alignés, - si(f) est une fonction de transfert d'un I-ième transducteur d'émission, - Di(f,0) est une fonction de directivité du 1-ième transducteur d'émission 30 dans un milieu d'émission des ondes ultrasonores, - e est la fonction exponentielle, - j est le nombre complexe tel que j2 = -1, 3060768 6 - k est un nombre d'onde défini par k = 2rrf/c où c est la célérité des ondes ultrasonores dans le milieu d'émission considéré, et - d est un pas inter-éléments, c'est-à-dire une largeur commune des transducteurs d'émission ajoutée à une distance entre deux transducteurs 5 voisins. De façon optionnelle également, le calcul de champ acoustique réalisé pour chacune des M' émissions initiales successives comporte en outre le calcul d'une valeur intégrée de champ A,',(0) à partir de chaque modèle simplifié de champ E,',(f,e) de la façon suivante : L fmax Imcl,m, - f.S1 ) - Di (f, e)e -jad sin 6+ 1=1 f min fmax Am,(0) Em,(f,O)df= fmin 10 où fmin et fmax sont des fréquences respectivement minimale et maximale d'une bande passante du signal temporel d'excitation commun e(t). De façon optionnelle également, le critère de sélection appliqué aux M' calculs de champs acoustiques comporte un seuil en amplitude en dessous duquel les contributions du champ acoustique sont considérées comme négligeables. 15 De façon optionnelle également, chaque colonne de la matrice initiale de codage [MC'] produisant une émission initiale dont le calcul de champ acoustique ne fournit pas de valeur supérieure ou égale au seuil en amplitude est supprimée. De façon optionnelle également, le critère de sélection appliqué aux M' calculs de champs acoustiques comporte en outre un seuillage angulaire consistant à écarter 20 toute contribution du champ acoustique en dehors d'un secteur angulaire prédéterminé. De façon optionnelle également, la matrice initiale de codage [MC'] est une matrice de Hadamard ou obtenue à partir d'une matrice de Hadamard. Il est également proposé un programme d'ordinateur téléchargeable depuis un 25 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 étapes d'un procédé d'acquisition de signaux selon l'invention, lorsque ledit programme est exécuté sur un ordinateur. Il est également proposé un dispositif de sondage à ultrasons comportant : 30 - une sonde comprenant une pluralité de transducteurs d'émission à ultrasons et une pluralité de transducteurs de réception à ultrasons, et - des moyens de commande des transducteurs et de traitement conçus pour mettre en oeuvre un procédé d'acquisition de signaux selon l'invention. 3060768 7 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 5 de sondage à ultrasons selon un mode de réalisation de l'invention, - les figures 2A, 2B et 2C illustrent trois exemples de champs acoustiques correspondant à respectivement trois émissions codées réalisées par un dispositif de sondage tel que celui de la figure 1, - les figures 3A, 3B et 3C illustrent trois exemples de calculs de champs 10 acoustiques correspondant à respectivement trois émissions codées réalisées par un dispositif de sondage tel que celui de la figure 1, - la figure 4 illustre un autre exemple de calcul de champ acoustique exploitant notamment les calculs des figures 3A, 3B et 3C, - la figure 5 illustre l'effet d'un seuillage réalisé sur le résultat du calcul de la 15 figure 4, - la figure 6 illustre l'effet d'une binarisation réalisée sur le résultat de la figure 5 en vue d'une suppression d'émissions codées, - la figure 7 représente schématiquement une installation expérimentale de tests comparatifs, 20 - la figure 8 est un diagramme illustrant des résultats de tests comparatifs obtenus avec ou sans dispositif de sondage selon l'invention sur l'installation expérimentale de la figure 7, - les figures 9A, 9B et 9C illustrent trois images obtenues par focalisation synthétique en tous points sur des signaux temporels obtenus par une 25 acquisition avec ou sans dispositif de sondage selon l'invention sur l'installation expérimentale de la figure 7, et - la figure 10 illustre les étapes d'un procédé d'acquisition de signaux ultrasonores mis en oeuvre par le dispositif de la figure 1. En référence à la figure 1, un dispositif de sondage 100 d'un objet 102 selon 30 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. L'objet 102 est par exemple une pièce mécanique que l'on souhaite examiner 35 par contrôle non destructif ou bien, dans un contexte médical, une partie de corps 3060768 8 humain 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 sépare. Mais dans un autre mode de réalisation équivalent, la sonde 104 pourrait être 5 au contact direct de l'objet 102. Les transducteurs 1081, ,108N sont conçus pour émettre des ondes ultrasonores en direction de l'objet 102 en réponse à des signaux de 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 10 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 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 N transducteurs 1081, 15 ... ,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 l'invention. De plus, le nombre L d'émetteurs pourrait tout à fait être différent du nombre N de récepteurs. 20 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 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 25 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 enregistré un programme d'ordinateur 118. Le programme d'ordinateur 118 comporte tout d'abord des instructions 120 30 pour définir M' émissions initiales successives à l'aide d'une matrice initiale de codage [MC'] de taille LxM', c'est-à-dire de taille NxM' dans l'exemple non limitatif considéré. Chaque coefficient de cette matrice représente un facteur multiplicatif appliqué à un signal temporel d'excitation e(t), commun à tous les transducteurs 1081, ,108N, pour son émission par le i-ième transducteur d'émission 35 à l'occasion de la j-ième émission. Ce facteur multiplicatif peut inclure un retard d'une 3060768 9 loi de retards appliquée à la j-ième émission initiale concernée. La matrice initiale de codage [MC'] peut être prédéterminée et enregistrée en mémoire, choisie à l'aide des instructions 120 parmi un ensemble de matrices initiales de codage enregistrées en mémoire, définie via une interface homme-machine à l'aide des instructions 120, etc. 5 Comme indiqué précédemment, la matrice initiale de codage [MC'] peut être la matrice identité pour M' = L = N émissions initiales successives conformes à une acquisition FMC. Elle peut aussi être une matrice dont les coefficients non nuls incluant des lois de retards sont situés autour de sa diagonale principale, pour M' REVENDICATIONS1. Procédé (200) d'acquisition de signaux par sondage ultrasonore comportant les étapes suivantes : - commande (208) de L transducteurs d'émission (1081, ..., 108N) pour M émissions successives d'ondes ultrasonores vers une zone d'intérêt (102), - commande (208) de N transducteurs de réception (1081, ..., 108N) de manière à recevoir simultanément et pendant une durée prédéterminée, pour chacune des M émissions successives, N signaux de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée dans la zone d'intérêt, - obtention (210) 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 que : - une matrice initiale [MC'], de taille LxM', de codage des émissions successives est préalablement définie (202) pour un nombre M' d'émissions initiales successives strictement supérieur à M, chaque coefficient de cette matrice représentant un facteur multiplicatif appliqué à un signal temporel d'excitation commun e(t) pour son émission par le i-ième transducteur d'émission à l'occasion de la jième émission, - un calcul de champ acoustique est réalisé (204) pour chacune des M' émissions initiales successives, - une matrice réduite de codage [MC], de taille LxM, est obtenue (206) à partir de la matrice initiale de codage [MC'] par suppression de M'- M colonne(s) correspondant à M'-M émissions initiale(s) supprimée(s) sur la base d'un critère de sélection appliqué aux M' calculs de champs acoustiques, et - la commande (208) des L transducteurs d'émission pour les M émissions successives d'ondes ultrasonores vers la zone d'intérêt est codée à l'aide de la matrice réduite de codage [MC] appliquée au signal temporel d'excitation commun e(t). 3060768 19 2. Procédé (200) d'acquisition de signaux selon la revendication 1, dans lequel la matrice [MR(t)] de signaux temporels ultrasonores est décodée (214) pour l'obtention d'une matrice décodée [MR'(t)] calculée par produit matriciel de la façon suivante : [Me (t)] = [MR(t)] - [M Cr - ([MC] - [M , 5 où « T » est le symbole de transposition matricielle. 3. Procédé (200) d'acquisition de signaux selon la revendication 1 ou 2, dans lequel le calcul de champ acoustique réalisé (204) pour chacune des M' émissions initiales successives comporte le calcul d'un modèle simplifié de champ Ei',(f,0) défini pour chaque colonne d'indice m' de la matrice initiale de codage [MC'] 10 de la façon suivante : L Em, (f, 0) = . (f) . e)e-fiad sin , 1=1 où : - f est une fréquence temporelle, - 0 est un angle par rapport à une normale à un plan ou axe principal des L transducteurs d'émission (1081, ..., 108N) lorsque ceux-ci sont 15 alignés, - si(f) est une fonction de transfert d'un l-ième transducteur d'émission, - DI(f,0) est une fonction de directivité du 1-ième transducteur d'émission dans un milieu d'émission des ondes ultrasonores, - e est la fonction exponentielle, 20 - j est le nombre complexe tel que j2 = -1, - k est un nombre d'onde défini par k = 2rrf/c où c est la célérité des ondes ultrasonores dans le milieu d'émission considéré, et - d est un pas inter-éléments, c'est-à-dire une largeur commune des transducteurs d'émission (1081, ..., 108N) ajoutée à une distance 25 entre deux transducteurs voisins. 4. Procédé (200) d'acquisition de signaux selon la revendication 3, dans lequel le calcul de champ acoustique réalisé (204) pour chacune des M' émissions initiales successives comporte en outre le calcul d'une valeur intégrée de champ A,',(0) à partir de chaque modèle simplifié de champ E,',(f,0) de la façon suivante : fmax MCI,m, 1rnaxfmax s1(f) - DI(f e)e- Pad sin e f min Am,(0) = Em,(f, , 0)df= 1=1 fmin 3060768 20 Où fmin et furax sont des fréquences respectivement minimale et maximale d'une bande passante du signal temporel d'excitation commun e(t). 5. Procédé (200) d'acquisition de signaux selon l'une quelconque des revendications 1 à 4, dans lequel le critère de sélection appliqué aux M' calculs de 5 champs acoustiques comporte un seuil en amplitude en dessous duquel les contributions du champ acoustique sont considérées comme négligeables. 6. Procédé (200) d'acquisition de signaux selon la revendication 5, dans lequel chaque colonne de la matrice initiale de codage [MC'] produisant une émission initiale dont le calcul de champ acoustique ne fournit pas de valeur supérieure ou 10 égale au seuil en amplitude est supprimée (206). 7. Procédé (200) d'acquisition de signaux selon l'une quelconque des revendications 1 à 6, dans lequel le critère de sélection appliqué aux M' calculs de champs acoustiques comporte en outre un seuillage angulaire consistant à écarter toute contribution du champ acoustique en dehors d'un secteur angulaire 15 prédéterminé. 8. Procédé (200) d'acquisition de signaux selon l'une quelconque des revendications 1 à 7, dans lequel la matrice initiale de codage [MC'] est une matrice de Hadamard ou obtenue à partir d'une matrice de Hadamard. 9. Programme d'ordinateur (118) téléchargeable depuis un réseau de 20 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 pour l'exécution des étapes d'un procédé (200) d'acquisition de signaux selon l'une quelconque des revendications 1 à 8, lorsque ledit programme est exécuté sur un ordinateur (112). 25 10. Dispositif de sondage à ultrasons (100) comportant : - une sonde (104) comprenant une pluralité de transducteurs (1081, ..., 108N) d'émission à ultrasons et une pluralité de transducteurs (1081, ..., 108N) de réception à ultrasons, et - des moyens (112) de commande des transducteurs et de traitement 30 conçus pour mettre en oeuvre un procédé (200) d'acquisition de signaux selon l'une quelconque des revendications 1 à 8.