FR 3041756 A1 20170331 FR 1559287 A 20150930 1 Procédé de détermination de la réflectance d'un objet et dispositif associé La présente invention concerne un procédé de détermination de la réflectance d'un objet et un dispositif de détermination de la réflectance d'un objet. Le document WO 2013/120956 Al décrit un procédé de mesure de la réflectance diffuse uniforme en au moins un point d'un objet à l'aide d'un dispositif comportant un moyen apte à émettre des illuminants de couleur exprimés sous la forme de flux lumineux et un capteur électronique d'images couleurs et un dispositif comportant un moyen apte à émettre des illuminants de couleur sous la forme de flux lumineux de couleurs et un capteur électronique d'images couleurs, de mesure de la réflectance diffuse uniforme en au moins un point d'un objet placé dans une zone se trouvant en face de et sensiblement perpendiculairement audit moyen apte à émettre des couleurs et se trouvant dans le champ de vision dudit capteur électronique d'images couleurs et étant soumis à un illuminant extérieur sous la forme d'un flux lumineux extérieur environnant constant et inconnu. Pour déterminer la réflectance d'un objet, il est connu d'utiliser des dispositifs spécialisés de haute précision comme des spectromètres à diffraction ou des spectromètres à double-photorécepteurs en parallèle. Cependant, de tels dispositifs sont coûteux et difficiles à utiliser pour des opérateurs non spécialisés. Il existe donc un besoin pour un procédé de détermination de la réflectance fiable et facile à mettre en oeuvre. Pour cela, il est proposé un procédé de détermination de la réflectance d'un objet, le procédé comprenant les étapes de l'éclairement de l'objet par un illuminant extérieur présentant un éclairement inconnu et variable, l'émission d'au moins un flash de lumière éclairant l'objet, chaque flash de lumière étant émis par une source et présentant un éclairement connu dans une gamme de longueurs d'onde, collection de l'onde réfléchie par l'objet pour former au moins une image sur un capteur, l'obtention d'une équation à plusieurs inconnues, l'équation étant obtenue à partir des images formées, la réflectance de l'objet et l'éclairement de l'illuminant extérieur étant deux inconnues de l'équation et la résolution de l'équation. L'étape de résolution de l'équation comprendt le calcul de points solution de l'équation, l'interpolation des points calculés par une fonction d'interpolation, et l'utilisation d'au moins une des approximations suivantes pour la résolution de l'équation : o une première approximation selon laquelle chaque image est issue de l'émission d'un flash de lumière distinct, o une deuxième approximation selon laquelle la fonction d'interpolation détermine les points de stabilité de l'équation. 3041756 2 Un tel procédé de détermination de la réflectance p d'un objet est facile à mettre en oeuvre et permet d'obtenir une modélisation fiable de la réflectance p réelle de l'objet même avec un illuminant extérieur variable. Une telle mise en oeuvre permet de réduire le temps de calcul, tout en conservant la précision du procédé. 5 Suivant des modes de réalisation particuliers, le procédé de détermination de la réflectance d'un objet comprend une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - la source et le capteur sont disposés sur un même appareil. - une pluralité de flashs de lumière sont émis, chaque flash présentant un maximum 10 d'éclairement en longueur d'onde, l'étape de collection étant mise en oeuvre pour chaque flash de lumière émis et au moins deux flashs de lumière présentent un maximum d'éclairement éloigné d'au moins 20 nanomètres. - l'étape de collection est mise en oeuvre plusieurs fois pour un même flash de lumière, l'équation obtenue étant un système d'équations surdéterminé, l'étape de 15 résolution étant mise en oeuvre pour une pluralité de systèmes d'équations déterminés en utilisant la première approximation pour obtenir une pluralité de fonctions de réflectance, le procédé comprenant, en outre, un calcul de la réflectance de l'objet en calculant une moyenne de la pluralité des fonctions de réflectance. 20 - la deuxième approximation est utilisée lors de l'étape de résolution de l'équation et dans lequel la fonction d'interpolation est une combinaison pondérée de fonctions de bases scellée par un nombre fini de points d'interpolation, notamment des splines cubiques, chaque point d'interpolation étant un point de stabilité de l'équation. 25 - une pluralité de flashs de lumière sont émis, chaque flash présentant un maximum d'éclairement en longueur d'onde, l'étape de collection étant mise en oeuvre pour chaque flash de lumière émis. Les points d'interpolations vérifient au moins la propriété suivante : le nombre de point d'interpolation est égal au nombre de flashs.-le procédé est un procédé de mesure de la réflectance diffuse uniforme RoeJ(X) n au 30 moins un point d'un objet à l'aide d'un dispositif comportant un moyen apte émettre des illuminants de couleur exprimés sous la forme de flux lumineux et un capteur électronique d'images couleurs, caractérisé en ce qu'il comporte les étapes suivantes : - placement dudit objet dans une zone se trouvant en face de et sensiblement 35 perpendiculairement audit moyen apte à émettre des illuminants de couleur sous la forme de flux lumineux de couleurs et se trouvant dans le champ de 3041756 3 vision dudit capteur électronique d'images couleurs, ledit objet étant également soumis à un illuminant extérieur sous la forme d'un flux lumineux extérieur environnant lext(X) inconnu et constant où X désigne la longueur d'onde, émission par ledit moyen d'une série de N illuminants S SOURCE(X)i (avec N entier naturel supérieur à un, i variant de 1 à N et X, la longueur d'onde), SSOURCE', (A) étant connu en fonction des paramètres d'entrée dudit moyen apte à émettre des flux lumineux de couleurs, capture par ledit capteur électronique d'images couleurs du flux lumineux réfléchi en au moins un point dudit objet et entrant dans le capteur, ledit flux lumineux étant noté Ecapteur,, (A) avec N entier naturel strictement supérieur à deux, i variant de 1 à N et X la longueur d'onde, et obtention de N équations « E, »: Ecapteur(x)i = ROBJ(x)*(lext(x)+sSOURCE(x)i) en raison du caractère additif de la lumière ondulatoire et par définition de la réflectance diffuse uniforme RoeJ(X) en au moins un point de l'objet (30) ; et 15 - détermination par ledit dispositif des deux fonctions inconnues continues Row(x) et lext(X) en résolvant le système des N équations E,: en intégrant chaque équation E, sur l'intersection des spectres source et capteur, en notant ID, chaque sensibilité dans la base colorimétrique choisie, chaque équation E, générant alors un ensemble d'équations 20 « E, intégrée » : Ecapteur (2)i * (2) * _ f RoBJ (2) * uext (2) sSOURCE (2)i) * x(2) * en calculant la valeur numérique correspondant au terme de gauche des équations E, intégrée à l'aide des paramètres de sortie du capteur numérique d'image ; et en exprimant les deux fonctions inconnues continues RoBJ(X) et lext(X) 25 à l'aide d'un nombre fini de points d'interpolation (X,, y,) reliés par au moins une fonction d'interpolation s(X) pour conserver le caractère continu desdites fonctions inconnues continues RoeJ(X) et lext(X), les xi étant des longueurs d'onde choisies dans l'intersection des spectres source et capteur et étant des paramètres d'entrée du procédé, 30 choisis pour minimiser le nombre de points d'interpolation à précision donnée ; et 5 10 3041756 4 en recherchant les paramètres yi des fonctions R°I3J(X) et lext(X) qui minimisent le système de moindres carrés IIA * X - BII2 résultant des équations E, intégrées. - le procédé comporte en outre une étape de détermination de la valeur de 5 l'illuminant extérieur lext(X). - le procédé comporte en outre une étape de transcription de la fonction RoeJ(X) de réflectance diffuse uniforme en au moins un point de l'objet en coordonnées CIE XYZ pour un illuminant donné. - le nombre de flashs est du même ordre de grandeur que le nombre de points 10 d'interpolation pour déterminer les valeurs de la réflectance diffuse uniforme RoeJ(X) en au moins un point de l'objet et de l'illuminant extérieur lext(X). - le procédé comporte une étape de détermination des valeurs de la réflectance diffuse uniforme leBJ(X) en au moins un point de l'objet et de l'illuminant extérieur lext(X) dans plusieurs bandes spectrales. 15 - ledit dispositif met en oeuvre un écran pour émettre les flashs de couleur et un capteur électronique d'image pour capter la lumière réfléchie par l'objet-cible. - ledit dispositif est un appareil photo ou une caméra avec flash intégré ou amovible. - ledit dispositif met en oeuvre des guides d'onde pour faire transiter l'émission et la réception des flashs de couleurs. 20 - ledit dispositif est mis en oeuvre pour réaliser des photographies spectrométriques d'objets et pour réaliser des adaptations chromatiques (balance des blancs) à volonté. - ledit dispositif est mis en oeuvre pour mesurer la couleur d'un élément compris dans le groupe suivant : matériaux, solides, liquides, gaz, peintures, tapisseries, 25 graphismes, textiles, plastiques, bois, métaux, sols, minéraux, végétaux et aliments. - ledit dispositif est mis en oeuvre pour la mesure des couleurs à but médical ou cosmétique sur l'homme et le vivant d'au moins un élément compris dans le groupe suivant : peau, boutons, grains de beauté, chevelure, pelage, maquillage, et dents. - ledit dispositif est mis en oeuvre pour l'utilisation de codes-barres en couleur, à une 30 ou plusieurs dimensions, et - ledit dispositif est mis en oeuvre dans un but d'assistance aux personnes daltoniennes et/ou aveugles. La présente description concerne aussi le dispositif de détermination de la réflectance d'un objet, l'objet étant éclairé par un illuminant extérieur présentant un 35 éclairement inconnu et variable, le dispositif comprenant une source, apte à émettre au moins un flash de lumière éclairant l'objet, chaque flash de lumière émis par la source 3041756 5 présentant un éclairement connu dans une gamme de longueurs d'onde, un capteur, apte à collecter l'onde réfléchie par l'objet pour former au moins une image. Une unité de traitement, propre mettre en oeuvre les étapes suivantes : obtention d'une équation à plusieurs inconnues, l'équation étant obtenue à partir des images formées, la réflectance 5 de l'objet et l'éclairement de l'illuminant extérieur étant deux inconnues de l'équation, et résolution de l'équation. L'étape de résolution de l'équation comprenant le calcul de points solution de l'équation, l'interpolation des points calculés par une fonction d'interpolation, et l'utilisation d'au moins une des approximations suivantes pour la résolution de l'équation : une 10 première approximation selon laquelle chaque image est issue de l'émission d'un flash de lumière distinct, une deuxième approximation selon laquelle la fonction d'interpolation détermine les points de stabilité de l'équation. Suivant des modes de réalisation particuliers, le dispositif de détermination de la réflectance d'un objet comprend une ou plusieurs des caractéristiques suivantes, prise(s) 15 isolément ou suivant toute combinaison techniquement possible : - le capteur et la source sont disposés sur un même appareil. - la source est un écran lumineux ou un ensemble de diodes électroluminescentes. - le capteur est choisi dans un groupe constitué d'un appareil photo, d'une caméra, d'un imageur multicanal et d'un imageur hyperspectral.- le dispositif comporte en 20 outre un moyen apte à émettre des illuminants de couleur sous la forme de flux lumineux de couleurs et un capteur électronique d'images couleurs, de mesure de la réflectance diffuse uniforme ROBJ(X) en au moins un point d'un objet placé dans une zone se trouvant en face de et sensiblement perpendiculairement audit moyen apte à émettre des couleurs et se trouvant dans le champ de vision dudit capteur 25 électronique d'images couleurs et étant soumis à un illuminant extérieur sous la forme d'un flux lumineux extérieur environnant noté lext(X) constant et inconnu. Le dispositif comporte également des moyens pour : - émettre une série de N illuminants SsOURCE,, (A) (avec N entier naturel supérieur à un, i variant de 1 à N et X, la longueur d'onde), SsOURCE,, (A) étant 30 connu en fonction des paramètres d'entrée dudit moyen apte à émettre des flux lumineux de couleurs, capture par ledit capteur électronique d'image couleur du flux lumineux réfléchi en au moins un point dudit objet et entrant dans le capteur, ledit flux lumineux étant noté Ecapteur) avec N entier naturel strictement supérieur à deux, i variant de 1 à N et X la longueur 35 d'onde, et obtention de N équations « E, »: Ecapteur(x)i = ROBJ(x)*(lext(x)+sSOURCE(x)i) 3041756 6 en raison du caractère additif de la lumière ondulatoire et par définition de la réflectance diffuse uniforme RoeJ(X) en au moins un point de l'objet ; et - déterminer les deux fonctions inconnues continues RoeJ(X) et leXt(X) en résolvant le système des N équations Ei: 5 en intégrant chaque équation E, sur l'intersection des spectres source et capteur, en notant ID, chaque sensibilité dans la base colorimétrique choisie, chaque équation E, générant alors un ensemble d'équations « E, intégrée » : f E capteur 7 Ut) * bj (2) * d2 = f R°13-I(2)(* ext (2) sSOURCE (2)i) * x(2) * dA en calculant la valeur numérique correspondant au terme de gauche 10 des équations E, intégrée à l'aide des paramètres de sortie du capteur numérique d'image ; et en exprimant les deux fonctions inconnues continues RoeJ(X) et leXt(X) à l'aide d'un nombre fini de points d'interpolation (XJ, y,) reliés par au moins une fonction d'interpolation s(X) pour conserver le caractère 15 continu desdites fonctions inconnues continues RoeJ(X) et lext(X) les xi étant des longueurs d'onde choisies dans l'intersection des spectres source et capteur et étant des paramètres d'entrée du procédé, choisis pour minimiser le nombre de points d'interpolation à précision donnée ; et 20 en recherchant les paramètres yi des fonctions RoeJ(X) et lext(X) qui minimisent le système de moindres carrés IIA * X - BII2 résultant des équations E, intégrées. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit de modes de réalisation de l'invention, donnée à titre d'exemple 25 uniquement et en référence aux dessins qui sont : - figure 1, un schéma d'un dispositif de détermination de la réflectance d'un objet, - figure 2, un ordinogramme d'un premier exemple de mise en oeuvre de procédé de détermination de la réflectance d'un objet, - figure 3, une représentation graphique d'un spectre pour plusieurs flashs de 30 lumière, - figure 4, un ordinogramme d'un deuxième exemple de mise en oeuvre de procédé de détermination de la réflectance d'un objet, - figure 5, une représentation graphique de l'erreur de la réflectance déterminée par rapport à une réflectance réelle, pour plusieurs déterminations, 3041756 7 - figure 6, une représentation graphique de l'erreur de la réflectance déterminée par rapport à la réflectance réelle, - figure 7, une représentation graphique de plusieurs réflectances déterminées et de la réflectance réelle, 5 - figure 8, une représentation graphique de la réflectance déterminée par la mise en oeuvre d'un procédé de l'état de la technique et de la réflectance réelle, - figure 9, une représentation graphique de la réflectance déterminée pour un exemple de mise en oeuvre de procédé de détermination et de la réflectance réelle. Un dispositif 1 de détermination de la réflectance d'un objet 4, un objet 4 et un 10 illuminant extérieur 6 sont représentés sur la figure 1. Le dispositif 1 de détermination de la réflectance d'un objet 4 comprend une source 10, un capteur 12 et une unité de traitement 14. Le dispositif 1 de détermination est propre à mettre en oeuvre un procédé de détermination de la réflectance d'un objet 4, dont un exemple de mise en oeuvre est 15 représenté sur la figure 2. La réflectance réelle d'un objet, notée réelle, est une fonction (À) de la longueur réelle, d'onde notée À. La réflectance donne une information sur la couleur de l'objet 4 au sens chromatique. La réflectance de l'objet 4 dépend de la matière de la surface de réflexion de 20 l'objet 4. La réflectance d'un objet 4 est définie comme le rapport entre la luminance reçue par l'objet 4 et la luminance réfléchie par l'objet 4. La luminance est une grandeur correspondant à la sensation visuelle de luminosité d'une surface. La luminance reçue par l'objet 4 se définit comme le quotient de l'intensité 25 lumineuse reçue par l'objet 4 par l'aire de la surface apparente de l'objet 4. La luminance réfléchie par l'objet 4 se définit comme le quotient de l'intensité lumineuse réfléchie par l'objet 4 par l'aire de la surface apparente de l'objet 4. La surface apparente est la projection de l'aire de l'objet 4 perpendiculairement à une direction d'observation. L'éclairement de l'objet 4 est connu à partir de la luminance reçue par l'objet 4 et de 30 la géométrie d'observation. Dans le procédé de détermination, la surface de l'objet 4 est considérée comme une surface Lambertienne. Une surface Lambertienne est une surface où la luminance est indépendante de la direction d'observation. L'éclairement de l'objet 4 correspond à un flux lumineux reçu par unité de surface. 35 La détermination de la réflectance consiste à trouver une réflectance déterminée notée p, la plus proche possible selon une norme, de la réflectance réelled'un objet 4 'réelle 3041756 8 sur une gamme de longueurs d'onde. La gamme de longueurs d'onde dépend de la source 10 et du capteur 12 du dispositif 1 de détermination. Par exemple, la réflectance est déterminée sur une gamme de longueurs d'onde dans le domaine du visible. La réflectance déterminée est notée p dans la suite de la description. 5 La réflectance déterminée p par le procédé de détermination est avantageusement proche de la réflectance réelle réelle sur une gamme de longueurs d'onde comme décrit par la suite. Une erreur F de la réflectance déterminée p par rapport à la réflectance réelle réelle est une fonction de la longueur d'onde définie à partir de l'écart-type entre la réflectance 10 réelle réelle et la réflectance déterminée p. L'erreur F est une fonction nulle si la réflectance déterminée p est égale à la réflectance réelle réelle pour toutes les longueurs d'ondes. Plus l'erreur F est grande, plus cela indique que l'écart entre la réflectance déterminée p et la réflectance réelle réelle est grand. 15 Il est entendu que le procédé de détermination de la réflectance p et le dispositif 1 de détermination de la réflectance p d'un objet 4, sont applicables pour déterminer la réflectance p de tout objet 4. Par exemple, l'objet 4 est une partie de la peau d'un patient, un code-barres en couleur, une peinture, un produit cosmétique comme un fond de teint ou autre. 20 L'objet 4 est disposé dans un environnement comprenant un nombre inconnu de sources d'éclairement de l'objet 4. L'ensemble des sources éclairant l'objet 4 peut varier au cours du procédé de détermination de la réflectance de l'objet 4. En outre, l'éclairement, provenant des différentes sources éclairant l'objet 4 et distinctes de la source 10 du dispositif 1 de détermination, peut fluctuer au cours du procédé de détermination de la 25 réflectance de l'objet 4. Par exemple, un objet 4, placé dans une vitrine éclairée, est éclairé par la lumière du jour traversant la vitrine et par les lampes intérieures du magasin, avec des flux lumineux inconnus et variables. L'ensemble des sources éclairant l'objet 4 et distinctes de la source 10 du dispositif 30 1 de détermination est représenté par un illuminant extérieur 6 présentant un éclairement inconnu et variable en fonction du temps, noté I. L'éclairement I de l'illuminant extérieur 6 dépend de l'ensemble des fluctuations des sources d'éclairement de l'objet 4 distinctes de la source du dispositif 10 éclairant l'objet 4. Dans des conditions particulières, l'éclairement I de l'illuminant extérieur 6 est fixe. 35 A chaque instant, l'objet 4 est éclairé par l'illuminant extérieur 6 et, éventuellement, par la source 10. L'éclairement reçu par l'objet 4 est la somme de l'éclairement provenant 3041756 9 de la source 10 du dispositif 1 de détermination avec l'éclairement de l'illuminant extérieur 6. A chaque instant t, l'objet 4 éclairé par l'illuminant extérieur 6 et, éventuellement, par la source 10, réfléchit une onde 20 dépendant de la réflectance réellede l'objet 4. 5 , réelle Dans le mode de réalisation du dispositif 1 de détermination représenté sur la figure 1, le capteur 12 et la source 10 sont disposés sur un même appareil 16. Par exemple, l'appareil 16 est une tablette tactile, un téléphone portable, un ordiphone, ou autre. La source 10 est apte à émettre au moins un flash de lumière 18 éclairant l'objet 4. 10 Par exemple, la source 10 est un écran lumineux ou un ensemble de lampes. Par exemple, l'ensemble de lampes est un ensemble de diodes électroluminescentes (désigné par l'acronyme « LED »). Un flash de lumière 18 est un flux lumineux émis pendant un court intervalle de temps. Par exemple, l'intervalle de temps d'émission est compris entre 1 ms (milliseconde) 15 et 2 s (seconde). L'intervalle de temps d'émission dépend des caractéristiques de la source 10 et du capteur 12. Le flux lumineux du flash de lumière 18 présente une intensité d'émission en fonction de la longueur d'onde dépendant de la source 10. La source 10 est apte à émettre chaque flash 18 dans le domaine du visible. Cela 20 signifie que pour chaque flash 18, l'intensité d'émission est supérieure à un seuil de perception de l'oeil humain pour au moins une longueur d'ondes comprises entre 380 nm (nanomètres) et 800 nm. En variante ou en complément, la source 10 est apte à émettre dans le domaine de l'infrarouge, notamment à une longueur d'onde comprise entre 800 nm et 1000 nm. 25 Chaque flash i de lumière 18 émis par la source 10 présente un éclairement de l'objet 4, noté E, et dépendant de la longueur d'onde. Pour chaque flash i, l'éclairement E, dans une gamme de longueurs d'onde, est une caractéristique de la source 10, connue. Les caractéristiques de la source 10 sont conservées dans une mémoire. Les caractéristiques de la source 10 sont déterminées avant la mise en oeuvre du 30 procédé de détermination. La gamme de longueurs d'onde est délimitée par une valeur de longueur d'onde minimum Àmin et une valeur de longueur d'onde maximum À'x. Chaque longueur d'onde de la gamme de longueurs d'onde est comprise entre la valeur de longueur d'onde minimum min et la valeur de longueur d'onde maximum À'x. La gamme de longueurs d'onde dépend 35 de la source 10 et du capteur 12 utilisé pour l'émission et la réception des flashs 18. 3041756 10 La source 10 est propre à émettre des flashs de lumière 18. L'éclairement E, d'un flash i présente un maximum d'éclairement dans la gamme de longueurs d'onde. Le maximum d'éclairement est un maximum global de l'éclairement E, en fonction de la longueur d'onde. Le maximum d'éclairement d'un flash i se trouve à une longueur d'onde 5 notée À. L'instant où la source 10 émet un flash est appelé instant d'émission. Chaque instant d'émission d'un flash i est noté t,. L'instant d'émission t, de chaque flash i est une donnée conservée dans une mémoire. La source 10 est propre à émettre plusieurs flashs 18 successivement. L'intervalle 10 de temps entre deux instants d'émission t,, t, de flashs i, j successifs est, par exemple, compris entre 1 ms et 2 s. La source 10 est propre à émettre plusieurs flashs 18 de couleurs différentes, c'est-à-dire présentant des spectres différents. La source 10 est propre à émettre au moins deux flashs i, j de lumière 18 émis présentant un maximum d'éclairement à des longueurs 15 d'ondes À, et À éloignées d'au moins 20 nanomètres (nm). Dans un exemple, la source 10 est propre à émettre quatre flashs : un flash bleu, un flash rouge, un flash vert et un flash blanc. Sur la figure 3, chaque courbe 110, 112, 114, 116 représente l'éclairement de l'objet 4 par la source 10 en fonction de la longueur d'onde, sans illuminant extérieur 6, pour un flash respectif. Sur la figure 3, le flash bleu « flash1 », 20 correspondant à la première courbe 110 présente un maximum d'éclairement à la valeur Àflash 1 , alors que le flash rouge « flash2 », correspondant à la deuxième courbe notée 112 présente un maximum d'éclairement à la valeur À -flash2. Le capteur 12 est apte à collecter l'onde réfléchie 20 par l'objet 4 pour former au moins une image. 25 Par exemple, le capteur 12 est un appareil photo ou une caméra. Le capteur 12 est propre à détecter des intensités lumineuses dans la gamme de longueurs d'onde d'émission de la source 10. La sensibilité du capteur 12 dans la gamme de longueurs d'onde est une caractéristique du capteur 12 conservée dans une mémoire. Les caractéristiques du 30 capteur 12 sont déterminées avant la mise en oeuvre du procédé de détermination. En outre, la partie sensible du capteur 12 n'est pas orientée vers la source 10. Ceci permet que la collection de l'onde réfléchie 20 ne soit pas perturbée par la lumière directe du flash de lumière 18 émis par la source 10. L'image formée contient des données colorimétriques. Pour chaque image k 35 formée, les données de l'image sont notées Bk. 3041756 11 L'instant où le capteur 12 forme une image est appelé instant de collection. Chaque instant de collection d'une image k est noté tk. De même, l'instant de collection tk de chaque image k est une donnée conservée dans une mémoire. Le capteur 12 est propre à collecter plusieurs images successivement. Le capteur 5 12 est rapide, c'est-à-dire que le capteur 12 est propre à former des images à des instants de collection rapprochées. L'intervalle de temps entre deux instant de collection tk, t, d'images de flashs k, I successifs est, par exemple, compris entre 1 ms et 2 s. .Dans le mode de réalisation du dispositif 1 représenté sur la figure 1, l'unité de traitement 14 est disposée sur l'appareil 16 où sont disposés la source 10 et le capteur 12. 10 L'unité de traitement 14 comprend, par exemple, des processeurs et des mémoires. L'unité de traitement 14 est apte à traiter des données. L'unité de traitement 14 est, de plus, propre à recevoir les données du capteur 12 relatives à chaque image k formée et à chaque instant de collection tk et les données de la source 10 relatives aux éclairements E, de chaque flash i émis et aux instants d'émission t,. 15 L'unité de traitement 14 est propre à obtenir une équation à plusieurs inconnues à partir des images formées. Dans la suite, l'équation obtenue est notée équation à résoudre (1). La réflectance p de l'objet 4 et l'éclairement I de l'illuminant extérieur 6 sont deux inconnues de l'équation à résoudre (1). 20 En outre, l'unité de traitement 14 est propre à assurer la résolution de l'équation à résoudre (1). L'obtention et la résolution de l'équation à résoudre (1) par l'unité de traitement 14 sont décrites dans la suite de la description. Le fonctionnement du dispositif 1 est maintenant décrit en référence à la figure 2 qui 25 est un ordinogramme d'un premier exemple de mise en oeuvre du procédé de détermination de la réflectance p de l'objet 4. Le procédé de détermination de la réflectance p comprend les cinq étapes suivantes : une étape d'éclairement 100, une étape d'émission 102, une étape de collection 104, une étape d'obtention 106 et une étape de résolution 108. 30 Lors de l'étape d'éclairement 100, à chaque instant noté t, l'objet 4 est éclairé par l'illuminant extérieur 6 présentant l'éclairement 1(t). L'étape d'émission 102 est mise en oeuvre par la source 10 du dispositif 1. Lors de l'étape d'émission 102, la source 10 émet à un instant d'émission t,, un flash de lumière i éclairant l'objet 4. 3041756 12 Dans un mode de réalisation préféré, lors de l'étape d'émission, la source 10 émet une pluralité de flashs de lumière 18, chaque flash i étant émis à différents instants d'émission t,. Les instants d'émission t, et les éclairements E, de chaque flash i de lumière 18 émis 5 sont transmis à l'unité de traitement 14. L'étape de collection 104 est mise en oeuvre par le capteur 12 du dispositif 1. Lors de l'étape de collection 104, l'onde réfléchie 20 par l'objet 4 est collectée pour former au moins une image à un instant de collection sur le capteur 12. Dans un mode de réalisation préféré, l'étape de collection 104 est mise en oeuvre 10 pour chaque flash de lumière 18 émis. Les données relatives aux instants de collection et aux images formées sont transmises à l'unité de traitement 14. L'unité de traitement 14 reçoit les données relatives aux éclairements des flashs, aux instants d'émission t,, aux instants de collection tk et aux images formées 15 respectivement de la source 10 et du capteur 12. Lors de l'étape d'obtention 106, il est obtenu l'équation à résoudre (1). L'étape d'obtention 106 est mise en oeuvre par l'unité de traitement 14. L'unité de traitement 14 convertit chaque image formée en une équation. Par exemple, pour chaque image k, l'unité de traitement 14 détermine si l'instant de 20 collection tk de l'image k a lieu à un instant d'émission d'un flash. Deux cas sont alors possibles. Dans un premier cas, si l'instant de collection t, d'une image i est un instant d'émission d'un flash, le flash est noté i, et l'unité de traitement 14 obtient, à partir de l'image k, une première équation (2) reliant les données Bi de l'image i formée à la 25 réflectance p de l'objet 4. La première équation (2) s'écrit sous la forme suivante dans le cas d'une surface Lambertienne: Bi =K1 EL + 1(t1)) * b (.1)a où : - K1 est une première constante, 30 - `"' désigne l'opération de multiplication, Â-ma - fminx f (À) désigne l'opération mathématique d'intégrer la fonction f sur la variable À dans l'intervalle FA À 1 L min I max., - 3041756 13 Dans un deuxième cas, si l'instant de collection tk d'une image k n'est pas un instant d'émission d'un flash, l'unité de traitement 14 obtient, à partir de l'image k, une deuxième équation (3) reliant les données Bk de l'image k formée à la réflectance p de l'objet 4. La deuxième équation (3) s'écrit sous la forme suivante : 2-max Bk = K1 p *I(ti) * b (À)dÀ 2-min 5 L'unité de traitement 14 extrait alors au moins une première équation (2) ou une deuxième équation (3) pour former l'équation à résoudre (1). Avantageusement, l'équation à résoudre (1) est un système d'équations formé à partir de premières équations (2) et /ou de deuxièmes équations (3) obtenues pour plusieurs images formées. 10 L'équation (1) comprend N1 premières équations (2) et N2 deuxièmes équations (3), où N1 est un nombre entier naturel non nul d'images formées à un instant de collection qui est un instant d'émission et N2 est un nombre entier d'images formées à un instant de collection sans émission de flash. La réflectance p de l'objet 4 et l'éclairement I de l'illuminant extérieur 6 sont deux 15 inconnues de l'équation à résoudre (1). Pour la suite, il est à noter que l'équation à résoudre (1) est aisément représentable sous la forme d'une équation matricielle, pour une surface Lambertienne. A l'issue de l'étape d'obtention 106, l'équation à résoudre (1) a été obtenue à l'aide de l'unité de traitement 14. 20 L'étape de résolution 108 est alors mise en oeuvre par l'unité de traitement 14. L'étape de résolution 108 vise à résoudre l'équation à résoudre (1). Plus précisément, lors de l'étape de résolution 108, l'objectif est de rechercher une réflectance p, solution de l'équation à résoudre (1). Dans le cas général, la solution de l'équation à résoudre (1) est fortement sensible 25 aux erreurs d'observation et/ou de modélisation. Pour limiter l'influence des erreurs sur la solution p, il est proposé de combiner trois sous-étapes lors de la résolution 108. Aussi, selon l'exemple de mise en oeuvre de la figure 2, l'étape de résolution 108 comprend trois sous-étapes. 30 Les trois sous-étapes 150, 152, 154 de la résolution 108 sont mises en oeuvre successivement ou en parallèle par l'unité de traitement 14. Lors de la première sous-étape 150, il est calculé les points solution de l'équation à résoudre (1). 3041756 14 L'unité de traitement 14 détermine le nombre N de points solution P à calculer. Le nombre N de points solution à calculer est un nombre entier naturel non nul. Chaque point solution P à calculer comporte deux coordonnées, une abscisse et une ordonnée. 5 Selon l'exemple proposé, les abscisses sont d'abord déterminées puis les ordonnées associées sont calculées. L'unité de traitement 14 détermine N longueurs d'onde de calcul Àp, chaque longueur d'onde de calcul étant l'abscisse d'un point solution. Typiquement, par exemple, le nombre N de points solutions P à calculer est compris 10 entre 4 et 10. Par exemple, les longueurs d'ondes de calcul À sont réparties uniformément sur la gamme de longueurs d'onde. Pour chaque point solution P, l'unité de traitement 14 calcule une valeur de réflectance pp associée à la longueur d'onde de calcul À vérifiant l'équation à résoudre (1). 15 A l'issue de la sous-étape de calcul 150 de points solution, l'unité de traitement 14 obtient une pluralité de points solution P. Chaque point solution P comprend une valeur de réflectance pp associée à une longueur d'onde de calcul À vérifiant l'équation à résoudre (1). Lors de la deuxième sous-étape 152, une interpolation des points solution P est 20 mise en oeuvre par l'unité de traitement 14 à l'aide d'une fonction d'interpolation. Pour mettre en oeuvre la deuxième sous-étape 152, il est utilisé les points solution P trouvés lors de la première sous-étape, l'équation à résoudre (1) et des critères d'interpolation. Les critères d'interpolation définissent le type de fonctions d'interpolation à 25 examiner. Selon un exemple, les critères d'interpolation délimitent un espace autour de l'équation à résoudre (1) par lequel les points de la fonction d'interpolation doit passer. Selon un autre exemple, les critères d'interpolation limitent les fonctions d'interpolation à utiliser. 30 Ainsi, dans un cas particulier, la fonction d'interpolation est écrite sous la forme d'une combinaison pondérée d'un nombre fini np de fonctions de bases it.k Par exemple, la fonction d'interpolation de la réflectance p s'écrit sous la forme suivante : np p(À) = 1 ak *(k(À) k=1 3041756 15 où les coefficients ak sont les poids associés aux fonctions de base (pk. Lors de la deuxième sous-étape 152, l'unité de traitement 14 détermine les valeurs des poids ak et les formes des fonctions de bases 4)k. Par exemple, selon un critère d'interpolation, l'unité de traitement 14 définit chaque 5 point solution P comme un point d'interpolation. En variante, d'autres critères d'interpolation sont utilisés lors de la deuxième sous-étape 152. Par exemple, selon un critère d'interpolation, les fonctions de base (pksont des splines cubiques scellées. Une spline cubique est un polynôme cubique défini par 10 morceaux. Chaque morceau de la fonction est une fonction polynômiale d'ordre trois sur chaque intervalle de longueurs d'ondes délimité par deux points d'interpolation. A l'issue de la deuxième sous-étape 152 des points calculés, l'unité de traitement 14 obtient une fonction d'interpolation vérifiant les critères d'interpolation. La demanderesse a constaté que cette manière de résoudre l'équation à résoudre 15 (1) conduit parfois à des solutions non optimales. Pour pallier ce problème, une troisième sous-étape 154 est mise en oeuvre en même temps que la première sous-étape 150 ou la deuxième sous-étape 152. Lors de cette troisième sous-étape 154, il est utilisé au moins une approximation parmi une première approximation et une deuxième approximation. 20 Selon la première approximation, chaque image est issue de l'émission d'un flash de lumière 18 distinct. Selon la deuxième approximation, la fonction d'interpolation détermine les points de stabilité de l'équation à résoudre (1). A l'issue de l'étape de résolution 108, l'unité de traitement 14 obtient une réflectance 25 déterminée p. Le choix des approximations résulte de tests de la demanderesse d'une pluralité d'approximations possibles, ces approximations ayant l'avantage de rendre le procédé de détermination de la réflectance fiable et facile à mettre en oeuvre. Dans chacun des cas, il en résulte une meilleure réflectance déterminée p à temps 30 de calcul égal. Des modes d'implémentation particuliers de chacune des approximations sont détaillés dans ce qui suit. La figure 4 est un ordinogramme illustrant une mise en oeuvre particulière du procédé de détermination de la réflectance p lorsque la première approximation est mise en 35 oeuvre. 3041756 16 Les mêmes étapes que pour la mise en oeuvre du procédé selon la figure 2 sont mises en oeuvre. L'étape de collection 104 est mise en oeuvre plusieurs fois lors d'un même flash de lumière. 5 En utilisant la deuxième approximation, il est considéré que chaque image provient d'un flash différent. Il en résulte qu'à l'étape d'obtention 106, l'équation à résoudre (1) comporte plus d'équations que d'inconnues. L'équation à résoudre (1) est donc surdéterminée. Il est alors extrait de l'équation à résoudre (1) une pluralité de sous-équations à 10 résoudre. Chaque sous-équation à résoudre forme un système surdéterminé. A l'étape de résolution 108, il est résolu chaque sous-équation. Il en résulte une pluralité de réflectances solution ,solution, solution, chacune vérifiant l'équation à résoudre (1). Pour déterminer la réflectance p, il est procédé au calcul de la moyenne de la 15 pluralité de réflectances solution ,solution solution lors d'une quatrième sous-étape 180 de l'étape de résolution 108. Par exemple, le calcul de moyenne est mis en oeuvre par un calcul de moyenne arithmétique. Selon un autre exemple, le calcul de moyenne est mis en oeuvre par un calcul de 20 moyenne quadratique. Un tel mode de réalisation est aisé à mettre en oeuvre puisqu'aucun flash supplémentaire n'est impliqué. En particulier, tout en améliorant la précision, un tel procédé est mis en oeuvre avec la même rapidité. Une mise en oeuvre particulière du procédé de détermination de la réflectance p est 25 maintenant décrite lorsque la deuxième approximation est mise en oeuvre. Dans cet exemple, il est proposé de réduire le nombre de points d'interpolation aux points de stabilité de l'équation. Selon la deuxième approximation, les points de stabilité de l'équation sont déterminés par les fonctions d'interpolation. 30 De préférence, les points d'interpolation sont répartis sur un intervalle de longueurs d'onde, l'intervalle de longueur d'onde étant, par exemple, compris entre 380 nanomètres et 780 nanomètres dans le cas d'un ordiphone. Cette deuxième approximation résulte de travaux de la demanderesse illustrés par les figures 5 à 7. 35 Une simulation numérique de l'erreur F a été réalisée pour plusieurs simulations de détermination de la réflectance p. 3041756 17 Dans chaque simulation, des valeurs de bruits différents sont ajoutées à chaque donnée d'images et à chaque éclairement de flash. Par exemple, le bruit ajouté est un bruit blanc gaussien. Le bruit est une modélisation des différents défauts à la source 10 ou au capteur 12. 5 La figure 5 représente la fonction d'erreur F en fonction de la longueur d'onde selon le nombre de points d'interpolation retenus pour la fonction d'interpolation de p (À). Chaque courbe 200, 202, 204, 206, 208 est obtenue respectivement pour neuf, huit, sept, six ou cinq points d'interpolation. L'analyse de la figure 5 montre que certaines longueurs d'onde sont plus sensibles 10 au bruit que d'autres et que plus le nombre de points d'interpolation est important, c'est-à- dire que plus la fonction d'interpolation comprend de fonctions de base, plus l'équation à résoudre (1) est sensible aux instabilités. La figure 6 représente la fonction d'erreur F en fonction de la longueur d'onde pour cinq points d'interpolation. L'échelle a été modifiée par rapport à l'échelle de la figure 5 pour 15 montrer les détails de la fonction d'erreur F. La fonction d'erreur F présente 4 minima 220, qui sont des points de stabilité de l'équation à résoudre (1). La réflectance p trouvée pour les longueurs d'ondes minimisant la fonction F est moins sensible au bruit. Dans l'exemple, le nombre de point d'interpolation est donc égal à 4. Avantageusement, le nombre de points d'interpolation est égal au nombre de flashs 20 18. Sur la figure 7, chaque courbe 230 en trait pleins représente la réflectance p en fonction de la longueur d'onde pour cent réflectances p calculées avec différentes simulation de bruit. La courbe 234 en pointillés représente la réflectance réelle, connue par d'autres moyens. Les points de stabilités 220 sont les des points autour desquels l'écart 25 entre la réflectance p réel et les réflectances p calculés est moins important qu'à l'écart des points de stabilité 220. Ces simulations illustrent l'intérêt de la deuxième approximation. Dans tous les modes de réalisation, le procédé de détermination de la réflectance p d'un objet permet d'obtenir une modélisation fiable de la réflectance p réelle de l'objet 30 même avec un illuminant extérieur variable. Par fiable, il est entendu que la réflectance p déterminée par l'interpolation dépend peu des erreurs liés aux bruits de la source 10 et du capteur 12. Les figures 8 et 9 illustrent l'avantage de la mise en oeuvre du procédé de détermination de la demanderesse. Sur les figures 8 et 9, la courbe 240 en pointillés 35 représente la réflectance réelle réelle connue par d'autres moyens. 3041756 18 Sur la figure 8, la courbe en trait plein 242 représente la réflectance déterminée avec la mise en oeuvre d'un procédé n'utilisant ni la première approximation, ni la deuxième approximation. Sur la figure 9, la courbe en trait plein 244 représente la réflectance déterminée avec la mise en oeuvre d'un procédé de détermination selon l'invention. 5 L'analyse des figures 8 et 9 montre que la réflectance déterminée à partir du procédé de détermination est plus proche de la réflectance réelle préelle que la réflectance déterminée à partir du procédé n'utilisant ni la première approximation, ni la deuxième approximation. Dans une variante, le capteur 12 et la source 10 sont disposés sur des appareils 16 différents. 10 Dans une variante, le capteur 12 et la source 10 sont disposés sur le même appareil 16 et l'unité de traitement 14 est située à l'écart de l'appareil 16. En outre, il est à noter que le procédé permet de déterminer la réflectance de la surface observée pour chaque point image du capteur 12. 15 REVENDICATIONS1.- Procédé de détermination de la réflectance d'un objet (4), le procédé comprenant les étapes de : éclairement de l'objet (4) par un illuminant extérieur (6) présentant un éclairement inconnu et variable , émission d'au moins un flash de lumière (18) éclairant l'objet (4), chaque flash de lumière (18) étant émis par une source (10) et présentant un éclairement connu dans une gamme de longueurs d'onde, collection de l'onde réfléchie (20) par l'objet (4) pour former au moins une image sur un capteur (12), obtention d'une équation à plusieurs inconnues, l'équation étant obtenue à partir des images formées, la réflectance de l'objet (4) et l'éclairement de l'illuminant extérieur (6) étant deux inconnues de l'équation, résolution de l'équation, l'étape de résolution de l'équation comprenant le calcul de points solution de l'équation, l'interpolation des points calculés par une fonction d'interpolation, et l'utilisation d'au moins une des approximations suivantes pour la résolution de l'équation : o une première approximation selon laquelle chaque image est issue de l'émission d'un flash de lumière distinct, o une deuxième approximation selon laquelle la fonction d'interpolation détermine les points de stabilité de l'équation. 25 2.- Procédé selon la revendication 1, dans lequel la source (10) et le capteur (12) sont disposés sur un même appareil (16). 3.- Procédé selon la revendication 1 ou 2, dans lequel une pluralité de flashs de 30 lumière (18) sont émis, chaque flash (18) présentant un maximum d'éclairement en longueur d'onde, l'étape de collection étant mise en oeuvre pour chaque flash de lumière (18) émis et au moins deux flashs de lumière (18) présentent un maximum d'éclairement éloigné d'au moins 20 nanomètres. 35 4.- Procédé selon l'une quelconque des revendications 1 à 3, dans lequel l'étape de collection est mise en oeuvre plusieurs fois pour un même flash de lumière (18), l'équation 3041756 20 obtenue étant un système d'équations surdéterminé, l'étape de résolution étant mise en oeuvre pour une pluralité de systèmes d'équations déterminés en utilisant la première approximation pour obtenir une pluralité de fonctions de réflectance, le procédé comprenant, en outre, un calcul de la réflectance de l'objet (4) en calculant une moyenne 5 de la pluralité des fonctions de réflectance. 5.- Procédé selon l'une quelconque des revendications 1 à 4, dans lequel la deuxième approximation est utilisée lors de l'étape de résolution de l'équation et dans lequel la fonction d'interpolation est une combinaison pondérée de fonctions de bases 10 scellée par un nombre fini de points d'interpolation, notamment des splines cubiques, chaque point d'interpolation étant un point de stabilité de l'équation. 6.- Procédé selon la revendication 5, dans lequel une pluralité de flashs de lumière (18) sont émis, chaque flash (18) présentant un maximum d'éclairement en longueur 15 d'onde, l'étape de collection étant mise en oeuvre pour chaque flash de lumière (18) émis, et les points d'interpolations vérifient au moins la propriété suivante : - le nombre de point d'interpolation est égal au nombre de flashs (18). 7.- Dispositif de détermination (1) de la réflectance d'un objet (4), l'objet (4) étant 20 éclairé par un illuminant extérieur (6) présentant un éclairement inconnu et variable, le dispositif (1) comprenant : une source (10), apte à émettre au moins un flash de lumière (18) éclairant l'objet (4), chaque flash de lumière (18) émis par la source (10) présentant un éclairement connu dans une gamme de longueurs d'onde, un capteur (12), apte à collecter l'onde réfléchie (20) par l'objet (4) pour former au moins une image, une unité de traitement (14), propre mettre en oeuvre les étapes suivantes : o obtention d'une équation à plusieurs inconnues, l'équation étant obtenue à partir des images formées, la réflectance de l'objet (4) et l'éclairement de l'illuminant extérieur (6) étant deux inconnues de l'équation, o résolution de l'équation, l'étape de résolution de l'équation comprenant : le calcul de points solution de l'équation, l'interpolation des points calculés par une fonction d'interpolation, et l'utilisation d'au moins une des approximations suivantes pour la résolution de l'équation : 3041756 21 o une première approximation selon laquelle chaque image est issue de l'émission d'un flash de lumière (18) distinct, o une deuxième approximation selon laquelle la fonction d'interpolation détermine les points de stabilité de l'équation. 5 8.- Dispositif selon la revendication 7, dans lequel le capteur (12) et la source (10) sont disposés sur un même appareil (16). 9.- Dispositif selon la revendication 7 ou 8, dans lequel la source (10) est un écran 10 lumineux ou un ensemble de diodes électroluminescentes. 10.- Dispositif selon l'une quelconque des revendications 7 à 9, dans lequel le capteur (12) est choisi dans un groupe constitué d'un appareil photo, d'une caméra, d'un imageur multicanal et d'un imageur hyperspectral.