FR 2968921 A1 20120622 FR 1060567 A 20101215 1 PROCEDE DE LOCALISATION D'UN MARQUEUR OPTIQUE DANS UN MILIEU DIFFUSANT DESCRIPTION DOMAINE TECHNIQUE ET ART ANTERIEUR L'invention concerne le domaine de l'imagerie optique, par exemple de fluorescence, dans les milieux diffusants et notamment les tissus biologiques. Elle met en oeuvre des méthodes optiques résolues en temps. Elle trouve une application au domaine médical, en particulier à la tomographie optique diffusive de fluorescence pour la détermination de la distribution de particules fluorescentes, par exemple des marqueurs. L'imagerie optique diffuse, et notamment l'imagerie de fluorescence, systèmes de diagnostic offre la perspective de non invasifs grâce a l'utilisation de rayonnements non-ionisants, faciles d'utilisation et peu coûteux. Dans les applications mettant en oeuvre l'imagerie de fluorescence, des marqueurs fluorescents sont injectés au sujet et se fixent sur certaines molécules spécifiques, par exemple des tumeurs cancéreuses. La zone d'intérêt est éclairée à la longueur d'onde d'excitation optimale du fluorophore (substance chimique d'une molécule capable d'émettre de la lumière de fluorescence après excitation) et le signal fluorescent est détecté. 2 On cherche ensuite, à partir de ce signal de fluorescence à obtenir des informations de localisation ou de concentration sur les marqueurs. L'imagerie optique de fluorescence (avec injection de fluorophore spécifique) est limitée à ce jour à des applications « petit animal » en raison du manque de marqueurs adaptés et injectables à l'humain, et du problème d'autofluorescence des tissus qui se pose pour la détection en profondeur. En effet, pour appliquer cette méthode au diagnostic de cancers chez l'homme, il est primordial que le signal spécifique situé plus profondément sous la peau que chez le petit animal puisse être détecté. Mais le signal spécifique à détecter faiblit avec la profondeur, principalement à cause de l'absorption et de la diffusion des tissus, et se confronte à un signal parasite qui perturbe la détection. Ce signal, appelé « auto-fluorescence », décrit la fluorescence de tissus auxquels aucune substance chimique ou fluorophore spécifique n'a été injecté : il s'agit de la fluorescence naturelle du tissu. Lorsque les marqueurs sont situés en profondeur (quelques cm), leur signal spécifique s'atténue. Alors, l'autofluorescence peut devenir d'amplitude non négligeable car elle est émise par tout le volume des tissus biologiques. De plus, le filtrage de la lumière d'excitation n'est jamais parfait et un peu de lumière d'excitation peut passer à travers le filtre de fluorescence censé la supprimer et venir s'ajouter au signal mesuré. D'autres éléments du système de mesure 3 (filtres par exemple) peuvent également avoir une faible fluorescence propre. Ces deux phénomènes, autofluorescence propre des tissus et signal parasite dû au système, sont des phénomènes qui viennent interférer et fausser la mesure de fluorescence. L'autofluorescence est un phénomène connu, mais est à ce jour rarement perçue comme un signal parasite. En cancérologie notamment, l'autofluorescence est utilisée pour distinguer des tissus cancéreux de tissus sains. Il ne s'agit alors pas d'injecter de marqueur spécifique, mais simplement d'observer l'autofluorescence de zones spécifiques et de comparer différentes zones d'un même individu. Par opposition, la spectroscopie optique de fluorescence utilise généralement des longueurs d'onde d'excitation en rouge ou proche infra rouge, qui assurent une absorption moindre du signal dans les tissus, et permettent ainsi une meilleure pénétration des tissus. L'autofluorescence des tissus est alors beaucoup plus faible et devient un signal à supprimer plutôt qu'à utiliser. Dans la suite on pourra regrouper les deux effets que sont l'autofluorescence propre des tissus et la production de signaux parasites dus au système sous le terme de signaux parasites. On connait différentes méthodes de soustraction d'images pour limiter l'effet de l'autofluorescence dans le cas d'une illumination du milieu par une source continue. 4 Selon une première de ces méthodes on réalise une soustraction simple du signal moyen d'une région de l'image qui ne contient pas de marqueurs fluorescents. Une limitation de ces techniques est la nécessité de trouver une zone de référence sans fluorescence représentative et de reproduire des conditions d'acquisition similaires. Selon un autre type de méthode on effectue une soustraction en utilisant des filtres « de fond » décalés vers le bleu : cette technique est mise en place dans les systèmes IVIS de Caliper Life Science. Une image est prise avec un filtre décalé vers le bleu dont le rôle est d'exciter l'autofluorescence, mais pas les fluorophores. Cette image de fond est ensuite soustraite - à un facteur d'échelle près - à l'image régulière avec le filtre d'excitation. Il n'est pas certain que cette fluorescence produite par un décalage vers le bleu soit représentative de l'autofluorescence à la longueur d'onde d'observation. On trouve aussi des méthodes à base de techniques spectrales, mais elles nécessitent évidemment de réaliser des acquisitions résolues en longueur d'onde. Ainsi une technique multispectrale avec des excitations à différentes longueurs d'onde a été décrite dans le document US 2005065440. D'autres travaux, voir notamment ceux de D.Wood et al. « Refining epifluorescence imaging and analysis with automated multiple-band flat-field correction", paru dans Nature Methods Application Notes (2008), ou S. Psycharakis et al., Autofluorescence removal from fluorescence tomography data using multispectral imaging (Proceedings Paper SPIE 66260I-66260I-7 (2007)) utilisent également une méthode multispectrale basée sur des modèles a priori de 5 spectres. Certains de ces travaux ont présenté des résultats en tomographie du petit animal. Enfin, il existe des méthodes de décomposition spectrale (mesure à une seule longueur d'onde) pour distinguer l'autofluorescence de la fluorescence. Un spectromètre permet de mesurer le spectre d'émission en chaque pixel de l'image. La décomposition sépare les deux spectres et permet d'obtenir deux images, une d'autofluorescence, une de fluorescence. Certaines de ces techniques effectuent une séparation de spectres par analyse en Composantes Principales (ACP), d'autres par Factorisation en Matrices Non Négatives. Là encore, la méthode ne s'applique qu'à des illuminations continues c'est-à- dire non pulsées, la factorisation en matrices non négatives permettant de reconstruire le spectre en longueur d'onde de différentes composantes constituant un signal de fluorescence détecté. Par spectre en longueur d'onde, on entend une distribution en énergie. (US 20090245611). D'autres techniques mettent en oeuvre une illumination du milieu par une lumière pulsée, également appelée fluorescence résolue en temps. Par lumière pulsée, on entend une impulsion lumineuse dont la durée est de quelques fs à quelques ns. Dans ce cas, la détection du signal émis par le tissu s'effectue 6 généralement entre deux impulsions consécutives de la source lumineuse. Ainsi, dans le document A.T. Kumar et al. « Feasibility of in vivo imaging of fluorescent proteins using lifetime contrast", Optics Letters 34, 2066-2068 (2009), on distingue la fluorescence de l'autofluorescence (AF) sur un petit animal en modélisant la distribution de l'AF et le déclin de fluorescence du fluorophore. Mais cette technique met en oeuvre une connaissance a priori des profils de déclin pour faire un ajustement sur les mesures temporelles. Cette technique permet de réaliser seulement une image en deux dimensions, mais ne permettent pas d'accéder à une information sur la localisation en profondeur des espèces d'intérêt. Enfin, dans la technique décrite dans l'article de A. Laidevant et al. « Analytical method for localizing a fluorescent inclusion in a turbid medium », Applied Optics, 46, 2131-2137 (2007), une mesure spécifique du signal parasite est effectuée sur le milieu en retirant l'inclusion fluorescente. Cette mesure est ensuite soustraite au signal de fluorescence brute. Cette technique est évidemment coûteuse en temps (elle nécessite de réaliser une mesure sans inclusion fluorescente puis une mesure avec inclusion fluorescente) et est délicate à mettre en oeuvre. D'une manière générale, il se pose donc le problème de trouver un nouveau procédé, ne nécessitant pas la mise en oeuvre d'une technique spectrale, pouvant être réalisée rapidement (en vue d'application in vivo) et permettant de différencier, les différentes 7 contributions de sources lumineuses, et notamment les contributions de sources fluorescentes et de sources parasites, en faisant appel à une technique de fluorescence résolue en temps. Il se pose également le problème de trouver un nouveau dispositif, permettant de mettre en oeuvre un tel procédé. EXPOSÉ DE L'INVENTION L'invention concerne d'abord un procédé de traitement de donnée relatives à au moins un marqueur optique dans un milieu diffusant, dans lequel : a) on réalise au moins une acquisition d'un signal optique résolu en temps, ce signal résultant d'une excitation du milieu à l'aide d'une source de rayonnement en impulsions, chaque acquisition comportant d'une part une composante d'intérêt, due audit marqueur, et d'autre part une composante parasite, b) on constitue une matrice X à partir de chaque acquisition réalisée à l'étape a). On peut ensuite traiter les données de ladite matrice X par factorisation de cette matrice en un produit de deux matrices non négatives A et S. Il est ensuite possible d'extraire, ou de déterminer ou de calculer ou d'identifier, la contribution d'une desdites composantes à au moins une desdites acquisitions, à partir des termes des matrices A et S. Selon un mode de réalisation, un tel 30 procédé peut en outre comporter une étape de calcul ou de détermination ou d'identification ou de localisation 8 de la position d'au moins un marqueur optique dans le milieu diffusant et/ou une représentation graphique de la répartition d'au moins un tel marqueur dans le milieu diffusant, à partir d'au moins une desdites composantes temporelles acquises. Selon ce mode de réalisation, l'invention permet donc une identification de la contribution au signal détecté, mais aussi une identification (même approchée) de la position, d'un ou de plusieurs marqueurs optiques dans le milieu diffusant étudié, en trois dimensions, sans avoir recours à une décomposition spectrale du rayonnement de fluorescence, et sans nécessiter de mesures temporelles censées ne représenter que le signal parasite (sans le marqueur). Un marqueur peut être une molécule présentant les propriétés optiques souhaitées. Selon un mode de réalisation, au moins un des marqueurs optiques est un marqueur fluorescent. Dans ce cas, au cours de l'étape a), on détecte un signal de fluorescence, dont une longueur d'onde correspond à au moins une longueur d'onde d'au moins un marqueur fluorescent, le milieu diffusant étant alors excité par une source de rayonnement en impulsions, la source émettant à une longueur d'onde d'excitation de la fluorescence du marqueur. La composante d'intérêt, due à au moins un des marqueurs fluorescents, est alors au moins une composante de fluorescence (mais il peut y en avoir plusieurs). 9 La composante parasite peut comprendre une composante d'autofluorescence du milieu, ainsi qu'une composante de diffusion de la lumière d'excitation. On peut donc ensuite extraire ou calculer ou déterminer ou identifier une contribution de fluorescence due à au moins l'un des marqueurs fluorescents à partir des termes des matrices A et S. On peut ensuite déterminer ou calculer la position de l'un au moins des marqueurs fluorescents, dans le milieu diffusant, et/ou une représentation graphique d'une répartition d'un ou de plusieurs marqueurs fluorescents, à partir de la contribution de fluorescence obtenue lors de l'étape d). Selon un autre mode de réalisation, au moins un des marqueurs optiques présente des propriétés d'absorption et/ou de diffusion différentes de celles du milieu dans lequel il est placé. Au moins une des composantes d'intérêt comporte alors un rayonnement absorbé et/ou diffusé. Dans ce cas, le rayonnement détecté est un rayonnement de diffusion émis par le milieu, en réponse à une excitation lumineuse impulsionnelle. Dans ce cas le rayonnement diffusé par le marqueur présente des caractéristiques temporelles distinctes de celles du rayonnement diffusé par le milieu. On peut ensuite déterminer ou calculer la position de l'un au moins des marqueurs, dans le milieu environnant, et/ou une représentation graphique d'une répartition d'un ou de plusieurs marqueurs, à partir des composantes d'intérêt. 10 Un procédé selon l'invention peut comporter, avant formation dudit tableau multidimensionnel, une étape de correction des données de signal optique de chaque acquisition par recentrage des données temporelles dans une plage temporelle autour d'un temps moyen du signal optique. Dans l'un ou l'autre des procédés ci-dessus, la première matrice A non négative du produit AS peut être une matrice dont les éléments ax,p sont des coefficients de pondération, ax,p étant la contribution, dans une acquisition indicée par l'indice x, d'une composante d'intérêt, par exemple de fluorescence, ou d'une composante de diffusion et/ou d'autofluorescence, ces composantes étant désignées par l'indice p. La seconde matrice non négative S peut comporter le profil temporel d'intérêt, par exemple de fluorescence, de chaque composante. Chaque ligne p de cette matrice peut représenter un profil temporel de signal caractéristique de la composante p (des valeurs numérisées de ce profil temporel). Le nombre de lignes du tableau S et le nombre de colonnes du tableau A correspondant alors au nombre de composantes à séparer. Le tableau X est formé en réalisant des acquisitions successives, une acquisition pouvant par exemple correspondre à une position donnée de la source et à une position donnée du détecteur, ou à la position d'une source donnée et à la position d'un détecteur donné (ici et dans le reste du texte, la position d'une source (respectivement : d'un détecteur) peut-être la position de l'extrémité distale d'une fibre qui transmet un rayonnement depuis ladite source vers le milieu étudié (ou depuis ce milieu vers le détecteur). Ainsi, la matrice X comprend les données temporelles mesurées, ces dernières étant également appelées acquisitions, indicées par l'indice x et notées X. Chaque acquisition correspond à un histogramme temporel du signal détecté par le détecteur. Les matrices X, A et S sont généralement de dimension 2. Lors de l'étape de traitement du tableau X comprenant les données temporelles mesurées, ces dernières étant également appelées acquisitions, on peut déterminer A et S par minimisation d'une fonction de coût F, cette fonction pouvant être : - le carré de la distance euclidienne entre la matrice X et le produit matriciel A.S : on a alors F - X -AS 2 ou bien encore la fonction : F -X-AS 2+a S-So 2 où So est la matrice initiale de S et a un coefficient réel dit de régularisation. La matrice X, rassemblant une ou plusieurs acquisitions, est traitée de préférence selon un procédé itératif. Par exemple, on réalise k itérations, les tableaux Ai et Si, obtenus lors de l'itération d'ordre i étant réactualisés à chaque itération. Le nombre d'itération peut être déterminé en fonction de l'évolution de la fonction de coût F entre des itérations successives, ou en fonction de la valeur de cette fonction de coût, en la comparant alors à un critère d'arrêt prédéterminé. Autrement dit, dans ce 12 dernier cas, le processus itératif s'arrête lorsque la fonction de coût atteint ledit critère d'arrêt. Ce nombre d'itérations peut également être déterminé empiriquement, en fonction de l'expérience de l'utilisateur. L'excitation du milieu peut être réalisée par une source d'excitation laser, qui peut être éventuellement focalisée à l'interface entre le milieu diffusant et le milieu extérieur. La lumière d'excitation va ensuite pénétrer dans le milieu diffusant, et y exciter des marqueurs ou des sources dans ce milieu, par exemple à 3 cm ou 5 cm en profondeur, c'est-à-dire à distance de l'interface, dans le milieu diffusant. Le rayonnement de fluorescence provient donc d'une zone en profondeur, par exemple comprise entre l'interface et environ 3 cm ou 5 cm à distance de l'interface, ou comprise entre 1 cm à distance de l'interface et 5 cm à distance de l'interface. L'excitation peut avoir lieu dans l'infra- rouge ou le proche infra rouge ou le rouge, par exemple à au moins une longueur d'onde comprise environ entre 600 et 900 nm. Selon un mode de réalisation, la fluorescence d'un marqueur fluorescent peut quant à elle être détectée à des longueurs d'onde supérieures à la longueur d'onde d'excitation du milieu. L'invention concerne également un dispositif de localisation d'un marqueur optique dans un milieu diffusant, ou de traitement de données relatives à un tel marqueur, comportant : 13 a) au moins une source de rayonnement en impulsions pour produire un faisceau d'excitation du milieu diffusant, b) des moyens pour réaliser au moins une acquisition d'un signal de diffusion optique résolu en temps d'au moins un marqueur optique dudit milieu, l'acquisition comportant d'une part une ou plusieurs composantes d'intérêt, et d'autre part une ou plusieurs composantes parasite, c) des moyens pour constituer une matrice X à partir de chaque acquisition, d) des moyens pour traiter les données de ladite matrice X par factorisation en deux matrices non négatives A et S, e) des moyens pour calculer ou déterminer ou identifier ou extraire la contribution d'une desdites composantes à au moins une desdites acquisitions à partir des termes des matrices A et S. Selon un mode de réalisation, un tel dispositif comprend en outre des moyens pour déterminer ou calculer ou identifier et/ou pour visualiser la répartition des intensités du ou des marqueurs optiques à partir d'au moins une contribution d'une desdites composantes. Le marqueur optique peut être un marqueur fluorescent. Dans ce cas, la composante d'intérêt associée à ce marqueur est une composante de fluorescence. Aussi, la source d'excitation produit des impulsions à la longueur d'onde d'excitation du marqueur fluorescent, et les moyens de détection sont 14 adaptés pour collecter le rayonnement de fluorescence du marqueur fluorescent. Selon ce mode, la source d'excitation permet d'éclairer le milieu, et de produire une zone, dite zone d'excitation dans ce dernier. En effet la lumière d'excitation pénètre dans le milieu, y diffuse, et excite les sources de fluorescence, marqueurs et tissus autofluorescents. Comme expliqué ci-dessus, les sources de fluorescence peuvent être localisées en profondeur, à distance sous l'interface. Un dispositif selon l'invention peut comporter en outre des moyens pour modifier la position de la source d'excitation et/ou des moyens de détection. Une acquisition du signal optique de diffusion résolu en temps, par exemple un signal de fluorescence, peut être réalisée pour chaque position de la source d'excitation et des moyens de détection. Ou bien un dispositif selon l'invention peut comporter plusieurs fibres optiques qui permettent d'amener le rayonnement de la source d'excitation en différents points du milieu et/ou diverses fibres optiques qui permettent de capter, en différents points du milieu ou en différents points de sa limite, un rayonnement optique, par exemple un rayonnement de fluorescence, qui provient de ce milieu, et de l'amener à un ou plusieurs détecteurs. Les moyens pour traiter la matrice d'acquisition X par factorisation en deux matrices non négatifs A et S, mettent en oeuvre un procédé selon l'invention, comme déjà décrit ci-dessus. 15 BRÈVE DESCRIPTION DES DESSINS - La figure 1 représente un dispositif pour mettre en oeuvre l'invention, - les figures 2A et 2B représentent un autre dispositif pour mettre en oeuvre l'invention, - les figures 3A-3B des acquisitions, d'une part de la fluorescence brute mesurée et, d'autre part, du signal parasite mesuré, pour deux couples source - détecteur, - la figure 4 est un exemple de superposition d'un ensemble d'acquisition de fluorescence brute, pour 6 sources et 4 détecteurs, - la figure 5 représente schématiquement la décomposition, - la figure 6 illustre comment est constituée une acquisition de fluorescence, - la figure 7 représente un milieu et des zones à propriétés optiques différentes du milieu environnant, - les figures 8A-8B et 9A-9B représentent des comparaisons entre des composantes de fluorescence mesurées et des composantes de fluorescence calculées, ainsi qu'entre des composantes parasites mesurées et calculés par la mise en oeuvre d'un procédé selon l'invention, - les figures 10A, 10B et 11A, 11B représentent des images de reconstruction, d'une part sans correction de la composante parasite (figure 9 A) puis avec correction de la composante parasite (figure 9 B) et, d'autre part, avec correction expérimentale de 16 la composante parasite (figure 10A) et avec correction selon la présente invention (figure 10B). EXPOSÉ DÉTAILLÉ DE MODES DE REALISATIONS DE L'INVENTION Dans cet exposé, on va tout d'abord considérer le cas où le (ou les) marqueurs optiques sont des marqueurs fluorescents. La figure 1 est un exemple de système expérimental permettant de mettre en oeuvre l'invention. L'illumination d'un milieu 20 à examiner est réalisée à l'aide d'une source 8 de rayonnement qui produit des impulsions de rayonnement, par exemple un rayonnement dans l'infra rouge ou même le proche infra rouge. Cette source de rayonnement est impulsionnelle, chaque impulsion de rayonnement pouvant avoir une largeur à mi-hauteur comprise entre, par exemple, quelques nanosecondes, par exemple 1 ns ou 10 ns, et quelques femtosecondes, par exemple 10 fs. Il peut par exemple s'agir d'un laser fonctionnant de manière impulsionnelle. Le faisceau de cette source de rayonnement peut être dirigé vers le milieu, par exemple avec une fibre optique 10. En variante, plusieurs fibres peuvent être utilisées pour envoyer des impulsions à divers endroits dans le milieu, définissant alors autant de sources de lumière, chaque extrémité de fibre pouvant être assimilée à une source de lumière s. Des moyens de commutation, par exemple une platine de translation, permettent alors la sélection d'une fibre parmi les différentes fibres. 17 Une platine permet de positionner la source de rayonnement devant telle ou telle fibre. Le rayonnement d'excitation passe de préférence par un filtre interférentiel 14 pour éliminer une grande partie de la lumière à une longueur d'onde supérieure à la longueur d'onde d'excitation. La lumière d'excitation diffuse dans le milieu diffusant 20 et va y exciter une ou plusieurs espèces fluorescentes. En réponse, le milieu émet un rayonnement comprenant : - une composante de fluorescence, due à la présence du ou des marqueurs fluorescents 22, cette composante étant parfois appelée fluorescence spécifique car elle provient des marqueurs fluorescents. C'est cette composante que l'on cherche généralement à déterminer, de façon à localiser la position des marqueurs fluorescents dans le milieu, - une composante parasite. Cette dernière peut notamment comporter une composante d'autofluorescence du milieu, et/ou une composante de diffusion du rayonnement d'excitation dans le milieu et/ou une ou plusieurs composantes parasites dues par exemple à la fluorescence propre de certains éléments matériels du système, par exemple des composants optiques. Une ou plusieurs fibres 12 collectent la lumière en provenance du milieu 20 étudié, dite lumière émise. Un filtre interférentiel 16 et/ou un filtre coloré peuvent être placés devant le détecteur 4 pour limiter la détection de lumière en dehors du spectre de 18 fluorescence des fluorophores 22 disposés dans le milieu 20 et optimiser l'élimination de la lumière d'excitation. Le filtre peut par exemple être un filtre passe haut, laissant passer les longueurs d'onde supérieures à 650 nm lorsque la source émet à la longueur d'onde 631 nm. Le rayonnement émis est donc envoyé vers des moyens de détection 4 et d'acquisition 24. Les moyens de détection 4 comprennent un détecteur de photons, relié à des moyens d'acquisition du signal détecté, ces derniers pouvant comporter des cartes d'acquisition d'un ordinateur. Ces moyens d'acquisition 24 permettent d'obtenir la distribution temporelle du signal détecté par les moyens de détection 4. Par distribution temporelle, on entend le temps, dit temps d'arrivée, séparant la détection d'un photon d'un temps initial to déterminé. Ce temps initial correspond généralement à l'instant timp où la source émet une impulsion lumineuse, éventuellement avec un léger décalage temporel. Une telle distribution est généralement discrétisée dans le temps et se présente sous la forme d'un histogramme communément dénommé par l'acronyme TPSF (Temporal Point Spread Function). D'une façon générale, les moyens de détection et d'acquisition permettent d'obtenir un signal de fluorescence résolu en temps, ce signal prenant généralement la forme de l'histogramme précédemment décrit. Eventuellement le rayonnement émis est transmis vers les moyens de détection 4 via une ou plusieurs fibres optiques 12, comme déjà expliqué 19 ci-dessus. Dans ce cas, l'extrémité de chaque fibre optique située du côté du milieu diffusant peut être assimilée à un détecteur. Les moyens de détection 4 et d'acquisition 24 comportent par exemple un photomultiplicateur couplé à une carte de comptage, cette dernière produisant une distribution temporelle des signaux détectés par le photomultiplicateur. Les moyens de détection 4 et d'acquisition 24 peuvent également comporter une caméra rapide intensifiée, cette caméra étant activée selon une porte temporelle dont l'ouverture est successivement décalée, ce qui permet de réaliser simultanément un histogramme de type TPSF en différents points de l'espace. Une ligne de synchronisation 11 permet de déclencher les moyens de détection 4 et/ou d'acquisition 24 à l'aide d'une impulsion engendrée par la source de rayonnement 8. Ainsi, ces moyens de détection 4 et d'acquisition 24 permettent d'obtenir, un signal de fluorescence résolu en temps, ou histogramme temporel du signal de fluorescence. Lorsque le milieu est excité par une source s et que le signal émis par le milieu est détecté par le détecteur d, cet histogramme est noté Xsd. Il peut être également être appelé acquisition résolue en temps correspondant au couple source-détecteur. Ainsi, une acquisition résolue en temps peut être référencée par XX ou Xsd, l'indice x étant un entier avec 1 Autrement dit, des fluorophores situés à une distance z de la frontière 20' du milieu comprise 21 entre 0 cm (donc situés très proches de la surface) et, par exemple 3 cm ou 5 cm vont pouvoir être détectés. Les moyens de détection 4 et d'acquisition 24, déterminent la distribution temporelle Xsd d'un rayonnement qui provient de la zone du milieu diffusant excitée par la source lumineuse s, qui traverse le milieu diffusant en direction de la frontière 20' entre le milieu diffusant et le milieu extérieur, puis qui atteint les moyens 4 de détection comprenant le détecteur d. Une analyse spectrale n'est pas nécessaire pour les besoins de la présente invention. Typiquement, le milieu étudié peut être un milieu vivant. Il peut s'agir par exemple d'une zone du corps humain ou animal. L'enveloppe corporelle constitue l'interface du milieu diffusant avec le milieu extérieur. Une source d'excitation est donc focalisée sur cette interface, ou dans le milieu. Des marqueurs préalablement injectés dans ce milieu diffusant permettent de localiser des zones telles des tumeurs. Comme déjà expliqué ci-dessus, il y a également excitation d'autres éléments du milieu, créant une composante parasite au signal détecté. Un autre exemple de dispositif selon l'invention est illustré en figure 2A, dans une application à un milieu qui est en fait un fantôme de prostate 20 contenant des marqueurs fluorescents 22. Un laser 8 fonctionne en impulsions, par exemple à une longueur d'onde proche de 775 nm. Les impulsions envoyées sont très brèves : elles sont par exemple de l'ordre de quelques dizaines de picosecondes (par 22 exemple de largeur temporelle à mi-hauteur comprise entre 10 ps et 50 ps). La sortie du laser 8 est fibrée et envoyée dans 6 fibres sources 101r 102, 103, 104, 105, 106. Ces fibres sont intégrées dans une coque de sonde endorectale 30 en contact avec le fantôme. Pour chaque fibre source (10i, i=1-6), la lumière se propage dans le milieu 20 jusqu'à la zone fluorescente. Les marqueurs fluorescents 22 ainsi excités émettent des photons à une longueur d'onde supérieure à celle d'excitation, ici aux environs de 800 nm. Mais il y a dans le signal d'autres contributions parasites, notamment de l'autofluorescence qui provient des tissus biologiques environnant les marqueurs fluorescents. En utilisant un fantôme tel que celui représenté sur la figure 2A, il est possible d'obtenir un signal d'autofluorescence en effectuant la soustraction de deux acquisitions réalisées respectivement avec et sans l'inclusion fluorescente. Cette sonde contient également 4 fibres de détection 12j (j=1-4) pour collecter le signal de fluorescence. Ces fibres acheminent la lumière détectée vers le détecteur d compris dans les moyens de détection 4, relié à une carte de comptage de photons, faisant partie des moyens d'acquisition 24, qui permet de construire l'histogramme des temps d'arrivée des photons ou distribution temporelle des photons. Les moyens de détection 4, et d'acquisition 24 ont déjà été décrits ci-dessus. Des filtres peuvent aussi être utilisés, comme déjà expliqué ci-dessus. 23 La figure 2B représente l'extrémité 30' de la sonde 30, avec les extrémités Si - S6 des fibres d'excitation 10i, i=1-6 et celles (D1 - D4) des fibres de détection 124 (j=1-4). À titre d'exemple, l'extrémité de la sonde 30 peut avoir une section sensiblement rectangulaire, par exemple de côté 20 mm et 18 mm. Les extrémités des fibres d'excitation et de détection sont réparties selon deux rangées, de manière alternée, sur chaque rangée, entre les fibres d'excitation et les fibres de détection, par exemple selon un pas d'environ 5 mm. On a donc, sur chaque rangée, les extrémités de 5 fibres, dont 3 fibres d'excitation et 2 fibres de détection. Pour chaque fibre d'excitation 10i, i=1-6 on a pu réaliser des mesures par chaque fibre de détection 124 (j=1-4). Ainsi, dans l'exemple de 6 fibres d'excitation et de 4 fibres de détection, on a pu réaliser 4 mesures (une pour chaque fibre de détection 124 (j=1-4)) avec une excitation par la fibre 10i, i=1-6 (i=1 - 6). Dans la suite de la description, le terme source désigne à la fois une source non fibrée, ou l'extrémité d'une fibre d'excitation lorsque la source est fibrée. De même, le terme détecteur désigne à la fois un détecteur, un groupe de pixels d'un détecteur pixelisé, ou l'extrémité d'une fibre de détection, optiquement couplée à un détecteur. On identifie dans la suite (ou sur les figures) par « source n, détecteur p » (ou sn, dp) une acquisition obtenue par excitation via source 24 d'excitation 10n et par détection via le détecteur 12p. Aussi, chaque acquisition est notée XsdP. Dans la suite de la description, XsdPest un vecteur-ligne, chaque ligne comportant des données représentant l'histogramme temporel de fluorescence, c'est-à-dire une distribution temporelle des photons détectés. XsdP pourrait également être représenté sous la forme d'un vecteur colonne. Chaque acquisition XsdP est constituée de termes notés XsdP,t, t représentant le pas de discrétisation temporel, avec 1 Pour chaque configuration, correspondant à un couple source détecteur, on mesure : - une acquisition dite de fluorescence brute, provenant du milieu avec l'inclusion fluorescente ; Le terme brut désigne le fait que le signal détecté comprend une composante de fluorescence, ou composante d'intérêt, et une composante parasite. Ce signal résolu en temps correspond à l'acquisition Xsd, - une acquisition constituée uniquement de ladite composante parasite, ce signal étant mesuré sans l'inclusion fluorescente ; notons que, en général, la 25 détermination d'une telle composante n'est pas possible, car il n'est pas possible de retirer les marqueurs fluorescents du milieu, - une acquisition Xdiff constituée de la diffusion du signal d'excitation dans le milieu (à la longueur d'onde d'excitation, sans filtre de fluorescence). Quelques exemples d'acquisitions obtenues avec ce dispositif sont illustrés en figures 3A - 3B : - pour la source 2 et le détecteur 2 (figure 3A), - pour la source 4 et le détecteur 2 (figure 3B). Chacune de ces figures montre simultanément le signal (résolu en temps) de fluorescence brut et le signal (résolu en temps) parasite. Ce dernier est parfois faible par rapport au signal de fluorescence brute (cas de la figure 3A) mais parfois important (cas de la figure 3B). On voit qu'il y a donc bien ici une perturbation des mesures de signal de fluorescence par les signaux parasites d'autofluorescence et/ou de diffusion. Pour traiter ce problème, on peut réaliser d'abord un recalage temporel des données d'acquisition temporelles des signaux de fluorescence brute Xsd, par rapport à un temps de réponse impulsionnel. Cette étape de correction est optionnelle, mais sera de préférence réalisée. Les inventeurs ont en effet constaté qu'un tel recalage temporel de chaque acquisition Xsd, réalisé préalablement à la décomposition de la matrice X, permettait d'améliorer les résultats de la factorisation de cette matrice. Un objectif de ce recalage temporel peut être de s'affranchir des variations des fonctions de réponse impulsionnelle du dispositif, ou de l'ensemble constitué par le dispositif et le milieu, entre les différents couples sources-détecteurs. En effet, dans une configuration telle que décrite dans la figure 2A, les sources et les détecteurs sont fibrés. Selon la longueur d'une fibre d'excitation 10i séparant le Laser 8 de l'extrémité de la sonde 30, selon la longueur d'une fibre d'émission 12i s'étendant entre l'extrémité de la sonde 30 des moyens de détection 4, et selon les positions des extrémités desdites fibres d'excitation 10i et d'émission 12i sur la sonde 30, une même impulsion du Laser ne sera pas détectée au même moment selon qu'on utilise un couple fibre d'émission - fibre de détection (ou couple source-détecteur) ou un autre couple. Il est alors utile de réduire les décalages temporels entre les différents couples sources-détecteurs lorsque le milieu est excité par une même impulsion. Afin de corriger un tel décalage peut par exemple calculer tout d'abord le temps de réponse impulsionnel T diffsd du signal de diffusion mesuré: J X diffsd (t) t dt T _diffsd = JXd~gsd(t)dt 0 Ainsi, à chaque couple source-détecteur correspond un temps moyen du signal de diffusion T diffsd, Ce temps moyen quantifie, pour une configuration source-détecteur donnée, la durée moyenne entre l'émission lumineuse et la détection des photons de diffusion par le détecteur. Lors d'une telle acquisition, on peut retirer le filtre 16 disposé en amont du détecteur 4, de telle sorte que le moyen de détection 4 détecte un signal à la longueur d'onde d'excitation, et non à la longueur d'onde de fluorescence comme ce sera le cas lors des acquisitions de fluorescence. On peut calculer ensuite le temps moyen théorique du signal diffusion, dans le milieu diffusant, ppelé de , prédit par le modèle de diffusion choisi (le « chapeau (^) » désigne le modèle, par opposition à la mesure). Par exemple, dans un milieu infini, ce temps moyen est déterminé en divisant la distance entre une source et un détecteur par la vitesse de la lumière dans le milieu. On rappelle que par source et détecteur, on entend ici une extrémité de la fibre d'excitation 10i et une extrémité de la fibre d'émission 12i , les deux extrémités étant fixées sur la sonde 30. Dans notre cas, on calcule T_diffsd de façon numérique. Pour ce faire, on peut résoudre numériquement l'équation suivante : G(» -O.(DDÔGs~~)+µQG,,» =8(i - is)+ s c Avec, Gs'= moment d'ordre 1 de la fonction de Green exprimant le transfert d'énergie entre la source et le détecteur. e moment d'ordre 0 de la fonction de Green exprimant le transfert d'énergie entre la source et le détecteur. 28 g, : coefficient d'absorption du milieu diffusant c : vitesse de la lumière D: coefficient de diffusion du milieu rs : coordonnées de la source, cette dernière pouvant être assimilée à l'extrémité de la fibre d'émission lorsque la source est fibrée. Finalement, on définit T impulsesd comme la difference des deux temps précédemment définis : T_impulsesd =T_dlfJsd -T -dl. sd T impulse sd représente le temps de réponse impulsionnel du dispositif, c'est-à-dire le temps moyen s'écoulant entre une impulsion de la source et la 15 détection des photons diffusés par le milieu diffusant. Ce temps moyen de réponse impulsionnel quantifie le délai entre l'émission d'une impulsion et la détection d'un signal de diffusion, corrigé de l'estimation du délai de diffusion dans le milieu. Autrement dit, 20 T impulsesd quantifie le temps de réponse impulsionnel du dispositif, tandis que T diffsd quantifie le temps de réponse impulsionnel de l'ensemble dispositif et milieu diffusant. D'une façon generale, on retiendra que T impulsesd et 25 T diffsd sont des temps de réponse impulsionnel, quantifiant un délai moyen entre l'émission d'une impulsion et la détection d'un signal de détection. Le recalage temporel des données mesurées peut être par rapport à un temps de réponse impulsionnel, qui peut10 29 être soit T diffsd, soit, et de préférence, 'impulse sd , tels que précédemment définis. L'objectif du recalage temporel est de permettre ensuite de représenter chaque acquisition selon un axe temporel, de telle sorte que, chaque temps de réponse impulsionnel considéré correspondant à la configuration sd a une même coordonnée selon l'axe des temps. Autrement dit, si X'sd désigne une acquisition Xsd temporellement recalée par rapport à un temps de réponse impulsionnel, alors : (t) = Xsd (t - T impulsesd) t) = Xsd (t - Tdiffsd) Autrement dit, chaque acquisition Xsd est recalée dans le temps en fonction du temps de réponse impulsionnel, déterminé pour chaque configuration source-détecteur. . On peut également parler de centrage de chaque acquisition Xsd selon le temps de réponse impulsionnel correspondant au couple source détecteur sd : en effet, selon l'équation précédente, pour chaque couple source détecteur, Xsd' (T impulse sd) = 0 (ou Xsd (T diffsd) = O. De cette façon, on limite les décalages temporels entre les signaux résolus en temps correspondant à chaque couple source-détecteur. Dans la suite de la description, on va décrire plus précisément la factorisation de la matrice ou 30 d'acquisition X, constituée d'une pluralité d'acquisitions Xsd. Le procédé décrit s'applique indifféremment à une matrice X constituée d'acquisitions Xsd recalées temporellement ou non. De préférence, les acquisitions Xsd sont recalées temporellement selon l'équation précédente. La figure 4 est un exemple d'une superposition d'un jeu de 24 mesures Xsd de fluorescence, pour 6 sources et 4 détecteurs. Les acquisitions sont recalées, comme expliqué ci-dessus, par rapport au temps moyen de la mesure de diffusion. Il reste des décalages entre les positions temporelles des différents pics, dus aux positions relatives des couples source-détecteurs et des marqueurs fluorescents : le trajet entre l'extrémité d'une fibre d'excitation et un marqueur n'est pas le même que celui entre l'extrémité d'une autre fibre d'excitation et le même marqueur, et le problème est le même avec les différentes positions des détecteurs (ou des extrémités des fibres de détection) par rapport à ce même marqueur. On explique maintenant comment on exploite ces données, de préférence corrigées ou recentrées, en vue de réaliser une image des marqueurs fluorescents dans le milieu. L'étape (optionnelle) précédente permet de recentrer les données sans s'affranchir de l'autofluorescence, ou, plus généralement, des signaux parasites, tandis que les étapes suivantes permettent de séparer ces composantes. Dans la suite de la description, on suppose que chaque acquisition est une acquisition corrigée en 31 temps. Mais une telle correction n'est pas toujours à mettre en oeuvre. Considérons une acquisition de fluorescence Xsd (c'est-à-dire un ensemble de données résultant d'un échantillonnage d'une courbe de fluorescence telle que l'une de celles des figures 3A- 3B) sur un milieu diffusant contenant un type de marqueur fluorescent. Cette mesure est la superposition de deux distributions temporelles distinctes : - si : distribution temporelle des photons d'intérêt émis par les marqueurs fluorescents, dite composante temporelle de fluorescence ou composante d'intérêt, - s2 : distribution temporelle de photons dits parasite, incluant l'autofluorescence et/ou de diffusion, dite composante temporelle parasite. Si on appelle al et a2 les quantités respectives de ces deux distributions, l'acquisition peut s'écrire comme la combinaison linéaire suivante des 2 distributions s1 et s2: Xsd = al s1 + a2 s2 = (al a =AS (1) où chaque terme Sp,t de la matrice S est la valeur du profil de la composante p correspondant à l'intervalle temporel t. Chaque intervalle temporel t correspond généralement à la fréquence de l'échantillonage temporel mis en oeuvre lors de l'acquisition30 32 Ainsi . - si,i désigne la ième valeur de la distribution si, - s2,i désigne la ième valeur de la distribution s2. Selon cette méthode on cherche à trouver des matrices A et S à coefficients non-négatifs, dont le produit est le plus proche possible de la matrice Xsd On peut généraliser cet exemple de traitement de données à une série X de NX acquisitions. X est alors une matrice non-négative XE RN»'T , dont chaque terme XX,t correspond à la valeur de l'acquisition XX dans le canal temporel t. Comme précédemment indiqué, l'indice x peut être remplacé par l'indice sd, ce dernier représentant la position d'une source s et d'un détecteur d. NX représente alors le nombre de couples source-détecteur considérés. On cherche alors à trouver les matrices non-négatives AE RNXXNP et SE RNPXN` telles que : XAS (2) Par matrice non-négative on entend une matrice dont tous les éléments sont non-négatifs, et Np représente le nombre de composantes temporelles, ou de distributions temporelles si, recherchées. A est dite matrice de poids. Un élément asd,P (> 0) de cette matrice représente la contribution, dans une acquisition de fluorescence Xsd, d'une composante temporelle de fluorescence ou d'une 33 composante temporelle de diffusion et/ou d'autofluorescence, ces composantes étant désignées par l'indice p (par exemple p = 1 désignant la composante de fluorescence et p = 2 désignant la composante parasite). S est la matrice des composantes de X à séparer les unes des autres (dans cet exemple, composante de fluorescence ou composante parasite), chaque ligne de cette matrice, notée Sp, comportant des données discrétisées représentant un profil temporel de la composante p. Si on considère uniquement 2 composantes temporelles (p = 2), cette équation peut être illustrée selon le schéma de la figure 5 pour le cas de 24 mesures. Le problème est ensuite de trouver des matrices particulières A et S pour satisfaire l'équation ci-dessus. Pour cela, on cherche à minimiser une fonction de coût F, par exemple le carré de la distance euclidienne entre X et AS, ou une adaptation de la divergence de Kullback-Leibler. On peut également utiliser une autre expression, comprenant le carré de ladite distance euclidienne, détaillée dans la suite de la description. Si on considère la distance euclidienne entre la matrice X et le produit des deux matrices A et S, alors on cherche à minimiser la fonction de coût F: OX - AS avec A>-0 et S>-0. 34 La solution est obtenue de manière approchée, par itérations. Différentes méthodes peuvent être utilisées pour minimiser la fonction de coût. On peut citer les lois de mise à jour multiplicatives introduites par Lee et Seung dans l'article « Learning the parts of objects by non negative matrix factorization», paru dans Nature, pages 788-791, 1999. Selon cette technique la procédure générale d'optimisation peut donc s'écrire - 1. initialisation de A et S par des valeurs aléatoires positives, - 2. mise à jour de A - 3. mise à jour de S - 4. recommencer les étapes 2 et 3 jusqu'à ce que la fonction de coût F, ou sa variation entre deux itérations successives, atteigne un critère de convergence, ou jusqu'à un certain nombre d'itérations. Dans le cas de Nx acquisitions (par exemple : les 24 mesures obtenues avec le dispositif des figures 2A et 2B), chacune étant discrétisée suivant Nt points dans le temps, on a donc : xNx 1 XNx,Nt ) aNx,1 aNx,2 Si l'indice p = 1 désigne la composante temporelle de fluorescence, et que l'indice p = 2 désigne la composante temporelle parasite, alors le produit ax,1 * S1 correspond à l'histogramme de la X= XII x1,Nt a11 a12 ( S11 S1,Nt S21 S2,Nt 35 composante temporelle de fluorescence correspondant à l'acquisition X. Cette composante temporelle de fluorescence peut être notée IX (ou Isd lorsque l'indice désigne un couple source-détecteur). Cette composante IX correspond à l'acquisition Xx corrigée de sa composante parasite. Autrement dit, IX représente la partie de Xx réellement due aux marqueurs fluorescents, d'où la dénomination composante d'intérêt. On peut également écrire que Ix représente le signal "net" de fluorescence, alors que Xx représente le signal "brut" (i-e non corrigé) de fluorescence. Le produit ax,2 * S2 correspond, lui, à la composante temporelle parasite de l'acquisition Xx On retrouve bien l'équation déjà présentée ci-dessus ax,1 * S1 + ax,2 * S2 XX. On voit donc qu'à partir des matrices A et S, on obtient une composante temporelle, désignée p, d'une acquisition XX, en multipliant les termes de la ligne Sp par le coefficient (ou poids) ax,p . Selon la valeur de p, cette composante peut être une composante d'intérêt (composante de fluorescence dans cet exemple) ou une composante parasite. D'une façon plus générale, cette composante p est obtenue en multipliant des termes de la matrice S, relatifs à cette même composante, par le terme ax,p, de la matrice A, représentant le poids de la contribution p dans l'histogramme X. On a donné en figure 5 l'exemple imagé du 30 produit d'une matrice S (pour une acquisition avec deux 36 sources fluorescentes) avec un tableau A afin d'obtenir le tableau X. On a donné ci-dessus un exemple de fonction de coût F à minimiser. Mais celle-ci peut avoir différentes formes. On peut aussi utiliser, de façon avantageuse le carré de la distance euclidienne entre X et AS, plus un terme de contrainte : F= X-AS 2+a S-So 2 (1) où So est la matrice S initiale, comme expliqué ci-dessous, et a un paramètre dit de régularisation. Ce dernier permet de contraindre plus ou moins la solution S à rester proche de l'initialisation choisie Sc. Le fait d'utiliser une telle fonction de coût permet d'accroître la robustesse de la méthode. Cela peut également permettre d'accroître la vitesse de convergence Un exemple d'algorithme pouvant être utilisé pour calculer le couple de matrices A et S dont la produit représente au mieux les données initiales X est alors, là encore, un algorithme itératif. Le nombre d'itérations N1t est déterminé soit de façon empirique, soit par un critère d'arrêt choisi, comme précédemment évoqué. L'algorithme comprend les étapes suivantes . 1. Initialisation des matrices A et S avec des valeurs strictement positives. On peut choisir effectuer l'initialisation de S avec des modèles de sources mesurés expérimentalement ou tirés de la littérature. 2. Choix du paramètre de régularisation a. Le choix se fait de façon empirique suivant la confiance que l'on donne à l'initialisation (plus a est grand, plus la matrice S calculée est contrainte à être proche de l'initialisation). Il y a ensuite mise à jour de A et S, à tour de rôle pendant Nit itérations. Plusieurs méthodes peuvent être mises en oeuvre pour cette mise à jour. Par exemple, on peut citer une descente de gradient alternativement sur A et S (appelée « ALS » pour « Alternative Least Square ») ou les lois de mises à jour multiplicatives indiquées par Lee et Seung dans l'article déjà ci-dessus. En s'inspirant de ces dernières, et en définissant XT la transposée d'une matrice X, si la fonction de coût à minimiser est F = X - AS 2 +OE S -So z , ce qui constitue la variante préférée, chaque itération, comprend les mises à jour suivantes pour A et S : (AZT X + OGSo ) SZ+i ~ SZ (AzT AzS +OESz ) puis . AZ+1(XSZ+IT ) Z (AZSZ+1SI+1T ) Ces lois de mise à jour particulières ont été développées par les inventeurs. On peut choisir une fonction de coût F plus classique, précédemment évoquée, telle que F = X - AS 2 Dans ce cas, à chaque itération, les lois de mise à jour sont celles décrites par Lee et Seung : (AZTX) (AZT AZS) Puis xs +iT AZ+1 AZ T (AZSZ+isl+i ) A partir des données des matrices A et S, on peut extraire, pour chaque acquisition XX, une composante de fluorescence, ou composante d'intérêt, due aux marqueurs, et notée Ix. On peut alors mettre en oeuvre une méthode de reconstruction basée sur le modèle de l'équation de la diffusion, par exemple telle que décrite dans J.Boutet et al. dans « Bimodal ultrasound and fluorescence approach for prostate cancer diagnosis", J. Biomed. Opt. 14, 064001 (2009). Pour chaque courbe temporelle obtenue, on peut calculer les moments d'ordre 0 et 1 du signal de fluorescence corrigé lx(t) : M° =J~Ix(t).dt 0 Mx = fo lx (t).t.dt 0 L'obtention de tels moments permet de déterminer la position des fluorophores, selon des techniques telles que par exemple décrites dans l'article de J.Boutet et al. cité ci-dessus, ou encore dans le document "Méthodes optiques résolues en temps pour la tomographie de fluorescence dans les milieux (2) (3) 39 diffusants",Thèse de Aurélie Laidevant soutenue le 5 octobre 2006 à l'Université Joseph Fourier - Grenoble - France ou encore dans la demande EP10179184 du 24 septembre 2010. Afin de tester les différentes étapes d'un procédé selon l'invention, on a pu réaliser différents tests. Pour chaque couple source-détecteur, les comparaisons entre les signaux temporels mesurés (moins le signal parasite de diffusion) et les signaux calculés par le procédé selon l'invention ont été réalisées en utilisant un fantôme de fluorescence. Des exemples de résultats sont donnés sur les figures 8A - 8B (figure 8A : pour le couple source 3, détecteur 3, figure 8B : pour le couple source 4, détecteur 2). Sur ces figures, le signal mesuré correspond à : - un premier signal mesuré lorsqu'une inclusion fluorescente est présente dans le fantôme, - auquel on a soustrait un second signal, ce dernier étant mesuré lorsque l'inclusion fluorescente est retirée et représentant donc la composante parasite. Sur ces figures, le signal mesuré est identifié par Smes et le signal calculé par Sea (composante fluorescente du premier signal décrit dans l'alinéa précédent, déterminée selon le procédé objet de l'invention). Les deux signaux sont pratiquement confondus La soustraction permet d'évaluer la composante de fluorescence du premier signal. 40 Dans toutes les comparaisons effectuées, dont deux exemples sont en figures 8A - 8B, on a une composante de fluorescence expérimentale (signal mesuré décrit précédemment) et une composante calculée (signal calculé décrit précédemment). Dans la plupart des cas, l'accord entre les valeurs de signal mesuré et le signal calculé est très bon. De même, pour chaque couple source- détecteur, des comparaisons ont été effectuées entre les signaux parasites mesurés et les signaux parasites calculés par le procédé selon l'invention. Deux exemples en sont donnés sur les figures 9A - 9B (figure 9A : pour le couple source 3, détecteur 3 ; figure 9B : pour le couple source 4, détecteur 2). Sur chacune de ces figures : - le signal mesuré Smes correspond au signal mesuré lorsque l'inclusion fluorescente est retirée. Il représente la composante parasite. - le signal calculé Scat est la composante parasite d'un signal mesuré lorsque l'inclusion fluorescente est placée dans le fantôme. Cette composante parasite est déterminée selon le procédé objet de l'invention. Dans certains cas, et en particulier sur les figures 9A - 9B, les deux signaux sont pratiquement confondus. Cela permet d'obtenir une composante parasite expérimentale (signal mesuré décrit précédemment) et calculée (signal calculé décrit précédemment). 41 On a pu constater, dans la plupart des cas un bon accord entre les mesures et le calcul, ce qui permet une validation expérimentale de la méthode. Une comparaison avec les résultats obtenus avec d'autres méthodes a été réalisée. Ainsi les figures 10A et 10B montrent respectivement une comparaison entre une localisation par un procédé sans correction puis par un procédé selon l'invention, donc avec correction. Dans cette expérience, la position réelle du marqueur fluorescent selon l'axe Z, par rapport à la limite du milieu diffusant est z = 2 cm. En appliquant une technique classique et sans correction de l'autofluorescence on obtient z = 1.7 cm (figure 10A). Avec correction expérimentale, on obtient : z = 2.12 cm (figure 10B). Lors de la correction expérimentale, on détermine expérimentalement la composante parasite (mesure sans inclusion fluorescente) et on la soustrait des mesures réalisées avec l'inclusion fluorescente La comparaison des figures 10A et 10B montre que, sans correction de l'autofluorescence, l'inclusion est reconstruite à une position trop proche de la sonde (comme sur les figures 11A-11B, les sources S sont identifiées sur ces figures 10 par les disques noirs de faible diamètre et les détecteurs D par des disques noirs de plus grand diamètre). Les figures 11A et 11B montrent respectivement une comparaison entre une comparaison dite « expérimentale » c'est-à-dire réalisée en soustrayant une mesure sur un fantôme réalisée sans l'inclusion fluorescente d'une mesure sur fantôme 42 réalisée avec l'inclusion fluorescente et une comparaison par un procédé selon l'invention. Plus précisément, la figure 11B montre un résultat obtenu avec correction par la méthode de séparation des composantes selon le procédé précédemment décrit en utilisant une fonction de coût basée sur la distance euclidienne pour laquelle on obtient z=2.19 cm. Les résultats de la correction expérimentale (figure 11A) et de la correction selon l'invention (figure 11B) sont similaires. Des étapes d'un procédé selon l'invention sont représentées en figure 6 : - dans une étape S1, on réalise une ou plusieurs acquisitions par excitation du milieu diffusant, par un faisceau laser fonctionnant en impulsions; il en résulte par exemple une ou plusieurs acquisitions, - lors de l'étape S2, on forme la matrice X, avec, pour chaque ligne, des données numérisées et échantillonnées d'une courbe de signal de fluorescence pour un couple (source, détecteur) ou un couple (fibre d'excitation, fibre de détection) donné, - dans une étape S3, optionnelle, on réalise une correction des données temporelles mesurées 25 précédemment pour le signal de diffusion, - l'équation X A.S peut ensuite être résolue, de manière itérative comme expliqué ci-dessus (étape S4), - on peut ensuite déterminer les 30 acquisitions I correspondant à la composante 43 recherchée, en l'occurrence la composante de fluorescence des marqueurs, - on peut alors procéder à une représentation graphique du ou des marqueurs fluorescents (étape S5). Lors d'une mise en oeuvre d'un procédé selon l'invention sur un patient, l'étape d'excitation par les sources de rayonnement peut être précédée par une étape d'injection d'un marqueur fluorescent dans le milieu dont la zone à étudier. On cherche à localiser ce ou ce(s) marqueur(s) fluorescent(s) dans ce milieu diffusant. Dans toute la description ci-dessus, les acquisitions Xx étaient représentées sous la forme de matrice ligne. Mais l'invention s'applique aussi lorsque ces acquisitions sont représentées sous forme de vecteurs-colonnes. Dans ce cas, chaque colonne x de la matrice X correspond à une acquisition, et on obtient la contribution de la composante p au vecteur Xx en effectuant le produit du terme vecteur Sp (pième colonne de la matrice S) par le coefficient ap,X. La factorisation de la matrice X en matrices non négatives permet alors de déterminer S et A telles que X = S A telles que : - X représente la matrice des acquisitions, de dimension (Nt, Nx), - S représente la matrice des profils des contributions, de dimension (Nt, Np), - A représente la matrice des poids, de dimension (Np, Nx). 44 Bien que dans cette description, on ait considéré deux composantes (une composante de fluorescente et une composante parasite), l'invention s'applique à un nombre plus important de composantes, par exemple plusieurs composantes de fluorescence et une composante parasite, ou seulement plusieurs composantes de fluorescence lorsque la composante parasite s'avère être négligeable. L'invention s'applique à d'autres cas que 10 celui de marqueurs fluorescents. En effet, l'injection de marqueurs fluorescents peut présenter certaines contraintes. Il est parfois préférable de s'en affranchir et de chercher à localiser, dans un milieu 15 120 (figure 7), non pas des marqueurs fluorescents, mais des zones 1221r 1222, 1223, 1224 présentant des propriétés optiques différentes de celles du milieu diffusant 120 environnant. En particulier, on peut rechercher des zones, dites zones d'intérêt, ces zones 20 d'intérêt étant alors considérées comme des marqueurs optiques et présentant : - un coefficient d'absorption et ou de diffusion différent de celui du milieu diffusant (dans ce cas on détecte un rayonnement diffusé à travers le 25 tissu, à la longueur d'onde d'excitation, - et/ou un coefficient de diffusion différent de celui du milieu diffusant environnant ; dans ce cas, on détecte un rayonnement diffusé différemment par les marqueurs et par le milieu 30 diffusant. 45 Le cas en de zones absorbantes, qui peuvent être considéré comme marqueurs, est celui de certaines tumeurs cancéreuses. Un procédé selon l'invention comprend alors des étapes identiques ou similaires à celui précédemment décrit dans le cadre de la fluorescence, à savoir : - l'éclairement d'un milieu diffusant par une lumière d'excitation en impulsions, de préférence dans l'infra-rouge ou dans le rouge, afin de limiter l'absorption dans les tissus, - la détection du signal diffusé résolu en temps par le milieu, à la longueur d'onde d'excitation, et la réalisation d'acquisitions prenant la forme d'histogrammes temporels XX du temps d'arrivée des photons, issus du milieu, sur le détecteur, - la constitution d'une matrice X rassemblant une ou plusieurs (NX) acquisition XX, - la décomposition de la matrice X en un produit de deux matrices non négatives A et S, - la détermination d'une composante d'intérêt, ou composante utile, en multipliant des termes de la matrice A par des termes de la matrice S, par exemple en multipliant un coefficient ax,p par les termes de la matrice S représentatifs de la contribution p. La composante utile IX à l'acquisition XX représente alors le signal diffusé par les zones d'intérêt à la longueur d'onde d'excitation Xex. On peut alors déterminer, par reconstruction, la position de ces zones dans le milieu. 46 Dans tous les cas un rayonnement d'excitation a une longueur d'onde inférieure à une ou plusieurs longueurs d'one détectées, quel que soit le type d'interaction entre le rayonnement d'excitation et le milieu, contenant un ou plusieurs marqueurs tels que défini ci-dessus : si le marqueur est un marqueur fluorescent, la longueur d'onde de fluorescence est supérieure à la longueur d'onde du rayonnement d'excitation, - si le marqueur a un coefficient d'absorption et/ou de diffusion différent de celui du milieu diffusant, alors le rayonnement détecté a une longueur d'onde supérieure ou égale au rayonnement d'excitation. REVENDICATIONS1. Procédé de localisation d'au moins un marqueur optique (22, 1221r 1222, 1223, 1224) dans un milieu diffusant, ledit marqueur ayant au moins une propriété optique différente du milieu diffusant, procédé dans lequel : a) on réalise au moins une acquisition d'un signal optique caractéristique du milieu (20, 120) et d'un ou plusieurs desdits marqueurs (22, 1221r 1222, 1223, 1224), par une interaction du milieu avec le rayonnement d'une source de rayonnement (8) en impulsions, chaque acquisition comportant d'une part une ou plusieurs composantes d'intérêt, due à un ou plusieurs marqueurs, d'autre part une composante parasite, due à une partie du milieu autre que les marqueurs, b) on forme, à partir des données de signal optique d'au moins une desdites acquisitions, un 20 tableau multidimensionnel X, c) on traite les données dudit tableau X par factorisation de ce tableau en un produit de, seulement, deux tableaux multidimensionnels non négatifs A et S, 25 d) on extrait la contribution, au signal optique, d'au moins une desdites composantes temporelles, à partir des données contenues dans les tableaux A et S. 30 48 2. Procédé selon la revendication 1, comportant en outre une étape de calcul de la position d'au moins un marqueur optique dans le milieu diffusant et/ou une représentation graphique de la répartition d'au moins un tel marqueur dans le milieu diffusant, à partir d'au moins une desdites composantes temporelles. 3. Procédé selon l'une des revendications précédentes, dans lequel le marqueur optique est un marqueur fluorescent, l'excitation du milieu étant alors réalisée à au moins une longueur d'onde d'excitation d'au moins un marqueur fluorescent, le signal optique résultant de l'excitation étant à au moins une longueur d'onde de fluorescence d'au moins un marqueur fluorescent, et dans lequel on extrait la contribution de fluorescence due à au moins un marqueur fluorescent à partir des termes des matrices A et S. 4. Procédé selon l'une des revendications 1 ou 2, dans lequel le marqueur optique est une zone présentant au moins un coefficient d'absorption et/ou de diffusion différent de celui du milieu diffusant environnant (20, 120), le signal optique résultant de la diffusion du signal d'excitation dans le milieu diffusant (20, 120) et dans lequel on extrait la contribution due à au moins une zone à coefficient d'absorption différent de celui du milieu diffusant. 5. Procédé selon l'une des revendications 1 à 4, dans lequel, avant formation dudit tableau multidimensionnel, on corrige les données de signaloptique de chaque acquisition par recentrage des données temporelles dans une plage temporelle par rapport à un temps de réponse impulsionnel. 6. Procédé selon l'une des revendications 1 à 5, dans lequel, lors de l'étape c), on détermine A et S par minimisation d'une fonction de coût. 7. Procédé selon la revendication 6, dans lequel, la fonction de coût est la distance OX- AS 2 entre l'image X et le produit A.S. 8. Procédé selon la revendication 6, dans lequel, la fonction de coût est la fonction : F- X-AS 2+a S-So 2 Où So est la matrice initiale et a un coefficient dit de régularisation. 9. Procédé selon l'une des revendications 1 à 8, dans lequel l'étape c) est réalisée par k itérations, des tableaux Ai+1 et Si+l, obtenus lors de l'itération d'ordre 1+1 étant déterminées à partir des tableaux Al et SI obtenus lors de l'itération d'ordre 1. 10. Procédé selon la revendication 9, dans lequel le nombre d'itération est déterminé en fonction des fluctuations des tableaux A et S, ou de façon automatique, en fonction de fluctuations de la fonction de coût au cours de 2 ou plusieurs itérations successives. 50 11. Procédé selon l'une des revendications 1 à 10, dans lequel le rayonnement d'excitation a un spectre au moins en partie dans l'infra rouge. 12. Procédé selon l'une des revendications 1 à 11, dans lequel le signal optique est détecté à au moins une longueur d'onde supérieure à 600 nm. 13. Dispositif d'identification d'au moins un marqueur optique dans un milieu diffusant, comportant : a) une source de rayonnement en impulsions (2, 4) pour produire un faisceau d'excitation du milieu et d'au moins un marqueur optique, b) des moyens (6, 8) pour réaliser au moins une acquisition d'un signal optique, caractéristique du milieu et d'un ou plusieurs marqueurs, chaque acquisition comportant d'une part une ou plusieurs composantes d'intérêt, due à un ou plusieurs marqueurs, d'autre part une composante parasite, due à une partie du milieu autre que les marqueurs, c) des moyens pour former, à partir des données de signal optique d'au moins une desdites acquisitions, un tableau multidimensionnel X, d) des moyens (24, 27, 28) pour traiter les données dudit tableau en un produit de, seulement, deux tableaux multidimensionnels non négatifs A et S, e) des moyens (24, 27, 28) pour extraire la contribution, au signal optique, d'au moins une desdites composantes temporelles, à partir des données contenues dans les tableaux A et S. 51 14. Dispositif selon la revendication 13, comportant en outre des moyens pour calculer la position d'au moins un marqueur optique dans le milieu diffusant et/ou une représentation graphique de la répartition d'au moins un tel marqueur dans le milieu diffusant, à partir d'au moins une desdites composantes temporelles. 15. Dispositif selon l'une des revendications 13 ou 14, comportant en outre des moyens, pour, avant formation dudit tableau multidimensionnel, corriger les données de signal optique de chaque acquisition par recentrage des données temporelles dans une plage temporelle autour d'un temps moyen du signal optique. 16. Dispositif selon l'une des revendications 13 à 15, comportant des moyens pour déterminer A et S par minimisation d'une fonction de coût. 17. Dispositif selon l'une des revendications 13 à 16, dans lequel la source de rayonnement en impulsions (2, 4) permet de produire un faisceau d'excitation au moins en partie dans l'infra rouge. 18. Dispositif selon l'une des revendications 13 à 17, dans lequel les moyens (6, 8) pour réaliser au moins une acquisition d'un signal optique comportant des moyens pour détecter un rayonnement à au moins une longueur d'onde supérieure à 600 nm.