FR 2990822 A1 20131122 FR 1254564 A 20120518 PROCEDE D'ETALONNAGE D'UN ECRAN EMISSIF ET D'UN CAPTEUR ELECTRONIQUE D'IMAGES EN COULEUR Domaine de l'invention La présente invention se rapporte au domaine de la mesure de la couleur. La présente invention se rapporte plus particulièrement à un procédé pour réaliser l'étalonnage colorimétrique d'un écran émissif et l'étalonnage colorimétrique et télémétrique d'un capteur électronique d'image couleur et l'étalonnage de la sonde de luminance extérieure (s'il en existe une), tous situés côte à côte, à l'aide d'une mire de couleur annexe (e.g. type 24 couleurs GretagMcBeth - Marque déposée /X-RiteColorCheker - Marque déposée) et d'une feuille de papier quadrillé annexe (e.g. 5x5 mm), et au dispositif associé. Etat de la technique Une couleur est un flux lumineux perçu par l'oeil humain. Il s'agit d'une distribution spectrale d'énergie sur une bande de longueur d'onde, en l'occurrence 380 à 780 nm pour le visible. Les spectrophotomètres numériques professionnels mesurent actuellement ces distributions avec un pas de 5 nm, soit 81 valeurs pour le visible. Dans la suite, nous désignerons par une lettre et un À entre parenthèses les quatre-vingt-une valeurs qui décrivent une distribution spectrale énergétique S(À) dans le visible. Cette courbe S(À) peut être interpolée avec un nombre réduit de points [Ài,S(À);], en l'occurrence dans la suite de ce document une fonction spline cubique scellée avec des pentes nulles à 380 et 780 nm. Du fait de la synthèse additive, toute couleur peut être exprimée par ses coordonnées trichronnatiques selon trois couleurs primaires. Il existe plusieurs systèmes de primaires, nous retiendrons dans la suite les deux normalisés que sont l'espace CIE XYZ (et sa variante CIE xyY avec plan de chromacité à luminance constante) et l'espace CIE L*a*b qui permet d'estimer correctement une différence de couleur avec une norme euclidienne appelée « deltaE ». Le « gamut » (limites) de ces deux espaces couvre l'intégralité du visible humain. Nous feront aussi référence à l'espace trichromatique sRGB correspondant au gamut de référence de la plupart des équipements électroniques actuels. Le gamut du sRGB ne couvre pas l'ensemble du visible humain, notamment dans les bleus-verts. La couleur d'un objet est le spectre C(À) égal au produit d'un illuminant 1(À) se réfléchissant sur une surface dont la réflectance énergétique est R(À). Autrement dit, la couleur d'un objet peut être caractérisée par ses coordonnées trichromatiques CIE XYZ ou L*a*b pour un illuminant donné (exemple : D65 ou D50) ou, mieux, par sa distribution spectrale de réflectance énergétique R(À). Cette deuxième méthode présente l'intérêt de pouvoir simuler la perception de la couleur sous différents illuminants (intérieur, extérieur). C'est la méthode des nuanciers de couleurs : l'utilisateur pose le nuancier à côté de l'objet et recherche par itération la couleur du nuancier qui se rapproche le plus de la couleur de l'objet (ex : Pantone Matching System - Marque déposée). Un écran émissif produit des couleurs au niveau de chaque pixel par synthèse additive de 3 primaires rouge, vert et bleue. Plus précisément, il existe 4 paramètres d'entrées : le niveau de luminance général de l'écran CECR RE (E[0 ;1]), valable pour tous les pixels (proportionnel à la luminance du rétro-éclairage dans le cas particulier d'un écran à cristaux liquides), un niveau de rouge CECR R (E[0 ;1]), un niveau de vert CEcR V (e[0 ;1]) et un niveau de bleu cECR B MO ;1]). Un capteur électronique d'image couleur est un colorimètre tri-stimulus pour chaque pixel, en ce sens qu'il contient 3 filtres de couleur (rouge, vert, bleu) et 3 photodétecteurs derrière qui convertissent la luminance en tension électrique. La plupart des capteurs d'image possèdent également en amont une lentille auto-focus fonctionnant sur le principe de la maximisation du contraste. Le paramètre de pilotage de la lentille est lié à sa dioptrie d. Les quatre paramètres de sortie du capteur électronique d'image pour chaque pixel sont cCEI R (e[o D, cCEI V (E[0 ;1]), cCEI B (E[o ]) et d, (dioptrie de la lentille autofocus - nombre réel). Actuellement, l'étalonnage colorimétrique d'un écran émissif nécessite une sonde électronique annexe d'un coût très supérieur à celui d'une mire de couleurs. Actuellement, l'étalonnage colorimétrique et télémétrique d'un capteur électronique d'image nécessite une sonde électronique annexe d'un coût très supérieur à celui d'une mire de couleurs. L'usage actuel d'une mire de couleur avec un capteur d'image (exemple : photographie) n'a pas pour fonction l'étalonnage colorimétrique du capteur d'image comme expliqué plus haut, simplement de réaliser l'adaptation chromatique des couleurs des objets photographiés, sous un illuminant donné (cf. balance des blancs). Actuellement, l'étalonnage de la sonde de luminance extérieure nécessite une sonde électronique annexe d'un coût très supérieur à celui d'une mire de couleurs. Le procédé et le dispositif de ce document découlent de la demande de brevet français N° FR 12 51484 (Franck HENNEBELLE) proposant un procédé de mesure de la couleur d'un objet par spectrométrie utilisant comme émetteur un écran émissif et comme récepteur un capteur électronique d'image situés côte à côte. L'art antérieur connaît pour ce dernier, par la demande de brevet PCT N° VVO 2004079314, un colorimètre, un boîtier de détection à colorimètre et un procédé de détermination de la couleur. Le fonctionnement de la mesure de couleur avec ce procédé de l'art antérieur est réalisé avec un boîtier développé spécialement et sans aucune mention du procédé d'étalonnage.30 Exposé de l'invention La présente invention concerne, dans son acception la plus générale, un procédé d'étalonnage d'un dispositif comportant un écran émissif et un capteur électronique d'images en couleur (tri stimulus), ledit capteur électronique d'images comportant une lentille autofocus, caractérisé en ce qu'il comporte les étapes suivantes : a) mesure de la surexposition dudit capteur électronique d'images relativement à la puissance lumineuse de l'écran ; b) mesure du seuil de détection dudit capteur électronique d'images ; c) calibration des fonctions gammas RN/B dudit capteur électronique d'images ; d) calibration de la fonction gamma du rétroéclairage de l'écran ; e) calibration des fonctions gammas RN/B de l'écran ; f) calcul d'une matrice de transfert colorinnétrique dudit capteur électronique d'images des CCEI RN/B linéairee [0;1] vers des coordonnées CIE XYZ ; g) calcul des spectres des trois primaires de l'écran rouge E(À)R MAX, vert E(À)v MAX et bleu E(À)B MAX; h) calcul d'ajustement entre la matrice de transfert colorimétrique de dudit capteur électronique d'images [mcEl] et les spectres des trois primaires de l'écran rouge E(À)R MAX, vert Eve MAX et bleu E(À)B MAX ; i) détermination de deux fonctions de transfert permettant de passer de la dioptrie d de ladite lentille autofocus à la mesure de distances dans le plan focal image, et à la distance entre le foyer de ladite lentille autofocus et le plan focal image ; et j) détermination de l'écartement entre le foyer de la lentille et le bord de l'écran. Le procédé selon la présente invention permet d'effectuer un étalonnage avec une mire de couleur et du papier quadrillé à un coût très réduit. De préférence, ledit écran et ledit capteur électronique d'images sont situés côte à côte. Avantageusement, ledit procédé comporte en outre une étape h') d'itération des étapes c) à h). De préférence, ledit procédé utilise une mire de couleur et du papier quadrillé. - Selon un mode de réalisation, les fonctions ,gammas de l'étape c) convertissent des CRIVIB linéairesE [0;1] captés en valeurs CRN/BE [0;1] afin de correspondre à la sensibilité de l'oeil humain. Selon un mode de réalisation, la fonction gamma de l'étape d) convertit le paramètre de pilotage du rétroéclairage CECR REE [0 ;1] en valeur linéaire cECR RE linéaire e [0 ;1] correspondant à la luminance effectivement produite par l'écran. Avantageusement, les fonctions gammas RN/B de l'écran de l'étape e) convertissent les CRN/BE[0 ;1] en CRA1113 linéairesE[0 ;1] afin de correspondre à la sensibilité de l'ceil humain. Selon un mode de réalisation, ladite matrice de transfert est fonction de filtres de couleurs placés devant des photodétecteurs. Avantageusement, ledit procédé comporte en outre une étape e') de calcul de la borne inférieure de saturation du capteur électronique d'image et de sa borne supérieure de détection sous la forme de valeurs de luminance relative à celle de l'écran, et d'étalonnage d'une sonde de luminance extérieure (si elle existe). Selon un mode de réalisation, ledit procédé est mis en oeuvre pour paramétrer un écran. Selon un mode de réalisation, ledit procédé est mis en oeuvre pour mesurer des distances au moyen du capteur électronique d'images. Selon un mode de réalisation, ledit procédé est mis en oeuvre pour corriger l'effet « fish eye » de la lentille. Selon un mode de réalisation, ledit procédé est mis en oeuvre pour mesurer des inclinaisons au moyen du capteur électronique d'images. Selon un mode de réalisation, ledit procédé est mis en oeuvre pour effectuer des mesures de couleur. Selon un mode de réalisation, ledit dispositif est un téléphone mobile. Selon un mode de réalisation, ledit dispositif est une tablette numérique. Selon un mode de réalisation, ledit dispositif est un écran d'ordinateur ou une télévision équipé(e) d'une caméra de visioconférence. Avantageusement, ledit capteur électronique d'images comporte un capteur de luminance. 30 Brève description des dessins On comprendra mieux l'invention à l'aide de la description, faite ci-après à titre purement explicatif, d'un mode de réalisation de l'invention, en référence aux Figures dans lesquelles : 20 25 - la Figure 1 illustre de façon schématique le dispositif mis en oeuvre par le procédé selon la présente invention ; - la Figure 2 représente les différentes étapes du procédé selon la présente invention ; - la Figure 3 illustre l'algorithme de génération de N flashs colorés répartis uniformément dans le gamut de l'écran (exemple pour un pas de 1/(21'p) = 1/4 avec p=2) ; - la Figure 4 représente un pic de gris inversé, qui permet de localiser une ligne de quadrillage ; - la Figure 5 illustre la détermination du point zéro provisoire, ainsi que la direction de la ligne, à travers les vecteurs U et V; et - la Figure 6 représente la détection de la première case du quadrillage. Description détaillée des modes de réalisation de l'invention La Figure 1 représente sur une coupe schématique le dispositif 10 mis en oeuvre par le procédé selon l'invention : un objet 30; l'écran émissif 11 ; la lumière environnante extérieure 40 ; le capteur d'image 12; l'application informatique 13 « colorimétriquement pré-étalonnée » qui pilote l'ensemble et interagit avec l'opérateur. Dans un mode de réalisation, l'application informatique 13 est partiellement ou totalement distribuée dans le réseau ou « nuage » (« cloud » en terminologie anglo-saxonne). La Figure 2 représente les différentes étapes du procédé selon la présente invention : a) mesure de la surexposition dudit capteur 12 électronique d'images relativement à la puissance lumineuse de l'écran 11 ; b) mesure du seuil de détection dudit capteur 12 électronique d'images ; c) calibration des fonctions gammas RN/B dudit capteur 12 électronique d'images ; d) calibration de la fonction gamma du rétroéclairage de l'écran 11 ; e) calibration des fonctions gammas RN/B de l'écran 11; f) calcul d'une matrice de transfert colorimétrique dudit capteur 12 électronique d'images des CCEI RAHE linéairee [0;1] vers des coordonnées CIE XYZ ; g) calcul des spectres des trois primaires de l'écran 11 rouge E(À)R mAx, vert E(À)v MAX et bleu E(À)B MAX ; h) calcul d'ajustement entre la matrice de transfert colorimétrique dudit capteur 12 électronique d'images j et les spectres des trois primaires de l'écran 11 rouge E(À)R MAX vert E(j)" MAX , et bleu E(À)B MAX ; i) détermination de deux fonctions de transfert permettant de passer de la dioptrie d de ladite lentille autofocus à la mesure de distances dans le plan focal image, et à la distance entre le foyer de ladite lentille autofocus et le plan focal image ; et D détermination de l'écartement entre le foyer de la lentille et le bord de l'écran 11. Dans des modes de réalisations particuliers, le procédé selon la présente invention comporte en outre les étapes suivantes : o e') calcul de la borne inférieure de saturation du capteur électronique d'image et de sa borne supérieure de détection sous la forme de valeurs de luminance relative à celle de l'écran (11), et étalonnage d'une sonde de luminance extérieure (si elle existe) ; et - h') itération des étapes c) à h). Dans le mode de réalisation présenté Figure 1, l'objet à mesurer 30 est placé sous l'écran 11, près du bord de l'écran 11 et dans le champ de vision du capteur d'image 12. A titre de remarque préliminaire, nous considérons que la description détaillée ci-après est basée sur l'utilisation d'une mire de couleur annexe de type 24 couleurs GretagMcBeth - Marque déposée / X-Rite ColorCheker - Marque déposée dont les cases sont numérotées de 01 (marron peau sombre) à #24 (noir) et dont le procédé connaît la distribution spectrale de réflectance énergétique Rcce(À) et l'écart-type %ETRcc#i(À) (exprimé en 13/0 de Rcc #i(À)). Le fonctionnement du dispositif n'est pas restreint à l'utilisation de la mire ColorChecker, il est même opportun de développer une mire alternative pour ne pas être sous le joug commercial de X-Rite (par exemple, teintes différentes-et nombre de cases différent, tout en ayant une échelle de gris pour l'étape c) et une bonne couverture du gannuts RGB pour les étapes d) à h)). A titre de remarque préliminaire, nous considérons ici que les notations sont : CEI = capteur électronique d'image couleur (la fonction vidéo n'est pas utilisée dans le cadre de l'étalonnage, seule est utilisée la prise d'images statiques), ECR=écran ; R/V/B = rouge ou vert ou bleu (équation valable pour chaque primaire) ; RE=rétroéclairage écran ; CC #i = case ColorChecker numéro i ; sensibilités de l'observateur standard CIE 1931 XYZ 2° xstdobs(À), ystdobs(À), zStdObs (À) pour ÀE[380nm ;780nm]. A titre de remarque préliminaire, nous considérons que la description détaillée ci-après est basée sur un écran émissif à cristaux liquides. Pour chaque pixel, la luminance de l'écran résulte de l'action combinée du rétroéclairage et des cristaux liquides RVB : yECR = yECR RE yECR RVB= yECR RE .(yECR R ±yECR V +yECR B). A titre de remarque préliminaire, nous considérons que l'ensemble des pixels de l'écran sont considérés comme homogènes dans leur fonctionnement. Idem pour les pixels de la caméra. A titre de remarque préliminaire, nous considérons qu'avant le lancement du cycle d'étalonnage, le procédé connaît les valeurs théoriques de sortie « usine » de : - la luminance maximale blanche de l'écran (yECR BL MAXen cd/m2); - le spectre lumineux blanc associé E(À)BL MAX obtenu pour RVB=(1,1,1) et cECR RE=1 ; - les valeurs pECR RN/B=yECR RN/BNECR BL (%). De plus, toutes les fonctions gamma ci-après sont exprimées sous la forme du formalisme défini par la norme de l'International ColorCommittee (ICC) version 4, à savoir : Y = (a.X+b)Agamma+c. En outre, comme indiqué- à l'étape h'), le mode de calcul de l'étalonnage est itératif pour les étapes c) à h) (mais une seule prise de mesures pour l'utilisateur). La première itération est réalisée avec les valeurs théoriques « usine ». De plus, l'ordre exact des opérations ci-dessous peut être modifié de manière à optimiser le temps passé par l'utilisateur, par exemple : - étapes a) à b) : l'utilisateur pointe la case ColorChecker #19 (blanche) ; - étapes c) à h') : l'utilisateur réalise des flashs de mesures pour les 24 cases (par exemple, 1 flash blanc suivis par 45 à 153 flashs de couleurs), puis le procédé réalise les calculs d'étalonnage en temps masqué ; - étape i) : l'utilisateur présente une feuille de papier quadrillé centrée sur l'axe de la caméra à différentes distances de la caméra pour chaque mesure, puis le procédé réalise les calculs d'étalonnage en temps masqué ; - étape j) : idem à i) mais dans une salle noire, à défaut dans un environnement très faiblement lumineux et stable. Enfin, l'utilisateur doit réaliser l'étalonnage dans une salle noire, à défaut dans un environnement très faiblement lumineux et stable. L'écran émissif fournit à l'utilisateur un éclairage d'appoint suffisant pour réaliser ses manipulations. Etape a) : calibration de la borne inférieure de saturation (luminance) du capteur d'image Le procédé allume l'écran à sa luminance maximale : CECR RE=1 et (RVB)EcR=(1,1,1). L'utilisateur est invité à pointer la case ColorChecker#19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses (RVB)cEl Si (RVB)cE1=(1,1,1) ce qui constitue la situation la plus courante, le procédé diminue le rétroéclairage (par dichotomie) jusqu'à atteindre le seuil où (RVB)cEl passe en dessous de (1,1,1). Le procédé retient le pas de rétroéclairage juste inférieur. Les valeurs trouvées constituent temporairement la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » { case ColorChecker#19 ; CECR RE =X1 ; (RVB)EcR =(1,1,1) }. Le procédé passe ensuite au pas juste supérieur de rétroéclairage et il diminue (RVB)EcR pour revenir au seuil où (RVB)cEl passe en dessous de (1,1,1). Les valeurs trouvées constituent temporairement la borne inférieure de saturation de la caméra dite «2 » ou « RVB » {case ColorChecker#19; CE CR RE=X2 ; (RVB)EcR .(x2,x2,x2)}. Si (RVB)cEl Nota : par la suite, le procédé sera capable d'exprimer la borne inférieure de saturation de la caméra sous la forme de valeurs de luminance relative à celle de l'écranE(À)BL MAX. Etape b) : calibration de la borne supérieure de détection (luminance) du capteur d'image Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; cECR RE=A', I ; (RVB)ECR=(1 ) L'utilisateur est invité à pointer la case ColorChecker#19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses (RVB)cEl Comme (RVB)cEI>(0,0,0), le procédé baisse (par dichotomie)la valeur du rétroéclairage pour parvenir au seuil où (RVB)cE1=(0,0,0).Le procédé retient le pas juste supérieur. Les valeurs trouvées constituent temporairement la borne supérieure de détection de la caméra dite « 1 »/« rétroéclairage » casecc=19 ; cECR RE=X3 ; (RVB)ECR=(1,1,1) }. Le procédé allume de nouveau l'écran, avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19; CECR RE=A^^-1 ; (RVB)ECR=(1 ,1 ) L'utilisateur est toujours invité à pointer la case ColorChecker#19 (blanche) dans une salle noire. Comme (RVB)cEl>(0,0,0), le procédé baisse (par dichotomie) les valeurs (RVB)EcR pour parvenir au seuil où (RVB)EcR=(0,0,0). Le procédé retient le pas juste supérieur. Les valeurs trouvées constituent temporairement la borne supérieure de détection de la caméra dite «2 »/« RVB » { caseColorChecker#19 ; CEcR RE=X1 ; (RVB)EcR=(x3,x3,x3) Nota : par la suite, le procédé sera capable d'exprimer la borne supérieure de détection de la caméra sous la forme de valeurs de luminance relative à celle de l'écranE(À)BL MAX. Etape c) : calibra tion des fonctions gammas R/V/B de la caméra (elles convertissent les CRI" linéairee [0;1] captés en valeurs CRAllBE (0,11 afin de coller à la sensibilité de l'oeil humain) Pour chaque primaire RN/B, le procédé va déterminer les paramètres suivants de la fonction gamma (conformément à la norme ICC v4): cCEI RN/B = (a10EI R/V/B-yCEI RN/B+bCEI R/V/B)AgammaCEI RN/B+cCEI R/V/B ; cCEI RN/Blinéaire=a2CEI R/V/B-yCEI R/V/B ; _ cCEI RN/B = (aCEI R/V/B-cCEI RN/Blinéaire+bCEI R/V/13)AgammaCEI R/V/B+cCEI RN/B Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; cECR RE=X1 ; (RVB)ECR = (1,1,1) 1. La luminance de l'écran est alors égale à YiEcR, que le procédé calcule par la connaissance de E(À)BL MAX et 0/0yECR RE (Xi) yi ECR= cyoyECR RE (Xi) v^) r7EID elM(E(À)BL MAX yStdObs.... dÀ ) - J380 ren cyoyECR RE (Xi) = ((al ECR RE-xi ±blECR RE)AgammaECR RE+cl ECR RE)NECR BL MAX; Nota : pour la première itération, cyoyECR RE (X1) = X1 E[0;1] ; Ensuite l'utilisateur est invité à pointer successivement les casesColorChecker #24 (noire) à #19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses (RVB) caméra. Comme le procédé connaît E(À)BL MAX et les Rcc #i (À), il calcule le taux de perte de luminance après réflexion Rcc #i, lié la case ColorChecker #i, ainsi que le poids de chaque mesure Wi associé à l'écart-type ETRoe : RCC #i .-790 mn (Ego BL MAX RCC #i (À) yStdObs.. ^- dÀ ) /I 7E'D E(À)BL MAX HG mn 380 ?lm yStdObs(À),dÀ ) ; - ETR CC #i =f7 BL MAX #i (E(À) (À) . %ETRcc #i (À). yStdObsv,,\ 7°0 > mi) dÀ ) /j" 380 nin 380 nm E(À)BL MAX yStdObs(À).dÀ ) ; - Wi = 1/ ((ETRcc #5^2) ; 25 30 En posant a2CEI RN/B =pECR RN/BNi ECR, le procédé obtient de la sorte 6 points de mesure (Xi, Yi) de poids Wi pour chaque primaire R/V/B où : = cCEI R/V/B linéaire i = YCEI R/V/B ECR = RCC #i; _ Yi = cCEI RN/Bi = (aCEI RN/B.cCEI RN/B linéaire i +pCEI RN/B)AgammaCEI RN/B cCEI R/V/B La solution de régression correspond au minimum du « Chi 2» : x2 CEI R/V/B=57 wi.(yi_((aCEI R/V/B;xi+pCEI RN/B)AgammaCEI RN/B+pCEI RN/B))A2 Le procédé détermine alors par l'algorithme des moindres carrés les inconnues gammaCEI RN/B;aCEI RN/B; pCEI RN/B et pCEI R/V/B; puis ai CEI RN/B=aCEI RN/B a2CEI RN/Bpour chaque primaire RVB. Un point important pour garantir la convergence est de procéder en deux passes : ' 1) une régression linéaire (sur base logarithme) avec b''''=ccEi=0 2) la résolution d'un système non linéaire surdéterminé avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton) pour minimiser davantage le « Chi 2» x2 CEI RN/Bavec bcEl et cCEI . Etape d) calibration de la fonction gamma du rétroéclairage de l'écran (elle convertit le paramètre de pilotage du rétroéclairage CECR RE E [0 ;1] en valeur linéaire CECR RE linéaire e [0 ; 1 ] correspondant à la luminance effectivement produite par l'écran) Le procédé va déterminer les paramètres suivants de la fonction gamma (conformément à la norme ICC v4): yECR RE = (ai ECR RE. cECR RE+piECR RE)Ag a m maECR ECR RE; - CECR RE linéaire=a2ECR RE.yECR RE; - CECR RE linéaire = (aECR RE.cECRRE+pECR RE ) A gammaECRRE+CECR RE Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; CECR RE=X1 et (RVB)écran =(1,1,1). La luminance de l'écran est alors égale à Y1 ECR, que le procédé calcule par la connaissance de E(À)BL MAX et oAyECR RE (Xi) ECR= cyoyECR RE (XI) CSO E(À BL MAX. yStdObs,. v\) dÀ ) ; jSE10 Mn, oAyECR RE (xl ) = ((al ECR RE.xi +bi ECRA ) gammaECR RE+cI ECR RE)/yECR BL MAX. Nota : pour la première itération, cyoyECR RE (X1) = X1 c[0;1] ; L'utilisateur est invité à pointer la case ColorChecker#19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses (RVB)cEl. Comme le procédé connaît E(À)BI- MAX et Rcc #19 (À), il calcule le taux de perte de luminance après réflexion Rcc #19 lié à la case ColorChecker #19: 'RO 72-7 E/À)BL MAX _ R tcc2n9 = fuo (Ewn MAX RCC #19,, dObs,, v^) dÀ ) /J" n v^) .y J330 nm '380 n. '^ StdOb Y s(À).dÀ ) ; 15 Ensuite, le procédé diminue le rétroéclairage jusqu'à la borne supérieure de détection de la caméra dite « 1 » ou « rétroéclairage » { caseCC=19 ; cECR RE=X3 et (RVB)écran =(1,1,1) }. La diminution du rétroéclairage s'opère en si possible N pas (N=10 à 16, à défaut 5, voire 2 au 20 minimum). En posant a2ECR RE=1/Y1ECR, et comme l'étape précédente a permis de (al CEI RN/BNCEI RN/13±bCEI R/V/B)AgammaCEI RN/B paramétrer CCEI RN/B = +êEl RN/B , le procédé obtient de la sorte N points de mesure (Xi, Yi) de poids 25 Wi=1 pour chaque primaire RVB où : _ = cECR RE i ; (oECR REI = entrée du dispositif) Yi = CECRRElinéairei = a2ECR RE yECR RE i = yECR RE i Y1 ECR avec yECR REi = yCEI RN/Bi iRCC #19 = ((cCEI RN/Bi RN/B )A /(gam maCEI RN/B)..bCEI RN/B)/(a CEI RN/B 30 R cc#19) ; (cCEI RN/B. I = sortie du dispositif) 10 La solution de régression correspond au minimum du « Chi 2» : x2 ECR RE _ Yi _((aECRRE.xi+bECR RE)Agam maECRRE±cECR RE))A2 Le procédé détermine alors par l'algorithme des moindres carrés les inconnues gammaECR RE aECR RE; bECR RE et cECR RE; puis ai ECR RE=aECR RE.a2ECR RE Un point important pour garantir la convergence est de procéder en deux passes : 1) une régression linéaire (sur base logarithme) avec bECR RE=cECR RE=0 ; 2) la résolution d'un système non linéaire surdéterminé avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton) pour minimiser davantage le « Chi 2 » X2 ECR REavec b ECR RE et cECR RE ^e. Etape e): calibration des fonctions gammas R/V/B de l'écran (elles convertissent les Cfmse ;1] en CR/v/BlinéaireE[O ;1] afin de coller à la sensibilité de l'oeil humain) Pour chaque primaire RN/B, le procédé va déterminer les paramètres suivants 20 de la fonction gamma (conformément à la norme ICC v4): yECR RN/B = (a1ECR RN/BcECR RN/B+biECR RN/Brg am maECR R/V/B+ci ECR R/V/B . cECR RN/B linéaire = a2ECR RN/B.yECR RAI/B . cECR RN/B linéaire = (aECR RN/B.cECR R/V/B+bECR R/V/B,) Agamma ECR RN/B+CECR R/V/B 25 Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; rétroéclairage =X1 et (RVB)écran =(1,1,1). La luminance de l'écran est alors égale à Y1 ECR, que le procédé calcule par la connaissance de 30 E(À)BL MAX et 0/0yECR RE po ) y1 ECR= 0/0yECR RE (Xi) 1-7SQ E(À)BL MAX yStdObs(À).dÀ ) ; 480 nna 0/0yECR RE (Xi) = ((ai ECR RE.xi +b1 ECR RE)A gammaECR RE ECR RE)NECR BL MAX . +C I Nota : pour la première itération, (voyECR RE (X1) = X1 E[0;1] ; L'utilisateur est invité à pointer la case ColorChecker#19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses (RVB)cEl. Comme le procédé connaît E(À)BL MAX et Rcc #19 (À), il calcule le taux de perte de luminance après réflexion RCC *19 lié à la case ColorChecker#19: _ Rcc #19 =7O nni (E(À)BL MAX ; RCC #19(À) ;yStdObsi, ; v\) dÀ CD) mn r E(À)BL MAX 3O :17n. MG 1,1771 ystdobs(À).dA ; Ensuite, le procédé diminue le rétroéclairage jusqu'à la borne supérieure de détection de la caméra dite « 2 »/« RVB » { case ColorChecker#19 ; cECR RE=X-I ; (RVB)écran =(x3,x3,x3) }. La diminution du rétroéclairage s'opère en si possible N pas (N=10 à 16, à défaut 5, voire 2 au minimum). 15 En posant a2ECR P,/V/B = pECR RN/B/yi ECR; et comme les étapes précédentes ont (ai cEiNcEl+bcE5A g ammaCEI+cCEI ; permis de paramétrer CcEl = le procédé obtient de la sorte N points de mesure (Xi, Yi) de poids Wi=1 pour chaque primaire R/V/B où : 20 _ cECR RN/B ; (oECR R/V/B. = = entrée du dispositif) - Yi = CEGR RIVIB linéaire i = a2ECR R/V/B yECR RN/B = yECR R/V/B ECR RN/B/yi ECR avec yECR RN/Bi = yCEI RN/Bi / RCC#19 . = ((Cc EI RN/Bi _cCEI RN/6 /(gammaCEI RN/B)4DCEI RN/B)/(ai CEI R/V/B R cce19) ; (CcEt RAI/Bi = sortie du dispositif) 25 Pour chaque primaire RN/B, la solution de régression correspond au minimum du « Chi 2 » : x2 ECR R/V/B = iy; it,ECR R/V/B v;+bECR R/V/BrgammaECR RN/B+CECR R/V/B., A harz--1 "-^^°R 'Al )) 2 30 Le procédé détermine par l'algorithme des moindres carrés les inconnues gammaECR RN/B; aECR RN/B, bECR R/V/B et cECR RN/B; puis al ECR RN/B=aECR R/V/B CR R/V/B a2E pour chaque primaire R/V/B. 10 Un point important pour garantir la convergence est de procéder en deux passes : R RN/= 1) une régression linéaire (sur base logarithme) avec bEC BcECR RN/B=0 ; 2) la résolution d'un système non linéaire surdéterminé avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton) pour minimiser davantage le « chi 2 » )(,2 ECR RN/B avec bECR RN/B et cECR RN/B Etape e') : calcul de la borne inférieure de saturation de la caméra et de sa borne supérieure de détection sous la forme de valeurs de luminance relative à celle de l'écran et, étalonnage de la sonde de luminance extérieure (s'il en existe une) qui consiste à déterminer la fonction gamma de la sonde qui convertit le paramètre de sortie CsLEE [0 ;1] en valeur de luminance effectivement reçue par la sonde (cd/m2) Pour chaque primaire, le procédé utilise les paramètres déterminés aux étapes d) et e) afin de calculer la borne inférieure de saturation de la caméra et de sa borne supérieure de détection sous la forme de valeurs de luminance relative à celle de l'écran(cd/m2) à l'aide des formules suivantes : yECR = yECR RE (yECR R yECR V + yECR B) yECR RE = (al ECR RE. cECR RE+b ECR RE\ A- ) yammaECR ECR RE yECR RN/B = (al ECR RN/B.cECR R/V/B+L..,1.1 ECR RN/B )AgammaECR RIV/B+C1 ECR RN/B Ensuite, le procédé étalonne la sonde de luminance extérieure (s'il en existe une). Cette opération consiste à déterminer la fonction gamma de la sonde qui convertit le paramètre de sortie CE [0 ;1] en valeur de luminance effectivement reçue par la sonde (cd/m2). La connaissance de cette fonction servira à éviter les sous-exposition/surexposition de la caméra lors des mesures en présence d'un illuminant extérieur. Le procédé va déterminer les paramètres suivants de la fonction gamma (conformément à la norme ICC v4): ySLE = (aSLE.cSLE+bSLE \ A- ) gaMMaSLE ±CSLE , Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : caseColorChecker #19 ; cECR RE =X1 et (RVB)écran =(1,1,1). L'utilisateur est invité à pointer la case ColorChecker #19 (blanche) dans une salle noire et le procédé mesure automatiquement les réponses de la sonde de luminance extérieure CSLE. Comme le procédé connaît E(À)BL MAX et RCC #19 (À), il calcule le taux de perte de luminance après réflexion Rcc #19 lié à la case ColorChecker #19: RCC #19 = (780 "1 (E(À)BL MAX . RCC #19(À) .yStdObs(A).dÀ r7S0 Mn. E(À)BL MAX -13SO -1380 mtni yStdObs(À).dÀ ) ; Ensuite, le procédé diminue le rétroéclairage jusqu'à la borne supérieure de détection de la caméra dite « 1 » ou « rétroéclairage » { caseCC=19 ; cECR RE=X3 et (RVB)écran =(1,1,1) }. La diminution du rétroéclairage s'opère en si possible N pas (N=10 à 16, à défaut 5, voire 2 au 20 minimum). Comme les étapes précédentes ont permis de déterminer la luminance produite par l'écran YEcR en fonction des paramètres d'entrée CECR RE et cECR RN/B, le procédé obtient de la sorte N points de mesure (Xi, Yi) de poids VVi=1 où : 25 _ = cSLE ; (cSLE = sortie du dispositif) Xi _ Yi = ySLE = yECR .R cc#19 avec : yECR = yECR RE. (yECR R yECR V + yECR B) yECR RE = (al ECR RE. cECR RE+bi ECR RE)AgammaECR RE+cl ECR RE 30 _ yECR R/V/B = (al ECR RN/B.cECR R/V/B+biECR R/V/B)AgammaECR RN/B+ciECR RN/B La solution de régression correspond au minimum du « Chi 2 » : ,x2 SLE = 1 wi Yi .....((aSLE.xi+bSLE)AgarnmasLE+csLE))A2 Le procédé détermine alors par l'algorithme des moindres carrés les inconnues gammaSLE, aSLE, bSLE et bSLE. Un point important pour garantir la convergence est de procéder en deux passes : sLE 1) une régression linéaire (sur base logarithme) avec bsLE=c=0 2) la résolution d'un système non linéaire surdéterminé avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton) pour minimiser davantage le « Chi 2» x,2 sLE avec bsLE et cSLE. Etape t) : Première passe (linéaire) de calcul de la matrice de transfert colorimétrique de la caméra des CCEI RI /B linéaire e [0;1] vers CIE XYZ (cette matrice est fonction des filtres de couleurs placés devant les photodétecteurs) Le procédé va déterminer la matrice de transfert colorimétrique de la caméra [mcE1], de dimension(3,3), des CcEl RN/BlinéaireE [0;1] vers CIE XYZ : XCEli CCEIRlinéaire 1 mc = v [YCE[t] ucEicEi é vi :afre i] ZCEli CCElBitnéatre Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; CECR RE=xi ; (RVB)écran =(1,1,1) }. Ensuite l'utilisateur est invité à pointer successivement les cases ColorChecker 01 à 24 la mire de couleurs dans une salle noire et le procédé mesure 30 automatiquement les réponses (RVB) caméra. Le procédé dispose de la sorte de 24 mesures (N=24). Pour chaque mesure i : XCE1 L CCEIR1i: [ t] = [McEI].[CCEIViinéa re : (1) ZeKli CCEISi!néaLre En posant : y CE11 yCE 12 CE13 p\AcEI, = [ yCE14 y CE15 CEM , yeE.17 yeE18 ye.E19 L'équation (1)devient : xCEli = y CEItcCEI R linéaire i yCE12.cCEI V linéaire i yCE13.cCEI B linéaire i ; yCEli = yCE14.cCEI R linéaire i yCE15.cCEI V linéaire i y CE16.0CEI B linéaire i ; zCEli = yCE17.cCEI R linéaire i yCE18.cCEI V linéaire i yCE19.cCEI B linéaire i ; Posons les variables « PhiCEI » (DcE1(i,k,XN/Z) telles que : = CCEI R linéaire i (DCE1(i,27x) =V u linéaire i (DcEi(i,3,x) = L, e,CEI B linéaire i aDcEl(i,k,X) = 0 pour k=4 à 9 (DCEIG. I K Y)= 0 pour k=1 à 3 (DcE1(i,4,y) =CCEI linéaire i ajcE1(,5,y) = CCEIV linéaire i oCE1(0,y) = f-sCEI B linéaire i OcEl(i,k,Y) = 0 pour k=7 à 9 CEI(iky) Y) = 0 pour k=1 à 6 (DCE1(i,7,y) = R linéaire i ocEi(j,8,z) = cCEI R linéaire i (DCE1(i,9,z) = cCEI V linéaire i L'équation (1) devient : xem SOMME(k = 11 9; yCEIk.PhiCEI(i,k,X)) [YCEII,]=[SOMME(k = 1 à 9; yCElk.PhiCE.1(i,k,Y):)] ZCEit SOMME(k = 1 à 9;),CElk.PWCE10,fr,Z)) Pour chaque mesure i, comme C(x)cEI = Rcc#im.Ec .A)BL VXG[380nm ;780nmi, l'équation (1) devient : _yoyEcR RE(x1). r72.13; (E(À)BL MAX.RCC #i(2\,) ( .xStdObs,, A) dÀ)=EZ-,1( yCElk.(DCEI(i;k;x)) / 380 Irti. '`- CEI 0- yCElk.d) ,k,y)) ,:yoyEcR RE(x i) 780 _ (E(À) - R (2^,).Y BL MAX CC #i StdObs(À). d À)=21, 380 mn (A) dA)=11,1( yc k.( ,KZ)) _cyoyEcR RE(xl). f 780 2-1m (E BL MAX. Roc #im.zStdObs ,- - . El DOEI(i 3eo >1»1 \ j avec : 0/0yECR RE (Xi) = ((al ECR RE.xi +biECR RE)A gammaECR RE+cl ECR RE)NECR BL MAX, Posons le vecteur 0 de dimension (9,1) tel que (ocE5T = (ycEli ,ycEi9). Posons le vecteur Y de dimension (3N,1) tel que : Y(1) = 00yECR RE (Xi) r730 Ja80 E(À)BL MAX.RCC #1(x).xStdObs,' v\) dÀ ) y(2) = yoyECR RE (Xi) f730 eun E(À)BL MAX.RCC #1(1).yStdObs,, çA) dÀ ) - 'a'ElD nen Y(3) = yoyECR RE (xi \ 1- J300 7Z0 21-un ( E(À)BL MAX.RCC #1(2\,).zStdObs ,' v\) dÀ ) / - nan Y(3.(N-1)+1) = yoyECR RE (XI) C30 mn ( E(mBL MAX. RCC #N(4.xStdObs(À).dÀ ) JSM nm Y(3.(N-1)+2) = yoyECR RE (Xi) r7N n7771( E(mBL MAX. ROC #N(x).yStdObs.% Y(3.(N-1)+3) = cyoyECR RE (Xi) y') dÀ ) M71 C60 S0 nen te E(À)BL MAX. RCC #Nm.zStdObs-^,, ' . v) dÀ ) 13mn L'équation (1) peut s'écrire sous la forme suivante avec comme seules inconnues lesOcElk : Yi = Eir:=1(0CEb K (DcEl(i,k,X/Y/Z)) (3N équations)25 Vu la forme linéaire vis-à-vis des OcEik de l'équation, le procédé va utiliser l'algorithme des moindres carrés linéaires pour les déterminer. Posons la matrice jacobienne J de dimension (3N,9) : Phi:CE:1 = >1' k = t,.7( P ET@ = 1. k, X.) .PiliCEI,(1 1',k = 9' ,PhirE.10', k.' x P hi:CE7 k J = [ ,P hiCE10,' = PitirE10' Phie 0',,k = PhirE1,(i', k, 27) .PhiCE10,, = 9, r) PhirE10 = N, k = Phi CE(i = Nk. Z) .PhirETO = N = 9, 2.,/) Posons la matrice 1/\/ de dimension (3N,3N) des poids de chaque mesure VVi XN/Z associé à l'écart-type ETRoeleiXN/Z (i=1 à N=24): - ETRoe XN/Z= r780 21:11 E(À)BL MAX. RCC #i (À) (YoETRcc #i (À)./x ykstdobs dÀ ) j3EO /rK7so E(À)Bi_ MAX x/y/zStdObs(À).dÀ ) ; J:3 - Wi XN/Z = 1/ ((ETRcc X/Y/Z)^2) ; wiX o 0 0 WlY 0 0 0 W,1 - 1/\/ = [ 0 0 WNX 0 0 0 0 0 IVNY 0 0 0 WN.Z. La solution de régression emin correspond au minimum du « Chi2 » (x2) : x2 = VVi XN/Z . (Yi - Ez.icocElk ; (DCEI(-, k,XN/Z)) )^2 ) (3N termes) Ce minimum estatteint pour :eminjr;vv; j)-1; ; T d'où les neuf valeurs de la matrice de transfert colorimétrique de la caméra. Etape g) Première passe (linéaire) de calcul des spectres des 3 primaires de l'écran rouge E(A)R MAX, vert E(A)V MAX et bleu E(A)B MAX Le procédé va interpoler le spectre de chaque primaire de l'écran à l'aide d'une fonction spline cubique scellée et de (NECR RN/B +1) points d'interpolation de coordonnées [Xk'ÀklYk=YECRR/V/BMAXk] pourk=0 àNECR RN/B - ÀO = 380 nm, ÀNECR RN/B =780 nm - pente nulle aux extrémités du spline : p0=pNECR RN/B=0 Les yECRRN/BMAX ksont les inconnues que le procédé va déterminer. Nota : la luminance est prise au maximum de chaque primaire, à l'instar du blanc : E(À)BL MAX = E(À)R MAX + E(À)V MAX + E(À)B MAX ; Nota : NECR fflii3 peut être réduit à une vingtaine de points environ via le choix d'une suite discrète de Àk choisis judicieusement par rapport aux points sensibles du spectre des primaires de l'écran, par exemple : 380, 410, 430, 440, 450, 460, 470, 485, 500, 515, 530, 550, 570, 585, 600, 615, 630, 670, 700, 740, 780 (nm). Cette tactique permet de réduire le nombre de flashs de couleurs à générer, donc le temps passé par l'utilisateur. Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; CECR RE=X1 ; (RVB)écran =(1,1,1) }. L'utilisateur est invité à pointer successivement les cases ColorChecker#01 à #24. Pour chaque case, le procédé génère N flashs colorés répartis uniformément dans le gamut de l'écran (typiquement N=45 ou 153 en fonction de la précision choisie, voir algorithme de génération par la suite) et le procédé mesure automatiquement les réponses (RVB)cEi. Au total, le procédé génère de la sorte 24.N mesures (i=1 à 24.N). Le procédé utilise un algorithme de stabilisation d'image. On utilise une image avec au moins un point de contraste. Nota : Le procédé génère des flashs de -50 ms, soit -50 flashs par seconde pour être compatible avec la réactivité des composants électroniques courants. L'algorithme pour générer les N flashs colorés est tel qu'illustré Figure 3. Nota : L'écran est capable de générer un grand nombre de couleurs, typiquement Ne [2 ; 224=16,7 millions] pour les systèmes RVB en 3x8 bits. Comme le spectre de chaque primaire de l'écran est une fonction spline scellée à pente nulle aux extrémités, le spectre de chaque primaire de l'écran peut s'écrire sous une forme linéaire grâce aux formules d'interpolation de Newton et d'Hermite : E(À)R/V/B = EoCIR RA,VE ECR RN/BMAXk. (DECR RN/B ) : (1) avec : -1 entre 1 et NECR RN/B tel que : 2',1_1 Posons la matrice jacobienne J de dimension (3.24.N,N*EcR) égale à: E C R = k = 1, AD Phi ECF 1,,K) [ Pk- r phiw-ECR(i, Pkf.ECRUk X) Phi E C Y) i4 ECR( .7) Phi ECR,X) Ph r = N 177 CR Y) Pht EC ECRiz,) = :Lk le-, X) Ph 'EGR« = 1,k = ECR.'rj Phi ECR = 24N,k= izj Phi :Ir 24N, *,,,,Z) Phi ECR(t =-7 24,k N F Posons la matrice W de dimension (3.24.N,3.24.N) des poids de chaque mesure égale à la matrice identité (le procédé ignore à ce stade l'écart type pour chaque flash de couleur). La solution de régressionemin correspond au minimum du « Chi2 » (x2) : x2 = X/Y/Z. (x/y/zcEti_zi,,,,:fr 03*EcRk*.a)*EcR(-, k*,X/Y/Z)) )^2)) ; Nota : somme de 24.3.N termes Ce minimum est atteint pour : Omin = (JT.J)-1.JT.Y , d'où les valeurs des points d'interpolation des 3 primaires de l'écran de coordonnées : [xk=Àk sk=yECRR/V/BMAXk] our k=0 àNECR RIV/B. Etape h) Seconde passe (non-linéaire) d'ajustement entre la matrice de transfert colorimétrique de la caméra TAICEI) et les spectres des 3 primaires MAX, vert E(A)1/ MAX de l'écran rouge E(A)R et bleu E(A)B MAX Le procédé utilise les 24.N mesures réalisées à l'étape g). Le procédé va procéder à la résolution d'un système non linéaire surdéterminé avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton) pour minimiser davantage le « Chi 2 »de l'étape g) : x2 = -Ex_2(Wi XN/Z. (X/Y/ZCEli ,""EcR re*ECRk*.(*ECR,, k --1 k*,X/Y/Z)) )^2)) ; Nota : somme de 24.3.N termes Le vecteur inconnu est la concaténation des OcElk et des 0*EcRk*, de dimension (n=9+N*EcR,1). Le vecteur initial pour l'algorithme de Levenberg-Marquardt est celui obtenu à l'issue de l'étape g). Nota : la dimension de l'espace vectoriel d'origine (n=9±N*EcR) doit être strictement inférieur à la dimension de l'espace vectoriel de destination (m=3.24. N). La matrice W des poids des mesures, de dimension (3.24.N) est définie avec les valeurs des splines des primaires de l'écran déterminées à l'issue de l'étape g) : - ETRoe X/Y/Z= f78() (E(À) i . Rcc #i(À) YoETRcc #1 (À). x/y/zstdobs(.^) .. d À ) 380 run ' ,rso m'un nm E (À ) x/y/zstdo(bs, s À) dÀ ) ; 380 _ %yECR RE (Xi) ( EwRi EmVi E(2)Bi) ; _ %yECR RE (Xi) = ((a.iECR RE.xi +blECR RE)AgammaECR RE+ci ECR RE)/yECR BL MAX ; EeN/B =cECR RiviBlinéaire N EocR,Rivis ECR R/V/BMAXk (DECR RN/B (1,LÇA ) ) ; cECR RN/B linéaire = (aECR RN/B.cECR RN/Bi+bECR RN/B) A gammaECR RN/B+C ECR RN/B - Wi X/Y/Z = 1/ ((ETRcc X/Y/Z)A2) ; 0 0 Wi Y 0 o 0 VII 1 27 - W = [ O 0 W24.N X 0 0 0 0 W24NY 0 0 0 Wa4N Z Etape h9: Itération des étapes c) à h) Le procédé utilise les mesures déjà réalisées par l'utilisateur, en temps masqué pour ce dernier. Avant de relancer les calculs des étapes c) à h), le procédé met à jour les valeurs théoriques de sortie « usine » suivantes avec les résultats de l'étape h) - la luminance maximale blanche de l'écran (yECR BL MAX en cd/m2) ; 20 - le spectre lumineux blanc associé E(À)BL MAX obtenu pour RVB=(1,1,1) et cECR RE=1 ; - les valeurs pECR RN/B = yECR RN/BNECR BL (%). L'itération est acceptée si l'ensemble des « Chi 2 » est réduit : 25 - Etape c) : x2 CE I RN/B=T.,:q w (aCEI R/V/B,xi+bCEI RN/B)AgammaCEI RN/B÷ cCEI Rtv/B))A2 (aECRRExi+bECR RE^ A - Etape d) : x2 cECR RE))A2 ECR RE (yi ) gammaECR RE + 15 - Etape e) : x2 ECR R/V/B wi.(yi (aECR R/V/B.xi+bECR RN/B)Ag am maECR RN/B+cECR RN/B))A2 Etape h) : x2 yis.z (VViVVIY/Z. (X/Y/ZcEli .,:f (0,,EcRk,,.(1),,EcR(j,k*,X/Y/Z)) )^2)) ; L'itération prend fin pour un nombre arbitraire d'itérations (cf. temps de calcul) et/ou lorsque l'ensemble des « Chi 2» sont inférieurs à un seuil de précision arbitraire (par exemple, la précision des calculs). Etape i) Détermination des 2 fonctions de transfert qui permettent de passer de la dioptrie d de la lentille autofocus à 1) la mesure de distances dans le plan focal image (ex le diamètre d'une tâche) ; 2) la distance entre le foyer de la lentille de la caméra et le plan focal image Le procédé éteint l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { caseColorChecker#19 ; CECR RE=X1 ; (RVB)écran =(1,1,1) }. Nota : réalisation des opérations dans une salle intérieure avec un éclairage non éblouissant, contrairement aux étapes précédentes réalisées dans une salle noire, à défaut dans un très environnement très faiblement éclairé et stable. L'utilisateur pose la feuille de papier quadrillé sur un support horizontal (par exemple une table bien plane) et il indique au procédé le pas du papier quadrillé utilisé (par exemple 5x5mm, si possible à fort contraste) Ensuite, l'utilisateur est invité à déplacer le dispositif la tête en bas, à différentes hauteurs par rapport au papier quadrillé, le plus à horizontal possible. A chaque pas, le procédé attend la mise au point de l'autofocus, puis il stocke le paramètre d'entrée de l'autofocus (dioptrie d) et il analyse l'image pour déterminer la position des croix du quadrillage, par exemple avec l'algorithme décrit à la fin de ce chapitre. Nota : le procédé tient compte de l'effet fish-eye éventuel de la lentille et stocke les coordonnées de tous les points d'intersection du quadrillage avec-sans effet fisheye pour corriger cet effet fisheye par la suite au moyen d'un procédé de morphing La position des croix du quadrillage permet de calculer : 1) la distance de bord à bord de l'image captée par la caméra dans le plan focal image (cf. nombre réel d'unités de quadrillages en abscisse et en ordonnée), d'où la capacité pour le dispositif, par simple homothétie, à transformer une distance entre 2 pixels à l'image en une distance métrée dans le plan focal image pour un paramètre d'entrée de l'autofocus - donné (dioptrie d). Nota : le procédé utilise une interpolation du nombre réel de quadrillages en abscisse et en ordonnée pour un paramètre d quelconque compris entre 2 prises d'étalonnage dl et d2 (par exemple, une fonction gamma comme à l'étape c) ou une fonction spline cubique scellée comme à l'étape g)) ; 2) la distance et l'inclinaison de la caméra par rapport à la feuille quadrillée, par exemple à l'aide de l'algorithme suivant issu du cours de Ernst Hairer, université de Genève, Suisse, section mathématiques: Les croix du quadrillage constituent N points dont les coordonnées dans le référentiel orthonormé 3D du papier quadrillé sont (xk, yk, zk=0) et dont les coordonnées dans le plan image de la caméra sont (uk, vk). Les 7 inconnues que le procédé va déterminer sont les suivantes. Elles forment par concaténation le vecteur x Cin=7: - (xCEI,yCEI,ZCEI) la position du foyer de la caméra dans le référentiel orthonormé 3D du papier quadrillé ; - le vecteur (a,b,c) qui donne la direction de vue avec norme correspondant à la distance du plan du film ; -Ol'angle de rotation de la caméra autour du vecteur (a,b,c). Considérons la base orthogonale dont les deux vecteurs et .e./ engendrent le plan image, riétant horizontal et vertical : -4 [ a ] ; oeac b2; o [ b -, V( b2)5Ca2-4-2)2+ c2) a2 b2 Alors, le vecteur dans l'espace qui montre du centre (XCEIsCEI,zCEI) de la lentille vers le point (uk, vk) sur le film est donné par : wkia H Alpha. k co.s e ta sin. the ta wk- = [ wk2]= [ b + cx,k i.-;+ 13k D# OÙ [ = [ Beta k -sint,ieta cos the ta 1vk v..rk3 -+ Les conditions à satisfaire sont que les vecteurs wket (xk-xcEisk-ycEl,zk-zcEl) 15 soient colinéaires, d'où les 3 conditions suivantes à minimiser pour chaque k: wkixk-NCEI wk2. (Ek - zCE1.) - wk.3.-(yk - yCEI) [ wk2 Ix [ yk, - yCEI = [ tvk 3. (Xlc - xCEI) - (zk - zCEI) ] : (1) wk3 zk - zCEI 14,7k1.(yk-- yCEI) - wk2.(xk- xCED Posons F : x C1rl=7-> F(x) E9r=3N tel que : F(x) 35(k- 1)) wk2.. (zik - zCEI) - wk3. (yk [ F (x') (35 - +1 j = [ wk3.(xk - xCEI) - - z.CED F(x) (3.(k- 1) , wk1.(.7,,k - yCEO .(x.k - .xCE.0 20 Pour N>2, m>n et il s'agit d'un système non linéaire surdéterminé. Le procédé va minimiser la norme euclidienne de F(x)avec l'algorithme de Levenberg-Marquardt (mélange entre l'algorithme du Gradient et celui de Gauss-Newton). De la solution xmin découlent les valeurs recherchées, à savoir la distance et l'inclinaison de la caméra par rapport à la feuille quadrillée. Nota : Un algorithme possible pour déterminer la position des croix du quadrillage dans le plan image est le suivant : A) Conversion des couleurs en niveaux de gris inversés . Pour analyser l'image, le procédé va t-convertir les couleurs de l'image en niveaux de gris en inversant noirs et blancs avec la formule : Col = 1 R . 0,3 + cv . 0,59 ± U-B . 0,11) Nota : CGI = 0 pour un blanc pur, CI = pour un noir pur Nota : le procédé effectue la conversion CGI à la demande par économie de temps de calcul B) Détermination des coordonnées (xImG,y11\4G) du point « Op » (p pour provisoire) et de la direction de la ligne du quadrillage rencontrée à « Op » : Le procédé analyse le CGI des pixels de l'image en partant du point central de l'image (xlIvIG=0, ylmG=0) et en remontant l'axe des ordonnées vers les y croissants (xlmG=0, ylmG>0) jusqu'à passer un pic de CGI (voir figure 1) qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « Op » d'ordonnée ylmG0p. La largeur du pic selon y1MG est appelée dylmG0p. Nota : si le pic a commencé à ylmG=0, le procédé balaie les CGI des pixels de part et d'autre de 0 pour localiser le milieu du pic. Nota : si la largeur du pic est très grande (par exemple un quart de la hauteur de l'image), cela signifie que la ligne est parallèle à l'axe des ordonnées yil\AG. Dans ce cas, le procédé passe directement à C). Ensuite, le procédé balaie les CI des pixels de part et d'autre du point « Op » selon l'axe des abscisses pour déterminer le milieu du pic de CGI. Le procédé obtient alors l'abscisse du point « Op» xlmG0p, ainsi que la largeur du pic selon XIMG appelé dx1mG0p. Ensuite, le procédé va déterminer la direction de la ligne rencontrée au point «Op » avec les formules suivantes (voir figure 2) : - tan(00p) = dylmG0p /dxImGp+ ; - = cos(eOp) + sin(00p) Ey ; - 1/0p = -sin(00p) Ex. + cos(00p) Ey ; Si -45°5_00p._ 45°, la ligne est dite horizontale ; elle est appelée « H+1/2 » Sinon, la ligne est dite verticale ; elle est appelée « V+1/2 » si 00p est négatif et « V-1/2 » si 00p est positif. Dans ce cas, on effectue une rotation de u-W.? et vop pour que UOpt soit le vecteur horizontal et1,77e0 le vecteur vertical. La largeur de la ligne vaut LOp = racine carrée ((sin(00p). delG0p)12 + (cos(00p). dylmG0p)1%2). C) Détermination des coordonnées (ximG,y1mG) des points +1/2pH, +1/2pB, - 25 1/2pH et -1/2pB Posons deux points « OpH » (haut) et « OpB » (bas) de part et d'autre de la ligne rencontrée à Op tels que (voir figure 3) : OpH = 0+ + 510p.vop ; 20 30 OpB = 0+ - 5.L0p.vor; ; En partant de OpH, le procédé balaie les CGI des pixels vers les Lioï..; croissants jusqu'à passer un pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2pH ». En partant de OpB, le procédé balaie les CGI des pixels vers les U0 croissants croissants jusqu'à passer un pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2pB ». En partant de OpH, le procédé balaie les CGI des pixels vers les LITOp décroissants jusqu'à passer un pic de Col qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2pH ». En partant de OpB, le procédé balaie les CGI des pixels vers les ii-e9 décroissants jusqu'à passer un pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2pB ». Nota : lorsque le procédé balaie les pixels parallèlement à une ligne du quadrillage, tous les N pixels (par exemple N=10), le procédé recale la position du milieu de la ligne du quadrillage et les vecteurs F., et 17 en déterminant le milieu du pic de CGI selon la direction TI; D) Détermination des coordonnées (xlmG,y1mG) des points +1/2pHD, +1/2pHG, puis de +1/2pD, +1/2ppD, puis de +1/2pG, +1/2ppG, puis de +1/2p, et +1/2pp Posons deux points « +1/2pHD » (droite) et « +1/2pHD » (gauche) de part et d'autre de la ligne rencontrée à +1/2pH tels que (voir figure 3) : (+1/2pHD) = (+1/2pH) + 5.L.0 ; (+1/2pHG) = (+1/2pH) - 5.L.0 1/2pii ; avec L la largeur de la ligne rencontrée au point +1/2pHD (cf. formule à la fin de l'étape B) En partant de +1/2pHD, le procédé balaie les CGI des pixels vers les V0p décroissants jusqu'à passer un pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2pD ». Le procédé continue de balayer les CGI des pixels vers les vo r décroissants jusqu'à passer un deuxième pic de 0G1 qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2ppD ». En partant de +1/2pHG, le procédé balaie les CGI des pixels vers les Fa décroissants jusqu'à passer un pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2pG ». Le procédé continue de balayer les CGI des pixels vers les IlOp> décroissants jusqu'à passer un deuxième pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « +1/2ppG ». Nota : lorsque le procédé balaie les pixels parallèlement à une ligne du quadrillage, tous les N pixels (par exemple N=10), le procédé recale la position du milieu de la ligne du quadrillage et les vecteurs .54 et t7en déterminant le milieu du pic de CGI selon la direction ri; Le procédé calcule alors les coordonnées de +1/2p = milieu de +1/2pD et +1/2pG. Le procédé calcule alors les coordonnées de +1/2pp = milieu de +1/2ppD et +1/2ppG.30 E) Détermination des coordonnées (xlmG,yelG) des points -1/2pHD, -1/2pHG, puis de -1/2pD, -1/2ppD, puis de -1/2pG, -1/2ppG, puis de -1/2p, et -1/2pp Posons deux points « -1/2pHD » (droite) et « -1/2pHD » (gauche) de part et d'autre de la ligne rencontrée à -1/2pH tels que (voir figure 3) : (-1/2pHD) = (-1/2pH) + 5.L.0 +1/2pH ; (-1/2pHG) = (-1/2pH) - 5.L.0 -F1/2zpH avec L la largeur de la ligne rencontrée au point -1/2pHD (cf. formule à la fin de l'étape B) En partant de -1/2pHD, le procédé balaie les CG1 des pixels vers les 7G décroissants jusqu'à passer un pic de CG1 qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2pD ». Le procédé continue de balayer les CG1 des pixels vers les décroissantsdécroissants jusqu'à passer un deuxième pic de CGI qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2ppD ». En partant de -1/2pHG, le procédé balaie les CG1 des pixels vers les ---e-)1/ décroissants jusqu'à passer un pic de CG1 qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2pG ». Le procédé continue de balayer les CG1 des pixels vers lesfrp' décroissants jusqu'à passer un deuxième pic de CG1 qui indique la présence d'une ligne du quadrillage. Le milieu de ce pic est appelé point « -1/2ppG ». Nota : lorsque le procédé balaie les pixels parallèlement à une ligne du quadrillage, tous les N pixels (par exemple N=10), le procédé recale la position du milieu de la ligne du quadrillage et les vecteurs if' et t! en déterminant le milieu du pic de CG1 selon la direction Fi"; Le procédé calcule alors les coordonnées de -1/2p = milieu de -1/2pD et -1/2pG. Le procédé calcule alors les coordonnées de -1/2pp = milieu de -1 /2ppD et -1 /2ppG. F) Détermination des coordonnées (xlmG,yimG) des points (0;0), (+0,5,0), (-0,5;0), (+0,5;0), (-0,5;0) et des pas horizontal et vertical du quadrillage (0;0) est le barycentre de +1/2p, +1/2pp, -1/2p, -1/2pp ; (+0,5;0) est le barycentre de +1/2p, +1/2pp ; (-0,5;0) est le barycentre de -1/2p, -1/2pp ; (0;+0,5) est le barycentre de +1/2p, -1/2p ; (0;-0,5) est le barycentre de +1/2pp, -1/2pp ; Le pas horizontal du quadrillage est égal à la distance entre les points (-0,5;0) et (+0,5 ;0). Le pas vertical du quadrillage est égal à la distance entre les points (0;-0,5) et (0;+0,5) G) Détermination des droites horizontales Hvi/2et verticales Vuj/2 La détermination d'une droite consiste àdéclarer son existence dans la liste des droites constituant le quadrillage et à instancier la liste de points d'intersection du quadrillage qui constituent ses segments. Le procédé part des points +1/2pH et +1 /2pB et balaie vers les croissants pour détecter un pic de 0G1 qui correspond à la ligne de quadrillage suivante. Pour optimiser le temps de calcul, le procédé utilise le pas horizontal du quadrillage déterminé à l'étape F) pour conjecturer la zone la plus probable de recherche, et ainsi de suite jusqu'à arriver au bord de l'image. Chaque ligne trouvée donne lieu à la déclaration d'une ligne verticale appelée « Vuj/2 » (exemple : V(+1/2)). Chaque point de la droite horizontale parcourue, situé à l'intersection avec une verticale, est ajouté à la liste de points déterminant la droite horizontale. Nota : Un point important pour contrer l'effet « fish-eye » de courbure de toute lentille consiste à recaler à chaque verticale trouvée la position de la ligne horizontale que le procédé parcourt et les vecteurs il et Ilen déterminant le milieu du pic de 0G1 selon la direction i7 ; Puis, le procédé part des points -1/2pH et -1/2pB et balaie vers les e décroissants pour détecter un pic de CG1 qui correspond à la ligne de quadrillage précédente. Pour optimiser le temps de calcul, le procédé utilise le pas horizontal du quadrillage déterminé à l'étape F) pour conjecturer la zone la plus probable de recherche, et ainsi de suite jusqu'à arriver au bord de l'image. Chaque ligne trouvée donne lieu à la déclaration d'une ligne verticale appelée « Vuj/2 » (exemple V(-1/2)). Chaque point de la droite horizontale parcourue, situé à l'intersection avec une verticale, est ajouté à la liste de points déterminant la droite horizontale. Nota : A nouveau, un point important pour contrer l'effet « fish-eye » de courbure de toute lentille consiste à recaler à chaque verticale trouvée la position de la ligne horizontale que le procédé parcourt et les vecteurs cf et 17 en déterminant le milieu du pic de CG1 selon la direction 7; Ensuite, de la même façon, le procédé va parcourir chaque ligne verticale trouvée pour déterminer les points la déterminant, ainsi que les intersections avec les lignes horizontales. Puis, de la même façon, le procédé va parcourir chaque ligne horizontale trouvée pour déterminer les points la déterminant, ainsi que les intersections avec les lignes verticales. Et ainsi de suite jusqu'à avoir toutes les lignes de l'image complètement déterminées par des listes de points d'intersection. H) Correction de l'effet fish-eye de la lentille (décourbure des droites) Le procédé parcourt chaque droite en partant du centre de l'image, vers la droite et la gauche pour les droites horizontales et vers le haut et le bas pour les droites verticales. Pour chaque droite horizontale, le procédé prend le point au centre, le premier point à droite et le premier point à gauche. Si les 3 points ne sont pas alignés, ils forment un triangle isocèle de hauteur h attestant d'une courbure de la droite due à l'effet fish-eye de la lentille. La correction à apporter est égale à la hauteur du triangle h pour les deux points de part et d'autre du point central (les coordonnées corrigées sont notées avec une étoile, la correction est effectuée dans la direction opposée du centre de l'image). Et ainsi de suite pour les points suivants de la droite. Le procédé est identique pour les droites verticales en partant toujours du point central. Nota : le procédé stocke les coordonnées de tous les points d'intersection du quadrillage avec-sans effet fisheyepour corriger cet effet fisheye par la suite au 30 moyen d'un procédé de morphing Etape j) : Détermination de l'écartement entre le foyer de la lentille de la caméra et le bord de l'écran (servira à guider la position de la cible lors des futures mesures de couleur, une fois l'étalonnage du dispositif terminé) Le procédé allume l'écran avec les valeurs de la borne inférieure de saturation de la caméra dite « 1 » ou « rétroéclairage » déterminées à l'étape a) : { case ColorChecker#19 ; cECR RE=X1 ; (RVB)écran =(1,1,1) }. Nota : retour à une salle noire, à défaut un environnement très faiblement lumineux et stable. L'écran émissif fournit à l'utilisateur un éclairage d'appoint suffisant pour réaliser ses manipulations. Comme à l'étape i), l'utilisateur pose la feuille de papier quadrillé sur un support horizontal (par exemple une table bien plane) et il indique au procédé le pas du papier quadrillé utilisé (par exemple 5x5mm, si possible à fort contraste) Ensuite, l'utilisateur est invité à déplacer le dispositif la tête en bas, à différentes hauteurs par rapport au papier quadrillé, le plus à horizontal possible. A chaque pas, le procédé attend la mise au point de l'autofocus, puis il stocke le paramètre d'entrée de l'autofocus (dioptrie d) et il analyse l'image pour déterminer la position de la ligne de délimitation entre la zone du plan focal image éclairé par l'écran et celle non éclairée (par exemple, via la technique du gradient de luminance des pixels). A chaque pas, la connaissance paramètre d'entrée de l'autofocus (dioptrie d) permet de connaître la distance du plan focal image à la lentille de la caméra et la distance de la ligne de démarcation de l'éclairage de l'écran par rapport à l'axe de la caméra (cf. étape i)). Avec au moins deux pas de mesure, le procédé en déduit le plan qui constitue la limite d'éclairage de l'écran, relativement vertical par rapport au plan formé par l'écran et la caméra côte à côte. Comme le système connaît la position de l'intersection de ce plan frontière dans un plan focal image de distance connue par rapport au foyer de la lentille de la caméra, le procédé est capable de déterminer l'intersection de ce plan frontière avec le plan relativement horizontal formé par l'écran et la caméra côte à côte et donc la distance entre le bord de l'écran et le foyer de la lentille de la caméra. L'invention est décrite dans ce qui précède à titre d'exemple. Il est entendu que l'homme du métier est à même de réaliser différentes variantes de l'invention sans pour autant sortir du cadre du brevet. 10 REVENDICATIONS1. Procédé d'étalonnage d'un dispositif (10) comportant un écran émissif (11) et un capteur (12) électronique d'images en couleur (tri stimulus), ledit capteur (12) électronique d'images comportant une lentille autofocus, caractérisé en ce qu'il comporte les étapes suivantes : a) mesure de la surexposition dudit capteur (12) électronique d'images relativement à la puissance lumineuse de l'écran (11) ; b) mesure du seuil de détection dudit capteur (12) électronique d'images ; c) calibration des fonctions gammas RN/B dudit capteur (12) électronique d'images ; d) calibration de la fonction gamma du rétroéclairage de l'écran (11) ; e) calibration des fonctions gammas RN/B de l'écran (11) ; f) calcul d'une matrice de transfert colorimétrique dudit capteur (12) électronique d'images des CCEI RAI/E3 linéairee [0;1] vers des coordonnées CIE XYZ ; g) calcul des spectres des trois primaires de l'écran (11) rouge Eve MAX, vert E(À)v MAX et bleu E(À)B MAX ; h) calcul d'ajustement entre la matrice de transfert colorimétrique dudit capteur (12) électronique d'images [McEl] et les spectres des trois primaires de l'écran (11) rouge E(À)R MAX, vert E(À)v MAX et bleu E i) détermination de deux fonctions de transfert permettant de passer de la dioptrie d de ladite lentille autofocus à la mesure de distances dans le plan focal image, et à la distance entre le foyer de ladite lentille autofocus et le plan focal image ; etj) détermination de l'écartement entre le foyer de la lentille et le bord de l'écran (11). 2. Procédé selon la revendication 1, caractérisé en ce que ledit écran (11) et ledit capteur (12) électronique d'images sont situés côte à côte. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'il comporte en outre une étape h') d'itération des étapes c) à h). 4. Procédé selon la revendication 1, 2 ou 3, caractérisé en ce qu'il utilise une mire de couleur et du papier quadrillé. 5. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce que les fonctions gammas de l'étape c) convertissent des CRAllB linéairesE [0;1] captés en valeurs CRiv/Be [0;1] afin de correspondre à la sensibilité de l'oeil humain. 6. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce que la fonction gamma de l'étape d) convertit le paramètre de pilotage du rétroéclairage CECR REe [0 ;1] en valeur linéaire cECR RE linéaire e [0 ;1] correspondant à la luminance effectivement produite par l'écran (11). 7. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce que les fonctions gammas RN/B de l'écran de l'étape e) convertissent les CRN/Be[ ;1] en CR/v/B linéairesE[0 ;1] afin de correspondre à la sensibilité de l'ceil humain. 8. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce que ladite matrice de transfert est fonction de filtres de couleurs placés devant des photodétecteurs. 9. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce qu'il comporte en outre une étape e') de calcul de la borne inférieure de saturation du capteur électronique d'image et de sa borne supérieure de détection sous la forme de valeurs de luminancerelative à celle de l'écran (11), et d'étalonnage d'une sonde de luminance extérieure (si elle existe). 10. Procédé selon l'une au moins des revendications précédentes, caractérisé en ce qu'il est mis en oeuvre pour paramétrer un écran. 11.Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce qu'il est mis en oeuvre pour mesurer des distances au moyen du capteur (12) électronique d'images. 12. Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce qu'il est mis en oeuvre pour corriger l'effet « fish eye » de la lentille. 13. Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce qu'il est mis en oeuvre pour mesurer des inclinaisons au moyen du capteur (12) électronique d'images. 14.Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce qu'il est mis en oeuvre pour effectuer des mesures de couleur. 15.Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce què ledit dispositif (10) est un téléphone mobile. 16.Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce que ledit dispositif (10) est une tablette numérique. 17. Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce que ledit dispositif (10) est un écran d'ordinateur ou une télévision équipé(e) d'une caméra de visioconférence. 18. Procédé selon l'une au moins des revendications 1 à 9, caractérisé en ce que ledit capteur (12) électronique d'images comporte un capteur de luminance.25