FR 2999292 A1 20140613 FR 1261802 A 20121207 DISPOSITIF ET PROCÉDÉ D'IMAGERIE PAR ULTRASONS AVEC FILTRAGE DES ARTEFACTS DUS AUX INTERFÉRENCES ENTRE MODES DE RECONSTRUCTION L'invention se situe dans le domaine de l'imagerie par ultrasons. Elle s'applique en particulier au contrôle non destructif de pièces. Le contrôle non destructif vise notamment à détecter la présence d'éventuels défauts dans une pièce, à les localiser et à les dimensionner. L'invention concerne plus particulièrement un procédé d'imagerie par ultrasons utilisant un traducteur multiélément et basé sur la technique de reconstruction appelée "focalisation en tous points". Dans le domaine du contrôle non destructif par ultrasons, l'utilisation de traducteurs multiéléments a considérablement augmenté les possibilités d'inspection de composants industriels. Par rapport à un traducteur monoélément, un traducteur multiélément permet de focaliser un faisceau d'ondes ultrasonores en un point déterminé d'une pièce en affectant des retards d'émission appropriés aux différents éléments du traducteur. La multiplication des lois de retard permet d'imager des zones très étendues de la pièce. Une focalisation multipoints peut également être effectuée de façon synthétique par post-traitement des signaux d'acquisition des différents éléments. L'une des méthodes de focalisation synthétique les plus développées à l'heure actuelle est appelée "méthode de reconstruction par focalisation en tous points" (FTP). Cette méthode dérive de la méthode dite SAR pour "Synthetic Aperture Radar", et s'applique notamment sur une méthode d'acquisition dite "matrice inter-élément" ou FMC selon l'expression anglosaxonne "Full Matrice Capture". La méthode d'acquisition FMC consiste à émettre successivement par M émetteurs à ultrasons une onde ultrasonore dans la pièce à imager, et à enregistrer dans N récepteurs ultrasonores les échos provenant de la propagation des ondes émises dans la pièce à inspecter. La méthode aboutit ainsi à la formation de MxN signaux élémentaires. La méthode de reconstruction FTP consiste à sommer de façon cohérente les signaux élémentaires de manière à obtenir des maxima d'amplitude, correspondant à des interférences constructives, pour les points de la pièce où sont effectivement localisés des défauts générant des échos des ondes ultrasonores émises. La cohérence des signaux est déterminée à partir des durées de parcours théoriques des ondes ultrasonores se propageant dans la pièce entre un émetteur et un récepteur donnés et passant par le point considéré. La méthode FTP a été développée initialement pour une focalisation synthétique en mode direct. Le mode direct fait référence à des trajets directs de l'onde ultrasonore entre un émetteur, un point de la zone à imager, et un récepteur. Autrement dit, le mode direct implique des trajets dépourvus de réflexion sur une surface de la pièce. La méthode FTP a ensuite été généralisée à d'autres modes de reconstruction. On parle alors de reconstruction multi-modale, en prenant en compte des trajets d'ondes ultrasonores faisant intervenir une ou plusieurs réflexions sur des surfaces de la pièce à imager. En particulier, il est possible de tenir compte de trajets comportant une réflexion de l'onde ultrasonore sur le fond de la pièce, c'est-à-dire la surface de la pièce opposée à celle insonifiée par le traducteur multiélément émetteur. On parle dans ce cas de mode "écho de coin". La réflexion peut se produire avant ou après une interaction avec le point à reconstruire. Certains trajets peuvent aussi comporter une réflexion avant l'interaction avec le point à reconstruire, et une réflexion après. On parle alors de mode "indirect". La focalisation FTP en mode écho de coin et celle en mode indirect permettent de compléter les informations fournies par la focalisation FTP en mode direct, car elles exploitent des contributions faisant intervenir d'autres phénomènes physiques et issues d'interactions avec des parties différentes de la pièce et du défaut. Elles permettent notamment d'imager d'autres parties du défaut. Par exemple, dans certains cas, la reconstruction en mode écho de coin présente l'avantage d'imager un défaut de la pièce sur toute sa longueur, et non uniquement ses extrémités. En revanche, dans certaines configurations, la focalisation FTP multimode présente l'inconvénient de générer des maxima d'amplitude n'ayant aucune origine physique dans la pièce. Ces maxima, appelés artefacts de reconstruction, sont dus au fait que, pour un point de la zone reconstruction, le temps de vol théorique du trajet calculé pour le mode courant peut être identique, ou sensiblement identique, au temps de vol d'une contribution réelle et correspondant à un mode de propagation différent dans la pièce. Par exemple, la durée de parcours théorique d'une onde en mode écho de coin peut être sensiblement égale à la durée de parcours effective en mode direct et interagissant avec le fond de la pièce. Tel est notamment le cas pour les points à reconstruire situés sous le traducteur multiélément. Pour certaines configurations, la méthode FTP en mode écho de coin peut donc interpréter un écho de fond en mode direct comme étant un écho en mode écho de coin provenant d'un défaut situé sous le capteur. La seule prise en compte des temps de vol ne permet donc pas de distinguer, pour tous les trajets possibles, un point correspondant à un écho sur le fond de la pièce d'un point correspondant à un éventuel défaut de la pièce. Un but de l'invention est notamment de remédier au problème des artefacts de reconstruction en mode écho de coin en fournissant un procédé d'imagerie par ultrasons filtrant ces artefacts. A cet effet, l'invention a pour objet un procédé d'imagerie par ultrasons pour imager une pièce. Le procédé comprend les étapes suivantes : - acquérir un ensemble de M x N signaux ultrasonores Si(t), où M est le nombre d'émetteurs pouvant émettre un signal ultrasonore dans la pièce à imager, N est le nombre de récepteurs pouvant recevoir l'un des signaux ultrasonores S1(t) provenant de la pièce à imager, l'indice i désigne l'un des émetteurs, l'indice j désigne l'un des récepteurs, et t est une variable 20 représentant le temps, - pour chaque point P d'une zone de la pièce à imager et pour chaque couple émetteur-récepteur ij, - déterminer un temps de vol T1(P) correspondant à une durée théorique nécessaire au signal ultrasonore S1(t) pour parcourir un trajet 25 depuis l'émetteur i jusqu'au récepteur] en passant par le point P, - déterminer un angle O1(P) formé entre une portion du trajet rejoignant le point P, et une portion du trajet quittant le point P, - pour chaque point P, déterminer une somme pondérée 1(P) des amplitudes extraites des signaux S1(t) aux temps de vol T1(P) 30 correspondant au point P considéré : N M 1(P) =f (0i1(P)).Si1(Ti1(P)) i=1 1=1 où f (t9ii(P)) est une fonction de pondération déterminée de manière à atténuer les signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de z par rapport aux signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de la valeur nulle. Selon une première forme de réalisation, la fonction de 5 pondération f (t9ii(P)) est une fonction seuil prenant la valeur nulle pour les angles O1(P) compris dans une plage de valeurs [7-t- z +61,'ii], et la valeur un pour les autres angles. Selon une deuxième forme de réalisation, la fonction de 10 pondération f (t9ii(P)) est une fonction affine prenant la valeur nulle pour un angle O1(P) sensiblement égal à ff, et la valeur un pour un angle O1(P) sensiblement égal à la valeur nulle. Selon une troisième forme de réalisation, la fonction de 15 pondération f (o(p)) est définie comme suit : f (61ii(P)) = cos2(61i1(P)/2) Chaque trajet entre un émetteur i et un récepteur] peut être formé de plusieurs portions, chaque portion étant définie entre deux points d'interaction successifs, les points d'interaction étant les émetteurs j, les 20 récepteurs j, les points P, et des surfaces de réflexion de la pièce à imager. Chaque suite de points d'interaction définit un mode. Une somme pondérée 1(P) peut alors être déterminée pour chaque mode. En outre, les temps de vol T1(P) peuvent être déterminés en fonction d'une polarisation du signal ultrasonore dans chaque portion du trajet entre l'émetteur i et le récepteur j 25 considérés. L'invention a également pour objet un dispositif d'imagerie par ultrasons permettant de mettre en oeuvre le procédé selon l'invention. En particulier, le dispositif comprend : 30 - M émetteurs aptes à émettre un signal ultrasonore dans la pièce à imager, - N récepteurs aptes à recevoir l'un des M x N signaux ultrasonores Si(t), où l'indice i désigne l'un des émetteurs, l'indice j désigne l'un des récepteurs, et t est une variable représentant le temps, - des premiers moyens de traitement permettant de déterminer, pour 5 chaque point P d'une zone de la pièce à imager et pour chaque couple émetteur-récepteur ij, un temps de vol T1(P) correspondant à une durée théorique nécessaire au signal ultrasonore S1(t) pour parcourir un trajet depuis l'émetteur i jusqu'au récepteur] en passant par le point P, - des deuxièmes moyens de traitement permettant de déterminer, pour 10 chaque point P d'une zone de la pièce à imager et pour chaque couple émetteur-récepteur ij, un angle O1(P) formé entre une portion du trajet rejoignant le point P, et une portion du trajet quittant le point P, - des troisièmes moyens de traitement permettant de déterminer, pour chaque point P, une somme pondérée 1(P) des amplitudes extraites des 15 signaux S1(t) aux temps de vol T1(P) correspondant au point P considéré : N M 1(P) =f (0i1(P)).Si1(Ti1(P)) i=1 1=1 où f (t9ii(P)) est une fonction de pondération déterminée de manière à atténuer les signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de Ir par rapport aux signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de la valeur nulle. 20 Les émetteurs et/ou les récepteurs peuvent être disposés de différentes manières, par exemple de façon linéaire, matricielle ou annulaire. Dans le cas d'une disposition linéaire ou matricielle, le ou les axes des récepteurs peuvent être parallèles ou perpendiculaires à le ou les axes des 25 émetteurs. Par ailleurs, les émetteurs et/ou les récepteurs peuvent chacun présenter une forme quelconque, notamment une forme annulaire. L'invention présente notamment l'avantage de filtrer les trajets en mode écho de coin ne correspondant pas à une réalité physique, tout en 30 conservant intacte l'image des défauts de la pièce. L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description qui va suivre, faite en regard de dessins annexés sur lesquels : - la figure 1 représente schématiquement un exemple de pièce à 5 imager sur laquelle est disposé un traducteur multiélément ; - la figure 2 illustre un premier mode de reconstruction, appelé mode direct ; - les figures 3A et 3B illustrent un deuxième mode de reconstruction, appelé mode écho de coin ; 10 - la figure 4 illustre un troisième mode de reconstruction, appelé mode indirect ; - la figure 5A représente un exemple de trajet en mode direct avec une polarisation transversale de l'onde ultrasonore et la figure 5B représente une image obtenue par ce mode de reconstruction ; 15 - la figure 6A représente un exemple de trajet en mode écho de coin avec une polarisation transversale de l'onde ultrasonore et la figure 6B représente une image obtenue par ce mode de reconstruction ; - la figure 7 représente, par un logigramme, un exemple de procédé d'imagerie basé sur une méthode de reconstruction par focalisation 20 en tous points sans filtrage des artefacts ; - la figure 8A montre un exemple d'image obtenue par le procédé de la figure 7 comportant des artefacts, et les figures 8B et 8C illustrent l'origine de ces artefacts ; - la figure 9 représente, par un logigramme, un exemple de 25 procédé d'imagerie selon l'invention ; - les figures 10A à 10C représentent des exemples d'images obtenues par le procédé de la figure 9 en mode écho de coin avec une fonction de pondération prenant trois valeurs seuils différentes ; - les figures 11A à 11C représentent des exemples d'images 30 obtenues par le procédé de la figure 7 en mode écho de coin pour différentes polarisations ; - les figures 12A à 12C représentent des images équivalentes aux images 11A à 11C mais obtenues par le procédé de la figure 9 avec une fonction de pondération sinusoïdale. 35 La figure 1 représente schématiquement un exemple de pièce à imager comportant un défaut et sur laquelle est disposé un traducteur multiélément. La pièce 1 comprend une entaille 2 s'étendant depuis l'une de ses surfaces 3. Un traducteur multiélément 4 est disposé sur une surface 5 de la pièce 1 opposée à la surface 3, afin de détecter la présence d'un éventuel défaut dans une zone de reconstruction 6. À titre d'exemple, on considère ici un traducteur linéaire comportant 48 éléments. Chaque élément comprend un émetteur apte à générer une onde ultrasonore à la surface 3 de la pièce 1, et un récepteur apte à délivrer un signal correspondant à un écho d'une onde ultrasonore générée par l'un des émetteurs. Les ondes ultrasonores ont par exemple une fréquence égale à 2 mégahertz. De manière générale, on considère un traducteur multiélément comprenant M émetteurs et N récepteurs, et on note S1(t) les signaux délivrés par les récepteurs, où l'indice i désigne l'un des émetteurs, l'indice j désigne l'un des récepteurs, et la variable t représente le temps. Les émetteurs émettent successivement une onde ultrasonore, de sorte que M x N signaux S1(t) peuvent être générés pour une acquisition donnée. Une telle méthode d'acquisition est appelée méthode de type "matrice inter-élément", ou "Full Matrice Capture" (FMC) en anglais. L'invention repose également sur une méthode de focalisation synthétique appelée "méthode de reconstruction par focalisation en tous points" (FTP). Cette méthode consiste à sommer de manière cohérente les signaux S1(t) de manière à obtenir des maxima d'amplitude aux endroits de la pièce où sont effectivement localisés des défauts. La cohérence des signaux pour un point P donné est déterminée à partir des durées de parcours théoriques des ondes ultrasonores le long de différents trajets entre les émetteurs j et les récepteurs j passant par le point P. Les trajets considérés peuvent être plus ou moins complexes, selon le nombre d'interactions que les ondes ultrasonores subissent sur des surfaces de la pièce. La méthode FTP peut prendre en compte différents types de trajets ; elle est alors qualifiée de multimode. Pour la suite de la description, on appelle points d'interaction les différents points où une onde ultrasonore est susceptible de subir une déviation. Il s'agit des points sur les surfaces de la pièce, ainsi que des points P de la zone de reconstruction où un défaut est 2 9992 92 8 supposé se trouver. Par extension, on considère également les émetteurs et les récepteurs du traducteur comme des points d'interaction. Un trajet d'onde ultrasonore peut ainsi être défini de manière univoque par une suite de points d'interaction. 5 La figure 2 illustre un premier mode de reconstruction, appelé mode direct. Dans ce mode, on considère que chaque onde ultrasonore suit un trajet direct, c'est-à-dire qu'elle est émise depuis l'un des émetteurs vers un point P de la zone de reconstruction, et qu'elle est réémise depuis le point 10 P vers l'un des récepteurs. Le trajet ne subit donc pas de réflexion sur l'une des surfaces de la pièce. Les figures 3A et 3B illustrent un deuxième mode de reconstruction, appelé mode écho de coin. Dans ce mode, on considère les 15 ondes ultrasonores subissant une seule réflexion sur l'une des surfaces de la pièce 1. La figure 3A représente un écho de coin se produisant avant une interaction avec un point P de la zone de reconstruction. Le trajet de l'onde ultrasonore comprend ainsi une première portion 31 entre l'un des émetteurs du traducteur 4 et une surface de la pièce 1, en l'occurrence la surface 3, 20 une deuxième portion 32 entre la surface 3 de la pièce 1 et le point P de la zone de reconstruction, et une troisième portion 33 entre le point P et l'un des récepteurs du traducteur 4. La figure 3B représente un écho de coin se produisant après une interaction de l'onde ultrasonore avec un point P de la zone de reconstruction. Le trajet de l'onde ultrasonore comprend ainsi une 25 première portion 35 entre l'un des émetteurs et le point P, une deuxième portion 36 entre le point P et la surface 3, et une troisième portion 37 entre la surface 3 et l'un des récepteurs. La figure 4 illustre un troisième mode de reconstruction, appelé 30 mode indirect. Dans ce mode, on considère les ondes ultrasonores subissant deux réflexions sur une surface de la pièce. En l'occurrence, sur la figure 4 est représenté un trajet comprenant une première portion 41 entre l'un des émetteurs du traducteur 4 et la surface 3, une deuxième portion 42 entre la surface 3 et un point P de la zone de reconstruction, une troisième portion 43 35 entre le point P et la surface 3, et une quatrième portion 44 entre la surface 3 et l'un des récepteurs. Plus généralement, le mode indirect peut aussi prendre en compte des trajets comprenant une première réflexion sur l'une des surfaces de la pièce, et une deuxième réflexion sur une autre surface de la pièce. En plus des différents modes de reconstruction, la méthode de reconstruction FTP peut prendre en considération la polarisation des ondes ultrasonores, longitudinale ou transversale. Plus précisément, elle peut tenir compte de la polarisation dans les différentes portions du trajet, sachant qu'une onde peut passer d'une polarisation longitudinale à transversale (ou inversement) à chaque point d'interaction par le processus dit de conversion de mode. La polarisation d'une onde influe sur sa vitesse de propagation, et donc sur la durée de parcours entre un émetteur et un récepteur. Les figures 5A et 5B illustrent un mode de reconstruction direct dans lequel l'onde ultrasonore a une polarisation transversale tout le long du trajet, à savoir sur une portion 51 entre un émetteur et un point P, et sur une portion 52 entre le point P et un récepteur. En notant L et T des portions de trajet pour lesquelles l'onde ultrasonore présente une polarisation longitudinale ou transversale, respectivement, et d l'interaction au point P, le mode de reconstruction de la figure 5A peut être appelé "mode direct TdT". La figure 5B représente l'amplitude de la somme des signaux S1(t) pour les différents points P de la zone de reconstruction considérée dans ce mode. Les figures 6A et 6B illustrent un mode de reconstruction écho de coin dans lequel l'onde ultrasonore a une polarisation transversale tout le long du trajet, à savoir sur une portion 61 entre un émetteur et la surface 3, sur une portion 62 entre la surface 3 et un point P, et sur une portion 63 entre le point P et un récepteur. Ce mode de reconstruction est appelé "mode écho de coin TTdT". La figure 6B représente l'amplitude de la somme des signaux S1(t) pour les différents points P de la zone de reconstruction considérée dans ce mode. Les figures 5B et 6B mettent par ailleurs en évidence le fait que le mode direct permet de visualiser uniquement les extrémités d'une entaille dans la pièce, alors que le mode écho de coin permet de visualiser l'entaille sur toute sa longueur. La figure 7 représente, par un logigramme, un exemple de procédé d'imagerie basé sur la méthode de reconstruction FTP. On considère à titre d'exemple la configuration représentée sur la figure 1, dans laquelle la pièce 1 comprend une entaille s'étendant depuis la surface 3, le traducteur multiélément 4 étant disposé sur la surface 5. Dans une étape 71, chacun des M émetteurs, indicés i, du traducteur génère successivement et individuellement une onde ultrasonore dans la pièce 1 depuis l'une de ses surfaces, en l'occurrence depuis la surface 5. Dans une étape 72, chacun des N récepteurs, indicés j, du traducteur reçoit un signal ultrasonore issu de la propagation de l'onde ultrasonore émise lors de l'étape 71 par un émetteur i. Un signal ultrasonore est reçu pour chaque émetteur i, donnant ainsi un signal S1(t) représentatif de l'écho d'onde reçu pour une émission par l'émetteur j et une réception par le récepteur]. Les étapes 71 et 72 sont représentées comme étant réalisées successivement. En réalité, ces étapes sont imbriquées : une nouvelle émission d'une onde ultrasonore par un émetteur i + 1 n'est réalisée qu'une fois que les signaux S1(t) ont été générés par les N récepteurs en réponse à l'onde émise par l'émetteur i. Les étapes 71 et 72 permettent ainsi de générer un ensemble de M x N signaux Sij(t). Ces étapes forment globalement une étape d'acquisition 712. Dans une étape 73, on détermine pour chaque point P de la zone de reconstruction et pour chaque couple émetteur-récepteur jj, une durée théorique nécessaire à l'onde ultrasonore émise par l'émetteur i pour rejoindre le récepteur j en passant par le point P. Cette étape peut être réalisée pour chaque mode de reconstruction (direct, écho de coin ou indirect), et pour chaque type de polarisation des ondes ultrasonores. Les durées théoriques sont nommées temps de vol Tii(P). Dans une étape 74, on détermine pour chaque point P, une somme 1(P) des amplitudes extraites des signaux S1(t) aux temps de vol t = Tii(P). Autrement dit, on somme les amplitudes des M x N signaux aux instants correspondant aux durées respectives nécessaires pour rejoindre l'un des récepteurs depuis l'un des émetteurs en passant par le point P considéré. La somme 1(P) peut être formulée par l'expression suivante : N M 1(P) = (Tii(P)) i=1 1=1 L'étape 74 est réalisée individuellement pour chaque mode de reconstruction et pour chaque type de polarisation considéré. Un exemple d"image obtenue par ce procédé est représenté sur la figure 5B pour le mode direct TdT, et sur la figure 6B pour le mode écho de coin TTdT. La figure 8A montre la présence d'artefacts dans un exemple d'image obtenue par le procédé de la figure 7 en mode écho de coin, et les figures 8B et 8C illustrent l'origine de ces artefacts. Sur la figure 8A, il apparaît que les artefacts se produisent sous la surface du traducteur 4. La figure 8B représente le trajet théoriquement suivi par une onde ultrasonore en mode écho de coin LLdL en considérant un émetteur j, un récepteur] et un point de reconstruction P situé sous le traducteur 4. Ce trajet comprend une première portion 81 entre l'émetteur j et la surface 3, une deuxième portion 82 entre la surface 3 et le point P, et une troisième portion 83 entre le point P et le récepteur]. La figure 8C représente le trajet suivi par une onde ultrasonore en mode direct LdL en considérant le même émetteur j et le même récepteur], et un point de reconstruction P2 situé sur la surface 3 au point d'interaction entre les portions 81 et 82 du trajet de la figure 8B. Ce trajet direct comprend une première portion 84 entre l'émetteur i et le point P2, et une deuxième portion 85 entre le point P2 et le récepteur]. À la lecture des figures 8B et 8C, il est clair que le trajet en mode écho de coin est identique au trajet en mode direct. En particulier, la portion 81 du trajet en mode écho de coin est identique à la portion 84 du mode direct, et la somme des portions 82 et 83 est identique à la portion 85. En conséquence, le temps de vol T1(P) du mode écho de coin est égal au temps de vol Ti1(P2) du mode direct. Lors de la reconstruction du point P en mode écho de coin, un maximum va ainsi se former du fait du temps de vol identique en mode direct. L'invention a pour objet de supprimer la contribution des trajets en mode écho de coin qui correspondent en réalité à des trajets en mode direct. Elle apporte une solution à ce problème en tenant compte de l'angle O1(P) formé entre la portion du trajet rejoignant le point de reconstruction P, et la portion quittant le point P. En l'occurrence, il s'agit de l'angle O1(P) formé entre la portion 82 et la portion 83. En fait, les trajets en mode écho de coin dont l'angle O1(P) est voisin de z n'ont pas de sens physique. Si la pièce 1 présente effectivement un défaut au point P, alors l'onde ultrasonore arrivant à ce point P subit nécessairement une déviation. La figure 9 représente, par un logigramme, un exemple de procédé d'imagerie selon l'invention. Ce procédé se base toujours sur la méthode de reconstruction FTP. Le procédé comprend les étapes 71, 72 et 73 du procédé de la figure 7. Il comprend en outre une étape 91 dans laquelle les angles O1(P) sont déterminés pour chaque point P de la zone à imager et pour chaque couple émetteur-récepteur ij. L'étape 74 du procédé de la figure 7 est substituée par une étape 92 dans laquelle les sommes 1(P) des amplitudes extraites des signaux S1(t) sont pondérées en fonction des angles O1(P) déterminés lors de l'étape 91. La pondération des signaux S1(t) peut être effectuée au moyen d'une fonction de pondération f (t9ii(P)). Cette fonction doit être déterminée de manière à atténuer les signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de z par rapport aux signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de la valeur nulle. La somme 1(P) peut alors être formulée par 20 l'expression suivante : N M 1(P) =f (0i1(P)).Si1(Ti1(P)) i=1 j=1 Il est à noter que les signaux S1(t) peuvent en outre subir différentes transformations. Par exemple, la transformée de Hilbert peut être appliquée à ces signaux. 25 Selon une forme particulière de réalisation, la fonction de pondération f (t9ii(P)) est une fonction seuil prenant la valeur nulle pour les angles O1(P) compris dans une plage de valeurs [7-t- - z + 61,'ii], et la valeur un pour les autres valeurs d'angle. Une telle fonction est qualifiée de seuillage dur. À titre d'exemple, la valeur seuil °seuil peut prendre la valeur 30 10 degrés, 30 degrés ou 60 degrés. Les figures 10A à 10C représentent des exemples d'images obtenues par le procédé de la figure 9 en mode écho de coin LLdL avec une fonction de pondération dont la valeur seuil °seuil prend respectivement les valeurs 10 degrés, 30 degrés et 60 degrés. La détermination de la valeur seuil °seuil permettant de supprimer les artefacts sans altérer la reconstruction des défauts dans la pièce peut s'avérer délicate pour certaines géométries de pièces. Pour cette raison, il est parfois préférable d'utiliser une fonction de pondération variant continument entre la valeur nulle et la valeur un. La fonction de pondération f (O(P)) est par exemple une fonction affine prenant la valeur nulle pour un angle O1(P) sensiblement égal à ff, et la valeur un pour un angle O1(P) sensiblement égal à la valeur nulle. Entre ces deux valeurs, la fonction de pondération f (t9i1(P)) prend une valeur intermédiaire. La fonction de pondération f (t9i1(P)) peut aussi être définie par une fonction sinusoïdale. Elle prend par exemple la forme suivante : f (61ii(P)) = cos2(61i1(P)/2) Les figures 11A à 11C représentent des exemples d'images obtenues par le procédé de la figure 7 en mode écho de coin LLdL, TLdT, et TTdT, respectivement, et les figures 12A à 12C représentent des exemples d'images obtenues par le procédé selon l'invention utilisant la fonction définie ci-dessus dans ces mêmes modes. La fonction sinusoïdale présente l'avantage de ne pas nécessiter de paramètre. L'invention a été décrite en considérant un traducteur multiélément linéaire. Dans un tel traducteur, les émetteurs et les récepteurs sont disposés linéairement. Généralement, l'axe sur lequel s'alignent les émetteurs est parallèle à l'axe sur lequel s'alignent les récepteurs. Les différents éléments du traducteur peuvent être utilisés à la fois en émission et en réception. De manière plus générale, l'invention s'applique à des reconstructions FTP en trois dimensions. Un traducteur à découpe multi-éléments en deux dimensions, par exemple matriciel, annulaire ou sectoriel est alors utilisé. Dans un traducteur matriciel, les éléments (émetteurs et/ou récepteurs) sont disposés le long de lignes et de colonnes. Ils ont généralement une forme rectangulaire ou circulaire. Néanmoins, ils peuvent présenter d'autres formes. En outre, les axes de la matrice d'émetteurs sont généralement parallèles aux axes de la matrice de récepteurs. Par traducteur annulaire, on entend un traducteur dans lequel chaque élément présente une forme annulaire. Les différents éléments forment par exemple des anneaux concentriques. Par traducteur sectoriel, on entend un traducteur dans lequel chaque élément présente une forme déterminée, par exemple une portion d'anneau. De manière générale, les émetteurs et les récepteurs peuvent prendre n'importe quelle forme, et être disposés de manière quelconque les uns par rapport aux autres, dans un plan ou dans l'espace. Par ailleurs, l'invention a été décrite en considérant que le traducteur multiélément est au contact direct de la pièce à imager. Cependant, l'invention s'applique tout aussi bien en utilisant des capteurs en immersion ou en interposant un sabot entre le traducteur multiélément et la pièce à imager. Les temps de vol sont alors déterminés en conséquence. L'invention porte aussi sur un dispositif d'imagerie permettant de mettre en oeuvre le procédé d'imagerie décrit ci-dessus. Le dispositif comprend alors de préférence un traducteur multiélément, et différents moyens permettant de réaliser les étapes 73, 91 et 92 décrites en référence à la figure 9. Ces moyens peuvent être mis en oeuvre sous forme matérielle, logicielle, ou par un mix entre ces deux formes. REVENDICATIONS1. Procédé d'imagerie par ultrasons pour imager une pièce (1), le procédé comprenant les étapes suivantes : - acquérir (712) un ensemble de M x N signaux ultrasonores Si(t), où M est le nombre d'émetteurs pouvant émettre un signal ultrasonore dans la 5 pièce à imager (1), N est le nombre de récepteurs pouvant recevoir l'un des signaux ultrasonores S1(t) provenant de la pièce à imager, l'indice i désigne l'un des émetteurs, l'indice j désigne l'un des récepteurs, et t est une variable représentant le temps, - pour chaque point P d'une zone (6) de la pièce à imager et pour 10 chaque couple émetteur-récepteur ij, - déterminer (73) un temps de vol T1(P) correspondant à une durée théorique nécessaire au signal ultrasonore S1(t) pour parcourir un trajet depuis l'émetteur i jusqu'au récepteur] en passant par le point P, - déterminer (91) un angle O1(P) formé entre une portion du trajet 15 (82) rejoignant le point P, et une portion du trajet (83) quittant le point P, - pour chaque point P, déterminer (92) une somme pondérée 1(P) des amplitudes extraites des signaux S1(t) aux temps de vol T1(P) correspondant au point P considéré : N M 1(P) =f (0i1(P)).Si1(Ti1(P)) i=1 1=1 où f (t9ii(P)) est une fonction de pondération déterminée de manière à 20 atténuer les signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de Ir par rapport aux signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de la valeur nulle. 2. Procédé selon la revendication 1, dans lequel la fonction de 25 pondération f (t9ii(P)) est une fonction seuil prenant la valeur nulle pour les angles O1(P) compris dans une plage de valeurs Ur - + 61,'ii], et la valeur un pour les autres angles. 3. Procédé selon la revendication 1, dans lequel la fonction de pondération f (t9i1(P)) est une fonction affine prenant la valeur nulle pour un angle O1(P) sensiblement égal à ff, et la valeur un pour un angle O1(P) sensiblement égal à la valeur nulle. 4. Procédé selon la revendication 1, dans lequel la fonction de pondération f (t9i1(P)) est définie comme suit : f (61ii(P)) = cos2(61i1(P)/2). 5. Procédé selon l'une des revendications précédentes, dans lequel chaque trajet entre un émetteur i et un récepteur j est formé de plusieurs portions (81, 82, 83), chaque portion étant définie entre deux points d'interaction successifs, les points d'interaction étant les émetteurs j, les récepteurs j, les points P, et des surfaces de réflexion (3, 5) de la pièce à imager (1), chaque suite de points d'interaction définissant un mode, une somme pondérée 1(P) étant déterminée pour chaque mode. 6. Procédé selon la revendication 5, dans lequel les temps de vol T1(P) sont déterminés en fonction d'une polarisation du signal ultrasonore 20 dans chaque portion du trajet (81, 82, 83) entre l'émetteur jet le récepteur] considérés. 7. Dispositif d'imagerie par ultrasons permettant d'imager une pièce, le dispositif comprenant : 25 - M émetteurs (4) aptes à émettre un signal ultrasonore dans la pièce à imager (1), - N récepteurs (4) aptes à recevoir l'un des M x N signaux ultrasonores Si(t), où l'indice j désigne l'un des émetteurs, l'indice j désigne l'un des récepteurs, et t est une variable représentant le temps, 30 - des premiers moyens de traitement permettant de déterminer, pour chaque point P d'une zone de la pièce à imager et pour chaque couple émetteur-récepteur ij, un temps de vol T1(P) correspondant à une durée théorique nécessaire au signal ultrasonore S1(t) pour parcourir un trajet depuis l'émetteur j jusqu'au récepteur] en passant par le point P,- des deuxièmes moyens de traitement permettant de déterminer, pour chaque point P d'une zone de la pièce à imager et pour chaque couple émetteur-récepteur ij, un angle O1(P) formé entre une portion du trajet (82) rejoignant le point P, et une portion du trajet (83) quittant le point P, - des troisièmes moyens de traitement permettant de déterminer, pour chaque point P, une somme pondérée 1(P) des amplitudes extraites des signaux S1(t) aux temps de vol T1(P) correspondant au point P considéré : N M 1(P) =f (0i1(P)).Si1 (Tii(P)) i=1 1=1 où f (t9ii(P)) est une fonction de pondération déterminée de manière à atténuer les signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) 10 est voisin de z par rapport aux signaux sij (Tii(P)) correspondant aux trajets dont l'angle O1(P) est voisin de la valeur nulle. 8. Dispositif selon la revendication 7, dans lequel les émetteurs sont disposés linéairement selon un premier axe. 15 9. Dispositif selon la revendication 7, dans lequel les émetteurs sont disposés sous forme matricielle selon un premier et un deuxième axe. 10. Dispositif selon l'une des revendications 7 à 9, dans lequel les 20 récepteurs sont disposés linéairement selon un troisième axe. 11. Dispositif selon l'une des revendications 7 à 9, dans lequel les récepteurs sont disposés sous forme matricielle selon un troisième et un quatrième axe. 25 12. Dispositif selon l'une des revendications 10 et 11, prise avec l'une des revendications 8 et 9, dans lequel le premier et le troisième axe sont parallèles. 30 13. Dispositif selon les revendications 9, 11 et 12, dans lequel le deuxième et le quatrième axe sont parallèles. 14. Dispositif selon la revendication 7, dans lequel les émetteurs et/ou les récepteurs présentent une forme annulaire.