Brevet FR3085096A1 — publié le 20200221

Procédé et dispositif d'imagerie amélioré pour imager un objet L'invention concerne le domaine de l'imagerie par ultrasons ou plus généralement par ondes élastiques. Elle s'applique en particulier au contrôle 5 non destructif d'objets ou 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 amélioré à sous-échantillonnage adaptatif qui permet de diminuer la complexité algorithmique, autrement dit le nombre d'opération et 10 le temps de calcul nécessaires pour obtenir une image par rapport à un procédé d'imagerie par ondes élastiques conventionnel. Le contrôle non destructif de structures ou de matériaux présente de nombreuses applications notamment, mais pas exclusivement, dans les domaines de l'industrie pétrolière, de l'industrie navale, de l'aéronautique, de 15 l'industrie automobile ou encore de la sidérurgie ou de l'industrie de l'énergie. Dans le domaine du contrôle non destructif d'objets, les techniques d'imagerie par capteurs ultrasons multiéléments sont couramment utilisées. En particulier, il existe une méthode d'imagerie dite « Total Focusing zo Method » ou TFM qui fournit une image d'une zone ou d'un volume d'une pièce à partir d'une collection de signaux ultrasonores obtenus par acquisition à l'aide d'un capteur ultrasons multiéléments. L'acquisition est par exemple du type « Full Matrix Capture » ou FMC. L'imagerie peut être obtenue en deux dimensions ou en trois dimensions selon que le capteur est 25 linéaire ou matriciel. Lors d'une acquisition du type FMC, chaque élément constitutif du capteur ou transducteur multiélément émet successivement une onde ultrasonore, tandis que tous les éléments sont actifs en réception. Pour un transducteur constitué de N éléments, N*N signaux sont donc acquis. 30 L'imagerie TFM consiste à afficher en chacun des pixels de la zone imagée la somme des N2 signaux retardés du temps de vol théorique associé à la 2 position du pixel dans la zone. Ce faisant, on obtient une image présentant des tâches de forte intensité la où sont présents les réflecteurs responsables des échos reçus. Ces réflecteurs sont le plus souvent associés a des discontinuités géométriques correspondant à des défauts ou encore des bords de pièce par exemple. Cette méthode permet donc d'imager précisément l'intérieur d'une structure, en mettant en exergue les réflecteurs des ondes ultrasonores (correspondants à des défauts tels que des trous ou des fissures...). Pour créer une image comprenant un nombre P de pixels avec la méthode d'imagerie TFM, il faut donc calculer les P x N2 temps de vol théoriques associés aux P pixels et aux N2 paires d'éléments (ou N2 signaux). Puis il faut calculer la somme des N2 signaux retardés de ces temps de vol théoriques, pour obtenir une valeur scalaire associée à chaque pixel représentant l'amplitude de réflexion des ondes par ce pixel. La juxtaposition de ces valeurs scalaires aux P pixels fournit l'image TFM permettant par exemple de diagnostiquer l'état de la structure. Le temps de calcul de l'image est donc proportionnel au coefficient P"N2. Ce temps de calcul limite la résolution de la grille (ou bien la taille de l'image) compatible avec la contrainte de calcul de l'image en temps réel, en particulier dans des appareils portables, qui sont naturellement limités en termes de masse et taille, utilisés pour réaliser des inspections sur site. Même si des appareils portables intégrant le procédé d'imagerie TFM en temps réel existent, ceux-ci sont onéreux et massifs notamment en raison du système de calcul nécessaire à l'imagerie. Cette limitation est aujourd'hui un frein au déploiement de la technique d'imagerie TFM. Ainsi, un problème à résoudre consiste à diminuer le nombre de calculs et/ou le temps de calcul nécessaire pour obtenir une image à une résolution donnée dans l'optique d'exploiter cette image pour détecter et caractériser des défauts dans une pièce à inspecter et ce, en temps réel. Un autre objectif de l'invention est de fournir une méthode d'imagerie dont les 3 calculs sont parallélisables afin de ne pas limiter la vitesse de calcul d'une image. Des méthodes permettant une mise en oeuvre temps réelle d'imagerie 5 dite TFM sont notamment décrites dans les références [1],[2],[3] et [4]. Les références [1] et [2] décrivent des implémentations optimisées pour des processeurs graphiques de type GPU (Graphics Processing Unit). Les références [3] et [4] décrivent des implémentations pour des circuits logiques programmables de type FPGA (Field Programmable Gate Array). 10 Toutes ces méthodes requièrent une puissance de calcul conséquente et dédiée afin de pouvoir effectuer une imagerie en temps réel. De ce fait, elles imposent un coût important de mise en oeuvre du dispositif de calcul de l'image. 15 L'invention propose une méthode d'imagerie mettant en oeuvre un sous-échantillonnage adaptatif et optimisé afin d'atteindre une haute résolution dans certaines zones d'intérêts de l'image qui correspondent aux zones de l'objet à imager dans lesquelles des réflecteurs sont présents, tout en limitant la résolution de l'image dans des zones qui présentent moins 20 d'intérêt du point de vue de l'inspection de l'objet. Ainsi, l'invention permet de limiter la puissance de calcul nécessaire à l'acquisition et le calcul d'une image pour diminuer le coût du dispositif de calcul associé. 25 L'invention a ainsi pour objet un procédé d'imagerie pour imager un objet à l'aide d'un dispositif d'imagerie comprenant une pluralité d'émetteurs aptes à émettre une onde élastique dans l'objet à imager et une pluralité de récepteurs aptes à recevoir une onde élastique provenant de l'objet à imager, ledit procédé comprenant les étapes de : 4 Déterminer une première image de l'objet, à un premier niveau d'échantillonnage spatial, à partir d'une pluralité de signaux reçus chacun par l'un desdits récepteurs, Définir, dans l'image, au moins une zone d'intérêt de l'objet à 5 imager, Déterminer une seconde image de l'objet à partir d'une pluralité de signaux reçus chacun par l'un desdits récepteurs, la seconde image présentant, dans chaque zone d'intérêt, un second niveau d'échantillonnage spatial supérieur au premier niveau 10 d'échantillonnage spatial. Selon un aspect particulier de l'invention, le premier niveau d'échantillonnage spatial est défini en décomposant l'objet en cellules et en sélectionnant un point dans chaque cellule. Selon un aspect particulier de l'invention, la sélection d'un point dans 15 chaque cellule suit une distribution aléatoire ou une distribution prédéfinie. Selon une variante particulière, le procédé d'imagerie selon l'invention comprend en outre une étape de sur-échantillonner la première image en affectant à chaque pixel d'un point d'une cellule la valeur du pixel calculé pour le point sélectionné de ladite cellule. 20 Selon une variante particulière, le procédé d'imagerie selon l'invention comprend en outre une étape de sur-échantillonner la première image en effectuant une interpolation des pixels de la première image. Selon un aspect particulier de l'invention, l'interpolation est polynômiale, linéaire ou réalisée par procédé Gaussien. 25 Selon un aspect particulier de l'invention, au moins une zone d'intérêt est élargie en intégrant tous les points situés à une distance prédéterminée de la zone d'intérêt. Selon un aspect particulier de l'invention, la distance prédéterminée dépend de la longueur d'onde de l'onde élastique. 30 Selon un aspect particulier de l'invention, une zone d'intérêt comprend au moins un pixel dont la valeur est un extrema local. 5 Selon un aspect particulier de l'invention, la définition, dans la première image, d'au moins une zone d'intérêt de l'objet à imager comprend les sous-étapes de : 5 - définir un seuil d'amplitude, - comparer chaque pixel de la première image au seuil d'amplitude, - définir si le pixel de la première image appartient ou non à une zone d'intérêt en fonction du résultat de la comparaison. Selon un aspect particulier de l'invention, l'image de l'objet est une 10 image en deux dimensions ou en trois dimensions. Selon une variante particulière, le procédé d'imagerie selon l'invention comprend en outre les étapes préalables de : Emettre, depuis chaque émetteur, une onde élastique dans l'objet à imager, 15 Recevoir, dans chaque récepteur, une onde élastique provenant de l'objet à imager. Selon un aspect particulier de l'invention, la détermination d'une image de l'objet comprend le calcul d'une pluralité de pixels correspondant chacun à un point de l'objet en exécutant les sous-étapes de : 20 Pour chaque couple émetteur-récepteur, déterminer un temps de vol correspondant à une durée théorique nécessaire à l'onde élastique pour parcourir un trajet depuis l'émetteur jusqu'au récepteur en passant par le point de l'objet, Déterminer une somme des amplitudes extraites d'un ensemble 25 d'ondes élastiques émises par lesdits émetteurs et reçus par lesdits récepteurs, aux temps de vol déterminés à l'étape précédente. Selon un aspect particulier de l'invention, l'onde élastique est une onde ultrasonore ou une onde guidée. 30 L'invention a aussi pour objet un dispositif d'imagerie pour imager un objet, le dispositif comprenant une interface pour recevoir une pluralité de 6 signaux reçus par une pluralité de récepteurs aptes à recevoir une onde élastique provenant d'un objet à imager et un imageur configuré pour mettre en oeuvre les étapes du procédé d'imagerie selon l'un quelconque des modes de réalisation de l'invention. 5 Dans une variante de réalisation, le dispositif d'imagerie selon l'invention comprend une pluralité d'émetteurs aptes à émettre une onde élastique dans l'objet à imager et une pluralité de récepteurs aptes à recevoir une onde élastique provenant de l'objet à imager. 10 D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : La figure 1, un schéma d'un dispositif d'imagerie selon un mode de réalisation de l'invention, 15 La figure 2, un organigramme détaillant les étapes de mise en oeuvre d'un procédé d'imagerie selon différents modes de réalisation de l'invention, La figure 3, un exemple illustrant l'imagerie d'une pièce, Les figures 3a,3a-bis,3b,3c,3d, cinq exemples d'images obtenues 20 après différentes étapes intermédiaires du procédé selon un mode de réalisation de l'invention appliqué à la pièce décrite à la figure 3, Les figures 4a,4b,4c, trois autres exemples d'images obtenues après différentes étapes intermédiaires du procédé selon un autre mode de réalisation de l'invention. 25 La figure 1 représente, sur un schéma, un exemple de dispositif d'imagerie par ondes élastiques selon un mode de réalisation de l'invention. Un tel dispositif comporte au moins un transducteur multiéléments TR comprenant une pluralité d'émetteurs/récepteurs d'ondes élastiques ChCi. 30 Chaque élément 0,2 du transducteur TR est apte à émettre une onde élastique dans un objet OBJ et à recevoir une onde élastique provenant d'un 7 objet. L'onde élastique est, par exemple, une onde ultrasonore ou une onde guidée. Lorsque l'onde élastique est une onde ultrasonore, le transducteur TR est, par exemple, un capteur piézoélectrique. Les émetteurs/récepteurs Ci,Ci sont agencés linéairement, sous forme de matrice ou distribués à 5 différentes positions sur ou à proximité de l'objet. Dans le cas d'un agencement compact linéaire, le dispositif est apte à produire une image en deux dimensions d'un plan de l'objet. Dans le cas d'un agencement compact matriciel, le dispositif est apte à produire une image en trois dimensions d'un volume d'une partie de l'objet. 10 Pour réaliser une acquisition de type TFM, c'est-à-dire une procédure d'émission et de réception successives d'ondes élastiques à travers l'objet à imager, un fluide couplant est, parfois, disposé entre le dispositif d'imagerie et l'objet à imager. Le fluide couplant est, par exemple, de l'eau. Le dispositif selon l'invention comporte en outre un imageur IMG 15 configuré pour déterminer une image de l'objet OBJ à partir des signaux reçus par les récepteurs du transducteur TR. L'imageur IMG est, par exemple, un processeur graphique, un processeur de signaux ou encore un circuit intégré, un circuit logique programmable, un circuit intégré propre à une application ou tout autre dispositif équivalent configurable pour exécuter 20 les étapes du procédé selon l'un quelconque des modes de réalisation de l'invention. Le dispositif selon l'invention peut aussi comporter un afficheur AFF, tel qu'un écran ou toute autre interface homme machine, pour restituer à un utilisateur une image de l'objet OBJ déterminée par l'imageur IMG. L'image 25 peut être en deux dimensions ou en trois dimensions. Dans un mode de réalisation particulier de l'invention, le dispositif comprend uniquement l'imageur IMG qui reçoit les signaux mesurés par les récepteurs du transducteur multi-éléments TR au moyen d'une interface adaptée. Dans ce scénario, le transducteur TR est distinct de l'imageur IMG 30 qui peut être, ou non, associé à un afficheur AFF. 8 Dans le cas où les ondes élastiques sont des ondes guidées, les différents capteurs du dispositif peuvent être positionnés en différents endroits de l'objet à imager. Un avantage à l'utilisation d'ondes guidées est qu'elles ont la capacité à être guidées dans une structure dans laquelle elles 5 se propagent sur une plus longue distance que les ondes ultrasonores. Ainsi, il est possible d'imager des objets de plus grande taille en utilisant des ondes guidées. La figure 2 détaille les étapes de réalisation d'un procédé d'imagerie 10 selon un mode de réalisation de l'invention, le procédé étant mis en oeuvre au moyen d'un dispositif du type décrit à la figure 1. Une première étape 201 du procédé consiste à déterminer une acquisition de signaux à l'aide d'un transducteur multiéléments du type décrit à la figure 1. Cette étape 201 peut être réalisée de différentes façons. Le 15 principe général d'une telle acquisition consiste à émettre successivement, au moyen de chaque émetteur du transducteur multiéléments, une onde élastique propre à se propager dans l'objet à imager, puis à acquérir successivement, au moyen de chaque récepteur du transducteur multiéléments, chacun des signaux émis par chaque émetteur. 20 Sur la figure 1 on a représenté le trajet d'une onde élastique émise par un émetteur Ci se propageant dans un objet OBJ, l'onde est réfléchie sur un point P de l'objet puis est reçue par un récepteur C. L'onde peut être une onde ultrasonore ou une onde guidée. Si le transducteur multiéléments comprend N émetteurs/récepteurs, le 25 nombre total de signaux acquis est égale à N2. La procédure d'acquisition, objet de l'étape 201 du procédé, suit, par exemple, une procédure du type cc Full Matrix Capture » ou FMC. Une deuxième étape 202 du procédé consiste à réaliser un pré-30 échantillonnage spatial de l'image, à faible densité. Autrement dit, lors de cette étape 202, on détermine, pour un premier ensemble de points 01 de 9 l'objet à imager, un pixel correspondant d'une image de l'objet. Le premier ensemble de points C), est choisi de sorte à ce que le nombre de points de cet ensemble est significativement inférieur au nombre de points nécessaire pour obtenir une image de l'objet à une résolution prédéterminée. Pour cela, 5 la résolution souhaitée de l'image ainsi que sa dimension sont, par exemple, fournis comme paramètres de l'étape 202, en fonction de l'application visée ou du matériau constitutif de l'objet. En général, le nombre de pixels souhaité pour obtenir une image permettant une inspection précise d'un objet dépend de la longueur d'onde de l'onde élastique utilisée. Par exemple, la distance 10 entre deux pixels est prise égale à un ratio de la longueur d'onde À, par exemple À/10 ou À/20. Ainsi, le premier ensemble de points 01 est, par exemple, choisi pour contenir un faible pourcentage, de l'ordre de 1%, du nombre de points normalement requis pour obtenir une image à la résolution souhaitée. 15 Les positions des points du premier ensemble de points 01 peuvent être choisies de différentes façons. Un premier choix consiste à sélectionner les points de façon aléatoire ou pseudo-aléatoire parmi l'ensemble des points de l'objet. Un autre choix consiste à sélectionner les points selon une distribution prédéterminée régulière, par exemple une distribution de 20 Poisson. La distance entre deux points sélectionnés est, par exemple, dépendante de la longueur d'onde de l'onde élastique émise à travers l'objet. Un mode de réalisation particulier de l'étape 202 consiste à segmenter l'objet à imager en cellules, de mêmes tailles ou de tailles différentes, puis à sélectionner, de façon aléatoire ou déterministe, un point dans chaque 25 cellule. Le nombre et les coordonnées des points du premier ensemble de points 01 peuvent aussi être prédéterminés a priori pour toute image ayant une taille donnée, ou pour une application donnée. 30 Ensuite, pour chaque point du premier ensemble de points 01, on détermine un pixel correspondant de l'image de l'objet à partir de l'ensemble 10 des signaux mesurés par les récepteurs à l'étape 201. Autrement dit, l'étape 201 consiste à imager l'objet avec un premier niveau d'échantillonnage spatial. Pour déterminer un pixel, différentes méthodes sont possibles. 5 Une méthode possible est la méthode dite TFM (Total Focusing Method). Cette méthode consiste à calculer, pour chaque point P de la zone de l'objet à imager, le pixel correspondant 1(P) comme étant égal à la somme des amplitudes extraites des signaux S1(t) reçus par les récepteurs du transducteur multiéléments, aux temps de vol théoriques t= T1(P)= t1(P) + ti(P) 10 correspondant aux chemins parcourus entre un émetteur 0 et un récepteur 0 en passant par le point P. Autrement dit, si N est le nombre de récepteurs et M le nombre d'émetteurs, on somme les amplitudes des M x N signaux reçus par les récepteurs aux instants correspondant aux durées respectives nécessaires 15 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 Al 1(P) = Su (Ti j (P)) i=1 j=1 Selon une variante de cette méthode, les termes de la somme peuvent être pondérés par des coefficients de pondération W;,1 afin de 20 calibrer les différents signaux pour prendre en compte des différences de gain entre les différents capteurs. Un exemple de réalisation de la méthode d'acquisition TFM est notamment décrit dans les documents [1] et [4]. D'autres méthodes sont possibles pour calculer la valeur d'un pixel 25 1(P) à partir des signaux reçus S1(t). On peut citer, sans être exhaustif, les méthodes suivantes : SAFT (Synthetic Aperture Focusing Technique) décrite dans [5], MUSIC (Multiple Signal Classification), Excitelet, DAS (Delay and Sum) décrite dans [6], ou encore les méthodes décrites dans les références [7] à [10] qui portent plus précisément sur l'utilisation d'ondes guidées. 11 De manière générale, l'invention est applicable à toute méthode d'imagerie permettant de déterminer un pixel pour un point d'un objet, en deux dimensions ou en trois dimensions, à partir des signaux reçus Sii(t). 5 La figure 3 illustre schématiquement la mise en oeuvre de l'invention pour imager un objet 300 au moyen d'un dispositif d'imagerie 301 apte à émettre des ondes élastiques 302 pour image l'objet 300. Sur la figure 3, on a représenté en superposition à l'objet 300, une image 303 obtenue via une méthode d'imagerie TFM selon l'art antérieur. Sur l'image, on peut visualiser 10 un défaut 310 ainsi que des échos 311 correspondant au fond de l'objet 300. Il convient de noter que les images suivantes illustrées sur les figures 3a, 3a-bis, 3b, 3c, 3d, 4a, 4b, 4c sont inversées selon l'axe vertical par rapport à l'image de la figure 3. 15 La figure 3a illustre le résultat obtenu après l'étape 202, pour l'exemple particulier décrit à la figure 3. Selon cet exemple, l'image est décomposée en cellules et un pixel (représenté par un point sur la figure 3a) est calculé dans chaque cellule. Sur la figure 3a-bis, on a représenté, pour illustration, en fond et en 20 superposition à l'image de la figure 3a, l'image de l'objet obtenue avec une méthode d'imagerie TFM selon l'art antérieur ainsi que les cellules. L'étape suivante 203 est une étape d'affinement de l'image calculée à l'étape 202. Il convient de noter que cette étape est optionnelle et que, dans 25 un mode de réalisation particulier de l'invention, elle peut être omise, auquel cas l'étape suivante 204 est directement appliquée au résultat de l'étape 202. Selon cette étape d'affinement 203, on effectue un sur-échantillonnage de l'image à partir des pixels calculés à l'étape 202. Ce sur-échantillonnage peut être réalisé de multiples façons. 30 Une première réalisation consiste à reprendre la décomposition en cellules effectuée à l'étape 202 et à affecter à tous les pixels d'une cellule, la 12 valeur du pixel de la cellule calculé à l'étape 202. La distance entre deux pixels calculés est, par exemple, prédéterminée et peut être fonction de la longueur d'onde de l'onde élastique. Une autre réalisation possible de l'étape 203 consiste à réaliser une 5 interpolation des pixels calculés à l'étape 202 pour affiner l'image. Différents types d'interpolation sont possibles, il peut s'agir d'une interpolation Gaussienne ou d'une interpolation linéaire ou encore d'une interpolation quadratique ou polynômiale ou bilinéaire polynômiale. L'interpolation bilinéaire peut s'interpréter comme une succession de 10 deux interpolations linéaires, une dans chaque direction. Contrairement à ce que son nom suggère, la fonction d'interpolation bilinéaire polynômiale n'est pas une forme linéaire mais quadratique, qui peut se mettre sous la forme : f(x,y)-ax+by-Ecxy+d f(x ,y) est la valeur interpolée au point de coordonnées ( x , y ) et a , b, 15 c et d sont des constantes déterminées à partir des quatre points voisins ( x1 ,y1 ), ( x2 , y1 ) , ( x1 , y2 ) , ( x2 , y2 ) du point ( x , y ) dont on cherche la valeur. Dans le cas où l'imagerie est en trois dimensions, on applique une zo interpolation en trois dimensions, par exemple une succession de trois interpolations linéaires, une dans chaque direction. Un objectif de l'interpolation est de permettre une identification, dans l'image, des zones d'amplitude les plus fortes correspondants aux défauts que l'on souhaite identifier. 25 L'interpolation Gaussienne ou régression Gaussienne ou interpolation par procédé Gaussien (en anglais Gaussian Process Regressor ou GPR) est une méthode dite de « machine learning » supervisée souvent présentée comme une alternative à l'interpolation classique. Il s'agit d'un procédé 30 itératif qui observe les points fournis par l'étape 202 du procédé un à un, en adaptant sa régression au fur et à mesure. La régression Gaussienne n'est 13 pas une régression paramétrique (c'est-à-dire qu'elle ne définit pas un nombre fixe de paramètres à estimer, par exemple la régression linéaire 'y= ax+b' a deux paramètres), mais elle essaie de réaliser la régression de manière probabiliste en cherchant des distributions gaussiennes permettant 5 de prédire au mieux la fonction à reconstruire. Un exemple de régression Gaussienne est décrit dans la référence [11]. En particulier, un avantage à l'utilisation d'une interpolation Gaussienne est qu'elle permet de mieux représenter les extrema de l'image 10 qui correspondent à des réflecteurs. Un autre moyen d'affiner la représentation des extrema de l'image consiste à augmenter le nombre de points sur lequel est réalisé l'interpolation. 15 Sur la figure 3b, on a représenté un exemple de résultat de l'étape 203 dans le cas d'une interpolation Gaussienne. On peut identifier sur cette image un indicateur de défaut aux alentours de la position (x,y)=(150,100) et une zone en haut de l'image correspondant au fond de la pièce. 20 Dans une étape suivante 204, on détermine, à partir de l'image obtenue après l'étape 203, ou après l'étape 202 si l'étape 203 est omise, un ensemble de zones d'intérêts de l'image qui correspondent aux zones de l'objet dans lesquelles se situe potentiellement un défaut ou plus généralement un élément que l'on souhaite visualiser avec précision. 25 Autrement dit, les zones d'intérêts de l'image correspondent aux zones de l'objet que l'on souhaite imager avec une résolution importante. De façon générale, pour définir une zone d'intérêt, on cherche à identifier les extrema locaux des valeurs des pixels de l'image qui correspondent, en moyenne, à des points remarquables de l'image que l'on 30 souhaite visualiser avec précision. 14 Une réalisation possible de l'étape 204 consiste à définir un seuil auquel sont comparées les valeurs de chaque pixel de l'image. En fonction du résultat de la comparaison, le pixel est classé ou non dans une zone d'intérêt. Si les valeurs des pixels présentent des valeurs positives et 5 négatives, la comparaison est réalisée entre la valeur absolue d'un pixel et le seuil. La valeur du seuil est un paramètre de l'invention qui peut être déterminée en fonction de la moyenne des valeurs des pixels et des valeurs absolues maximales des pixels. 10 Le seuil peut être, par exemple déterminé comme un pourcentage de la valeur moyenne de l'ensemble des pixels calculés à l'étape précédente. Ce pourcentage est, par exemple, égal à 100%,50% ou 200% ou toute autre valeur paramétrée par l'utilisateur. Il s'agit d'un seuil en proportion de 15 l'énergie totale. Les pixels de l'image qui présentent une valeur absolue supérieure à ce seuil sont sélectionnés pour faire partie d'une zone d'intérêt. D'autres possibilités existent pour fixer la valeur du seuil. Elle peut être fixe mais aussi réglable en fonction de l'application visée ou déterminé à partir des premières mesures sur l'objet à imager. 20 Une autre alternative est de définir une zone d'intérêt en se basant sur un nombre de points (et non pas sur l'énergie). Par exemple, une zone d'intérêt peut être définie en sélectionnant un nombre prédéterminé de pixels dont la valeur est la plus grande dans l'image. Une autre alternative est de calculer les maxima et/ou minima locaux 25 de la cartographie estimée et de sélectionner un nombre prédéterminé de pixels aux alentours de ces extrema pour définir les zones d'intérêts. La figure 3c donne un exemple de zones d'intérêts identifiées à partir de l'image obtenue sur la figure 3b. Les zones d'intérêts sont détectées 30 d'une part autour d'un défaut situé approximativement aux coordonnées 15 (x=150, y=100) et d'autre part en haut de l'image, qui correspond au fond de la pièce imagée. Dans une dernière étape 205, on défini un second ensemble 02 de 5 points de l'objet, situés exclusivement dans les zones d'intérêts définies à l'étape 204 précédente, puis on calcule les pixels associés à ces points. Un objectif de l'étape 205 est de réaliser un échantillonnage spatial de l'image uniquement dans les zones d'intérêts définies, c'est-à-dire pour les zones de l'objet que l'on souhaite imager avec une résolution importante. Le 10 niveau d'échantillonnage choisi est supérieur au niveau d'échantillonnage utilisé pour l'étape de pré-échantillonnage 202. Ainsi, on défini dans les zones d'intérêts une grille d'échantillonnage plus dense que celle choisie à l'étape 202 pour définir le premier ensemble Qi de points. Dans une variante particulière de l'étape 205, l'échantillonnage spatial 15 est appliqué dans les zones d'intérêts élargies aux points situés à une distance prédéfinie de ces zones. La distance dépend, par exemple, de la longueur d'onde de l'onde élastique. Cette variante permet de s'assurer que des points importants, situés entre deux zones d'intérêts par exemple, sont bien sélectionnés. 20 Le nombre de points ou la densité des points du second ensemble 02 est choisie en fonction de la résolution souhaitée. Ces paramètres peuvent dépendre de la longueur d'onde de l'onde élastique. Le second ensemble Q2 peut comprendre des pixels déjà calculés à l'étape 202, c'est-à-dire pour des points qui appartiennent aussi au premier 25 ensemble Qi. Selon l'implémentation choisie, ces pixels sont recalculés ou non. La figure 3d illustre un exemple de sur-échantillonnage réalisé à l'étape 205 pour les zones d'intérêts définies à l'étape 204. La figure 3d correspond à l'image finale qui est composée des pixels calculés à l'étape 30 202 de pré-échantillonnage et des pixels calculés à l'étape 205 de sur-échantillonnage. 16 Pour une image de même taille et même résolution, le gain de réduction du nombre de pixels calculés pour l'image de la figure 3d par rapport à une image calculée avec une méthode TFM classique est de l'ordre de 80%. L'erreur moyenne estimée entre l'image de la figure 3d et une image 5 calculée avec une méthode TFM classique est de l'ordre de 4%. Les figures 4a,4b,4c illustrent les résultats obtenus par le procédé selon l'invention lorsque l'étape 202 ne comporte pas d'interpolation mais un simple sur-échantillonnage selon lequel chaque pixel d'une cellule prend la 10 valeur du pixel de la cellule calculé à l'étape 201. La figure 4a illustre l'image obtenue après l'étape 201. La figure 4b illustre les zones d'intérêts déterminées dans l'image après l'étape 203. La figure 4c illustre l'image finale obtenue après l'étape 204. Le gain de réduction du nombre de pixels calculés pour l'image de la 15 figure 4c par rapport à une image calculée avec une méthode TFM est de l'ordre de 89%. L'erreur moyenne entre les deux images est de l'ordre de 30/0. Dans une variante de réalisation de l'invention, l'étape 202 ou la 20 combinaison des étapes 202 et 203 est/sont itérée(s) dans chaque zone d'intérêt définie à l'étape 204 afin de permettre de mettre à jour la définition des zones d'intérêts à chaque itération et ainsi d'améliorer la précision de la définition de ces zones. 17 Références [1] J. Lambert, A. Pedron, G. Gens, F. Bimbard, L. Lacassagne et E. lakovleva, «Performance evaluation of total focusing method on GPP and 5 GPU,» Design and Architectures for Signal and Image Processing (DASIP), 2012 Conference on, 2012. [2] G. Rougeron, J. Lambert, E. lakovleva, L. Lacassagne et N. Dominguez, «Implementation of a GPU accelerated total focusing reconstruction method within CIVA software,» AIP Conference Proceedings, 2014. 10 [3] J. Guo, X. Li, X. Gao, Z. Wang et Q. Zhao, «Implementation of total focusing method for phased array ultrasonic imaging on FPGA,» Ninth International Symposium on Precision Engineering Measurement and Instrumentation, 2015. [4] M. Njiki, A. Elouardi, S. Bouaziz, O. Casula et O. Roy, «A real-time 15 implementation of the total focusing method for rapid and precise diagnostic in non destructive evaluation,» Application-Specific Systems, Architectures and Processors (ASAP), 2013 IEEE 24th International Conference on. [5] SAFT : Mayer, K., et al. "Three-dimensional imaging system based on Fourier transform synthetic aperture focusing technique." Ultrasonics 28.4 20 (1990): 241-255. [6] MUSIC : Schmidt, Ralph. "Multiple emitter location and signal parameter estimation." IEEE transactions on antennas and propagation 34.3 (1986): 276-280. [7] Excitelet : Quaegebeur, N.; Masson, P.; Langlois-Demers, D. & Micheau, 25 P. Dispersion-based imaging for structural health monitoring using sparse and compact arrays Smart Materials and Structures, 10P Publishing, 2011, 20, 025005 [8] DAS : Michaels, J. E. Detection, localization and characterization of damage in plates with an in situ array of spatially distributed ultrasonic 30 sensors Smart Materials and Structures, 10P Publishing, 2008, 17, 035035 18 [9]MVDR : Hall, J. S. & Michaels, J. E. Minimum variance ultrasonic imaging applied to an in situ sparse guided wave array IEEE transactions on ultrasonics, ferroelectrics, and frequency control, IEEE, 2010, 57, 2311-2323 [10] RAPID : Zhao, X.; Gao, H.; Zhang, G.; Ayhan, B.; Yan, F.; Kwan, C. & 5 Rose, J. L. Active health monitoring of an aircraft wing with embedded piezoelectric sensor/actuator network: I. Defect detection, localization and growth monitoring Smart materials and structures, 10P Publishing, 2007, 16, 1208, [11] M.Ebden, "Gaussian Processes for Regression: a quick introduction", REVENDICATIONS1. Procédé d'imagerie pour imager un objet (OBJ) à l'aide d'un dispositif d'imagerie comprenant une pluralité d'émetteurs (C1) aptes à émettre une onde élastique dans l'objet à imager et une pluralité de récepteurs (C1) aptes à recevoir une onde élastique provenant de l'objet à imager, ledit procédé comprenant les étapes de : Déterminer (202) une première image de l'objet, à un premier niveau d'échantillonnage spatial, à partir d'une pluralité de signaux reçus chacun par l'un desdits récepteurs, Définir (204), dans l'image, au moins une zone d'intérêt de l'objet à imager, Déterminer (205) une seconde image de l'objet à partir d'une pluralité de signaux reçus chacun par l'un desdits récepteurs, la seconde image présentant, dans chaque zone d'intérêt, un second niveau d'échantillonnage spatial supérieur au premier niveau d'échantillonnage spatial. 2. Procédé d'imagerie selon la revendication 1 dans lequel le premier niveau d'échantillonnage spatial est défini (204) en décomposant l'objet en cellules et en sélectionnant un point dans chaque cellule. 3. Procédé d'imagerie selon la revendication 2 dans lequel la sélection d'un point dans chaque cellule suit une distribution aléatoire ou une distribution prédéfinie. 4. Procédé d'imagerie selon l'une des revendications 2 ou 3 comprenant en outre une étape (203) de sur-échantillonner la première image en affectant à chaque pixel d'un point d'une cellule la valeur du pixel calculé pour le point sélectionné de ladite cellule. 20 5. Procédé d'imagerie selon l'une des revendications 1,2 ou 3 comprenant en outre une étape (203) de sur-échantillonner la première image en effectuant une interpolation des pixels de la première image. 5 6. Procédé d'imagerie selon la revendication 5 dans lequel l'interpolation est polynômiale, linéaire ou réalisée par procédé Gaussien. 7. Procédé d'imagerie selon l'une des revendications précédentes dans lequel au moins une zone d'intérêt est élargie en intégrant tous les points 10 situés à une distance prédéterminée de la zone d'intérêt. 8. Procédé d'imagerie selon la revendication 7 dans lequel la distance prédéterminée dépend de la longueur d'onde de l'onde élastique. 15 9. Procédé d'imagerie selon l'une des revendications précédentes dans lequel une zone d'intérêt comprend au moins un pixel dont la valeur est un extrema local. 10. Procédé d'imagerie selon l'une des revendications précédentes dans 20 lequel la définition (204), dans la première image, d'au moins une zone d'intérêt de l'objet à imager comprend les sous-étapes de : définir un seuil d'amplitude, comparer chaque pixel de la première image au seuil d'amplitude, définir si le pixel de la première image appartient ou non à une 25 zone d'intérêt en fonction du résultat de la comparaison. 11. Procédé d'imagerie selon l'une des revendications précédentes dans lequel l'image de l'objet est une image en deux dimensions ou en trois dimensions. 30 12. Procédé d'imagerie selon l'une des revendications précédentes comprenant en outre les étapes préalables (201) de : - Emettre, depuis chaque émetteur, une onde élastique dans l'objet à imager, 21 - Recevoir, dans chaque récepteur, une onde élastique provenant de l'objet à imager. 13. Procédé d'imagerie selon l'une des revendications précédentes dans 5 lequel la détermination (202,205) d'une image de l'objet comprend le calcul d'une pluralité de pixels correspondant chacun à un point (P) de l'objet en exécutant les sous-étapes de : Pour chaque couple émetteur-récepteur (Ci,0), déterminer un temps de vol correspondant à une durée théorique nécessaire à 10 l'onde élastique pour parcourir un trajet depuis l'émetteur (Ci) jusqu'au récepteur (Ci) en passant par le point (P) de l'objet, Déterminer une somme des amplitudes extraites d'un ensemble d'ondes élastiques émises par lesdits émetteurs et reçus par lesdits récepteurs, aux temps de vol déterminés à l'étape 15 précédente. 14. Procédé d'imagerie selon l'une des revendications précédentes dans lequel l'onde élastique est une onde ultrasonore ou une onde guidée. 20 15. Dispositif d'imagerie pour imager un objet, le dispositif comprenant une interface pour recevoir une pluralité de signaux reçus par une pluralité de récepteurs aptes à recevoir une onde élastique provenant d'un objet (OBJ) à imager et un imageur (IMG) configuré pour mettre en oeuvre les étapes du procédé d'imagerie selon l'une quelconque des revendications 25 précédentes. 16. Dispositif d'imagerie selon la revendication 15, comprenant une pluralité d'émetteurs (Ci) aptes à émettre une onde élastique dans l'objet à imager et une pluralité de récepteurs (Ci) aptes à recevoir une onde élastique 30 provenant de l'objet à imager.