FR 3026493 A1 20160401 FR 1459136 A 20140926 Procédé et dispositif d'imagerie acoustique. DOMAINE DE L'INVENTION La présente invention est relative aux procédés et dispositifs d'imagerie acoustique. ARRIERE-PLAN DE L'INVENTION Le document EP2101191 décrit un procédé d'imagerie acoustique, dans lequel on fait successivement émettre dans un milieu à imager, par un réseau de transducteurs acoustiques, une série de n ondes incidentes E±(t) acoustiques planes ou divergentes, les n onde incidentes étant propagées dans le milieu à imager et générant respectivement n ondes réverbérées, et on fait capter par le réseau de transducteurs des signaux représentatifs des ondes réverbérées. RESUME DE L'INVENTION Les procédés de ce type donnent satisfaction, mais présentent l'inconvénient que la sensibilité des images obtenues diminue rapidement avec la profondeur du fait de l'atténuation des ondes ultrasonores dans le milieu à imager, notamment des tissus biologiques. Une solution à ce problème est d'augmenter l'amplitude des ondes ultrasonores émises. Or, cette amplitude ne peut pas être augmentée au-delà de certains seuils, dépendant des transducteurs ultrasonores et / ou de l'électronique qui commande les transducteurs et / ou des normes de sécurité en vigueur pour ce qui concerne les applications médicales. D'autres solutions basées sur un codage temporel ou spatial des signaux émis ont été proposées, mais ces solutions nécessitent l'utilisation d'électroniques d'émissions complexes et couteuses, qui ne peuvent en pratique pas être utilisés sur des appareils d'imagerie échographique classiques. La présente invention a notamment pour but de 35 pallier ces inconvénients et de proposer un procédé d'imagerie acoustique permettant une sensibilité améliorée sans perte de cadence d'imagerie, sans augmenter notablement la complexité ou le coût des appareils d'imagerie acoustique sur lesquels il est mis en oeuvre. A cet effet, la présente invention propose un procédé d'imagerie acoustique comprenant les étapes suivantes : (a) une étape d'émission / réception au cours de laquelle on fait successivement émettre dans un milieu à imager, par un réseau de transducteurs acoustiques, une série de n ondes acoustiques incidentes E±(t) et on fait capter par le réseau de transducteurs acoustiques, n ondes réverbérées Ri(t) qui sont réverbérées respectivement par le milieu à imager suite à l'émission des n ondes incidentes, n étant un entier naturel au moins égal à 2, lesdites ondes acoustiques incidentes étant obtenues en combinant linéairement n ondes incidentes élémentaires E0±(t) de façon correspondante à la formule : É(0=1-/c.É0(t) (1) où : H, est une matrice carrée d'ordre n (dite matrice de codage), Effl=[Ei(01 est un vecteur ayant pour composantes les n ondes incidentes, EID(04E0i(01 est un vecteur ayant pour composantes les n ondes incidentes élémentaires, lesdites ondes incidentes élémentaires étant soit des ondes planes, soit des ondes divergentes, soit des ondes focalisées en différents points de l'espace présentant respectivement n fronts d'ondes différents, (b) une étape de décodage au cours de laquelle on détermine n ondes réverbérées élémentaires R0(t) en combinant linéairement les n ondes réverbérées Ri(t) captées, de façon correspondante à la formule : liO(t) =2.Hd.ii(t) (2) où : À est une constante non nulle, Hd est une matrice carrée d'ordre n (dite matrice de décodage), les matrices H, et Hd étant telles que Hc.Hd=D, où D est une matrice diagonale d'ordre n dont tous les éléments diagonaux dii sont non nuls, est un vecteur ayant pour composantes les n ondes 10 réverbérées, ilD(t)=W0i(01 est un vecteur ayant pour composantes les n ondes réverbérées élémentaires, (c) une étape de construction d'image, au cours de laquelle on réalise une image synthétique à partir des n 15 ondes réverbérées élémentaires ROi(t). L'invention permet ainsi d'augmenter notablement la sensibilité des images échographiques obtenues, sans avoir besoin d'augmenter le nombre de tirs d'ondes incidentes (et donc sans baisser la cadence d'images), sans avoir besoin 20 d'augmenter la complexité de l'appareil d'échographie, et sans avoir besoin d'augmenter l'amplitude des ondes acoustiques incidentes. En effet, les n ondes réverbérées élémentaires R0±(t) sont représentatives respectivement des signaux qui 25 seraient obtenus en émettant séparément les n ondes incidentes élémentaires E0±(t) puis en captant les ondes réverbérées correspondantes, mais lesdites n ondes réverbérées élémentaires R0±(t) ont une amplitude supérieure à l'amplitude des signaux qui seraient captés en 30 émettant séparément les n ondes incidentes élémentaires E0i(t), ce qui induit une amélioration de la sensibilité. Ce gain en amplitude et donc en sensibilité, est dû au fait que chaque onde réverbérée élémentaire R0(t) combine des informations venant de n tirs d'ondes incidentes, donc est plus riche en information. L'invention est utilisable dans tous les procédés où l'image échographique est synthétisée à partir de n tirs d'ondes planes ou divergentes, notamment tel que décrit dans le document EP2101191, par exemple pour l'imagerie des ondes de cisaillement, l'imagerie des réseaux sanguins (imagerie Doppler ultrasensible ou autre, notamment pour l'imagerie fonctionnelle du cerveau), etc. 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 : - les matrices H, et Hd sont telles que Hc.Hd=D, où D est une matrice diagonale d'ordre n dont tous les éléments diagonaux dii sont supérieurs à 1 ; - n est un entier naturel pair non nul ; - la matrice Hc est une matrice de Hadamard d'ordre n et la matrice Hd est la transposée de la matrice Hc ; - chacune desdites ondes acoustiques est représentée par une pluralité de signaux temporels, chaque signal temporel représentant ladite onde acoustique au niveau d'un transducteur appartenant audit réseau de transducteurs ; - le procédé comporte en outre, après l'étape de 25 décodage, une étape de construction d'image, au cours de laquelle on réalise une image synthétique à partir des ondes réverbérées élémentaires ROi(t). Par ailleurs, l'invention a également pour objet un dispositif pour la mise en oeuvre d'un procédé tel que 30 défini ci-dessus, comprenant: (a) des moyens d'émission adaptés pour faire successivement émettre dans un milieu à imager, par un réseau de transducteurs acoustiques, une série de n ondes acoustiques incidentes Ei(t), lesdites ondes 35 acoustiques incidentes étant obtenues en combinant linéairement n ondes incidentes élémentaires E0±(t) de façon correspondante à la formule : É(0=1--/c.É0(t) (1) où : H, est une matrice carrée d'ordre n, Effl=[Ei(01 est un vecteur ayant pour composantes les n ondes incidentes, É0(04E0i(01 est un vecteur ayant pour composantes les n ondes incidentes élémentaires, lesdites ondes incidentes élémentaires étant soit des ondes planes, soit des ondes divergentes, soit des ondes focalisées en différents points de l'espace présentant respectivement n fronts d'ondes différents, (b) des moyens de réception adaptés pour capter par le réseau de transducteurs acoustiques, successivement n ondes réverbérées Ri(t) qui sont réverbérées respectivement par le milieu à imager suite à l'émission des n ondes incidentes, (c) des moyens de décodage adaptés pour déterminer n ondes réverbérées élémentaires R0±(t) en combinant linéairement les n ondes réverbérées Ri(t) captées, de façon correspondante à la formule : Ii0(t)=Hdli(t) (2) où : Hd est une matrice carrée d'ordre n, les matrices H, et Hd étant telles que Hc.Hd=D où D est une matrice diagonale d'ordre n dont tous les éléments diagonaux dii sont non nuls, est un vecteur ayant pour composantes les n ondes 30 réverbérées, Ii0(04R0i(01 est un vecteur ayant pour composantes les n ondes réverbérées élémentaires, (d) des moyens de construction d'image, adaptés pour réaliser une image synthétique à partir des n ondes réverbérées élémentaires ROi(t). 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 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, - la figure 2 est un schéma bloc d'une partie du dispositif de la figure 1, - la figure 3 est une série de chronogrammes représentant les ondes incidentes élémentaires composant les différentes ondes acoustiques incidentes émises successivement pour réaliser une image synthétique, - les figues 4a a 4d sont des images Doppler ultrarapides montrant la vascularisation d'un cerveau, les figures 4a et 4c étant réalisées dans deux plans différents par imagerie synthétique sans la présente invention, et les figures 4b et 4d étant réalisées respectivement dans les mêmes plans par imagerie synthétique avec la présente invention, et - les figures 5a et 5b sont des images élastographiques montrant la propagation d'une onde de cisaillement dans des tissus humains, obtenues par imagerie synthétique respectivement avec et sans l'invention. DESCRIPTION DETAILLEE Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. Les figures 1 et 2 montrent un exemple de dispositif adapté pour réaliser une imagerie échographique 35 synthétique d'un milieu 1 à imager, par exemple des tissus d'un sujet humain ou animal à étudier, par émission et réception d'ondes ultrasonores de compression, par exemple de fréquences comprises entre 2 et 40 MHz, notamment entre 2 et 10 MHz. L'émission et réception d'ondes ultrasonores peut se faire au moyen d'un réseau 2 de p transducteurs ultrasonores 2a (Ti-T,), comprenant par exemple quelques centaines de transducteurs 2a, ce réseau 2 pouvant être par exemple une barrette de transducteurs (réseau 1D) adapté pour réaliser une image bidimensionnelle (2D) ou encore un réseau à deux dimensions adapté pour réaliser une image tridimensionnelle (3D). Le dispositif d'imagerie comporte également, par exemple : - un circuit électronique 3 commandant le réseau de transducteurs 2 et adapté pour faire émettre des ondes ultrasonores par le réseau de transducteurs et acquérir les signaux captés par ce réseau de transducteurs ; - un ordinateur 4 ou similaire pour commander le 20 circuit électronique 3 et visualiser les images ultrasonores obtenues à partir desdits signaux captés. Comme représenté sur la figure 2, le circuit électronique 3 peut comprendre par exemple : p convertisseurs analogique / digital 11 (A/Di) 25 connectés individuellement aux p transducteurs Ii-In du réseau 2 de transducteurs, p mémoires tampon 12 (Bi) respectivement connectées aux p convertisseurs analogique / digital 11, - une unité centrale 13 (CPU) communiquant avec 30 les mémoires tampon 12 et l'ordinateur 4 et adapté au traitement des signaux envoyés vers le réseau de transducteurs 2 et reçus dudit réseau de transducteurs, - une mémoire 14 (MEM) connectée à l'unité centrale 13. 35 Le dispositif d'imagerie représenté sur les figures 1 et 2 est adapté pour réaliser une imagerie ultrasonore synthétique du milieu 1 comme décrit notamment dans le document EP2101191. On peut distinguer plusieurs étapes principales 5 dans cette méthode d'imagerie : (a) une étape d'émission / réception (b) une étape de décodage, (c) une étape de construction d'image. Ces différentes étapes vont maintenant être 10 détaillées. Les étapes d'émission et de décodage font appel à des matrices de codage H, et décodage Hd, qui seront décrites préalablement. Ces matrices de codage H, et décodage Hd sont matrices carrées d'ordre n, où n est un entier naturel 15 au moins égale à 2, préférentiellement pair. n peut par exemple être compris entre 2 et 100, avantageusement entre 4 et 20. Les matrices de codage H, et de décodage Hd sont telles que Hc.Hd=D où D est une matrice diagonale 20 d'ordre n dont tous les éléments diagonaux dil sont non nuls. Avantageusement, D peut être une matrice diagonale d'ordre n dont tous les éléments diagonaux dil sont positifs, avantageusement supérieurs à 1. 25 Plus avantageusement encore, la matrice Hc peut être proportionnelle à une matrice de Hadamard Hn d'ordre n et la matrice Hd est la transposée de la matrice Hc. Autrement dit : 30 H, = 2.H' (0) et Hd=t1-1, (0'), Où À est une constante non nulle. Des exemples de matrices de Hadamard Hn sont donnés ci-après, pour n=2 et n=4. 35 Pour n=2 : 9 H2= Pour n=4 : H4= Dans le cas où la matrice de codage est une matrice de Hadamard (Hc=Hfl), la matrice diagonale D susmentionnée est égale à n fois la matrice identité (autrement dit tous les éléments diagonaux dii de D sont égaux à n). (a) étape d'émission / réception (al) Emissions : Au cours de l'étape d'émission / réception, l'unité centrale 13 fait successivement émettre dans un milieu 1, par le réseau 2 de transducteurs, une série de n ondes acoustiques incidentes Ei(t), notamment des ondes ultrasonores comme décrit précédemment. Les ondes acoustiques 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. Chacune des ondes acoustiques incidentes est obtenue en combinant linéairement n ondes incidentes élémentaires E0±(t) qui sont soit des ondes planes, soit des ondes divergentes, soit des ondes focalisées en différents points de l'espace présentant respectivement n fronts d'ondes différents. Cette combinaison linéaire se fait de façon correspondante à la formule : É(t) = I I c.É0(t) (1) où : Effl=[Ei(01 est un vecteur ayant pour composantes les n ondes incidentes, a)(0=WOM est un vecteur ayant pour composantes les n ondes incidentes élémentaires. Les ondes acoustiques incidentes peuvent être calculées à l'avance et mémorisées par exemple dans la mémoire 14. Avantageusement, chacune des n ondes acoustiques incidentes E±(t) peut être représentée par une pluralité de signaux temporels E±(t) représentant ladite onde acoustique incidente telle qu'émise respectivement par les différents transducteurs Tj dudit réseau 2 de transducteurs. Ainsi, E±(t) peut être représentée par un vecteur Ei(t)=[Eij(t)]. De même, chacune des n ondes incidentes élémentaires E0(t) peut être représentée par une pluralité de signaux temporels E0±(t) représentant ladite onde incidente élémentaire respectivement au niveau des différents transducteurs Tj dudit réseau 2 de transducteurs. Ainsi, E0(t) peut être représentée par un vecteur E0i(t)=[E0ij(t)]. La figure 3 illustre un cas particulier avec n = 4, où les ondes acoustiques incidentes E±(t) sont des combinaisons linéaires de quatre ondes incidentes élémentaires qui sont des ondes planes E0±(t) d'inclinaisons différentes, dont les fronts d'ondes sont représentés selon l'axe X parallèle au réseau de transducteurs et en fonction du temps t. Dans cet exemple, la matrice de codage utilisée est la matrice de Hadamard H4 susmentionnée. Les ondes planes affectées d'un coefficient linéaire -1 (-E0i(t)) sont représentées en traits mixtes et les ondes planes affectées d'un coefficient linéaire +1 (+E0i(t)) sont représentées en traits pleins. Comme représenté sur les vues de détail en haut de la figure 3 pour les ondes E1(t) et E2(t), les formes d'ondes des ondes incidentes élémentaires E0±(t) peuvent être des sinusoïdes amorties durant par exemple environ 3 cycles, les ondes planes affectées d'un coefficient linéaire -1 (-E0i(t)) ayant simplement une forme d'onde opposée à E0i(t). (a2) Réceptions : Chacune des ondes acoustiques incidentes E±(t) se propage dans le milieu 1 qui renvoie une onde réverbérée correspondante Ri(t). Avant l'émission de l'onde acoustique incidente suivante, le réseau 2 de transducteurs capte l'onde réverbérée Ri(t) qui est transmise à l'unité centrale 13 et mémorisée dans la mémoire 14. L'onde réverbérée Ri(t) est mémorisée sous la forme d'une pluralité de signaux temporels R±(t) représentant ladite onde réverbérée telle que captée respectivement par les différents transducteurs Ti du réseau 2 de transducteurs. Ainsi, Ri(t) peut être représentée par un vecteur Ri(t)=[Rii(t)]. (b) étape de décodage Après émission des n ondes acoustiques incidentes E±(t) et réception des n ondes réverbérées Ri(t), l'unité centrale 13 procède à une étape de décodage au cours de laquelle on détermine n ondes réverbérées élémentaires R0(t) en combinant linéairement les n ondes réverbérées Ri(t) captées, de façon correspondante à la formule : Ii0(t)=Hdli(t) (2) où : est un vecteur ayant pour composantes les n ondes réverbérées, Ii0(t)=W0i(01 est un vecteur ayant pour composantes les n ondes réverbérées élémentaires. Les ondes réverbérées élémentaires RO(t) sont déterminées sous la forme d'une pluralité de signaux temporels R0±(t) représentant ladite onde réverbérée élémentaire au niveau des différents transducteurs Ti du réseau 2 de transducteurs. Ainsi, RO(t) peut être représentée par un vecteur ROi(t)=[R0ii(t)], ou plus couramment par une matrice exactement ROi(t)=[ROijk(tk)] après échantillonnage temporel en K instants tk. (c) étape de construction d'image Après l'étape de décodage (b), l'unité centrale 13 procède à une étape de construction d'image, au cours de laquelle on réalise une image synthétique à partir des n ondes réverbérées élémentaires ROi(t), par exemple comme enseigné dans le document EP2101191 ou dans l'article de 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). Les ondes réverbérées élémentaires ROi(t), généralement présentées sous forme matricielle ROi(t)=[ROijk(tk)] comme indiqué précédemment, constituent les matrices de données brutes de départ (ou « données brutes RF ») des procédés décrits dans les documents susmentionnés. Plus particulièrement, à partir des n matrices R0±(t) de données brutes RF, l'unité centrale 13 calcule des matrices de données cohérentes synthétiques RF, respectivement en M points Pm(x, z) du milieu à imager 1 (dits points de focalisation virtuels), m étant un entier compris entre 1 et M, x l'abscisse du point Pm sur un axe X parallèle au réseau 2 de transducteurs dans le cas d'un réseau 2 linéaire de transducteurs (x étant remplacé par x,y dans le cas d'un réseau 2 bidimensionnel) et z est la profondeur dans le milieu 1, selon un axe Z perpendiculaire à l'axe X. Chacune de ces M matrices de données cohérentes synthétiques RF comporte p signaux temporels RFcoherentm,j(t) correspondant aux signaux qui seraient captés respectivement par les p transducteurs Tj si les transducteurs émettaient une onde incidente focalisée au point Pm. Les matrices de données RF cohérentes peuvent être 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 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 correspondant à chaque one incidente élémentaire étant connue, et la vitesse de propagation c connue, l'unité centrale 13 peut calculer pour point Pm le temps de propagation Tec(i,m) de plane E0 (t) étant chaque l'onde incidente élémentaire i jusqu'au point Pm, et le temps de propagation (i, m, j) de l'onde réverbérée élémentaire depuis le point Pm vers le transducteur Tj, donc le temps de trajet total aller-retour 2(i,m,j) = i-ec(irPra) re c(r Pm r ) - Le signal spatialement cohérent pour le transducteur Tj, correspondant au point de focalisation virtuel Pm, est alors calculé selon la formule : RFcoherentn1(t)=IB(i)RFraw11(4-(i,m,j)) (3) où B(i) est une fonction de pondération pour la contribution de chaque onde incidente élémentaire (éventuellement, les coefficients B(i) peuvent être tous égaux à 1) Les matrices de données cohérentes REcoherentmi peuvent ensuite être éventuellement affinées en corrigeant les effets d'aberrations dans le milieu 1, par exemple comme expliqué dans les documents susmentionnés EP2101191 ou Montaldo et al. L'unité centrale 13 peut alors procéder à une formation de voie en réception classique pour obtenir une image ultrasonore comme expliqué dans les documents susmentionnés EP2101191 ou Montaldo et al., et / ou procéder à des traitements plus élaborés (par exemple calcul d'une image Doppler ou Doppler couleur), connus en soi. Le procédé de construction d'image est similaire lorsque le réseau 2 de transducteurs est bidimensionnel, ou lorsque les ondes incidentes élémentaires et les ondes réverbérées élémentaires sont divergentes. Dans tous les cas, l'image obtenue présente une sensibilité nettement accrue par rapport à une image obtenue par imagerie synthétique classique par émission d'ondes planes (telle que décrite par les documents susmentionnés EP2101191 ou Montaldo et al). Ainsi, les figues 4a et 4c des images Doppler ultrarapides montrant la vascularisation d'un cerveau de rat, réalisées dans deux plans différents par imagerie synthétique avec émission d'ondes planes (sans la présente invention), tandis que les figures 4b et 4d sont réalisées respectivement dans les mêmes plans par imagerie synthétique avec la présente invention : la sensibilité nettement accrue des images des figues 4b et 4d est visible à De même, les figures 5a et 5b sont des images élastographiques montrant la propagation d'une onde de cisaillement dans des tissus humains, obtenues par imagerie synthétique respectivement avec et sans l'invention : la sensibilité supérieure de l'image de la figure 5b est visible à l'oeil.5 REVENDICATIONS1. Procédé d'imagerie acoustique comprenant les étapes suivantes : (a) une étape d'émission / réception au cours de laquelle on fait successivement émettre dans un milieu (1) à imager, par un réseau (2) de transducteurs acoustiques, une série de n ondes acoustiques incidentes E±(t) et on fait capter par le réseau de transducteurs acoustiques, n ondes réverbérées Ri(t) qui sont réverbérées respectivement par le milieu à imager suite à l'émission des n ondes incidentes, n étant un entier naturel au moins égal à 2, lesdites ondes acoustiques incidentes étant obtenues en combinant linéairement n ondes incidentes élémentaires E0(t) de façon correspondante à la formule : É(0=1--/c.É0(t) (1) où : H, est une matrice carrée d'ordre n, Effl=[Ei(0] est un vecteur ayant pour composantes les n ondes incidentes, E0(04E0i(01 est un vecteur ayant pour composantes les n ondes incidentes élémentaires, lesdites ondes incidentes élémentaires étant soit des ondes planes, soit des ondes divergentes, soit des ondes focalisées en différents points de l'espace, présentant respectivement n fronts d'ondes différents, (b) une étape de décodage au cours de laquelle on détermine n ondes réverbérées élémentaires R0(t) en combinant linéairement les n ondes réverbérées Ri(t) captées, de façon correspondante à la formule : RO(t) = Hd.R(t) ( 2 ) où : Hd est une matrice carrée d'ordre n, les matrices H, et Hdétant telles que Hc.Hd=D où D est une matrice diagonale d'ordre n dont tous les éléments diagonaux dil sont non nuls, est un vecteur ayant pour composantes les n ondes 5 réverbérées, i0(0=[ROiffl] est un vecteur ayant pour composantes les n ondes réverbérées élémentaires, (c) une étape de construction d'image, au cours de laquelle on réalise une image synthétique à partir des n 10 ondes réverbérées élémentaires ROi(t). 2. Procédé selon la revendication 1, dans lequel les matrices H, et Hd sont telles que Hc.Hd=D, où D est une matrice diagonale d'ordre n dont tous les éléments diagonaux dii sont supérieurs à 1. 15 3. Procédé selon la revendication 1 ou la revendication 2, dans lequel n est un entier naturel pair non nul. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel la matrice Hc est 20 une matrice de Hadamard d'ordre n et la matrice Hd est la transposée de la matrice Hc. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel chacune desdites ondes acoustiques est représentée par une pluralité de 25 signaux temporels, chaque signal temporel représentant ladite onde acoustique au niveau d'un transducteur appartenant audit réseau de transducteurs. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le milieu (1) à 30 imager comprend un tissu biologique vivant humain ou animal. 7. Dispositif d'imagerie acoustique, comprenant: (a) des moyens d'émission (13, 2) adaptés pour faire successivement émettre dans un milieu (1) à imager, parun réseau (2) de transducteurs acoustiques, une série de n ondes acoustiques incidentes Ei(t), lesdites ondes acoustiques incidentes étant obtenues en combinant linéairement n ondes incidentes élémentaires E0(t) de façon correspondante à la formule : É(t) = .E0(t) (1) où : H, est une matrice carrée d'ordre n, Effl=[Ei(0] est un vecteur ayant pour composantes les n ondes incidentes, É0(04E0i(01 est un vecteur ayant pour composantes les n ondes incidentes élémentaires, lesdites ondes incidentes élémentaires étant soit des ondes planes, soit des ondes divergentes, soit des ondes focalisées en différents points de l'espace présentant respectivement n fronts d'ondes différents, (b) des moyens de réception (2, 13) adaptés pour capter par le réseau (2) de transducteurs acoustiques, successivement n ondes réverbérées Ri(t) qui sont réverbérées respectivement par le milieu à imager suite à l'émission des n ondes incidentes, (c) des moyens de décodage (13) adaptés pour déterminer n ondes réverbérées élémentaires R0(t) en combinant linéairement les n ondes réverbérées Ri(t) captées, de façon correspondante à la formule : Ii0(t)=Hdli(t) (2) où : Hd est une matrice carrée d'ordre n, les matrices H, et Hd étant telles que Hc.Hd=D où D est une matrice diagonale 30 d'ordre n dont tous les éléments diagonaux dil sont non nuls, est un vecteur ayant pour composantes les n ondes réverbérées,ilD(t)=W0i(01 est un vecteur ayant pour composantes les n ondes réverbérées élémentaires, (d) des moyens de construction d'image (13), adaptés pour réaliser une image synthétique à partir des n ondes réverbérées élémentaires ROi(t).