FR 2984502 A1 20130621 FR 1161835 A 20111216 Procédé de reconstruction des propriétés optiques d'un milieu avec calcul d'un signal corrigé en fonction d'une première fonction de modélisation pour un milieu de référence et d'une deuxième distribution pour un milieu à caractériser, et système de reconstruction associé La présente invention concerne un procédé de reconstruction de propriétés optiques d'un milieu à l'aide d'un système de reconstruction comportant au moins une source de rayonnement propre à éclairer un milieu et au moins un détecteur propre à recevoir un signal émis par le milieu. L'invention concerne également un tel système de reconstruction. L'invention se situe notamment dans le domaine de l'imagerie médicale, et en particulier dans le domaine de l'imagerie de diffusion sur des tissus biologiques. L'invention s'applique à la détection de différentes pathologies, telles que des cancers du sein, de la prostate ou encore du tube digestif, le cerveau.... L'invention s'applique également à la dermatologie, notamment à la caractérisation d'une réaction d'inflammation cutanée. On connaît du document EP 2 302 362 Al un procédé et un système du type précité. Le procédé vise à localiser un ou plusieurs fluorophores contenus dans un milieu diffusant en éclairant le milieu diffusant par une source de rayonnement et en détectant, pour au moins un détecteur, le signal émis par le milieu à la longueur d'onde de fluorescence. Une distribution temporelle du signal reçu par le détecteur est alors réalisée pour chaque couple source-détecteur. Cette distribution temporelle dépend notamment des valeurs du temps de réponse moyen de la fonction source Ts et du temps de réponse moyen du détecteur -rd pour chaque couple source-détecteur, et il est alors nécessaire de connaître lesdites valeurs afin de localiser le ou les fluorophores. Les valeurs du temps moyen de la fonction source Ts et du temps de réponse moyen du détecteur -rd sont mesurées lors d'une opération de calibrage au cours de laquelle la source est placée face au détecteur en l'absence du milieu diffusant. Plus généralement, la distribution temporelle mesurée à l'aide du système de reconstruction correspond à la vraie distribution temporelle convoluée par une réponse instrument, également appelée IRF (de l'anglais Instrument Response Function). La réponse instrument représente la répartition temporelle de l'impulsion générée par la source de rayonnement et détectée par le détecteur en l'absence de milieu diffusant. Il est alors nécessaire de mesurer la réponse instrument afin de pouvoir déterminer ultérieurement la vraie distribution temporelle. Toutefois, la détermination de la réponse instrument pour chaque couple source-détecteur nécessite des mesures expérimentales longues et fastidieuses. Le but de l'invention est donc de proposer un procédé de reconstruction et un système de reconstruction permettant de reconstruire les propriétés optiques souhaitées du milieu tout en ne nécessitant pas une détermination expérimentale de la réponse instrument. A cet effet, l'invention a pour objet un procédé de reconstruction des propriétés optique d'un milieu, caractérisé en ce que: - l'élaboration, pour au moins un couple source-détecteur, d'une première distribution d'un signal reçu par le détecteur correspondant pour un milieu de référence, le signal reçu étant émis par le milieu de référence suite à l'éclairement dudit milieu par la source correspondante, - la détermination, pour ledit au moins un couple source-détecteur, d'une première fonction de modélisation d'un signal de diffusion de la lumière entre la source et le détecteur dans le milieu de référence, - l'élaboration, pour ledit au moins un couple source-détecteur, d'une deuxième distribution d'un signal reçu par le détecteur correspondant pour un milieu à caractériser, le signal reçu étant émis par le milieu à caractériser suite à l'éclairement dudit milieu par la source correspondante, et - le calcul d'un signal corrigé en fonction de la première fonction de modélisation et de la deuxième distribution. Par reconstruction des propriétés optiques, on entend par exemple : - la reconstruction de propriétés d'absorption, ces dernières étant notamment caractérisées par le coefficient d'absorption, noté pl,. - la reconstruction de propriétés de diffusion, ces dernières étant notamment caractérisées par le coefficient de diffusion réduit il', ou par le coefficient de diffusion D. - la reconstruction de propriétés de fluorescence, ces dernières étant notamment caractérisées par une fonction de réponse F d'un fluorophore, ou par une concentration c d'un fluorophore, ou encore par toute autre grandeur exprimant une quantité q d'un fluorophore, ce dernier étant par exemple endogène ou exogène. La reconstruction des propriétés optique permet de déterminer la distribution spatiale de coefficients optiques tels que ceux décrits ci-dessus, bien connus dans le domaine de l'imagerie optique de diffusion. Il s'agit d'une reconstruction en deux ou en trois dimensions. Suivant d'autres aspects avantageux de l'invention, le procédé de reconstruction comporte une ou plusieurs des caractéristiques suivantes, prises isolément ou suivant toutes les combinaisons techniquement possibles : - le calcul du signal corrigé est réalisé également en fonction de la première distribution, - le procédé comprend en outre la détermination, pour ledit au moins un couple source-détecteur, d'une deuxième fonction de modélisation d'un signal de diffusion de la lumière entre la source et le détecteur dans le milieu à caractériser, et le signal corrigé est calculé également en fonction de la deuxième fonction de modélisation, - l'étape de calcul comporte une opération de comparaison du produit de la première distribution et de la deuxième fonction de modélisation avec le produit de la deuxième distribution et de la première fonction de modélisation, - l'opération de comparaison comporte une opération arithmétique du produit de la première distribution et de la deuxième fonction de modélisation avec le produit de la deuxième distribution et de la première fonction de modélisation, - l'opération arithmétique est une soustraction du produit de la première distribution et de la deuxième fonction de modélisation du produit de la deuxième distribution et de la première fonction de modélisation, - l'opération arithmétique est un ratio entre le produit de la deuxième distribution temporelle et de la première fonction de modélisation et le produit de la première distribution temporelle et de la deuxième fonction de modélisation, - les première et deuxième distributions sont des distributions temporelles, et le produit de la première distribution et de la deuxième fonction de modélisation et le produit de la deuxième distribution et de la première fonction de modélisation sont des produits de convolution, - les première et deuxième distributions sont des distributions continues, et le produit de la première distribution et de la deuxième fonction de modélisation et le produit de la deuxième distribution et de la première fonction de modélisation sont des multiplications, - les première et deuxième distributions sont des distributions fréquentielles, et le produit de la première distribution et de la deuxième fonction de modélisation et le produit de la deuxième distribution et de la première fonction de modélisation sont des multiplications, - la première distribution est élaborée pour une longueur d'onde d'excitation, la deuxième distribution est élaborée pour une longueur d'onde d'émission, et le signal corrigé est exprimé en fonction du produit d'un scalaire positif et de la première distribution temporelle, le produit de la deuxième distribution et de la première fonction de modélisation étant fonction du produit du scalaire positif et de la première distribution temporelle, le scalaire positif étant égal au ratio entre une réponse instrument correspondant à la détection réalisée à la longueur d'onde d'émission et une réponse instrument correspondant à la détection réalisée à la longueur d'onde d'excitation, et - les propriétés optiques comprennent au moins un élément parmi le groupe consistant en : + les propriétés d'absorption de la lumière, caractérisées notamment par le coefficient d'absorption, + les propriétés de diffusion, caractérisées notamment par le coefficient de diffusion réduit ou par le coefficient de diffusion, et + les propriétés de fluorescence, caractérisées notamment par une fonction de réponse d'un fluorophore, ou par une concentration du fluorophore, ou encore par une grandeur fonction d'une quantité du fluorophore. L'invention a également pour objet un système de reconstruction comprenant: - au moins une source de rayonnement propre à éclairer un milieu, - au moins un détecteur propre à recevoir un signal émis par le milieu, caractérisé en ce qu'il comprend en outre : - des moyens d'élaboration, pour au moins un couple source-détecteur, d'une première distribution d'un signal reçu par le détecteur correspondant pour un milieu de référence, le signal reçu étant émis par le milieu de référence suite à l'éclairement dudit milieu par la source correspondante, - des premiers moyens de détermination, pour ledit au moins un couple source-détecteur, d'une première fonction de modélisation d'un signal de diffusion de la lumière entre la source et le détecteur dans le milieu de référence, - des moyens d'élaboration, pour ledit au moins un couple source-détecteur, d'une deuxième distribution d'un signal reçu par le détecteur correspondant pour un milieu à caractériser, le signal reçu étant émis par le milieu à caractériser suite à l'éclairement dudit milieu par la source correspondante, et - des moyens de calcul d'un signal corrigé en fonction de la première fonction de modélisation et de la deuxième distribution. Suivant d'autres aspects avantageux de l'invention, le système de reconstruction comprend une ou plusieurs des caractéristiques suivantes, prises isolément ou suivant toutes les combinaisons techniquement possibles : - les moyens de calcul sont propres à calculer le signal corrigé également en fonction de la première distribution, - le système comprend en outre des deuxièmes moyens de détermination, pour ledit au moins un couple source-détecteur, d'une deuxième fonction de modélisation d'un signal de diffusion de la lumière entre la source et le détecteur dans le milieu à caractériser, et le signal corrigé est calculé également en fonction de la deuxième fonction de modélisation. Ces caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple, et faite en référence aux dessins annexés sur lesquels : - la figure 1 est une représentation schématique d'un système de reconstruction selon l'invention, et - la figure 2 est un organigramme d'un procédé de reconstruction selon l'invention. Sur la figure 1, un système de reconstruction 10 est destiné à réaliser un examen d'un milieu 12 via l'acquisition d'une image du milieu 12, puis le traitement de cette image. Le système de reconstruction 10 comprend une pluralité de sources de rayonnement 14 aptes à éclairer le milieu 12 et une pluralité de détecteurs 16 propres à recevoir un signal émis par le milieu 12. Les sources de rayonnement 14 et les détecteurs 16, sont, par exemple, agencés dans une sonde, non représentée. Le système 10 comprend une unité de traitement d'informations 18 et un écran 20 d'affichage d'une image du milieu. Dans le mode de réalisation décrit, le système 10 est un système temporel, et chaque source de rayonnement 14 est une source de rayonnement pulsée. Chaque détecteur 16 est un détecteur résolu en temps, c'est-à-dire un détecteur permettant de mesurer la distribution temporelle du temps d'arrivée des photons. Le détecteur 16 est, par exemple, un dispositif de type TCSPC (de l'anglais Time-Correlated Single Photon Counting). Selon l'exemple représenté, le milieu 12 présente une surface d'observation 22, en forme d'un plan parallèle à un axe longitudinal X et un axe transversal Y, et au contact de laquelle les sources 14 et les détecteurs 16 sont propres à être appliqués. Le milieu 12 présente une direction d'observation 24 s'étendant selon un axe vertical Z et sensiblement perpendiculairement à la surface d'observation 22. Lors d'une première phase, le milieu 12 est un milieu de référence, noté MR et non représenté, également appelé fantôme. Le milieu de référence MR est connu en soi et est, par exemple, décrit dans le document FR 2 950 241 Al. Le fantôme MR a des caractéristiques optiques connues, de préférence proches à quelques dizaines de pourcents, par exemple à 50% près, voire à 30% près, de celles d'un milieu à caractériser, noté MC, que l'on souhaite examiner Le fantôme MR est, par exemple, un fantôme ayant des coefficients d'absorption pa° et de diffusion D° homogènes de valeurs médianes pour les tissus humains, par exemple égales à 0,2 cm-1, et respectivement à 10 cm-1. Lors d'une deuxième phase, le milieu 12 est le milieu à caractériser MC, visible sur la figure 1. Le milieu à caractériser MC comporte des tissus biologiques d'un animal ou de l'homme. Le milieu à caractériser MC est, par exemple, une zone d'un organe tel que le cerveau, un sein, la prostate, le tube digestif ou encore d'autres organes dans lesquels des fluorophores sont susceptibles d'être injectés. Plus généralement, le milieu à caractériser MC est un milieu biologique, par exemple un organe, dont on souhaite déterminer des propriétés optiques de diffusion et notamment la distribution spatiale de coefficients tels que le coefficient d'absorption ou le coefficient de diffusion. Dans l'exemple de réalisation de la figure 1, le milieu à caractériser MC est un milieu diffusant contenant des inclusions 26, dont les propriétés optiques d'absorption ou de diffusion sont différentes de celles du milieu de référence MR. Une seule inclusion 26 est représentée sur la figure 1 pour la clarté du dessin. Ainsi, en réalisant une cartographie du coefficient d'absorption pl, et du coefficient de diffusion D, on localise leurs hétérogénéités, ce qui permet de localiser les inclusions. Chaque source de rayonnement pulsée 14 comporte une source lumineuse pulsée 28 et une fibre optique d'excitation 30 reliée à la source pulsée 28 pour la transmission de l'impulsion lumineuse jusqu'au milieu 12. Lorsqu'une telle source fibrée est utilisée, on assimile l'extrémité libre de la fibre d'excitation 30 à la source s. En variante non représentée, chaque source de rayonnement pulsée 14 comporte une fibre optique d'excitation 30 reliée à une unique source lumineuse pulsée commune à la pluralité de sources de rayonnement 14. Selon cette variante, le système 10 comprend en outre un commutateur optique ou un multiplexeur pour sélectionner la fibre d'excitation dans laquelle le faisceau lumineux est envoyé. En variante encore, l'ensemble des sources de rayonnement pulsée 14 est formé d'une unique source lumineuse pulsée et d'un dispositif à miroir de type MEMS (de l'anglais MicroElectroMechanical Systems), non représenté, pour balayer le milieu 12 30 avec la lumière issue de la source lumineuse pulsée. En variante encore, la source 14 est une source multiplexée, délivrant des impulsions lumineuses à différentes longueur d'onde. Dans le cas d'un examen de la surface 22 du milieu ou d'un examen à une faible profondeur du milieu 12, c'est-à-dire une profondeur de quelques millimètres, la longueur d'ondes de chaque source de rayonnement pulsée 14 est, de préférence, dans le proche infrarouge ou le visible, c'est-à-dire comprise entre 500 nm et 1300 nm. Le taux de répétition est d'environ 50 MHz. Les impulsions émises par chaque source de rayonnement pulsée 14 présentent une durée comprise entre 500 picosecondes et 50 femtosecondes, chaque source lumineuse pulsée 28 étant propre à générer une impulsion de largeur temporelle de durée comprise entre 500 picosecondes et 50 femtosecondes. Chaque détecteur résolu en temps 16 comporte une fibre optique de détection 32 reliée à un module de détection résolu en temps 34. L'extrémité libre de cette fibre optique de détection 32 est assimilée au détecteur d. Dans l'exemple de réalisation de la figure 1, le module de détection 34 comporte un organe de détection 36 pour chaque détecteur 16. En variante non représentée, le module de détection 34 comporte un organe de détection commun à plusieurs détecteurs 16, en particulier un unique organe de détection pour l'ensemble des détecteurs 16. L'organe de détection 36 est, par exemple, un photomultiplicateur, ou une photodiode à avalanche, également appelée APD (de l'anglais Avalanche Photodiode), ou une photodiode à avalanche à photon unique, également SPAD (de l'anglais Single-Photon Avalanche Diode) ou encore un tube intensificateur d'images ayant une ou plusieurs anodes (de l'anglais Multi-Channel Plate). La sonde est une sonde compacte pour le diagnostic de certains cancers, tels qu'une sonde portative pour le diagnostic du cancer du sein, une sonde endorectale pour le diagnostic du cancer de la prostate, ou encore une sonde de dermatologie. En variante, la sonde est une sonde d'endoscope telle qu'une sonde souple pour le diagnostic du cancer digestif. L'unité de traitement de l'information 18 comporte un processeur de données 38 et une mémoire 40 associée au processeur 38. L'unité de traitement 18 comporte des moyens 42 d'élaboration, pour au moins un couple source 14 - détecteur 16, d'une première distribution Asd d'un signal reçu par le détecteur 16 correspondant pour le milieu de référence MR lors de la première phase, le signal reçu étant émis par le milieu de référence MR suite à l'éclairement dudit milieu MR par la source 14 correspondante, et d'une deuxième distribution Bsd d'un signal reçu par le détecteur 16 correspondant pour le milieu à caractériser MC lors de la deuxième phase, le signal reçu étant émis par le milieu à caractériser MC suite à l'éclairement dudit milieu MC par la source 14 correspondante. Dans le mode de réalisation décrit, les première et deuxième distributions sont des distributions temporelles respectivement notées Asd(t) et Bsd(t). Les moyens d'élaboration 42 sont, de préférence, réalisés sous forme d'une ou plusieurs cartes électroniques connectées au détecteur 16 et permettant de mesurer la distribution temporelle du temps d'arrivée des photons. Chaque source de lumière pulsée 28 comporte un laser impulsionnel. En variante, chaque source lumineuse pulsée 28 comporte une diode laser. En variante encore, chaque source lumineuse pulsée 28 comporte une source de lumière constante dont l'intensité lumineuse est modulée en des impulsions de durée équivalente par un dispositif à obturation rapide. Ainsi, selon ce premier mode de réalisation, la source est apte à émettre un rayonnement lumineux prenant la forme d'une impulsion temporelle. L'extrémité de chaque fibre d'excitation 30 s'étend perpendiculairement à la surface d'observation 22 c'est-à-dire selon l'axe Z, ou encore de manière oblique par rapport à l'axe Z afin d'émettre la lumière en biais par rapport à la surface d'observation 22. La mémoire 40 est apte à stocker un premier logiciel 44 de détermination, pour ledit au moins un couple source 14 - détecteur 16, d'une première fonction Gsd° de modélisation d'un signal de diffusion de la lumière entre la source 14 et le détecteur 16 dans le milieu de référence MR. La mémoire 40 est apte à stocker un logiciel 46 de calcul d'un signal corrigé en fonction de la première distribution, de la première fonction de modélisation et de la deuxième distribution. La mémoire 40 est apte à stocker un deuxième logiciel 48 de détermination, pour ledit au moins un couple source 14 - détecteur 16, d'une deuxième fonction Gsd de modélisation d'un signal de diffusion de la lumière entre la source 14 et le détecteur 16 dans le milieu à caractériser MC. Le logiciel de calcul 46 est alors propre à calculer le signal corrigé également en fonction de la deuxième fonction de modélisation. En variante, les premiers moyens de détermination 44, les moyens de calcul 46 et les deuxièmes moyens de détermination 48 sont réalisés sous forme de composants logiques programmables ou encore sous forme de circuits intégrés dédiés. Les première et deuxième fonctions de modélisation Gsd°, Gsd sont, par exemple, des fonctions de Green, bien connues dans le domaine de l'imagerie médicale. Chaque fonction de Green Gsd°, Gsd représente la densité de photons au niveau du détecteur 16 lorsque le milieu 12, c'est-à-dire le milieu de référence MR pour la première fonction de modélisation et respectivement le milieu à caractériser MC pour la deuxième fonction de modélisation, est éclairé par la source 14, où s et d sont des indices désignant respectivement la source 14 et le détecteur 16. La première distribution temporelle Asd(t) vérifie alors l'équation suivante : Asd (t) = IRF sd (t)* G:1 (t) (1) où IRFsd(t) représente la réponse instrument, c'est-à-dire l'influence de la source s et du détecteur d sur la première distribution élaborée. La deuxième distribution temporelle Bsd(t) vérifie l'équation suivante : Bsd (t) = IRF sd (t)* Gsd (t) (2) Dans le mode de réalisation décrit, les première et deuxième fonctions de modélisation et le signal corrigé sont respectivement notés Gsd°(t), Gsd(t) et Ysd(t). Le logiciel de calcul 46 est propre à effectuer une opération de comparaison du produit de la première distribution Asd et de la deuxième fonction de modélisation Gsd avec le produit de la deuxième distribution Bsd et de la première fonction de modélisation Gsd°. L'opération de comparaison comporte une opération arithmétique du produit de la première distribution Asd et de la deuxième fonction de modélisation Gsd avec le produit de la deuxième distribution Bsd et de la première fonction de modélisation Gsd°. L'opération arithmétique est, par exemple, une soustraction du produit de la première distribution Asd et de la deuxième fonction de modélisation Gsd du produit de la deuxième distribution Bsd et de la première fonction de modélisation Gsd°. Dans le mode de réalisation décrit, le produit de la première distribution Asd(t) et de la deuxième fonction de modélisation Gsd(t) et le produit de la deuxième distribution Bsd(t) et de la première fonction de modélisation Gsd°(t) sont des produits de convolution, les première Asd(t) et deuxième Bsd(t) distributions étant des distributions temporelles. Le signal corrig é Ysd(t) vérifie l'équation suivante : Ysd (t) - Bsd (t) * (t) Asd (t) * Gsd (t) (3) où Bsd(t) représente respective la deuxième distribution temporelle, _ Gsdom = Gso (rd,t) = G° (rs, rd, t) est une fonction de Green représentant la densité de photon à l'instant t, au niveau du détecteur situé en rd, lorsque le milieu de référence est éclairé par une source s située en rs., cette fonction étant appelée première fonction de modélisation, - Gsd (t) = Gs (rd,t) = G (rs, rd, t) est une fonction de Green représentant la densité de photon à l'instant t, au niveau du détecteur situé en rd, lorsque le milieu examiné est éclairé par une source s située en rs, cette fonction étant appelée deuxième fonction de modélisation, et - Asd(t) représente la première distribution temporelle. En variante, l'opération arithmétique est un ratio entre le produit de la deuxième distribution temporelle et de la première fonction de modélisation et le produit de la première distribution temporelle et de la deuxième fonction de modélisation. Selon cette variante, le signal corrig é Ysd(t) vérifie alors l'équation suivante : B d (t) * G(s)d (t) Y sd (t) = s (4) Asd (t) * G sd (t) Le fonctionnement du système de reconstruction 10 selon l'invention est décrit ci- après à l'aide de la figure 2 représentant un organigramme du procédé de traitement d'image selon l'invention. Initialement, lors de l'étape 100, le milieu de référence MR, également appelé fantôme, est éclairé par la ou les sources 14 correspondantes. En cas d'une pluralité de sources 14, l'éclairement du milieu 12 est effectué successivement pour chacune des sources 14. Lorsqu'il est éclairé, le milieu de référence MR émet en retour un signal lumineux à destination du ou de chaque détecteur 16, et la première distribution du signal reçu Asd(t) est élaborée pour chaque couple source 14 - détecteur 16, également noté (s, d), à l'aide des moyens d'élaboration 42. Lors de l'étape 110, le système de reconstruction 10 détermine ensuite, à l'aide des premiers moyens de détermination 44 et pour chaque couple source-détecteur (s, d), la première fonction Gsd°(t) de modélisation du signal de diffusion de la lumière entre la source 14 et le détecteur 16 dans le milieu de référence MR. Les étapes 100 et 110 correspondent à la première phase réalisée sur le milieu de référence MR. Un opérateur change alors le milieu 12, et remplace le milieu de référence MR par le milieu à caractériser MC, afin que le système de reconstruction 10 effectue les étapes suivantes du procédé de traitement d'image. Le milieu à caractériser MC est d'abord éclairé par la ou les sources 14 correspondantes (étape 120). Lorsqu'il est éclairé, le milieu à caractériser MC émet en retour un signal lumineux à destination du ou de chaque détecteur 16, et la deuxième distribution du signal reçu Bsd(t) est élaborée pour chaque couple source-détecteur (s, d), à l'aide des moyens d'élaboration 42. Lors de l'étape 130, le système de reconstruction 10 détermine ensuite, à l'aide des deuxièmes moyens de détermination 48 et pour chaque couple source-détecteur (s, d), la deuxième fonction Gsd(t) de modélisation du signal de diffusion de la lumière entre la source 14 et le détecteur 16 dans le milieu à caractériser MC. Le signal corrigé Ysd(t) est alors calculé, lors de l'étape 140 et à l'aide des moyens de calcul 46, en fonction de la première distribution Asd(t), de la première fonction de modélisation Gsd°(t), de la deuxième distribution 13,d(t) et de la deuxième fonction de modélisation Gsd(t). Le signal corrigé Ysd(t) vérifie, par exemple, l'équation (3) précédente. En complément, la deuxième fonction de modélisation Gsd(t) est déterminée de manière approchée à l'aide de la fonction de Green Gsdi(t) pour le milieu à caractériser. Le processus de reconstruction des propriétés optiques est généralement itératif, de telle sorte qu'au cours de chaque itération, on obtient des valeurs pal(r) et D1(r), auxquelles correspondent des fonctions de Green G si d (t) pour chaque couple source-détecteur (s, d). Le processus itératif vérifie alors les équations suivantes : Y sd (t) = B sd (t) * G (s)d (t) - Asd (t) * Gisd (t) IRF sd (t)*Gsd (t)* G(s)d (t) - IRF sd (t)* G (s)d (t)* Gisd (t) IRF sd (t)* G (s)d (t)* sd (t) Gisd (t)) (5) Asd (t) * (Gsd (t) Gisd (t)) - Asd (t)* (.1. G si ,t) * , t) * Gd (r, t)cli + r> Gis ,t) * , t) * Gid , t)cli) :--"" où pl, et D représentent respectivement le coefficient d'absorption et le coefficient de diffusion du milieu à caractériser MC, - IRFsd représente la réponse instrument correspondant au couple source-détecteur utilisé, - Gs(r,t) = G(rs,r,t) est la fonction de Green représentant la densité de photons à un endroit r du milieu MC lorsque le milieu MC est éclairé par la source s située en rs, - Gd(r,t) = G(rd,r,t) est la fonction de Green représentant la densité de photons à l'endroit r du milieu MC lorsque le milieu MC est éclairé par la source s située en rd. Cette fonction de Green s'écrit également Gd(r,t) = G(r,rd,t). Ainsi, cette fonction de Green représente également la densité de photons au niveau du détecteur (rd) lorsque le milieu est éclairé par une source située en r. Autrement dit, Gsd(t) correspond aux vraies valeurs pa(r, t) et D(r, t), c'est-à-dire aux valeurs recherchées, et Gsdl(t) correspond à des valeurs approchées pal(r, t) et D1(r, t), avec pal(r, t) = pa(r, t) + Eva(r, t) et D1(r, t) = D(r, t) + ÔD(r, t). L'objectif du procédé de reconstruction est alors de minimiser Eva(r, t) et ÔD(r, t), ce qui permet de faire tendre pl (r, t) et D1(r, t) respectivement vers pa(r, t) et D(r, t). Les coefficients pl, et D sont constants dans le temps, et les valeurs pa(r, t) et D(r,t) s'écrivent pa(r, t) = pa(r)8(t) et D(r,t) = D(r)8(t), avec am représentant une distribution de Dirac à l'instant t. Le calcul du signal corrigé Ysd(t), également appelé distribution temporelle corrigée, s'affranchit donc de la connaissance précise de la réponse instrument IRFsd(t) d'après les équations (3), (4) ou encore (5). Ceci permet d'éviter les mesures expérimentales longues et fastidieuses mises en oeuvre dans le procédé de l'état de la technique afin de déterminer la réponse instrument du système de reconstruction 10. Le signal corrigé Ysd(t) est ensuite mis en oeuvre dans des méthodes classiques de reconstruction du signal optique de diffusion. Ces méthodes mettent, par exemple, en oeuvre des transformées de la distribution temporelle Ysd(t), telles que des transformées de Mellin ou des transformées de Mellin-Laplace. Après une étape de maillage du milieu à caractériser MC en voxels, on obtient un système linéaire reliant : - les grandeurs mesurées, comprenant par exemple des transformées de Mellin ou de Mellin-Laplace de Ysd(t) - des grandeurs modélisées, à savoir les transformées des fonctions de Green, et - les inconnues, 4u a(r'' t) et bD(i-'' t) , l'indice m représentant le voxel de centre rm. Les transformées de la distribution corrigée Ysd, telles que les transformées de Mellin-Laplace Ysd(P'n), sont calculées pour différentes valeurs d'un ordre n comprises entre 0 et nmax, de préférence pour toutes les valeurs successives comprises entre 0 et nmax, suivant l'équation suivante : = (Bs(r1).Gs°,(/' - A( .G(r -1)) (6) + j=n Or, la distribution corrigée Ysd(t) vérifie l'équation suivante : Asd (t)* Gis * , t) * Gd (i- ,t)cli- + V>" Gis (i- ,t) * e(r-,t) * V>rGid 04(7)aYsd (t) = Les équations (6) et (7) permettent alors d'obtenir l'équation suivante, étant donné que la transformée de Mellin-Laplace d'un produit de convolution de deux termes est égale au produit des transformées de Mellin-Laplace de chacun de ces deux termes : YI") = - A(i:ik).(fGsi(p,i)(i)xypa(i-).G d1( P -1) + f V" G 51( P (i- ) (i- ) :Q" G dl(") (i- ) di-) (8) j+k=n L'équation (8) permet ainsi de calculer la transformée de Mellin-Laplace à l'ordre n de la distribution corrigée Ysd(P'n) des transformées de Mellin-Laplace de la première distribution temporelle Asd(t) et des fonctions de Green Gs1(r), G2(r) approchant la deuxième fonction de modélisation Gsd(t). Le milieu 12 est discrétisé en une pluralité M de voxels référencés m, avec m compris entre 1 et M, afin de calculer les transformées de Mellin-Laplace des fonctions de Green Gsl(r), Gd1(r). Les fonctions de Green Gsl(r), Gd1(r) discrétisées pour chaque voxel m sont notées Gsl(r,), Gdi (rm). La reconstruction de propriétés optiques du milieu 12 est alors effectuée, afin de déterminer par exemple les vecteurs pi et b> , dont les termes pa(m) et D(m), discrétisés pour chaque voxel m, c'est-à-dire µa = ) et D (i- ), constituent les cartes des propriétés optiques recherchées. L'étape de reconstruction vise alors à résoudre le système matriciel suivant : Y = WX , (9) où Y représente un vecteur des observations, W représente une matrice de passage, et X représente un vecteur des inconnues. Le vecteur des observations Y contient les Nn transformées de Mellin-Laplace de la distribution corrigée Ysd(P'n), avec Nn égal à n,'+1, calculées pour toutes les valeurs successives de l'ordre n comprises entre 0 et n,' et pour les couples source-détecteur (s, d) considérés. On définit, pour chaque triplet (s, d, n) des indices s, d et de l'ordre n, un indice I de la manière suivante : / = (s - 1) x Nd x Nn + (d -1)x Nn + n (10) où Nd est le nombre de détecteurs 16 considérés. La valeur maximale de l'indice I est égale à Imax, tel que : = Ns x Nd x Nn (11) Le vecteur des observations Y comporte alors Imax lignes. La matrice de passage W comprend une première partie WP comportant des premiers termes notés WP (/,m) et une deuxième partie W D comportant des deuxièmes termes notés W D (1 ,M) . Cette matrice est également appelée matrice de sensibilité des mesures Y aux propriétés optiques X. Les premiers termes WP (/,m) vérifient l'équation suivante : W" (I m) sl(p,i) ).Gd(p,j)(.m)Vm I+ j+k=n et les deuxièmes termes W D (1 ,M) vérifient l'équation suivante WD (/,m) =- ).'"Gd1(")(i-',).V', (13) i+ j+k=n OÙ V, représente un élément de volume entourant la maille m, tel que le volume de la cellule de Voronoi. Autrement dit, V, représente le volume du voxel m. Selon ce mode de réalisation, au cours de chaque itération, les transformées de Mellin Laplace, de différents ordres, de Gsl(t) et de Gd1(t) sont déterminées, notamment en utilisant l'expression suivante : ( 5(r) si n = 0 -V D(i)VGs1(")(i)+ ma (i)+ G si(p,n) (i.)P (14) Gl(p,n1) si n > 0 Ainsi, lors de chaque itération, on détermine Gs1("-0)puis on détermine les transformées d'ordre supérieur de façon itérative, en utilisant l'expression ci-dessus lorsque n > 0. Lors de la première itération, on utilise D = D° et pl, = Ita°, et lors des itérations suivantes, on (12) C c s utilise D = D1 et pl, = pal. Le même raisonnement s'applique naturellement pour la détermination des transformées de Mellin Laplace de la fonctionGd1(t). La matrice de passage W comporte alors Imax lignes et 2M colonnes. Chaque terme WP (/,m) , WD (/,m) de la matrice de passage W est déterminé par modélisation, à chaque itération, en fonction des propriétés optiques déterminées lors de l'itération précédente, ou lors de la première itération, en fonction de propriétés optiques initialisées par l'opérateur. Cette initialisation est réalisée en considérant par exemple un milieu 12 homogène. Le vecteur des inconnues X comporte 2M lignes et une colonne, et contient les inconnues pa(m) et D(m) pour chacun des M voxels. Les M premières lignes correspondent aux inconnues pa(m) et les M lignes suivantes du vecteur X correspondent aux inconnues D(m). L'inversion de la matrice de passage W est réalisée selon des algorithmes d'inversion bien connus de l'homme du métier, tels qu'un algorithme de descente de gradient, une technique de reconstruction algébrique (ART), une décomposition en valeurs singulières (SVD), ou encore une méthode des gradients conjugués. Le processus s'arrête lorsqu'un critère de convergence est atteint, par exemple lorsque la distance entre deux vecteurs X successifs est inférieure à un seuil prédéterminé. On conçoit ainsi que le procédé de reconstruction et le système de reconstruction 10 selon l'invention permettent de reconstruire les propriétés optiques souhaitées tout en ne nécessitant pas une détermination expérimentale de la réponse instrument. Selon une variante de ce mode de réalisation, on reconstruit la position de fluorophores dans le milieu à caractériser MC. Pour cela, on effectue une mesure Asd(t) correspondant à l'excitation du milieu de référence MR par une longueur d'onde, dite longueur d'onde de référence et notée Xref. Cette longueur d'onde de référence correspond à la longueur d'onde d'excitation du fluorophore Xex. Le signal Asd(t) est donc détecté à la longueur d'onde d'excitation. On fait alors l'hypothèse que la réponse du détecteur Dd(t) évolue de façon négligeable entre la longueur d'onde d'excitation Xex et la longueur d'onde de référence Xref, ou que cette évolution est maîtrisée, par exemple par une relation de proportionnalité. On présume alors que la réponse instrument IRF,,Âr (t) correspondant à la détection réalisée à la longueur d'onde d'excitation est reliée à la réponse instrument IRF: (t) correspondant à la détection réalisée à la longueur d'onde d'émission du fluorophore, par une relation de type /RF2/- (t) = alRF,7ï (t), a étant un scalaire positif. Cette mesure est réalisée sur le milieu de référence MR, dont les propriétés optiques sont connues. De préférence, ces dernières sont proches à quelques dizaines de pourcents, par exemple à 50% près, de celles du milieu à caractérisé MC. Asd (t) = IRF jex (t)* G(s)d (t) , (15) G s% (t) étant la fonction de Green correspondant à la densité de photon en i = id , dans le milieu de référence MR, lorsque ce dernier est éclairé par une source impulsionnelle de type Dirac en i = is Chaque mesure 13,d(t) correspond au signal de fluorescence, à la longueur d'onde d'émission X', mesuré par un détecteur d situé en rd, lorsque le milieu examiné est excité par une source d'excitation s, située en r = rs, et délivrant une impulsion lumineuse à la longueur d'onde d'excitation Xex- Chaque mesure 13,d(t) est modélisée de la façon suivante : B sd (t) = IRF sd2em (t)* f G s2- (i- , t) * F (i- , t) * G d2'n' (i- , t)di- (16) avec : - G54 (i- ,t) : fonction de Green correspondant à la densité de photons en i lorsque le milieu est éclairé par une source impulsionnelle Dirac située en is , - Gd2- (i,t): fonction de Green correspondant à la densité de photons en id lorsque le milieu est éclairé par une source impulsionnelle Dirac située en i , et - Mi, t): réponse du fluorophore localisé dans le volume di- , F (i- ,t) = F (i- ) exp(- -t) t désignant le temps de vie du fluorophore. F(r) est déterminé I" selon l'équation suivante : F (i- ) -,u,(i)17 ,17 étant le rendement de fluorescence. On établit alors la grandeur corrigée Ysd(t), telle que Ysd (t) = B sd (t) * G s% (t) = MFsdllem (t)* f G 54 (i- , t) * F (i- , t) * G d4 (i- , t)di- * G s% (t) 25(17) = crAsd (t) * I G 54 (i- , t) * F (i- , t) * Gd (i- , t)d i- On relie alors la grandeur corrigée Ysd(t) à des grandeurs mesurées ou connues (a, A sd (t)), à des grandeurs modélisées G 54 (i- ,t) , Gd- (i- ,t) , ainsi qu'aux inconnues F (i- , t) . En disposant de plusieurs grandeurs corrigées, correspondant à différents couples sources-détecteurs, on reconstruit les grandeurs F (i- ,t) selon des algorithmes d'inversion connus, tels que ceux mentionnés précédemment. Selon un deuxième mode de réalisation, pour lequel les éléments identiques au premier mode de réalisation, décrit précédemment, sont repérés par des références identiques, et ne sont pas décrits à nouveau, le système de reconstruction 10 est un système continu. Chaque source de rayonnement 14 est une source de rayonnement continue également appelée source de lumière continue, telle qu'une source blanche filtrée, ou encore telle qu'une source laser monochromatique. Chaque détecteur 16 est un détecteur continu propre à mesurer l'atténuation de l'intensité lumineuse. Chaque détecteur continu 16 est un capteur d'image, tel qu'un capteur CCD (de l'anglais Charge-Coupled Device) ou encore un capteur CMOS (de l'anglais Complementary Metal-Oxide Semiconductor). En variante, chaque détecteur continu 16 est un détecteur ponctuel, tel qu'un photomultiplicateur, une photodiode ou une photodiode à avalanche. Selon le deuxième mode de réalisation, les première et deuxième distributions sont des distributions continues respectivement notées Asd et Bsd. Les moyens d'élaboration 42 sont, de préférence, réalisés sous forme d'un logiciel apte à être stocké dans la mémoire 40. Selon le deuxième mode de réalisation, les première et deuxième fonctions de modélisation et le signal corrigé sont respectivement notés Gsd°, Gsd et Ysd. Selon le deuxième mode de réalisation, le produit de la première distribution Asd et de la deuxième fonction de modélisation Gsd et le produit de la deuxième distribution Bsd et de la première fonction de modélisation Gsd° sont des multiplications, les première Asd et deuxième Bsd distributions étant des distributions continues. L'opération arithmétique est, par exemple, une soustraction du produit de la première distribution continue Asd et de la deuxième fonction de modélisation Gsd du produit de la deuxième distribution continue Bsd et de la première fonction de modélisation Gsd o- Le signal corrigé Ysd vérifie l'équation suivante : Ysd - Bsd - Gsd Asd - Gsd (18) OÙ Bsd représente respective la deuxième distribution continue, Gsd° représente la première fonction de modélisation, Asd représente la première distribution continue, et Gsd représente la deuxième fonction de modélisation. En variante, l'opération arithmétique est un ratio entre le produit de la deuxième distribution continue Bsd et de la première fonction de modélisation Gsd° et le produit de la première distribution continue Asd et de la deuxième fonction de modélisation Gsd. Selon cette variante, le signal corrigé Ysd vérifie alors l'équation suivante : B - ° y sd sd Gsd (19) Asd - Gsd En complément, la deuxième fonction de modélisation Gsd est déterminée de manière approchée à l'aide de la fonction de Green Gsd' pour le milieu à caractériser, et on obtient alors les équations suivantes : Y sd B sd - G:d Asd - Gisd IRFsd - Gsd - Gsd - IRF sd - G(s)d - Gisd IRFsd - Gsd - Gsd Gsd) (20) Asd - Gsd Gisd) - Asd - Gis ) - (Sti ) - G id + f >Gis ) - 5D(r) - Gid(r)cli) où Gs(r) = G(rd,r) est la fonction de Green représentant la densité de photons à un endroit r du milieu MC lorsque le milieu MC est éclairé par la source s située en rs, et Gd(r) = G(rd,r) est la fonction de Green représentant la densité de photons à l'endroit r du milieu MC lorsque le milieu MC est éclairé par la source s située en rd. Autrement dit, Gsd correspond aux vraies valeurs pa(r) et D(r), et Gsd' correspond à des valeurs approchées pal(r) et Dl(r), avec pal(r) = pa(r) + 151.1,(r) et Dl(r) = D(r) + ÔD(r). Le fonctionnement de ce deuxième mode de réalisation est pour le reste analogue à celui du premier mode de réalisation décrit précédemment, et n'est pas décrit à nouveau. Les avantages de ce deuxième mode de réalisation sont analogues à ceux du premier mode de réalisation décrit précédemment, et ne sont pas décrits à nouveau. La reconstruction des propriétés optiques pa(r) ou D(r) est ensuite effectuée en supposant soit que Spa = 0, soit que 5D (r) = 0, selon des algorithmes d'inversion connus, tels que ceux mentionnés précédemment. Selon une variante de ce deuxième mode de réalisation, on reconstruit la position de fluorophores dans le milieu à caractériser MC. Pour cela, on effectue une mesure Asa correspondant à l'excitation du milieu de référence MR par la longueur d'onde de référence Xref. Cette longueur d'onde de référence correspond à la longueur d'onde d'excitation du fluorophore Xex. Le signal Asd est donc détecté à la longueur d'onde d'excitation Xex. On fait alors l'hypothèse que la réponse du détecteur Dd évolue de façon négligeable entre la longueur d'onde d'excitation Xex et la longueur d'onde de référence Xref, ou que cette évolution est maîtrisée, par exemple par une relation de proportionnalité. On présume alors que la réponse instrument /RF,7/ correspondant à la détection réalisée à la longueur d'onde d'excitation est reliée à la réponse instrument /RFs7/- correspondant à la détection réalisée à la longueur d'onde d'émission, par une relation de type /RF2/- = alRFs7i' , a étant un scalaire positif. Cette mesure est réalisée sur le milieu de référence MR, dont les propriétés optiques sont connues. De préférence, ces dernières sont proches à quelques dizaines de pourcents, par exemple à 50% près, de celles du milieu à caractériser MC. Asd IRFs2dex .Gs% , (21) G s% étant la fonction de Green correspondant à la densité de photon en r = rd , dans le milieu de référence MR, lorsque ce dernier est éclairé par une source ponctuelle en r = rs Chaque mesure Bsd correspond au signal de fluorescence, à la longueur d'onde d'émission Xem, mesuré par un détecteur d situé en rd, lorsque le milieu à caractériser MC est excité par une source d'excitation s, située en r = rs, et excitant le milieu à la longueur d'onde d'excitation Xex. Chaque mesure Bsd est modélisée de la façon suivante : Bsd = IRFs2dem .1 G52- (i).F (i).G (12- (i)di- (22) Avec : G 52- (i) : fonction de Green correspondant à la densité de photons en r lorsque le milieu est éclairé par une source ponctuelle située en is Gd4' (i): fonction de Green correspondant à la densité de photons en id lorsque le milieu est éclairé par une source ponctuelle située en F(i): réponse du fluorophore localisé dans le volume di- , cette fonction ayant été précédemment définie. On établit alors la grandeur corrigée Ysd, telle que IRFs2dem f G 52- (i- ).F (i).G (12- (i- )di- .G Ysd B sd *G s°d (23) sd .f G(i).F (i- ).G (12- ceA On relie alors la grandeur corrigée Ysd(t) à des grandeurs mesurées ou connues (a, Asd(t)), à des grandeurs modélisées G 54 (i) , Gd4' (i) , ainsi qu'aux inconnues F (i- ) En disposant de plusieurs grandeurs corrigées, correspondant à différents couples sources-détecteurs, on reconstruit les grandeurs F (i- ) selon des algorithmes d'inversion connus, tels que ceux mentionnés précédemment. Selon un troisième mode de réalisation, pour lequel les éléments identiques au premier mode de réalisation, décrit précédemment, sont repérés par des références identiques, et ne sont pas décrits à nouveau, le système de reconstruction 10 est un système comportant une source de rayonnement fréquentielle. Chaque source de rayonnement 14 est une source de rayonnement fréquentielle, c'est-à-dire une source modulée en intensité à une fréquence donnée. Selon cette variante, chaque source de rayonnement 14 comporte, par exemple, une source laser modulée en intensité à des fréquences radio, c'est-à-dire à des fréquences de l'ordre de quelques centaines de MHz. Chaque détecteur 16 est un détecteur fréquentiel et est propre à mesurer l'amplitude de l'intensité lumineuse et sa phase par rapport à la source 14. Selon le troisième mode de réalisation, les première et deuxième distributions sont des distributions fréquentielles complexes respectivement notées Asd(w) et Bsd(w). Les moyens d'élaboration 42 sont, de préférence, réalisés sous forme d'un logiciel apte à être stocké dans la mémoire 40. Selon le troisième mode de réalisation, les première et deuxième fonctions de modélisation et le signal corrigé sont respectivement notés Gsd°(W), Gsd(w) et Ysd(w). Ces trois fonctions sont complexes. Selon le troisième mode de réalisation, le produit de la première distribution Asd(w) et de la deuxième fonction de modélisation Gsd(w) et le produit de la deuxième distribution Bsd(w) et de la première fonction de modélisation Gsd°(w) sont des multiplications, les première Asd(w) et deuxième Bsd(w) distributions étant des distributions fréquentielles. L'opération arithmétique est, par exemple, une soustraction du produit de la première distribution fréquentielle Asd(w) et de la deuxième fonction de modélisation et de la première Gsd(w) du produit de la deuxième distribution fréquentielle Bsd(w) fonction de modélisation Gs2(W). Le signal corrigé Ysd(w) vérifie l'équation suivante : Ysd (co) - B sd (W) - G(s)d (W) - A sd (CO) - G sd (CO) OÙ Bsd(w) représente respective la deuxième distribution fréquentielle, Gsd°(w) représente la première fonction de modélisation, Asd(w) représente la première distribution fréquentielle, et Gsd(w) représente la deuxième fonction de modélisation. (24) Le fonctionnement de ce troisième mode de réalisation est pour le reste analogue à celui du premier mode de réalisation décrit précédemment, et n'est pas décrit à nouveau. Les avantages de ce troisième mode de réalisation sont analogues à ceux du premier mode de réalisation décrit précédemment, et ne sont pas décrits à nouveau. On conçoit ainsi que le procédé de reconstruction et le système de reconstruction selon l'invention permettent de reconstruire les propriétés optiques souhaitées tout en ne nécessitant pas une détermination expérimentale de la réponse instrument. REVENDICATIONS1.- Procédé de reconstruction de propriétés optiques d'un milieu, à l'aide d'un système de reconstruction (10) comportant au moins une source de rayonnement (14, s) propre à éclairer un milieu (12, MR, MC) et au moins un détecteur (16, d) propre à recevoir un signal émis par le milieu (12, MR, MC), le procédé étant caractérisé en ce qu'il comprend les étapes suivantes : - l'élaboration (100), pour au moins un couple source-détecteur (14, s, 16, d), d'une première distribution (Asd(t) ; Asa ; Asd(W)) d'un signal reçu par le détecteur (16, d) correspondant pour un milieu de référence (MR), le signal reçu étant émis par le milieu de référence (MR) suite à l'éclairement dudit milieu (MR) par la source (14, s) correspondante, - la détermination (110), pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une première fonction (Gsd°(t) ; Gsa° ; Gsd°(w)) de modélisation d'un signal de diffusion de la lumière entre la source (14, s) et le détecteur (16, d) dans le milieu de référence (MR), - l'élaboration (120), pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une deuxième distribution (Bsd(t) ; Bsd ; Bsd(W)) d'un signal reçu par le détecteur (16, d) correspondant pour un milieu à caractériser (MC), le signal reçu étant émis par le milieu à caractériser (MC) suite à l'éclairement dudit milieu (MC) par la source (14, s) correspondante, et - le calcul (140) d'un signal corrigé (Ysa(t) ; Ysa ; Ysd(w)) en fonction de la première fonction de modélisation (Gsd°(t) ; Gsa° ; Gsa°(w)) et de la deuxième distribution (Bsa(t) ; Bsd ; Bsa(w))- 2.- Procédé selon la revendication 1, dans lequel le calcul du signal corrigé est réalisé également en fonction de la première distribution (Asd(t) ; Asa ; Asa(w))- 3.- Procédé selon la revendication 1 ou 2, dans lequel le procédé comprend en outre la détermination (130), pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une deuxième fonction (Gsa(t)) de modélisation d'un signal de diffusion de la lumière entre la source (14, s) et le détecteur (16, d) dans le milieu à caractériser (MC), et dans lequel le signal corrigé (Ysd(t) ; Ysa ; Ysa(w)) est calculé également en fonction de la deuxième fonction de modélisation (Gsd(t) ; Gsd ; Gsa(w))- 4.- Procédé selon la revendication 3, dans lequel l'étape de calcul (140) comporte une opération de comparaison du produit de la première distribution (Asd(t) ; Asa ; Asd(W)) et de la deuxième fonction de modélisation (Gsd(t) ; Gsd ; Gsd(W)) avec le produit de la deuxième distribution (Bsd(t) ; Bsd ; Bsd(W)) et de la première fonction de modélisation (Gsd°(t) ; Gsa° ; Gsa°(w))- 5.- Procédé selon la revendication 4, dans lequel l'opération de comparaison comporte une opération arithmétique du produit de la première distribution (Asa(t) ; Asa ; Asa(w)) et de la deuxième fonction de modélisation (Gsd(t) ; Gsa ; Gsd(w)) avec le produit de la deuxième distribution (Bsd(t) ; Bsa ; Bsa(w)) et de la première fonction de modélisation (Gsd°(1) ; Gsd° ; Gsa°(w))- 6.- Procédé selon la revendication 5, dans lequel l'opération arithmétique est une soustraction du produit de la première distribution (Asa(t) ; Asa ; Asa(w)) et de la deuxième fonction de modélisation (Gsd(t) ; Gsa ; Gsd(W)) du produit de la deuxième distribution (Bsd(t) ; Bsa ; Bsd(W)) et de la première fonction de modélisation (Gsd°(t) ; Gsd° ; Gsa°(w))- 7.- Procédé selon la revendication 5, dans lequel l'opération arithmétique est un ratio entre le produit de la deuxième distribution temporelle (Bsd(t) ; Bsa ; Bsd(W)) et de la première fonction de modélisation (Gsd°(t) ; Gsa° ; Gsa°(w)) et le produit de la première distribution temporelle (Asa(t) ; Asa ; Asa(w)) et de la deuxième fonction de modélisation (Gsd(t) ; Gsa ; Gsa(w))- 8.- Procédé selon l'une quelconque des revendications 4 à 7, dans lequel les première et deuxième distributions sont des distributions temporelles (Asa(t), Bsd(t)), et dans lequel le produit de la première distribution (Asa(t)) et de la deuxième fonction de modélisation (Gsd(t)) et le produit de la deuxième distribution (Bsa(t)) et de la première fonction de modélisation (Gsa°(t)) sont des produits de convolution. 9.- Procédé selon l'une quelconque des revendications 4 à 7, dans lequel les première et deuxième distributions sont des distributions continues (Asa, Bsa), et dans lequel le produit de la première distribution (Asa) et de la deuxième fonction de modélisation (Gsa) et le produit de la deuxième distribution (Bsa) et de la première fonction de modélisation (Gsd°) sont des multiplications. 10.- Procédé selon l'une quelconque des revendications 4 à 7, dans lequel les première et deuxième distributions sont des distributions fréquentielles (Asa(w), Bsd(w)), et dans lequel le produit de la première distribution (Asa(w)) et de la deuxième fonction de modélisation (Gsd(w)) et le produit de la deuxième distribution (Bsd(w)) et de la première fonction de modélisation (Gsd°(w)) sont des multiplications. 11.- Procédé selon la revendication 1, dans lequel la première distribution (Asa(t) ; Asa ; Asa(w)) est élaborée pour une longueur d'onde d'excitation (Xex), la deuxième distribution (Bsd(t) ; Bsa ; Bsd(W)) est élaborée pour une longueur d'onde d'émission (X'), et dans lequel le signal corrigé (Ysa(t) ; Ysa ; Ysa(w)) est exprimé en fonction du produit d'un scalaire positif (a) et de la première distribution temporelle (Asa(t) ; Asa ; Asa(w)), le produit de la deuxième distribution (Bsd(t) ; Bsa ; Bsd(w)) et de la première fonction de modélisation (Gsd°(t) ; Gsd° ; Gsd°(w)) étant fonction du produit du scalaire positif (a) et dela première distribution temporelle (Asd(t) ; Asa ; Asd(W)), le scalaire positif (a) étant égal au ratio entre une réponse instrument correspondant à la détection réalisée à la longueur d'onde d'émission (IRFsdÀem(t)) et une réponse instrument correspondant à la détection réalisée à la longueur d'onde d'excitation (IRFsd"e((t)). 12.- Procédé selon l'une quelconque des revendications précédentes, dans lequel les propriétés optiques comprennent au moins un élément parmi le groupe consistant en : - les propriétés d'absorption de la lumière, caractérisées notamment par le coefficient d'absorption (14), - les propriétés de diffusion, caractérisées notamment par le coefficient de diffusion réduit (i.ts) ou par le coefficient de diffusion (D), et - les propriétés de fluorescence, caractérisées notamment par une fonction de réponse (F) d'un fluorophore, ou par une concentration (c) du fluorophore, ou encore par une grandeur fonction d'une quantité (q) du fluorophore. 13.- Système de reconstruction (10), comprenant : - au moins une source (14, s) de rayonnement propre à éclairer un milieu (12, MR, MC), - au moins un détecteur (16, d) propre à recevoir un signal émis par le milieu (12, MR, MC), caractérisé en ce qu'il comprend en outre : - des moyens (42) d'élaboration, pour au moins un couple source-détecteur (14, s, 16, d), d'une première distribution (Asd(t) ; Asd ; Asd(W)) d'un signal reçu par le détecteur (16, d) correspondant pour un milieu de référence (MR), le signal reçu étant émis par le milieu de référence (MR) suite à l'éclairement dudit milieu (MR) par la source (14, s) correspondante, - des premiers moyens (44) de détermination, pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une première fonction (Gsd°(t) ; Gsd° ; Gsd°(w)) de modélisation d'un signal de diffusion de la lumière entre la source (14, s) et le détecteur (16, d) dans le milieu de référence (MR), - des moyens (42) d'élaboration, pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une deuxième distribution (Bsd(t) ; Bsd ; Bsd(W)) d'un signal reçu par le détecteur (16, d) correspondant pour un milieu à caractériser (MC), le signal reçu étant émis par le milieu à caractériser (MC) suite à l'éclairement dudit milieu (MC) par la source (14, s) correspondante, et - des moyens (46) de calcul d'un signal corrigé (Ysd(t) ; Ysd ; Ysd(w)) en fonction de la première fonction de modélisation (Gsd°(t) ; Gsd° ; Gsd°(w)) et de la deuxième distribution (Bsd(t) ; Bsd ; Bsd(w))- 14.- Système (10) selon la revendication 13, dans lequel les moyens de calcul (46) sont propres à calculer le signal corrigé (Ysd(f) ; Ysa ; Ysd(w)) également en fonction de la première distribution (Asd(t) ; Asa ; Asa(w))- 15.- Système (10) selon la revendication 13 ou 14, dans lequel le système (10) comprend en outre des deuxièmes moyens (48) de détermination, pour ledit au moins un couple source-détecteur (14, s, 16, d), d'une deuxième fonction (Gsd(t) ; Gsd ; Gsd(W)) de modélisation d'un signal de diffusion de la lumière entre la source (14, s) et le détecteur (16, d) dans le milieu à caractériser (MC), et dans lequel le signal corrigé (Ysa(t) ; Ysa ; Ysa(w)) est calculé également en fonction de la deuxième fonction de modélisation (Gsd(t) ; Gsd ; Gsa(w))-