La présente invention concerne un procédé de localisation d'au moins une cible dans un milieu électromagnétiquement absorbant. Elle concerne également un programme d'ordinateur et un dispositif correspondants, ainsi qu'une application de ce procédé pour la détection et la localisation de tumeurs dans des tissus biologiques.

L'invention s'applique plus particulièrement à un procédé de localisation d'au moins une cible dans un milieu électromagnétiquement absorbant, comportant les étapes suivantes :

- émission d'au moins un signal excitateur électromagnétique à partir d'au moins une source,

- réception, par au moins un capteur acoustique, d'un signal acoustique résultant de l'émission de ce signal excitateur,

- détection, dans le signal acoustique reçu par le capteur acoustique, d'un instant de réception d'une réponse au signal excitateur, cette réponse résultant d'une perturbation acoustique due à une hétérogénéité électromagnétique de la cible dans le milieu, et

- estimation d'une distance entre la cible et le capteur acoustique à l'aide de cet instant de réception.

Le contexte de l'invention est donc celui de l'imagerie optoacoustique, dite également photoacoustique, ou thermoacoustique, appliquée à la localisation d'hétérogénéités (ou cibles) enfouies dans un milieu caractérisé par des propriétés électromagnétiques et acoustiques différentes de ces hétérogénéités. C'est par exemple le cas d'hétérogénéités telles que des tumeurs dans des tissus biologiques. L'imagerie optoacoustique ou thermoacoustique est ainsi en plein essor actuellement dans le domaine biomédical, en particulier parce qu'elle est réputée non invasive et non ionisante. Son principe est par exemple décrit dans l'article de C. Li et L. Wang, intitulé « Photoacoustic tomography and sensing in biomedicine », publié dans Physics in Medicine and Biology, vol. 54, pages R59-R97, 2009.

Les techniques actuelles d'imagerie optoacoustique ou thermoacoustique exploitent la génération instantanée d'une perturbation acoustique locale due à une hétérogénéité électromagnétique d'une cible dans un milieu électromagnétiquement absorbant et diffusant (si on se place en imagerie optoacoustique) ou diffractant (si on se place en imagerie thermoacoustique), lorsqu'un signal excitateur électromagnétique est émis depuis une source. Plus précisément, l'énergie électromagnétique émise, en étant partiellement absorbée dans le milieu observé, se dissipe thermiquement, ce qui crée des dilatations dans le milieu. En particulier, l'hétérogénéité électromagnétique de la cible crée localement une hétérogénéité de dilatation qui engendre à son tour la perturbation acoustique locale révélatrice de cette hétérogénéité.

Cette perturbation acoustique locale est détectée au niveau d'un capteur acoustique sous la forme d'une réponse avec un retard, par rapport à l'émission du signal excitateur, significatif de la distance entre la cible et le capteur acoustique. Cette distance peut donc être estimée à l'aide de l'instant de réception de la réponse précitée par le capteur acoustique. En effet, l'imagerie optoacoustique ou thermoacoustique utilise la propriété selon laquelle la génération de la perturbation acoustique au niveau de la cible se fait depuis la source à une vitesse voisine de celle de la lumière (vitesse de propagation électromagnétique), ou en tout cas très grande devant celle du son dans le milieu, alors que la propagation de la perturbation acoustique depuis la cible vers le capteur acoustique se fait à la vitesse du son. Par conséquent, en synchronisant la réception avec l'émission et en connaissant la vitesse du son dans le milieu, on déduit simplement la distance entre la cible et le capteur acoustique de la durée s'étant écoulée entre l'émission du signal excitateur et la réception de la réponse à ce signal.

Il est ainsi possible de localiser la cible dans l'espace tridimensionnel sur une sphère dont le centre est le capteur acoustique et dont le rayon est la distance estimée. Mais plusieurs émissions-réceptions sont nécessaires pour localiser plus précisément la cible. Ainsi les technologies actuelles impliquent une certaine complexité de capture et de traitement de signaux pour permettre une localisation de la cible. Comme en outre ces technologies sont généralement intégrées dans un principe plus global d'imagerie tomographique, leur complexité engendre également une plus grande complexité de reconstruction tomographique.

Une solution est de combiner la technologie optoacoustique ou thermoacoustique avec par exemple une technologie purement ultrasonore, pour limiter le nombre d'émissions/réceptions nécessaires. Une telle solution est proposée dans l'article de J. Zalev et al, intitulé « Clinical feasibility study of combined optoacoustic and ultrasonic imaging modality providing coregistered functional and anatomical maps of breast tumors », publié dans Photons Plus Ultrasound: Imaging and Sensing 2012, Proceedings of SPIE, vol. 8223, 2012. Elle reste complexe techniquement puisqu'elle combine deux technologies proches mais distinctes.

Par ailleurs, les technologies actuelles d'imagerie optoacoustique ou thermoacoustique limitent l'exploitation du signal acoustique résultant de l'émission du signal excitateur à l'inhomogénéité électromagnétique de la cible. Pour se concentrer sur cette propriété, en particulier quand le milieu dans lequel se trouve la cible est lui-même acoustiquement inhomogène, il est même enseigné de s'affranchir de ce problème en limitant la bande passante du capteur acoustique ou en effectuant des mesures acoustiques dans une configuration géométrique perpendiculaire à l'illumination du milieu. C'est ce qui est divulgué dans l'article de S. Ermilov et al, intitulé « Laser optoacoustic imaging System for détection of breast cancer », publié dans Journal of Biomédical Optics, vol. 14(2), pages 024007-1 à 024007-14, mars/avril 2009. Pourtant d'autres réponses au signal excitateur sont présentes dans le signal reçu par le capteur acoustique et mériteraient sans doute d'être exploitées plutôt que d'avoir recours à des combinaisons de technologies différentes. Ainsi par exemple, une inhomogénéité électromagnétique existe entre la source d'émission et le milieu, mais elle n'est pas exploitée dans les documents précités de l'art antérieur.

En revanche, dans la demande de brevet internationale WO 201 1/096198 A1 , l'inhomogénéité électromagnétique de la source par rapport au milieu est exploitée, plus précisément son inhomogénéité optique. Elle conduit à détecter, outre la réponse résultant d'une perturbation acoustique due à une hétérogénéité électromagnétique de la cible dans le milieu, une autre réponse résultant d'une perturbation acoustique due à l'hétérogénéité optique de la source par rapport au milieu. Cette autre réponse se propage à la vitesse du son dans le milieu, depuis la source vers le capteur acoustique. En connaissant la distance qui sépare la source du capteur acoustique, on en déduit une estimation de la vitesse du son dans le milieu. Cette solution permet donc de se passer d'une connaissance a priori de la vitesse du son dans le milieu pour estimer la distance séparant la cible du capteur acoustique. Mais là encore, plusieurs émissions-réceptions sont nécessaires pour localiser plus précisément la cible ce qui rend cette technologie complexe.

Il peut ainsi être souhaité de prévoir un procédé de localisation d'au moins une cible dans un milieu qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités.

Il est donc proposé un procédé de localisation d'au moins une cible dans un milieu électromagnétiquement absorbant, comportant les étapes suivantes : - émission d'au moins un signal excitateur électromagnétique à partir d'au moins une source,

- réception, par au moins un capteur acoustique, d'un signal acoustique résultant de l'émission de ce signal excitateur,

- détection, dans le signal acoustique reçu par le capteur acoustique, d'un premier instant de réception d'une première réponse au signal excitateur, cette première réponse résultant d'une perturbation acoustique due à une hétérogénéité électromagnétique de la cible dans le milieu,

- estimation d'une première distance entre la cible et le capteur acoustique à l'aide de ce premier instant de réception,

comportant en outre les caractéristiques suivantes :

- détection, dans le même signal acoustique reçu par le capteur acoustique, d'un second instant de réception d'une seconde réponse au signal excitateur, cette seconde réponse résultant d'une perturbation acoustique due à une hétérogénéité acoustique de la cible dans le milieu,

- estimation d'une seconde distance entre la source et la cible à l'aide de ce second instant de réception, et

- obtention d'une localisation de la cible à partir des première et seconde distances estimées.

En exploitant l'hétérogénéité acoustique de la cible dans le milieu, un procédé de localisation selon l'invention permet d'améliorer très significativement les technologies existantes. En effet, cela conduit à détecter, outre la première réponse résultant d'une première perturbation acoustique due à l'hétérogénéité électromagnétique de la cible dans le milieu, une seconde réponse résultant d'une seconde perturbation acoustique due à l'hétérogénéité acoustique de la cible. La seconde perturbation acoustique est générée au niveau de la source et se propage, depuis la source jusqu'à la cible, à la vitesse du son dans le milieu, puis continue à se propager, depuis la cible jusqu'au capteur acoustique, également à la vitesse du son dans le milieu. Par conséquent, en synchronisant la réception avec l'émission et en connaissant la vitesse du son dans le milieu, on déduit simplement la longueur du trajet source-cible-capteur de la durée s'étant écoulée entre l'émission du signal excitateur et la réception de la seconde réponse précitée. Connaissant en outre la distance entre la cible et le capteur grâce à la première réponse, on en déduit la distance entre la source et la cible à l'aide de l'instant de réception de la seconde réponse. On peut aussi déduire la distance entre la source et la cible directement à partir de la durée écoulée entre les premier et second instants de réception. Ainsi, avec une seule émission, on dispose de deux informations, la distance source-cible et la distance cible-capteur, pour localiser plus précisément la cible dans le milieu. A nombre d'émissions-réceptions constant, on obtient donc une meilleure localisation. En variante, pour une précision de localisation souhaitée de la cible, il devient possible de réduire le nombre d'émissions-réceptions.

De façon optionnelle :

- l'estimation de la distance entre la cible et le capteur acoustique se fait par multiplication d'un premier temps de vol de la première réponse entre la cible et le capteur acoustique avec une valeur prédéterminée de vitesse d'onde acoustique dans le milieu,

- le premier temps de vol est estimé en tant que temps écoulé entre l'instant d'émission du signal excitateur et le premier instant de réception.

De façon optionnelle également :

- l'estimation de la distance entre la source et la cible se fait par multiplication d'un second temps de vol de la seconde réponse entre la source et la cible avec la valeur prédéterminée de vitesse d'onde acoustique dans le milieu,

- le second temps de vol est estimé en tant que différence entre, d'une part, le temps écoulé entre l'instant d'émission du signal excitateur et le second instant de réception, et, d'autre part, le premier temps de vol.

De façon optionnelle également, la localisation de la cible dans l'espace tridimensionnel est obtenue par intersection de deux sphères, la première sphère ayant pour centre le capteur acoustique et pour rayon la première distance, la seconde sphère ayant pour centre la source et pour rayon la seconde distance.

De façon optionnelle également, le signal excitateur émis est un signal optique issu d'une source lumineuse de fréquence et/ou d'intensité modulées et les deux réponses au signal excitateur sont des réponses impulsionnelles résultant de deux perturbations acoustiques dues à une double hétérogénéité optique et acoustique de la cible dans le milieu.

De façon optionnelle également, un procédé de localisation d'au moins une cible selon l'invention peut comporter les étapes suivantes :

- émission de signaux excitateurs électromagnétiques à partir de plusieurs sources réparties en périphérie du milieu dans lequel est située la cible, - réception, par plusieurs capteurs acoustiques répartis en périphérie du milieu dans lequel est située la cible, de signaux acoustiques résultant des émissions de ces signaux excitateurs, et

- reconstruction tomographique d'une image de la cible dans le milieu dans lequel elle est située, à l'aide des localisations obtenues de la cible.

De façon optionnelle également, les sources et les capteurs acoustiques sont répartis régulièrement sur un cercle autour du milieu dans lequel est située la cible.

Il est également proposé une application d'un procédé de localisation d'au moins une cible selon l'invention pour la détection et la localisation de tumeurs dans des tissus biologiques.

Il est également proposé un programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions pour l'exécution des étapes d'un procédé de localisation d'au moins une cible selon l'invention, lorsque ledit programme est exécuté sur un ordinateur.

Il est également proposé un dispositif de localisation d'au moins une cible dans un milieu électromagnétiquement absorbant, comportant :

- au moins une source d'émission d'au moins un signal excitateur électromagnétique,

- au moins un capteur acoustique d'un signal acoustique résultant de l'émission de ce signal excitateur,

- un calculateur conçu pour détecter, dans le signal acoustique reçu par le capteur acoustique, un premier instant de réception d'une première réponse au signal excitateur, cette première réponse résultant d'une perturbation acoustique due à une hétérogénéité électromagnétique de la cible dans le milieu, et pour estimer une première distance entre la cible et le capteur acoustique à l'aide de ce premier instant de réception, dans lequel le calculateur est en outre conçu pour :

- détecter, dans le même signal acoustique reçu par le capteur acoustique, un second instant de réception d'une seconde réponse au signal excitateur, cette seconde réponse résultant d'une perturbation acoustique due à une hétérogénéité acoustique de la cible dans le milieu,

- estimer une seconde distance entre la source et la cible à l'aide de ce second instant de réception, et - fournir une localisation de la cible à partir des première et seconde distances estimées.

L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels :

- la figure 1 représente schématiquement la structure générale d'un dispositif de localisation de cible dans un milieu électromagnétiquement absorbant, selon un mode de réalisation de l'invention,

- la figure 2 illustre les étapes successives d'un procédé mis en œuvre par le dispositif de la figure 1 ,

- la figure 3 illustre un exemple d'application du procédé de la figure 2,

- la figure 4 représente schématiquement la structure générale d'un exemple de système d'émission/réception pour le dispositif de la figure 1 , et

- la figure 5 illustre l'intégration du procédé de la figure 2 dans un schéma plus global de reconstruction tomographique par imagerie optoacoustique ou thermoacoustique.

Dans ce qui suit, on se place dans le domaine de l'optoacoustique (ou photoacoustique) de façon purement illustrative et non limitative.

Le dispositif 10 représenté schématiquement sur la figure 1 est ainsi un dispositif d'imagerie optoacoustique permettant la localisation d'au moins une cible T située dans un milieu 12 optiquement absorbant et diffusant (puisque l'on se place en imagerie optoacoustique). La cible T, par exemple une tumeur dans un milieu observé 12 constitué d'un tissu biologique, présente une double hétérogénéité optique et acoustique dans le milieu 12 qui est par ailleurs homogène du point de vue de ces deux caractéristiques. Cette double hétérogénéité de la cible T est exploitée par le dispositif 10 pour permettre sa localisation, et éventuellement sa visualisation.

A cette fin, le dispositif de localisation 10 comporte un émetteur 14 comportant au moins une source d'émission d'au moins un signal excitateur électromagnétique E, plus précisément optique dans cet exemple. Cet émetteur 14 est conçu pour illuminer optiquement le milieu 12 et engendrer des perturbations acoustiques dues à la double hétérogénéité optique et acoustique de la cible T dans son milieu 12 et à l'hétérogénéité optique de la source d'émission par rapport au milieu 12. Ainsi, de façon connue en soi et selon les principes généraux de l'imagerie optoacoustique, l'émetteur 14 est une source lumineuse de fréquence et/ou d'intensité modulées, par exemple un émetteur laser puisé de type Nd :Yag à 1064 nm, de largeur temporelle d'impulsion laser d'environ 3 ns, présentant un taux de répétition de 10 Hz et une puissance inférieure à 20 mJ/cm ² si le milieu 12 est un tissu biologique. Compte tenu des propriétés optiques du milieu 12, il se crée instantanément, c'est-à-dire à une vitesse proche de celle de la lumière ou en tout cas très grande devant celle du son dans le milieu 12, une carte CPinit de distribution de pression initiale dans le milieu 12 révélatrice des dilatations mentionnées précédemment. Cette distribution de pression initiale engendre ensuite une onde de pression acoustique qui se propage à la vitesse du son selon les propriétés acoustiques du milieu 12. Si des éléments dans le milieu, tels que la cible T, sont optiquement et acoustiquement inhomogènes, des perturbations acoustiques correspondant à ces hétérogénéités sont détectables par réceptions de signaux acoustiques résultant de l'émission du signal excitateur E.

Le dispositif de localisation 10 comporte donc en outre un récepteur 16 comportant au moins un capteur acoustique d'un signal acoustique R résultant de l'émission du signal excitateur E. Le ou les capteurs acoustiques du récepteur 16 sont bien connus de l'homme du métier et ne seront pas détaillés. On peut notamment se reporter aux documents de l'état de la technique cités précédemment.

Le dispositif de localisation 10 comporte en outre une carte d'acquisition 18 permettant un prétraitement du signal reçu R, par exemple un filtrage et une numérisation, ainsi qu'une synchronisation du signal prétraité à l'aide d'un signal « sync » reçu de l'émetteur 14. Le signal reçu R comporte en effet un certain nombre de réponses impulsionnelles acoustiques à chaque impulsion optique émise par l'émetteur 14, chacune de ces réponses impulsionnelles étant, comme cela sera détaillé par la suite, due à une hétérogénéité électromagnétique, plus précisément optique dans cet exemple, dans le milieu 12. Il est donc important de synchroniser le signal reçu R par rapport au signal excitateur E pour permettre une estimation des instants de réception des réponses impulsionnelles acoustiques détectables dans le signal reçu R par rapport à l'instant d'émission du signal excitateur E.

Le dispositif de localisation 10 comporte en outre un calculateur 20 programmé pour commander l'émetteur 14 à l'aide d'un signal C en fonction d'une illumination souhaitée du milieu 12 et pour analyser le signal reçu par le récepteur 16 et prétraité par la carte d'acquisition 18. Ce signal reçu et prétraité est noté Rs en sortie de la carte d'acquisition 18.

Le calculateur 20 tel que représenté schématiquement sur la figure 1 , comporte une unité de traitement 22 associée de façon classique à une mémoire 24 (par exemple une mémoire RAM). L'unité de traitement 22 peut par exemple être mise en œuvre dans un dispositif informatique tel qu'un ordinateur classique comportant un processeur associé à une ou plusieurs mémoires pour le stockage de fichiers de données et de programmes d'ordinateurs. L'unité de traitement 22 peut alors elle-même être considérée comme formée de ce processeur associé à la mémoire 24 agissant en tant que mémoire de stockage des instructions que le processeur exécute sous forme de programmes d'ordinateurs.

La mémoire 24 telle qu'illustrée sur la figure 1 comporte ainsi fonctionnellement deux programmes d'ordinateurs ou deux fonctions d'un même programme d'ordinateur 26 et 28. On notera en effet que les programmes d'ordinateurs 26 et 28 sont présentés comme distincts, mais cette distinction est purement fonctionnelle. Ils pourraient tout aussi bien être regroupés selon toutes les combinaisons possibles en un ou plusieurs logiciels. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique mettant en œuvre le calculateur 20 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d'ordinateur) pour la réalisation des mêmes actions.

Le premier programme d'ordinateur 26 comporte des instructions pour l'exécution d'une génération du signal C par l'unité de traitement 22. Ce type de programme d'ordinateur est bien connu et ne sera pas détaillé.

Le second programme d'ordinateur 28 comporte des instructions pour l'exécution d'une analyse du signal Rs par l'unité de traitement 22. Plus précisément, ces instructions sont conçues pour détecter, par une analyse du signal Rs :

- un premier instant ti de réception d'une première réponse au signal excitateur E, cette première réponse résultant d'une première perturbation acoustique due à l'hétérogénéité optique de la cible T dans le milieu 12,

- un deuxième instant t ₂ de réception d'une deuxième réponse au signal excitateur E, cette deuxième réponse résultant d'une deuxième perturbation acoustique due à l'hétérogénéité acoustique de la cible T dans le milieu 12, et éventuellement

- un troisième instant t ₃ de réception d'une troisième réponse au signal excitateur E, cette troisième réponse résultant d'une troisième perturbation acoustique due à l'hétérogénéité optique de la source d'émission du signal excitateur E par rapport au milieu 12. La détection optionnelle du troisième instant t ₃ est par exemple réalisée conformément à l'enseignement du document WO 201 1/096198 A1 précité. Elle permet de déterminer, sur connaissance a priori de la distance notée SD entre la source d'émission du signal excitateur E et le capteur du signal acoustique R, quelle est la vitesse moyenne v de propagation d'une onde acoustique dans le milieu 12. En notant t ₀ l'instant d'émission du signal excitateur E, la vitesse v s'estime à l'aide de la relation suivante :

SD

v = .

f3 ^{— f} 0

En variante, la vitesse v peut être considérée comme connue a priori, de sorte que son estimation n'est pas nécessaire.

La détection du premier instant ti est réalisée conformément à l'enseignement du document WO 201 1/096198 A1 ou des autres documents de l'état de la technique cités précédemment. La propriété, selon laquelle la génération de la première perturbation acoustique au niveau de la cible se fait depuis la source à une vitesse voisine de celle de la lumière alors que la propagation de la perturbation acoustique depuis la cible vers le capteur acoustique se fait à la vitesse du son, est alors exploitée pour effectuer le calcul suivant :

SD

TD = v. (t ₁ - t ₀ ) , ou TD = - où TD est la distance entre la cible T et le capteur du signal acoustique R. Il est en effet considéré avec une très bonne approximation que la première perturbation acoustique est générée au niveau de la cible T à l'instant t ₀ alors que la réponse acoustique qu'elle engendre est reçue par le capteur du signal acoustique R à l'instant ti .

La détection du deuxième instant t ₂ est réalisée en contradiction avec l'enseignement des documents de l'état de la technique cités précédemment. Aucun de ces documents n'exploite la propriété selon laquelle une perturbation (i.e. la deuxième perturbation précitée) est générée au niveau de la cible en raison de son hétérogénéité acoustique dans le milieu 12. Pire, certains documents de l'état de la technique enseignent même de filtrer le signal acoustique R de manière à supprimer, entre autre, cette deuxième perturbation, pour supprimer l'effet des inhomogénéités acoustiques dans le milieu 12 perçues d'une façon générale comme un artéfact à éliminer. Au contraire et conformément à l'invention, la propriété, selon laquelle la génération de la deuxième perturbation acoustique au niveau de la cible se fait depuis la source à la vitesse du son et se propage ensuite de la cible jusqu'au capteur également à la vitesse du son, est exploitée pour effectuer le calcul suivant :

ST + TD = v. (t ₂ - t„) ,

où ST est la distance entre la source d'émission du signal excitateur E et la cible T. La deuxième perturbation acoustique est en effet générée au niveau de la source à l'instant t ₀ alors que la réponse acoustique qu'elle engendre est reçue par le capteur du signal acoustique R à l'instant t ₂ .

On déduit des deux calculs précédents une expression de la distance ST :

SD

ST = v. t ₂ - t ₁ ) , ou ST =——. t ₂ - t ₁ ) .

t3 ^t 0

Par conséquent, après avoir détecté les instants ti , t ₂ et éventuellement t ₃ , les instructions du second programme d'ordinateur 28 sont conçues pour :

- estimer la distance TD entre la cible et le capteur récepteur du signal acoustique R à l'aide du premier instant de réception ti , et éventuellement du troisième instant de réception t ₃ si la vitesse v n'est pas connue, selon le calcul détaillé précédemment,

- estimer la distance ST entre la source d'émission du signal excitateur E et la cible à l'aide du premier instant de réception ti , du deuxième instant de réception t ₂ , et éventuellement du troisième instant de réception t ₃ si la vitesse v n'est pas connue, selon le calcul détaillé précédemment,

- fournir une localisation de la cible à partir des distances ST et TD précédemment estimées.

La fourniture de la localisation se fait d'une façon générale dans l'espace tridimensionnel par intersection de deux sphères, la première sphère ayant pour centre le capteur récepteur du signal acoustique R et pour rayon la distance TD, la seconde sphère ayant pour centre la source d'émission du signal excitateur E et pour rayon la distance ST. Elle peut se faire plus précisément et avantageusement dans un plan d'illumination du milieu 12, l'intersection des deux sphères devenant ainsi une intersection de deux cercles. Ces deux cercles ayant a priori deux points d'intersection, il subsiste un doute sur la localisation de la cible T. Ce doute peut être levé à l'aide d'une autre émission/réception ou plus simplement à l'aide de considérations géométriques, par exemple si l'un des deux points se situe à l'extérieur du milieu observé 12. Ainsi, dans un plan, la localisation de la cible peut être résolue sans ambiguïté à partir d'une seule émission de signal excitateur E. En revanche, dans l'espace tridimensionnel, une émission supplémentaire peut être nécessaire pour résoudre la localisation sans ambiguïté par intersection de trois sphères. On notera qu'en outre les milieux à hétérogénéités acoustiques induisent une incohérence des positions estimées : l'incertitude sur une position globale peut cependant être minimisée en multipliant les angles d'examen du milieu observé 12.

Un procédé de localisation de la cible T dans le milieu 12 mis en œuvre par le dispositif 10 précédemment décrit va maintenant être détaillé en référence à la figure 2.

Au cours d'une première étape d'émission 100, le signal excitateur optique E est émis à l'instant t ₀ à partir d'une source S de l'émetteur 14 sur réception par ce dernier du signal de commande C émis par l'unité de traitement 22 du calculateur 20.

Au cours d'une étape suivante de réception 102, le signal acoustique R résultant de cette émission est reçu, à partir de l'instant t ₀ et au moins jusqu'aux instants ti , t ₂ et éventuellement t ₃ , par un capteur D du récepteur 16 puis prétraité par la carte d'acquisition 18 qui transmet le signal acoustique prétraité Rs à l'unité de traitement 22 du calculateur 20.

Au cours d'une étape optionnelle de détection 104 suivant l'étape 102, le troisième instant t ₃ est détecté par exécution du second programme d'ordinateur 28. Le principe de cette détection sera détaillé en référence à la figure 3.

Au cours d'une étape de détection 106 suivant l'étape 102, le premier instant ti est détecté par exécution du second programme d'ordinateur 28. Le principe de cette détection est le même que celui de l'étape optionnelle 104.

Au cours d'une étape de détection 108 suivant l'étape 102, le deuxième instant t ₂ est détecté par exécution du second programme d'ordinateur 28. Le principe de cette détection est le même que celui de l'étape optionnelle 104.

Au cours d'une étape d'estimation 1 10 suivant l'étape 106, et éventuellement l'étape optionnelle 104 le cas échéant, la distance TD entre la cible T et le capteur acoustique D est estimée par exécution du second programme d'ordinateur 28 conformément au calcul détaillé précédemment.

Au cours d'une étape d'estimation 1 12 suivant l'étape 108, et éventuellement l'étape optionnelle 104 le cas échéant, la longueur ST+TD est estimée par exécution du second programme d'ordinateur 28 conformément au calcul détaillé précédemment.

Au cours d'une étape d'estimation 1 14 suivant les étapes 1 10 et 1 12, la distance ST entre la source S et la cible T est estimée par exécution du second programme d'ordinateur 28 par soustraction de TD à la longueur ST+TD. En variante, au lieu de procéder au calcul en deux étapes 1 12 et 1 14, on notera que la distance ST peut être directement calculée à l'aide des instants ti et t ₂ (et éventuellement t ₃ ) conformément au calcul détaillé précédemment.

Enfin, au cours d'une dernière étape de localisation 1 16, une localisation géométrique de la cible T par intersection de sphères ou de cercles est obtenue par exécution du second programme d'ordinateur 28, à partir des distances estimées ST et TD. De façon optionnelle et en variante, d'autres algorithmes de reconstruction tomographique peuvent être utilisés tels que par exemple une méthode de sommation de rétroprojections filtrées ou toute méthode algébrique connue.

Dans l'exemple illustré sur la figure 3, le procédé de la figure 2 mis en œuvre par le dispositif 10 de la figure 1 est appliqué avec une seule source S d'émission de signal d'excitation E et un seul capteur acoustique D par souci de simplicité.

En partie supérieure gauche de la figure 3, les étapes 100, 102 et 104 sont appliquées au milieu 12 sans la cible T. Le signal Rs obtenu présente une réponse de forme impulsionnelle dans laquelle l'instant t ₃ est détecté de la façon suivante. Cet instant est choisi comme correspondant au passage à zéro de la réponse impulsionnelle, celle-ci présentant un pic supérieur suivi d'un pic inférieur. Dans l'exemple particulier et expérimental illustré sur la figure 3, l'instant t ₃ est détecté à 5,39.10 ^"5 s, l'origine des temps étant choisie à t ₀ . Pour une vitesse de son v = 1485 m/s, on en déduit SD ~ 8 cm. En pratique, c'est plutôt le calcul inverse qui est effectué. Connaissant par exemple a priori la distance SD de 8 cm entre la source S et le capteur D, on en déduit la valeur à 1485 m/s de la vitesse du son dans le milieu 12 grâce à la détection de t ₃ = 5,39.10 ^"5 s.

En partie supérieure droite de la figure 3, les étapes 100, 102 et 104 sont appliquées au milieu 12 avec la cible T. Le signal Rs obtenu présente toujours une réponse de forme impulsionnelle dans laquelle l'instant t ₃ peut être détecté de la même façon que précédemment. Mais, dans l'exemple illustré, les réponses impulsionnelles dues aux inhomogénéités optique et acoustique de la cible T ne sont pas visibles, celles-ci étant d'amplitudes très largement inférieures à celle due à l'inhomogénéité optique de la source.

Il est donc avantageux dans cet exemple de procéder à une soustraction entre les deux signaux Rs obtenus en partie supérieure gauche et droite de la figure 3. Concrètement, dans une application de détection et localisation d'une tumeur dans un tissu biologique, il suffit de disposer de signaux Rs de référence relatifs à un tissu sain qu'il convient de soustraire aux signaux captés pour un tissu biologique à tester. Dans certaines applications et selon les contrastes de propriétés optiques et acoustiques entre les milieux, les trois réponses impulsionnelles peuvent être visibles sur le même signal acoustique de sorte qu'une soustraction avec un signal de référence n'est pas toujours nécessaire. Mais c'est en revanche généralement nécessaire pour des applications d'observation de tissus biologiques épais.

Dans l'exemple illustré, après la soustraction précitée, on obtient le signal Rs illustré en partie inférieure de la figure 3. Deux réponses impulsionnelles correspondant respectivement aux instants ti et t ₂ sont maintenant bien visibles. La première, correspondant à l'instant ti , présente un pic inférieur suivi d'un pic supérieur entre lesquels le passage à zéro détermine l'instant ti . Dans l'exemple particulier et expérimental illustré sur la figure 3, l'instant ti est détecté à 3,39.10 ^~5 s, l'origine des temps étant choisie à t ₀ . Pour une vitesse de son v = 1485 m/s, on en déduit TD ~ 5 cm. La deuxième réponse impulsionnelle, correspondant à l'instant t ₂ , présente un pic supérieur suivi d'un pic inférieur entre lesquels le passage à zéro détermine l'instant t ₂ . Dans l'exemple particulier et expérimental illustré sur la figure 3, l'instant t ₂ est détecté à 6,07.10 ^~5 s, l'origine des temps étant choisie à t ₀ . Pour une vitesse de son v = 1485 m/s, on en déduit ST ~ 4 cm.

Enfin, par reconstruction géométrique (intersection de deux sphères ramenée à une intersection de deux cercles si l'on se place dans un plan), la localisation de la cible T est identifiée à une inconnue près éventuellement entre deux localisations T et T'. Comme indiqué précédemment, cette inconnue peut être résolue par exemple si T' se trouve géométriquement à l'extérieur du milieu 12.

Comme illustré sur la figure 4, pour des applications d'imagerie par reconstruction tomographique, l'ensemble émetteur/récepteur 14, 16 peut comporter plusieurs sources et plusieurs capteurs acoustiques. Le milieu 12 observé est par exemple en forme de disque et plusieurs sources S1 , S2, S3, S4, S5, S6, S7, S8 sont réparties régulièrement sur un cercle en périphérie de ce disque. De même, plusieurs capteurs acoustiques D1 , D2, D3, D4, D5, D6, D7, D8 sont répartis régulièrement sur un cercle en périphérie de ce disque, intercalés entre les sources. Grâce à cette configuration de l'ensemble émetteur/récepteur 14, 16, il peut être procédé à une reconstruction tomographique optoacoustique d'une image de la cible T dans le milieu 12 dans lequel elle est située, à l'aide des localisations obtenues de la cible T par émissions/réceptions successives réalisées à l'aide de ces multiples sources et capteurs. Un exemple de procédé d'imagerie par reconstruction tomographique intégrant un procédé de localisation selon l'invention est illustré sur la figure 5.

Au cours d'une étape 200, un ensemble de mesures par émissions/réceptions successives est réalisé à l'aide des sources optiques et des capteurs acoustiques de l'émetteur 14 et du récepteur 16.

Au cours d'une étape 202 suivant l'étape 200, les étapes 100 à 1 16 du procédé de la figure 2 sont exécutées pour chaque émission/réception de manière à fournir une localisation précise de la cible T dans le milieu 12.

Au cours d'une étape 204 suivant l'étape 202, la carte CPinit de distribution de pression initiale dans le milieu 12 est établie par simulation de l'illumination du milieu 12 par les sources optiques de l'émetteur 14, cette simulation intégrant, conformément à l'invention, la localisation de la cible T telle qu'estimée à l'étape 202.

A l'aide d'informations géométriques G sur les sources optiques de l'émetteur 14, sur les capteurs acoustiques du récepteur 16 et sur les contours extérieurs du milieu 12 illuminé, un modèle complet de ce milieu 12 incluant la cible T est établi lors d'une étape 206. De façon optionnelle et plus généralement, toute information a priori (morphologie, paramètres physiques ou physiologiques moyens, ...) peut être exploitée pour établir ce modèle.

Les mesures réalisées lors de l'étape 200 sont simulées sur le modèle établi à l'étape 206 au cours d'une étape suivante 208.

Les mesures simulées à l'étape 208 sont comparées aux mesures réelles de l'étape 200 au cours d'une étape 210. Cette comparaison se fait de façon connue en soi par minimisation d'une fonction d'erreur, par exemple une erreur quadratique.

Au cours d'une étape de test suivante 212, la fonction d'erreur est comparée à un seuil ε prédéterminé. Si elle reste supérieure à ce seuil, le procédé passe à une étape 214 de révision de la carte CPinit de distribution de pression initiale. Sinon, le procédé passe à une étape 216 d'enregistrement, et éventuellement d'affichage, d'une image représentant le milieu 12 avec la cible T correctement localisée, cette image étant définie sur la base du modèle établi à l'étape 206.

Au cours de l'étape 214, la révision de la carte CPinit de distribution de pression initiale se fait de façon connue en soi par optimisation de la fonction d'erreur. A l'issue de cette étape, le procédé reprend à l'étape 206 de manière à mettre à jour le modèle établi précédemment en fonction de la carte CPinit de distribution de pression initiale telle que révisée à l'étape 214. Il apparaît clairement qu'un dispositif et un procédé de localisation tels que ceux décrits précédemment permettent de localiser plus efficacement une cible dans un milieu illuminé par une technologie optoacoustique. Notamment, les principes de détection proposés permettent de s'affranchir d'une combinaison de technologies différentes pour parvenir à une localisation, puis à une représentation, précise de cibles dans un milieu. Ils exploitent de façon optimale les signaux acoustiques reçus en réponse aux excitations optiques émises, en tirant parti des propriétés d'inhomogénéité optiques mais également acoustiques des cibles à détecter.

Les résultats sont particulièrement convaincants dans le domaine de l'imagerie médicale, pour la détection de tumeurs dans des tissus biologiques. En outre, en imagerie médicale, l'application d'un procédé de localisation selon l'invention est compatible avec l'utilisation d'agents de contrastes comme dans les autres technologies conventionnelles.

On notera par ailleurs que l'invention n'est pas limitée aux modes de réalisation décrits précédemment.

En particulier, la description a porté plus précisément sur la mise en œuvre de l'invention dans le contexte d'une technologie d'acquisition optoacoustique. Mais il est bien évident que ses principes s'adaptent simplement et de façon connue en soi par l'homme du métier à une technologie d'acquisition thermoacoustique, notamment parce que l'hypothèse de la propagation bien plus rapide des ondes électromagnétiques dans le milieu comparée à celle du son est également vérifiée en technologie thermoacoustique.

En fait, la technologie thermoacoustique, comme l'optoacoustique, utilise des émetteurs d'ondes électromagnétiques et des récepteurs acoustiques. Ce qui distingue principalement ces deux technologies, ce sont les propriétés, de diffusion ou diffraction, exploitées du milieu observé en fonction des longueurs d'ondes émises. Lorsque les longueurs d'ondes émises engendrent un phénomène de diffusion dans le milieu, c'est l'équation de transfert radiatif (dite équation ETR) qui s'applique et qui est exploitée par la technologie optoacoustique. Lorsque les longueurs d'ondes émises engendrent un phénomène de diffraction dans le milieu, ce sont les équations de Maxwell qui s'appliquent et qui sont exploitées par la technologie thermoacoustique. Le fait d'engendrer un phénomène de diffusion ou de diffraction est une question de taille des constituants du milieu observé par rapport aux longueurs d'ondes émises. Ainsi par exemple, lorsque l'on souhaite observer un milieu constitué d'un tissu biologique, la technologie optoacoustique s'applique si l'émetteur émet des ondes dans le spectre de l'optique visible ou du proche infrarouge, car alors le milieu est diffusant, et la technologie thermoacoustique s'applique si l'émetteur émet des micro-ondes, car alors le milieu est diffractant. Lorsque l'on souhaite observer un milieu constitué d'un matériau manufacturé transparent, la technologie thermoacoustique s'applique si l'émetteur émet des ondes dans le spectre de l'optique visible, car alors le milieu est diffractant.

Plus précisément, la transposition de l'enseignement détaillé précédemment en technologie thermoacoustique se fait sur la base des considérations suivantes :

- le milieu à observer est diffractant et non pas diffusant,

- en ce qui concerne l'émission de signaux d'excitation, il convient de se placer dans la bonne gamme de longueurs d'ondes en fonction de la nature du milieu observé de manière à engendrer le phénomène de diffraction exploitable en technologie thermoacoustique, par exemple la gamme des micro-ondes pour l'examen de tissus biologiques et la gamme des ondes optiques visibles pour l'examen de matériaux transparents,

- les contrastes électromagnétiques sont à variation locale des propriétés diélectriques, par variation des parties réelle et imaginaire des constantes de permittivité et conductivité diélectrique du milieu observé, de sorte que seules les équations de propagation changent,

- les équations régissant la carte de distribution initiale de pression sont les équations de Maxwell et pas l'équation de transfert radiatif ETR. Il apparaîtra plus généralement à l'homme de l'art que diverses modifications peuvent être apportées aux modes de réalisation décrits ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. Dans les revendications qui suivent, les termes utilisés ne doivent pas être interprétés comme limitant les revendications aux modes de réalisation exposés dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents que les revendications visent à couvrir du fait de leur formulation et dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en œuvre de l'enseignement qui vient de lui être divulgué.