FR 3008801 A1 20150123 FR 1357036 A 20130717 PROCÉDÉ ET DISPOSITIF D'IMAGERIE PAR ULTRASONS AVEC PREDICTION DES ARTEFACTS INDUITS ENTRE MODES DE RECONSTRUCTION DESCRIPTION DOMAINE TECHNIQUE Le domaine de l'invention est celui de l'imagerie par ultrasons. L'invention s'applique en particulier au contrôle non destructif de pièces ou d'installations industrielles pour y détecter la présence d'éventuels défauts, 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 les technique de focalisations synthétiques. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Dans le domaine du contrôle non destructif par ultrasons, l'une des méthodes de focalisation synthétique les plus développées, appelée "méthode de reconstruction par focalisation en tous points" (FTP), consiste en une focalisation synthétique des champs émis et reçus par sommation cohérente des données issues d'une acquisition du type Full Matrix Capture (FMC) ou matrice inter-éléments. 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. La méthode aboutit ainsi à la formation de M*N signaux élémentaires Su(t) avec 1 Pour les trajets de type échos de coin et échos indirects, et comme illustré par les figures lb et lc, les trajets pris en compte pour le calcul des temps de vol incluent une ou plusieurs réflexions sur le fond au cours du trajet entre l'émetteur et le point de la zone de reconstruction et/ou au cours du trajet entre le point de reconstruction et le récepteur. Plus précisément, le mode écho de coin prend en considération une réflexion de l'onde ultrasonore sur le fond de la pièce, avant ou après une interaction avec le point à reconstruire, tandis que le mode écho indirect prend en compte à la fois une réflexion avant l'interaction avec le point à reconstruire et une réflexion après. La focalisation FTP en mode écho de coin et celle en mode indirect permettent typiquement 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. Outre le type de trajet, un mode de reconstruction est également caractérisé par le mode de propagation du signal ultrasonore (longitudinal, transversal) sur chaque portion du trajet. On utilise conventionnellement les lettres L et T pour désigner respectivement une propagation longitudinale, respectivement transversale. Dans le cas des trajets directs, la première lettre correspond au trajet entre l'émission et le lieu de reconstruction dans la pièce, la seconde au trajet entre ce lieu et la réception. Dans le cas des trajets en écho de coin, la première lettre représente le trajet entre l'émission et le fond de pièce, la deuxième le trajet entre le fond de pièce et le lieu de reconstruction, la troisième, le trajet entre ce lieu et la réception. Un exemple de résultat de reconstruction par FTP à partir de données expérimentales est présenté sur les figures 2a et 2b. Ce résultat montre que l'exploitation combinée de différents modes de reconstruction présente un intérêt pour la caractérisation des défauts. En effet, dans le cas d'une reconstruction selon le mode direct LL, les deux échos de diffraction (entourés en pointillés sur l'image à droite de la figure 2a) sont correctement repositionnés aux deux extrémités, hautes et basses, de l'entaille, tandis que pour une reconstruction en mode écho de coin TTL, ce dernier est correctement repositionné sur toute la longueur de l'entaille (cf. image à droite de la figure 2b). Les deux échos de diffraction permettent alors de déterminer la hauteur du défaut et l'écho de coin permet de déterminer son orientation. Cependant, les reconstructions FTP multi-modes réalisées sur différentes géométries ont mis en avant la présence d'indications n'ayant aucune origine physique dans le composant inspecté. L'analyse de telles indications, aussi appelées artefacts de reconstruction, montre que leur présence est due à des temps de parcours équivalents pour des modes de reconstruction différents. Les échos de défauts ou de géométrie détectés par certains modes peuvent être reconstruits pour d'autres modes dès lors qu'il y a coïncidence des temps de vols. Un exemple de ce genre d'artefact provenant de l'écho de géométrie est présenté sur les figures 3a-3b. Cet exemple illustre le fait qu'une indication réellement détectée dans une reconstruction en mode direct LL engendre une contribution dans la reconstruction en mode direct LT dans les zones où les temps de vol coïncident. Pour la reconstruction obtenue en mode LL (figure 3a), l'écho de géométrie EF, correspondant à la réflexion de l'onde ultrasonore sur le fond de la pièce, est correctement repositionnée dans la pièce. Pour la reconstruction en mode LT (figure 3b), on constate la présence d'artefacts A de très forte amplitude sous l'ouverture géométrique du capteur. Ces artefacts de reconstruction proviennent du fait que, pour un couple émetteur-récepteur, pour un point de reconstruction P* suivant le mode LT, le temps de vol correspondant est identique à celui du point P (toujours pour le même couple émetteur-récepteur) reconstruit à une position différente dans la pièce suivant le mode LL. Il en résulte que l'amplitude attribuée au point P lors de la reconstruction LL est attribuée au point P* lors de la reconstruction LT. Dans cet exemple, si l'amplitude de P est lié à l'écho provenant du fond de la pièce, celle-ci sera attribuée de manière impropre au point P* induisant alors l'apparition d'un artefact dans le mode LT. Comme le montre cet exemple, le risque de confusion entre une indication associée à un défaut et un artefact de reconstruction est non négligeable. Néanmoins, si la présence des artefacts dans les reconstructions obtenues par des techniques de focalisation synthétiques ou des techniques de focalisation conventionnelles est un fait bien connu, à ce jour, les procédures d'identification de tels artefacts sont complexes et reposent sur la connaissance de la géométrie de la pièce et du défaut. Elles sont en outre limitées à l'identification d'artefacts provenant de l'écho de géométrie. Il existe ainsi un besoin pour des outils qui permettent d'analyser les artefacts présents dans les reconstructions par focalisation synthétique multi-modes, en permettant notamment de distinguer des échos dus à des défauts de ceux dus à la géométrie de la pièce, et de filtrer les échos indésirables, tout en restant simples à mettre en oeuvre. EXPOSÉ DE L'INVENTION L'invention vise à répondre à ces besoins et propose pour ce faire, selon un premier aspect, un procédé d'imagerie par ultrasons pour imager une pièce, caractérisé par la mise en oeuvre des étapes suivantes : sélection d'une première sous-région de la pièce à partir d'une première image d'une région de la pièce, détermination, pour chaque point de la première sous-région sélectionnée, des temps de vol correspondant aux parcours selon un premier mode de reconstruction passant par le point depuis un émetteur i jusqu'à un récepteur j pour un ensemble de M*N couples émetteur-récepteur d'un signal ultrasonore; détermination d'une deuxième sous-région de la pièce, un point de la région appartenant à la deuxième sous-région lorsqu'un temps de vol du parcours selon un deuxième mode de reconstruction passant par le point depuis un émetteur i jusqu'à un récepteur j dudit ensemble de M*N couples émetteur-récepteur coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée. Certains aspects préférés mais non limitatifs de ce procédé sont les suivants : - il comprend une reconstruction d'une deuxième image de la région de la pièce par une focalisation synthétique suivant le deuxième mode de reconstruction selon laquelle, pour chaque point de la région, et pour chaque couple émetteur-récepteur i-j, on exploite un signal ultrasonore considéré à un temps de vol du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point de la région, uniquement si ledit temps de vol selon le deuxième mode de reconstruction coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée ; - il comprend une reconstruction d'une deuxième image de la région de la pièce par une focalisation synthétique suivant le deuxième mode de reconstruction selon laquelle, pour chaque point de la région, et pour chaque couple émetteur-récepteur i-j, on exploite un signal ultrasonore considéré à un temps de vol du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point de la région, uniquement si ledit temps de vol selon le deuxième mode de reconstruction diffère des temps de vol des parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par les points de la première sous-région sélectionnée. - un mode de reconstruction est caractérisé par un type de trajet comportant plusieurs portions et un mode de propagation du signal ultrasonore sur chaque portion du trajet ; - le premier mode de reconstruction est un mode selon lequel chaque signal ultrasonore suit un trajet direct à deux portions, le mode de propagation du signal étant latéral sur chacune des portions ; - le deuxième mode de reconstruction est un mode selon lequel chaque signal ultrasonore suit un trajet direct à deux portions et le mode de propagation du signal est différent sur chacune des portions, ou un mode selon lequel chaque signal ultrasonore suit un trajet indirect à trois portions au moins. BRÈVE DESCRIPTION DES DESSINS D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en référence aux dessins annexés sur lesquels, outres les figures la-1c, 2a-2b et 3a-3b déjà discutées ci-dessus : - les figures 4a et 4b illustrent un premier mode de reconstruction de type direct LL et un deuxième mode de reconstruction de type direct LT ; la figure 5 illustre différentes indications présentes dans l'image d'une pièce reconstruite par un premier mode de reconstruction de type direct LL; les figures 6a-6c illustrent l'application d'un procédé conforme à l'invention avec un premier mode de reconstruction de type direct LL et un deuxième mode de reconstruction de type direct LT ; les figures 7a-7c illustrent le filtrage des artefacts dans une reconstruction selon un mode de type direct LT ; les figures 8a-8c et 9a-9c illustrent l'application d'un procédé conforme à l'invention avec un premier mode de reconstruction de type direct LL et un deuxième mode de reconstruction de type écho de coin LLL, avec sélection de deux régions différentes dans une image reconstruite en considérant un premier mode de reconstruction de type direct LL; les figures 10a-10d illustrent le filtrage des artefacts dans une reconstruction selon un mode de type écho de coin LLL. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention porte sur un procédé d'imagerie par ultrasons pour imager une pièce, mettant en oeuvre une reconstruction par focalisation synthétique, par exemple une focalisation en tous points. En référence à la figure 4a, le procédé exploite un traducteur multiélément 1 disposé sur la surface d'une pièce 2 à imager, comportant M émetteurs d'un signal ultrasonore, indicés i, et N récepteurs d'un signal ultrasonore, indicés j, fournissant des signaux Su(t) représentatifs de l'écho de l'onde ultrasonore reçu pour une émission par un émetteur i et une réception par un récepteur j. Toujours en référence à la figure 4a, le procédé comporte une première étape de sélection d'une première sous-région Z de la pièce 2 à partir d'une première image IA(Z) d'une région Z (ic Z) de la pièce 2. Cette sélection peut être réalisée manuellement par un opérateur qui visualise et sélectionne une sous-région qui contient un écho dont l'origine est identifiable (écho de fond de pièce, ou écho de défaut par exemple). En variante, la sélection de la sous-région peut être automatisée, par exemple au moyen d'outils d'analyse automatique d'images aptes à dessiner les contours de sous-régions présentant les plus fortes amplitudes (typiquement par seuillage de l'image et conservation des points dont les amplitudes sont supérieures au seuil). Le procédé comporte ensuite une étape de détermination, pour chaque point P de la première sous-région sélectionnée Z, des temps de vol des parcours selon un premier mode de reconstruction passant par le point P depuis un émetteur i jusqu'à un récepteur j pour l'ensemble de M*N couples émetteur-récepteur. En d'autres termes, on vient extraire les temps de vol selon le premier mode de reconstruction associés à chaque point de la sous-région sélectionné pour un mode de reconstruction donné. On forme alors, pour chacun des M*N couples émetteur-récepteur, un sous-ensemble de temps de vol TuA (Z) comprenant l'ensemble des temps de vol selon le premier mode de reconstruction (mode A) depuis l'émetteur i jusqu'à un récepteur j passant par l'un des points de la sous-région sélectionnée Z . Dans un mode de réalisation possible de l'invention, le procédé comporte une étape de reconstruction de ladite première image IA(Z) de la région Z par une focalisation synthétique suivant le premier mode de reconstruction d'ondes ultrasonores (mode A) exploitant, pour chaque point P de la région, les M*N signaux ultrasonores Su(t) considérés, pour chaque couple émetteur-récepteur i-j, à un temps de vol TuA(P) correspondant au parcours selon ledit premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point P de la région. L'ensemble de ces temps de vol, calculés pour tous les couples émetteur-récepteur i-j et pour tous les points P de la région reconstruite Z, est dénoté TA (Z) =A (Z)1 Dans le cadre de J m " ce mode de réalisation les différents sous-ensembles de temps de vol TuA(i) sont alors connus de la reconstruction de la première image suivant le premier mode A. Le procédé comporte par ailleurs une étape de détermination d'une deuxième sous-région de la pièce, un point P de la région Z appartenant à la deuxième sous-région lorsqu'un temps de vol TuB(P) du parcours selon un deuxième mode de reconstruction (mode B) passant par le point P depuis un émetteur i jusqu'à un récepteur j dudit ensemble de M*N couples émetteur-récepteur coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction passant par l'un des points de la première sous-région sélectionnée Z depuis l'émetteur i jusqu'au récepteur j (T 1,7 (13)e TuA ) - A titre d'exemple, les échos présents dans la sous-région sélectionnée de la figure 4a induisent des artefacts dans une image reconstruite selon le mode B dans la zone entourée en pointillé sur la figure 4b qui comprend les points pour lesquels le temps de vol en mode B correspond à un temps de vol en mode A passant par un point de la sous-région sélectionnée i . On comprend donc que l'invention réalise une transformation en mode B de la sous-région sélectionnée Z . Dans un mode de réalisation possible de l'invention, le procédé comporte par ailleurs une étape de reconstruction d'une deuxième image IB(Z) de la région Z de la pièce par une focalisation synthétique suivant le deuxième mode de reconstruction (mode B). Mais à la différence d'une reconstruction conventionnelle exploitant l'ensemble des temps de vol pour tous les couples émetteur-récepteur et tous les points de la région, la reconstruction de la deuxième image selon l'invention met un oeuvre une isolation, ou a contrario un filtrage, d'artefacts de reconstruction. Le procédé peut ainsi comprendre une reconstruction d'une deuxième image IB(Z) de la région Z de la pièce suivant le deuxième mode de reconstruction (mode B) exploitant, pour chaque point P de la région Z, et pour chaque couple émetteur- récepteur i-j, un signal ultrasonore Su(t) considéré à un temps de vol TuB(P) du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point P de la région, uniquement si ledit temps de vol TuB(P) coïncide (isolation d'artefacts), respectivement diffère (filtrage d'artefacts), avec un temps de vol d'un parcours selon le premier mode de reconstruction passant par l'un des points de la sous-région sélectionnée depuis l'émetteur i jusqu'au récepteur j. En d'autres termes, la reconstruction de la deuxième image suivant le mode B exploite, lors d'une isolation d'artefacts, et pour un couple émetteur-récepteur i-j donné, uniquement les temps de vol préalablement extraits formant le sous-ensemble A ZT ( ). A contrario, la reconstruction de la deuxième image suivant le mode B exploite, lors d'un filtrage d'artefacts, et pour un couple émetteur-récepteur i-j donné, uniquement les temps de vol n'appartement pas au sous-ensemble préalablement extrait TuA . Dans le cas d'une isolation, les indications présentes dans la sous-région sélectionnée Z à partir de la première image IA(Z) ont la même origine que celles présentes dans la reconstruction de la deuxième image. Dans l'exemple des figures 6a et 6b, il est alors possible de déduire que les indications liées à l'écho de fond de pièce et à l'écho de défaut détecté (zonez) ) suivant le premier mode A induisent les artefacts entourées en pointillé dans la deuxième image IB(Z). Pour réaliser l'isolation d'artefacts de reconstruction dans la deuxième image, on détermine, pour chaque point P de la région Z, une somme d'amplitudes des signaux acoustiques limitée aux signaux acoustiques reçus à un récepteur j depuis un émetteur i lorsque le temps de vol TuB(P) du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point P coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par l'un des points de la sous-région sélectionnée ( 7;7 (P)e TjA (Z) ). Ainsi, lors de la reconstruction suivant le mode B, on ne considère que les contributions Su(T,,B(P)) des points pour lesquels les temps de vols coïncident en mode B avec ceux de la sous-région sélectionnée dans le mode A. Pour l'isolation, la somme d'amplitudes des signaux acoustiques M N déterminées en chaque point P peut ainsi s'exprimer selon /B(P)=ELAB(P) où i=1 j =1 ^ S (TB (P)) Si { l'amplitude A(P) est donnée par AijB(P)= ''' ''' 0 sin on Pour réaliser le filtrage d'artefacts de reconstruction dans la deuxième image, on détermine, pour chaque point P de la région Z, une somme d'amplitudes des signaux acoustiques limitée aux signaux acoustiques reçus à un récepteur j depuis un émetteur i lorsque le temps de vol TijB( P ) du parcours selon le deuxième mode de TB (P)E TY A (Z) Y - reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point P ne coïncide pas avec les temps de vol des parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par les points de la sous-région sélectionnée ( T,,B (P)e TiA (i) ). Ainsi, lors de la reconstruction suivant le mode B, on ne considère pas les contributions Su(T,,B(P)) des points pour lesquels les temps de vols coïncident en mode B avec ceux de la sous-région sélectionnée dans le mode A. Pour le filtrage, la somme d'amplitudes des signaux acoustiques M N déterminées en chaque point P peut ainsi s'exprimer selon /B(P)=ELAB(P)' où l'amplitude A (P) est donnée par A (13)= sy y (TB (ID)) 0 i=1 j =1 { si TB (13)e T A (i) Y Y sin on Compte tenu de la différence de célérité existant entre les ondes L et T, l'analyse des images est préférentiellement réalisée en considérant tout d'abord les modes de propagation les plus rapides (direct LL par exemple) puis les plus lents (direct LT et écho de coin LLL par exemple). En effet, dans le cas du mode direct LL, les risques de présence d'artefacts de reconstruction provenant d'un mode de propagation plus lent sont réduits. En particulier, pour la sous-région considérée, ce mode ne pourra pas contenir d'indication liée à un écho de géométrie détecté suivant les modes direct LT ou écho de coin LLL, ce qui rend la nature des indications observables plus facilement interprétable. La figure 5 donne à cet égard un exemple d'image reconstruite selon le mode direct LL sur laquelle on observe la présence d'un écho de fond de pièce EF à l'aplomb du traducteur, d'un écho de géométrie EB lié au bourrelet de soudure, d'indications liées à l'écho d'émission EE et un écho ED lié à l'entaille débouchante. Le procédé précédemment décrit peut être appliquée sur ces différentes indications afin de prédire les artefacts qui leur seront associés dans d'autres modes de reconstruction. Les figures 6a-6c donnent un exemple de mise en oeuvre du procédé lors d'une isolation d'artefacts dans une reconstruction selon le mode direct LT. La figure 6a est une image reconstruite selon le mode LL (mode A) dans laquelle des échos sont observables dans la sous-région sélectionnée i . Ces échos induisent des artefacts dans une image reconstruite selon le mode LT (mode B) dans la zone entourée en pointillés sur la figure 6b qui comprend les points pour lesquels le temps de vol du parcours en mode LT correspond à un temps de vol du parcours en mode LL passant par un point de la sous- région sélectionnée Z (transformation en mode B de la sous-région i ). La figure 6c illustre une reconstruction en mode LT standard. La comparaison des figures 6b et 6c permet alors d'interpréter la nature d'un grand nombre d'indications, en particulier la présence d'artefacts A. Les figures 7a-7c illustrent quant à elle un exemple de mise en oeuvre du procédé lors d'un filtrage d'artefacts dans une reconstruction selon le mode LT. La figure 7a correspond à la figure 6c et illustre une reconstruction en mode LT standard. La figure 7b correspond à la figure 6b et illustre une transformation en mode B de la sous-région sélectionnée i . La figure 7c illustre quant à elle le résultat du filtrage conforme à l'invention, où la deuxième image de la région est débarrassée des artefacts A de la figure 6c induits par les échos observés dans la sous-région Z de la première image. Les figures 8a-8c, respectivement 9a-9c, donnent un exemple de mise en oeuvre du procédé lors d'une isolation d'artefacts dans une reconstruction selon le mode écho de coin LLL. La figure 8a, respectivement 9a, est une image reconstruite selon le mode direct LL (mode A) dans laquelle des échos sont observables dans une sous- région sélectionnée Z 1, respectivement Z 2. Ces échos induisent des artefacts dans une image reconstruite selon le mode LLL (mode B) dans la zone entourée en pointillés sur la figure 8b, respectivement 9b, qui correspond à la transformation en mode B de la sous- région i . La figure 8c, respectivement 9c, illustre une reconstruction en mode LLL standard. La comparaison des figures 8b et 8c, respectivement 9b et 9c, permet alors d'interpréter la nature d'un grand nombre d'indications, en particulier la présence d'artefacts. Les figures 10a-10d illustrent quant à elles un exemple de mise en oeuvre du procédé lors d'un filtrage d'artefacts dans une reconstruction selon le mode LLL. La figure 10a correspond aux figures 8c et 9c et illustre une reconstruction en mode LLL standard. La figure 10b correspond à la figure 8b et illustre transformation en mode B de la sous-région Z 2. La figure 10c correspond à la figure 9b et illustre transformation en mode B de la sous-région Z 1. La figure 10d illustre quant à elle le résultat du filtrage conforme à l'invention, où la deuxième image de la région est débarrassée des artefacts induits par les échos observés dans les sous-régions Z 1 et Z 2 de la première image. L'invention n'est pas limitée au procédé précédemment décrit, mais s'étend également à un dispositif d'imagerie par ultrasons pour imager une pièce 1 apte à mettre en oeuvre ce procédé, notamment un dispositif comprenant M émetteurs d'un signal ultrasonore dans la pièce à imager 1 et N récepteurs d'un signal ultrasonores Su(t) provenant de la pièce à imager, l'indice i désignant l'un des émetteurs, l'indice j désignant l'un des récepteurs, et t représentant le temps , caractérisé en ce qu'il comporte des moyens de traitement configurés pour : - sélectionner une première sous-région de la pièce à partir d'une première image d'une région de la pièce 1, déterminer, pour chaque point P de la première sous-région sélectionnée, des temps de vol des parcours selon un premier mode de reconstruction passant par le point P depuis un émetteur i jusqu'à un récepteur j pour un ensemble de M*N couples émetteur-récepteur d'un signal ultrasonore; déterminer une deuxième sous-région de la pièce, un point P* de la région appartenant à la deuxième sous-région lorsqu'un temps de vol du parcours selon un deuxième mode de reconstruction passant par le point P* depuis un émetteur i jusqu'à un récepteur j dudit ensemble de M*N couples émetteur-récepteur coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée. Et l'invention s'étend également à un produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé précédemment décrit, lorsque ledit programme est exécuté sur un ordinateur. On aura compris de la présente description que l'invention permet d'analyser l'origine d'indications présentes dans des reconstructions par focalisation synthétique multimodale et ainsi de distinguer des échos dus à des défauts de ceux dus à la géométrie de la pièce. L'invention présente en outre l'avantage de s'appliquer à toutes les reconstructions multimodales. Elle utilise en outre les mêmes paramètres qu'une focalisation synthétique standard, et donc, ne nécessite pas de calculs ou de mesures supplémentaires. Elle peut de plus être utilisée non seulement pour une identification des artefacts, que pour un filtrage de ces derniers. REVENDICATIONS1. Procédé d'imagerie par ultrasons pour imager une pièce (1), caractérisé par la mise en oeuvre des étapes suivantes : sélection d'une première sous-région (Z) de la pièce à partir d'une première image (IA(Z)) d'une région (Z) de la pièce (1), détermination, pour chaque point de la première sous-région sélectionnée (i ), des temps de vol ( TuA(i)) correspondant aux parcours selon un premier mode de reconstruction (A) passant par le point depuis un émetteur i jusqu'à un récepteur j pour un ensemble de M*N couples émetteur-récepteur d'un signal ultrasonore; détermination d'une deuxième sous-région de la pièce, un point (P) de la région appartenant à la deuxième sous-région lorsqu'un temps de vol (T,,B(P)) du parcours selon un deuxième mode de reconstruction (B) passant par le point (P) depuis un émetteur i jusqu'à un récepteur j dudit ensemble de M*N couples émetteur- récepteur coïncide avec un temps de vol ( TA(i)) d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée. 2. Procédé selon la revendication 1, comprenant une reconstruction d'une deuxième image (IB(Z)) de la région de la pièce par une focalisation synthétique suivant le deuxième mode de reconstruction selon laquelle, pour chaque point (P) de la région, et pour chaque couple émetteur-récepteur i-j, on exploite un signal ultrasonore Su(t) considéré à un temps de vol (T,,B(P)) du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point (P) de la région, uniquement si ledit temps de vol selon le deuxième mode de reconstruction coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée. 3. Procédé selon la revendication 2, dans lequel la reconstruction de la deuxième image (IB(Z)) comprend, pour chaque point (P) de la région, la détermination d'une somme d'amplitudes de signaux acoustiques limitée aux signaux acoustiques reçus à un récepteur j depuis un émetteur i lorsque le temps de vol du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point (P) coïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par l'un des points de la première sous-région sélectionnée. 4. Procédé selon la revendication 1, comprenant une reconstruction d'une deuxième image (IB(Z)) de la région de la pièce par une focalisation synthétique suivant le deuxième mode de reconstruction selon laquelle, pour chaque point (P) de la région, et pour chaque couple émetteur-récepteur i-j, on exploite un signal ultrasonore Su(t) considéré à un temps de vol (T,,B(P)) du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point de la région, uniquement si ledit temps de vol selon le deuxième mode de reconstruction diffère des temps de vol des parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par les points de la première sous-région sélectionnée. 5. Procédé selon la revendication 4, dans lequel la reconstruction de la deuxième image (IB(Z)) comprend, pour chaque point de la région, la détermination d'une somme d'amplitudes des signaux acoustiques limitée aux signaux acoustiques reçus à un récepteur j depuis un émetteur i lorsque le temps de vol du parcours selon le deuxième mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par le point ne coïncide pas avec l'un des temps de vol des parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j en passant par les points de la première sous-région sélectionnée. 6. Procédé selon l'une des revendications 1 à 5, dans lequel un mode de reconstruction est caractérisé par un type de trajet comportant plusieurs portions et un mode de propagation du signal ultrasonore sur chaque portion du trajet. 7. Procédé selon la revendication 6, dans lequel le premier mode de reconstruction est un mode selon lequel chaque signal ultrasonore suit un trajet direct à deux portions, le mode de propagation du signal étant latéral sur chacune des portions. 8. Procédé selon la revendication 7, dans lequel le deuxième mode de reconstruction est un mode selon lequel chaque signal ultrasonore suit un trajet direct à deux portions et le mode de propagation du signal est différent sur chacune des portions, ou un mode selon lequel chaque signal ultrasonore suit un trajet indirect à trois portions au moins. 9. Dispositif d'imagerie par ultrasons pour imager une pièce (1) mettant en oeuvre une reconstruction par focalisation synthétique, comprenant M émetteurs d'un signal ultrasonore dans la pièce à imager (1) et N récepteurs d'un signal ultrasonores Su(t) provenant de la pièce à imager, l'indice i désignant l'un des émetteurs, l'indice j désignant l'un des récepteurs, et t représentant le temps , caractérisé en ce qu'il comporte des moyens de traitement configurés pour : - sélectionner une première sous-région de la pièce à partir d'une première image d'une région de la pièce (1), déterminer, pour chaque point de la première sous-région sélectionnée, des temps de vol correspondant aux parcours selon un premier mode de reconstruction passant par le point depuis un émetteur i jusqu'à un récepteur j pour un ensemble de M*N couples émetteur-récepteur d'un signal ultrasonore; déterminer une deuxième sous-région de la pièce, un point (P) de la région appartenant à la deuxième sous-région lorsqu'un temps de vol du parcours selon un deuxième mode de reconstruction passant par le point depuis (P) un émetteur i jusqu'à un récepteur j dudit ensemble de M*N couples émetteur-récepteurcoïncide avec un temps de vol d'un parcours selon le premier mode de reconstruction depuis l'émetteur i jusqu'au récepteur j passant par l'un des points de la première sous-région sélectionnée. 10. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 8 lorsque ledit programme est exécuté sur un ordinateur.