FR 3050528 A1 20171027 FR 1653587 A 20160422 1 Procédé et dispositif d'estimation de propriétés optiques d'un échantillon Description DOMAINE TECHNIQUE L'invention se situe dans le domaine de la caractérisation de propriétés optiques d'un échantillon, par exemple un échantillon de peau. ART ANTERIEUR Les mesures optiques, permettant de caractériser des propriétés optiques d'échantillons, sont répandues. On peut notamment citer des mesures basées sur la détection d'un signal rétrodiffusé ou réfléchi par un échantillon illuminé par un faisceau lumineux. Il peut par exemple s'agir de la spectrométrie Raman, de l'imagerie de fluorescence ou de la spectrométrie de réflectance diffuse. La spectrométrie de réflectance diffuse, usuellement désignée par l'acronyme anglosaxon DRS (Diffuse Reflectance Spectroscopy), consiste à exploiter la lumière rétrodiffusée par un objet diffusant soumis à une illumination ponctuelle. Cette technique a été décrite dans le document EP2762064 ainsi que dans la publication Sorgato V. "Non-contact quantitative diffuse reflectance spectroscopy", Proceedings of SPIE Vol. 9538, 95380U (2015), désignée par la suite par « Sorgato 2015 ». Ces documents décrivent une sonde comportant une fibre optique centrale, reliée à une source de lumière et destinée à diriger un faisceau lumineux fin vers un échantillon, en particulier un échantillon de peau. Des fibres optiques disposées autour de la fibre optique centrale, dites fibres de détection, recueillent un rayonnement optique rétrodiffusé par l'échantillon. Des moyens d'analyse spectrale du rayonnement rétrodiffusé permettent d'estimer des propriétés optiques de l'échantillon, en particulier des propriétés optiques de diffusion et des propriétés optiques d'absorption. La sonde peut être appliquée au contact ou à distance de l'échantillon à caractériser. Cette technique permet une caractérisation quantitative de ces propriétés optiques, mais il s'agit d'une caractérisation locale. Les propriétés optiques sont déterminées dans un volume élémentaire de l'échantillon, ce volume élémentaire dépendant de la géométrie de la sonde. Il est possible d'effectuer des mesures de type DRS en différentes positions discrètes sur un échantillon, selon un maillage spatial, de façon à obtenir une caractérisation bidimensionnelle des propriétés optiques de l'échantillon. Une telle démarche est présentée dans la publication Nichols BS et al. "A Quantitative Diffuse Reflectance Imaging (QDRI) System for Comprehensive 3050528 2 Surveillance of the Morphological Landscape in Breast Tumor Margins", PLoS ONE 10, 2015, désignée par la suite par "Nichols 2015". Le dispositif et la méthode décrits dans cette publication permettent une caractérisation des propriétés optiques selon un maillage bidimensionnel discret comprenant 49 points répartis régulièrement selon un arrangement 5 matriciel. De façon alternative, à la DRS, l'imagerie spectrale de réflectance, désignée par l'acronyme anglosaxon MSI « Multi spectral Imaging » permet une caractérisation bidimensionnelle d'une propriété optique d'un échantillon. La publication Bjorgan A, "Estimation of skin optical parameters for real-time hyperspectral imaging applications", J. Biomed. Opt. 19(6), 066003, Jun 10 04, 2014, désignée par la suite par « Bjorgan 2014 », en est un exemple. Cette publication décrit l'utilisation d'images d'un échantillon, formées dans différentes bandes spectrales, pour évaluer une distribution spatiale bidimensionnelle du coefficient d'absorption de l'échantillon, et cela dans les différentes bandes spectrales. L'évaluation se base sur une valeur du coefficient d'absorption réduit préalablement déterminé en se basant sur des données disponibles dans la 15 littérature. Étant basée sur des images, cette méthode peut être appliquée selon un large champ d'observation, et présente à la fois des bonnes résolutions spatiale et spectrale. Cependant, comme décrit par la suite, les inventeurs ont constaté qu'une telle méthode peut aboutir à une erreur de mesure importante lorsqu'on cherche à obtenir une quantification précise du coefficient d'absorption. 20 Les inventeurs ont recherché une méthode permettant une caractérisation quantitative des propriétés optiques de diffusion et d'absorption d'un échantillon, conférant un champ d'observation étendu, de façon à obtenir rapidement une distribution spatiale bidimensionnelle de ces propriétés optiques. Cette méthode met en oeuvre des dispositifs simples et peu coûteux, et est suffisamment robuste pour pouvoir être appliquée au-delà d'un laboratoire de recherche, 25 par exemple dans une salle d'opération ou dans un laboratoire d'analyses. EXPOSE DE L'INVENTION Un objet de l'invention est un procédé d'estimation de propriétés optiques d'un échantillon, comportant les étapes suivantes : a) illumination de la surface de l'échantillon à l'aide d'un premier faisceau lumineux 30 produit par une première source de lumière ; b) acquisition d'une première image d'un rayonnement lumineux rétrodiffusé par l'échantillon ainsi illuminé, à l'aide d'un premier capteur d'image ; 3050528 3 c) à l'aide de la première image, détermination, en plusieurs régions d'intérêt de la surface de l'échantillon, d'une première grandeur d'intérêt, représentative d'une quantité dudit rayonnement lumineux rétrodiffusé par l'échantillon ; d) estimation d'une première propriété optique, en chacune desdites régions d'intérêt, à 5 l'aide de ladite première grandeur d'intérêt ; le procédé étant caractérisé en ce que préalablement à l'étape d), le procédé comporte une étape de mesure d'une propriété optique, dite auxiliaire, de l'échantillon, différente de la première propriété optique, chaque estimation réalisée au cours de l'étape d) prenant en compte la propriété optique auxiliaire ainsi mesurée sur l'échantillon. 10 Lors de l'étape c), une grandeur d'intérêt est obtenue pour chaque région d'intérêt. Ainsi, la propriété optique auxiliaire, mesurée sur l'échantillon, est utilisée pour déterminer la première propriété optique dudit échantillon. Le premier faisceau lumineux est émis, par la première source de lumière, dans une première bande spectrale. De préférence, la première propriété optique et la propriété optique auxiliaire 15 sont déterminées dans une même bande spectrale. Selon un mode de réalisation, la première propriété optique est une propriété d'absorption de la lumière de l'échantillon, la propriété optique auxiliaire étant une propriété de diffusion de la lumière dans l'échantillon. Selon un autre mode de réalisation, la première propriété optique est une propriété de diffusion de la lumière de l'échantillon, la propriété optique auxiliaire étant 20 une propriété d'absorption de la lumière dans l'échantillon. De préférence, l'étape de mesure de la propriété optique auxiliaire comporte les sous-étapes suivantes : i) illumination de la surface de l'échantillon à l'aide d'un faisceau lumineux auxiliaire, produit par une source de lumière auxiliaire, de façon à former, à la surface de 25 l'échantillon, une zone élémentaire d'illumination, correspondant à la partie de la surface illuminée par ledit faisceau lumineux auxiliaire ; ii) mesure, à l'aide d'un photodétecteur auxiliaire, d'un signal de rétrodiffusion, représentatif d'un rayonnement rétrodiffusé par l'échantillon, ledit rayonnement rétrodiffusé émanant d'une zone élémentaire de rétrodiffusion, séparée de la zone 30 élémentaire d'illumination et située à une distance, dite distance de rétrodiffusion, de cette dernière ; iii) estimation de la propriété optique auxiliaire de l'échantillon en fonction du signal de rétrodiffusion mesuré lors de la sous-étape ii). 3050528 4 Le faisceau lumineux auxiliaire est émis, par la source de lumière auxiliaire, dans une deuxième bande spectrale. La première bande spectrale et la deuxième bande spectrale se recouvrent ou sont confondues. La sous-étape ii) peut comporter l'obtention, à partir du signal de rétrodiffusion, d'une grandeur 5 d'intérêt dite auxiliaire, représentative d'une quantité du rayonnement rétrodiffusé par l'échantillon, à une distance de rétrodiffusion, suite à son illumination par le faisceau lumineux auxiliaire; la sous-étape iii) comporte alors une comparaison de ladite grandeur d'intérêt auxiliaire avec une pluralité d'estimations de cette grandeur d'intérêt auxiliaire, chaque estimation étant réalisée en considérant une valeur prédéterminée de ladite propriété 10 optique auxiliaire. La propriété optique auxiliaire est alors obtenue en déterminant la propriété optique auxiliaire minimisant ladite comparaison. De préférence, la sous-étape ii) est effectuée en considérant différentes distances de rétrodiffusion, l'étape iii) prenant en compte les signaux de rétrodiffusion respectivement mesurés à chacune de ces distances de rétrodiffusion. 15 Selon un mode de réalisation, lors de la sous-étape i), l'échantillon est illuminé successivement ou simultanément par une pluralité de faisceaux lumineux auxiliaires de façon à former, à la surface de ce dernier, une pluralité de zones élémentaires d'illumination espacées les unes des autres ; les sous-étapes ii) et iii) sont alors mises en oeuvre en considérant un signal rétrodiffusé émis à au moins une distance de rétrodiffusion de chaque zone élémentaire d'illumination, de 20 façon à estimer la propriété optique auxiliaire dans autant d'éléments de surface que de faisceaux lumineux auxiliaires Selon un mode de réalisation, l'étape d) comporte une comparaison de la première grandeur d'intérêt déterminée lors de l'étape c) avec une pluralité d'estimations de cette grandeur d'intérêt, chaque estimation étant réalisée en considérant une valeur prédéterminée de la 25 première propriété optique, et en prenant en compte la propriété optique auxiliaire préalablement mesurée sur l'échantillon. La première propriété optique peut notamment être obtenue en déterminant la propriété optique minimisant ladite comparaison. Selon un mode de réalisation, l'étape c) comporte une segmentation de la première image en différentes régions d'intérêts, une première grandeur d'intérêt étant déterminée en chaque 30 région d'intérêt. La segmentation peut être réalisée en fonction de l'intensité de pixels de l'image, ou en se basant sur une estimation a priori de la composition de l'échantillon au niveau de chaque région d'intérêt. Chaque région d'intérêt de l'image est optiquement couplée à une région d'intérêt de la surface de l'échantillon. Une région d'intérêt de l'image peut comprendre 3050528 5 un pixel ou une pluralité de pixels de ladite image. L'étape c) peut comprendre les sous-étapes suivantes : i) détermination d'une intensité représentative de chaque région d'intérêt de la première image ; 5 application d'un premier facteur de calibration à ladite intensité représentative de chaque région d'intérêt, de façon à obtenir ladite première grandeur d'intérêt. Le premier facteur de calibration peut être obtenu à partir d'un échantillon de calibration dont les propriétés optiques sont connues. Il peut représenter un ratio entre une quantité de rayonnement rétrodiffusé par ledit échantillon de calibration, cette quantité étant modélisée, 10 sur une intensité mesurée en plaçant ledit échantillon de calibration à la place de l'échantillon examiné. Selon un mode de réalisation, les étapes c) à d) sont mises en oeuvre de façon itérative en considérant, lors de la première itération, un premier facteur de calibration correspondant à un échantillon de calibration choisi arbitrairement. Au cours d'une itération suivante, lors de l'étape 15 c), la première grandeur d'intérêt est déterminée en appliquant un premier facteur de calibration correspondant à la première propriété optique estimée lors de l'étape d) d'une itération précédente. Le premier facteur de calibration peut également être choisi de façon à correspondre à la propriété auxiliaire mesurée sur l'échantillon. 20 Selon un mode de réalisation, la propriété optique auxiliaire de l'échantillon peut être mesurée en une pluralité d'éléments de surface (5i) de l'échantillon. Le procédé peut alors comporter une interpolation de ladite propriété optique auxiliaire à partir d'au moins deux propriétés optiques auxiliaires préalablement mesurées, de façon à établir une distribution spatiale, notamment bidimensionnelle, de ladite propriété optique auxiliaire. 25 Selon un mode de réalisation, la propriété optique auxiliaire de l'échantillon est mesurée selon au moins un élément de surface, ledit élément de surface étant déterminé à partir de la première image acquise lors de l'étape b). Un autre objet de l'invention est un dispositif pour estimer des propriétés optiques d'un échantillon comportant : 30 une première source de lumière, apte à émettre un premier faisceau lumineux se propageant vers l'échantillon ; un premier capteur d'image, apte à acquérir une première image d'un rayonnement lumineux rétrodiffusé par l'échantillon, dit premier rayonnement rétrodiffusé ; 3050528 6 une source de lumière auxiliaire, apte à émettre un faisceau lumineux auxiliaire se propageant vers l'échantillon, de façon à former, à la surface de ce dernier, une zone élémentaire d'illumination, correspondant à la partie d'une surface de l'échantillon illuminé par ledit faisceau lumineux auxiliaire ; 5 un photodétecteur auxiliaire apte à mesurer un rayonnement rétrodiffusé par l'échantillon, dit deuxième rayonnement rétrodiffusé, à une distance de rétrodiffusion de la zone élémentaire d'illumination, suite à l'illumination de l'échantillon par la source de lumière auxiliaire; un microprocesseur, configuré pour déterminer une propriété optique auxiliaire de 10 l'échantillon à partir du rayonnement rétrodiffusé mesuré par le photodétecteur auxiliaire, le microprocesseur étant également configuré pour déterminer une première propriété optique différente de la propriété optique auxiliaire, en différentes régions d'intérêt de la surface de l'échantillon, à partir de la première image acquise par le premier capteur d'image, en prenant en compte ladite propriété optique auxiliaire. 15 Le photodétecteur auxiliaire peut notamment être un capteur d'image. Il peut s'agir du premier capteur d'image. Selon un mode de réalisation, la source de lumière auxiliaire est apte à former simultanément une pluralité de zones élémentaires d'illumination, espacées les unes des autres, à la surface de l'échantillon. Le dispositif peut alors comporter un masque, interposé entre la source de lumière 20 auxiliaire et l'échantillon, le masque comportant des portions transparentes, permettant le passage du faisceau lumineux auxiliaire, et des portions opaques, bloquant le faisceau lumineux auxiliaire, de telle sorte que chaque zone élémentaire d'illumination est formée par une projection d'une portion transparente sur la surface de l'échantillon. Selon un mode de réalisation, la source de lumière auxiliaire est apte à former successivement 25 une pluralité de zones élémentaires d'illumination, espacées les unes des autres, à la surface de l'échantillon. Selon un mode de réalisation, le microprocesseur est apte à mettre en oeuvre les étapes c) et d) du procédé tel que décrit ci-dessus. Un autre objet de l'invention est un procédé d'estimation de propriétés optiques d'un 30 échantillon, comportant les étapes suivantes : a) illumination de la surface de l'échantillon à l'aide d'un premier faisceau lumineux, produit par une première source de lumière ; 3050528 7 b) acquisition d'une première image d'un rayonnement lumineux rétrodiffusé par l'échantillon, ainsi illuminé, à l'aide d'un premier capteur d'image; c) sélection d'une région d'intérêt dans la première image ; d) détermination d'une pluralité de premiers facteurs de calibration, chaque premier 5 facteur de calibration prenant la forme d'un ratio entre : une modélisation d'une quantité de rayonnement rétrodiffusé par l'échantillon de calibration dont la première propriété optique est connue, et une intensité représentative d'une image de calibration acquise, par ledit premier capteur d'image, lorsqu'un échantillon de calibration est illuminé par la 10 première source de lumière, de telle sorte qu'on dispose d'une pluralité de facteurs de calibration correspondant respectivement à des échantillons de calibration ayant des premières propriétés optiques de différentes valeurs, chaque premier facteur de calibration étant associé à une valeur de ladite première propriété optique ; 15 e) détermination d'une première grandeur d'intérêt, représentative d'une quantité de rayonnement rétrodiffusé par l'échantillon, dans ladite région d'intérêt, en multipliant une intensité représentative de la région d'intérêt par un premier facteur de calibration ; f) estimation d'une première propriété optique en différentes régions d'intérêt de la 20 surface de l'échantillon, à l'aide de ladite première grandeur d'intérêt ainsi déterminée; g) réitération des étapes e) et f) en considérant, lors de l'étape e), un premier facteur de calibration correspondant à valeur de la première propriété optique estimée lors de l'étape f) de l'itération précédente, la réitération étant effectuée jusqu'à l'atteinte d'un critère de convergence. 25 Les étapes c) à g) peuvent être effectuées pour différentes régions d'intérêt de la première image, chaque région d'intérêt de l'image étant optiquement couplée à une région d'intérêt à la surface de l'échantillon. Selon un mode de réalisation, le procédé comprend une mesure d'une propriété optique auxiliaire, telle que précédemment définie, sur un ou plusieurs éléments de surface de 30 l'échantillon. Lors de l'étape e) la première grandeur d'intérêt peut être déterminée en considérant un premier facteur de calibration correspondant à ladite propriété optique auxiliaire. 3050528 8 Par premier facteur de calibration correspondant à une propriété optique, on entend un premier facteur de calibration représentatif d'un échantillon de calibration présentant ladite propriété optique. Un tel premier facteur de calibration peut être mesuré ou interpolé en se basant sur des facteurs de calibration mesurés. 5 D'autres avantages et caractéristiques ressortiront plus clairement de la description qui va suivre de modes particuliers de réalisation de l'invention, donnés à titre d'exemples non limitatifs, et représentés sur les figures listées ci-dessous. FIGURES La figure 1A représente une première composante d'un dispositif selon l'invention, pour estimer 10 une propriété optique d'un échantillon à partir d'une première image, selon une première modalité. La figure 1B représente une deuxième composante d'un dispositif selon l'invention, pour estimer une propriété optique auxiliaire d'un échantillon, selon une deuxième modalité. La figure 2 illustre les étapes d'un procédé selon l'invention. 15 La figure 3A montre une image d'échantillons tests obtenue à l'aide de la configuration représentée sur la figure 1A. La figure 3B représente une comparaison entre une première grandeur d'intérêt, mesurée à l'aide de la configuration représentée sur la figure 1A, et des modélisations de ladite première grandeur d'intérêt, chaque modélisation étant effectuée en fonction d'une valeur d'un 20 coefficient d'absorption et d'un coefficient de diffusion réduit de l'échantillon. La figure 3C montre des coefficients d'absorption d'un échantillon déterminés en mettant en oeuvre les étapes décrites en lien avec la figure 2 dans trois longueurs d'onde différentes. La figure 4A montre des estimations, à différentes longueurs d'onde, du coefficient d'absorption de quatre échantillons test en se basant sur une prise en compte arbitraire d'un coefficient de 25 diffusion réduit. La figure 4B illustre les erreurs de mesure associées aux estimations présentées sur la figure 4A. La figure 4C montre des estimations, à différentes longueurs d'onde, du coefficient d'absorption de quatre échantillons test en se basant sur une mesure d'un coefficient de diffusion réduit effectuée sur chaque échantillon. La figure 4D illustre les erreurs de mesure associées aux 30 estimations présentées sur la figure 4C. La figure 5A représente un mode de réalisation de la deuxième composante d'un dispositif selon l'invention. 3050528 9 Les figures 5B et 5C schématisent respectivement deux modes de réalisation différents d'un dispositif selon l'invention. La figure 5D montre un masque pouvant être appliqué dans le mode de réalisation exposé sur la figure 5C. La figure 6 montre trois images pouvant être obtenues en mettant en oeuvre l'invention. 5 EXPOSE DE MODES DE RÉALISATION PARTICULIERS La figure lA représente une première composante la d'un dispositif 1 permettant de former une image d'un échantillon et d'en déduire une propriété optique, ou une distribution spatiale d'une telle propriété optique. Le terme propriété optique désigne une propriété gouvernant l'absorption ou la diffusion de la 10 lumière dans l'échantillon. Il peut s'agir d'une propriété optique d'absorption, par exemple le coefficient d'absorption lia, ou une propriété optique de diffusion par exemple un coefficient de diffusion réduit us' ou un coefficient de diffusion us ou un coefficient d'anisotropie de diffusion. Ces coefficients sont connus de l'homme du métier. La première composante la permet d'appliquer une première modalité, dite d'imagerie 15 spectrale de diffusion, à l'échantillon. D'une façon générale, selon cette première modalité, l'échantillon 5 est illuminé par une première source de lumière 10 et on forme une première image /1 (À) de ce dernier à l'aide d'un premier capteur d'image 18. La source de lumière peut être une source de lumière blanche, auquel cas le premier capteur d'image 18 est spectralement résolu, permettant l'acquisition d'images de l'échantillon selon différentes bandes spectrales À. 20 Par bande spectrale, on entend une longueur d'onde ou une plage de longueurs d'onde. De façon alternative, la première source de lumière 10 peut émettre un premier faisceau lumineux 12, incident à l'échantillon, successivement selon différentes bandes spectrales À, auquel cas le premier capteur d'image 18 acquiert autant d'images de l'échantillon que de bandes spectrales. 25 Selon l'exemple représenté, la première source de lumière 10 est une source de lumière blanche, associée à une roue à filtres 15, cette dernière comportant une pluralité de filtres dont la bande passante définit la bande spectrale X du premier faisceau lumineux 12, émis par la source et se propageant jusqu'à l'échantillon 5. La première source de lumière 10 peut être couplée à un diffuseur 16, permettant d'augmenter un angle solide Q d'illumination de 30 l'échantillon 5. De préférence, une part significative de la surface 5s de l'échantillon est illuminée par le premier faisceau lumineux 12, l'aire de la surface illuminée étant de préférence supérieure à 4 cm2, voire supérieure à 10 cm2. Cette surface est typiquement de l'ordre de 100 cm2. 3050528 10 Sous l'effet de l'illumination par le premier faisceau lumineux 12, l'échantillon rétrodiffuse un rayonnement 14, dit premier rayonnement rétrodiffusé, ce dernier étant détecté par le premier capteur d'image 18, à travers un objectif 17. Le dispositif comporte également un processeur 30, par exemple un microprocesseur, relié à 5 une mémoire 32 et apte à exécuter des instructions stockées dans cette dernière, de façon à mettre en oeuvre certaines étapes du procédé décrit ci-après. Le processeur peut être relié à un écran 34. Le processeur 30 effectue lesdites étapes à partir des premières images acquises par le capteur d'image 18 et/ou à partir de signaux de rétrodiffusion détectés par un photodétecteur auxiliaire 28 décrit en lien avec la figure 1B. 10 À partir de chaque première image 11(À) acquise par le capteur d'image 18, dans une bande spectrale X, on peut estimer des propriétés optiques de l'échantillon 5, dans ladite bande spectrale, comme décrit ci-après, en lien avec la figure 2. La figure 1B représente une deuxième composante lb du dispositif 1 permettant la mise en oeuvre de l'invention selon une deuxième modalité dite de spectrométrie de réflectance diffuse 15 (DRS). Comme décrit dans les documents relatifs à la DRS cités dans l'art antérieur, une source de lumière auxiliaire 20, par exemple une source de lumière blanche, est couplée à une fibre optique 21, dite fibre optique d'illumination, de façon à former un faisceau lumineux auxiliaire 22 se propageant vers l'échantillon 5. La fibre optique d'illumination s'étend entre une extrémité proximale 21p disposée face à la source de lumière auxiliaire 20, collectant la lumière émise par 20 cette dernière, et une extrémité distale 21a disposée face à l'échantillon 5. L'extrémité distale 21a peut être appliquée contre l'échantillon 5, ou être maintenue à distance de l'échantillon 5, comme représenté sur la figure 1B. Dans ce dernier cas, un système optique 26, par exemple un objectif ou une lentille, permet une focalisation du faisceau lumineux auxiliaire 22 sur la surface 5s de l'échantillon. Le système optique 26 peut présenter un facteur de grandissement inférieur 25 ou supérieur à 1. Dans l'exemple représenté, le facteur de grandissement est égal à 1. Le faisceau lumineux auxiliaire 22 forme, à la surface 5s de l'échantillon, une zone élémentaire d'illumination 23, correspondant à l'intersection du faisceau lumineux auxiliaire 22 avec ladite surface. Le terme zone élémentaire désigne une zone de forme délimitée et comportant une dimension inférieure ou égale à 1 cm, et de préférence inférieure ou égale à 1 mm. Par 30 dimension, on entend un diamètre, une diagonale ou une largeur. L'aire d'une telle zone élémentaire est de préférence inférieure à 1 cm2, voire inférieure à 1 mm2. Sur l'exemple représenté sur la figure 1B, la zone élémentaire d'illumination 23 est située dans un plan focal 3050528 11 image du système optique 26, tandis que l'extrémité distale 21d de la fibre optique d'illumination 21 est située dans le plan focal objet de ce système optique. Le diamètre de la fibre optique d'illumination 21 est généralement compris entre 50 um et 1 mm, et est par exemple égal à 500 um. 5 Le dispositif comporte une pluralité de fibres optiques de détection 27, chacune s'étendant entre une extrémité distale 27d, faisant face à l'échantillon, et une extrémité proximale 27g. Chaque fibre de détection 27 est apte à collecter un rayonnement rétrodiffusé 24, dit deuxième rayonnement rétrodiffusé, par l'échantillon 5 lorsque ce dernier est illuminé par le faisceau lumineux auxiliaire 22. Les extrémités distales 27d de chaque fibre de détection 27 sont 10 disposées à des distances différentes de l'extrémité distale 21d de la fibre d'illumination 21, de façon à ce que chaque fibre de détection soit apte à collecter un rayonnement rétrodiffusé par l'échantillon, émanant de la surface de l'échantillon à une distance Dn, de la zone d'illumination 23, dite distance de rétrodiffusion, différente l'une de l'autre. Le principe de cette modalité est de déterminer une quantité de photons rétrodiffusés par 15 l'échantillon selon une ou, de préférence plusieurs distances de rétrodiffusion Dn. Pour cela, chaque fibre optique de détection 27 collecte un rayonnement rétrodiffusé 24 émanant d'une zone élémentaire de détection 25, cette dernière s'étendant à une distance de rétrodiffusion Dn, de la zone élémentaire d'illumination 23. Le photodétecteur auxiliaire 28 détecte un signal rétrodiffusé Sn dont l'intensité dépend d'une quantité de photons rétrodiffusés à la distance de 20 rétrodiffusion Dn. Sur la figure 1B, on a représenté trois fibres de détection 27 : une première fibre de détection 271, apte à collecter un rayonnement de diffusion 241 émanant de la surface de l'échantillon au niveau d'une zone élémentaire de détection 251 s'étendant à une première distance de rétrodiffusion D1 par rapport à la zone 25 élémentaire d'illumination 23. Une deuxième fibre de détection 272, apte à collecter un rayonnement de diffusion 242 émanant de la surface de l'échantillon au niveau d'une zone élémentaire de détection 252 s'étendant à une deuxième distance de rétrodiffusion D2 par rapport à la zone élémentaire d'illumination 23. 30 Une troisième fibre de détection 273, apte à collecter un rayonnement de diffusion 243 émanant de la surface de l'échantillon au niveau d'une zone élémentaire de détection 3050528 12 253 s'étendant à une troisième distance de rétrodiffusion D3 par rapport à la zone élémentaire d'illumination 23. Les première, deuxième et troisième distances de rétrodiffusion sont non nulles et différentes les unes des autres. Selon le facteur de grandissement du système optique 26, elles peuvent 5 s'étendre de quelques dizaines de um à quelques mm, voire cm de la zone élémentaire d'illumination 23. Les zones élémentaires de rétrodiffusion 25 sont de préférence séparées de la zone élémentaire d'illumination 23, le terme séparé signifiant qu'il n'y a pas d'intersection entre une zone élémentaire de rétrodiffusion 25 et la zone élémentaire d'illumination 23. Plusieurs fibres optiques de détection peuvent être conjuguées à des zones élémentaires de 10 détection s'étendant à une même distance de rétrodiffusion de la zone élémentaire d'illumination. Dans ce cas, ces fibres peuvent être disposées selon un cercle autour de la zone élémentaire d'illumination. Cela permet d'augmenter la quantité de rayonnement de rétrodiffusion détecté, correspondant à une même distance de rétrodiffusion. On améliore ainsi la sensibilité de la mesure. Cela permet également d'obtenir d'effectuer une moyenne des 15 rayonnements rétrodiffusés à chaque distance de rétrodiffusion. De telles configurations, dites concentriques, sont décrites dans le document EP2762064 ou dans la publication « Sorgato 2015 » précédemment citée. L'extrémité proximale 27p de chaque fibre optique de détection est couplée à un commutateur optique 29, permettant un couplage sélectif de chaque fibre optique avec un photodétecteur 20 auxiliaire 28, ce dernier étant, dans cet exemple, un spectrophotomètre. Le spectrophotomètre est apte à former un spectre en longueur d'onde de chaque rayonnement rétrodiffusé selon une distance de rétrodiffusion. Le processeur 30 est apte à traiter les spectres de façon à estimer une propriété optique auxiliaire de l'échantillon, comme expliqué en lien avec la figure 2. Une différence notable entre la première composante la du dispositif, représentée sur la figure 25 lA et la deuxième composante lb du dispositif, représentée sur la figure 1B, concerne les extensions spatiales du premier faisceau lumineux 12 et du faisceau lumineux auxiliaire 22. Comme précédemment évoqué, la surface de l'échantillon éclairée par le premier faisceau lumineux 12 est supérieure à quelques cm2, et le plus souvent supérieur 5 cm2, voire 10 ou 100 cm2, tandis que la surface de l'échantillon éclairée par le faisceau lumineux auxiliaire 22 est 30 inférieure à 1 cm2. La première modalité permet d'effectuer une détermination d'une propriété optique simultanément sur une surface importante de l'échantillon, c'est-à-dire sur différentes régions d'intérêt de la surface 5s de l'échantillon, tandis que la deuxième modalité permet une 3050528 13 mesure quantitative d'une propriété optique locale de l'échantillon. Ainsi, la première modalité permet d'obtenir, à partir d'une même image /1(À), une distribution spatiale bidimensionnelle d'une propriété optique, ce que ne permet pas la deuxième modalité, sauf à multiplier le nombre de fibres optiques d'illumination et de fibres optiques de détection, comme décrit dans 5 la publication « Nichols 2015» précédemment citée. On va à présent décrire les principales étapes d'un procédé de détermination de propriétés optiques d'un échantillon, en lien avec la figure 2. Les étapes 100 à 150 concernent la première modalité, en l'occurrence l'imagerie spectrale de diffusion, tandis que les étapes 200 à 240 sont concernent la deuxième modalité, en l'occurrence la spectrométrie de réflectance diffuse. 10 Étape 100: illumination de l'échantillon avec la première source de lumière. Au cours de cette étape, l'échantillon est illuminé avec la source de lumière blanche 10, émettant un premier faisceau lumineux 12 se propageant jusqu'à l'échantillon dans une bande spectrale X ajustée par la roue à filtres 15. Dans cet exemple, la source de lumière utilisée est une source halogène Schott de référence KL 2500 LCD. Elle est associée à une roue à filtre 15 15 comportant trois filtres définissant des bandes spectrales X respectivement centrées sur les longueurs d'onde 500 nm, 600 nm et 700 nm. La largeur de bande de chaque filtre, correspondant à la largeur à mi-hauteur du pic d'émission, est d'environ 10 nm. Le diffuseur 16 est un anneau au niveau duquel aboutissent des guides de lumière provenant de la première source de lumière 10. Il permet de répartir la lumière, de façon la plus uniforme possible, selon 20 un angle solide d'émission SI Les étapes 110 et 150 ci-dessous sont décrites en référence à une bande spectrale, sachant qu'elles peuvent être mises en oeuvre successivement selon chaque bande spectrale du premier faisceau 12 émis par la première source de lumière et atteignant l'échantillon. Étape 110 : Acquisition d'une image. 25 Le premier capteur d'image acquiert une première image /1(À) d'un rayonnement 14 formé par une rétrodiffusion, par l'échantillon 5 du premier faisceau lumineux 12, dans une longueur d'onde ou une bande spectrale À. La figure 3A représente une telle image, l'échantillon étant formé de différents échantillons test, présentant des propriétés de diffusion optique différentes, et décrits en lien avec les figures 4A à 4D. 30 L'image acquise 11(À) est de préférence corrigée d'une image de noir 'black' obtenue dans le noir, selon une même durée d'acquisition, la première source de lumière étant éteinte. L'image de noir' 'black est représentative du bruit du premier capteur d'image et d'un niveau d'obscurité 3050528 14 du dispositif. Elle est également de préférence normalisée par une image dite de blanc /I.latfield, également désignée par le terme anglosaxon « flat field », destinée à corriger les hétérogénéités spatiales de l'illumination de l'échantillon et de l'efficacité de détection de l'image. En effet, l'éclairement de l'échantillon par la première source de lumière peut ne pas être homogène. De 5 plus, en particulier du fait de composants optiques dans l'objectif 17, l'efficacité de détection n'est pas homogène, notamment en raison d'effets de vignettage à la périphérie de l'image. L'image acquise est normalisée par l'image de blanc Iflatfield, cette dernière étant obtenue en effectuant une image d'un échantillon dont la rétrodiffusion est considérée comme homogène, et en normalisant l'image ainsi obtenue par la valeur d'intensité du pixel le plus intense de ladite 10 image. Le recours à une telle image de blanc est une pratique courante. Étape 120: Segmentation de la première image en régions d'intérêt et sélection d'une région d'intérêt. Au cours de cette étape, on sélectionne une ou plusieurs régions d'intérêt ROI dans la première image /1(À) obtenue au cours de l'étape 110. Une région d'intérêt peut correspondre à 15 différents pixels adjacents de l'image présentant un même niveau d'intensité, ou un niveau d'intensité comparable. Chaque région d'intérêt de l'image est associée à une région d'intérêt de l'échantillon, à laquelle elle est couplée par l'objectif 17. Dans l'image de la figure 3A, on a représenté 4 régions d'intérêt, chacune étant délimitée par un rectangle. Chaque région d'intérêt correspond à un échantillon test décrit en lien avec les figures 4A à 4D. Une telle 20 sélection permet d'éviter des artéfacts dus notamment à la réflexion spéculaire. Une région d'intérêt peut ne comprendre qu'un seul pixel. Elle peut également regrouper des pixels non adjacents mais présentant une même intensité ou une intensité comparable. Une région d'intérêt peut également correspondre à une région de l'échantillon dans laquelle on considère, de façon arbitraire ou sur la base de mesures préalables, qu'une propriété optique 25 de l'échantillon est homogène, ou qu'une composition de l'échantillon est homogène. Étape 130: Détermination d'une intensité moyenne dans chaque région d'intérêt. Au cours de cette étape, on détermine une grandeur d'intérêt représentative de chaque région d'intérêt. Si /1,Ro/(.1) désigne l'image de la région d'intérêt ROI, on effectue une moyenne de l'intensité des pixels dans cette région d'intérêt, de façon à obtenir une intensité moyenne 30 / oI (À). 3050528 15 Cette intensité moyenne peut être normalisée par une durée d'acquisition ôt de l'image, de telle Rst01 (A- sorte que /M0'7(.1) = 11, ), (1) l'exposant norm désignant le fait que l'image est normalisée par la durée d'acquisition. Étape 140 : Prise en compte de facteurs correctifs. 5 Cette étape consiste à appliquer un facteur de calibration f (À), dit premier facteur de calibration, à l'intensité moyenne. Ce premier facteur de calibration prend en compte l'intensité du premier faisceau lumineux 12 incident à l'échantillon, ainsi que la sensibilité du premier capteur d'image 18. Il est déterminé au cours d'une phase de calibration, mettant en oeuvre un échantillon de calibration, ou fantôme de calibration, dont les propriétés optiques sont connues. 10 Le premier facteur de calibration f (À) correspond à un ratio entre une réflectance modélisée p odel I\ -norm de l'échantillon de calibration et l'intensité moyenne d'une image Icalib (À) de "crnalib l'échantillon acquise par le dispositif, en conservant les mêmes paramètres expérimentaux que lors de l'analyse de l'échantillon : mêmes réglages de la source de lumière 10 et du capteur d'image 18, même distance de l'échantillon de calibration vis-à-vis de la source et du capteur 15 d'image, même niveau obscurité dans lequel est placé l'échantillon. RMicle!(A.)~ormLe premier facteur de calibration est tel que f (À) = m (2) Icalib Par ailleurs, certains paramètres du dispositif peuvent évoluer entre l'instant to ou la calibration précédemment décrite a été réalisée et l'instant t de la mesure sur l'échantillon examiné. Il peut par exemple s'agir de la variation de l'intensité de la source de lumière, de la réponse des guides 20 de lumière, ou une variation de la sensibilité du capteur d'image, et, d'une façon plus générale, de différents composants du dispositif. Afin de prendre en compte cette variabilité, un premier facteur de correction k(À) est appliqué, de telle sorte que k(À) = It-re f(A-) , (3) où It_ref (À) et lto-re f (À.) It0-re f (Al) fantôme de référence respectivement à l'instant de la mesure t et à un instant initial t0. Le 25 fantôme de référence peut être un échantillon de calibration ou un fantôme comportant une surface dont la réflectance est constante. Il peut par exemple s'agir d'une pastille appliquée sur un support que l'on dispose face au capteur d'image 18 avant d'effectuer les mesures. Les mesures du fantôme de référence, à l'instant initial to et à l'instant de la mesure t sont effectuées selon les mêmes paramètres expérimentaux. 30 L'application du premier facteur de calibration f (À) et du premier facteur de correction k(À) permet d'obtenir une réflectance R(.1) de l'échantillon analysé 5, de telle sorte que : désignent respectivement une intensité, mesurée par le capteur d'image 18, d'un 3050528 16 R(À) = TnroT(À) x k (À) x f (À) (4) Ainsi, on obtient, à partir de la première image, une première grandeur d'intérêt R(.1) représentative de la réflectance de l'échantillon, dans la région d'intérêt considérée, cette réflectance représentant une quantité de rayonnement rétrodiffusé par ladite région d'intérêt 5 normalisée par une intensité du premier faisceau lumineux 12. Étape 150: estimation d'une première propriété optique Au cours de l'étape 150, on cherche à estimer une première propriété optique pi de l'échantillon, par exemple le coefficient d'absorption ma(À), dans la bande spectrale À du premier faisceau lumineux 12. 10 Pour cela, on utilise une modélisation, dans la bande spectrale considérée, de la réflectance de l'échantillon, cette modélisation prenant en compte une pluralité de propriétés optiques p, parmi lesquelles le coefficient d'absorption lia, le coefficient de diffusion réduit Vs. On effectue des modélisations en faisant varier le coefficient d'absorption et le coefficient de diffusion réduit respectivement entre une valeur minimale et une valeur maximale. On obtient ainsi une 15 réflectance modélisées R inwderi \(/1) pour chaque valeur de lia et de Vs considérées. La 0-441 s) modélisation peut être réalisée sur la base d'un modèle analytique ou en utilisant une méthode stochastique de type Monte Carlo. On peut effectuer une comparaison, par exemple sous la forme d'un écart quadratique entre la ri, owdpe st ) , réflectance mesurée R(À) et chaque réflectance modélisée R v . ) selon l'expression 20 A(P-a(À), Fils (À)) = j (R (À) - Ridlicjac,ilersi)(À))2 (5) avec 0.1,,(À), Vs(À)) = argmin Gl(R(À) - R(µadaf,s1) (À))2) (5') 1-ta(A),Plis(1) On a représenté, sur la figure 3B, une telle comparaison 3,(1.1,,(À), Vs(À)) en fonction de (lia (À) et Vs(À) en se basant sur une image mesurée sur un échantillon test dans une bande spectrale centrée sur 600 nm. Les propriétés optiques (u.a(À), Vs(À)) de l'échantillon sont celles 25 minimisant la différence 3,(1.1,,(À), Vs(À)). Cependant, cette différence ne prend pas une valeur minimum pour un point de l'espace (1.1,,(À), Vs(À)), mais selon une multitude de points, formant une ligne. Ainsi, le problème est indéterminé, puisqu'il ne comporte pas une solution, c'est-à-dire un couple (lia (À), Vs(À)) mais une pluralité, voire une infinité de solutions. Afin 3050528 17 d'estimer une de ces deux propriétés optiques, par exemple u.a,(À), la valeur de l'autre propriété optique, en l'occurrence Vs(.1), doit être préalablement établie. Selon l'art antérieur, pour corriger cette indétermination, on considère une valeur prédéterminée du coefficient de diffusion réduit Vs. Dans la publication Bjorgan 2014 , le 5 coefficient Vs est fixé selon un modèle de diffusion analytique, basé sur des échantillons considérés comme représentatifs de l'échantillon analysé. Les inventeurs ont constaté qu'une telle démarche conduit à une erreur de mesure pouvant être rédhibitoire, comme décrit ci-après en lien avec les figures 4A et 4B. Ils en ont conclu que pour lever cette indétermination, il est préférable non pas d'utiliser une valeur prédéterminée du 10 coefficient de diffusion réduit, censée être représentative de l'échantillon examiné, mais de mesurer cette valeur directement sur l'échantillon examiné. Pour cela, une deuxième modalité est utilisée, visant à estimer, au moins localement, une propriété optique auxiliaire de l'échantillon 5, en l'occurrence le coefficient de diffusion réduit Vs. 15 Autrement dit, la méthode proposée vise à estimer une première propriété optique de l'échantillon, par exemple lia, dans une bande spectrale, à l'aide d'une première modalité, cette estimation étant basée sur une mesure d'une propriété optique auxiliaire de l'échantillon, par exemple Vs, dans ladite bande spectrale, par le biais d'une deuxième modalité de mesure indépendante de la première modalité. 20 Dans cet exemple, la propriété optique auxiliaire est mesurée par une méthode de type DRS, visant à former une zone élémentaire d'illumination 23 sur la surface 5s de l'échantillon 5, et à détecter au moins un rayonnement diffusé et de préférence au moins deux rayonnements rétrodiffusés 241, 242, 243 respectivement selon deux distances de rétrodiffusion D1, D2, D3 non nulles et différentes l'une de l'autre. Cette méthode est appliquée localement, sur au moins un 25 élément de la surface 5i de la surface 5s, comme représenté sur la figure 1B.Le dispositif mis en oeuvre peut être celui décrit en lien avec la figure 1B, et appliqué selon la méthode décrite dans la publication « Sorgato 2015 » dont les principales étapes sont les suivantes : Étape 200: Illumination Au cours de cette étape, on illumine l'échantillon en activant la source de lumière auxiliaire 20, 30 de façon à diriger un faisceau lumineux auxiliaire, ou deuxième faisceau lumineux 22 vers la surface 5s de l'échantillon 5, l'intersection du faisceau lumineux auxiliaire 22 avec la surface 5s 3050528 18 de l'échantillon formant une zone élémentaire d'illumination 23. La source de lumière peut être par exemple une lampe tungstène halogène de référence HL 2000 de Ocean Optics. Étape 210: Détection de signaux rétrodiffusés selon différentes distances de rétrodiffusion. Au cours de cette étape, le photodétecteur auxiliaire 28 détecte au moins deux signaux 5 rétrodiffusés 241, 242 selon deux distances de rétrodiffusion D1, D2 différentes. Dans cet exemple, le photodétecteur auxiliaire 28 est un spectromètre de type MAYA2000 Pro de Ocean Optics. On forme alors, à chaque distance de rétrodiffusion Dn, un signal détecté, dit signal de rétrodiffusion, Sn(À) dans la bande spectrale À , l'indice n désignant la distance de rétrodiffusion. Chaque signal de rétrodiffusion Sn(À) est de préférence corrigé d'un signal de 10 noir S mesuré dans le noir, selon une même durée d'acquisition, la source de lumière auxiliaire étant éteinte. Le signal de noir S black est représentatif du bruit du photodétecteur auxiliaire 28 et d'un niveau d'obscurité dans lequel l'échantillon est placé. De façon alternative, l'échantillon peut être successivement illuminé, par la source de lumière auxiliaire 20, selon différentes bandes spectrales X, le photodétecteur auxiliaire 28 n'étant pas 15 résolu spectralement. Il peut en particulier s'agir d'un capteur d'image, comme décrit en lien avec les figures 5A à 5C. Étape 220: Normalisation et application de facteurs correctifs pour obtenir une grandeur d'intérêt auxiliaire, ou deuxième grandeur d'intérêt. Au cours de cette étape, chaque signal de rétrodiffusion détecté Sn(À), corrigé par le signal de 20 noir S -black, est normalisé par une durée d'acquisition (St et par une quantité de lumière Ssource(À) émise par la source de lumière auxiliaire 20. Une telle quantité de lumière peut être mesurée en disposant une fibre optique, dite de retour d'excitation 21e, entre la source de lumière auxiliaire 20 et le photodétecteur auxiliaire 28. On applique également un facteur de calibration, dit deuxième facteur de calibration, fn(À) à 25 chaque signal de rétrodiffusion Sn(À) ainsi corrigé et normalisé. Ce deuxième facteur de calibration, dépendant de la distance de rétrodiffusion Dn prend en compte l'effet de différents composants de la fibre d'illumination 21 et de la fibre de détection 27n sur le signal de rétrodiffusion. De façon analogue au premier facteur de calibration f (À) décrit dans l'étape 140, chaque deuxième facteur de calibration est déterminé au cours d'une phase de calibration, 30 mettant en oeuvre un échantillon de calibration dont les propriétés optiques sont connues. Ce facteur de calibration correspond à un ratio entre une réflectance modélisée &n'ent() de l'échantillon de calibration à la distance de rétrodiffusion Dn et une réflectance Rcalib,n(À) de 3050528 19 l'échantillon de calibration mesurée par le dispositif, à la même longueur d'onde et à la même distance de rétrodiffusion, en conservant les mêmes paramètres expérimentaux qu'avec l'échantillon analysé : mêmes réglages de la source 20 et du photodétecteur 28, même distance de l'échantillon de calibration vis-à-vis de la source et du photodétecteur. Dmodel 5 Le deuxième facteur de calibration est tel que fn(À) - '-calib'nl'`.1 (6) Rcalib,n(X) _ scaub,n(x)-sbiacka) (61, avec : Rc ali b ,n - SsourceMX8tcalib Scalib,n(À) désignant le signal de rétrodiffusion, à une distance de rétrodiffusion Dn, mesuré par le photodétecteur auxiliaire 28 lorsque l'échantillon mesuré est l'échantillon de calibration, atcallb désignant la durée de l'acquisition du signal de rétrodiffusion Scalib,n(2)- 10 Par ailleurs, certains paramètres du dispositif peuvent évoluer entre l'instant to auquel la calibration précédemment décrite a été réalisée et l'instant t de la mesure. Il peut par exemple prendre en compte une variation de la sensibilité du photodétecteur auxiliaire 28 ou le vieillissement des fibres optiques de détection 27. Afin de prendre en compte cette variabilité, un facteur de correction, dit deuxième facteur de 15 correction, kn (À) peut être appliqué, de telle sorte que kn(À) = It'i(À) (7) où ft,n(À) et fto,n(A-)' from. (À) représentent respectivement des deuxièmes facteurs de calibration mesurés aux instants t et to. Cette étape permet d'obtenir une réflectance de l'échantillon Rn(À), à la distance de rétrodiffusion Dn et à la longueur d'onde À, selon l'expression : 20 (À) x kn (À) (8) Su Rn(À) = Snso00rce-Sbiackx Sta) fn Cette réflectance correspond à une grandeur d'intérêt auxiliaire, représentative d'une intensité d'un rayonnement rétrodiffusé à une distance de rétrodiffusion Dn, non nulle, par rapport à la zone d'illumination élémentaire 23. Ce rayonnement rétrodiffusé peut également être désigné par le terme « deuxième rayonnement rétrodiffusé ». 25 Étape 230: Détermination des propriétés optiques locales de l'échantillon Pour au moins une longueur d'onde X, et en considérant au moins autant de distances de rétrodiffusion Dn différentes qu'il y a de propriétés optiques à estimer, on détermine les propriétés optiques p présentant le moindre écart entre la réflectance Rn (À), déterminée lors 3050528 20 de l'étape 220, à la longueur d'onde X, et une réflectance modélisée Rnmpode1V1..)71,, 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 Dn. Cette détermination peut être réalisée par la minimisation d'un écart quadratique, et par exemple selon l'expression : 5 p = argmin(Eirv,=i(R N désigne le nombre de distances de rétrodiffusion prises en compte, Reel est une réflectance modélisée, à la distance de rétrodiffusion Dn, 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. 10 Dans cet exemple, les propriétés optiques recherchées sont un(X) et µ',(À). Ainsi, le couple (1.1,,(X), Vs(À)) recherché est celui présentant le moindre écart entre la réflectance mesurée Rn (À), à la longueur d'onde X, et une réflectance modélisée Rnin,loadells,(À) pour différentes valeurs de p.a(X) et de p.',(X), à ladite distance de rétrodiffusion Dn. Cette détermination peut être réalisée selon l'expression : 15 (un (X), Vs (À)) = argminwup., sm)(Eirv,=ifeael,ls, (À) - Rn (À)) 2) (9'), Dmelsod (A1) désignant une réflectance modélisée, à la distance de rétrodiffusion Dn, en "n,pa,p considérant différentes valeurs de lia(X) et µ,'(À). Les différentes valeurs de réflectance modélisée Rn7ioad ells, 71, o.) sont obtenues, pour une pluralité de couples de valeurs lia(a), Vs(À) au cours d'une phase de paramétrage, par simulation 20 numérique mettant en oeuvre une méthode de type Monte Carlo ou par un modèle analytique. Étape 240 : transmission de la propriété optique auxiliaire ou réitération. Au cours de cette étape, la propriété optique auxiliaire, en l'occurrence le coefficient de diffusion réduit Vs(À) , déterminée lors de l'étape 230, est transmise pour être utilisée pour le calcul de la première propriété optique. Les propriétés optiques déterminées lors de l'étape 230 25 peuvent également être utilisées pour sélectionner un deuxième facteur de calibration, utilisé lors d'une nouvelle itération de l'étape 220, comme décrit ci-après. De préférence, lors de l'étape 220, on dispose d'une pluralité de deuxièmes facteurs de calibration fn(À), réalisés en utilisant des échantillons de calibration dont les propriétés optiques p, en particulier lia (X), et/ou Vs(À), sont différentes. On peut également calculer des 30 facteurs de calibration par interpolation entre deux facteurs de calibration respectivement imeei (À) - Rn (À)) 2) (9) 3050528 21 obtenus avec deux échantillons de calibration présentant des propriétés optiques différentes. On dispose alors de deuxièmes facteurs de calibration, notés fif (À), associés à chaque distance de rétrodiffusion Dn, et à différentes valeurs des propriétés optiques p. Lors d'une première itération de l'étape 220, on utilise un deuxième facteur de calibration fn7)(À) correspondant à 5 des propriétés optiques définies de façon arbitraire. Suite à l'étape 240, l'étape 220 peut être réitérée, en prenant en compte un deuxième facteur de calibration fn7)(À) correspondant aux propriétés optiques p déterminées lors de l'étape 230. Le procédé peut être réitéré plusieurs fois, mais les inventeurs estiment que deux itérations sont généralement suffisantes. Le deuxième facteur de correction Icn peut n'être déterminé que sur un seul échantillon de 10 calibration et utilisé quel que soit deuxième le facteur de calibration fn7)(À) appliqué. L'adaptation du deuxième facteur de calibration à des propriétés optiques de l'échantillon déterminées au cours d'une itération précédente, méthode dite du facteur de calibration adaptatif, a été décrite dans la publication Sorgato 2015. Cela permet de limiter les erreurs de mesures dues à la prise en compte d'un facteur de calibration non représentatif des propriétés 15 optiques réelles de l'échantillon examiné. Ainsi, l'étape 230 permet de mesurer, localement, une propriété optique auxiliaire de diffusion, en l'occurrence Vs(À), cette propriété étant utilisée, dans l'étape 150, pour estimer la première propriété optique, en l'occurrence 1.1%,(À) à l'aide de la première modalité. Il est préférable que la bande spectrale À, à laquelle la propriété optique auxiliaire est mesurée, soit proche ou incluse 20 dans la bande spectrale À à laquelle la première propriété optique est estimée. La mesure de la propriété optique auxiliaire µ'S(À) permet de lever l'indétermination évoquée en lien avec l'étape 150, et d'estimer, au cours cette étape, la valeur µa(À) minimisant la fonction d'écart A(p.a, V) mentionnée en lien avec l'équation (5). Les étapes 100 à 150 sont de préférence mises en oeuvre pour différentes zones d'intérêt ROI 25 de la première image, ce qui permet d'estimer la première propriété optique en chacune de ces zones d'intérêt. On obtient alors une estimation de la première propriété optique en différents points de la surface de l'échantillon. Une région d'intérêt peut correspondre à un pixel ou à un groupe de pixels de la première image. Un groupe de pixels peut être composé de pixels de même intensité, ou d'intensité située dans une même plage de valeurs. Un groupe de pixels 30 peut également résulter d'un maillage arbitraire de la surface de l'échantillon en différentes régions d'intérêt ROI. 3050528 22 La figure 3C illustre des fonctions d'écart A (1.1,,(À), u.',(À)), telles qu'explicitées selon l'équation (5), obtenues en utilisant un premier échantillon de test IL1-1,5%, dans trois longueurs d'onde : À = 500 nm, À =600 nm et À =700 nm, le coefficient µ'S(À) de l'échantillon de test ayant été mesuré, à ces trois longueurs d'onde. Le premier échantillon de test IL1-1.5% est un échantillon 5 comportant une concentration d'encre de Chine permettant d'obtenir un coefficient d'absorption µa(À) égal à 1 cm-1 à À=600 nm, et comportant une concentration volumique de 1.5% d'intralipide. Cette fonction d'écart ne comporte qu'une inconnue, à savoir u.a(À). Pour chaque longueur d'onde, la valeur µa(À) minimisant la fonction d'écart est le coefficient d'absorption de l'échantillon. Sur la figure 3C, on a représenté les valeurs estimées, 10 représentées par un cercle, ces dernières étant proches des valeurs réelles, représentées par un astérisque *. Dans ce qui précède, on estime une propriété d'absorption µa(À) (première propriété optique) à l'aide de la modalité d'imagerie spectrale de réflectance (MSI -1" modalité), tandis qu'on estime une propriété de diffusion µ'S(À) (propriété otique auxiliaire) à l'aide de spectrométrie 15 de réflectance diffuse (DRS - 2' modalité) appliquée sur l'échantillon. De façon alternative, on pourrait estimer une propriété de diffusion µ'S(À) (première propriété optique) par imagerie spectrale de réflectance et mesurer une propriété d'absorption µa(À) (propriété optique auxiliaire) par spectrométrie de réflectance diffuse. Cependant, la profondeur de champ de la DRS est plus élevée pour la mesure d'un coefficient de diffusion que pour la mesure d'un 20 coefficient d'absorption, cette profondeur de champ étant d'environ 2mm pour la mesure de µ'S(11) et d'environ 500p.m pour la mesure de lia (À), avec le dispositif de mesure par DRS décrit ci-avant. Ainsi, la mesure d'un coefficient d'absorption par DRS suppose une mise au point particulièrement précise du dispositif représenté sur la figure 1B. Cette contrainte de mise au point est moindre lorsqu'on utilise la DRS pour déterminer un coefficient de diffusion u.'s(À). 25 Par ailleurs, la profondeur de champ de l'imagerie spectrale de réflectance, appliquée à la détermination de µa(À) est plus élevée que la profondeur de champ de la DRS, avec le dispositif décrit ci-avant. Pour ces raisons de profondeur de champ, il est préférable de recourir à la DRS pour mesurer une propriété optique de diffusion u.'s(À), et à la MSI pour mesurer une propriété optique d'absorption lia (À). 30 Par ailleurs, quel que soit le mode de réalisation, la première propriété optique estimée par la première modalité, en différentes régions d'intérêt de l'échantillon, peut être recalée, dans au 3050528 23 moins une région d'intérêt, en utilisant une première propriété optique auxiliaire mesurée par la deuxième modalité. L'avantage de l'imagerie de réflectance est d'obtenir une estimation de la première propriété optique de l'échantillon en différentes zones d'intérêt ROI de la première image /1(À), chaque 5 région d'intérêt ROI étant définie par un pixel ou une pluralité de pixels (x, y) de la première image. Ainsi, lorsque la première propriété optique est un coefficient d'absorption, on obtient une distribution spatiale u.a_zy(.1) du coefficient d'absorption à la surface 5s de l'échantillon. La deuxième modalité permet de mesurer une propriété optique auxiliaire de l'échantillon, par exemple un coefficient de diffusion réduit u.',(.1) en une ou plusieurs éléments de surface 5, de 10 l'échantillon, comme décrit par la suite, en particulier en lien avec les figures 5B et 5C. Selon une variante, à l'issue de l'étape 150, l'étape 140 peut être réitérée, en utilisant, dans l'expression (4), un premier facteur de calibration tenant compte de la première propriété optique résultant de l'étape 150 pour la région d'intérêt considérée, et éventuellement de la propriété optique auxiliaire. Un tel facteur de calibration peut être obtenu à partir d'un 15 échantillon de calibration dont la première propriété optique correspond à la première propriété optique résultant de l'étape 150, ou par interpolation entre des facteurs de calibration obtenus respectivement à partir de différents échantillons, l'interpolation visant à estimer un facteur de calibration correspondant à la première propriété optique préalablement estimée, et/ou la propriété optique auxiliaire mesurée sur l'échantillon. Les étapes 150 et 140 peuvent 20 ainsi être réitérées jusqu'à l'atteinte d'un critère de convergence, par exemple un nombre d'itérations prédéterminé ou un faible écart entre deux premières propriétés optiques respectivement estimées au cours de deux itérations successives. Autrement dit, au cours de cette variante, on applique, lors de la première modalité, un premier facteur de calibration adapté à la première propriété optique de l'échantillon, et éventuellement à la propriété 25 optique auxiliaire de l'échantillon. Cela permet d'obtenir une meilleure estimation de la première propriété optique. Précisons qu'un tel procédé itératif ne nécessite pas une mesure de la propriété optique auxiliaire, mais peut se baser sur une estimation de cette dernière, sachant qu'il demeure préférable de procéder à une mesure de la propriété optique auxiliaire sur l'échantillon examiné. 30 Les inventeurs ont comparé l'estimation du coefficient d'absorption 14(.1), par imagerie spectrale de réflectance, en considérant soit une valeur arbitraire d'un coefficient de diffusion réduit, comme dans l'art antérieur, soit une valeur de ce coefficient mesurée par DRS, sur 3050528 24 l'échantillon examiné. Ces estimations ont été réalisées selon trois longueurs d'onde, respectivement 500 nm, 600 nm et 700 nm. La figure 4A représente les estimations de µa(À) des quatre échantillons test 1%IL, 1,5%IL, 2%IL, 3%IL, dont la concentration volumique en intralipide est respectivement de 1%, 1,5%, 2% et 3%. 5 Ces échantillons test présentent un même coefficient d'absorption µ,(À), ce dernier étant représenté en pointillés sur les figures 4A et 4C. Sa valeur s'élève à 0.4 cm-1 à 600 nm. Chaque échantillon a son propre coefficient de diffusion réduit, de telle sorte que le pourcentage d'intralipide correspond à une valeur d'un coefficient de diffusion réduit. Ces estimations sont effectuées en se basant sur des images de chaque échantillon de test, telles que représentées 10 sur la figure 3A. Sur cette dernière, on a représenté des régions d'intérêt identifiées par un pourcentage, chaque région d'intérêt correspondant à un échantillon test. Chaque région d'intérêt a été traitée en considérant une mesure du coefficient de diffusion réduit µ'S(À) effectuée sur l'échantillon de test 1.5%IL, ce dernier étant, dans cet exemple, comme représentatif des autres échantillons de test. On constate que seule l'estimation du coefficient 15 d'absorption µa(À) de l'échantillon 1.5%IL est fiable. Les estimations liées aux autres échantillons sont affectées d'une erreur de mesure importante, bien qu'étant basées sur le coefficient de diffusion réduit 1.1's (À) de l'échantillon de test 1.5%IL, dont les compositions sont pourtant proches des autres échantillons de test. La figure 4B illustre l'erreur relative moyenne des estimations du coefficient d'absorption µa(À) 20 dans les trois longueurs d'onde, représentées sur la figure 4A, pour chaque échantillon de test. Seule l'estimation concernant l'échantillon 1.5%IL est satisfaisante. Afin d'obtenir une estimation fiable du coefficient d'absorption, il est nécessaire de se baser sur une mesure du coefficient de diffusion réduit µ'S(W) sur l'échantillon examiné, et non sur d'autres estimations, même lorsque ces dernières sont obtenues à partir d'échantillons dont les 25 compositions comme proches de l'échantillon examiné. La figure 4C représente les estimations de µa(À) obtenues sur chaque échantillon test 1%IL, 1,5%IL, 2%IL, 3%IL en considérant 4 valeurs différentes du coefficient de diffusion réduit Vs(À), ces valeurs ayant été respectivement mesurées sur chacun de ces échantillons de test. Le coefficient de diffusion réduit µ'S(À) de chaque échantillon test a été estimé par DRS aux trois longueurs d'onde mentionnées ci-dessus. 30 On a ensuite estimé le coefficient d'absorption µa(À) de chaque échantillon de test par imagerie spectrale de réflectance, en prenant en compte le coefficient de diffusion réduit mesuré sur 3050528 25 chacun d'entre eux. On observe que les valeurs estimées de µa(À) sont nettement plus proches de la valeur réelle. La figure 4D illustre l'erreur moyenne de l'estimation du coefficient d'absorption µa(11.), dans les trois longueurs d'onde. On observe que ces erreurs sont moindres que celles représentées sur 5 la figure 4B. La mesure du coefficient d'absorption 11,(11.) d'un échantillon, basée sur une mesure du coefficient de diffusion réduit 11;(g) mesuré sur ledit échantillon, est donc plus fiable. Les inventeurs ont réalisé une autre série de mesures en mesurant la propriété optique auxiliaire, c'est-à-dire le coefficient de diffusion réduit, de façon itérative, le deuxième facteur de calibration prenant en compte ledit coefficient mesuré lors d'une itération précédente. Cette 10 méthode d'adaptation du facteur de calibration, précédemment décrite, permet d'adapter le deuxième facteur de calibration aux propriétés optiques de l'échantillon, et résulte en une meilleure estimation du coefficient de diffusion réduit. De ce fait, l'erreur de mesure de la première propriété optique, c'est-à-dire le coefficient d'absorption, est systématiquement inférieure à 5%. La précision de la mesure est donc améliorée. 15 Un des points clefs de l'approche proposée par les inventeurs est d'effectuer une mesure locale d'une propriété optique auxiliaire afin de mieux estimer une première propriété optique à partir de la première image /1(À). La propriété optique auxiliaire est estimée ponctuellement, au niveau d'un élément de surface 5i de la surface 5s de l'échantillon. La localisation de cet élément de surface 5i peut être déterminée à partir de la première image /1(À). Cette première image 20 peut par exemple être segmentée en zones présentant des pixels Ii (xy) de même intensité, ou d'intensité comparable, la propriété optique auxiliaire étant estimée sur une ou plusieurs zones issues de cette segmentation. Ainsi, la première image /1(.1) peut permettre de définir la localisation de l'élément de surface 5i ou un maillage d'éléments de surface 5i sur lesquels la propriété optique auxiliaire est mesurée. La propriété optique auxiliaire peut être mesurée 25 successivement, sur chaque élément de surface ainsi défini, en appliquant les étapes 200 à 240 précédemment décrites, le dispositif représenté sur la figure 1B effectuant un balayage à la surface de l'échantillon. Selon une variante, les éléments de surface 5i sont déterminés arbitrairement, par exemple sous la forme d'un maillage, régulier ou non, établi à la surface de l'échantillon. La propriété optique 30 auxiliaire est alors mesurée en chaque point de ce maillage. L'application de la deuxième modalité de mesure sur différents éléments de surface 5i permet une mesure discrète de la propriété optique auxiliaire, selon chaque élément de surface réparti 3050528 26 discrètement à la surface de l'échantillon, selon un maillage. Il est alors possible d'interpoler les différentes mesures, de façon à estimer la propriété optique auxiliaire entre chaque maille. L'interpolation peut être linéaire ou polynomiale. Ainsi, la première propriété optique, au niveau d'une région d'intérêt ROI de la première image 11(À) peut être déterminée à partir d'une 5 propriété optique auxiliaire, soit mesurée, au niveau cette région d'intérêt, lorsque cette dernière coïncide avec un élément de surface 5i, soit interpolée à partir de propriétés optiques auxiliaires respectivement mesurées en différents éléments de surface de l'échantillon, en particulier adjacents de ladite région d'intérêt. Les figures 5A à 5C représentent des variantes du dispositif 1. Sur la figure 5A, on a représenté 10 la deuxième composante lb permettant d'effectuer des mesures par DRS. Sur cette figure, par rapport à la figure 1B, on a supprimé les fibres optiques de détection 27 ainsi que le commutateur optique 29. Le photodétecteur auxiliaire 28 est un capteur d'image auxiliaire, ou deuxième capteur d'image. Selon cette configuration, le photodétecteur auxiliaire n'est pas spectralement résolu, la bande spectrale du faisceau lumineux auxiliaire 22 étant modulable. La 15 source de lumière auxiliaire 20 peut être associée à des filtres optiques, par exemple sous la forme d'une roue à filtre 35. De façon alternative, la source de lumière auxiliaire comporte des sources élémentaires de lumières, par exemple des diodes laser ou des diodes électroluminescentes, émettant dans différentes bandes spectrales. Le capteur d'image auxiliaire 28 est optiquement couplé à la surface 5s de l'échantillon, de telle sorte que les zones 20 élémentaires de rétrodiffusion 25 soient distinctes les unes des autres, et distinctes de la zone élémentaire d'illumination 23. La source de lumière auxiliaire peut émettre un faisceau lumineux auxiliaire collimaté, formant la zone d'illumination 23, auquel cas la fibre optique d'illumination 21 et son optique 26 peut être supprimée. De même que dans la figure 1B, une fibre de retour d'excitation 21e peut guider le faisceau lumineux auxiliaire 22 émis vers le 25 photodétecteur auxiliaire 28. La figure 5B représente un dispositif 1 comportant les premières et deuxièmes composantes la et lb. Sur ce dispositif, le premier capteur d'image 18 et le capteur d'image auxiliaire 28 sont confondus et forment un même capteur d'image, permettant d'acquérir : une première image 11(À) de l'échantillon 5 lorsque ce dernier est illuminé par un 30 premier faisceau lumineux 12, ce dernier couvrant une surface étendue de l'échantillon ; 3050528 27 une deuxième image 12(À) de l'échantillon 5 lorsque ce dernier est illuminé par un faisceau lumineux auxiliaire 22, ce dernier formant une zone élémentaire d'illumination 23 sur la surface de l'échantillon. Le premier faisceau 12 et le faisceau auxiliaire 12, 22 sont respectivement émis par la première 5 source de lumière 10 et par la source de lumière auxiliaire 20 dans une bande spectrale, et éventuellement successivement dans différentes bandes spectrales. L'optique 17 peut être adaptée entre l'acquisition de la première image et de la deuxième image. De même que dans les figures 1B et 5A, une fibre de retour d'excitation 21e peut guider le faisceau lumineux auxiliaire 22 émis vers le photodétecteur auxiliaire 28. 10 La première image /1(À) est traitée par le microprocesseur 30, de façon à estimer une première propriété optique, selon les étapes 100 à 150 préalablement décrites. La deuxième image 12(À) est traitée par le microprocesseur 30 de façon à déterminer, sur cette image, l'intensité de signaux rétrodiffusés Sn(A) à différentes distances de rétrodiffusion On puis à mettre en oeuvre le procédé décrit en lien avec les étapes 200 à 240 pour mesurer une propriété optique 15 auxiliaire de l'échantillon, cette dernière étant utilisée, lors de l'étape 150, pour estimer la première propriété optique. Dans les modes de réalisation précédemment décrits, la première image /1(À) permet une estimation d'une propriété optique à partir de différentes régions d'intérêt ROI de la première image, chaque région d'intérêt correspondant à un pixel /1,x3,(À), ou à un groupe de pixels, de 20 ladite image. La deuxième modalité permet d'estimer ponctuellement, c'est-à-dire au niveau d'un élément de surface 5i, la propriété optique auxiliaire. Le faisceau lumineux auxiliaire 22 peut être successivement dirigé vers différents éléments de surface 5i, de façon à estimer successivement la propriété optique auxiliaire de l'échantillon au niveau de chaque élément de surface 5i. 25 Le mode de réalisation présenté sur la figure 5C permet une formation simultanée de différentes zones élémentaires d'illumination 23 réparties à la surface de l'échantillon. Pour cela, on a disposé un masque 36 entre la source de lumière auxiliaire 20 et la surface de l'échantillon 5. Un tel masque est constitué de portions transparentes 361 et de portions opaques 362, ces portions étant projetées à la surface de l'échantillon 5. La figure 5D montre un exemple d'un tel 30 masque. On obtient alors autant de zones d'illumination élémentaires 231,232, 233 que de portions transparentes projetées. On peut disposer un système optique 26 pour focaliser le faisceau lumineux auxiliaire 22 sur la surface de l'échantillon. Une fibre de retour d'excitation 3050528 28 21e peut relier la source de lumière auxiliaire 20 au capteur d'image. Le capteur d'image auxiliaire 28, qui est ici confondu avec le premier capteur d'image 18, forme une deuxième image 12(À), sur laquelle apparaissent différentes zones élémentaires d'illumination 23 et, les différentes zones de rétrodiffusion élémentaires 25 associées à chaque zone élémentaire 5 d'illumination. Les étapes 200 à 240 sont mises en oeuvre à partir de l'intensité des zones de rétrodiffusion élémentaires 25 associées à une même zone élémentaire d'illumination 23. La propriété optique auxiliaire est estimée, en différentes localisations discrètes 5i, à partir de la même image 12 (À). Ce dispositif évite un balayage du faisceau lumineux auxiliaire 20 à la surface de l'échantillon, puisque le masque permet une formation simultanée de plusieurs faisceaux 10 auxiliaires 221, 222, 223 Le masque peut être constitué d'un arrangement régulier de portions opaques et de portions transparentes. Il peut s'agir d'une grille ou d'une plaque de verre sur laquelle un motif opaque a été formé. Il peut également s'agir d'un écran à cristaux liquides, ce qui permet de former un masque dont le motif opaque est modulable. La figure 6 schématise un mode de réalisation de l'invention, appliqué sur la surface d'un 15 échantillon. Une première image 11(À) est formée, à partir de laquelle on détermine des éléments de surface sur lesquels on mesure, de façon discrète, la propriété optique auxiliaire, en l'occurrence le coefficient de diffusion réduit. L'interpolation des mesures permet d'obtenir une répartition spatiale de la propriété optique auxiliaire us'(11). Cette répartition spatiale est utilisée pour estimer le coefficient d'absorption u.a(11), et d'obtenir une répartition 20 spatiale u.a,xy(À)de ce dernier. On dispose alors d'une estimation quantitative de la distribution spatiale des propriétés optiques d'absorption et de diffusion à la surface de l'échantillon examiné. La méthode et le dispositif objets de l'invention pourront être mis en oeuvre sur des échantillons biologiques, par exemple des tissus biologiques, à des fins d'aide au diagnostic, notamment dans 25 le cas de pathologies de la peau. Elle peut également être mise en oeuvre dans le suivi de la vascularisation de la peau, par exemple suite à des greffes, ou à l'analyse d'échantillons de peau prélevés au cours d'une opération, en particulier pour délimiter des zones tumorales de zones saines. Au-delà des échantillons biologiques, la méthode pourra être mise en oeuvre pour contrôler des 30 échantillons, par exemple dans l'industrie agroalimentaire ou le contrôle de l'environnement. REVENDICATIONS1. Procédé d'estimation de propriétés optiques d'un échantillon, comportant les étapes suivantes : a) illumination de la surface de l'échantillon à l'aide d'un premier faisceau lumineux (12), produit par une première source de lumière (10) ; b) acquisition, à l'aide d'un premier capteur d'image (18), d'une première image (11(.1)) d'un rayonnement lumineux (14) rétrodiffusé par l'échantillon ainsi illuminé ; c) à l'aide de ladite première image, détermination, en plusieurs régions d'intérêts (ROI) de la surface de l'échantillon (5s), d'une première grandeur d'intérêt (R(.1)), représentative d'une quantité dudit rayonnement lumineux rétrodiffusé (14) par l'échantillon ; d) estimation d'une première propriété optique (14(.1)) en chacune desdites régions d'intérêt (ROI), à l'aide de ladite première grandeur d'intérêt ; le procédé étant caractérisé en ce que préalablement à l'étape d), le procédé comporte une étape de mesure d'une propriété optique dite auxiliaire (u.'s(.1)) de l'échantillon, différente de ladite première propriété optique, chaque estimation réalisée lors de l'étape d) prenant en compte la propriété optique auxiliaire ainsi mesurée sur l'échantillon. 2. Procédé selon la revendication 1, dans lequel : soit la première propriété optique est une propriété d'absorption de la lumière (u.a(11)) de l'échantillon, la propriété optique auxiliaire étant une propriété de diffusion de la lumière (Vs(.1)) dans l'échantillon ; soit la première propriété optique est une propriété de diffusion de la lumière (Vs(.1)) de l'échantillon, la propriété optique auxiliaire étant une propriété d'absorption de la lumière (14(.1)) dans l'échantillon. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape de mesure de la propriété optique auxiliaire (u.'s(.1)) comporte les sous-étapes suivantes : i) illumination de la surface de l'échantillon à l'aide d'un faisceau lumineux auxiliaire (22), produit par une source de lumière auxiliaire (20), de façon à former, à la surface de l'échantillon, une zone élémentaire d'illumination (23), correspondant à la partie de la surface illuminée par ledit faisceau de lumineux auxiliaire ; 3050528 30 ii) mesure, à l'aide d'un photodétecteur auxiliaire (28), d'un signal de rétrodiffusion (4(.1)), représentatif d'un rayonnement rétrodiffusé (24) par l'échantillon, ledit rayonnement rétrodiffusé émanant d'une zone élémentaire de rétrodiffusion (25), séparée de la zone élémentaire d'illumination (23) et située à une distance (Dn), dite 5 distance de rétrodiffusion, de cette dernière ; iii) estimation de la propriété optique auxiliaire de l'échantillon en fonction du signal de rétrodiffusion (S,,(À)) mesuré lors de la sous-étape ii). 4. Procédé selon la revendication précédente dans lequel : la sous-étape ii) comporte l'obtention, à partir du signal de rétrodiffusion (Sn(À)), d'une 10 grandeur d'intérêt dite auxiliaire (Rn(À)), représentative d'une quantité du rayonnement (24) rétrodiffusé par l'échantillon, à une distance de rétrodiffusion (Dn), suite à son illumination par le faisceau lumineux auxiliaire; la sous-étape iii) comporte une comparaison de ladite grandeur d'intérêt auxiliaire (Rn(À)) avec une pluralité d'estimations (RZades,(À)) de cette grandeur d'intérêt, 15 chaque estimation étant réalisée en considérant une valeur prédéterminée de ladite propriété optique auxiliaire (1.1',(À)). 5. Procédé selon l'une quelconque des revendications 3 ou 4 dans lequel : lors de la sous-étape i), l'échantillon est illuminé successivement ou simultanément par une pluralité de faisceaux lumineux auxiliaires (221, 222, 223) de façon à former, à la 20 surface de ce dernier, une pluralité de zones élémentaires d'illumination (231, 232, 233) espacées les unes des autres ; les sous-étapes ii) et iii) sont mises en oeuvre en considérant un signal rétrodiffusé émis à une distance de rétrodiffusion (Dn) de chaque zone élémentaire d'illumination, de façon à estimer la propriété optique auxiliaire dans autant d'éléments de surface (5i) 25 que de faisceaux lumineux auxiliaires. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape d) comporte, dans chaque région d'intérêt (ROI), une comparaison de la première grandeur d'intérêt (R(À)) déterminée lors de l'étape c) avec une pluralité d'estimations (R(T,i'adaf,ls)(À)) de cette grandeur d'intérêt, chaque estimation étant réalisée en considérant une valeur 30 prédéterminée de la première propriété optique (ua(À)), et en prenant en compte la propriété optique auxiliaire (u.',(À)) préalablement mesurée sur l'échantillon. 3050528 31 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel la propriété optique auxiliaire (12,(.1)) de l'échantillon est mesurée en une pluralité d'éléments de surface (5i) de l'échantillon. 8. Procédé selon la revendication 7, comportant une interpolation de ladite propriété optique 5 auxiliaire (12,(.1)) à partir d'au moins deux propriétés optiques auxiliaires préalablement mesurées, de façon à établir une distribution spatiale (11'sxy (11)) de ladite propriété optique auxiliaire. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel la propriété optique auxiliaire (u.'s(.1)) de l'échantillon est mesurée selon au moins un élément de surface 10 (5i), ledit élément de surface étant déterminé à partir de la première image (4(.1)) acquise lors de l'étape b). 10. Dispositif (1) pour estimer des propriétés optiques d'un échantillon (5) comportant : une première source de lumière (10), apte à émettre un premier faisceau lumineux (12) se propageant vers l'échantillon ; 15 un premier capteur d'image (18), apte à acquérir une première image (4(.1)) d'un premier rayonnement lumineux (14) rétrodiffusé par l'échantillon ; une source de lumière auxiliaire (20), apte à émettre un faisceau lumineux auxiliaire (22) se propageant vers l'échantillon, de façon à former, à la surface de ce dernier, une zone élémentaire d'illumination (23), correspondant à la partie d'une surface (5s) de 20 l'échantillon illuminée par ledit faisceau lumineux auxiliaire ; un photodétecteur auxiliaire (28) apte à mesurer un rayonnement rétrodiffusé (24) par l'échantillon, à une distance de rétrodiffusion (Dn) de la zone élémentaire d'illumination (23), suite à l'illumination de l'échantillon par la source de lumière auxiliaire (20) ; un microprocesseur (30), configuré pour déterminer une propriété optique auxiliaire 25 (u.'s(.1)) de l'échantillon à partir du rayonnement rétrodiffusé (24) mesuré par le photodétecteur auxiliaire (28), le microprocesseur (30) étant également configuré pour déterminer une première propriété optique (u.a(11)) différente de la propriété optique auxiliaire, en différents points de la surface (5s) de l'échantillon, à partir de la première image acquise par le premier capteur d'image (18), en prenant en compte ladite 30 propriété optique auxiliaire. 11. Dispositif selon la revendication 10, dans lequel le photodétecteur auxiliaire (28) est un capteur d'image. 3050528 32 12. Dispositif selon la revendication 11, dans lequel le photodétecteur auxiliaire (28) est le premier capteur d'image (18). 13. Dispositif selon l'une quelconque des revendications 10 à 12, dans lequel la source de lumière auxiliaire (20) est apte à former simultanément une pluralité de zones élémentaires 5 d'illumination (231, 232, 233), espacées les unes des autres, à la surface de l'échantillon (5). 14. Dispositif selon la revendication 13, comportant un masque (36), interposé entre la source de lumière auxiliaire (20) et l'échantillon, le masque comportant des portions transparentes (36a), permettant le passage du faisceau lumineux auxiliaire (22), et des portions opaques (36b), bloquant le faisceau lumineux auxiliaire, de telle sorte que chaque zone élémentaire 10 d'illumination (231, 232, 233) est formée par une projection d'une portion transparente (36a) sur la surface de l'échantillon. 15. Dispositif selon l'une quelconque des revendications 10 à 14, dans lequel le microprocesseur est apte à mettre en oeuvre les étapes c) et d) du procédé objet de l'une quelconque des revendications 1 à 9. 15