FR 3008806 A1 20150123 FR 1357178 A 20130722 Procédé et dispositif d'imagerie acousto-électrique. DOMAINE DE L'INVENTION La présente invention est relative aux procédés et 5 dispositifs d'imagerie acousto-électrique. Des organes comme le coeur, les muscles squelettiques et le cerveau sont parcourus en permanence par des impulsions électriques qui transportent les informations dans les neurones, ou qui déclenchent les contractions 10 musculaires ou du myocarde. Pouvoir imager la propagation de ces impulsions est extrêmement important pour diagnostiquer de nombreuses pathologies et comprendre les mécanismes cérébraux par l'exploration fonctionnelle du cerveau. 15 L'imagerie acousto-électrique exploite l'interaction entre les ultrasons et les courants électriques pour déterminer la valeur du courant électrique aux points d'interaction entre les ultrasons et les tissus, classiquement au niveau de la tache focale d'une onde 20 ultrasonore focalisée. ARRIERE-PLAN DE L'INVENTION Le document US8057390 décrit un exemple de procédé d'imagerie acousto-électrique, dans lequel on émet des ondes ultrasonores focalisées pour former, ligne par ligne, 25 une image du courant. Ce processus d'acquisition est lent, et ce d'autant plus que, les signaux électriques obtenus étant très faibles, un niveau de moyennage important est nécessaire. On obtient donc de faibles débits d'images. RESUME DE L'INVENTION 30 La présente invention a notamment pour but de pallier cet inconvénient. A cet effet, l'invention propose un procédé d'imagerie acousto-électrique, comprenant : (a) une étape de mesure au cours de laquelle on 35 fait émettre par un ensemble de transducteurs Ti, dans un champ d'observation d'un milieu à imager, un nombre N au moins égal à 2 d'ondes ultrasonores incidentes 1 non focalisées dans le champ d'observation et ayant des fronts d'ondes différents, et on fait capter par au moins un capteur électrique en contact avec le milieu à imager, des signaux électriques bruts Erawl(t) respectivement pendant la propagation des ondes incidentes 1, (b) une étape de formation d'image, au cours de laquelle on détermine à partir des signaux électriques bruts Erawl(t) obtenus à l'étape (a), une image du milieu comprenant une cartographie de courants électriques (c'est-à-dire une cartographie de valeurs électriques représentatives des densités de courant locales en chaque point du milieu). Grâce à ces dispositions, on peut réaliser une imagerie ultrarapide des impulsions électriques dans le milieu observé, et éventuellement filmer la propagation d'impulsions électrique en profondeur dans les tissus, en temps-réel et avec une résolution millimétrique. Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - au cours de l'étape b) on détermine, au moins à partir des N signaux électriques bruts Erawl(t), pour un 25 nombre M de points de focalisation fictifs Pk dans champ d'observation, des valeurs électriques Ecoherentk correspondant chacune au signal électrique qui aurait été capté si une onde ultrasonore focalisée au point Pk avait été émise par lesdits transducteurs ; 30 - au cours de l'étape (b), on applique aux signaux électriques bruts Erawl(t) une transformée en ondelettes inverse WT 1, puis une transformée de Radon inverse R (les signaux électriques bruts Erawl(t) peuvent bien entendu subir un traitement préliminaire avant la transformée de 35 Radon inverse R1); - au cours de l'étape (b), on superpose à la cartographie de courants électriques, une image échographique du milieu réalisée avec l'ensemble de transducteurs ; - au cours de l'étape (a), on fait capter par les transducteurs Ti des signaux acoustiques RFrawl,i(t) représentatifs d'ondes ultrasonores réverbérées par le milieu respectivement à partir des ondes incidentes 1, au cours de l'étape (b), on détermine à partir des N 10 ensembles de signaux RFrawl,i(t) captés, M signaux acoustiques cohérents RFcoherentk,i(t) correspondant aux signaux acoustiques qui auraient été reçus par les transducteurs Ti si une onde ultrasonore focalisée au point Pk avait été émise par lesdits transducteurs, et on calcule 15 l'image échographique du milieu à partir des signaux acoustiques cohérents ; au cours de l'étape (b), on détermine l'image échographique par formation de voies à partir des signaux acoustiques cohérents ; 20 le milieu à imager est un tissu humain ou animal. Par ailleurs, l'invention a également pour objet un dispositif pour la mise en oeuvre d'un procédé d'imagerie acousto-électrique, comprenant un ensemble de transducteurs Ti, moins un capteur électrique, et des moyens de commande 25 et traitement adaptés pour : (a) faire émettre par un ensemble de transducteurs Ti, dans un milieu à imager, un nombre N d'ondes ultrasonores incidentes non focalisées 1 ayant des fronts d'ondes différents, et faire capter par au moins un capteur 30 électrique en contact avec le milieu à imager, des signaux électriques bruts Erawl(t) respectivement pendant la propagation des ondes incidentes 1, (b) déterminer à partir des signaux électriques bruts Erawl(t), une image du milieu comprenant une 35 cartographie de courants électriques. BREVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre 5 d'exemple non limitatif, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue schématique d'un dispositif pour la mise en oeuvre d'un procédé selon un mode de réalisation de l'invention, et 10 - la figure 2 est un schéma bloc d'une partie du dispositif de la figure 1. DESCRIPTION DETAILLEE Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. 15 La figure 1 montre un exemple de dispositif d'imagerie acousto-électrique adapté pour imager un milieu 1 par émission et réception d'ondes ultrasonores de compression (par exemple de fréquences comprises entre 0.2 et 40 MHz), avec mesure simultanée de valeurs électriques. 20 Le milieu 1 à imager peut consister notamment en des tissus d'un patient ou d'un animal, notamment un muscle (myocarde ou autre) ou un cerveau. Le dispositif d'imagerie comporte par exemple : - un réseau 2 de n transducteurs ultrasonores, 25 comprenant par exemple quelques centaines de transducteurs et adapté pour réaliser une image bidimensionnelle (2D) d'un champ d'observation (zone d'intérêt, balayée par les ondes ultrasonores) dans le milieu 1 à imager ; - une baie électronique 3 ou similaire commandant 30 le réseau 2 de transducteurs et adaptée pour acquérir les signaux captés par ce réseau de transducteurs ; - un ordinateur 4 ou similaire pour commander la baie électronique 3 et visualiser les images ultrasonores obtenues à partir desdits signaux captés. 35 Le réseau 2 de transducteurs peut par exemple être un réseau linéaire formé par une barrette de transducteurs juxtaposés selon un axe X, l'axe Z perpendiculaire à l'axe X désignant la direction de la profondeur dans le champ d'observation. Dans ce qui suit, les transducteurs seront notés Ti, i étant un indice désignant le rang de chaque transducteur selon l'axe X. La description qui suit sera faite en prenant pour exemple ce type de réseau 2 de transducteurs, mais d'autres formes de de réseau de transducteurs sont également possibles dans le cadre de la présente invention, notamment des réseaux bidimensionnels. Le dispositif comporte en outre au moins un capteur électrique El (figure 2), constitué par exemple par deux électrodes mesurant une différence de potentiel électrique. Ce capteur électrique peut avantageusement être fixé au le réseau 2 de transducteurs et adapté pour entrer en contact avec le milieu 1 à imager en même temps que les transducteurs du réseau 2. Comme représenté sur la figure 2, la baie électronique 3 peut comprendre par exemple : - n+1 convertisseurs analogique / digital 5 (A/Di - A/De) connectés individuellement aux n transducteurs Ti du réseau 2 de transducteurs et au capteur électrique El, - n+1 mémoires tampon 6 (Bi-Be) respectivement connectées aux n convertisseurs analogique / digital 5, - une unité centrale 8 (CPU) communiquant avec les mémoires tampon 6 et l'ordinateur 4, - une mémoire 9 (MEM) connectée à l'unité centrale 8, - un processeur numérique de signal 10 (DSP) 30 connecté à l'unité centrale 8. On notera que les n+1 convertisseurs analogique / digital 5 (A/Di - A/De) peuvent être identiques, de même que les n+1 mémoires tampon 6 (Bi-Be), de sorte que le dispositif utilisé peut être simplement un dispositif 35 utilisé classiquement en imagerie acoustique ultrarapide. Ce dispositif permet de mettre en oeuvre un procédé d'imagerie acousto-électrique du milieu 1, qui inclut notamment les étapes suivantes, mises en oeuvre par l'unité centrale 8 assistée du processeur 8 et du processeur numérique de signal 10 : a) Mesure (émission / réception et enregistrement des données brutes), b) détermination d'une image du milieu comprenant une cartographie de valeurs électriques. Etape (a) : Mesure (émission / réception et enregistrement des données brutes) : Le réseau 2 de transducteurs et le capteur électrique El sont mis en contact avec le milieu 1 et un nombre N d'ondes incidentes ultrasonores est émis dans le milieu 1 par les transducteurs Ti (N peut être compris par exemple entre 2 et 100, notamment entre 5 et 10). Les ondes incidentes en question sont non focalisées (plus précisément, non focalisées dans le champ d'observation) et ont des respectivement fronts d'onde différents, c'est-à- dire des fronts d'onde de formes différentes et / ou d'orientation différentes. Avantageusement, les ondes incidentes peuvent être des ondes planes dont les fronts d'onde F respectifs (le front d'onde F d'une seule onde est représenté sur la figure 1) ont des inclinaisons toutes différentes, caractérisées par leurs angles d'inclinaison respectifs e mesurés entre leur direction de propagation V et l'axe Z, ou encore des ondes divergentes émises comme si elles provenaient de différents points de l'espace. L'exemple des ondes planes sera considéré dans ce qui suit. Les ondes incidentes sont généralement des impulsions de moins d'une microseconde, typiquement environ 1 à 10 cycles de l'onde ultrasonore à la fréquence centrale. Les tirs d'ondes incidentes peuvent être séparés les uns des autres par exemple d'environ 50 à 200 microsecondes. Chaque onde incidente rencontre dans le milieu 1 des diffuseurs qui réverbèrent l'onde incidente. L'onde ultrasonore réverbérée est captée par les transducteurs Ti du réseau. Le signal ainsi capté par chaque transducteur Ti provient de l'ensemble du milieu 1, puisque l'onde incidente n'est pas focalisée en émission. De même, le capteur électrique El capte un signal électrique E(t) pendant la propagation de l'onde ultrasonore incidente, et ce signal électrique résulte de l'interaction entre l'onde incidente et le milieu 1 à imager sur toute la ligne représentée par le front d'onde, à chaque instant de mesure. Les signaux réverbérés captés par les n transducteurs Ti sont alors numérisés par les convertisseurs analogique- digital correspondants A/Di et mémorisés dans les mémoires tampon correspondantes Bi, tandis que le signal électrique est numérisé par le convertisseur analogique-digital A/De et mémorisé dans la mémoire tampon correspondante Be. Les signaux ainsi mémorisés dans les mémoires tampon après chaque tir incidents seront appelées ci-après données brutes. Ces données brutes consistent en n+1 signaux temporels bruts RFrawl,i(t) et Erawl(t) captés respectivement par les transducteurs Ti et le capteur électrique El après le tir 1 d'onde ultrasonore incidente. Après chaque tir 1 d'onde incidente, les signaux mémorisés dans les mémoires tampons Bi - Be sont transférés dans la mémoire 9 du processeur de signal 10 aux fins de traitement par ce processeur. A la fin de l'étape (a), la mémoire 9 contient donc N matrices (vecteurs) de n+1 signaux bruts. Etape (b): détermination d'une image du milieu comprenant une cartographie de valeurs électriques : Deux méthodes seront expliquées ci-après pour réaliser cette étape (b). bl) Première méthode : Synthèse de données cohérentes : A partir des N matrices de données brutes, un nombre M de matrices (vecteurs) de données cohérentes synthétiques 5 est calculé par le processeur 8, respectivement en M points Pk(x, z) du champ d'observation (k étant un entier compris entre 1 et M et x, z étant les coordonnées du point Pk sur les axes X, Z). Chacune de ces M vecteurs de données cohérentes synthétiques comporte n signaux temporels 10 RFcoherentk,i(t) correspondant aux signaux qui seraient captés respectivement par les transducteurs Ti si les transducteurs émettaient une onde incidente focalisée au point Pk. Les matrices de données cohérentes peuvent être 15 obtenues par exemple en supposant une vitesse de propagation c homogène dans tout le milieu 1 pour les ondes de compression ultrasonores, selon le principe expliqué notamment dans le document EP2101191 ou dans l'article de Montaldo et al. « Coherent plane-wave compounding for very 20 high frame rate ultrasonography and transient elastography » (IEEE Trans Ultrason Ferroelectr Freq Control 2009 Mar ; 56(3): 489-506). La direction de propagation de l'onde plane correspondant à chaque tir 1 étant connue, et la vitesse de 25 propagation c étant connue, le processeur 8 peut calculer pour chaque point Pk le temps de propagation rec(1,k) de l'onde incidente 1 jusqu'au point Pk, et le temps de propagation r - rec (1,k, i) de l'onde réverbérée depuis le point Pk vers le transducteur Ti, donc le temps de trajet 30 total aller-retour r(1,k, i) = rec(1,PJ i). Le signal acoustique spatialement cohérent pour le transducteur Ti, correspondant au point de focalisation virtuel Pk, est alors calculé selon la formule : RFcoherenC, =IB(1)RFrawhi(z-(1, k, i, j)) ( 1 ) où B(1) est une fonction de pondération pour la contribution de chaque tir 1 d'onde incidente (dans les cas courants, les valeurs B(1) peuvent être toutes égales à 1). Ce signal RFcoherentkii présente une seule valeur pour chaque point Pk. De la même façon, on peut calculer un signal électrique cohérent Ecoherentk : Ecoherentk(t)=IB(1)Erawi(r(1,k,i,j)) (ibis) Cette valeur électrique est celle qui serait mesurée par le capteur électrique El si une onde ultrasonore incidente focalisée en Pk avait émise, particulièrement si on émet un nombre suffisant d'ondes incidentes pour obtenir une image acousto-électrique, par exemple 40 à 100 ondes incidentes pour obtenir une grande résolution d'image. Ces valeurs Ecoherentk sont représentatives des courants électriques aux points Pk, de la même façon que les valeurs électriques captées dans les méthodes connues d'imagerie acousto-électrique susmentionnées, et donnent donc une cartographie des courants électriques dans le champ d'observation. Les matrices de données cohérentes RFcoherentk et éventuellement les valeurs Ecoherentk peuvent ensuite être éventuellement affinées en corrigeant les effets d'aberrations dans le milieu 1, par exemple comme expliqué par exemple dans les documents EP2101191 ou Montaldo et al « Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography » (IEEE Trans Ultrason Ferroelectr Freq Control 2009 Mar ; 56(3): 489- 506). La cartographie des courants électriques peut être présentée à l'écran de l'ordinateur 4, éventuellement en superposition avec une image échographique B mode du milieu 1 ou une autre image dudit milieu 1, notamment une image échographique obtenue à partir des matrices Ecoherentk par formation de voie en réception, comme expliqué par exemple dans le document EP2101191 susmentionné. b2) Deuxième méthode : Transformée de Radon et 5 ondelettes : A partir des signaux électriques Erawk(t), on peut également remonter directement aux valeurs locales de courants électriques aux points Pk, comme il sera expliqué ci-après. 10 Le signal électrique brut Erawk(t) peut être modélisé comme suit : Eraw k = f K pJ (x, y, z)4P(x, y, z)dxdydz (2 ) volume Où : K est une constante d'interaction de l'ordre de 10-9 Pa-1, 15 p est la résistivité du milieu, OP est la variation de pression y est une coordonnée selon un axe Y perpendiculaire au plan (X, Z) et J est la distribution de densité de courant détectée, 20 c'est-à-dire, le produit scalaire du vecteur densité de courant par le vecteur de sensitivité des électrodes du capteur électrique El. L'onde ultrasonore émise étant une onde plane impulsionnelle, AP(x,y,z) peut être paramétrisé en fonction 25 de l'angle d'émission 0 et le temps t. En ignorant la direction Y, nous avons: AP(x,z) = OP(-qsin0 + ctcose, qcos0 + ctsin0), où q et et sont des coordonnées respectivement selon la direction du front d'onde F et selon la direction de 30 propagation V. En considérant l'onde ultrasonore émise comme une impulsion de Dirac, c'est-à-dire une impulsion infiniment courte, le signal acousto-électrique devient: Eraw k ±sc J(-q. sin â + ct. cos 6 , q. cos + ct. sin 8)dq ( 3 ) Kp Ou de façon équivalente : Erawk = RI (B, ct) = J (x , z).S (x. sin B + z. cos B - ct)dxdz ( 4) où R[J] est la transformée de Radon. En pratique, l'onde incidente n'est pas une impulsion de Dirac mais un signal impulsionnel de bande de fréquence finie, ce qui résultera en une convolution par rapport à la variable ct de la transformée de Radon : Eraw k Kp (8 ,t) =W (ct) ® RI (8 ,ct) (5) où W(ct) est la forme d'onde émise et ® est le produit de convolution. Par exemple, une émission ultrasonore typique produit le noyau (kernel) de convolution suivant: (6) où n et m peuvent être ajustés à l'intérieur de la bande de fréquence du transducteur. Ce noyau de convolution est équivalent à une transformée en ridgelettes (« ridgelet transform ») [E. J. Candes, "Ridgelets: theory and applications," Stanford University, 1998] de la 20 distribution de la densité de courant. En pratique, m = n et ce noyau de convolution devient une décomposition en ridgelettes avec les paramètres suivant: a = nÀ, b = ct et O. La décomposition en ridgelettes possède plusieurs 25 propriétés mathématiques comme une relation de Parseval- Plancherel, une formule de reconstruction, une représentation parcimonieuse d'objets à variation lente loin de discontinuités linéaires, et peut être exprimée comme une composition d'une transformation en ondelettes et 30 de la transformée de Radon. Plus spécifiquement, en notant les transformées en Kp ondelettes et en ridgelettes par WT[.], et RT[.] , respectivement, il peut être démontré que Erawk = WT[R(J)] (7) Les inversions de la transformée en ondelettes et de 5 la transformée de Radon sont des problèmes bien connus. En effet, il existe des inversions exactes, respectivement WT-1 et R-1 pour ces deux transformées et nous avons donc: 10 1 J=R-d'abord, inverser courant les ondes Eraw ( 8 ) . 15 ce après WT -1( k) l'inversion se produit en deux étapes: la transformée en ondelettes WT, puis de Radon R. ainsi une cartographie de la densité de le champ d'observation (zone balayée par à l'intérieur du milieu à imager, et très rapide, ce qui permet de une acquisition En On Kp pratique, inverser la transformée obtient dans tout incidentes) suivre en temps réel des phénomènes électriques très rapides, en obtenant un véritable film de la propagation des impulsions électriques. Par ailleurs, il est souhaitable de maximiser le 20 rapport signal à bruit (SNR), la résolution et le débit d'images. Une approche consiste à émettre les ondes incidentes sous formes d'impulsions les plus courtes possible, ce qui optimise la résolution. Cependant, cela correspond à une 25 faible énergie émise et donc, à un faible SNR. Autrement, il est aussi possible de diviser la bande de fréquence en sous-bandes correspondant à des émissions plus longues (et donc plus énergétiques). Il en résulte, en théorie, une augmentation du SNR mais une diminution du 30 débit d'images (puisqu'il faut plusieurs émissions pour former une image). Finalement, une troisième approche consiste en l'émission d'un « chirp » qui peut être utilisé pour faire Kp de la compression d'impulsion. Cette approche permet de maximiser le SNR tout en conservant le débit d'images. Le SNR peut aussi être amélioré en limitant l'effet du bruit. Puisque la transformée en ridgelettes est une base parcimonieuse qui représentera la distribution de densité de courant avec un petit nombre de grands coefficients et un grand nombre de petits coefficients, le débruitage peut être obtenu simplement en appliquant un seuillage sur les signaux obtenus. Une première approche consiste en un seuillage éliminant les 'petits' coefficients. Autrement, il est aussi possible d'utiliser la physique du problème. Par exemple, les coefficients contenant principalement du bruit peuvent être identifiés en performant l'intercorrélation sur des fenêtres de signaux reçus pour deux émissions de polarités opposées. De plus ces signaux peuvent être soustrais pour éliminer les artefacts systémiques. Plusieurs techniques existent pour l'inversion de la transformée de Radon. La plus commune est probablement la rétroprojection filtrée, qui consiste en l'application d'un filtre de type rampe avant la rétroprojection (correspondant à la formation de voies). Pour éviter cette étape qui augmente le niveau de bruit, il est aussi possible d'émettre les ondes ultrasonores incidentes sous forme d'une impulsion adaptée qui inclura ce filtre. D'autres stratégies comme l'acquisition comprimée (« compressed sensing ») sont aussi bien adaptées. Par ailleurs, on peut calculer les matrices RFcoherentk comme expliqué dans la méthode bl) susmentionnée pour former en outre une image échographique bidimensionnelle (B mode) du champ d'observation, par formation de voie en réception, comme expliqué par exemple dans le document EP2101191 susmentionné. Cette image échographique B mode (ou une autre image 35 échographique ou non) du champ d'observation peut être éventuellement superposée à la cartographie des valeurs électriques déterminée précédemment, et on peut afficher sur l'écran de l'ordinateur 4 à la fois l'image échographique du milieu et la cartographie des courants électriques. REVENDICATIONS1. Procédé d'imagerie acousto-électrique, comprenant : (a) une étape de mesure au cours de laquelle on fait émettre par un ensemble (2) de transducteurs Ti, dans un champ d'observation d'un milieu (1) à imager, un nombre N au moins égal à 2 d'ondes ultrasonores incidentes 1 non focalisées dans le champ d'observation et ayant des fronts d'ondes différents, et on fait capter par au moins un capteur électrique (El) en contact avec le milieu à imager, des signaux électriques bruts Erawl(t) respectivement pendant la propagation des ondes incidentes 1, b) une étape de formation d'image, au cours de 15 laquelle on détermine à partir des signaux électriques bruts Erawl(t) obtenus à l'étape (a), une image du milieu comprenant une cartographie de courants électriques. 2. Procédé selon la revendication 1, dans lequel au cours de l'étape b) on détermine, au moins à partir des 20 N signaux électriques bruts Erawl(t), pour un nombre M de points de focalisation fictifs Pk dans champ d'observation, des valeurs électriques Ecoherentk correspondant chacune au signal électrique qui aurait été capté si une onde ultrasonore focalisée au point Pk avait été émise par 25 lesdits transducteurs. 3. Procédé selon la revendication 1, dans lequel au cours de l'étape (b), on applique aux signaux électriques bruts Erawl(t) une transformée en ondelettes inverse WT-1, 30 puis une transformée de Radon inverse R-1. 4. Procédé selon la revendication 1 ou la revendication 2 ou la revendication 3, dans lequel au cours de l'étape (b), on superpose à la cartographie de courants électriques, une image échographique du milieu réalisée 35 avec l'ensemble (2) de transducteurs. . Procédé selon la revendication 4, dans lequel au cours de l'étape (a), on fait capter par les transducteurs Ti des signaux acoustiques RFrawl,i(t) représentatifs d'ondes ultrasonores réverbérées par le milieu 5 respectivement à partir des ondes incidentes 1, au cours de l'étape (b), on détermine à partir des N ensembles de signaux RFrawl,i(t) captés, M signaux acoustiques cohérents RFcoherentk,i(t) correspondant aux signaux acoustiques qui auraient été reçus par les 10 transducteurs Ti si une onde ultrasonore focalisée au point Pk avait été émise par lesdits transducteurs, et on calcule l'image échographique du milieu à partir des signaux acoustiques cohérents. 6. Procédé selon la revendication 5, dans lequel au 15 cours de l'étape (b), on détermine l'image échographique par formation de voies à partir des signaux acoustiques cohérents. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le milieu à imager est un tissu 20 humain ou animal. 8. Dispositif pour la mise en oeuvre d'un procédé d'imagerie acousto-électrique selon l'une quelconque des revendications précédentes, comprenant un ensemble de transducteurs Ti, moins un capteur électrique (El), et des 25 moyens de commande et traitement (8, 10, 4) adaptés pour : (a) faire émettre par un ensemble (2) de transducteurs Ti, dans un milieu (1) à imager, un nombre N d'ondes ultrasonores incidentes non focalisées 1 ayant des fronts d'ondes différents, et faire capter par au moins un 30 capteur électrique en contact avec le milieu à imager, des signaux électriques bruts Erawl(t) respectivement pendant la propagation des ondes incidentes 1, (b) déterminer à partir des signaux électriques bruts Erawl(t), une image du milieu comprenant une 35 cartographie de courants électriques.