Brevet FR2890203A1 — publié le 20070302

Procédé de reconstruction de la distribution d'éléments fluorescents dans un milieu diffusant Domaine technique de l'invention L'invention concerne un procédé de reconstruction de la distribution d'éléments fluorescents dans un milieu diffusant ayant sensiblement une forme cylindrique finie, comportant une étape de formulation de fonctions de transfert de l'énergie io dans le milieu entre au moins une source de lumière d'excitation ponctuelle et au moins un détecteur. État de la technique En imagerie optique et en tomographie optique diffusive, l'objet examiné, ou au moins la zone examinée, présente souvent une forme cylindrique. Il peut s'agir d'un tube cylindrique dans lequel est disposée une souris, par exemple, ou une partie du corps humain. La tomographie optique diffusive classique comporte la reconstruction de contrastes d'absorption et/ou de diffusion, tandis que la tomographie optique diffusive de fluorescence comporte la reconstruction de la concentration en fluorophores et/ou du temps de vie de molécules fluorescentes. Dans ce dernier cas, on cherche à déterminer la distribution de marqueurs fluorescents dans un milieu diffusant à partir de données de mesure. On trouve classiquement deux approches pour le traitement des données de mesure. Dans une première approche, analytique, un cylindre infini est 30 considéré afin de faire une approximation bidimensionnelle, ce qui ne permet pas d'obtenir de résultats satisfaisants dans le cas où le cylindre n'est pas infini. Dans une seconde approche, qui est numérique et tridimensionnelle, un cylindre de taille finie est considéré. Dans ce cas, on utilise par exemple la méthode des éléments finis, la méthode des différences finies ou la méthode des intégrales de frontières. Cependant, ces méthodes sont coûteuses en temps de calcul. Dans l'article The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis de S.R. Arridge et al. (Phys. Med. Biol., 1992, Vol. 37, No 7, 1531-1560), la propagation de la lumière dans des tissus biologiques confinés dans un cylindre fini est donnée par analogie avec la conduction de la chaleur dans des solides, sans cependant tenir compte de la fluorescence. L'article de S.R.Arridge fait référence à l'ouvrage Conduction of heat in solids de H.S. Carslaw et al. (1986, Oxford at the Clarendon Press), qui décrit une solution de l'équation de conduction de la chaleur dans un cylindre fini par l'intermédiaire de fonctions de Green. L'expression est limitée aux points de la surface du cylindre. Objet de l'invention L'invention a pour but de remédier à ces inconvénients et, en particulier, de proposer un procédé de reconstruction de la distribution d'éléments fluorescents dans un milieu diffusant ayant sensiblement la forme d'un cylindre fini, ce procédé permettant d'utiliser une approche analytique et tridimensionnelle. Selon l'invention, ce but est atteint par le fait que l'étape de formulation comporte la formulation d'une pluralité de premières fonctions de transfert d'énergie représentatives respectivement du transfert de l'énergie entre la source de lumière d'excitation ponctuelle et les éléments fluorescents et la formulation d'une pluralité de secondes fonctions de transfert d'énergie représentatives respectivement du transfert de l'énergie entre les éléments fluorescents et le détecteur. Selon un développement de l'invention, les premières et secondes fonctions de transfert sont des fonctions de Green résolvant l'équation de diffusion et correspondant à un volume cylindrique fini, les fonctions de Green étant exprimées en fonction des fonctions de Bessel modifiées. io Description sommaire des dessins D'autres avantages et caractéristiques ressortiront plus clairement de la description qui va suivre de modes particuliers de réalisation de l'invention donnés à titre d'exemples non limitatifs et représentés aux dessins annexés, dans lesquels: La figure 1 illustre la propagation de la lumière dans un milieu diffusant ayant sensiblement une forme cylindrique finie. La figure 2 représente l'organigramme d'un mode de réalisation particulier du 20 procédé selon l'invention. Description de modes particuliers de réalisation Sur la figure 1, une source S de lumière d'excitation ponctuelle est placée à un point rs et émet une lumière ayant une première longueur d'onde XSF et une amplitude Q. La lumière émise se propage dans un volume V d'un milieu diffusant ayant sensiblement une forme cylindrique finie, dans lequel sont répartis des éléments fluorescents F selon une distribution que l'on cherche à déterminer. Sur la figure 1, une première onde diffusive LSF émise par la source S excite un élément fluorescent FI qui émet ensuite un rayonnement ayant une seconde longueur d'onde ?FD: l'intensité de cette seconde onde diffusive LFD est mesurée par le détecteur D. La distribution des éléments fluorescents doit alors être déterminée à partir des mesures effectuées par le détecteur D. Ces mesures sont effectuées à partir d'un ensemble de positions de source. Dans le cas de plusieurs détecteurs, les mesures effectuées par tous les détecteurs seront prises en compte. Sur la figure 1, deux autres éléments fluorescents F2 et F3 sont représentés. Les éléments fluorescents FI, F2 et F3 étant répartis dans le volume V, le détecteur io reçoit une densité de photons (/)M mesurée composée par toutes les secondes ondes LFD émises par l'ensemble des éléments fluorescents FI, F2 et F3. Pour des raisons de clarté, les ondes correspondant aux éléments fluorescents F2 et F3 ne sont pas représentées. On définit une première fonction de transfert d'énergie G(ÂsF,rs,FF) représentative du transfert de l'énergie entre la source S de lumière d'excitation ponctuelle et les éléments fluorescents F, et une seconde fonction de transfert d'énergie G(AFD,rF,FD) représentative du transfert de l'énergie entre les éléments fluorescents F et le détecteur D. On peut montrer que la densité théorique de photons OM mesurée par le détecteur D peut être exprimée, de façon approximative, par une intégrale comportant le produit des /première et seconde fonctions derr transfert G: 0M(rs,rD) = Q(rs) f G(xsF+rs,FF)N(rF)G(xFD,rF,rD)u,F (1), v où le paramètre S lrF)r (rF) __ 1-icvt(rF) (2) dépend du rendement quantique n de l'élément fluorescent, de l'absorption locale 8u(rF) due aux éléments fluorescents et de l'amortissement (1-icor(rF))"1 lié au temps de vie des éléments fluorescents F et à la pulsation w imposée extérieurement au rayonnement, la lumière étant émise de façon pulsée. Les première et seconde fonctions de transfert G peuvent ainsi être traitées comme des fonctions de Green déterminées en résolvant des équations de diffusion pour le volume cylindrique V fini et pour une longueur d'onde A donnée: V2G(A,r,ro)+ kxG(A,r,r'o) = 8(r - r0) (3), S étant la fonction de Dirac et k2 = 't "-+ tw le nombre d'onde. Dans ces Dx cfDA io expressions, en est la vitesse de la lumière dans le milieu, , est le coefficient d'absorption du milieu, DA est le coefficient de diffusion et r et ro sont les variables spatiales de la fonction de Green. Les paramètres a et D,, sont évalués aux longueurs d'onde correspondantes ? SF et 2 FD. De plus, les solutions des équations de diffusion doivent satisfaire à des conditions aux limites sur la surface délimitant le volume cylindrique, par exemple les conditions de Dirichlet ou les conditions de Neumann. Le problème mathématique est ainsi analogue au problème de la conduction de la chaleur dans un milieu condensé présentant une forme cylindrique finie. Pour un volume cylindrique fini, les fonctions de Green résolvant les équations de diffusion et satisfaisant aux conditions aux limites peuvent être exprimées en fonction des fonctions de Bessel modifiées: i 'ç p rz p7tzo +c. I,,(rkp)Fn(akp,rokp) G(,r,ro)= sin sin jJ cos n(9 0o) si rsro, Dl p= l l n=_oo I,,(akp) 1sin pr sin p o In(rokp)Fn(akp,rkp) = cosn(0 00) D,,7rl p_1 l l n__oo In (akp) si r > ro Dans ces équations (4), 1 est la hauteur du cylindre, a est le rayon du cylindre, r = (r,6,z) et ro = (ro, 60,zo) sont les variables spatiales de la fonction de Green en coordonnées cylindriques, F,,(u,v) = In(u)K,,(v)-Kn(u)In(v), In et Kn sont les fonctions de Bessel modifiées respectivement de première et de deuxième espèce d'ordre n. De plus, on définit kr =.k,2 + p27r2 /12. Les équations (4) sont ensuite introduites dans l'équation (1). Pour un calcul discret, l'espace du volume cylindrique V est divisé en N volumes élémentaires selon un maillage approprié. L'intégrale de l'équation (1) est alors remplacée par 1 o une somme sur les N volumes élémentaires dv j constituant le volume cylindrique V: N (/ OM(rS+rD) _.(rS)IG/ IASF,rS,rj) Nlrj)G(AFDrj+rD)dvj (5). j=1 Selon l'invention, le procédé de reconstruction comporte la formulation d'une 15 pluralité de N premières fonctions de transfert d'énergie G(AsF,rs,Fj) représentatives respectivement du transfert de l'énergie entre la source S de lumière d'excitation ponctuelle et les éléments fluorescents F et la formulation d'une pluralité de N secondes fonctions de transfert d'énergie G(2FD,rj, FD) représentatives respectivement du transfert de l'énergie entre les éléments 20 fluorescents et le détecteur. Ainsi, on associe à chaque volume élémentaire dv une première fonction de transfert d'énergie G(Â.SF, rs,rj) et une seconde fonction de transfert d'énergie G(ÂFD,rj,FD). Lorsqu'on considère une pluralité de Ns sources S de lumière d'excitation 25 ponctuelle et une pluralité de ND détecteurs D, une représentation matricielle peut être choisie: 1OM _ NDxNs r I GFDl NDxNNxNLGSF]NxNs LQ] Ns (6). Dans cette équation, chaque colonne de [gM]NoxNs représente la mesure sur ND détecteurs pour une source S donnée. Pour chaque source S et pour chaque détecteur D, une matrice de conversion J 5 peut être calculée, ce qui permet de relier de manière linéaire la densité de photons OM mesurée au paramètre (3: N YSD E JSDjN j j=1 L'ensemble des combinaisons source-détecteur permet de construire l'équation matricielle qui est ensuite résolue dans un algorithme de reconstruction, soit en calculant l'erreur entre les mesures expérimentales et cette équation théorique matricielle (par exemple par l'intermédiaire de la méthode de rétroprojection d'erreurs ou d'algorithmes de type ART (Algebraic Reconstruction Technique)), soit en inversant directement la matrice J (par exemple par l'intermédiaire d'algorithmes de type SVD (Single Value Decomposition)). Dans le cas du calcul d'erreur entre la mesure expérimentale et la théorie, on utilise une distribution théorique des éléments fluorescents qui est ajustée au fur et à mesure de manière à réduire cette erreur. Dans le cas de l'inversion de la matrice, on obtient, dans une première étape, le paramètre (3 qui dépend de la distribution des éléments fluorescents par l'intermédiaire de l'absorption locale 8 (FFF) due aux éléments fluorescents et par l'intermédiaire de l'amortissement (1icwr(rF))-' lié au temps de vie i des éléments fluorescents F. La connaissance du paramètre (3 permet ainsi de déterminer la distribution et la concentration locale des éléments fluorescents. Comme représenté à la figure 2, les paramètres et les variables sont définies (fonction F1 sur la figure 2) en début du procédé de reconstruction. On définit (7). ainsi la géométrie du cylindre (la hauteur l et le rayon a), les positions des sources (rs) et détecteurs (7D), le maillage du milieu et les paramètres constitutifs du milieu comme les coefficients de diffusion (DA) et d'absorption (u) et les longueurs d'ondes (XSF, 2FD). Ensuite les fonctions de Green G sont déterminées (fonction F2 sur la figure 2) conformément aux équations (4). La matrice de conversion J peut alors être déterminée. D'autre part, la mesure expérimentale (fonction F3 sur la figure 2) et un traitement des données expérimentales (fonction F4 sur la figure 2) permettent de déterminer les densités de photons 4)M mesurées par les différents détecteurs D. Les paramètres liés aux éléments fluorescents F sont ensuite reconstruits (fonction F5 sur la figure 2). Un traitement supplémentaire (fonction F6 sur la figure 2) permet ensuite de représenter sous forme d'image la distribution des éléments fluorescents F ainsi que de leur concentration. Le calcul des fonctions de Green conformément aux équations (4) peut être optimisé en temps de calcul en utilisant une représentation synthétique des fonctions de Green: G = Eapup (8) p In (rkp)Fn (akp, r kp) cos n(0 B ) . II n= co In (akp) où ap = sin p' sin pO up = Egn p et gn p = est plus judicieux de calculer d'abord, et une fois pour toutes, les fonctions ap (ne dépendant que de z) et up (ne dépendant que de r et 0) séparément, et d'effectuer ensuite leur produit scalaire. Le choix de l'ordre de troncature pour les indices n et p peut être effectué conformément au théorème de Shannon. 25 Pour n, 100 est un ordre suffisant et pour p, 20 est un ordre suffisant. Revendications 1. Procédé de reconstruction de la distribution d'éléments fluorescents (F) dans un milieu diffusant (V) ayant sensiblement une forme cylindrique finie, comportant une étape de formulation de fonctions de transfert de l'énergie dans le milieu entre au moins une source (S) de lumière d'excitation ponctuelle et au moins un détecteur (D), caractérisé en ce que l'étape de formulation comporte la formulation d'une pluralité de premières fonctions (G(A,sF,FFs,rJ)) de transfert d'énergie représentatives respectivement du transfert de l'énergie entre la source (S) de lumière d'excitation ponctuelle et les éléments fluorescents (F) et la formulation d'une pluralité de secondes fonctions (G(ÂFD,ri,FD)) de transfert d'énergie représentatives respectivement du transfert de l'énergie entre les éléments fluorescents (F) et le détecteur (D). 2. Procédé selon la revendication 1, caractérisé en ce que les premières et secondes fonctions de transfert (G) sont des fonctions de Green résolvant l'équation de diffusion et correspondant à un volume cylindrique fini, les fonctions de Green (G) étant exprimées en fonction des fonctions de Bessel modifiées (Fn, In, Kn).