FR 2938181 A1 20100514 FR 0857699 A 20081113 1 PROCEDE ET DISPOSITIF DE LOCALISATION DE FLUOROPHORES OU D'ABSORBEURS DANS UN MILIEU ENVIRONNANT DESCRIPTION DOMAINE TECHNIQUE ET ART ANTÉRIEUR L'invention concerne le domaine de l'imagerie optique diffuse sur les tissus biologiques par des méthodes résolues en temps. Elle s'applique à la localisation de fluorophores ou d'absorbeurs dans un milieu diffusant tel que, par exemple, un organe d'un animal ou d'un être humain (cerveau, sein, tout organe où des fluorophores peuvent être injectés). En particulier, un intérêt d'une molécule qui absorbe est le suivant. Certaines tumeurs cancéreuses sont visibles par le simple fait qu'elles présentent un coefficient d'atténuation plus important que les tissus sains. Dans ce cas, il est important de pouvoir identifier ces tumeurs. Dans ce contexte, on peut également être amené à utiliser un fluorophore, qui est un marqueur fluorescent spécifique. Celui-ci vient se fixer de façon préférentielle sur les cellules cibles d'intérêt (par exemple des cellules cancéreuses). Un tel fluorophore offre un meilleur contraste de détection qu'un marqueur non spécifique. Les techniques optiques d'imagerie moléculaire de fluorescence ont pour but de localiser spatialement ces marqueurs fluorescents. Les systèmes de tomographie optique utilisent diverses sources de lumière. Il existe donc 2 des appareils qui fonctionnent en mode continu, d'autres en mode fréquentiel (qui utilisent des lasers modulés en fréquence) et enfin des appareils fonctionnant en mode temporel, qui utilisent des lasers pulsés. Les données temporelles sont celles qui contiennent le plus de contenu informationnel sur le tissu traversé, mais pour lesquelles les techniques de reconstruction sont les plus complexes. La mesure en chaque point d'acquisition est en effet une fonction dépendante du temps (appelée TPSF pour � Temporal Point Spread Function �, ou fonction d'étalement temporel). On cherche à extraire des paramètres simples de la TPSF, dont on connaît par ailleurs l'expression théorique. Ensuite, la résolution du problème inverse permet de retrouver la distribution des marqueurs fluorescents. Dans certains cas, on cherche à localiser un fluorophore en vue d'une intervention ultérieure. A l'heure actuelle, la localisation de cellules malades peut se faire à l'aide des techniques de biopsie. Mais, afin de trouver la position, même approximative, de ces cellules, on doit effectuer plusieurs prélèvements. Cette technique invasive est évidemment délicate à mettre en oeuvre et est consommatrice de temps. Il se pose donc le problème de la détermination de la localisation, même de manière approximative, d'un fluorophore ou d'un absorbeur fixé sur une zone d'un milieu, par exemple un milieu biologique, cela avec un nombre limité de mesures, et une durée de dépouillement relativement faible, voire 3 proche du temps réel. L'invention trouve toute son utilité pour réaliser une telle détermination dans des milieux présentant une accessibilité limitée, par exemple lors du dépistage de tumeurs cancéreuses de la prostate à partir de mesures endoscopiques, ou lors du dépistage de cancers du sein, ou d'ostéosarchomes.... EXPOSÉ DE L'INVENTION L'invention concerne d'abord un procédé de localisation d'au moins un fluorophore ou d'au moins un élément (molécule notamment) absorbant - encore appelé absorbeur - dans un milieu diffusant, à l'aide d'une source de rayonnement apte à émettre un rayonnement d'excitation de ce fluorophore et des moyens de détection apte à mesurer : - un signal de fluorescence (a)f1uo) émis par ce fluorophore, - ou un signal réémis à la même longueur d'onde que celle d'excitation (issue de la source de rayonnement) par l'élément absorbant qui absorbe et diffuse le rayonnement absorbé. L'élément absorbant peut faire partie d'une zone d'un milieu qui présente un coefficient d'absorption supérieur à celui du milieu diffusant. On a déjà donné l'exemple de certaines tumeurs cancéreuses. Ce procédé de localisation comporte les étapes suivantes: a) pour au moins 3 couples différents de positions de la source de rayonnement et des moyens de détection, au moins une excitation par un rayonnement provenant de la source de rayonnement, et au moins une 4 détection par les moyens de détection du signal de fluorescence émis par ce fluorophore après cette excitation, ou du signal d'émission émis par l'absorbeur ou par l'élément absorbant, b) pour chacun de ces couples, l'identification d'une surface sur laquelle le fluorophore (ou l'élément absorbant) est situé, ou d'un volume comprenant cette surface et dans lequel le fluorophore ou l'absorbeur est situé, c) une estimation de la localisation du fluorophore ou de la localisation de l'élément absorbant dans son milieu diffusant, par calcul de l'intersection des trois surfaces, ou éventuellement dans un volume autour de cette intersection. Un tel procédé peut être mis en oeuvre quelles que soient les positions respectives du fluorophore ou de l'absorbeur, de la source et du détecteur. La précision du procédé ne sera pas la même suivant que l'on prend des positions de la source et des moyens de détection rapprochées ou éloignées du fluorophore ou de l'absorbeur. Mais le procédé permet dans tous les cas une localisation, même grossière, avec un faible nombre de mesures, avec un dépouillement rapide, voire en temps réel. Un tel procédé permet donc d'orienter la réalisation de biopsies ou d'examens non invasifs complémentaires. En outre ce procédé est de mise en oeuvre très rapide, on peut réaliser une localisation d'un fluorophore ou d'absorbeur en moins d'une minute, ou même en moins de 15 s. Pour chaque couple des au moins 3 couples différents, on peut définir un histogramme représentant le nombre de photons détectés en fonction du temps ; cet histogramme peut être appelé fonction temporelle de 5 fluorescence. La source de rayonnement ou de lumière et les moyens de détection peuvent être respectivement l'extrémité d'une fibre optique qui amène la lumière d'excitation dans le milieu et l'extrémité d'une fibre optique qui prélève une partie de la lumière d'émission. Pour la source, il peut également s'agir d'une source lumineuse de type Laser ou LED. En ce qui concerne le détecteur, il peut s'agir d'un tube photomultiplicateur, ou d'un capteur d'image Selon l'invention l'élément visé par la détection peut être un fluorophore unique, ou une pluralité ou une agrégation de fluorophores groupés sensiblement en un même lieu. Il peut également s'agir d'un absorbeur ou d'une pluralité d'absorbeurs groupés sensiblement en un même lieu (agrégation). Dans la suite du texte, un absorbeur sera assimilé à un fluorophore dont la longueur d'onde d'émission est égale à la longueur d'onde d'excitation, et dont la durée de vie est nulle. Autrement dit un milieu absorbant, ou un absorbeur, est un fluorophore dont la longueur d'onde d'émission est égale à la longueur d'onde d'excitation, et dont la constante de temps 'L est nulle. On n'emploiera donc que l'expression � fluorophore �, sauf lorsque des spécificités des absorbeurs sont expliquées. 6 Le volume autour de l'intersection des trois surface calculées peut résulter de l'intersection de trois volumes, chaque volume étant associé à l'une des surfaces et comportant l'ensemble des points de la surface et l'ensemble des points qui satisfont à l'équation de la surface avec une marge d'erreur, par exemple + 10% ou + 5%. Le milieu diffusant environnant le fluorophore peut avoir toute forme de géométrie : par exemple, il est de type infini. Il peut aussi être de type semi-infini, limité par une paroi. Un milieu semi infini, pour lequel chacune des extrémités de fibres d'amenée du signal d'excitation et de prélèvement du signal de fluorescence est disposée à plus de 1 cm, ou de 1,5 cm ou de 2 cm de la paroi qui limite ce milieu peut être traité, du point de vue du procédé selon l'invention, comme un milieu infini. Le milieu diffusant peut encore former une tranche (géométrie dite � slab �), limitée par deux surface parallèles. L'invention s'applique aussi à un milieu de forme quelconque, dont la surface extérieure est décomposée en une série de plans. Dans certains modes de réalisation, les surfaces à l'intersection desquelles se trouve le fluorophore, peuvent avoir la forme d'ellipsoïdes, mais il peut également s'agir de surfaces 3D n'ayant pas cette forme. Selon un mode de réalisation, l'étape b) comporte pour chaque couple de position, un calcul d'un moment temporel normé, d'ordre n (n> 1), d'une fonction temporelle de fluorescence. Un tel moment peut être obtenu par une dérivée d'ordre n d'une transformée de ladite fonction temporelle de fluorescence, par exemple une transformée de Fourier, ou une transformée de Laplace ou une transformée de Mellin. On rappelle qu'un moment normé d'ordre n (n> 1), noté mn d'une distribution g(t) est défini par ftng(t)dt m = n +n, g(t)dt Dans certains cas, il peut être avantageux de réaliser une excitation par une source de lumière, et une détection par les moyens de détection, pour 4 couples différents de positions de la source de rayonnement et des moyens de détection. On peut ensuite : - sélectionner parmi les quatre couples, le couple de points, dit quatrième couple de points, pour lequel la durée moyenne d'arrivée des photons, de l'émetteur au récepteur, est la plus courte, et calculer, pour ce quatrième couple de points, l'équation d'une surface sur laquelle le fluorophore est situé ; - et, pour chacun des trois autres couples de points, après calcul d'une première équation d'une surface sur laquelle le fluorophore est situé, on peut effectuer une différence entre cette équation et celle associée au quatrième couple, pour obtenir une deuxième équation indépendante de la durée de vie du fluorophore. 8 On peut donc, selon un mode de réalisation, lors de l'étape b) pour chaque couple de position, calculer une équation de chaque surface, indépendante de la durée de vie du fluorophore, qui peut résulter de la différence entre le temps moyen calculé pour chaque couple de position et le temps moyen calculé pour un autre couple de position. Cet autre couple de position peut être celui ayant un temps moyen minimal ou ayant un temps moyen calculé minimal. Un procédé selon l'invention peut en outre comporter une représentation visuelle ou graphique de la position ou de la répartition du ou des fluorophores. L'invention concerne également un procédé pour localiser une zone dans un milieu, comportant l'introduction dans ce milieu d'au moins un fluorophore ou d'un absorbeur, et la localisation de ce fluorophore ou de cet absorbeur dans ledit milieu par un procédé selon l'invention. La zone peut être par exemple un milieu constitué d'un organe humain ou animal, le fluorophore ou l'absorbeur permettant de marquer cette zone. Le but peut être de localiser des cellules malades dans un tel milieu. Comme déjà indiqué ci- dessus, l'invention permet de réaliser cette localisation en un temps rapide, ou même en temps réel. Un procédé selon l'invention peut être mis en oeuvre de manière invasive ou non invasive. L'invention concerne également un dispositif pour mettre en oeuvre un procédé selon l'invention, et permettant la localisation d'un fluorophore dans un milieu diffusant. 9 L'invention concerne donc également un dispositif pour la localisation d'un fluorophore dans un milieu diffusant, qui peut être de l'un des types évoqués plus haut, comportant : - au moins une source de lumière ou de rayonnement, apte à émettre un rayonnement d'excitation de ce fluorophore, - au moins un moyen de détection apte à mesurer un signal de fluorescence ((Dfiuo ) émis par ce fluorophore, - des moyens de calcul, pour calculer : a) pour chaque couple de points parmi au moins 3 couples différents de positions d'une source de rayonnement et d'un moyen de détection, la détermination, après une excitation réalisée par un rayonnement provenant de la source de rayonnement, et une détection par les moyens de détection des signaux de fluorescence émis par ce fluorophore après cette excitation, d'une surface sur laquelle le fluorophore est situé, ou un volume comprenant cette surface et dans lequel le fluorophore est situé, b) une estimation de la localisation du fluorophore dans son milieu environnant, par calcul de l'intersection des trois surfaces ou d'un volume autour de cette intersection. Les moyens de calcul permettent de calculer, pour chaque couple de position, un moment temporel normé, d'ordre n (n> 1), d'une fonction temporelle de fluorescence. Ce moment peut être obtenu à partir d'une fonction fréquentielle, déduite de 10 ladite fonction temporelle de fluorescence par exemple par transformée de Fourier Selon un mode avantageux, les moyens de calcul permettent, pour chaque couple de position, un calcul d'une équation de chaque surface, indépendante de la durée de vie du fluorophore. A cette fin, les moyens de calcul permettent de réaliser la différence entre le temps moyen calculé pour chaque couple de position et le temps moyen calculé pour un autre couple de position. BRÈVE DESCRIPTION DES DESSINS - Les figures 1A et 1B représentent chacune un exemple de dispositif expérimental pour la mise en oeuvre de l'invention, - les figures 2A et 2B représentent respectivement une série d'impulsions laser et de photons uniques émis, et une courbe de fluorescence obtenue à partir des données relatives aux photons uniques, - la figure 3 représente une configuration théorique, en deux dimensions, - les figures 4A-4F représentent le dispositif expérimental d'une mesure réalisée, le positionnement des points source et de détection, des courbes de fluorescence mesurées, et la position calculée et la position réelle d'un fluorophore; - la figure 5A représente des couples de positions d'une source et d'un détecteur, dans un plan XY, en vue d'une mesure selon l'invention, en configuration semi-infinie, 11 - les figures 5B-5D représentent des portions d'ellipsoïdes utilisés dans la cadre d'une mesure selon l'invention, en configuration semiinfinie, - les figures 6A-6C représentent chacune une localisation d'un fluorophore obtenue avec un procédé selon l'invention, en configuration semiinfinie, - les figures 7A - 7D représentent des surfaces mises en oeuvre dans le cadre d'un autre mode de réalisation de l'invention, - les figures 8A-8D représentent des portions de surfaces utilisées dans la cadre d'une mesure selon l'invention, en configuration de tranche. - la figure 9 représente un schéma de déroulement d'un procédé selon l'invention. les figures 10A-10D représentent diverses dispositions, par rapport à un milieu étudié, de moyens d'excitation et de détection, ou de fibres optiques amenant à ce milieu un rayonnement d'excitation et transmettant un rayonnement à détecter. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS La figure 1A est un exemple de système expérimental 2 pour mettre en oeuvre un procédé selon l'invention. Des représentations schématiques d'un milieu et de moyens d'excitation et de détection sont donnés plus loin en liaison avec les figures 10A-10D. Le dispositif de la figure 1A utilise comme détecteur 4 un photomultiplicateur et une carte TCSPC (rime Correlated Single Photon Counting), en fait 12 intégrée dans un ensemble de moyens 24 de traitement de données. La lumière est émise par une source 8 de rayonnement en impulsions, est envoyée et collectée par des fibres 10, 12 qui peuvent être déplacées. Les deux fibres peuvent être montées sur des moyens de déplacement en translation verticaux et horizontaux (suivant les axes X, Y et Z de la figure 1A). Pour cela on peut mettre en oeuvre une ou des platines montées en translation pour déplacer les extrémités des fibres optiques. On peut aussi utiliser un laser, lui aussi monté sur une platine montée en translation. Un détecteur peut lui aussi être monté sur une telle platine. Dans tous les cas la translation peut être pilotée par ordinateur. La source 8 d'impulsions de rayonnement peut également être utilisée comme moyen de déclenchement de la carte TCSPC (voir la liaison 9 entre la source 8 et les moyens 24). Il est également possible de travailler avec des impulsions dans le domaine de la femto seconde, à condition de disposer de la source de rayonnement adéquate, c'est-à-dire d'une source laser 8 dont chaque impulsion a une largeur temporelle également dans le domaine de la femto seconde. Selon un mode particulier de réalisation, la source 8 est une diode laser pulsée, par exemple à une longueur d'onde voisine de 630 nm et avec un taux de répétition d'environ 50 MHz. 13 La lumière laser passe de préférence par un filtre interférentiel 14 pour éliminer d'éventuels modes secondaires. La fibre 12 collecte la lumière en provenance du milieu 20 étudié, en particulier la lumière diffusée par ce milieu. Un filtre interférentiel 16 et/ou un filtre coloré absorbant les faibles longueurs d'onde peuvent être placés devant le détecteur 4 pour sélectionner la lumière de fluorescence (par exemple : A > 650 nm, la source étant à une longueur d'onde de, par exemple, 631 nm) d'un fluorophore 22 disposé dans le milieu 20 et optimiser l'élimination de la lumière d'excitation. Dans le texte de la présente demande, quel que soit le mode de réalisation de l'invention, il sera notamment question de la position de la source de rayonnement, et/ou de la position des moyens de détection. Lorsque des fibres sont utilisées, ces positions s'entendent le plus souvent comme celles des extrémités des fibres qui amènent le rayonnement dans le milieu diffusant 20, ou à l'interface de ce milieu, et /ou celles qui collectent le rayonnement diffusé, ces dernières étant placées dans le milieu ou à l'interface de ce milieu, mais ne s'entendent alors pas comme étant celles de la source 8 proprement dite ou du détecteur 4 proprement dit. En variante, on peut utiliser une source de lumière délivrant un faisceau lumineux, par exemple un laser en émission, et une caméra munie d'un objectif en sortie, ce qui permet un couplage optique entre la face sensible de la caméra et la face de sortie de l'objet 14 observé (ce qui est le cas du dispositif de la figure 1B). En application de type � endoscopie �, on utilise préférentiellement un système à fibres, les fibres pouvant être par exemple intégrées dans une sonde endo-rectale, au moins une fibre transmettant la source de lumière pulsée d'excitation à la zone à examiner. Cette source de lumière pulsée peut être externe à la sonde. Au moins une autre fibre transmet le signal optique d'émission de la zone à examiner vers un photodétecteur, qui peut lui aussi être externe à la sonde. Selon la technique TCSPC (pour � Time Correlated Single Photon Counting �, ou comptage de photons uniques correlés en temps) on détecte, à l'aide du photomultiplicateur, un photon émis par le fluorophore après une impulsion de la source de rayonnement. Le système permet donc une détection résolue en temps des impulsions de fluorescence. Il permet de récupérer des photons de fluorescence. La figure 2A représente une série d'impulsions laser Li (i = 1 - 4) et une série de photons uniques pi (i = 1 - 4) correspondants, ces derniers étant détectés par un système de type TCSPC (Time Correlated Single Photon Counting). Chaque photon est en fait détecté par rapport au départ de l'impulsion correspondante : sur la figure 2A, ti représente la durée écoulée entre chaque impulsion laser Li et l'instant de détection de chaque photon pi. Il est donc ensuite possible d'établir une répartition statistique ou histogramme du temps d'arrivée des photons, comme illustré sur la figure 2B, 15 qui représente le nombre de photons de fluorescence détectés, en fonction du temps écoulé t par rapport à chaque impulsion laser. On peut déterminer, à partir d'un tel histogramme, des paramètres statistiques, notamment une moyenne du temps d'arrivée, ou temps moyen (il s'agit en fait de la moyenne des abscisses pondérée par les ordonnées de l'histogramme). Une telle courbe (Df1uo (t) qui, on le voit (voir également l'exemple de la figure 4C), permet d'utiliser toute l'information sur une large fenêtre temporelle, de part et d'autre du point d'intensité maximum (et pas seulement dans la partie montante du signal) peut donc ensuite être traitée pour en tirer des informations caractéristiques, tel que le temps d'arrivée, ou temps moyen comme on va l'expliquer ci-dessous. Des moyens électroniques 24 tels qu'un micro ordinateur ou un calculateur sont programmés pour mémoriser et traiter les données de la carte TCSPC. Plus précisément une unité centrale 26 est programmée pour mettre en oeuvre un procédé de traitement selon l'invention. Des moyens 27 d'affichage ou de visualisation permettent, après traitement, de représenter le positionnement du ou des fluorophores dans le milieu 20 examiné. D'autres techniques de détection peuvent être employées, par exemple avec une caméra intensifiée, et par exemple une caméra intensifiée ultra-rapide de type � gated camera � ; dans ce cas la caméra s'ouvre sur une porte temporelle, de largeur par 16 exemple environ 200 ps, puis cette porte est décalée, par exemple de 25 ps en 25 ps. La figure 1B est un exemple d'un autre système expérimental 2 utilisant comme détecteur 32 une caméra rapide : la position relative de la source (ou le détecteur) et l'objet peut aisément être réalisée. Un faisceau 30 d'excitation de la fluorescence d'un milieu 20, contenant un ou plusieurs fluorophores 22, est émis par une source 37 de rayonnement (non représentée sur la figure) qui peut être du même type que celle présentée ci-dessus en liaison avec la figure 1A. Un photodétecteur 36 permet de commander des moyens 40 formant une ligne à retard. La référence 24 désigne, comme sur la figure 1A, des moyens électroniques de traitement de données de type micro-ordinateur ou calculateur, programmés pour mémoriser et traiter les données de la caméra 32. Une unité centrale de ces moyens 24 est programmée pour mettre en oeuvre un procédé de traitement selon l'invention. Là encore, des moyens d'affichage ou de visualisation permettent, après traitement, de représenter le positionnement ou la localisation d'un ou de plusieurs fluorophores dans le milieu 20 examiné. La distance entre le faisceau lumineux émis par la source, et le détecteur, ou la partie du milieu optiquement couplée au détecteur, est par exemple comprise entre quelques mm et quelques cm, par exemple entre 1 mm et 20 cm dans le cas de milieux biologiques. C'est le cas notamment dans le cas d'une intégration dans une sonde endo-rectale. Le positionnement relatif de la source et du détecteur par rapport à un organe à 17 examiner peut être déterminé par des moyens de type optique ou ultrasonore, par exemple une sonde d'echographie. En variante des dispositifs ci dessus, on peut placer un nombre n de sources lumineuses et un nombre p de détecteurs, avec par exemple 1 AX avec une durée de vie T. Cette durée de vie correspond à la durée moyenne pendant laquelle les électrons excités restent dans cet état avant de retourner à leur état initial. Les développements qui vont suivre montrent que le cas d'un absorbeur, et donc d'une émission à longueur d'onde sensiblement égale à celle d'excitation avec une durée de vie nulle, peut être déduit de cas d'un fluorophore. L'équation de propagation de la lumière d'excitation dans un milieu diffusant produite par la source pulsée de rayonnement s'écrit : 1 Y'x (r't) +v.(ùDx (r)vOx (r,t))+pax(r)Y'x (r,t) _ Sx (rs ,t) où - rs désigne la position de la source, r désigne une position quelconque dans le milieu, la fonction a)X (r, t) désigne le débit de fluence du rayonnement d'excitation au point r l'instant t. le coefficient DX(r) est le coefficient de diffusion au point r,à la longueur d'onde d'excitation Xx SX(rs, t) désigne le flux de rayonnement émis, à l'instant t, par la source au point rs. (1) cn at30 19 Celle de la lumière d'émission s'écrit : 1 aom(r,t) +o.(-Dm(r)VOm(r,t))+�a. (r�m(r,t) _ g�ax-->m fo t exp(-(t-t'))~x(r,t')dt' n at r r (2) avec : la fonction cm(r, t) qui désigne le débit de fluence du rayonnement d'émission au point r l'instant t, le coefficient Dm(r), qui est le coefficient de diffusion au point r à la longueur d'onde d'émission Xm. Les conditions aux limites étant les suivantes: (r,t)+2AD,(r)~~x ,t=0 (3) 0,n(r,t)+2ADä,(r) an (r,t)=0 (4) et ce pour tout point r de la frontière 6S2 du milieu diffusant. 20 Ces équations sont aussi développées été expliquées dans la thèse de d'Aurélie Laidevant, � Méthodes optiques résolues en temps pour la tomographie de fluorescence dans les milieux 25 diffusants �, Université Joseph Fourier, Grenoble, 5 octobre 2006, p.51 et suivantes. Dans ces équations . - Pax et Jam sont les coefficients d'absorption du milieu diffusant à la longueur d'onde 30 d'excitation XXet d'émission Xm respectivement,15 20 ù faxm est le coefficient d'absorption associé à la concentration locale de fluorophores à la longueur d'onde d'excitation Xx. - n est le rendement quantique du fluorophore, - a)(r, t) est le débit de fluence associé à la densité de photons, en photons/m2s, ou, de manière équivalente, en Watt/m2, - A est égal à (1-Reff) / (1+Reff) , où Reff est le coefficient de réflexion dû à la différence d'indice entre l'air et le milieu 22 environnant le fluorophore. La condition aux limites exprime que, à la frontière du milieu, le débit de fluence diffusé est nul. La fonction théorique a)m peut-être considérée comme proportionnelle à la TPSF, cette dernière étant mesurée expérimentalement (figure 2B). Ainsi, les moments temporels normés de la TPSF sont égaux aux moments temporels normés de la fonction Ï. La résolution du système d'équations couplées obtenu n'est pas aisée, et l'exploitation du volume de données que contient le signal temporel est difficile. On s'intéresse en général à l'exploitation de tel ou tel intervalle de temps extrait de ces fonctions dépendantes du temps (appelées TPSF pour � Temporal Spread Function �, ou fonction d'étalement temporel)- pour réaliser une étape d'ajustement (`fitting') avec les données expérimentales. Pour exploiter l'information temporelle par le biais d'autres grandeurs, extraites des TPSF mesurées, on peut définir d'autres mesures telles que 21 la mesure des moments - transformées de Mellin de la TPSF, ou des Transformées de Fourier (et donc travailler en mode fréquentiel), ou des Transformée de Laplace ou par paramétrisation de la TPSF (c'est-à-dire la détermination de paramètres caractéristiques tels que les coordonnées du maximum, et/ou le coefficient de pente montante et/ou de pente descendante de la TPSF. L'intérêt de telles grandeurs est de réduire le volume de données en supprimant les éventuelles redondances existant entre elles Dans la suite de l'invention, on s'intéresse principalement au calcul de moments d'ordre n de l'histogramme, un tel calcul pouvant être réalisé par transformée de Fourier. En particulier, on peut, partant de la fonction cm, en calculer tout moment temporel normé, et notamment le moment d'ordre 1 normé. Une telle grandeur correspond au temps moyen d'arrivée des photons. Les moments temporels normés de la fonction cm, sont définis comme : +c kOm(r,t)dt JOm(r,t)dt (6) Le moment d'ordre 0 est: Ç1 (r, t)dt 22 On va maintenant traiter le cas d'un milieu infini. On peut montrer que le moment temporel normé d'ordre 1 de la fonction cm, définissant le temps moyen d'arrivée des photons w, en milieu infini s'écrit : tom(r,t)dt/ ~ om(r,t)dt= oo= x + m + déclin fluo -w (7) + +z 2cn .jpax Dx 2cn .j,u. Dm Ou X, m et déclin fluo désignent respectivement la durée moyenne du trajet sourcefluorophore, la durée du trajet fluorophore -détecteur, et la durée de vie moyenne de l'espèce excitée (cette dernière étant égale à 0 en présence d'un absorbeur).On rappelle que, en outre : rsr = rs - r rrd = r-rd m1 est le résultat de la mesure On voit qu'on a là l'équation d'une surface 3D, ayant ici la forme d'un ellipsoïde. Plus généralement, le moment normé d'ordre 20 k peut s'exprimer sous la forme : k~(w) mk = (l)k k x aw w=0 (w)w0 Où (')(w) est la transformée de Fourier de a)m: ~(w) _ f O(t)e-zwcdt a~(~) _ - tO(t)e-'tdt 2~c _~ aw 27c rsY (8) 25 23 Par exemple, on peut montrer que la variance de la distribution temporelle s, écritV = )2 >=m2 -mi , soit en milieu infini V= rsr + rrd +r2 2c2 pw3 cDx ZCn /Ja3mDm On observe qu'on obtient donc, également dans le cas de la variance, la définition d'une surface tridimensionnelle, ici, celle d'un ellipsoïde. Cette grandeur (la variance) peut être calculée dans le cadre de la présente invention à partir des moments d'ordre 2 et 1. Plus généralement encore on peut, selon l'invention, calculer tout moment d'ordre supérieur à 1, et établir une relation liant au moins un de ces moments aux distances respectives source-fluorophore et fluorophore-détecteur. Si l'on considère des moments d'ordre supérieur, l'équation reliant ces derniers, ou une combinaison de ces derniers, aux distances respectives source-fluorophore et fluorophore- détecteur, est également celle d'une surface tridimensionnelle. Préférentiellement, on utilise le moment d'ordre 1 normé. D'une manière générale, on peut voir un fluorophore comme capable d'absorber une certaine quantité de lumière (dans son état excité) et d'en réémettre une partie à une longueur d'onde supérieure (lors de sa désexcitation), avec un retard correspondant à la durée dans l'état excité. A partir des expressions données ci dessus on peut donc passer du cas du fluorophore à l'absorbeur simple en écrivant /1ax = ,uam , Dx =Dm et z = 0 . Ceci vaut pour tous les cas exposés dans la présente demande. Les considérations ci-dessus valent donc également pour le cas d'un objet absorbant. Le milieu infini est typiquement celui que l'on rencontre lorsque des mesures invasives sont réalisées, avec des fibres 10, 12 pour amener la lumière dans un milieu 21 (cas de la figure 10A) et pour prélever le signal de fluorescence émis par le fluorophore F, les extrémités des fibres étant positionnées dans le milieu diffusant, chacune à une distance d10 et d12 d'au moins 1 cm ou 1,5 cm d'une paroi 21 qui délimite le milieu. On peut plus généralement considérer qu'on se trouve dans un milieu infini lorsque la source et le détecteur sont placés à une profondeur suffisante dans le milieu diffusant. En pratique, cela correspond à des mesures invasives, réalisées en disposant par exemple des fibres optiques, ou une partie de ces fibres, à l'intérieur d'un milieu diffusant. On considère une première hypothèse, celle d'un milieu où les propriétés optiques sont identiques aux 2 longueurs d'onde d'excitation et de mesure, "x et à ce qui est par exemple sensiblement le cas lorsque ces deux longueurs d'onde sont distantes de par exemple de 10 nm ou de quelques dizaines de nm. Rappelons aussi que = Xm lorsqu'on n'a qu'un absorbeur non fluorescent. ; Les propriétés optiques du milieu aux deux longueurs d'ondes peuvent donc être 25 considérées comme identiques. On a alors une simplification de l'équation (7) ci-dessus : Ceci constitue la formule d'un ellipsoïde de foyers S et D, soit le lieu des points r tels que la distance source-fluorophore rsr plus la distance fluorophore-détecteur rrd est égale au produit de la vitesse apparente dans le milieu diffusant 2cn ~aD multiplié par le temps moyen d'arrivée des photons corrigé du retard introduit par le temps de vie T. Ainsi pour trouver la position 3D du fluorophore, il suffit d'effectuer 3 mesures à trois positions de sources et 3 positions de détecteurs. Le fluorophore est donc à l'intersection de 3 ellipsoïdes définis dans l'espace par 3 équations différentes. De préférence, on sélectionne des couples de positions définissant des directions sensiblement perpendiculaires les unes aux autres. Pour chaque couple, on procède à l'établissement d'une TPSF ou d'un histogramme comme expliqué ci-dessus. Dans la deuxième hypothèse, le cas où les propriétés optiques ne peuvent pas être considérées identiques aux 2 longueurs d'onde, l'équation (7) devient: (9) (10) Ceci est toujours l'équation d'un ellipsoïde mais il traduit maintenant le lieu des points r tels que la distance source-fluorophore rsr soit � pondérée � par la vitesse apparente vX (d'où le terme rsr/vx), et la distance fluorophore-détecteur rrd soit � pondérée � par la vitesse apparente vm (d'où le terme rrd/vm) Là encore, pour trouver la position 3D du fluorophore, on effectue 3 mesures à trois positions de sources et 3 positions de détecteurs. Le fluorophore est donc à l'intersection de 3 ellipsoïdes définis dans l'espace par 3 équations différentes. De préférence, on sélectionne des couples de positions définissant des directions sensiblement perpendiculaires les unes aux autres. Pour chaque couple de positions (source, détecteur), on procède à l'établissement d'une TPSF ou d'un histogramme comme expliqué ci-dessus. Egalement de préférence, lors de chaque acquisition, on tentera d'augmenter le plus possible la distance source- détecteur. Les expressions des temps moyens peuvent être données analytiquement pour des géométries autres que celles du milieu infini, plus complexes : milieu semi-infini, avec une interface plane ; milieu plan à faces parallèles... De manière générale, si l'on connait l'expression, analytique ou numérique, de la transformée de Fourier ii)(w) de la fonction alors : mi = =a(D(w)'w=0x 1 , aw d~(w)~w = 0 27 où (D(w) = 1 $O(t)e-zwtdt , &D(w) _ - i $tO(t)e-zwtdt 27r _oo aw 2~c La surface décrite n'est plus un ellipsoïde, mais reste une surface 3D ; Donc, là encore, on localise un fluorophore par recherche des intersections entre 3 surfaces, chacune obtenue pour un couple de point (position source, position mesure) différent. D'une manière générale, dans tous les cas concernés par l'invention (ceux-ci-dessus et les autres exposés ci-dessous), pour trouver la position 3D du fluorophore, on effectue 3 mesures avec trois positions, de préférence différentes entre elles, de sources et 3 positions, de préférence différentes entre elles, de détecteurs (autrement dit, une mesure pour chacun des 3 couples (rsi, rdi), i= 1,2,3, avec rsi # rs . pour i j et rdi rd , pour i j). Le fluorophore est donc à l'intersection de 3 ellipsoïdes définis dans l'espace par 3 équations différentes. De préférence, on sélectionne 3 couples de positions définissant des directions sensiblement perpendiculaires les unes aux autres. Dans certains cas, dont on verra un exemple plus loin, il peut être utile d'effectuer 4 mesures pour 4 couples de positions (source, détecteur). Pour chaque couple de positions (source, détecteur), on procède à l'établissement d'une TPSF ou d'un histogramme comme expliqué ci-dessus. Un autre exemple est le cas d'un milieu semi-infini. C'est le cas de tout organe sondé en surface, par exemple prostate, sein, cerveau, dont 28 l'épaisseur est suffisamment grande devant les distances caractéristiques source-détecteur. Ce cas est celui que l'on rencontre lorsque des mesures sont réalisées, avec des fibres 10, 12 pour amener la lumière et pour prélever le signal de fluorescence, les extrémités des fibres étant positionnées à la surface d'une paroi 21 qui délimite le milieu (figure 10B). Plus généralement, il s'agit du cas où source et détecteur sont positionnés à la surface de cette paroi 21 qui délimite le milieu diffusant dans lequel le fluorophore F est disposé C'est notamment le cas où on souhaite être non invasif, source et détecteur étant positionnés à proximité de la surface du milieu. Ainsi, dans le cas d'un diagnostic d'une tumeur de la prostate, source et détecteur seront placés en périphérie d'une sonde endorectale. Il s'agit dans ce cas de vérifier non seulement l'équation de propagation mais également les conditions aux limites associées, déjà indiquées ci-dessus . O(r, t) + 2AD(r) ~ (r, t) = 0 pour tout point r de la frontière 5Q du milieu diffusant, et cela pour les longueurs d'onde d'excitation et d'émission. On montre que cette condition est équivalente à écrire, de manière plus simple, ( ~ ( r , t) = 0 dr E a extrapolée où la frontière extrapolée est située à une distance dextrapolée = 2AD de la frontière physique, comme expliqué dans l'article de A.Laidevant et al. � Effects of the surface boundary on the determination of the optical properties of a turbid medium with time resolved reflectance", Applied Optics, 45, 19, 4756 (2006). Une façon de satisfaire à cette condition est d'utiliser la méthode des sources images, comme mentionné dans l'article ci-dessus, mais aussi comme expliqué dans la thèse d'Aurélie Laidevant, � Méthodes optiques résolues en temps pour la tomographie de fluorescence dans les milieux diffusants �, Université Joseph Fourier, Grenoble, 5 octobre 2006, p.55-56: - on connait la solution en milieu infini pour une source placée en r+. Une façon de satisfaire 0(r,t)=0 er ea'ext"polee est de placer une source virtuelle, de contribution négative, à une position r-, symétrique de r+ par rapport à la frontière extrapolée aS2extrapozee Et tout ceci est valable également dans le domaine fréquentiel, et l'on a: 1 e`k(')T+ e'k(w)Y _ 1 v(w)=((w)ù((w)=4TCC DE r r I 4-zen]) [G (k(w)'r+)ùG (k(w),r )l n avec : w k~(w) _ - a +Z D c. D et . : tranformation de Fourier du débit de fluence émis par la source positive dans l'hypothèse d'un milieu infini : tranformation de Fourier du débit de fluence émis par la source négative dans l'hypothèse d'un milieu infini r+ = Irs+ - rl distance à la source excitatrice positive, la source excitatrice positive étant la source excitatrice r_ =lr, rl: distance à la source excitatrice négative, cette source négative étant factice, et définie afin d'obtenir q5(r,t) = 0 Vr e a-1 extrapolée rs+ : coordonnées de la source excitatrice positive rs_ : coordonnées de la source excitatrice négative G�k(w),r) : fonction de Green satisfaisant l'équation k(r,o)r e r D'autre part, on a vu que mi = =ioI'(w)lw=0x _ 1 aw D(w)l w = 0 On en déduit, pour le milieu semi-infini de diffusion en milieu infini. G�k(w),r) = x1/2oo ù 1 eik(0) etk(0)r /laxDx eik(0) + eik(0)r r+ r_ r,+ -r rs_-r (0)) 1 r+G° (k(0), r+) ù rG�k (0), r ) 2c,2V,uaD G�k(0),r+)ùG�k(0),r) x,1/2oo _20 De même, on peut montrer que uam Ird+ r 1 eik(0)r'+ ù eik(0)r'- 1 e Dm .f (rd+ -r , rd_ -r , km(0)) 2cn JJtam1lm f am Ird r m1/20 ù JtamDm eik(0)r'+ eik(0)r'_ n JUamDm ù.V Dm Ird+ rl d- rl e r'+ rd+ -r rd_ ùr Avec : rd+ : position de la source positive, qui est dans ce cas la source fluorescente rd position de la source négative de la source fluorescente, source factice permettant d'obtenir Yom (r, t) = 0 Vr e extrapolée r'+ = rd+ - rl= distance à la source fluorescente positive r'- = rd- - rl= distance à la source fluorescente négative L'utilisation d'une source image en milieu semi-infini est également décrite dans l'article de M.S.Paterson et al. � time resolved reflectance and transmittance for the non invasive measurement of tissue optical properties �, Applied Optics, vol.28, p.2331 - 2336, 1989. Ce document traite également d'une géométrie en forme de � tranche � (� slab �). On peut montrer que pour le signal de fluorescence, on a toujours : 1/2c= x,1/2,o + m,1/2c +1- f (Irs+ - rIrs- -rkx(0)) + f (Ir rd+l,lr km (0)) Vx Vm25 32 Ceci définit une fonction implicite et une surface dans l'espace en trois dimensions. Là encore, et comme déjà expliqué ci dessus, pour trouver la position 3D du fluorophore, on effectue 3 mesures à trois positions de sources et 3 positions de détecteurs. Le fluorophore est donc à l'intersection de 3 surfaces définies dans l'espace par 3 équations différentes. Pour chaque couple de positions, on procède à l'établissement d'une TPSF ou d'un histogramme comme expliqué ci-dessus. De préférence, on sélectionne des couples de positions définissant des directions sensiblement perpendiculaires les unes aux autres. En pratique, les conditions que définit la géométrie en milieu semi-infini ne sont pas très différentes de celles du milieu infini, dès lors que l'on positionne la source et le détecteur (dans le cas de la figure 1A, les extrémités des fibres 10 et 12) loin de l'interface, à plus de 1 cm ou 2 cm de celle - ci (comme en figure 10A). Dans ce cas, en effet, tout se passe comme si le milieu était infini, du point de vue du fluorophore. On peut vérifier la validité de l'équation 1/2 _ x,1/2w + m,1/2w +2 f(rs+ùr,rsùr,kx(0))+.f (rùrd+,rùrd ,km(0))+ Vx Vm en se plaçant dans le cas d'une géométrie infinie. Il n'y a pas de source fictive et on a donc : rs+ = rs, r- = I rs_ - r l = +oo, 5 10 33 r'- = r - rd = +oo, On retrouve alors l'expression correspondant au milieu infini. On va traiter maintenant le cas d'un milieu à faces parallèles (géométrie en forme de � tranche � ou � slab �). Il est déjà question de ce type de géométrie dans l'article de M.S.Paterson et al. déjà cité ci-dessus. Dans ce cas, on montre que, pour le signal d'excitation : -'slab (D, (w) _ fq) s+,l,w)ùtex (r_,i'w) 1 ik(w)Y e e 4,zcD f x[ r+ 15 La sommation dans cette expression résulte du fait que l'on est en présence de deux plans parallèles qui limitent le milieu. Il y a donc contribution de deux sources virtuelles négatives, mais chacun des deux plans multiplie les effets de la source 20 virtuelle négative associée à l'autre plan. Il faut donc ajouter la contribution de ces différents effets, sous la forme d'une somme à l'infini indicée par i. De façon analogue, pour le signal 25 d'émission : 1 e`k(w)r+; e`k(w)r- T1 [ 4~zcnDm r '+i r' i slab (co) = On appelle slab le moment normé d'ordre 1 de gslab(t), pouvant être obtenu par la transformée de Fourier Ît.siab(w) et m1 = = i ffl)slab(w) = OX 1 aw (Dslab(w) w = 0 r.G`°(k(0),r )-r �k(0),r ) f(r- r,r- r,k(0)) G�k(0),r+,l)-G�k(0),r) - 2cn~uaxDx et 1 r'+,zG�k(0),r'+,z)-r's,zG�k(0),r'_z) -~.f(rd+,z -r,rd -r,k(0)) 2cn.J,umnDm G`°(k(0),r'+,l)-G°(k(0),r',l) 2c,,V,uamD Or, slab slab slab 1/2o, = x,1/2. + m,1/2. +i Soit : Dans la pratique cependant, pour des milieux d'épaisseur suffisante, la somme converge assez rapidement, par exemple lorsque i est égal à 10 ou 20 voisin de 10, par exemple i = 8, 9, 11 ou 12. La surface décrite est encore une surface 3D ; Donc, là encore, on localise un fluorophore par recherche des intersections entre 3 surfaces, chacune obtenue pour un couple de point (position source, 25 position mesure) différent. Pour trouver la position 3D du fluorophore, on effectue 3 mesures avec trois positions différentes de sources et 3 positions différentes de détecteurs, comme déjà expliqué ci-dessus. .f(lrsi+ ù rllrs1_ ùrl,kx(0)) + f (Ir ùrai+l,lr - ratù1, km(0)) Vx Vm slab +r slab slab m 15 35 On va traiter maintenant le cas d'un milieu de forme quelconque : ce milieu est limité par une enveloppe, mais la forme de celle-ci est a priori indéterminée La méthode la plus générale pour traiter ce cas est la suivante. On procède d'abord à un maillage du milieu. Les équations sont ensuite résolues de manière numérique (par exemple par la méthode des volumes finis, ou la méthode des éléments finis, etc...). On peut exprimer le signal de fluorescence de la manière suivante : d'm(r,Jd,r,t) x r,t) * F(t) *Om(r,rd,t) , où : F(t) _ _/ e z r 20 25 1 a0x(r,t)+o.(-Dx(r)oçx(r,t))+� (r)0x(r,t)=8(rs -r)8(t), avec la condition aux limites (3) ci-dessus. Et : 1 (r ,t) +V.(ùDm(r)Vgm(r,t))+pan, (r)(5m(r,t) = 8(rd ùr)8(t) avec la condition aux limites (4) ci-dessus. On obtient donc ainsi, en tout point rla valeur des fonctions gx(rs,r,t) et Om(r,rd,t), dont on peut calculer les valeurs des moments, également en tout point r, en appliquant la relation (11) ci-dessus. n at at On a donc la relation : signal mesuré ù x + m + déclin flua = f'(Ir ùr1)+ f'(Ir ùrd1)+z où les fonctions f sont donc calculées numériquement. On cherche ensuite les points r tels que cette relation soit vraie. On peut inclure les paramètres D et ü dans la fonction f' et ainsi obtenir l'expression : signal mesuré = fX' (Dx, }pax, abs (rs-r)) + fm' (Dm, pain, l r-rd l) + i, avec : fX'(Dx, pax,lrs-rl) : solution de l'équation d'excitation fm'(Dm, pain, lr-rdl) solution de l'équation d'émission. Par résolution numérique, on trouve des points r qui satisfont à cette équation et on peut définir une surface passant par tous ces points, par exemple par interpolation d'un point à l'autre. Cette surface n'est pas nécessairement un ellipsoïde. Enfin, on peut traiter le cas où le milieu n'est pas homogène du point de vue de ses propriétés optiques, qui deviennent alors dépendantes de r. Dans ce cas, comme dans le cas d'un milieu homogène de forme quelconque, on procède à la résolution numérique de 3 équations, et on identifie le couple de points satisfaisant aux trois équations. On détermine alors numériquement l'expression de d3(w) pour ce milieu, satisfaisant 37 l'équation de diffusion écrite dans le domaine fréquentiel : \7.(ùD(r)VCD(r, w)) + r Pa (r) ù \ Cn / D(r, co) = S(rs ) Et on en déduit ensuite =ia`l)(0)lw-0x ti 1 aw ^:D(w)lw=0 Quelle que soit la géométrie envisagée, il est possible de s'affranchir de la dépendance par rapport à la durée de vie T, en effectuant une 4ème mesure avec un 4eme couple de positions (position de la source, position du détecteur) différent de chacun des 3 précédents. Pour de couple supplémentaire de positions (source, détecteur), on procède à l'établissement d'une TPSF ou d'un histogramme comme expliqué ci-dessus. On obtient alors l'équation d'une 4ème surface. Comme indiqué ci-dessus dans chacun des cas traités on obtient en fait un calcul du temps moyen à partir de chaque signal de fluorescence. Il est décrit dans le document EP 1 884 765 que l'on peut réaliser le calcul d'une variable indépendante de z en effectuant la différence entre le temps moyen calculé pour chaque signal de fluorescence et le temps moyen calculé pour un signal de fluorescence particulier, par exemple le signal de fluorescence particulier est celui ayant un temps moyen minimal ou ayant un temps moyen calculé minimal. Parmi les 4 mesures effectuées pour les 4 couples de positions, on pourra sélectionner celle qui correspond à un temps moyen mesuré minimal. L'équation 38 de la surface qui correspond à cette mesure minimale sera soustraite de chacune des 3 autres équations. On obtiendra alors 3 nouvelles équations de 3 nouvelles surfaces, qui peuvent ne pas être des ellipsoïdes, mais à l'intersection desquelles sera localisé le fluorophore, ou bien le fluorophore sera localisé dans un volume qui contient cette intersection. Quelle que soit la géométrie envisagée, le calcul de l'intersection de 3 ellipsoïdes, ou 3 surfaces, définis, ou définies, dans l'espace par 3 équations différentes est effectué à l'aide d'un traitement numérique mis en oeuvre par un calculateur, par exemple les moyens 24 des figures 1A et 1B. Ce traitement conduit à une localisation de l'intersection des 3 surfaces avec une certaine précision. Autrement dit, la résolution du système d'équations que constituent les définitions des 3 surfaces est effectuée avec une certaine précision, ce qui donnera une solution approchée : le fluorophore n'est alors pas localisé exactement à l'intersection de 3 surfaces, mais au voisinage de cette intersection. Après avoir trouvé une solution avec une première précision, on peut effectuer de nouveau le calcul pour trouver une autre solution avec une deuxième précision. Chaque surface peut elle-même être définie avec une certaine précision. Ceci définit non plus une surface au sens strict, mais un volume qui s'appuie sur ladite surface. Dans le cas d'un ellipsoïde, il peut s'agir d'une couronne ellipsoïdale Ce volume est 39 délimité par deux surfaces, sensiblement parallèles à la surface définie strictement par l'équation correspondante et proches de celle-ci la proximité étant définie par la précision associée à la surface, qui devra être de préférence inférieure à 25%, et encore de préférence inférieure à + 10% .Autrement dit, on cherche alors non pas l'intersection de 3 surfaces, mais de 3 volumes. Il en résulte non pas un point unique, mais l'identification d'un volume, en général suffisamment petit pour être compatible avec la localisation avec une certaine incertitude; ce volume contient l'intersection des 3 surfaces telles que définies strictement par les 3 équations initiales. D'une manière générale, on peut distinguer différents types d'analyse. Il y a d'abord le cas d'une analyse invasive pour lequel source et détecteur sont placés à l'intérieur du milieu diffusant. Si la profondeur de la source et du détecteur par rapport à la frontière du milieu est suffisante, on considèrera qu'on se trouve dans les conditions d'un milieu infini (cas de la figure 10A). Un premier cas d'analyse non invasive est celui pour lequel source et détecteur sont placés au contact d'une frontière du milieu diffusant (cas de la figure 10B). Un deuxième cas d'analyse non invasive, est celui pour lequel source et détecteur ne sont pas en contact physique avec une frontière 21 du milieu diffusant 20, mais sont en contact optique avec le 40 milieu. Dans ce cas, la source est par exemple un laser 8 focalisé à la surface 21 du milieu et le détecteur est un photodétecteur 4 en contact optique avec la surface du milieu (cas de la figure 10C). Le premier cas non invasif peut-être combiné au second cas non invasif, la source 8 étant située à distance (comme le laser 8 de la figure 10D) et la détection s'effectuant au contact, par exemple à l'aide d'une fibre optique 12 qui amène le rayonnement à détecter vers le détecteur 4 (cas de la figure 10D). Par ailleurs, l'excitation peut s'effectuer au contact, par exemple au moyen de fibres optiques Généralement, en mode non invasif, on ne considère pas qu'on se trouve dans une géométrie infinie, et on choisit alors un autre type de géométrie: la géométrie semi-infinie, ou la géométrie de type tranche, ou la géométrie de forme quelconque. L'intérêt de se rapprocher de la géométrie infinie, semi-infinie ou de type tranche est de permettre d'utiliser une relation analytique reliant une grandeur mesurée (par exemple le temps moyen d'arrivée des photons) aux distances respectives source-fluorophore (ou agrégation ou accumulation locale de fluorophore) et détecteur -flurophore (ou agrégation ou accumulation locale de fluorophore). Ainsi, pour 3 acquisitions différentes, c'est-à-dire réalisée avec 3 couples différents (source détecteur), le problème revient à déterminer les solutions satisfaisant à trois équations analytiques, ce qui peut-être entreprit rapidement à l'aide de 41 moyens de calculs actuels. En revanche, comme décrit ci-dessous, le fait de se placer en géométrie quelconque revient à résoudre le problème de façon numérique, au détriment du temps de calcul. Dans les différents cas expliqués ci-dessus en liaison avec les figures 10A-10D, la référence 20 désigne un milieu étudié, dont la limite est la paroi 21. Ce milieu 20 peut être par exemple un organe d'un animal ou d'un être humain, par exemple le cerveau, ou un sein (à titre d'exemples de fluorophores pour ces différents milieux, on peut citer le vert d'indocyanine, ou la fluorescéine). Sur ces figures, le point F désigne un fluorophore ou une agrégation de fluorophores. Mais on rappelle qu'il peut s'agir aussi d'un absorbeur (ou encore d'un fluorophore ré-émettant à la même énergie que celle d'excitation, et avec une durée 2 = 0) ou d'une agrégation d'absorbeurs. Dans ce cas, le milieu 20 peut être encore un organe, la zone F identifiant une tumeur cancéreuse, visibles par le simple fait qu'elle présente un coefficient d'atténuation plus important que les tissus sains environnants. Des exemples vont maintenant être décrits. Exemples . Exemple I- Cet exemple est un calcul en deux dimensions. Il ne concerne donc pas une mesure réelle, mais illustre la méthode dans un cas théorique, dans un plan. Le milieu est supposé infini, à 2 dimensions, avec les propriétés optiques suivantes. Le coefficient de diffusion réduit est p's = 10 cm-1 ce qui conduit à : D=1/3p', (cette définition de D est la plus répandue, mais il existe d'autres définitions, par exemple D=1/3 (p's+ p,), ce qui ne change rien si Pa �p', qui définit en général le cadre de validité de l'approximation de la diffusion). Le coefficient d'absorption est Pa = 0, lcm 1 . On suppose en outre que les propriétés optiques sont les mêmes aux deux longueurs d'ondes, d'indice de réfraction n=1.0, ce qui donne une vitesse de propagation dans le milieu en = 3.108/n. Il en résulte une vitesse apparente : vapp = 2 .J,uQD = 3.4469x10' m.s-1. n Le fluorophore choisi a un temps de vie z=1.5 ns, ordre de grandeur classique des fluorophores que l'on utilise en imagerie moléculaire optique. On suppose que l'on réalise : - une première mesure ( La figure 3 représente l'intersection (zone III) des deux ellipses I et II ainsi définies, cette intersection donne la position du fluorophore. En deux dimensions, le procédé selon l'invention donne donc des résultats satisfaisants. Des exemples en 3 dimensions vont maintenant être donnés. Exemple II- Le montage expérimental (figure 4A) est sensiblement celui de la figure 1A avec une diode laser 8 de longueur d'onde d'émission à environ 635 nm (diode laser de H&B) en tant que source d'excitation et un photomultiplicateur 4 couplé à une carte TCSPC. La cuve 21 est remplie d'un milieu diffusant liquide 20 composé d'eau et de billes d'Acronal� (résine acrylique, BASF), ce qui conduit aux paramètres optiques suivants: - pa = 0.01 cm-1, - 11's = 8 . 7 cm-1 On considère que les propriétés optiques sont identiques à la longueur d'onde d'excitation et à celle d'émission, du fait de la proximité des longueurs d'onde d'excitation et de diffusion, suffisante pour considérer ces coefficients constants sur cette plage de longueur d'onde. Le milieu fluorescent considéré est du Cy5 de concentration 10pMol. La durée de vie z est de 0,96 ns, et l'indice de réfraction no du milieu est de 1,33. 44 La figure 4A représente de manière plus détaillée le positionnement de l'échantillon dans la cuve 21. Celui-ci est placé dans un tube capillaire 23 en verre, de diamètre interne 1 mm, sur une hauteur d'environ 2 mm. On retrouve également sur cette figure la portion des fibres 10 et 12 qui sont plongées dans le milieu examiné. L'extrémité de chaque fibre est située à 3,5 cm sous le niveau du liquide, ce qui suffit à assurer la condition de milieu infini. La puissance mesurée en sortie de fibre d'excitation 10 est de 100 pW. Sont également représentées sur cette figure les positions respectives de l'échantillon et des extrémités de fibres. Du point de vue de la modélisation, la source n'étant pas ponctuelle, tout se passe comme si elle devenait ponctuelle à une distance d'un libre parcours moyen de transport soit à une distance de 1/p's (ce qui est bien le cas dans la configuration envisagée, puisque la distance entre les extrémités de fibres et l'échantillon est de 0,5 cm, tandis que 1/p's mesl) sont les suivantes (on rappelle que mesl est le moment d'ordre 25 1 mesuré): - Pour (S1, Dl): 2, 2700 cm ; - Pour (S2, D2) : 2,0059 cm ; - Pour (S3, D3) : 2,7284 cm ; - Pour (S4, D4) : 2,3931 cm ; 30 - Pour (S5, D5) : 2,5258 cm. 46 Partant des mesures effectuées pour 3 des 5 couples de points ((Si, D1), (S3, D3) et (S4, D4)), on a calculé trois surfaces tridimensionnelles, en fait trois ellipsoïdes, conformément à l'enseignement de l'invention. L'intersection de ces trois surfaces a également été calculée. Les 3 surfaces ont été élaborées en considérant un maillage : En x : de 50 points entre x -0,5 cm et x=2 cm; En y : de 30 points entre y = -0,5 cm et y=1 cm; En z : de 50 points entre z = -6 cm et z=-3,5 cm. Le fluorophore a été localisé en position (x1r y1r z1) suivante . - 0,041 On rappelle que la position réelle est (0,0, -4). Les figures 4D et 4E représentent dans chacun des plans XY et ZY (donc respectivement en vue de dessus est en vue de côté de la cuve 21) la position p1 mesurée et la position Po calculée, avec un léger décalage entre ces deux positions. La figure 4F représente ces deux positions en trois dimensions, avec, également, un léger décalage entre ces deux positions. Malgré ce décalage, on constate que l'information obtenue par le calcul est tout à fait satisfaisante pour bon nombre de cas, par exemple pour une localisation approchée dans un organe du corps humain. 47 Le temps nécessaire pour ces calculs, pour le maillage choisi (avec le maillage retenu) avec un logiciel � Matlab �, sur processeur Intel Core2, à 2.13GHz, avec 1Go de RAM est inférieur à 3 secondes: - durée du calcul des distances di pour les trois couples (Si, Di) sélectionnés: 2,6 s - durée du calcul des équations des 3 ellipsoïdes (1, 2, et 3) : 0,12 s - durée du calcul de la détermination de l'intersection : 0,04 s. La durée totale de calcul avec un matériel informatique standard est tout à fait compatible avec une mesure sur ou dans un organe du corps humain, par exemple lors d'une intervention chirurgicale ou dans une unité d'analyse. Le calcul a ensuite été refait, en prenant en compte l'ensemble des 5 mesures : on a donc calculé l'intersection des 5 surfaces. Le fluorophore a alors été localisé en position (x2r y2r z2) = (0,06 ; 0,12 ; 3,91). Cette nouvelle position est assez proche de la précédente : la localisation effectuée avec seulement trois surfaces est donc suffisante. Exemple III- Cet exemple est en trois dimensions, et en milieu semi-infini. On considère 3 couples de positions de la source et du détecteur (rs, Par exemple, on a les 3 positions suivantes de la source (toutes les coordonnées sont en cm): 48 rs1=[0.0, 0, 0], rsz= [0.5, 0, 0], rs3= [0.75, -1, 0]; Et les 3 positions suivantes du détecteur (toutes coordonnées également en cm): rd1=[1.0, 0, 0], rdz= [1.5, 0, 0], rd3= [0.75, 0,5, 0] Ces 3 couples sont représentés en figure 5A. Les propriétés optiques sont les mêmes que dans l'exemple précédent. On suppose que l'on réalise : - une première mesure cm, pour le couple de positions (rs1, une deuxième mesure ( ,oùr) X Vapp = 4.6000 15 cm, pour le couple de positions (rs3, Alors la position du fluorophore appartient à l'intersection des ellipsoïdes, qui ont pour équations suivantes: f(Irsz+ rl,lrsz_ ùrl,k(0)) + f (Ir ùraz+I,Ir ùrdz_I,k(0)) =Vapp x ( (t),nesi ùr) = 4.6 moment a'ordre 1 mesuré= ml 73 .f (Irs3+ -rI,Ir0 - -rI,k(0))+ f (Ir ùrd3+I,Ir ùrd3_I,k(0)) =Vapp x ( (t)mes1 ùr) =5 moment d'ordre 1 mesurxl .f (fr ù r',ksi _ ùri, k(0))+.f (Ir ùrd1+I,Ir ùrd1_I,k(0)) = v,pp x �1') mesl ùr)=4.61 moment d'ordre 1 mesuré=ml Sur chacune des figures 5B, 5C, 5D est représentée une partie de chacun des ellipsoïdes correspondants. ( , -z-) x Vapp = d1) r 4.6000 20 25 49 La zone d'intersection de ces 3 ellipsoïdes est identifiée par la zone I en figure 6A, sur laquelle sont en outre représentées les trois positions S1, S2, S3 de la source et les trois positions Dl, D2, D3 du détecteur utilisées pour les trois mesures. La zone P est la zone de localisation du fluorophore. Les figures ont été élaborées en considérant un maillage : - En x : de 101 points entre x = -2 et x=3; - En y : de 101 points entre x = -2 et x=2; - En z : de 101 points entre x = -3 et x=0; La solution de la figure 6A est une solution approchée, à 10% près. Il est possible d'imposer une précision plus importante. Ainsi, en figures 6B et 6C, sont représentées des solutions (zone P) plus précises : En figure 6B : solution à 2% près, En figure 6C : solution à 1,4% près. On voit que l'accroissement de précision apporte peu d'information supplémentaire pour une localisation approchée du fluorophore. Il est intéressant d'évaluer le temps nécessaire pour ces calculs, avec un logiciel 25 � Matlab �, sur processeur Intel Core2, à 2.13GHz, avec 1Go de RAM: - pour le calcul des moments pour les 3 couples source-détecteur en chaque point du maillage : 5.6835 s, 30 - pour le calcul des équations des 3 ellipsoïdes (1, 2, et 3): 2.086609 s, 15 20 - pour le calcul de la détermination de l'intersection : 2.253164 s. Il en résulte une durée totale d'environ 10 s pour le maillage choisi. Ce résultat illustre bien le caractère � temps réel � ou rapide du procédé. Dans le cadre d'une application chirurgicale, un praticien peut réaliser une analyse selon l'invention alors même qu'il est en présence d'un patient. Une durée de 10 s, ou même une durée légèrement supérieure, est tout à fait compatible avec le maintien d'un patient sur une table d'opération ou dans une salle d'analyse. La détermination de la position du fluorophore dépend de son temps de vie. Pour s'en affranchir, on peut considérer des mesures différentielles, comme expliqué dans la demande de brevet EP 1 884 765. Par exemple : Avec les configurations de sources-détecteurs ci-dessus, en milieu semi-infini, on a les ellipsoïdes suivants : i f (Irs,+ rl,lrsl_ rl,k(~)) + Î (Ir ùrdi+I,Ir ùrdi_I, k(0)) = v vv x �tL ùr)=4.7 ù,r mome td ordre 1mesuré =mi z : f (Irsz+ ùrhlrsz_ -rl,k(0))+ f(Ir ùraz+I,Ir ùrd2_ k(0)) = v x ( (t)mesl -r)=5 ,pp moment d'ordre 1 mesuré=ml Le Sème correspond à une mesure dont le trajet optique est le plus court des trois. C'est cette mesure qui, comme expliqué dans la demande EP 1 884 765, peut être utilisée comme mesure de référence. : f (Irs3+ ùrl, Ir 3- -Ilk(0)) + f (Ir ùrd3+I, Ir -rd3_I, k(0)) = v pP x ( (t) mesl ù a) = 4.6 ù,,. moment d'ordre 1 mesuré=ml 25 51 Si l'on considère les quantités 13 = 0.1 et z3 =0.4, on obtient des relations indépendantes du temps de vie. Ces quantités ne définissent plus des ellipsoïdes mais des surfaces plus complexes. Les figures 7A et 7B représentent, suivant deux orientations différentes, l'intersection entre les deux surfaces 3D S1_S3 (zone I sur ces figures) et S2_S3 (zone II sur ces figures). En 3D, une troisième surface est nécessaire pour identifier une localisation possible du fluorophore. On considère donc un ensemble, non pas de 3 mesures, mais de 4 mesures. Par exemple, on a les 4 positions suivantes de la source (toutes les coordonnées sont en cm): rs1=[0.0, 0, 0], rsz= [0.5, 0, 0], rs3= [0.75, -1, 0]; rs4= [-0.75 ; -1.5 ; 0]; Et les 4 positions suivantes du détecteur (toutes coordonnées également en cm): rd1=[1.0, 0, 0], rdz= [1.5, 0, 0], rd3= [0.75, 0, 5, 0] ; rd4= [0.75, 1, 0 ] . On obtient donc, outre les 3 surfaces E1r z, 3 déjà définies, une surface E4 (qui définit elle aussi un ellipsoïde) . : f (Irs4+ ùr1>is4- -rI,k(0)) + f (ir -rd4+1>Ir -rd4_i> k(0)) = vnPP X ( ~t~mesl ù,,.ù moment d'ordre 1 mesuré=ml ùr) =5.030 52 Cette surface E4 va elle aussi être normalisée. On effectue donc la différence U-U qui définit elle-même une surface dont l'intersection avec les surfaces U_ U et S2-S3 va permettre de localiser le fluorophore. La figure 7C représente les 3 surfaces 1- 3, (zone I) , S2-S3 (zone II) et S4-S3 (zone III), tandis que la figure 7D représente une coupe de l'ensemble dans un plan XZ La zone d'intersection des 3 surfaces, er donc la zone de localisation du fluorophore est identifiée par la zone P sur cette figure 7D. Exemple IV- Cet exemple est en trois dimensions, et en géométrie de type � slab �. La tranche considérée est d'épaisseur L = 1, 5 cm. On considère 3 couples de positions de la source et du détecteur (rs, rd) Par exemple, on a les 3 positions suivantes de la source (toutes les coordonnées sont en cm): rs1=[0.0, 0, 0], rs2= [0.5, 0, 0], rs3= [0.75, -1, 0]; Et les 3 positions suivantes du détecteur (toutes coordonnées également en cm): rd1=[1.0, 0, 0], rd2= [1.5, 0, 0], rd3= [0.75, 0, 5, 0] Comme pour les exemples précédents, on obtient la définition de trois ellipsoïdes, pour les trois valeurs suivantes de vapp (Ct>mes) : 1,84 ; 1,93 ; 1,81. 53 Sur chacun des figures 8 A, 8B, 8C est représentée une partie de chacun des ellipsoïdes correspondants. Sur ces figures, la tranche constituant le milieu est délimitée par les deux plans parallèles: Z = 0 et Z = - 1,5 cm. La zone d'intersection de ces 3 ellipsoïdes est identifiée par la zone I en figure 8D, sur laquelle sont en outre représentées les trois positions de la source et les trois positions du détecteur, utilisées pour les trois mesures. De manière générale, un exemple de procédé selon l'invention peut être mis en oeuvre selon le schéma représenté en figure 9. Les étapes de cet exemple, et d'ailleurs de toute autre mise en oeuvre d'un procédé selon l'invention, peuvent être précédées d'une étape - non représentée - d'injection dans le milieu étudié d'au moins un fluorophore ou absorbeur destiné à venir se fixer sur une zone d'intérêt. En fait lorsqu'on injecte un fluorophore, on en injecte souvent plusieurs. On peut d'ailleurs assimiler une zone fluorescente à une agglomération de fluorophores. Il n'y a généralement pas qu'une seule molécule fluorescente, mais une concentration locale de molécules fluorescentes. De la même façon, lorsqu'on travaille en absorption, il n'y a pas qu'une seule molécule absorbante, mais une concentration de molécules absorbantes. Cette concentration peut être considérée comme ponctuelle si la distance la séparant du détecteur est suffisamment grande devant ses dimensions. Typiquement, une 54 concentration de rayon r est considérée comme ponctuelle, à une distance d'observation d donnée, si d > 10 r. Certains absorbeurs peuvent déjà être dans le milieu, par exemple lorsqu'il s'agit d'une zone présentant une absorption différente de son environnement. On positionne ensuite les moyens d'émission de rayonnement et de détection par rapport au milieu. L'une des configurations déjà décrites peut être sélectionnée. En particulier on peut réaliser un examen invasif u non invasif. Au cours d'une première étape S1 de ce procédé, les données, ou paramètres, d'entrée sont fournis: il s'agit des propriétés optiques du milieu diffusant, aux deux longueurs d'onde de d'excitation et d'émission (dans certains cas, à ces deux longueurs d'onde, les propriétés optiques sont identiques), la géométrie du milieu, les positions des sources et du détecteur (au sens déjà expliqué ci-dessus), et le choix du maillage pour le calcul des surfaces. Au cours d'une deuxième étape (S2), les fichiers de mesures (TPSF) pour chaque couple source-détecteur sont lus. Un prétraitement des données peut être réalisé. A partir des données introduites lors de la première étape, une expression analytique de chaque surface va pouvoir être déterminée ou calculée pour chaque point du maillage (étape S3). Les données issues de l'étape S2 vont permettre un calcul des temps moyens pour chacune 30 des mesures effectuées (étape S4). 55 Les étapes S1 et S3 d'une part, S2 et S34 d'autre part, sont représentées en figure 9 comme sensiblement réalisées simultanément, mais ce n'est pas nécessaire, quoique souhaitable en vue d'un déroulement le plus rapide possible du procédé. A partir des données de calcul des étapes S3 et S4, on recherche les points du maillage tels que Si = , pour tout i, à x% près (étape S5), ce qui signifie 1 Au besoin, les calculs peuvent être réitérés en fixant une valeur de X supérieure (auquel cas on réduit la précision) ou inférieure (auquel cas on accroît la précision) à la valeur précédente (étape S7).20 REVENDICATIONS1. Procédé de localisation d'au moins un fluorophore (22) ou d'au moins un absorbeur dans un milieu diffusant (20), à l'aide d'une source (8, 10) de rayonnement apte à émettre un rayonnement d'excitation de ce fluorophore ou de cette molécule et des moyens (4, 12) de détection apte à mesurer un signal de fluorescence (q)f1uo) émis par ce fluorophore (22) ou un signal d'émission émis par cet absorbeur, comportant: a)- pour au moins 3 couples différents de positions de la source de rayonnement et des moyens de détection, au moins une excitation par un rayonnement provenant de la source (8) de rayonnement, et au moins une détection par les moyens (4) de détection du signal de fluorescence émis par ce fluorophore après cette excitation, ou du signal d'émission émis par cet absorbeur, b) pour chacun de ces couples, l'identification d'une surface sur laquelle le fluorophore ou l'absorbeur est situé, ou d'un volume comprenant cette surface et dans lequel le fluorophore ou l'absorbeur est situé, c) une estimation de la localisation du fluorophore ou l'absorbeur dans son milieu diffusant, par calcul de l'intersection des trois surfaces, et éventuellement d'un volume autour de cette intersection.30 57 2. Procédé selon la revendication 1, dans lequel le milieu environnant le fluorophore ou l'absorbeur est de type infini ou de type semi-infini ou forme une tranche, limitée par deux surfaces parallèles, ou est de forme quelconque, sa surface extérieure étant discrétisée en une série de plans. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel le milieu diffusant est de type infini, les surfaces ayant alors la forme d'ellipsoïdes. 4. Procédé selon l'une des revendications 1 à 3, l'étape b), comportant, pour chaque couple de position, un calcul d'un moment temporel normé, d'ordre n (n> 1), d'une fonction temporelle du signal de fluorescence ou du signal d'émission par l'absorbeur, ou d'un moment normé d'ordre n (n> 1) d'une fonction fréquentielle, déduite de ladite fonction temporelle du signal de fluorescence ou du signal d'émission par l'absorbeur. 5. Procédé selon la revendication 4, la fonction fréquentielle étant déduite de ladite fonction temporelle par transformée de Fourier ou par transformée de Laplace ou par transformée de Mellin 6. Procédé selon la revendication 4 ou 5, dans lequel le moment temporel normé d'ordre n de ladite fonction temporelle ou le moment normé d'ordre n est le moment d'ordre 1. 58 7. Procédé selon l'une des revendications 1 à 6, dans lequel on réalise une excitation par la source (8) de rayonnement, et une détection par les moyens (4) de détection pour 4 couples différents de positions de la source de rayonnement et des moyens de détection. 8. Procédé selon la revendication 7, dans lequel, pour la localisation d'au moins un fluorophore 10 (22) : on sélectionne parmi les quatre couples de points, le couple de points, dit quatrième couple de points, pour lequel la durée moyenne d'arrivée des photons, depuis l'instant d'émission par l'émetteur à 15 l'instant de réception par le récepteur, est la plus courte, et on calcule, pour ce couple de points, l'équation d'une surface sur laquelle le fluorophore est situé, - et, pour chacun des trois autres couples 20 de points, on calcule une première équation d'une première surface sur laquelle le fluorophore est situé, et on effectue une différence entre cette équation et celle associée au quatrième couple, pour obtenir une deuxième équation d'une deuxième surface indépendante 25 de la durée de vie du fluorophore 9. Procédé selon l'une des revendications 1 à 7, dans lequel, pour la localisation d'au moins un fluorophore (22), l'étape b), comporte, pour chaque 30 couple de position, un calcul d'une équation de chaque surface, indépendante de la durée de vie du fluorophore 59 10. Procédé selon la revendication 9, ladite équation indépendante de z résultant de la différence entre le temps moyen mesuré ou calculé pour chaque couple de position et le temps moyen mesuré ou calculé pour un quatrième couple de position. 11. Procédé selon la revendication 10, le quatrième couple de position étant un couple pour lequel le signal de fluorescence a un temps moyen mesuré ou un temps moyen calculé inférieur à celui des 3 couples de position. 12. Procédé selon l'une des revendications 1 à 11, comportant en outre une représentation visuelle ou graphique de la position ou de la répartition du ou des fluorophores ou de l'absorbeur. 13. Procédé selon l'une des revendications 1 à 12, les signaux étant détectés par technique TCSPC ou par caméra. 14. Procédé selon l'une des revendications 1 à 13, la source de rayonnement et les moyens de détection étant respectivement l'extrémité d'une fibre optique (10) qui amène la lumière d'excitation dans le milieu et l'extrémité d'une fibre optique (12) qui prélève une partie de la lumière d'émission. 15. Dispositif pour la localisation d'au 30 moins un fluorophore ou d'au moins un absorbeur dans un milieu diffusant, comportant :60 - une source (8, 10) de rayonnement apte à émettre un rayonnement d'excitation, - des moyens (4, 12) de détection apte à mesurer un signal de fluorescence (Q)fluo) émis par ce 5 fluorophore (22) ou un signal d'émission émis par cet absorbeur, - des moyens (24) de calcul, pour : a) pour chaque couple de points parmi au moins 3 couples différents de positions de la source de 10 rayonnement et des moyens de détection, calculer ou déterminer, après une excitation réalisée par un rayonnement provenant de la source (8) de rayonnement, et une détection par les moyens (4) de détection des signaux émis par le fluorophore ou l'absorbeur après 15 cette excitation, une surface sur laquelle le fluorophore ou l'absorbeur est situé, ou un volume comprenant cette surface et dans lequel le fluorophore ou l'absorbeur est situé, b) pour estimer la localisation du 20 fluorophore ou de l'absorbeur dans son milieu environnant, par calcul de l'intersection des trois surfaces ou d'un volume autour de cette intersection. 16. Dispositif selon la revendication 15, 25 les moyens de calcul étant programmés pour calculer, pour chaque couple de position, un moment temporel normé, d'ordre n (n> 1), d'une fonction temporelle, ou un moment normé d'ordre n (n> 1) d'une fonction fréquentielle, déduite de ladite fonction temporelle 30 d'émission.61 17. Dispositif selon l'une des revendications 15 ou 16, les moyens de calcul étant programmés pour réaliser un calcul d'une équation d'une surface, indépendante de la durée de vie du fluorophore 5 18. Dispositif selon la revendication 17, les moyens de calcul étant programmés pour réaliser un calcul d'une équation d'une surface, indépendante de la durée de vie, par différence entre le temps moyen 10 mesuré ou calculé pour chaque couple de position et le temps moyen mesuré ou calculé pour un quatrième couple de position. 19. Dispositif selon la revendication 17 15 ou 18, le quatrième couple de position étant un couple pour lequel le signal de fluorescence a un temps moyen mesuré ou un temps moyen calculé inférieur à celui des 3 couples de position. 20 20. Dispositif selon l'une des revendications 15 à 19, lesdits moyens de détection comportant des moyens de type TCSPC ou des moyens de type caméra. 25 21. Dispositif selon l'une des revendications 15 à 20, comportant en outre des moyens (27) de représentation visuelle ou graphique de la localisation du fluorophore ou de la molécule absorbante. 30