RELFEXION CALIBRE

Description

Domaine technique

L'invention se situe dans le domaine de la caractérisation des échantillons, et notamment des échantillons biologiques, et plus particulièrement la peau.

Description de l'art antérieur

Les mesures optiques, permettant de caractériser les propriétés optiques d'échantillons, sont répandues. On peut notamment citer les mesures basées sur la détection d'un signal rétrodiffusé par un échantillon illuminé par un faisceau lumineux. Il s'agit en particulier de la spectrométrie aman, de l'imagerie de fluorescence ou de la spectrométrie de réflectance diffuse.

La spectrométrie de réflectance diffuse, souvent désignée par l'acronyme anglais DRS (Diffuse Réflectance Spectroscopy), consiste à exploiter la lumière rétrodiffusée par un objet diffusant soumis à un éclairement, généralement ponctuel. Cette technique s'avère performante pour caractériser des propriétés optiques superficielles d'échantillons, en particulier les propriétés de diffusion ou d'absorption.

Mise en œuvre sur la peau, cette technique permet par exemple de caractériser le derme ou l'épiderme, comme décrit dans le document EP 2762064. Ce document décrit une sonde de mesure destinée à être appliquée contre la peau. Cette sonde comporte une fibre optique centrale, dite fibre d'émission, reliée à une source de lumière, et destinée à diriger un faisceau lumineux vers un échantillon de peau. Des fibres optiques, disposées autour de la fibre centrale, dites fibres de détection, recueillent un signal optique rétrodiffusé par le derme, ce signal optique étant ensuite détecté par un photodétecteur. Des moyens d'analyse spectrale du signal optique, couplés à des algorithmes de calcul, permettent d'estimer des paramètres du derme, en particulier la concentration de certains chromophores, par exemple l'oxyhémoglobine ou la déoxyhémoglobine et/ou des propriétés optiques gouvernant le parcours des photons dans le derme, notamment le coefficient de diffusion réduit μ ₃ ' ainsi que le coefficient d'absorption μ ₃ .

Ainsi, la sonde comporte une ligne d'illumination, destinée à illuminer l'échantillon, comportant la source de lumière et la fibre optique d'émission. La sonde comporte également une ligne de détection, destinée à détecter une lumière rétrodiffusée par l'échantillon, comportant les fibres optiques de détection et le photodétecteur. Les propriétés de la ligne d'illumination et de la ligne de détection sont prises en compte grâce à une étape de calibration, permettant la détermination d'un facteur de calibration. Ce dernier est obtenu en effectuant une mesure sur un échantillon de calibration, dont les propriétés optiques sont connues. Ce facteur de calibration, désigné par le terme M ^std dans cette demande, est ensuite appliqué au signal mesuré par le photodétecteur, désigné par le terme M ^skin .

Le document Qin J, "Hyperspectral diffuse réflectance imaging for rapid, non contact measurement of the optical properties of turbid materials"Applied Optics Vol 45 N°32, 10 november 2006, décrit un procédé de détermination de propriétés optiques d'un échantillon par spectrométrie de réflectance diffuse. Ce procédé comporte une détection, par un photodétecteur spectrométrique, d'un rayonnement rétrodiffusé par un échantillon à plusieurs distances de rétrodiffusion. Le signal ainsi détecté est multiplié par un facteur de calibration. L'estimation des propriétés optiques est réalisée par un ajustement à l'aide d'un modèle de diffusion représentant l'évolution du coefficient de diffusion réduit dans l'échantillon de calibration. Ainsi, la détermination des propriétés optiques se base sur une connaissance a priori de l'échantillon et d'un modèle de diffusion de la lumière dans l'échantillon analysé, ce modèle décrivant l'évolution de la valeur du coefficient de diffusion en fonction de la longueur d'onde. La détermination des propriétés optiques n'est quantitative que pour un échantillon analysé dont le coefficient de diffusion réduit le même modèle que les échantillons de calibration. On comprend que la nécessité de se baser sur un modèle constitue une limitation contraignante. Un tel procédé ne convient pas à un échantillon complexe, dont le modèle de diffusion a priori n'est pas connu. Par ailleurs, la prise en compte de ce modèle entraîne que différentes estimations d'une propriété optique, à différentes longueurs d'onde, ne sont pas indépendantes les unes des autres, puisque liées par le modèle.

Les inventeurs ont constaté que les procédés précédemment cités ne sont pas optimaux. Un objectif de la présente invention est d'améliorer les procédés de l'art antérieur, de façon à déterminer les propriétés optiques d'un échantillon avec une précision accrue.

Exposé de l'invention

Un objet de l'invention est un procédé de détermination d'une propriété optique d'un échantillon, comportant les étapes suivantes :

i) illumination d'une surface de l'échantillon, à l'aide d'un faisceau lumineux produit par une source de lumière, de façon à former, à la surface dudit échantillon, une zone élémentaire d'illumination, correspondant à la partie de ladite surface éclairée par ledit faisceau ;

ii) acquisition, à l'aide d'un photodétecteur, d'un signal de rétrodiffusion, représentatif d'un rayonnement rétrodiffusé, par l'échantillon, à une distance, dite distance de rétrodiffusion, de ladite zone élémentaire d'illumination ;

iii) sélection d'un facteur de calibration ;

iv) application dudit facteur de calibration à chaque signal de rétrodiffusion, de façon à obtenir une grandeur d'intérêt, associée à ladite distance de rétrodiffusion ;

v) répétition des étapes iv) à v), en mettant à jour le facteur de calibration, en fonction de ladite propriété optique ainsi déterminée, jusqu'à l'atteinte d'un critère d'arrêt ou un nombre d'itérations prédéterminé ;

Les étapes iv) à vi) peuvent être répétées jusqu'à l'atteinte d'un critère d'arrêt ou d'un nombre d'itérations prédéterminé.

Selon un mode de réalisation, au moins un facteur de calibration est un facteur de calibration mesuré, en effectuant un ratio entre :

^■ une estimation d'une grandeur représentative d'un rayonnement rétrodiffusé émanant d'une surface d'un échantillon de calibration, à une distance de rétrodiffusion d'une zone élémentaire d'illumination dudit échantillon de calibration, lorsque l'échantillon de calibration est illuminé par ledit faisceau lumineux ;

■ une mesure de ladite grandeur, à l'aide d'un signal de rétrodiffusion détecté par ledit photodétecteur, l'échantillon de calibration étant illuminé par ledit faisceau lumineux.

Un facteur de calibration peut être déterminé par interpolation à partir de deux facteurs de calibration mesurés, lesdits facteurs de calibration mesurés étant obtenus en utilisant respectivement deux échantillons de calibration dont les propriétés optiques sont différentes. Selon un mode de réalisation, lors de l'étape vi), le facteur de calibration peut être remplacé par un facteur de calibration déterminé en fonction de la propriété optique définie lors de l'étape v) précédent ladite étape vi). Cette propriété optique, considérée pour la mise à jour du facteur de calibration, peut notamment être une propriété optique de diffusion. Il peut par exemple s'agir d'un coefficient de diffusion ou d'un coefficient de diffusion réduit.

Le procédé peut comporter au moins une des caractéristiques suivantes, prises isolément ou selon les combinaisons techniquement réalisables : ^■ L'étape iv) comporte l'application d'un facteur de réactualisation, correspondant à la distance de rétrodiffusion, à chaque facteur de calibration, le facteur de réactualisation ayant été préalablement déterminé, en :

mesurant, à un instant t, un signal de rétrodiffusion, représentant un rayonnement rétrodiffusé émanant de la surface d'un échantillon de calibration, à une distance de rétrodiffusion d'une zone élémentaire d'illumination dudit échantillon de calibration, l'échantillon de calibration étant illuminé par ledit faisceau lumineux ;

comparant ledit signal de rétrodiffusion mesuré à l'instant t à un signal mesuré, dans des mêmes conditions, à un instant to, antérieur à l'instant t.

^■ Lors de l'étape v) la détermination de ladite propriété optique comporte une comparaison entre :

une grandeur d'intérêt,

une pluralité d'estimations de ladite grandeur d'intérêt, chaque estimation étant réalisée en considérant une valeur prédéterminée de ladite propriété optique.

^■ Chaque signal de rétrodiffusion est acquis à une pluralité de longueurs d'onde, de telle sorte que la grandeur d'intérêt et le facteur de calibration peuvent prendre la forme de fonctions spectrales, définies sur ladite pluralité de longueurs d'onde.

■ Ladite grandeur d'intérêt, associée à une distance de rétrodiffusion, est obtenue par l'application d'un ratio entre l'intensité d'un signal de rétrodiffusion, correspondant à ladite distance de rétrodiffusion, par l'intensité dudit faisceau lumineux, mesurée par le photodétecteur, auquel cas ladite grandeur d'intérêt est une réflectance.

^■ Lors de l'étape v), on considère différentes valeurs de ladite grandeur d'intérêt, chaque valeur correspondant à une distance de rétrodiffusion différente.

^■ L'échantillon examiné est un tissu humain, animal, végétal, ou une denrée alimentaire. Un autre objet de l'invention est un support d'enregistrement d'informations, lisible par un processeur, comportant des instructions pour l'exécution d'un procédé décrit ci-dessus, ces instructions étant aptes à être exécutées par un processeur.

Un autre objet de l'invention est un dispositif de mesure d'un signal optique produit par un échantillon comportant :

une source de lumière apte à émettre un faisceau lumineux vers une surface dudit échantillon, de façon à former, sur ladite surface, une zone élémentaire d'illumination un photodétecteur, apte à acquérir, un signal de rétrodiffusion, représentatif d'un rayonnement rétrodiffusé par l'échantillon à une distance, dite distance de rétrodiffusion, de ladite zone élémentaire d'illumination ;

le dispositif étant caractérisé en ce qu'il comporte également :

- un processeur, apte à mettre en œuvre du procédé décrit ci-dessus.

Ce dispositif peut notamment comporter un système optique, configuré pour assurer un couplage optique entre le photodétecteur et une zone élémentaire de détection située à la surface de l'échantillon, de laquelle émane ledit rayonnement rétrodiffusé.

Figures

La figure 1 représente un dispositif permettant l'application de l'invention.

La figure 2 est une vue en coupe de ce dispositif, selon un plan perpendiculaire à l'axe Z et passant par l'extrémité distale des fibres de détection.

La figure 3A représente une configuration de mesure dite « à distance », selon laquelle chaque fibre optique de détection est placée à distance de l'échantillon analysé.

La figure 3B représente une configuration de mesure dite « au contact », selon laquelle chaque fibre optique de détection est placée au contact de l'échantillon analysé.

La figure 4A représente les principales étapes d'un procédé selon l'invention.

La figure 4B représente les principales étapes d'une variante du procédé représenté sur la figure 4A.

La figure 5 représente différents facteurs de calibration pour 3 longueurs d'onde (λ = 470 nm, λ = 607 nm et λ =741 nm), et pour trois échantillons de calibration différents, chaque facteur de calibration ayant un coefficient de diffusion réduit μ différent l'un de l'autre.

La figure 6 représente une modélisation de la réflectance, à une distance de rétrodiffusion donnée, pour différentes valeurs du coefficient d'absorption μ ₃ et du coefficient de diffusion réduit

Les figures 7A, 7B, 7C et 7D représentent les résultats d'essais expérimentaux, montrant l'influence du choix du facteur de calibration sur l'estimation du coefficient de diffusion réduit (figures 7A et 7C) ou du coefficient d'absorption (figure 7B et figure 7D).

Les figures 8A, 8B, 8C et 8D représentent les résultats d'essais expérimentaux comparatifs, montrant l'influence de la mise en œuvre d'un procédé selon l'invention sur l'estimation du coefficient de diffusion réduit (figures 8A et 8C) ou du coefficient d'absorption (figure 8B et figure 8D), selon une configuration de mesure dite au contact, schématisée sur la figure 3B. Les figures 9A, 9B, 9C et 9D représentent les résultats d'essais expérimentaux comparatifs, montrant l'influence de la mise en œuvre d'un procédé selon l'invention sur l'estimation du coefficient de diffusion réduit (figures 9A et 9C) ou du coefficient d'absorption (figure 9B et figure 9D), selon une configuration de mesure dite à distance, schématisée sur la figure 3A.

Exposé de modes de réalisation particuliers

La figure 1 représente un premier mode de réalisation d'un dispositif 1 selon l'invention. Il comporte une source de lumière 10, qui, dans cet exemple, est une source de lumière blanche commercialisée par Océan Optics sous la référence HL2000. La source de lumière 10 comprend, dans cet exemple, une fibre optique d'émission 12, s' étendant entre une extrémité proximale 14 et une extrémité distale 16. La fibre optique d'émission 12 est apte à collecter la lumière par une extrémité proximale 14 et à émettre un faisceau lumineux 20 vers l'échantillon par une extrémité distale 16, ledit faisceau lumineux étant alors dirigé vers la surface d'un échantillon 50. Dans une telle configuration, la source de lumière 10 est dite fibrée.

Le diamètre de la fibre optique d'émission 12 est compris entre 100 μιη et 1 mm, et est par exemple égal à 500 μιη.

Le dispositif comporte également une pluralité de fibres optiques de détection 22i, 222, 22 ₃ ,...22 _f ...22 _F , l'indice f étant compris entre 1 et F, F désignant le nombre de fibres optiques de détection dans le dispositif. F est un entier naturel compris généralement compris entre 1 et 100, et préférentiellement compris entre 5 et 50. Chaque fibre de détection 22i, 22 ₂ , 22 ₃ ,.. 22f,...22 _F s'étend entre une extrémité proximale 24i, 242, 243 ,.. 24f...24 _F et une extrémité distale 26i, 262, ...26 _f ...26 _F . Sur la figure 1, les références 22, 24 et 26 désignent respectivement l'ensemble des fibres de détection, l'ensemble des extrémités proximales des fibres de détection et l'ensemble des extrémités distales des fibres de détection. Le diamètre de chaque fibre optique de détection 22 est compris entre 50 μιη et 1 mm, et est par exemple égal à 100 μιη. L'extrémité proximale 24 de chaque fibre optique de détection 22 est apte à être optiquement couplée à un photodétecteur 40. L'extrémité distale 26i, 26 ₂ ,...26 _F de chaque fibre optique de détection 22 est apte à collecter respectivement un rayonnement 52i, 52 ₂ , ...52 _F rétrodiffusé par l'échantillon 50, lorsque ce dernier est exposé au faisceau lumineux 20.

Le photodétecteur 40 est apte à détecter chaque rayonnement rétrodiffusé 52i, 522, ...52 _F de façon à former un signal, dit signal de rétrodiffusion (Si,S2, ...SN) comme décrit par la suite. Le photodétecteur 40 peut-être un photodétecteur non résolu spectralement, par exemple une photodiode ou un photodétecteur matriciel de type CCD ou CMOS. Dans cet exemple, il s'agit d'un spectrophotomètre, apte à établir le spectre en longueur d'onde du rayonnement collecté par une fibre optique de détection 22 auquel il est couplé. L'homme du métier choisira un photodétecteur spectrométrique lorsqu'il privilégie une bonne résolution spectrale ou un photodétecteur matriciel lorsqu'il privilégie une bonne résolution spatiale.

Le photodétecteur 40 est apte à être connecté à un processeur 48, ce dernier étant relié à une mémoire 49 comportant des instructions, ces dernières pouvant être exécutées par le processeur 48, pour mettre en œuvre le procédé représenté sur les figures 4A ou 4B, et décrit ci-après. Ces instructions peuvent être sauvegardées sur un support d'enregistrement, lisible par un processeur, de type disque dur, CD OM ou autre type de mémoire.

Les fibres optiques de détection 22 s'étendent parallèlement les unes aux autres, parallèlement à un axe longitudinal Z autour de la fibre optique d'émission 12. Elles sont maintenues fixes les unes par rapport aux autres par un élément de maintien 42. Leurs extrémités distales 26 sont coplanaires, et définissent, dans cet exemple, un plan de détection 44.

Une fibre optique 13, dite fibre de retour d'excitation, relie la source de lumière 10 au photodétecteur 40. Cette fibre optique est utile pour réaliser une mesure S _SO urce représentant l'intensité de la source, détaillée ultérieurement.

La figure 2 représente une vue en coupe du dispositif, dans le plan de détection 44, formé par l'ensemble des extrémités distales 26 des F fibres de détection. Dans cet exemple, F est égal à 36. Comme on peut le voir, les fibres optiques de détection sont réparties selon :

un premier groupe Gi de six fibres optiques de détection 21 _\ ...21 _& , disposées régulièrement le long d'un cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 26i...26 ₆ de chaque fibre de ce groupe est distante de l'extrémité distale 16 de la fibre optique d'émission 12 d'une première distance di égale à 300 μιη ;

un deuxième groupe G ₂ de six fibres optiques de détection 22 ₇ ...22i ₂ , disposées régulièrement le long d'un cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 26 ₇ ...26i ₂ de chaque fibre de ce groupe est distante de l'extrémité distale 16 de la fibre optique d'émission 12 d'une deuxième distance d ₂ égale à 700 μιη ; un troisième groupe G3 de six fibres optiques de détection 22i3...22i8, disposées régulièrement le long d'un cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 26i ₃ ...26i ₈ de chaque fibre de ce groupe est distante de l'extrémité distale 16 de la fibre optique d'émission 12 d'une troisième distance d ₃ égale à 1,1 mm ;

un quatrième groupe G4 de six fibres optiques de détection 22ig...2224, disposées régulièrement le long d'un cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 2619...2624 de chaque fibre de ce groupe est distante de l'extrémité d istale 16 de la fibre optique d'émission 12 d'une quatrième distance d4 égale à 1,5 mm ;

un cinquième groupe G ₅ de six fibres optiques de détection 2225...2230, disposées régulièrement le long d'un cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 2625...2630 de chaque fibre de ce groupe est distante de l'extrémité distale 16 de la fibre optique d'émission 12 d'une cinquième distance d ₅ égale à 2 mm ;

un sixième groupe G ₆ de six fibres optiques de détection 22 ₃ i...22 ₃₆ , disposées régulièrement le long d'u n cercle centré sur la fibre optique d'émission 12, de telle sorte que l'extrémité distale 2631...2636 de chaque fibre de ce groupe est distante de l'extrémité distale 16 de la fibre optique d'émission 12 d'une sixième distance d6 égale à 2,5 μιη.

Lorsqu'on évoque une distance entre deux fibres, ou entre une fibre ou un faisceau lumineux, on entend une distance centre à centre.

Ainsi, chaque extrémité distale 26f d'une fibre optique de détection 22f est placée, dans un plan perpendiculaire à l'axe longitudinal Z selon lequel s'étendent ces fibres, à une distance df de la source de lumière 10 (c'est-à-dire de l'extrémité distale 16 de la fibre d'émission 12), et, par conséquent, à une distance df du faisceau lumineux 20 d irigé vers l'échantillon 50.

Selon une variante, les extrémités distales 26f de chaque fibre optique de détection définissent une su rface courbe, s'adaptant par exemple à la courbure de la surface de l'échantillon 50.

Comme précédemment indiqué, le d ispositif comporte un photodétecteur 40, apte à être couplé à l'extrémité proximale 24 _f de chaque fibre optique de détection 22 _f . Dans cet exemple, le photodétecteur est un spectrophotomètre, apte à déterminer le spectre d'un rayonnement 52i...52 _F rétrodiffusé par l'échantillon lorsque ce dernier est exposé au faisceau lumineux 20. Pour cela, les extrémités proximales 24 de chaque groupe de fibres optiques de détection, décrits ci-dessus, sont regroupées et sont, groupe par groupe, successivement couplées au photodétecteur 40 au moyen d'un commutateur optique 41. Sur la figure 1, la référence 52 désigne un rayonnement rétrodiffusé par l'échantillon. Dans cet exemple, le dispositif comporte également un système optique 30, présentant un facteur de grandissement G et un axe optique Z'. Dans cet exemple, l'axe optique Z' est confondu avec l'axe longitudinal Z selon lequel s'étendent les fibres optiques de détection, ce qui constitue une configuration préférée.

D'une façon générale, le système optique 30 permet de former une image de la surface de l'échantillon 50 sur le plan de détection 44 formé par les extrémités distales 26 de chaque fibre optique de détection 22, avec un facteur de grandissement G donné. Ainsi, chaque extrémité distale 26i, 26 ₂ ...26 _F est respectivement conjuguée à une zone élémentaire de détection 28i, 28 ₂ , 28 _f ...28 _F de la surface de l'échantillon. De cette façon, chaque fibre optique de détection 22i, 222,...22 _F est apte à collecter respectivement un rayonnement élémentaire 52i, 522, 52f...52 _F rétrodiffusé par l'échantillon, chaque rayonnement élémentaire 52i, 52 _2, ...52 _F émanant d'une zone élémentaire de détection 28i, 28 ₂ ...28 _F , à la surface de l'échantillon.

Ainsi, chacune desdites extrémités distales 26i, 26 ₂ ...26 _F peut être située dans un plan focal image du système optique 30, et conjuguée à une zone élémentaire de détection 28i, 28 ₂ ...28 _F située dans le plan focal objet dudit système optique, à la surface de l'échantillon. De même, l'extrémité distale 16 de la fibre d'émission 12 est conjuguée à une zone élémentaire d'illumination 18 à la surface de l'échantillon. D'une façon générale, la zone élémentaire d'illumination constitue le point d'impact du faisceau lumineux 20 à la surface de l'échantillon 50.

D'une façon générale, quel que soit le mode de réalisation, le terme zone élémentaire désigne une zone de forme délimitée à la surface de l'échantillon. Une telle zone élémentaire est de préférence ponctuelle, c'est-à-dire que son diamètre ou sa diagonale sont inférieurs à 1 cm, et de préférence inférieurs à 1 mm, voire inférieure à 500 μιη.

Une zone élémentaire de détection peut également prendre une forme annulaire, centrée sur la zone élémentaire d'illumination, en définissant un anneau, ou un arc d'anneau, circulaire ou polygonal. L'épaisseur de l'anneau est alors de préférence inférieure à 1 cm. Une zone élémentaire de détection 28 peut avoir une forme quelconque, sous réserve que cette zone élémentaire soit délimitée par un contour, et distante d'une zone élémentaire d'illumination 18, cette dernière pouvant également avoir une forme quelconque, mais délimitée et distincte d'une zone élémentaire de détection 28.

Une zone élémentaire d'illumination 18 est traversée par le faisceau lumineux 20, se propageant en direction de l'échantillon 50, tandis qu'une zone élémentaire de détection 28 _f est traversée par un rayonnement rétrodiffusé 52 _f , ce rayonnement étant produit par la rétrodiffusion, dans l'échantillon, du faisceau lumineux 20. Le couplage optique, réalisé par le système optique 30, permet à chaque fibre de détection 22 _f de collecter le rayonnement rétrodiffusé élémentaire 52 _f , ce dernier correspondant au rayonnement rétrodiffusé traversant la zone élémentaire 28 _f . L'élément de maintien 42 peut assurer une liaison rigide entre les fibres optiques de détection 22 et le système optique 30, de façon à maintenir le plan de détection 44, formé par les extrémités distales 26 des fibres optiques de détection à une distance fixe du système optique 30.

En lien avec la figure 3A, si d _f est la distance entre l'extrémité distale 26 _f d'une fibre de détection 22 _f et l'extrémité distale 16 de la fibre d'émission 12, distance calculée dans un plan perpendiculaire à l'axe optique Z', la distance D _f entre la zone élémentaire de détection 28 _f , conjuguée de ladite extrémité distale 26 _f , et la zone élémentaire d'illumination 18, conjuguée de ladite extrémité distale 16, est telle que :

La distance D _f est appelée distance de rétrodiffusion, car elle correspond à la distance, à partir de la zone élémentaire d'illumination 18, à laquelle on collecte les photons rétrodiffusés. Cela correspond à la distance entre la zone élémentaire d'illumination 18 et une zone élémentaire de détection 28f.

Ainsi, comme représenté sur la figure 3A, la présence du système optique 30 permet de placer les fibres de détection 22 à distance de l'échantillon à caractériser, selon une configuration dite « à distance ». Cette distance est typiquement de quelques cm, par exemple entre 1 et 30 cm.

Selon une variante, le dispositif est similaire à celui représenté sur la figure 1, mais il ne comporte pas de système optique 30. Cette variante correspond à une configuration de mesure dite configuration « au contact ». Selon cette variante, les fibres de détection 22 sont, de préférence, appliquées directement au contact de l'échantillon (50), du fait de l'absence de système optique entre l'échantillon (50) et l'extrémité distale de chaque fibre optique de détection 22. Dans ce cas, pour chaque fibre optique de détection, D _f = d _f .

D'une manière générale, chaque fibre optique de détection 22 _f est apte à collecter un rayonnement de rétrodiffusion 52 _f issu d'une zone élémentaire de détection 28 _f , cette dernière étant située à une distance de rétrodiffusion D _f de la zone élémentaire d'illumination 18. Dans cet exemple, du fait de l'arrangement concentrique des fibres de détection autour de la fibre d'illumination, décrit en lien avec la figure 2, le dispositif permet de définir N distances de rétrodiffusion distinctes Di..D _n ...D _N , N étant ici égal à 6. A chaque distance de rétrodiffusion D _n correspond une pluralité de rayonnements rétrodiffusés, provenant de zones élémentaires de détection différentes. Par exemple, à la distance de rétrodiffusion Di correspondent les signaux optiques de rétrodiffusion 52i, 522, 52 ₃ , 52 ₄ , 52 ₅ et 52ε. Comme précédemment décrit, les fibres de détection correspondant à un même groupe, c'est-à-dire à une même distance de rétrodiffusion, sont couplées. De ce fait, chaque rayonnement rétrodiffusé correspondant à une même distance de rétrodiffusion D _n (l≤n<6) est adressé simultanément vers le photodétecteur 40, ce dernier produisant un signal S _n , dit signal de rétrodiffusion, représentatif d'un ou plusieurs rayonnements rétrodiffusés à ladite distance de rétrodiffusion D _n .

Quelle que soit la configuration du dispositif représenté sur la figure 1, au contact ou à distance, ce dispositif comporte :

une ligne d'illumination, apte à illuminer la surface de l'échantillon, cette ligne comprenant la source de lumière 10 et la fibre optique d'émission 12 ;

une ligne de détection, apte à détecter un rayonnement rétrodiffusé 52 _f par l'échantillon de façon à former un signal de rétrodiffusion S _n ; cette ligne comprend le photodétecteur 40 ainsi que l'éventuel système optique 30 lorsque le dispositif en est équipé. On décrit à présent, en lien avec la figure 4A, les principales étapes d'un procédé itératif pouvant être mis en œuvre par le dispositif précédemment décrit, afin d'estimer une ou plusieurs propriétés optiques p de l'échantillon examiné. Ce procédé itératif s'applique à la configuration « au contact » ou à la configuration « à distance » précédemment décrites.

Le terme propriété optique p désigne par exemple un ou plusieurs facteurs gouvernant l'absorption et/ou la diffusion des photons dans l'échantillon examiné, en particulier un coefficient d'absorption, un coefficient de diffusion, un coefficient de diffusion réduit, un coefficient d'anisotropie de la diffusion. Dans cet exemple, les propriétés optiques déterminées sont le coefficient d'absorption μ ₃ et le coefficient de diffusion réduit μ . l ^ère étape 110 : application du dispositif préalablement décrit, en regard de l'échantillon 50.

2 ^ième étape 120 : illumination de l'échantillon en dirigeant un faisceau lumineux 20 contre la surface de l'échantillon, la partie de la surface de l'échantillon illuminée constituant la zone élémentaire d'illumination 18.

g ⁱ eme étape 13Q _: collection d'un rayonnement 52i, 52 ₂ ...52 _F rétrodiffusé par l'échantillon, émanant respectivement de chaque zone élémentaire de détection 28i, 28 ₂ , ... 28 _F , par la fibre optique 22i, 222,...22 _F de détection dont l'extrémité distale 26i, 26 ₂ ,... 26 _F est respectivement conjuguée à ladite zone élémentaire 28i, 28 ₂ , ... 28 _F .

4' ^eme étape 140 : mesure, à l'aide d'un photodétecteur 40, d'un signal de rétrodiffusion S _n représentatif de la rétrodiffusion à chaque distance de rétrodiffusion D _n . Comme précédemment indiqué, le signal détecté S _n est, dans cet exemple, établi en cumulant les signaux optiques collectés par les fibres optiques de détection d'un même groupe G _n , c'est-à-dire correspondant à une même distance de rétrodiffusion. Le signal de rétrodiffusion S _n cumule alors plusieurs rayonnements rétrodiffusés 52 _n , chacun d'entre eux étant émis selon une même distance de rétrodiffusion D _n .

Lorsque le détecteur est un détecteur spectrométrique, il génère le spectre du signal détecté S _n , noté Sp(Sn), duquel on peut extraire des composantes spectrales S _n ( ), représentant le signal rétrodiffusé à la distance D _n , et à la longueur d'onde λ.

g ⁱ eme étape 15Q ; A l'aide de chaque signal S _n ( ), associé à une distance de rétrodiffusion D _n , détermination d'une grandeur d'intérêt R _n ( ), sur la base de laquelle seront déterminées les propriétés optiques p de l'échantillon examiné. Dans cet exemple, la grandeur d'intérêt R _n ( ) est une réflectance de l'échantillon. D'une façon générale, le terme réflectance représente l'intensité d'un rayonnement rétrodiffusé par l'échantillon, normalisée par l'intensité du faisceau incident à l'échantillon. Sa valeur dépend de la longueur d'onde λ, du fait de l'évolution des propriétés optiques du milieu diffusant examiné en fonction de la longueur d'onde.

Dans cet exemple, la réflectance R _n ( ) dépend du signal rétrodiffusé S _n ( ) à la distance D _n , normalisé par une quantité de lumière S _SO urce^) émise par la source, à la longueur d'onde λ, du temps d'acquisition du signal rétrodiffusé S _n , ainsi que d'un facteur de calibration. La réflectance R _n ( ) peut être définie selon l'expression : ou :

S _n (X) est le signal de rétrodiffusion détecté, correspondant à la distance de rétrodiffusion D _n ;

5 - S _re f (X) est un signal de référence, représentatif de signaux parasites, tel le bruit du détecteur 40 ou des réflexions parasites de l'éventuel système optique 30, obtenu en activant la source de lu mière, mais sans échantillon, ce dernier étant par exemple remplacé par un écran a bsorbant de type écran noir ;

S _source (X) est le signal produit par la source de lumière. S _source (X) peut notamment 10 être éta bli en couplant la source de lumière au photodétecteur, par exemple au moyen d'une fibre optique dite de retour d'excitation 13, représentée sur la figure 1 ; Dans ce cas, le photodétecteur acquiert un signal S _source _ _direct , duquel on peut soustraire un signal S _re f_ _source représentatif du bruit du détecteur. Si t _source désigne le temps d'acquisition du signal S _{source→irect} , S _source peut être tel que :

^ _p Ssource-direct( —S _re f- _source (A^

-L- ³ ^source (A) ^— I ·

^L source

fn W ^{est u n} facteur de calibration, correspondant à la distance de rétrodiffusion D _n et à la longueur d'onde λ. L'exposant i désigne le rang de l'itération, tandis que l'indice n désigne la distance de rétrodiffusion D _n . Ce facteur prend en compte l'effet de différents composants de la ligne d'illumination et de la ligne de détection sur le signal de

20 rétrodiffusion. Il s'agit de tenir compte, par exemple, de l'efficacité de collection par les fibres de détection 22, de la sensibilité du photodétecteur 40, de l'inhomogénéité du faisceau d'illumination 20 ou, le cas échéant, de l'efficacité de collection de la lumière rétrod iffusée par le système optique 30. Ce facteur de calibration est déterminé au cours d'une phase de calibration. Cette phase de calibration, mettant en œuvre des

25 échantillons de calibration, est réalisée préalablement ou postérieurement à la mesure sur l'échantillon, et est décrite par la suite.

t est le temps d'acquisition du signal de rétrodiffusion S _n .

La phase de calibration évoquée ci-dessus a pour but d'établir un facteur de calibration fn pcaiib W ^en appliquant le dispositif décrit ci-dessus en regard d'un échantillon de calibration,

30 dont les propriétés optiques p _ca iib sont connues. Par exemple, f _niPca ub peut être tel que :

Jn.pcalib W ^~ s _calib _ _n Ç.)-s _ref(7i) ^~ _{R (λ)} ^Δ Ι

ssource( ) ^xt: calib Scaiib-nQ ^est un signal de rétrodiffusion détecté, correspondant à la distance de rétrodiffusion D _n en utilisant le même dispositif que celui mis en œuvre pour acquérir le signal de rétrodiffusion S _n (X), le dispositif étant utilisé dans la même configuration : même source, même positionnement par rapport à l'échantillon ;

- S _re f (X) est le signal de référence décrit en lien avec l'expression (1) ;

S _source (X) est le signal représentant l'intensité du faisceau d'illumination produit par la source de lumière, décrit en lien avec l'expression (1) ;

- tcaiw est le temps d'acquisition du signal S _calib _ _n (X) ;

R _C aiib-n ^est ' ^a réflectance de l'échantillon de calibration, associée à une distance de rétrodiffusion D„. Dans cet exemple,

^{est u ne} estimation, de la réflectance R _ca iib-nQ , cette estimation pouvant être réalisée par modélisation du parcours de la lumière dans l'échantillon de calibration, en particulier au moyen de codes de calcul de type Monte-Carlo ou par un modèle analytique. Ainsi, le facteur de calibration f _niPC aiib W ^{est u ne} comparaison entre une grandeur d'intérêt, en l'occurrence une réflectance, modélisée et cette même grandeur d'intérêt mesurée par le dispositif, sur un échantillon de calibration. Cette comparaison prend généralement la forme d'un ratio.

Cependant, dans l'art antérieur, ce facteur de calibration est obtenu sur un échantillon, dont les propriétés optiques p _œ iib sont connues, mais ne sont pas forcément représentatives des propriétés optiques de l'échantillon examiné. Or, les inventeurs ont déterminé que la valeur de ce facteur de calibration peut évoluer, selon les propriétés optiques de l'échantillon. Par exemple, la figure 5 représente différentes valeurs de ce facteur de calibration, obtenu en utilisant différents échantillons de calibration, aux longueurs d'onde λ = 470 nm, λ = 607 nm et λ = 741 nm. Pour réaliser ces essais, on a constitué des échantillons d'eau, dont les propriétés de diffusion et d'absorption sont respectivement modifiées par incorporation d'intralipide et d'encre de Chine. Les échantillons de calibration utilisés comportent une concentration d'intralipide %IL respectivement égale à 1%, 2% et 3% ce qui leur confère des propriétés de diffusion différentes, le coefficient d'absorption étant égal à 0,4 cm ¹ à 600 nm ^"1 . Les facteurs de calibration représentés ont été déterminés en considérant une distance de rétrodiffusion de 1.1 mm, le dispositif 1 étant placé à distance de chaque échantillon de calibration, la distance entre l'échantillon et les fibres optiques de détection s'élevant à 20 cm. En se basant sur des essais expérimentaux, décrits en lien avec les figures 7A à 7D, les inventeurs ont estimé qu'il était préférable d'utiliser un facteur de calibration qui soit le plus représentatif possible des propriétés optiques de l'échantillon examiné. Or, lors de la première itération, ces propriétés ne sont pas connues. Aussi, lors de la première itération (i = 1), on utilise un facteur de calibration initial, noté déterminé arbitrairement, par exemple en se basant sur un a priori quant aux propriétés optiques de l'échantillon examiné. gième étape 160 : pour au moins une longueur d'onde λ, et en considérant au moins autant de distances D _n de rétrodiffusion différentes qu'il y a de propriétés optiques à estimer, détermination des propriétés optiques (p) présentant le moindre écart entre la réflectance R _n (X), déterminée lors de l'étape précédente, à la longueur d'onde λ, et une réflectance modélisée Rn ^del (À), cette réflectance étant modélisée en considérant une pluralité de valeurs desdites propriétés optiques p, à ladite distance de rétrodiffusion D _n . Cette détermination peut être réalisée par la minimisation d'un écart quadratique, et par exemple selon l'expression :

p = argmin(∑ _{= 1} (¾ ^odei (A) - R _n (A)) ² ) (3)

p

- N désigne le nombre de distances de rétrodiffusion prises en compte,

Rn,p ^{del est u ne} réflectance modélisée, à la distance de rétrodiffusion D _n , en prenant en compte des valeurs prédéterminées d'au moins une propriété optique p. Le paramètre p peut correspondre à une propriété optique, ou un ensemble de propriétés optiques.

Dans cet exemple, les propriétés optiques recherchées sont μ ₃ (λ) et μ^ίλ). Ainsi, le couple (μ ₃ (λ), μ^λ)) recherché est celui présentant le moindre écart entre la réflectance mesurée R _n (A), à la longueur d'onde λ, et une réflectance modélisée pour différentes valeurs de μ ₃ (λ) et de μ^ίλ), à ladite distance de rétrodiffusion D _n . Cette détermination peut être réalisée selon l'expression

μ _α (λ), μ' ₅ (λ)) = ar _g mm^ , _sW ^ _=i R^^{ ) - R _n (λ)) ² ) (3'), désignant une réflectance modélisée, à la distance de rétrodiffusion D _n , en considérant différentes valeurs de μ ₃ et μ$ .

Les différentes valeurs de réflectance modélisée sont obtenues, pour une pluralité de couples de valeurs μ ₃ , μ$' au cours d'une phase de paramétrage, par simulation numérique mettant en œuvre une méthode de type Monte-Carlo ou par un modèle analytique. Un modèle analytique peut être utilisé, de préférence, qu'au-delà d'une certaine distance de rétrodiffusion. Pour une distance de rétrodiffusion D _n donnée, on peut établir une pluralité de réflectances 6μ _5! (λ) modélisées en fonction de μ ₃ et de μ ₃ '. La figure 6 donne un exemple de représentation de telles réflectances modélisées, en considérant une distance de rétrodiffusion D _n égale à 700 μιη et en prenant en compte des valeurs du coefficient d'absorption comprises entre 0 et 10 cm ^"1 , ainsi que des valeurs du coefficient de diffusion réduit comprises entre 0 et 80 cm ^"1 . Les étapes 150 et 160 sont mises en œuvre par le processeur 48, préalablement programmé à cet effet, et dont les données d'entrée sont les mesures des signaux de rétrodiffusion S _n (A) réalisées par le photodétecteur 40. Chaque facteur de calibration, ainsi que chaque valeur R™ ^del (A), peuvent être stockés dans une mémoire, par exemple la mémoire 49, reliée au processeur 48. ieme étape 170 _{: m} j _se à jour du facteur de calibration.

La mise en œuvre de cette étape suppose que différents facteurs de calibration fn _:P (X) correspondant à des échantillons de calibration de propriétés optiques p connues, à une distance de rétrodiffusion D _n et à une longueur d'onde, λ aient été préalablement déterminés. D'une façon générale, la notation f _niP (X) correspond à un facteur de calibration correspondant aux propriétés optiques p, à la distance de rétrodiffusion D _n , pour la longueur d'onde λ. Ce facteur de calibration peut être obtenu à l'aide d'une mesure sur un échantillon de calibration, auquel cas on peut également le désigner par f _niCa ii _bp W _> l'indice calibp se référant à l'échantillon de calibration de propriétés optiques p. Il peut également être déterminé par calcul d'interpolation, comme décrit ci-après.

Ces différents facteurs de calibration peuvent être obtenus expérimentalement, en utilisant des échantillons de calibration de propriétés optiques p connues, comme décrit en lien avec la figure 5 ou l'équation (2). Lorsque plusieurs mesures expérimentales ont été effectuées, on peut déterminer des facteurs de calibration interpolés, entre deux facteurs de calibrations fn _,p W, f _{n p} ' W correspondant respectivement à des échantillons de propriétés optiques p et p'. L'interpolation peut être une interpolation linéaire.

On peut alors disposer d'une bibliothèque de facteurs de calibration f _np {X) correspondant à différentes distances de rétrodiffusion D _n , à différentes propriétés optiques p et à différentes longueurs d'onde λ. Ces facteurs de calibration sont stockés dans une mémoire, par exemple la mémoire 49 reliée au processeur 48. Les inventeurs estiment qu'il est suffisant qu'entre deux itérations, les facteurs de calibration soient mis à jour en fonction d'une propriété de diffusion de l'échantillon, une mise à jour en fonction d'une propriété d'absorption pouvant être omise.

L'étape 170 consiste à mettre à jour le facteur de calibration mis en œuvre dans le procédé, en remplaçant chaque facteur de calibration fn (À), associé à une distance de rétrodiffusion D _n lors de l'itération courante / ^' par un facteur de calibration correspondant aux propriétés optiques p déterminées lors de l'étape 160, ou par un facteur de calibration associé à une propriété optique la plus proche de la propriété optique p déterminée lors de l'étape 160 précédente. Aussi, lors de l'étape 170, fn ^{+ 1} (X) = f _n, W, le paramètre p étant le paramètre optique déterminé lors de l'étape 160. Ce facteur de calibration fn ^{+ 1} (X) est alors utilisé lors de l'étape 150 de l'itération suivante i+1.

Le processus itératif est stoppé au bout d'un nombre d'itérations prédéterminé, ou lorsque l'écart entre propriétés optiques p ¹ , p ^l+1 déterminées au cours de l'étape 160 de deux itérations successives i et i+1 est inférieur à un seuil prédéterminé. On passe alors à l'étape 180 de sortie de l'algorithme. La figure 4B représente une variante de ce procédé, dans lesquelles les étapes 110 à 180 sont analogues à celles exposées en lien avec la figure 4A. Cependant, préalablement à la mise en œuvre de ce procédé, on place le dispositif 1 face à un échantillon de calibration dont les propriétés optiques p _ca ii _b Sont connues. En effet, les inventeurs ont constaté qu'un facteur de calibration f _niP (À) n'est pas stable dans le temps, cela du fait de l'évolution des propriétés des composants constituant la ligne d'illumination et la ligne de détection. Cela peut provenir d'une évolution normale de ces composants, par exemple l'usure d'une fibre optique, ou le vieillissement de la source, ou encore un léger déplacement du système optique. De ce fait, un facteur de calibration f _n>p>t0 W déterminé à un instant t ₀ peut-être différent d'un facteur de calibration déterminé à un instant t, et cela d'autant plus que l'intervalle temporel ΔΤ = t- to est important. Afin de prendre en compte une telle dérive, les inventeurs ont mis en œuvre une réactualisation des facteurs de calibrations f _n>v>t oW déterminés à un instant to et stockés en mémoire. C'est l'objet des étapes 100 à 106.

Les étapes 100, 101, 102 et 103 sont respectivement analogues aux étapes 110, 120, 130 et 140, la seule différence étant que l'échantillon analysé est l'échantillon de calibration. Les fibres optiques de détection 22 collectent une pluralité de rayonnements rétrodiffusés l'exposant * désignant le fait qu'on utilise un échantillon de calibration. Le photodétecteur 40 forme alors autant de signaux de rétrodiffusion S _œ iib-i(A)...S _œ iib-N(A) qu'il y a de distances de rétrodiffusion DI... DN différentes.

Comme indiqué en lien avec l'équation (2), l'étape 104 permet l'obtention d'une réflectance Rcaiw-nO , à chaque distance de rétrodiffusion D _n et pour chaque longueur d'onde λ considérée. Lors de l'étape 105, on détermine un facteur de calibration f _{n pca} iib* t (λ) correspondant de l'échantillon de calibration, selon l'équation 2.

Lors de l'étape 106, on détermine un facteur de réactualisation k _{n t} Qv), associé à une distance de rétrodiffusion D _n et à une longueur d'onde λ, k _t (X) étant tel que :

fn _p caii _b * _t W désignant le facteu r de calibration, à la longueur d'onde λ, réalisé à l'instant courant (instant t), sur un échantillon de calibration de propriétés optiques Pcaiib*, et correspondant à la distance de rétrodiffusion D _n ;

fn _p caiw* _t o W désigna nt le facteur de calibration, à la longueur d'onde λ, réalisé à un instant t ₀ , préala blement à l'instant courant, sur le même échantillon de calibration, de propriétés optiques p _ca iib*, et correspondant à la distance de rétrodiffusion D _n .

En combinant les équations (4) et (2), on obtient :

s calib -n,to ⁽ ^> ^~S ref,to k) _^

k Q-v _ _ fn.pcalib'.tW _ S _sourcef0 (X)' ^{^ t} callb*,t

n ^,t f (X) ^s calib* -n,t ^v >- ^s ref,t( _i ) ^ * '

fn.pcalib'.toW — Xt _C alib*,to

^• >source,tv )

Les indices t et t ₀ se réfèrent respectivement aux instants de mesure t et t ₀ . L'exposant * représente une mesure réalisée sur un échantillon de calibration utilisé pour déterminer le facteur de réactualisation. Le facteur de réactualisation k _{n t} ( ) est essentiellement gouverné par l'évolution des signaux rétrodiffusés S _ca i _ib *_ _{n t0} (X) et S _ca i _ib *_ _{n t} (λ).

Aussi, d'une façon plus générale, le facteur de réactualisation k _{n t} (X) est déterminé en comparant :

un signal de rétrodiffusion S _caUb * _ _{n t} Q ) , représentant un rayonnement rétrodiffusé émanant de la surface d'un échantillon de calibration, à une distance de rétrodiffusion (D _n ) d'une zone d'illumination dudit échantillon de calibration, ce dernier étant illuminé par ledit faisceau lu mineux (20) ;

et un signal de rétrodiffusion S _caUb *_ _{n tQ} Q), mesuré, dans les mêmes conditions, à un instant to, antérieur à l'instant t. L'échantillon de calibration utilisé pour la détermination du facteur de réactualisation, selon les équations (4) et (4'), peut être n'importe quel échantillon de calibration. De préférence, il s'agit d'un échantillon de calibration aisément transportable, dont les propriétés optiques sont particulièrement stables, en particulier entre les instants t et t ₀ . Il peut par exemple s'agir d'un échantillon réalisé à l'aide d'une résine solide, dont les propriétés optiques d'absorption et de diffusion sont respectivement ajustées par adjonction d'encre de Chine et de particules diffusantes d'oxyde de Titane (T1O2).

Les inventeurs ont estimé qu'un tel facteur de réactualisation peut s'appliquer à l'ensemble des facteurs de calibration préalablement calculés, qu'ils soient issus d'autres échantillons de calibration, moins stables ou moins transportables, ou de calculs d'interpolation. Ainsi, chaque facteur de calibration f _n>v>t o W, après avoir été déterminé à un instant to, antérieur à l'instant t, et stocké dans la mémoire 49, peut-être simplement réactualisé par la formule de mise à jour :

fn.p.tW = kn.tOO ^x fn.p.toW ( ⁵ )

où :

- f _n ,p,tW désigne le facteur de calibration, correspondant aux propriétés optiques p, et à la distance de rétrodiffusion D _n , réactualisé à l'instant courant t ;

fn.p.to W désigne le facteur de calibration, correspondant aux propriétés optiques p, et à la distance de rétrodiffusion D _n , déterminé à l'instant to et stocké dans la mémoire 49. On remarque qu'un seul échantillon de calibration peut suffire pour déterminer le facteur de réactualisation k _{n t} (À), et permettre la réactualisation de l'ensemble des facteurs de calibration fn.p.to W établis auparavant, correspondant à la distance de rétrodiffusion D _n à laquelle est associé le facteur de réactualisation, et stockés dans la mémoire 49.

Le facteur de réactualisation k _{n t} (À) est ensuite mis en œuvre, sous forme de terme multiplicatif, dans l'étape 150 de détermination d'une grandeur d'intérêt R _n ( ) à partir du signal rétrodiffusé S _n ^). L'expression (1) peut alors être remplacée par l'expression ( ) :

* (λ) = x ti W x k _n>t (A) (ΐ')

On va à présent décrire des essais expérimentaux mettant en œuvre le dispositif représenté sur les figures 1 et 2, dans une configuration à distance ou dans une configuration au contact. Le dispositif est disposé en regard d'échantillons tests comportant une base d'eau, et dont les propriétés optiques d'absorption et de diffusion sont ajustées respectivement par ajout d'encre de Chine et d'Intralipide. Durant ces essais, les propriétés optiques p recherchées sont le coefficient d'absorption μ ₃ et le coefficient de diffusion réduit

Les figures 7A et 7B représentent des estimations respectives du coefficient de diffusion réduit et du coefficient d'absorption en fonction de la longueur d'onde, lors d'essais réalisés au contact sur un échantillon de test dont les propriétés optiques sont connues : son coefficient d'absorption est égal à 1 cm ¹ à 600nm, tandis que son coefficient de diffusion réduit, équivalent à une concentration de 1,5% d'intralipide, s'élève à 22 cm ¹ à 600 nm.

Sur la figure 7 A, chaque courbe en pointillés correspond à la valeur théorique du coefficient de diffusion réduit de 4 échantillons de calibration, en fonction de la longueur d'onde. Chaque échantillon de calibration a un même coefficient d'absorption μ ₃ = 0.4 cm ¹ à 600 nm et une concentration d'intralipide respectivement égale à 1%, 1.5%, 2% et 3%. Leurs coefficients de diffusion réduits μ$ , à 600 nm, sont respectivement 13,5 cm , 20,3 cm ^"1 , 27 cm ¹ et 40,6 cm ^"1 . Ces quatre échantillons de calibration sont respectivement désignés par « CF1% », « CF1.5% », « CF2% » et « CF3% ». Ces échantillons de calibration sont utilisés pour établir un facteur de calibration f _n>v ( ), associé aux propriétés optiques p de chaque échantillon.

L'échantillon test a fait l'objet d'une illumination par la source de lumière 10, durant laquelle on a détecté le signal rétrodiffusé $2(λ)...$6(λ) correspondant respectivement à 5 distances de rétrodiffusion D2... D6. Le spectre en longueur d'onde de chacun de ces signaux détectés a été réalisé, dans une bande spectrale comprise entre 470 nm et 880 nm. La réflectance de l'échantillon test aux différentes distances de rétrodiffusion a été déterminée en utilisant l'expression (1).

A chaque longueur d'onde λ, on a successivement considéré chaque facteur de calibration f _n . _p W, associé à chaque échantillon de calibration, de façon à calculer 4 mesures de la réflectance. On a alors estimé le coefficient de diffusion réduit μ ₅ '(λ) (cf. courbes en traits pleins de la figure 7A) et le coefficient d'absorption μ ₃ (λ) (cf. courbes en traits pleins de la figure 7B), chaque estimation étant respectivement associée à la prise en compte d'un facteur de calibration établi à l'aide d'un échantillon de calibration, comme indiqué dans la légende de ces courbes.

Sur la figure 7A, la courbe en pointillés désignée <^ _s - _t est » correspond au coefficient de diffusion réduit μ ₅ (λ) de l'échantillon de test. Elle représente la valeur exacte du coefficient de diffusion réduit en fonction de la longueur d'onde. Les courbes en traits pleins correspondent aux estimations de ce coefficient, à différentes longueurs d'onde. Il apparaît que la prise en compte d'un facteur de calibration basé sur l'échantillon de calibration CF3% aboutit à une estimation erronée de μ ₅ (λ). Les estimations utilisant un facteur de calibration établi avec les autres échantillons de calibration (CF1.5%, CF1% et CF2%) sont plus conformes à la réalité, la meilleure estimation étant obtenue avec le facteur de calibration établi avec l'échantillon de calibration CF1.5%. Cela confirme l'hypothèse de base de cette invention, selon laquelle les propriétés optiques d'un échantillon sont d'autant mieux estimées qu'elles sont calculées en se basant sur un facteur de calibration f _n>v (X) représentatif des propriétés optiques dudit échantillon.

Sur la figure 7B, la courbe en pointillés désignée <^ _a - _t es _t » correspond au coefficient d'absorption μ ₃ (λ) de l'échantillon de test. Elle représente la valeur exacte du coefficient d'absorption en fonction de la longueur d'onde. Les courbes en traits pleins correspondent à des estimations de ce coefficient, à différentes longueurs d'onde, chaque estimation étant réalisée en considérant un facteur de calibration f _niP ( ) déterminé respectivement avec chaque échantillon de calibration. De même que sur la figure 7 A, les estimations basées sur l'échantillon CF3% mènent à des résultats erronés, la meilleure estimation étant celle prenant en compte l'échantillon de calibration CF1.5%.

Les figures 7C et 7D représentent respectivement des résultats similaires aux figures 7A et 7B, le dispositif utilisé comportant un système optique 30, permettant son utilisation à distance de l'échantillon. Dans cette configuration, l'extrémité distale des fibres optiques de détection est placée à 20 cm de la surface de l'échantillon. On constate que la prise en compte d'un facteur de calibration représentatif des propriétés optiques de l'échantillon examiné améliore très significativement l'estimation du coefficient d'absorption. Plus précisément, on remarque que la prise en compte d'un facteur de calibration basé sur des propriétés optiques différentes de l'échantillon examiné mène à des erreurs significatives dans l'estimation du coefficient d'absorption, comme le montre la courbe CF3% de la figure 7D.

Les figures 8A, 8B et 8C, 8D représentent des estimations des propriétés optiques (μ _ϊ ( ) et μ ₃ (λ)) respectivement selon l'art antérieur et en mettant en œuvre l'invention, le dispositif étant appliqué au contact de quatre échantillons tests. Les valeurs réelles du coefficient de diffusion réduit de chaque échantillon de test sont représentées en pointillés sur les figures 8A et 8C. Les valeurs réelles du coefficient d'absorption de chaque échantillon de test sont représentées en pointillés sur les figures 8B et 8D. Ces 4 échantillons tests, désignés IL1%, IL1.5%, IL2%, IL3% sont respectivement identiques aux échantillons de calibration CF1%, CF1.5%, CF2% et CF3% précédemment décrits. Sur chaque figure, les courbes en traits pleins correspondent à des estimations des coefficients μ ₅ (λ) ou μ ₃ (λ) de chaque échantillon de test.

La figure 8A représente des estimations du coefficient de diffusion réduit de chaque échantillon test. Sur cette figure, pour chaque estimation, on a utilisé le même facteur de calibration, établi à l'aide de l'échantillon de calibration CF1%. Le coefficient de diffusion réduit est correctement estimé pour l'échantillon IL1% , puisque le procédé utilise un facteur de calibration établi avec ce même échantillon. Le coefficient de diffusion réduit de l'échantillon ILl.5% est également déterminé correctement. En revanche, les coefficients de diffusion réduits des échantillons IL2% et IL3% ne sont pas estimés avec une précision satisfaisante. La figure 8C représente des mesures similaires, en mettant en œuvre l'algorithme précédemment décrit, avec, lors de la première itération, l'utilisation d'un facteur de calibration établi à l'aide de l'échantillon de calibration IL1%. Contrairement aux résultats présentés sur la figure 8A, le coefficient de diffusion réduit de chaque échantillon de test a été correctement estimé. La figure 8B représente des estimations du coefficient d'absorption de chaque échantillon test. Sur cette figure, pour chaque échantillon test, on a utilisé le même facteur de calibration, établi à l'aide de l'échantillon de calibration CF1%. La figure 8D représente des mesures similaires, en mettant en œuvre l'algorithme précédemment décrit, avec, lors de la première itération, l'utilisation d'un facteur de calibration établi à l'aide de l'échantillon de calibration CF1%. La précision de l'estimation est satisfaisante dans les deux cas, mais la mise en œuvre de l'algorithme accroît cette précision.

Pour chacune de ces figures, on a estimé la racine carrée de l'erreur quadratique moyenne, notée ε, normalisée, estimée selon l'expression :

avec :

p X) = valeur réelle de la propriété optique p à la longueur d'onde λ, la propriété optique étant soit le coefficient de diffusion réduit μ soit le coefficient d'absorption μ ₃ . p X) = estimation de la propriété optique p à la longueur d'onde λ, la propriété optique étant soit le coefficient de diffusion réduit μ soit le coefficient d'absorption μ ₃ .

Les résultats correspondant aux différentes figures 8A (estimation de sans mise en œuvre de l'invention), 8B (estimation de μ ₃ sans mise en œuvre de l'invention), 8C (estimation de μs' avec mise en œuvre de l'invention), et 8D (estimation de μ ₃ 3νθε mise en œuvre de l'invention), sont reportés sur le tableau 1 ci-dessous :

Tableau 1

Les figures 9A, 9B, 9C et 9D représentent des essais similaires à ceux reportés respectivement reportés sur les figures 8A, 8B, 8C et 8D, la seule différence étant que le dispositif est utilisé selon une configuration « à distance », et en mettant en œuvre un système optique de focalisation, la distance entre l'extrémité de chaque fibre de détection et l'échantillon étant de 20cm.

Pour chacune de ces figures, on a également estimé la racine carrée de l'erreur quadratique moyenne, ε, telle que précédemment définie. Les résultats sont reportés sur le tableau 2 ci- dessous.

Tableau 2

La mise en œuvre d'un algorithme selon l'invention permet d'améliorer significativement la précision des estimations de propriétés optiques de l'échantillon. Bien que les essais décrits aient été réalisés en mettant en œuvre une source de lumière blanche et un photodétecteur 40 spectrométrique, des configurations basées sur une source de lumière monochromatique, ou une pluralité de sources de lumière émettant dans des bandes spectrales différentes, et/ou la détection d'un signal de rétrodiffusion à l'aide d'un photodétecteur non spectrométrique sont envisageables. En particulier, la source de lumière blanche peut être remplacée par différentes sources de lumière émettant dans des bandes spectrales différentes λι, λ2...λι_. Ainsi, le faisceau d'illumination 20 peut comporter, simultanément ou successivement, différentes bandes spectrales λι, λ2...λι_. Le dispositif peut également comporter une source de lumière, comprenant une pluralité de filtres optiques passe-bande, aptes à être successivement placés en regard de la source. De cette façon, le faisceau d'illumination 20 comporte successivement, différentes bandes spectrales λι, λ2...λι_.

D'une façon générale, la source de lumière, quelle qu'elle soit, est apte à former, à la surface de l'échantillon, une zone élémentaire d'illumination telle que précédemment définie. Le recours à des fibres optiques pour former le faisceau d'illumination n'est pas nécessaire. Une source de lumière pourrait être une source Laser, ou une autre source de lumière, par exemple une diode électroluminescente. La source de lumière peut être couplée à un système optique de mise en forme, permettant la formation du faisceau lumineux 20 et sa projection sur la surface de l'échantillon afin de définir la zone d'élémentaire d'illumination 18. De même, le photodétecteur peut être une photodiode, ou un photodétecteur matriciel de type CCD ou CMOS. Chaque pixel du photodétecteur est alors couplé à une zone élémentaire de détection soit par le système optique 30, soit en étant placé au contact de la surface de l'échantillon, ou éventuellement par l'intermédiaire de fibres optiques. L'utilisation d'un tel photodétecteur permet d'obtenir un grand nombre de distances de rétrodiffusion différentes. Elle est à privilégier dans les applications requérant une bonne résolution spatiale. Lorsque la source de lumière est apte à former un faisceau d'illumination 20, successivement, dans différentes bandes spectrales λι, λ ₂ ...λι_, un tel photodétecteur peut détecter mesurer un signal de rétrodiffusion S _n ( ), successivement, dans chacune des bandes spectrales. De préférence, la propriété optique est alors déterminée dans chaque bande spectrale, indépendamment l'une de l'autre, en mettant en œuvre les étapes décrites ci-avant. Comme on peut le voir en lien avec les exemples décrits ci-avant, la largeur d'une bande spectrale peut être inférieure à 10 nm, de façon à avoir une estimation précise de l'évolution de la propriété optique considérée en fonction de la longueur d'onde.

Le nombre N de distances de rétrodiffusion peut également varier. D'une façon générale, ce nombre doit être supérieur ou égal au nombre de propriétés optiques à déterminer.

L'invention peut être mise en œuvre pour caractériser les propriétés optiques superficielles d'un échantillon. Appliquée à la peau, elle permet par exemple la détection précoce de pathologies, le contrôle de la vascularisation ou de la perfusion d'un principe actif. Elle peut être appliquée à toute application, de type contrôle non destructif, permettant d'estimer ou de suivre l'évolution d'une propriété optique à proximité de la surface d'un échantillon. Il peut par exemple s'agir d'applications dans le domaine de l'agroalimentaire, afin de contrôler la qualité ou la composition d'aliments.