Brevet FR3118408A1 — publié le 20220701

L’invention est un procédé de détection d’une occlusion artérielle ou veineuse à partir de mesures de photons rétrodiffusés, les mesures étant résolues en temps. Les mesures permettent d’obtenir une distribution temporelle, correspondant à un nombre de photons rétrodiffusés mesurés en fonction du temps. Différentes intensités temporelles sont déterminées à partir de la distribution temporelle, à chaque intensité temporelle correspondant une profondeur dans l’échantillon. A partir des intensités temporelles, le procédé comporte une détermination d’une occlusion artérielle ou veineuse. Figure 7B. Procédé d’analyse d’un échantillon par mesure résolue en temps de l’intensité de photons rétrodiffusés Le domaine technique de l’invention est l’application de mesures optiques au suivi de la vascularisation d’un lambeau sous-cutané, ayant été enfoui sous la peau d’une personne dans le cadre d’une chirurgie réparatrice. ART ANTERIEUR Dans le domaine de la chirurgie plastique reconstructrice, lorsqu’une perte de substance cutanée n’est pas suturable, ou ne peut pas cicatriser, ou nécessite un apport de matière, une technique chirurgicale possible est d’utiliser un lambeau pédiculé. Un tel lambeau est un fragment de peau et de tissu cellulaire sous-cutané, prélevé sur une autre partie du corps, et présentant une vascularisation autonome, passant par un pédicule. Au cours de la chirurgie de reconstruction et dans le cas d’une mastectomie par exemple, le lambeau est enfoui sous la peau d’origine pour reconstruire le volume du sein. Le pédicule est alors raccordé à une veine et une artère, de façon à permettre une irrigation du lambeau. Ce type d’intervention suppose un suivi régulier de l’oxygénation du lambeau en post opératoire. Une vascularisation défaillante du lambeau enfoui, par exemple lors d’une thrombose, peut en effet avoir des conséquences graves, par exemple la survenue d’une nécrose suite à un manque d’oxygénation. Une telle situation nécessite généralement une intervention chirurgicale dans les plus brefs délais. Cependant la surveillance de la bonne oxygénation du lambeau est difficile du fait de sa non accessibilité. On sait que les mesures optiques peuvent permettre, dans certaines conditions, un suivi de l’oxygénation de tissus superficiels. La spectroscopie de réflectance diffuse, usuellement désignée par l’acronyme DRS (Diffuse Reflectance Spectroscopy), est une technique d’analyse non destructive permettant d’estimer des propriétés de propagation de la lumière dans un milieu analysé. Cette technique est par exemple décrite dans EP2762064 ou EP3054282 ou EP3054281 ou EP3311138. Elle consiste à illuminer le milieu par un faisceau de lumière incident, et à détecter des photons rétrodiffusés par le milieu analysé, à distance du faisceau incident. Il est fréquent que la détection soit effectuée à différentes longueurs d’onde et ou à plusieurs distances du faisceau incident, de façon à obtenir des propriétés spectrales de propagation de la lumière dans le milieu analysé. Les propriétés de propagation de la lumière comprennent généralement des propriétés d’absorption et/ou des propriétés de diffusion de la lumière. Il s’agit notamment de coefficients d’absorption ou de diffusion, ces derniers représentant respectivement des probabilités d’absorption et de diffusion d’un photon par unité de longueur. L’estimation des propriétés de propagation de la lumière, à certaines longueurs d’onde, permet une estimation d’une concentration d’analytes dans le milieu analysé. Ainsi, la DRS peut être utilisée pour estimer un taux d’oxygénation de tissus superficiels, par une comparaison entre une quantité d’oxyhémoglobine et une quantité de hémoglobine totale (oxyhémogobline + désoxyhémoglobine). La spectroscopie de réflectance résolue en temps (TRS : Time-Resolved Reflectance Spectroscopy) prend en compte un aspect temporel de la mesure. Le milieu analysé est illuminé par un faisceau lumineux impulsionnel, de faible durée, typiquement comprise entre 10 et 100 ps. On détecte, à distance de ce dernier, les photons rétrodiffusés par le milieu, et on en détermine une distribution temporelle des photons détectés. La distribution temporelle permet de discriminer les photons détectés entre des photons prompts, correspondant aux premiers instants de la distribution temporelle, et des photons retardés, correspondant aux derniers instants de la distribution temporelle. Lors de l’interprétation des mesures, on prend en compte le fait que plus un photon est retardé, plus il est probable qu’il ait traversé des couches profondes du milieu. Une des applications de la TRS est l’analyse de milieux dans lesquels les propriétés de propagation de la lumière sont susceptibles de varier selon la profondeur, et en particulier entre des couches superficielles et des couches profondes. Des exemples d’application sont donnés dans les publications suivantes : Selb. J et al. «Improved sensitivity to cerebral hemodynamics during brain activation with a time-gated optical system: Analytical model and experimental validation », Journal of Biomerical Optics, 10(1), 2007. Contini, D et al. "Novel method for depth-resolved brain functional imaging by time-domain NIRS", Proc. SPIE 6629, Diffuse Optical Imaging of Tissue, 662908. 12 July 2007. Zucchelli L, Contini D, Re R, Torricelli A, Spinelli L. « Method for the discrimination of superficial and deep absorption variations by time domain fNIRS. » Biomed Opt Express . 2013;4(12):2893‐2910. Published 2013 Nov 20. doi:10.1364/BOE.4.002893. Ces publications décrivent des exemples d’application de la TRS à des milieux diffusants formés de différentes couches superposées, présentant des propriétés de propagation de la lumière différentes les unes des autres. Lorsque le milieu n’est pas homogène, et comporte plusieurs couches superposées, les photons retardés traversent chaque couche superposée. Le signal résultant des photons retardés dépend à la fois des propriétés optiques des couches superficielles et des couches profondes. Afin de caractériser correctement les couches profondes, la contribution des couches superficielles au signal détecté doit être prise en compte. Les inventeurs proposent une méthode, basée sur la TRS, permettant un suivi post-opératoire de la vascularisation d’un lambeau. La méthode est non invasive. Elle peut permettre une estimation de la concentration de l’hémoglobine ou de l’oxyhémoglobine dans le lambeau. Un premier objet de l’invention est un procédé de détection d’une occlusion veineuse ou artérielle dans un tissu biologique entre au moins un instant de mesure et un instant de référence, le tissu biologique étant délimité par une surface, le tissu biologique comportant une couche superficielle et une couche profonde, la couche superficielle s’étendant entre la surface du tissu et la couche profonde, le procédé comportant : a) illumination du tissu biologique par une source de lumière impulsionnelle, la source de lumière étant disposée en regard de la surface du tissu biologique, la source de lumière émettant un faisceau d’illumination formant une zone d’illumination à la surface du tissu biologique, le faisceau d’illumination étant émis simultanément ou successivement dans deux longueurs d’onde d’illumination différentes l’une de l’autre, chaque longueur d’onde s’étendant dans une bande spectrale d’absorption de l’oxyhémoglobine et/ou de la désoxyhémoglobine ; b) à chaque longueur d’onde, détection de photons rétrodiffusés par le tissu biologique, après s’être propagés dans le tissu biologique, par un photodétecteur, les photons rétrodiffusés détectés émanant d’une zone de détection à la surface du tissu biologique ; c ) à chaque longueur d’onde, obtention d’une distribution temporelle, représentant un nombre de photons détectés par le photodétecteur en fonction du temps, la distribution temporelle s’étendant selon une durée, à partir d’un instant initial ; les étapes a) à c) étant effectuées à l’instant de référence et à l’instant de mesure ou chaque instant de mesure ; le procédé étant caractérisé en ce qu’il comporte, à l’instant de mesure ou chaque instant de mesure et à l’instant de référence : d) segmentation de la durée de chaque distribution temporelle en intervalles de détection, à chaque intervalle de détection étant associé un rang, le rang étant d’autant plus élevé que l’intervalle de détection est éloigné de l’instant initial de la distribution temporelle, chaque intervalle de détection étant associé à une couche élémentaire s’étendant à une profondeur dans le tissu biologique ; e) pour chaque couche élémentaire, à partir de chaque distribution temporelle résultant de c), établies à l’instant de mesure et à l’instant de référence, détermination d’intensités temporelles, correspondant à un nombre de photons détectés durant chaque intervalle de détection; et en ce qu’il comporte, à l’instant de mesure ou à chaque instant de mesure f) pour chaque intervalle de détection, calcul d’un ratio entre l’intensité temporelle à l’instant de mesure et l’intensité temporelle à l’instant de référence ; g) à partir des ratios résultant de f), établis pour chaque longueur d’onde, calcul d’une grandeur d’intérêt pour chaque couche élémentaire respectivement associée à chaque intervalle de détection, la grandeur d’intérêt étant représentative d’une variation de la concentration d’hémoglobine, dans ladite couche élémentaire, entre l’instant de mesure et l’instant de référence ; h) pour chaque couche élémentaire, comparaison de la grandeur d’intérêt, calculée lors de g), avec une grandeur d’intérêt de base, calculée pour une couche élémentaire de base, la couche élémentaire de base étant une des couches élémentaires du tissu biologique, de façon à obtenir, pour chaque couche élémentaire, une grandeur d’intérêt différentielle ; i) en fonction d’une évolution, selon la profondeur du tissu biologique, des grandeurs d’intérêt différentielles, détermination de la survenue d’une occlusion veineuse ou artérielle dans la couche profonde. Avantageusement les étapes a) à h) sont mises en œuvre à différents instants de mesure ; l’étape i) comporte une formation d’un profil, pour chaque couche de élémentaire, chaque profil représentant une variation de la grandeur d’intérêt différentielle, dans ladite couche élémentaire, en fonction de l’instant de mesure ; l’étape i) comporte un calcul d’un critère de corrélation entre différents profils, correspondant respectivement à différentes couches élémentaires, le critère de corrélation traduisant une corrélation entre les différents profils; de telle sorte que la survenue de l’occlusion veineuse ou artérielle est détectée lorsque de part et d’autre d’une couche élémentaire, le critère de corrélation franchit une valeur seuil prédéterminée. Le procédé peut comporter, suite à i) : j) détermination d’une couche élémentaire d’interface, choisie parmi les différentes couches élémentaires, la couche élémentaire d’interface étant une couche de part et d’autre de laquelle le critère de corrélation franchit la valeur seuil, la couche élémentaire d’interface étant considérée comme formant une interface entre la couche superficielle et la couche profonde, la couche élémentaire d’interface étant associée à un intervalle temporel d’interface ; k) à l’instant de référence et à un instant de mesure, à partir de chaque distribution temporelle résultant de c) établie à l’instant de mesure, calcul : d’une première intensité temporelle, comportant une contribution, à la distribution temporelle, de photons détectés antérieurement à l’intervalle temporel d’interface, lesdits photons étant considérés comme s’étant propagés uniquement dans la couche superficielle ; d’une deuxième intensité temporelle ; comportant une contribution à la distribution de photons détectés postérieurement à l’intervalle temporel d’interface, lesdits photons étant considérés comme s’étant propagés dans la couche superficielle et dans la couche profonde l) calcul d’un premier ratio, à partir des premières intensités temporelles respectivement déterminées à l’instant de référence et à l’instant de mesure et calcul d’un deuxième ratio, à partir des deuxièmes intensités temporelles respectivement déterminées à l’instant de référence et à l’instant de mesure ; m) estimation des variations des concentrations d’oxyhémoglobine et de désoxyhémoglobine ou de grandeurs proportionnelles auxdites variations, dans la couche profonde, entre l’instant de mesure et l’instant de référence, à partir des ratios résultant de l) ; n ) à partir des estimation résultant de m), identification de l’occlusion comme étant une occlusion veineuse ou une occlusion artérielle. Le procédé peut comporter : o) estimation des variations des concentrations d’oxyhémoglobine et de désoxyhémoglobine, ou de grandeurs proportionnelles auxdites variations, dans la couche superficielle entre l’instant de mesure et l’instant de référence, à partir des ratios résultant de m). Selon une possibilité, dans l’étape h), la couche de base est : une couche élémentaire prédéterminée ; ou une couche élémentaire de rang différent, par exemple inférieur, de la couche élémentaire considéré, l’écart entre le rang de la couche élémentaire de base et la couche élémentaire considéré étant constant. Selon une possibilité, la source de lumière est une source laser. Un deuxième objet de l’invention est un dispositif, destiné à détecter une survenue d’une occlusion veineuse ou artérielle dans un tissu biologique, entre au moins un instant de mesure et un instant de référence, le tissu biologique étant délimité par une surface, le dispositif comportant : une source de lumière impulsionnelle, la source de lumière étant disposée en regard de la surface du tissu biologique, la source de lumière étant configurée pour émettre un faisceau d’illumination, formant une zone d’illumination, à la surface du tissu biologique, la source de lumière étant configurée pour émettre simultanément ou successivement le faisceau d’illumination dans deux longueurs d’onde d’illumination différentes l’une de l’autre, chaque longueur d’onde s’étendant dans une bande spectrale d’absorption de l’oxyhémoglobine et/ou de la désoxyhémoglobine ; un photodétecteur, configuré pour détecter des photons rétrodiffusés par le tissu biologique, après s’être propagés dans le tissu biologique, les photons rétrodiffusés émanant d’une zone de détection à la surface du tissu biologique ; une carte de comptage, destinée à établir une distribution temporelle, représentant un nombre de photons détectés par le photodétecteur en fonction du temps, à chaque longueur d’onde ; une unité de traitement, programmée pour mettre en œuvre les étapes d) à i) d’un procédé selon le premier objet de l’invention à partir des distributions temporelles établies par la carte de comptage à chaque longueur d’onde, à l’instant de mesure ou à chaque instant de mesure. L'invention sera mieux comprise à la lecture de l'exposé des exemples de réalisation présentés, dans la suite de la description, en lien avec les figures listées ci-dessous. FIGURES La schématise un dispositif permettant une mise en œuvre de l’invention. La montre une zone d’illumination et une zone de détection formées à la surface d’un échantillon. La schématise une impulsion d’illumination et une distribution temporelle de photons rétrodiffusés par l’échantillon, et détectés par un photodétecteur, les photons émanant de l’échantillon à travers la zone de détection. La illustre une configuration selon laquelle le dispositif comporte plusieurs fibres d’illumination et de détection. La illustre une distribution temporelle de photons détectés segmentée en différents intervalles temporels. La schématise une subdivision de l’échantillon en différentes couches, chaque couche correspondant à une profondeur atteinte par des photons détectés durant un intervalle temporel. La montre les principales étapes permettant une mise en œuvre de l’invention. La montre une segmentation de l’échantillon en deux couches, séparées par une interface s’étendant à une profondeur de séparation par rapport à la surface de l’échantillon. La montre une distribution temporelle de photons détectés segmentée en deux intervalles temporels, chaque intervalle temporel correspondant respectivement à une première couche et à une deuxième couche de l’échantillon. La schématise un échantillon utilisé lors d’essais expérimentaux. La représente un profil, montrant une évolution, à plusieurs instants de mesure successifs, d’une grandeur d’intérêt, proportionnelle à une concentration différentielle apparente dans différentes profondeurs de l’échantillon. Une concentration différentielle apparente est une comparaison entre des concentrations apparentes respectivement à un instant de mesure et à un instant de référence. La correspond à l’évolution de la quantité d’oxyhémoglobine, au cours d’une occlusion artérielle. La représente un profil, montrant une évolution, à plusieurs instants de mesure successifs, de la grandeur d’intérêt décrite en lien avec la . La correspond à l’évolution de la quantité de désoxyhémoglobine, au cours de la même occlusion artérielle. La figure 5D représente un profil montrant une évolution, à plusieurs instants de mesure successifs, d’une comparaison entre les grandeurs d’intérêts représentées sur la 5B entre deux couches élémentaires de rang et -3, pour différents ordres w successifs La figure 5E représente un profil montrant une évolution, à plusieurs instants de mesure successifs, d’une comparaison entre les grandeurs d’intérêts représentées 5C entre deux couches élémentaires de rang et -3, pour différents ordres w successifs La montre des coefficients de corrélation entre les profils représentés sur la . La montre des coefficients de corrélation entre les profils représentés sur la . La résulte du produit des coefficients de corrélation obtenues en 6A et 6B. La représente des estimations, selon le temps, de concentrations d’oxyhémoglobine et de désoxyhémoglobine résultant respectivement d’une mise en œuvre d’une méthode de spectrométrie de réflectance diffuse résolue en temps sans prise en compte de la couche superficielle (courbes a et b respectivement) et d’une mise en œuvre d’une méthode de spectrométrie de réflectance diffuse non résolue en temps, selon l’art antérieur (courbes c et d). La est représentative d’une occlusion artérielle. La représente des concentrations d’oxyhémoglobine (courbe a) et de désoxyhémoglobine (courbe b) illustrées dans la , corrigée selon la mise en œuvre de l’invention ainsi qu’un lissage temporel des concentrations estimées à différents instants de mesure. La montre une estimation des concentrations d’oxyhémoglobine (courbe a) et de désoxyhémoglobine (courbe b) au cours d’une occlusion artérielle selon la mise en œuvre de l’invention, ainsi que le résultat de mesures résultant d’une sonde Licox. (courbe c). La est équivalente à la , dans le cas d’une occlusion veineuse. EXPOSE DE MODES DE REALISATION PARTICULIERS La représente un dispositif 1 permettant un suivi de la vascularisation d’un lambeau enfoui sous un une couche de peau. L’objectif est de déterminer une survenue d’une occlusion veineuse ou artérielle, en particulier durant les premières heures ou les premiers jours après la chirurgie. Le dispositif comporte une source de lumière 10 impulsionnelle. La source de lumière 10 est impulsionnelle. Elle est configurée pour émettre un faisceau lumineux 11 impulsionnel, la durée de chaque impulsion étant typiquement entre 1 ps et 100 ps. La fréquence d’impulsion peut être par exemple comprise entre 1 et 100 MHz. Le faisceau lumineux 11, dit faisceau d’illumination, se propage en direction d’un échantillon de tissu biologique 20 à analyser. L’échantillon 20 est délimité par une surface 21. L’intersection du faisceau d’illumination 11 et de la surface 21 de l’échantillon forme une zone d’illumination 12, spatialement délimitée. La zone d’illumination est de préférence ponctuelle : elle est par exemple inscrite dans un cercle de diamètre inférieur à 1mm, voire inférieur à quelques centaines de µm, par exemple 100 µm ou 1 mm. La zone élémentaire d’illumination 12 est représentée sur la . La source de lumière peut être disposée au contact du tissu biologique 20 ou à distance de ce dernier. De préférence, la source de lumière est une source de lumière laser femtoseconde. Dans l’exemple représenté, la source de lumière est disposée à distance de l’échantillon. Le faisceau lumineux 11 est transporté jusqu’à la surface 21 de l’échantillon par une fibre optique d’illumination 10’. Les photons formant le faisceau d’illumination 11 se propagent dans le tissu biologique à analyser. Le tissu biologique 20 est formé d’un milieu diffusant, susceptible d’absorber les photons, les propriétés de propagation des photons dans le milieu dépendant notamment de propriétés d’absorption ou de diffusion dans le milieu. De façon usuelle, les propriétés d’absorption peuvent être quantifiées par un coefficient d’absorption linéaire . De façon connue, le coefficient d’absorption linéaire quantifie une probabilité d’absorption par le milieu par unité de longueur, à la longueur d’onde λ. Il est usuellement exprimé en cm -1 . Les propriétés de diffusion peuvent être quantifiées par un coefficient de diffusion ou un coefficient de diffusion réduit , qui quantifient une probabilité de diffusion par le milieu par unité de longueur, à la longueur d’onde λ. Il est usuellement exprimé en cm -1 . Le dispositif comporte un photodétecteur 16. Il peut par exemple s’agir d’un tube photomultiplicateur, de type tube photomultiplicateur hybride, usuellement désigné par l’acronyme HPM (Hybrid Photo Multiplier – Photomultiplicateur Hybride). Les photons détectés par le photodétecteur 16 émanent d’une zone élémentaire de détection 14 à la surface 21 de l’échantillon 20. La zone élémentaire de détection 14 est de préférence ponctuelle, en étant, à l’instar de la zone élémentaire d’illumination 12, inscrite dans un diamètre inférieur à 1 mm voire à 250 µm. La zone élémentaire de détection 14 est confondue ou séparée de la zone élémentaire d’illumination 12. La distance entre la zone élémentaire d’illumination 12 et la zone élémentaire de détection 14 est une distance de rétrodiffusion . Elle peut être de quelques centimètres ou de l’ordre du centimètre, voire peut être inférieure à 1 cm. Les photons 13 émanant de l’échantillon, à travers la zone élémentaire de détection 14, sont des photons rétrodiffusés par l’échantillon selon la distance de rétrodiffusion . Le photodétecteur 16 peut être disposé au contact ou à distance de l’échantillon 20. Dans l’exemple représenté, une fibre optique de détection 16’ transporte les photons rétrodiffusés 13 entre la zone élémentaire de détection 14 et le photodétecteur 16. La source de lumière 10 peut être monochromatique ou polychromatique. Dans l’exemple représenté, la source de lumière est un laser femtoseconde de type super continuum (supercontinu) émettant des impulsions lumineuses selon une fréquence de 80 MHz. Une sélection de la longueur d’onde d’illumination est effectuée par un filtre optique 10 f . Dans l’exemple décrit, le filtre 10 f est une roue à filtres comportant différents filtres interférentiels passe-bande permettant une sélection de la longueur d’onde d’illumination selon l’une des valeurs suivantes : 750 nm, 800 nm ou 850 nm. Le filtre 10 f peut également comporter un filtrage passe-haut, en complément des filtres interférentielles passe-bande. Dans le mode de réalisation de l’invention, on utilise deux longueurs d’onde d’illumination et qui correspondent à des longueurs d’onde d’absorption de la désoxyhémoglobine et de l’oxyhémoglobine. Par exemple, 750 nm et 850 nm. La notation désigne indifféremment une longueur d’onde d’illumination, qu’il s’agisse de ou de . Le photodétecteur 16 est relié à une carte électronique de comptage 17 de type TCSPC (Time-Correlated Single Photon Counting – Comptage de photon unique résolu en temps). Ce type de carte électronique de comptage est bien connu dans le domaine des mesures optiques résolues en temps. La carte de comptage établit une distribution temporelle des photons détectés par le photodétecteur 16. La distribution temporelle est établie statistiquement suite à une multitude d’impulsions lumineuses émises successivement par la source de lumière 10. La distribution temporelle correspond à un histogramme représentant un nombre de photons rétrodiffusés 13 détectés (axe des ordonnées) en fonction du temps (axe des abscisses). La carte de comptage 17 est reliée à une unité de traitement 18. L’unité de traitement 18 comporte un microprocesseur programmé pour mettre en œuvre les étapes décrites par la suite, en lien avec la figure 3, à partir des distributions temporelles établies par la carte de comptage 17. Le dispositif peut comporter un ou plusieurs atténuateurs optiques, chaque atténuateur optique étant usuellement désigné « densité optique ». L’homme du métier adapte l’atténuation en fonction de la puissance de la source de lumière et de la sensibilité du photodétecteur. La carte de comptage 17 est synchronisée par une ligne de synchronisation 15, qui détecte chaque impulsion lumineuse émise par la source de lumière 10. Cela permet de déterminer un intervalle temporel dt entre l’émission et la détection de chaque photon détecté par le photodétecteur 16. Sur la , on a représenté une impulsion lumineuse de photons incidents 11, survenant à un instant d’illumination t 1 0 , et une distribution temporelle de photons rétrodiffusés 13, cette dernière s’étendant selon un intervalle temporel Δt 13 . En multipliant le nombre de zones élémentaires d’illumination 12 et de zones élémentaires de détection 14, il est possible d’adresser séquentiellement différentes distances de rétrodiffusion . La représente une configuration selon laquelle la source de lumière 10 est reliée à cinq fibres optiques d’illumination 10’, par l’intermédiaire d’un commutateur optique (switch optique). Les fibres optiques d’illumination 10’ sont réparties respectivement au centre et au milieu des côtés d’un carré de côté 12 mm. Selon cette configuration, le dispositif 1 comporte quatre fibres optiques de détection 16’, respectivement reliées à quatre photodétecteurs 16 différents. Les fibres optiques de détection 16’ sont respectivement disposées sur les quatre coins du carré de côté 12 mm précédemment décrit. Une telle configuration permet d’adresser des distances de rétrodiffusion respectivement égales à 6 mm, 8.5 mm et 13.4 mm. La puissance de la source de lumière incidente est monitorée par un prélèvement de la lumière d’excitation assuré par une fibre optique 15’. Sur la figure 2A, on a schématisé une segmentation d’une distribution temporelle entre différentes fenêtres temporelles centrées correspondant à des intervalles temporels . Chaque intervalle temporel s’étend autour d’un instant de détection . La segmentation de la distribution temporelle permet d’estimer un nombre de photons détectés durant chaque intervalle temporel , ce qui correspond une intensité temporelle pour chaque intervalle temporel . L’indice est un entier correspondant au rang de chaque intervalle temporel avec , étant le nombre d’intervalles temporels considérés. La durée de chaque intervalle temporel peut être inférieure à une nanoseconde, en étant par exemple de l’ordre de quelques dizaines de ps . Cette segmentation temporelle peut être faite, à partir de la distribution temporelle, par application d’une fenêtre temporelle en créneau de largeur fixe, les fenêtres temporelles étant décalées les unes par rapport aux autres. Les fenêtres temporelles peuvent se superposer l’une à l’autre. Par exemple, la largeur de chaque fenêtre temporelle peut être de 0.2 ns, le décalage de deux fenêtres temporelles consécutives s’élevant à 0.1 ns. Selon une autre possibilité, les fenêtres temporelles peuvent être définies en appliquant des transformées de Mellin Laplace de différents ordres à la distribution temporelle , comme décrit dans EP2605002 ou dans les publications : L. Hervé et al. « Time-domaine diffuse optical tomography processing by suing the Mellin-Laplace transform », Applied Optics 51(25), 5978-5988 (2012) ; J. Zaouaoui « Chromophore decomposition in multispectral time-resolved diffuse optical tomography », Biomedical optics Express 8(10), 4772 (2017). De façon connue, à chaque intervalle temporel correspond une profondeur susceptible d’être atteinte par les photons dans chaque couche de l’échantillon 20. Plus le rang de l’instant est élevé, plus les photons détectés peuvent être considérés comme ayant atteint une profondeur importante dans l’échantillon. Comme on peut le voir sur la figure 2B, La profondeur est déterminée par rapport à la surface 21 de l’échantillon, comportant la zone élémentaire d’illumination 12. Ainsi, à chaque intervalle temporel centré sur l’instant de détection peut être associé à une couche élémentaire de rang s’étendant à une profondeur , où est l’épaisseur de la couche élémentaire, qui dépend de la largeur de l’intervalle temporel considéré. Comme précédemment évoqué, le dispositif est configuré pour permettre une détection d’une occlusion veineuse ou artérielle dans un lambeau enfoui sous une couche de peau. Le dispositif est appliqué contre la surface 21 d’un tissu biologique 20, ce dernier comportant une couche superficielle 20 1 , de peau, et une couche profonde 20 2 . La couche superficielle est interposée entre la surface du tissu et la couche profonde. Sous l’effet d’une occlusion artérielle ou veineuse, la concentration d’hémoglobine dans la couche profonde peut subir une variation progressive par rapport à la couche superficielle, entre l’instant de référence et l’instant de mesure. Par hémoglobine, il est entendu les deux états de l’hémoglobine, à savoir l’oxyhémoglobine et la désoxyhémoglobine. Deux longueurs d’onde d‘illumination sont successivement utilisées qui correspondent à des longueurs d’onde absorbée par l’oxyhémoglobine et la désoxyhémoglobine. Le terme couche est à distinguer du terme « couche élémentaire ». Une couche correspond à une partie macroscopique de l’échantillon dans laquelle les propriétés optiques sont considérées comme homogènes. Une couche élémentaire correspond à une « microcouche », dont l’épaisseur et la profondeur dépendent de la segmentation temporelle effectuée sur la distribution temporelle. Comme représenté sur la figure 1A, la couche superficielle et la couche profonde sont séparées l’une de l’autre par une interface 20 1/2 . L’interface est située à une profondeur d’interface par rapport à la surface 21 de l’échantillon. Cependant, la profondeur d’interface n’est pas connue. La profondeur d’interface d 1/2 peut être située à quelques mm, par exemple 3 mm ou 5 mm, à quelques cm, par exemple 1 ou 2 cm, par rapport à la surface 21 de l’échantillon. Un aspect important de l’invention est d’estimer expérimentalement des instants de la distribution temporelle, respectivement associés à la couche superficielle et à la couche profonde. Il s’agit de discriminer précisément : les photons prompts, considérés comme s’étant propagés à travers la couche superficielle 20 1 , sans s’être propagés à travers la couche profonde 20 2 , ou de façon marginale : sur la , la trajectoire des photons prompts est matérialisée par des lignes en tirets. les photons retardés, considérés comme s’étant propagés à travers la couche superficielle 20 1 et à travers la couche profonde 20 2 : sur la , la trajectoire des photons retardés est matérialisée par des lignes en tirets mixtes. La figure 3 schématise les principales étapes du procédé, dans un cas d’application pour lequel on cherche à déterminer la survenue d’une occlusion veineuse ou artérielle dans la couche profonde du tissu biologique examiné. Il s’agit des étapes 100 à 130. Le procédé peut également permettre une identification d’une occlusion détectée : il s’agit de savoir si l’occlusion est veineuse ou artérielle. Cela correspond aux étapes 140 et 150. Au cours des étapes 140 et 150, on cherche à estimer une variation de la concentration d’hémoglobine 22 k , par rapport à une concentration de référence . L’hémoglobine peut être de l’oxyhémoglobine (k = 1) ou la désoxyhémoglobine (k=2). La concentration de référence correspond à une concentration de l’hémoglobine à un instant de référence . On se place dans un cas de figure selon lequel la concentration de l’hémoglobine est susceptible de varier dans le temps, et de manière différente dans les deux couches. Un aspect important de l’invention est de pouvoir déterminer une variation de concentration d’oxyhémoglobine et de désoxyhémoglobine (ou au moins une grandeur proportionnelle à chaque variation) dans la couche profonde, sans qu’il soit nécessaire de déterminer la variation de concentration d’oxyhémoglobine et de désoxyhémoglobine dans la couche superficielle. Le dispositif permet ainsi d’effectuer des mesures de réflectance optique résolues en temps peuvent permettre un suivi de l’oxygénation du lambeau, de façon non invasive. Le principe de l’invention est schématisé sur les figures 4A et 4B. Sur la , on a représenté l’échantillon, avec les photons respectivement issus de la couche superficielle 20 1 et de la couche profonde 20 2 . Lorsqu’une occlusion veineuse ou artérielle survient dans la couche profonde, la réponse de la couche profonde à l’égard du faisceau d’illumination varie. Cela se traduit par une variation du signal rétrodiffusé. Ce signal est sensible aux variations de la couche profonde, mais aussi de la couche superficielle : si la couche superficielle évolue (du fait de la chirurgie par exemple), alors le signal rétrodiffusé va traduire l’évolution des deux couches ; pour être à même de détecter et identifier correctement l’occlusion survenue en couche profonde, il faut corriger ce signal de la contribution de la couche superficielle : c’est l’objet de l’invention. Un premier objectif de l’invention est décrit en lien avec les étapes 100 à 130 décrites ci-après : il s’agit de détecter la survenue d’une occlusion veineuse ou artérielle dans la couche profonde 20 2 , sachant que la profondeur de l’interface 20 1/2 entre les deux couches n’est pas connue. Un deuxième objectif, complémentaire au premier, et optionnel, est d’effectuer un suivi de la concentration d’oxyhémoglobine ou d’oxyhémoglobine dans la couche profonde 20 2 . Cela correspond aux étapes 140 et 150. Pour cela, on détermine expérimentalement un premier intervalle temporel et un deuxième intervalle temporel , qui correspondent respectivement aux contributions des photons « prompts » et des photons retardés à la distribution temporelle T. Détection de la survenue d’une occlusion Etape 100 : acquisition Au cours de cette étape, l’échantillon 20 est illuminé de façon impulsionnelle. A partir de photons rétrodiffusés 13 détectés par le photodétecteur 16 après chaque impulsion, on constitue une distribution temporelle à chaque longueur d’onde d’émission . L’étape 100 peut être effectuée séquentiellement ou simultanément pour plusieurs longueurs d’onde , de façon à obtenir, pour chaque longueur d’onde, une distribution temporelle . L’indice désigne le rang de chaque longueur d’onde. La distribution temporelle représente une évolution d’une intensité des photons rétrodiffusés détectés par le photodétecteur 16., en fonction du temps Selon un mode de réalisation, l’intensité peut être normalisée par l’intensité émise par la source de lumière 10. Dans ce cas, une fibre, dite de retour d’excitation, s’étend entre la source de lumière 10 et un photodétecteur 16. Une telle normalisation permet de prendre en compte les fluctuations de la source de lumière. Dans la suite, le terme intensité désigne indifféremment une intensité ou une intensité normalisée par l’intensité du faisceau d’illumination. Dans cet exemple, l’étape 100 est mise en œuvre à différentes instants de mesure , de façon à permettre un contrôle régulier de l’échantillon examiné. Elle est notamment mise en œuvre à un instant de référence , lors de laquelle l’échantillon examiné est considéré comme dans un état de référence. Lorsque l’échantillon est un tissu cutané ayant fait l’objet d’une reconstruction chirurgicale par un lambeau, l’instant de référence peut correspondre à une mesure effectuée immédiatement suite à l’intervention chirurgicale, l’échantillon étant considéré comme correctement oxygéné. Au cours de l’instant de référence , on obtient une distribution temporelle de référence à chaque longueur d’onde . Au cours de chaque instant de mesure , la distribution temporelle d’intensité est notée . Etape 110 : détermination d’intensités temporelles. Au cours de cette étape, chaque distribution temporelle , fait l’objet d’une décomposition en fenêtres temporelles, de façon à obtenir une intensité temporelle de photons rétrodiffusés détectés dans l’intervalle temporel à la longueur d’onde . Le terme intensité temporelle désigne une intensité détectée dans chaque intervalle temporel d . Comme précédemment indiqué, la décomposition temporelle peut être effectuée par une application d’une transformée de Mellin-Laplace à la distribution temporelle , ou par l’application de fenêtres temporelles en créneaux (ou d’autres types de fenêtres temporelles). A l’issue de cette étape, on dispose d’intensités temporelles de la distribution temporelle chaque intensité temporelle étant associée à un intervalle temporel centré sur un instant de détection . Le référentiel temporel des mesures est donné par les instants de mesure (espacés de quelques minutes à plusieurs heures). Chaque instant de mesure est décomposé en intensité temporelles pour des instants successifs , rapprochés d’une durée de quelques dizaines de ps ou centaines de picosecondes (ps) ( par exemple 200ps). A chaque instant de détection correspond à un intervalle temporel de détection s’étendant de part et d’autre de l’instant de détection . Etape 120 : Explicitation de la concentration différentielle. Au cours de cette étape, on détermine, à partir de chaque intensité temporelle résultant de l’étape 110, une concentration de l’hémoglobine 22 k dans chaque couche élémentaire traversée par les photons détectés dans la fenêtre temporelle de rang , s’étendant autour de l’instant de détection . Les inventeurs considèrent que selon une première approximation, correspond à une variation du coefficient d’absorption, à la longueur d’onde , déterminée d’après la fenêtre temporelle de rang , entre l’instant de mesure et l’instant de référence : et sont respectivement un coefficient d’absorption à l’instant de mesure et à l’instant de référence , pour la longueur d’onde . est le libre parcours moyen des photons sélectionnés sur . En première approximation, est indépendant de la longueur d’onde ; il peut être estimé par étant un coefficient de proportionnalité. Ainsi, l’étape 120 étape comporte une comparaison, dans chaque fenêtre temporelle de rang , entre une intensité et une intensité de référence , établies, pour la même fenêtre temporelle, respectivement à partir des distributions temporelles et . A partir de (1), on peut écrire : Or est la concentration de l’hémoglobine (k = 1 : oxyhémoglobine ; k=2 désoxyhémoglobine k =2) lors de l’instant de mesure ; est un coefficient d’absorption linéaire l’hémoglobine pour une concentration unitaire (unité cm²/g lorsque la concentration est exprimée en g/cm 3 ), également désigné coefficient d’extinction. est une donnée connue, disponible dans des bases de données. Lors de l’instant de référence, est la concentration de l’hémoglobine lors de l’instant de référence ; La combinaison de (2), (3) et (4) donne : avec : correspond à une variation de la concentration l’hémoglobine , correspondant à la fenêtre temporelle de rang , entre l’instant de mesure et l’instant de référence . Par la suite, est désignée par le terme «concentration différentielle apparente». Lorsque deux longueurs d’onde , sont mises en œuvre, l’expression (5) peut être présentée sous forme matricielle : où : Il est à noter que le recours aux expressions (1), (2), (5) et (7) repose sur une hypothèse simplificatrice, selon laquelle la concentration d’hémoglobine, à une couche élémentaire , d’ordre , est estimée sans prendre en compte des inhomogénéités du milieu, à des couches élémentaires d’ordre inférieur . Or, comme précédemment évoqué, l’échantillon peut présenter des inhomogénéités de la concentration d’hémoglobine en fonction de la profondeur. L’hypothèse simplificatrice utilisée au cours de cette étape revient à considérer que la concentration de l’hémoglobine est homogène entre la surface de l’échantillon et la couche élémentaire Ainsi, chaque valeur peut être assimilée à une concentration différentielle apparente, le terme « apparente » désignant le fait qu’elle est estimée sur la base d’une hypothèse d’homogénéité de la concentration de l’hémoglobine dans l’échantillon, jusqu’à la profondeur . Une telle simplification est justifiée par le fait que chaque concentration différentielle est utilisée, dans les étapes 130 et 140, pour effectuer une segmentation de l’échantillon entre une première couche 20 1 et une deuxième couche 20 2 , la concentration de l’hémoglobine dans chaque couche ainsi segmentée étant considérée comme homogène. On considère à présent, pour chaque couche élémentaire une grandeur d’intérêt , telle que : La grandeur d’intérêt est proportionnelle à une variation de la concentration différentielle apparente de la couche élémentaire entre l’instant de référence et l’instant de mesure. La grandeur peut être considérée comme une concentration surfacique apparente. Elle correspond à la concentration apparente multipliée par le libre parcours moyen des photons . La grandeur est équivalente à une concentration surfacique apparente, obtenue comme la contribution du milieu entre la surface à la couche élémentaire d’ordre . Chaque grandeur d’intérêt ( pour l’oxyhémoglobine et pour la désoxyhémoglobine) peut être calculée directement à partir de et de en résolvant l’équation (10), qui résulte de l’équation (7) : Etape 130 : comparaison des grandeurs d’intérêt résultant de l’étape 120. L’invention consiste à comparer ces grandeurs d’intérêt obtenues en considérant respectivement deux couches élémentaires de rangs différents Pour chaque couche élémentaire , on obtient une grandeur d’intérêt différentielle , correspondant à une différence entre : une grandeur d’intérêt , obtenue pour une couche élémentaire de base , de rang (avec ) une grandeur d’intérêt, obtenue pour la couche élémentaire de rang . Dans cet exemple, on établit une comparaison telle que La couche élémentaire de base peut être de rang fixe, par exemple = 1, correspondant aux photons les plus prompts de la distribution, auquel cas : La couche élémentaire de base peut être telle que , auquel cas Lorsque les intensités temporelles , sont obtenues en appliquant une transformée de Mellin Laplace, la couche temporelle de base peut être tel que . Il est préférable que le rang de la couche temporelle de base soit fixe, ou obtenu par un décalage fixe par rapport au rang de la couche élémentaire considérée, par exemple ou Pour chaque couche élémentaire , la quantité représente une variation, en fonction de , de la grandeur d’intérêt entre la couche élémentaire et la couche élémentaire de base . A partir des intensités temporelles résultant de l’étape 110, on peut calculer, pour différents instants de mesure T, une grandeur d’intérêt différentielle pour différentes couches élémentaires d’ordre . L’évolution, en fonction de , de permet d’identifier d’éventuelles discontinuités survenant, au cours du temps, dans l’échantillon. Lors de la survenue d’une occlusion, veineuse ou artérielle, dans la couche profonde, une discontinuité apparaît entre , pour les rangs correspondant à la couche superficielle, et pour les rangs correspondant à la couche profonde. Autrement dit, lorsqu’une occlusion survient, on observe une décorrélation entre les valeurs de obtenues pour les couches élémentaires situées respectivement dans la couche superficielle 20 1 et dans la couche profonde 20 2 . Lorsque les étapes 100 à 130 sont mises en œuvre à différents instants de mesure, on peut obtenir, pour chaque couche élémentaire de rang , un profil traçant l’évolution temporelle de la grandeur d’intérêt différentielle . En cas de survenue d’une occlusion artérielle ou veineuse, la décorrélation entre la couche superficielle et la couche profonde est plus aisée à détecter sur la base d’un tel profil. La décorrélation entre les valeurs de obtenues respectivement entre la couche superficielle et la couche profonde peut être déterminée visuellement, par un utilisateur du dispositif. De préférence, un indicateur de corrélation est calculé, qui permet de quantifier une corrélation entre deux valeurs de pour deux couches de rang différents. Il peut s’agir d’un coefficient de corrélation r². Plus la valeur de l’indicateur de corrélation est élevée, plus les valeurs de , pour deux couches de rang différents, sont proches. En cas de survenue d’une occlusion, l’indicateur de corrélation entre deux couches élémentaires de rangs différents franchit un seuil prédéterminé, en dessous duquel on considère qu’il y a une décorrélation. Essais expérimentaux On a mis en œuvre le procédé précédemment décrit, en utilisant un échantillon comportant un lambeau de cochon implanté en dessous d’une couche de peau, tel que représenté sur la . La couche de peau formait couche superficielle 20 1 tandis que le lambeau formait une couche profonde 20 2 . Le lambeau était relié à une artère et à une veine. Un clamp permettait de bloquer l’irrigation soit de l’artère (occlusion artérielle), soit de la veine (occlusion veineuse). Le dispositif mettait en œuvre une source de lumière et des photodétecteurs fibrés, comme représenté sur la . Les paramètres expérimentaux étaient les suivants : épaisseur de la couche superficielle (peau + graisse) : environ 1 cm ; épaisseur du lambeau : environ 1 cm ; source de lumière utilisée : laser supercontinuum fibré - Fianium Wl-SC-480-10 fréquence d’impulsion 80 MHz ; photodétecteurs : photomultiplicateur hybride - Becker-Hickl HPM 100-50. Un filtre passe haut – fréquence de coupure à 690 nm – était associé à chaque photodétecteur. Référence filtre : HLP690 - Melles Griot. carte électronique TCSPC : Becker-Hickl TCSPC-SPCM 134 sélection de la longueur d’onde : premier filtrage par deux filtres passe haut Thorlabs FEL 1000 et FGL 570, couplés à une roue à filtres comportant des filtres centrés sur les longueurs d’onde de 750 nm, 800 nm et 850 nm – Largeur à mi-hauteur 10 nm. Un atténuateur de lumière était disposé entre chaque fibre d’illumination 10’ et la source de lumière : densité optique variable entre 0 et 4. Distance de rétrodiffusion : 13 mm Le protocole d’acquisition a été décomposé comme suit entre T ref et T ref + 20 min (minutes) : le lambeau est irrigué normalement par la veine et l’artère, t 0 correspondant à une série d’acquisitions ; entre T ref + 20 min et T ref + 40 min : clampage de l’artère (occlusion artérielle) ou de la veine (occlusion veineuse) : le lambeau n’est plus irrigué correctement; entre T ref + 40 min et T ref + 60 min : le clampage est arrêté. Irrigation normale du lambeau. Entre T ref et T ref + 60 min, le procédé décrit dans les étapes 100 à 130 a été suivi, en différents instants de mesure , chaque instant de mesure étant espacée l’un de l’autre de 45 secondes. Lors de l’étape 120, on a extrait différentes intensités temporelles en appliquant des fonctions de segmentation temporelles en appliquant une transformée de Mellin Laplace. Lors de cette étape, on a considéré un nombre total d’intervalles temporels égal à 20. Dans l’exemple donné par les figures 5B à 5E, il s’agit d’une occlusion artérielle. Les figures 5B et 5D représentent respectivement, à chaque instant de mesure (axe des abscisses) des valeurs de et de pour les valeurs de comprises entre 4 et 20. Lorsque k=1, il est rappelé qu’on adresse l’oxyhémoglobine. Les figures 5C et 5E représentent respectivement, à chaque instant de mesure (axe des abscisses) des valeurs de et de pour les valeurs de comprises entre 4 et 20. Lorsque k=2, il est rappelé qu’on adresse la désoxyhémoglobine. Les courbes 5B à 5E correspondent à des profils temporels traçant les évolutions temporelles respective de , , , . Sur les figure 5B et 5C, avant le clampage (T ref ref + 20 min), les courbes sont proches les unes des autres, ce qui traduit une certaine homogénéité. Durant le clampage, (T ref + 20 ref + 40 min), les courbes se séparent les unes des autres, ce qui traduit l’apparition d’une inhomogénéité de la concentration d’oxyhémoglobine (courbe 5B) et de désoxyhémoglobine (courbe 5C). Sans surprise, les grandeurs d’intérêt , qui sont proportionnelles à une variation de la concentration d’oxyhémoglobine (ou de désoxyhémoglobine) entre l’instant de mesure et l’instant de référence, traduisent, au cours du clampage artériel dans cet exemple, une diminution de la concentration d’oxyhémoglobine et une légère augmentation de la concentration de désoxyhémoglobine. Pour définir s’il y a ou non occlusion, on cherche à mettre en évidence une différence d’évolution dans le tissu d’intérêt. Pour cela, on s’intéresse en particulier aux grandeurs d’intérêt différentielles . Pour chaque analyte (oxyhémoglobine et désoxyhémoglobine), les valeurs sont faiblement dispersées avant l’occlusion. L’occlusion engendre une dispersion plus importante. Sur la base des figures 5D et 5E, on note un comportement similaire pour les profils lorsque Les profils temporels représentés sur les figures 5B à 5E permettent une détection de la survenue d’une occlusion artérielle à partir des grandeurs d’intérêt différentielles . La détection d’une survenue d’une occlusion peut être effectuée visuellement, ou en calculant un indicateur de corrélation entre les valeurs de entre des couches élémentaires d’ordres différents. Chaque grandeur d’intérêt différentielle est représentative d’une variation spatiotemporelle d’une concentration d’oxyhémoglobine ou de désoxyhémoglobine. La variation temporelle correspond au terme , qui représente une variation de entre deux instants différents (instant de mesure et instant de référence). La variation spatiale provient de la différence entre les grandeurs entre la couche de rang et la couche de rang ’. Les figures 6A à 6C représentent respectivement des coefficients de corrélation calculés respectivement : à partir des profils de la (oxyhémoglobine) ; à partir des profils de la (desoxyhémoglobine) ; en combinant les profils de la et 5E. La figure 6A montre un indicateur de corrélation entre des profils de la figure 5D correspondant à des rangs différents. Sur cette figure, l’axe des abscisses et l’axe des ordonnées correspondent à un rang . Le niveau de gris correspond au coefficient de corrélation. De même la figure 6B montre un indicateur de corrélation entre des profils de la figure 5E correspondant à des rangs différents. La est une synthèse figures 6A et 6B : elle représente le produit des coefficients de corrélation présentés en et 6B. Sur cette figure, on observe une bonne corrélation des couches élémentaires 1 à 11, ainsi que des couches élémentaires 13 à 20. Estimation de variations de concentrations d’hémoglobine entre l’instant de mesure et l’instant de référence . On va décrire à présent des étapes optionnelles, qui permettent une analyse plus approfondie des tissus examinés. Après avoir détecté la survenue d’une inclusion, l’objectif est d’identifier le type d’inclusion, afin de déterminer s’il s’agit d’une occlusion veineuse ou artérielle. Etape 1 4 0 : segmentation temporelle La décorrélation entre les grandeurs d’intérêt ou les grandeurs d’intérêt différentielles, décrites dans l’étape 130 et sur les figures 5B à 5E, ainsi que 6A à 6C, permet de définir une couche élémentaire de rang 1/2 , de part et d’autre de laquelle les grandeurs d’intérêt, ou les grandeurs d’intérêt différentielles, ne sont pas corrélées. Une telle couche élémentaire correspond à une couche d’interface . Elle est considérée comme formant une interface entre la couche superficielle 20 1 d’une part, et la couche profonde 20 2 d’autre part. La couche d’interface est associée à un intervalle temporel d’interface : cf. figures 4A et 4B A partir de l’intervalle temporel d’interface , il est possible d’estimer : un premier intervalle temporel , centré sur un instant de détection correspondant aux photons prompts considérés comme n’ayant pas traversé la couche profonde 20 2 ; un deuxième intervalle temporel , centré sur un instant de détection , correspondant aux photons retardés considérés comme ayant traversé la couche superficielle 20 1 ainsi que la couche profonde 20 2 . Etape 150 : estimation d’une variation de concentration d’hémoglobine dans la couche profonde 20 2 , entre l’instant de mesure et l’instant de référence. Dans cette étape, on considère que l’on a deux couches : la couche 20 1 (couche superficielle) et la couche 20 2 (couche profonde). Sur les distributions temporelles et , les contributions , de la couche 20 1 et , de la couche 20 2 correspondent respectivement : aux photons détectés durant le premier intervalle temporel , aux photons détectés durant l’intervalle temporel ; A partir des mesures , , , on peut écrire : (15) correspond au libre parcours moyen dans la couche 20 1 : correspond à une moyenne des libres parcours moyens des photons contribuant à et . est représenté sur la figure 4A. De même, correspond au libre parcours moyen dans la couche 20 2 . correspond à une moyenne des libres parcours moyens des photons contribuant à et En utilisant le fait que et il vient (17) Sachant que, la matrice est connue et que les grandeurs et , et résultent des distributions temporelles mesurées à l’instant de mesure et à l’instant de référence, à chaque longueur d’onde, on obtient , c’est-à-dire des grandeurs et , proportionnelles aux variations de concentrations et respectivement de l’oxyhémoglobine et de la désoxyhémoglobine dans la deuxième couche 20 2 . Cela permet d’identifier le type d’occlusion, c’est-à-dire de déterminer s’il s’agit d’une occlusion veineuse ou artérielle. Si l’on souhaite estimer les concentrations et , le libre parcours moyen peut être estimé. Il est également possible d’estimer et , dans la première couche, à partir de l’expression (18) A partir d’une estimation de , on peut obtenir et . Essais expérimentaux A partir des figures 5B à 5E et 6A à 6C, on a estimé la couche d’interface de rang w 1/2 entre la première couche 20 1 et la deuxième couche 20 2 . On a mis en œuvre les étapes 140 et 150 de façon à estimer, à chaque instant de mesure, des grandeurs et respectivement proportionnelles aux variations de concentrations d’oxyhémoglobine et de désoxyhémoglobine entre chaque instant de mesure et l’instant de référence. La représente des estimations, selon le temps, de concentrations d’oxyhémoglobine et de désoxyhémoglobine résultant respectivement d’une mise en œuvre d’une méthode de spectrométrie de réflectance diffuse résolue en temps sans prise en compte de la couche superficielle (courbes a et b respectivement) et d’une mise en œuvre d’une méthode de spectrométrie de réflectance diffuse non résolue en temps, selon l’art antérieur (courbes c et d). La est représentative d’une occlusion artérielle. La figure 7B représente des variations de concentrations d’oxyhémoglobine (courbe a) et de désoxyhémoglobine (courbe b) illustrées dans la figure 7A, corrigée selon la mise en œuvre de l’invention ainsi qu’un lissage temporel des concentrations estimées à différents instants de mesure. Sur la figure 7B, l’axe des ordonnées correspond aux variations de concentrations d’oxyhémoglobine et de désoxyhémoglobine, dans la couche 20 2 , multipliées par : il s’agit de et . Lors des essais expérimentaux, des sondes Licox (marque déposée) ont été mises en œuvre dans la couche superficielle et dans la couche profonde. Une sonde Licox, connue de l’homme du métier, permet de mesurer une pression partielle d’oxygène dans les tissus. La figure 8A représente une estimation de la variation de concentration d’oxyhémoglobine (courbe a) et de désoxyhémoglobine (courbe b), chaque variation étant multipliée par résultant de la mise en œuvre des étapes 100 à 150, lors d’une occlusion artérielle. On a également représenté une courbe résultant de la sonde Licox (courbe c). La figure 8B représente une estimation de la variation de concentration d’oxyhémoglobine (courbe a) et de désoxyhémoglobine (courbe b), chaque variation étant multipliée par résultant de la mise en œuvre des étapes 100 à 150, lors d’une occlusion veineuse. On a également représenté une courbe résultant de la sonde Licox (courbe c). Une sonde Licox ne permet pas de remonter à une estimation d’une variation concentration. De plus, il s’agit d’une mesure invasive. La sonde Licox permet de détecter la survenue d’une occlusion, qu’il s’agisse d’une occlusion artérielle ou veineuse, sans pouvoir discriminer le type d’occlusion. On observe que les mesures selon l’invention, non invasives, sont cohérentes avec les mesures Licox pour détecter la survenue d’une occlusion. Elles permettent une estimation d’une variation de la concentration d’oxyhémogline et de désoxyhémoglobine par rapport à un instant de référence. On peut ainsi identifier le type d’occlusion : artérielle ou veineuse. Les mesures telles qu’illustrées sur les figures 8A et 8B ont été effectuées sur quinze cochons différents : les occlusions artérielles ou veineuses ont été correctement identifiées dans chaque essai. Procédé de suivi de la vascularisation d’un tissu biologique entre au moins un instant de mesure ( ) et un instant de référence ( ), le tissu biologique étant délimité par une surface (21), le tissu biologique comportant une couche superficielle (20 1 ) et une couche profonde (20 2 ), la couche superficielle s’étendant entre la surface du tissu et la couche profonde, le procédé comportant : a) illumination du tissu biologique par une source de lumière impulsionnelle (10), la source de lumière étant disposée en regard de la surface du tissu biologique, la source de lumière émettant un faisceau d’illumination (11) formant une zone d’illumination (12) à la surface du tissu biologique, le faisceau d’illumination étant émis simultanément ou successivement dans deux longueurs d’onde d’illumination différentes l’une de l’autre, chaque longueur d’onde s’étendant dans une bande spectrale d’absorption de l’oxyhémoglobine et/ou de la désoxyhémoglobine ; b) à chaque longueur d’onde, détection de photons rétrodiffusés par le tissu biologique, après s’être propagés dans le tissu biologique, par un photodétecteur (16), les photons rétrodiffusés détectés émanant d’une zone de détection (14) à la surface du tissu biologique ; c ) à chaque longueur d’onde, obtention d’une distribution temporelle ( ), représentant un nombre de photons détectés par le photodétecteur en fonction du temps, la distribution temporelle s’étendant selon une durée, à partir d’un instant initial ; les étapes a) à c) étant effectuées à l’instant de référence et à l’instant de mesure ou chaque instant de mesure ; le procédé étant caractérisé en ce qu’il comporte, à l’instant de mesure ou chaque instant de mesure et à l’instant de référence : d) segmentation de la durée de chaque distribution temporelle en intervalles de détection ( ), à chaque intervalle de détection étant associé un rang ( ), le rang étant d’autant plus élevé que l’intervalle de détection est éloigné de l’instant initial de la distribution temporelle, chaque intervalle de détection étant associé à une couche élémentaire ( ) s’étendant à une profondeur ( w ) dans le tissu biologique ; e) pour chaque couche élémentaire ( ), à partir de chaque distribution temporelle résultant de c), établies à l’instant de mesure et à l’instant de référence, détermination d’intensités temporelles ( , ), correspondant à un nombre de photons détectés durant chaque intervalle de détection ( ) ; et en ce qu’il comporte, à l’instant de mesure ou à chaque instant de mesure f) pour chaque intervalle de détection ( ), calcul d’un ratio entre l’intensité temporelle à l’instant de mesure et l’intensité temporelle à l’instant de référence ; g) à partir des ratios résultant de f), établis pour chaque longueur d’onde, calcul d’une grandeur d’intérêt ( ) pour chaque couche élémentaire respectivement associée à chaque intervalle de détection, la grandeur d’intérêt étant représentative d’une variation de la concentration d’hémoglobine, dans ladite couche élémentaire, entre l’instant de mesure et l’instant de référence ; h) pour chaque couche élémentaire ( ), comparaison de la grandeur d’intérêt, calculée lors de g), avec une grandeur d’intérêt de base, calculée pour une couche élémentaire de base, la couche élémentaire de base étant une des couches élémentaires du tissu biologique, de façon à obtenir, pour chaque couche élémentaire, une grandeur d’intérêt différentielle ; i) détermination d’une évolution, selon la profondeur du tissu biologique, des grandeurs d’intérêt différentielles . Procédé selon l’une quelconque des revendications précédentes, dans lequel : les étapes a) à h) sont mises en œuvre à différents instants de mesure ; l’étape i) comporte une formation d’un profil, pour chaque couche de élémentaire ( ), chaque profil représentant une variation de la grandeur d’intérêt différentielle , dans ladite couche élémentaire, en fonction de l’instant de mesure ( ) ; l’étape i) comporte un calcul d’un critère de corrélation entre différents profils, correspondant respectivement à différentes couches élémentaires, le critère de corrélation traduisant une corrélation entre les différents profils. Procédé selon la revendication 2, comportant, suite à i) : j) détermination d’une couche élémentaire d’interface ( ), choisie parmi les différentes couches élémentaires, la couche élémentaire d’interface étant une couche de part et d’autre de laquelle le critère de corrélation franchit la valeur seuil, la couche élémentaire d’interface étant considérée comme formant une interface entre la couche superficielle (20 1 ) et la couche profonde (20 2 ), la couche élémentaire d’interface étant associée à un intervalle temporel d’interface ( ); k) à l’instant de référence et à un instant de mesure, à partir de chaque distribution temporelle résultant de c) établie à l’instant de mesure, calcul : d’une première intensité temporelle ( ), comportant une contribution, à la distribution temporelle, de photons détectés antérieurement à l’intervalle temporel d’interface, lesdits photons étant considérés comme s’étant propagés uniquement dans la couche superficielle ; d’une deuxième intensité temporelle ( ), comportant une contribution, à la distribution de photons détectés postérieurement à l’intervalle temporel d’interface, lesdits photons étant considérés comme s’étant propagés dans la couche superficielle et dans la couche profonde l) calcul d’un premier ratio , à partir des premières intensités temporelles respectivement déterminées à l’instant de référence et à l’instant de mesure et calcul d’un deuxième ratio , à partir des deuxièmes intensités temporelles respectivement déterminées à l’instant de référence et à l’instant de mesure ; m) estimation des variations des concentrations d’oxyhémoglobine et de désoxyhémoglobine ( ) ou de grandeurs proportionnelles auxdites variations ( , ), dans la couche profonde (20 2 ), entre l’instant de mesure et l’instant de référence, à partir des ratios résultant de l) ; Procédé selon la revendication 3, comportant : o) estimation des variations des concentrations d’oxyhémoglobine et de désoxyhémoglobine ( , ), ou de grandeurs proportionnelles auxdites variations ( , ), dans la couche superficielle (20 1 ) entre l’instant de mesure et l’instant de référence, à partir des ratios résultant de m). Procédé selon l’une quelconque des revendications précédentes, dans lequel dans l’étape h), la couche de base est : une couche élémentaire prédéterminée ; ou une couche élémentaire de rang différent de la couche élémentaire considéré, l’écart entre le rang de la couche élémentaire de base et la couche élémentaire considéré étant constant. Procédé selon l’une quelconque des revendications précédentes, dans lequel la source de lumière est une source laser. Dispositif, destiné à détecter une survenue d’une occlusion veineuse ou artérielle dans un tissu biologique (20), entre au moins un instant de mesure ( ) et un instant de référence ( ), le tissu biologique étant délimité par une surface (21), le dispositif comportant : une source de lumière impulsionnelle (10), la source de lumière étant disposée en regard de la surface du tissu biologique, la source de lumière étant configurée pour émettre un faisceau d’illumination (11), formant une zone d’illumination (12), à la surface du tissu biologique, la source de lumière étant configurée pour émettre simultanément ou successivement le faisceau d’illumination dans deux longueurs d’onde d’illumination différentes l’une de l’autre, chaque longueur d’onde s’étendant dans une bande spectrale d’absorption de l’oxyhémoglobine et/ou de la désoxyhémoglobine ; un photodétecteur (16), configuré pour détecter des photons rétrodiffusés par le tissu biologique, après s’être propagés dans le tissu biologique, les photons rétrodiffusés émanant d’une zone de détection (14) à la surface du tissu biologique ; une carte de comptage (17), destinée à établir une distribution temporelle ( ), représentant un nombre de photons détectés par le photodétecteur en fonction du temps, à chaque longueur d’onde ; une unité de traitement (18), programmée pour mettre en œuvre les étapes d) à i) d’un procédé selon l’une quelconque des revendications précédentes, à partir des distributions temporelles établies par la carte de comptage à chaque longueur d’onde, à l’instant de mesure ou à chaque instant de mesure.