FR 3156215 A1 20250606 FR 2313582 A 20231205 Méthode d’ajustement d’un modèle géométrique d’un robot multiaxe équipé d’une caméra , application pour un robot d’application de produit de revêtement et robot configuré pour mettre en œuvre une telle méthode Un modèle géométrique d’un robot multiaxe équipé d’une caméra (30) comprend une première matrice ( et une deuxième matrice de passage Une méthode d’ajustement de ce modèle comprend : viser (1004), avec la caméra disposée dans un premier emplacement (k1), au moins deux points ; déterminer les coordonnées ( de chaque point visé; déplacer (1005) la caméra ; d) viser (1006), avec la caméra disposée dans un deuxième emplacement (k2), les mêmes points (P l ) ; déterminer les coordonnées ( de chaque point visé ; calculer (1008), pour chaque point visé, une différence ( entre ses coordonnées nt dans un repère de base (BF), exprimées en utilisant les première et deuxièmes matrice de passage ( ; calculer (1010, 1018) une valeur globale de différence (F, G). Des variables (X 1 -X 6 , ) de coefficients des première et deuxième matrices de passage ( sont déterminées en minimisant (1012, 1018 ; 1021) la valeur globale calculée à l’étape g). Figure pour l'abrégé : Fig 4 Méthode d’ajustement d’un modèle géométrique d’un robot multiaxe équipé d’une caméra, application pour un robot d’application de produit de revêtement et robot configuré pour mettre en œuvre une telle méthode La présente invention concerne une méthode d’ajustement d’un modèle géométrique d’un robot multiaxe disposé dans un local et équipé d’une caméra. Dans le domaine de l’application de produit de revêtement au moyen d’un robot multiaxe équipé d’un organe d’application, tel qu’une tête d’impression, il importe de pouvoir commander le fonctionnement de cet organe d’application de façon précise, notamment en tenant compte du positionnement de cet organe d’application par rapport à une surface à revêtir. Pour ce faire, il est connu d’équiper un robot multiaxe d’une caméra qui permet de repérer l’environnement du robot, notamment une surface à revêtir, en regard de laquelle doit être positionné l’organe d’application. À ce jour, la position de la caméra dans le local, c’est-à-dire la position et/ou l’orientation d’un repère lié à la caméra dans un repère lié au local dans lequel est disposé le robot, est définie de façon relativement imprécise. En effet, si on utilise un modèle géométrique du robot basé sur l’hypothèse que la liaison entre la caméra et le poignet du robot est connue et rigide et que les axes du robot multiaxe, qui sont généralement au nombre de six, fonctionnent de façon parfaite, alors ce modèle utilisé ne correspond pas à la réalité physique du robot équipé de la caméra, car il existe des tolérance de montage de la caméra sur le poignet du robot, car la position de la caméra sur le poignet peut avoir changé et car le robot lui-même n’est pas parfait. Des problèmes comparables se posent avec des robots équipés d’une caméra et d’un outil qui peut être différent d’un organe d’application de produit de revêtement. Dans le domaine de la chirurgie, il est connu de CN115431278A de calibrer la position du centre d’outil (ou « tool center point » – TCP) d’un repère lié à une caméra embarquée sur un robot multiaxe en partant du principe que ce robot est parfait. Si une telle approche peut être envisagée dans le cas d’un robot de chirurgie, elle n’est pas transposable à un robot industriel, tel qu’un robot d’application d’un produit de revêtement, fonctionnant dans un environnement potentiellement perturbé et à grande vitesse, au point qu’il peut ne pas correspondre à son modèle théorique sur le plan géométrique, cinématique et/ou dynamique. Il est également connu de CN115533893A de déterminer le TCP d’une caméra en imposant à un robot de toucher, au moyen d’une pointe, un profil en forme de sphère. La mise en œuvre d’une telle méthode est complexe et consommatrice de temps. D’autre part, EP1555508A1 enseigne d’utiliser un système de mesure en prenant des photographies d’un objet, tout en faisant tourner une caméra autour d’un axe de son repère, ce qui permet de réduire des erreurs de mesure. Cette méthode est également complexe et consommatrice de temps. CN115741720A enseigne par ailleurs de calibrer la position angulaire d’un capteur en utilisant des calculs basés sur la méthode de Levenberg-Marquardt. Le calibrage ainsi réalisé est limité à l’optimisation du point d’origine de la mesure angulaire des capteurs et fait également l’hypothèse que le robot utilisé est parfait. C’est à ces inconvénients qu’entend plus particulièrement remédier l’invention en proposant une nouvelle méthode d’ajustement d’un modèle géométrique d’un robot multiaxe, en tenant compte de la réalité de l’assemblage entre la caméra et le poignet du robot et de la réalité de la structure du robot multiaxe. À cet effet, l’invention concerne une méthode d’ajustement d’un modèle géométrique d’un robot multiaxe disposé dans un local et équipé d’une caméra, ce modèle comprenant une première matrice de passage entre un repère lié à la caméra et un repère lié à un poignet du robot, et une deuxième matrice, de passage entre le repère lié au poignet du robot et un repère de base lié au local, caractérisée en ce que cette méthode comprend au moins les étapes suivantes consistant à : viser, avec la caméra disposée dans un premier emplacement de prise de vue situé dans le local, au moins deux points d’une cible fixe dans le repère de base ; déterminer les coordonnéesde chaque point visé à l’étape a) dans le repère lié à la caméra qui est dans le premier emplacement de prise de vue ; déplacer la caméra vers un deuxième emplacement de prise de vue situé dans le local, différent du premier emplacement de prise de vue ; viser, avec la caméra disposée dans le deuxième emplacement de prise de vue, les points de la cible fixe déjà visés à l’étape a) ; déterminer les coordonnéesde chaque point visé à l’étape d) dans le repère lié à la caméra qui est dans la deuxième position ; calculer, pour chaque point visé par la caméra aux étapes a) et d), une différenceentre les coordonnées de ce point dans le repère de base, exprimées en fonction des coordonnées déterminées aux étapes b) et e) en utilisant les première et deuxièmes matrice de passage ; calculer au moins une valeur globale de différence au moyen d’une fonction ayant comme variable une différence calculée à l’étape f) ; et en ce que des variables de coefficients de la première matrice de passageet des variablesde coefficients de la deuxième matrice de passagesont déterminées en minimisant la valeur globale calculée à l’étape g). Grâce à l’invention, la détermination des variables des coefficients de la première matrice et des variables des coefficients de la deuxième matrice permet une adaptation efficace du modèle géométrique du robot multiaxe, y compris en ce qui concerne la caméra embarquée sur le poignet de ce robot. Le calcul des variables des coefficients de la première matrice et de la deuxième matrice permet de tenir compte non seulement de l’assemblage entre la caméra et le poignet du robot mais également du dimensionnement géométrique du robot lui-même, notamment grâce à une optimisation des paramètres de Denavit-Hartenberg modifiés. Selon des aspects avantageux mais non obligatoires de l’invention, une telle méthode peut incorporer une ou plusieurs des caractéristiques suivantes prise(s) selon toutes combinaisons techniquement admissibles : - Une première valeur globale de différence est calculée à l’étape g) et en ce que la méthode comprend des étapes postérieures à l’étape g) et consistant à : déterminer les variables des coefficients de la première matrice de passageen minimisant la première valeur globale calculée à l’étape g) ; mettre à jour la première matrice de passageavec des coefficients qui intègrent les variables déterminées à l’étape h) calculer, pour chaque point visé par la caméra aux étapes a) et d), une différenceentre les coordonnées de ce point dans le repère de base, exprimées en fonction des coordonnées déterminées aux étapes b) et e) en utilisant la première matrice de passage mise à jour à l’étape i) ; calculer une deuxième valeur globale de différence au moyen d’une fonction ayant comme variable les différences calculées à l’étape i) ; déterminer les variablesdes coefficients de la deuxième matrice de passage en minimisant la deuxième valeur globale calculée à l’étape j) ; et mettre à jour la deuxième matrice de passage avec des coefficient qui intègrent les variablesdéterminées à l’étape l). - La méthode comprend une étape h’) postérieure à l’étape g) et consistant à déterminer les variables des coefficients de la première matrice et les variablesdes coefficients de la deuxième matrice déterminés, en minimisant une valeur globale unique calculée à l’étape g). - La méthode comprend des étapes, postérieures à l’étape h) ou à l’étape h’) et consistant à : calculer, pour chaque point visé par la caméra aux étapes a) et d), une différenceentre les coordonnées de ce point dans le repère de base, exprimées en fonction des coordonnéesdéterminées aux étapes b) et e) en utilisant la première matrice de passage mise à jour avec les dernières variables de coefficients déterminées pour cette première matrice et la deuxième matrice de passage mise à jour avec les dernières variablesdes coefficients déterminées pour cette deuxième matrice. calculer une valeur globale de différence au moyen d’une fonction ayant comme variable les différences calculées à l’étape n) ; comparer la valeur globale calculée à l’étape o) à une valeur de seuil ; si la comparaison de l’étape p) montre que la valeur globale calculée à l’étape o) est supérieure à la valeur de seuil, mettre à nouveau en œuvre les étapes f) et suivantes ; si la comparaison de l’étape p) montre que la valeur globale calculée à l’étape o) est inférieure à la valeur de seuil, figer les première et deuxième matrices de passage avec les dernières variablesdéterminées pour les utiliser dans le modèle géométrique du robot multiaxe. - La différence calculée à l’étape f), j) ou n) s’exprime sous la forme où est la différence entre les coordonnées du même point numéroté vu par la caméra à partir de positions k1 et k2 ; est un vecteur représentant la position du point l vu par la caméra à partir de la position k1 ; est un vecteur représentant la position du point l vu par la caméra à partir de la position k2 ; est la première matrice de passage; ou est la deuxième matrice de passage. - La valeur globale de différence calculée à l’étape g), k) ou o) est la somme des carrés des valeurs des différences calculées à l’étape f), i) ou l) et s’exprime sous la forme où est la différence entre les coordonnées du même point numéroté l vu par la caméra à partir d’emplacements k1 et k2 ; l est le numéro d’ordre d’un des points visés par la caméra aux étapes a) et d), compris entre 1 et L k1 est le numéro d’ordre du premier emplacement, compris entre 1 et K ; k2 est le numéro d’ordre du deuxième emplacement, compris entre 1 et K ; K est le nombre d’emplacements que peut prendre la caméra pour viser les points de la cible. - La première matrice de passage s’exprime sous la forme où les grandeurs Xi, pour i entier naturel entre 1 et 6, sont les variables déterminées en minimisant la valeur globale calculée à l’étape g) et correspondent pour i entier naturel entre 1 et 3, à des distances entre un centre du repère lié à la caméra et le centre du repère lié au poignet (21); pour i entier naturel entre 4 et 6, à des angles représentatifs d’une orientation du repère lié à la caméra dans le repère lié au poignet. - La deuxième matrice s’exprime sous la forme où les grandeurs sont des angles mesurés sur l’axe A i du robot multiaxe, pour le k-ième point visé par la caméra au cours des étapes a) et d) ; les grandeurs sont les paramètres théoriques de Denavit-Hartenberg modifiés pour l’axe A i du robot multiaxe ; et les grandeurs sont des variables pour les coefficients de la deuxième matrice de passage, déterminées en minimisant la valeur globale calculée à l’étape g). - La minimisation de la valeur globale de différence est effectuée par la méthode des moindres carrés avec résolution d’un système de six équations à six inconnues lors de la détermination des variables de coefficient de la première matrice de passageet résolution d’un système de vingt-quatre équations à vingt-quatre inconnues lors de la détermination des variablesde coefficient de la deuxième matrice de passage, dans le cas d’une méthode telle que décrite précédemment, ou résolution d’un système de trente équations à trente inconnues lors de la détermination des variables de coefficient des première et deuxième matrices de passage, dans le cas d’une autre méthode telle que décrite précédemment. Selon un deuxième aspect, l’invention concerne l’application d’une méthode telle que mentionnée ci-dessus pour l’ajustement du modèle géométrique d’un robot d’application de produit de revêtement équipé d’une tête d’impression ou d’un projecteur de produit de revêtement. Selon un troisième aspect, l’invention concerne un robot équipé d’un outil, d’une caméra et d’une unité électronique de contrôle programmée pour mettre en œuvre automatiquement une méthode telle que mentionnée ci-dessus. L’invention sera mieux comprise et d’autres avantages de celle-ci apparaitront plus clairement à la lumière de la description qui va suivre, donnée uniquement à titre d’exemple et faite en référence aux dessins annexés dans lesquels : FIG. 1 La FIG. 1 est une représentation schématique de principe, sur trois inserts A), B) et C), d’une cible utilisée dans la méthode de l’invention ; FIG. 2 La FIG. 2 est une représentation schématique de principe d’un robot multiaxe conforme à l’invention en cours d’utilisation avec la cible représentée à l’insert B) de la FIG. 1 ; FIG. 3 La FIG. 3 est une représentation schématique de principe de différents repères et matrices de passage utilisés dans la méthode de l’invention ; FIG. 4 La FIG. 4 est un schéma blocs d’une méthode conforme à un premier mode de réalisation de l’invention ; et FIG. 5 La FIG. 5 est un schéma blocs, analogue à la FIG. 4 , pour une méthode conforme à un deuxième mode de réalisation de l’invention. La FIG. 1 représente, sur trois inserts A), B) et C), trois cibles C1, C2 et C3 qui peuvent être utilisées pour mettre en œuvre la méthode de l’invention, avec le robot multiaxe 20 représenté à la FIG. 2 . Chacune des cibles C1, C2 et C3 comprend un certain nombre de points remarquables P 1 à P 4 pour la cible C1, P 1 à P 8 pour la cible C2, P 1 à P 7 pour la cible C3. On note L le nombre de points remarquables d’une cible, ce nombre L étant égal à 8 dans l’exemple de la FIG. 2 . On note l le numéro d’ordre d’un point remarquable P l de la cible ; avec l un entier naturel compris entre 1 et L. Les exemples de cible représentés à la FIG. 1 ne sont pas limitatifs et toute forme de cible est envisageable, pour autant qu’elle permet de repérer au moins deux points remarquables P l . Les points remarquables P l peuvent être les angles d’une figure géométrique ou des tâches dessinés sur une plaque bidimensionnelle, comme représenté sur les inserts A) et B) ou les angles d’une structure tridimensionnelle, comme représenté sur l’insert C) de la FIG. 1 . De façon générale, un point remarquable P l est un point d’une cible qui peut être repéré par une caméra, notamment parce qu’il tranche avec son environnement. À titre d’exemple non limitatif, la FIG. 2 montre l’utilisation de la cible C2 avec le robot multiaxe 20. Le robot multiaxe 20 comprend un bras formé de tronçons articulés entre eux autour de six axes A 1 à A 6 , ainsi qu’un poignet 21 qui forme l’extrémité distale de ce bras. Selon une variante non représentée de l’invention, le poignet peut être articulé par rapport à l’extrémité du bras du robot multiaxe 20, autour d’un septième axe. Une tête d’impression 10, comprenant un corps 12 équipé de buses 14, est montée sur le poignet 21. Elle est destinée à appliquer un produit de revêtement, tel qu’une peinture ou un vernis, sur des objets non représentés, par exemple des carrosseries de véhicule automobile. Les six axes d’articulation A 1 à A 6 permettent de déformer le bras du robot multiaxe 20 pour déplacer la tête d’impression 10 dans l’espace, en regard des surfaces des objets à revêtir. Une caméra 30 est montée sur le poignet 21, en étant fixée sur celui-ci par des moyens mécaniques appropriés, par exemple des vis ou des organes de clipsage. La fixation de la caméra 30 sur le poignet 21 doit être suffisamment solide pour résister aux accélérations subies par la caméra en cours de déplacement par le robot. Avantageusement, la caméra est une caméra de type CCD ou de type Laser avec une ou deux caméra(s), de type binoculaire et/ou profilomètre. Une problématique qui se pose avec un robot multiaxe du type de celui représenté à la FIG. 2 est de connaître la position de la caméra 30, dans un repère de base fixe BF, lié à un local LO dans lequel est disposé le robot multiaxe 20. En effet, un positionnement correct de la caméra 30 dans le local LO est nécessaire pour repérer précisément l’emplacement des objets visés par la caméra. La position de la caméra dans le local LO dépend du positionnement réel de la caméra 30 par rapport au poignet 21 du robot et du positionnement réel du poignet dans le local LO. Le positionnement de la caméra 30 par rapport au poignet 21 peut, en première analyse, être approximé en faisant l’hypothèse que la caméra est immobilisée sur le poignet dans une position théorique correspondant à un plan de montage. Le positionnement du poignet 21 dans le local LO peut, également en première analyse, être approximé en utilisant un modèle théorique de la géométrie du robot multiaxe 20, notamment basé sur des paramètres de Denavit-Hartenberg modifiés ou DHM. Toutefois, les approximations mentionnées ci-dessus conduisent à des erreurs de positionnement, qui peuvent induire des défauts dans l’application de produit de revêtement au moyen de la tête d’impression 10, si celle-ci n’est pas correctement disposée et orientée par rapport à une surface en cours de revêtement. La méthode de l’invention vise à réduire, voire supprimer, les erreurs de positionnement. Elle est mise en œuvre de façon itérative, jusqu’à réduire la valeur d’une fonction représentative d’une erreur de positionnement entre des points remarquables P l observés sur une cible telle que la cible C2. La méthode d’ajustement du modèle de l’invention est mise en œuvre dans un calculateur 40 qui est représenté à la FIG. 2 sous la forme d’un ordinateur et qui communique avec un contrôleur 24 du robot multiaxe 20 disposé dans une base 22 de ce robot. Le calculateur 40 est programmé pour mettre en œuvre automatiquement la méthode de l’invention. Le contrôleur 24 et le calculateur 40 forment ensemble une unité électronique de contrôle du robot multiaxe 20. En variante, les parties 24 et 40 de cette unité de contrôle sont une seule entité physique, qui peut être intégrée à la base 22. On note R Caméra un repère lié à la caméra et R PG un repère lié au poignet 21. La méthode de l’invention prévoit de relever la position d’au moins points remarquables P l de la cible C2 en amenant la caméra successivement dans deux emplacements distincts de prise de vue dans le local LO où est situé le robot multiaxe 20, c’est-à-dire dans deux emplacements distincts par rapport à un repère de base BF fixe par rapport à ce local LO. On note L le nombre de points remarquables P l de la cible C2, ce nombre L étant un nombre entier supérieur ou égal à 2. On note k avec k égal à 1 ou à 2 le numéro de l’emplacement de la caméra lors du repérage d’un point remarquable P l . Un emplacement est défini par un lieu où se trouve la caméra et une orientation de cette caméra lors d’une prise de vue, c’est-à-dire lors du repérage d’un point remarquable. Un emplacement de prise de vue peut également être dénommé position de prise de vue ou emplacement de visée. Pour chaque emplacement de prise de vue k et chaque point remarquable P l , les coordonnées de ce point dans le repère R Caméra s’expriment sous la forme d’un vecteur comprenant une abscisse, une ordonnée et une hauteur de ce point dans ce repère. On note le vecteur exprimant les coordonnées d’un point remarquable P l pour l’emplacement k, exprimée dans le repère caméra R Caméra . En variante, le nombre d’emplacements de prise de vue distincts utilisé est supérieur ou égal à trois. Pour la suite de la présente description, on utilise k1 et k2 pour repérer deux emplacements distincts de la caméra parmi les K emplacements que peut prendre la caméra pour viser les points P l de la cible en étant déplacée par le poignet 21, avec K un entier naturel supérieur ou égal à 2. Avantageusement, le produit des nombres L et K diminué de 6, soit L * K – 6, est supérieur ou égal au nombre de variables du modèle géométrique du robot. Ainsi, la position d’un point remarquable P l dans le repère caméra peut être exprimée sous la forme lorsque la caméra est dans le premier emplacement k1 et sous la forme lorsque la caméra est dans le deuxième emplacement k2. Ces positions des points remarquables P l peuvent également être exprimées dans le repère R PG lié au poignet sous la forme (équation 1) (équation 2) où est une première matrice de passage du repère R Caméra lié à la caméra 30 au repère R PG lié au poignet 21. Cette première matrice de passage peut s’exprimer sous la forme (équation 3) où les seize coefficients de la première matrice de passage sont exprimés en fonction de six variables X 1 à X 6 . Les variables X 1 , X 2 , X 3 correspondant à une translation du centre du repère R Caméra lié à la caméra par rapport au centre du repère R PG lié au poignet 21, alors que les variables X 4 , X 5 et X 6 correspondent à des angles de rotation des axes des deux repères les uns par rapport aux autres. Dans le cas de la matrice présentée en équation 3, il s’agit des angles de la convention de Roulis-Tangage-Lacet. Une autre représentation de ces angles peut être employée. La position des points remarquables P l n’est pas connue dans le repère de base BF mais peut s’exprimer dans ce repère, pour chaque emplacement de prise de vue k, en fonction de leur position dans le repère R PG lié poignet 21 en appliquant à cette position une deuxième matrice de passage du repère du poignet au repère de base. On note ou cette deuxième matrice de passage. Elle est définie pour chaque emplacement de prise de vue k de la caméra, donc du poignet 21. En pratique, la position du repère lié au poignet dans le repère de base dépend de l’emplacement de prise de vue dans lequel se situe la caméra 30. Ce repère lié au poignet lorsque la caméra est dans l’emplacement de prise de vue k peut donc s’exprimer sous la forme . Il existe donc, pour chaque emplacement de prise de vue k, une matrice de passage, du repère lié au poignet au repère de base BF. Cette matrice de passage est notée pour le premier emplacement de prise de vue et pour le deuxième emplacement prise de vue. Ainsi, comme représenté à la figure 3, la position d’un point remarquable P l repérée dans le repère caméra s’exprime dans le repère du poignet en utilisant la première matrice de passage et s’exprime, pour chaque emplacement k1 et k2 dans le repère de base BF en utilisant la deuxième matrice de passage , avec k égal à k1 ou k2. De façon avantageuse, la deuxième matrice de passage s’exprime en utilisant les paramètres de Denavit-Hartenberg modifiés, aussi nommés DHM, ainsi que des corrections associées, ce qui permet de modéliser le robot de sa base 22 à son poignet 21. Les paramètres de Denavit-Hartenberg modifiés sont connus en soi. Ils sont également parfois nommés paramètres de Khalil-Kleinfinger. A titre d’exemple, ces paramètres de Denavit-Hartenberg modifiés peuvent être exprimés, pour chaque axe de rotation A i du robot multiaxe 20 sous la forme : Angle théorique entre les axes et obtenu par la rotation autour de l’axe du repère associé à l’axe Distance théorique entre les axes et obtenu le long de l’axe Angle théorique entre les axes et obtenu par la rotation autour de l’axe Distance théorique entre l’axe et obtenu le long de l’axe On définit également les différences suivantes : Variation d’angle Variation de distance Variation d’angle Variation de distance Dans ces conditions, la deuxième matrice de passage qui permet de passer d’un point du repère R PG lié au poignet dans le repère de base peut s’exprimer sous la forme (équation 4) Cette expression est valable en chaque emplacement de prise de vue k, avec k égal à k1 ou k2. représente la valeur d’un angle de rotation de deux parties du robot multiaxe articulé autour d’un axe A i lorsque le robot multiaxe 20 est dans l’emplacement de prise de vue k. Les coefficients de cette deuxième matrice de passage, avec s et w des entiers naturels compris entre 1 et 4, peuvent s’exprimer comme suit : : matrice de rotation pour et : matrice de translation pour pour Les paramètres de Denavit-Hartenberg modifiés constituent les valeurs théoriques du robot, tandis que les grandeurs constituent des variables pour les coefficients de la deuxième matrice de passage . On note Mod le modèle géométrique du robot multiaxe 20 équipé de la caméra 30. Comme visible à la FIG. 4 , la méthode d’ajustement du modèle géométrique Mod comprend une première étape 1000 d’initialisation dans laquelle le calculateur 40 est démarré. Dans une étape 1002, le modèle géométrique Mod est défini comme un modèle initial Mod 0 qui est lui-même défini par une première matrice de transfert initiale et par une deuxième matrice de transfert initiale . La première matrice de transfert initiale peut être construite par le calcul à partir de la position théorique de la caméra 30 par rapport au poignet 21. La deuxième matrice de transfert initiale peut être construite à partir d’un modèle géométrique théorique du robot multiaxe 20. Le modèle initial Mod 0 est utilisé dans la première partie de la méthode, entre les étapes 1004 et 1012 définies ci-après. Dans une étape suivante 1004, le robot est amené dans le premier emplacement de prise de vue k1 où la caméra 30 vise les points remarquables P l et les repère dans le repère R Caméra lié à la caméra. Dans cette étape 1004, pour chaque point P l , ses coordonnées sont exprimées dans le repère R Caméra sous la forme Ensuite, dans une étape 1005, la caméra 30 est déplacée vers le deuxième emplacement de prise de vue. Lorsque la caméra est dans le deuxième emplacement de prise de vue k2, lors d’une étape 1006 qui suit l’étape 1005, chaque point remarquable P l est visé par la caméra, donc repéré, et ses coordonnées sont exprimées dans le repère R Caméra lié à la caméra sous la forme En pratique, chaque point P l garde la même position, qu’il soit repéré par la caméra à partir de son emplacement de prise de vue k1 ou qu’il soit repéré par la caméra à partir de son emplacement de prise de vue k2. Ainsi, la position de chaque point remarquable P l dans le repère de base BF n’est pas connue mais on sait qu’elle est indépendante de la façon dont elle est observée par la caméra. Pour chaque point remarquable P l, on définit une différence entre ses positions exprimées dans le repère de base BF à partir du repérage effectué dans le premier emplacement de prise de vue k1 et du repérage effectué dans le deuxième emplacement de prise de vue, sous la forme (é quation 5) est donc la différence des coordonnées d’un point remarquable P l exprimées dans le repère de bases BF, à partir de ses coordonnées détectées par la caméra 30 dans le repère R Caméra lié à la caméra, à partir du premier emplacement de prise de vue k1 et à partir du deuxième emplacement de prise de vue k2. Cette différence est calculée par le calculateur 40, pour chaque point remarquable P l , lors d’une étape 1008 de la méthode qui fait suite aux étapes 1006 et 1008. Théoriquement, cette différence devrait être égale à zéro, puisque les points remarquables sont fixes dans le repère de base BF. En pratique, si cette différence est non nulle, on peut faire l’hypothèse que les coefficients de la première matrice de passage initiale ne sont pas exactement représentatifs de la position réelle de la caméra 30 par rapport au poignet 21, notamment à cause des tolérance de fabrication de ces composants et des ajustements réalisés au montage de ces pièces l’une sur l’autre. Lors d’une étape subséquente 1010, une fonction d’erreur F est définie comme étant la somme des carrés et des différences déterminées à l’étape 1008, sous la forme (é quation 6) La valeur de la fonction d’erreur F est une valeur globale de différence entre les coordonnées des points remarquables P l déterminées à partir des deux emplacements de prise de vue. Lors d’une étape 1012, également mise en œuvre par le calculateur 40, la fonction F est minimisée en jouant sur les variables X 1 à X 6 de la première matrice de passage . Cette optimisation est effectuée par la résolution d’un système de six équations non-linéaires, par exemple par la méthode des moindres carrés. Ce système d’équations comporte six inconnues, à savoir les variables X 1 à X 6 . En variante, la résolution du système d’équations non-linéaires est effectuée au moyen de la méthode de Levenberg-Marquardt ou une autre méthode. Dans une étape suivante 1014, la première matrice de passage est mise à jour sous la forme d’une version optimisée de la première matrice de passage , dont les coefficients intègrent les variables X 1 à X 6 déterminées à l’étape 1012 est intégrée au modèle Mod à la place de la première matrice de passage initial En d’autres termes, à partir de l’étape 1014, les calculs effectués par le calculateur 40 tiennent compte d’une version optimisée de la première matrice de passage, à savoir . Dans une étape suivante 1016, une deuxième différence est calculée entre les coordonnées des points remarquables P l sous la forme (é quation 7) Théoriquement, cette différence devrait également être égale à zéro. En pratique, si cette différence est non nulle, et supérieur à la précision de la répétabilité du robot, on peut faire l’hypothèse que les coefficients de la deuxième matrice de passage initiale ne sont pas exactement représentatifs de la structure et du fonctionnement du robot multiaxe 20, notamment à cause des tolérances de fabrication et de l’usure de ses articulations. Lors d’une étape suivante 1018, une fonction d’erreur G est définie comme étant la somme des carrés et des différences déterminées à l’étape 1016, sous la forme : (é quation 8) La valeur de la fonction d’erreur G est une autre valeur globale de différence entre les coordonnées des points remarquables P l déterminées à partir des deux emplacements de prise de vue. Lors d’une étape 1020, également mise en œuvre par le calculateur 40, la fonction G est minimisée en jouant sur les variables 1 de la deuxième matrice de passage , pour chaque emplacement de prise de vue. Cette optimisation est effectuée par la résolution d’un système de vingt-quatre équations non-linéaires, par exemple par la méthode des moindres carrés. Ce système d’équations comporte vingt-quatre inconnues, à savoir les quatre variables . pour chaque axe A i , avec i un entier naturel entre 1 et 6. En variante, la résolution du système d’équations non-linéaires est effectuée au moyen de la méthode de Levenberg-Marquardt ou une autre méthode. La résolution de ce système de 24 équations à 24 inconnues permet de créer une deuxième matrice de passage optimisée Dans une étape suivante 1022, la version optimisée de la deuxième matrice de passage , dont les coefficients intègrent les variables déterminées à l’étape 1020, est intégrée au modèle Mod à la place de la première matrice de passage initiale La méthode de l’invention comprend une étape 1024, mise en œuvre après les étapes 1014 et 1022, de calcul des différences en utilisant le modèle Mod qui intègre les premières et deuxièmes matrices optimisées, c’est-à-dire mises à jour avec les dernières variables X 1 -X 6 , déterminées aux étapes 1012 et 1020. L’étape 1024 consiste également à calculer les valeurs globales de différence à partir des fonctions d’erreur F et G et à comparer ces valeurs globales de différence à deux valeurs de seuil F0 et G0. Si les valeurs globales de différence des fonctions F et G sont inférieures ou égales aux valeurs de seuil F0 et G0, alors on considère que le modèle Mod est correctement ajusté et le nouveau modèle géométrique du robot multiaxe 20 équipé de la caméra 30 est confirmé, c’est-à-dire figé, dans une étape 1026 sur la base des première et deuxième matrices de passage définies aux dernières étapes 1014 et 1022 mises en œuvre. Dans le cas contraire, les étapes 1008 et suivantes sont à nouveau mises en œuvre, sur la base des première et deuxième matrices de passage déjà partiellement optimisées, définies aux dernières étapes 1014 et 1022 mises en œuvre. Les étapes 1008 à 1024 sont alors répétées par itération, jusqu’à obtenir des valeurs des fonctions d’erreur F et G inférieures aux valeurs de seuil F0 et G0. Dans le deuxième mode de réalisation de l’invention représenté à la FIG. 5 , les étapes analogues à celles du premier mode de réalisation portent les mêmes références et ne sont pas décrites en détail. Les étapes 1000 à 1008 de cette deuxième méthode sont identiques aux étapes 1000 à 1008 de la première méthode. La méthode du deuxième mode de réalisation diffère de celle du premier mode de réalisation en ce que la première matrice de passage n’est pas optimisée par recherche d’un minimum de la fonction F considérée seule. Les étapes 1016 et 1018 suivent directement l’étape 1012. Les fonctions d’erreur F et G sont calculées sous la forme d’une valeur globale unique qui est minimisée dans une même étape 1021, laquelle consiste en la résolution d’un système de trente équations non-linéaires à trente inconnues, à savoir les variables X 1 à X 6 , avec i un entier naturel entre 1 et 6. Comme les pour les étapes 1014 et 1020 de la première méthode de la FIG. 4 , la résolution du système d’équations peut être effectuée par la méthode des moindres carrés, par la méthode de Levenberg-Marquardt ou par une autre méthode. Ensuite, une étape 1023 correspondant à la réunion des étapes 1014 et 1022 de la méthode de la FIG. 4 est mise en œuvre, avant qu’une étape de comparaison 1024 soit effectuée, avec le même fonctionnement que pour la première méthode. Lorsque les valeurs des deux fonctions F et G sont inférieures ou égales aux deux valeurs de seuil F0 et G0, le modèle est considéré comme ajusté et le nouveau modèle géométrique du robot multiaxe 20 équipé de la caméra 30 est confirmé, c’est-à-dire figé, dans une étape 1026 sur la base des première et deuxième matrices de passage définies à la dernière étape 1023 mise en œuvre. Dans le cas contraire, les étapes 1008 et suivantes sont à nouveau mises en œuvre, sur la base des première et deuxième matrices de passage déjà partiellement optimisées, définies à la dernière étape 1023 mise en œuvre. La méthode peut être mise en œuvre pour autant qu’au moins deux points remarquable P l sont visés et repérés aux étapes 1004 et 1006. En d’autres termes le nombre de points remarquables utilisés est supérieur ou égal à deux. Il n’est pas nécessairement égal au nombre de points remarquables de la cible. Lorsque plus de deux emplacement de prise de vue sont utilisés pour déterminer la position des points remarquable, la définition des différences et des erreurs est adaptée. En variante, l’une au moins des fonctions d’erreur F et G est construite sans faire intervenir le carré des différences individuelles ou . Ces fonctions peuvent être égales, par exemple, à la somme des valeurs absolues de ces différences ou à une autre valeur calculée à partir de ces différences. En variante non représentée de l’invention, le robot multiaxe 20 peut être équipé d’un organe d’application autre qu’une tête d’impression, par exemple un projecteur pneumatique ou rotatif de produit de revêtement, éventuellement de type électrostatique. Selon une autre variante, le robot multiaxe 20 peut être équipé d’un outil autre qu’un organe d’application d’un produit de revêtement, par exemple un outil d’usinage, un outil de soudage ou un outil de préhension. Toute caractéristique décrite pour un mode de réalisation ou une variante dans ce qui précède peut être mise en œuvre pour les autres modes de réalisation et variantes décrits précédemment, pour autant que techniquement faisable. Méthode d’ajustement d’un modèle géométrique d’un robot multiaxe (20) disposé dans un local (LO) et équipé d’une caméra (30), ce modèle comprenant une première matrice ( de passage entre un repère (R Caméra ) lié à la caméra et un repère (R PG ) lié à un poignet du robot, et une deuxième matrice ( , de passage entre le repère lié au poignet du robot et un repère de base (BF) lié au local, caractérisée en ce que cette méthode comprend au moins les étapes suivantes consistant à : viser (1004), avec la caméra disposée dans un premier emplacement (k1) de prise de vue situé dans le local, au moins deux points (P l ) d’une cible (C1, C2, C3) fixe dans le repère de base (BF); déterminer les coordonnées ( de chaque point visé à l’étape a) dans le repère (R Caméra ) lié à la caméra (30) qui est dans le premier emplacement (k1) de prise de vue ; déplacer (1005) la caméra vers un deuxième emplacement (k2) de prise de vue situé dans le local, différent du premier emplacement (k1) de prise de vue ; viser (1006), avec la caméra disposée dans le deuxième emplacement (k2) de prise de vue, les points (P l ) de la cible fixe déjà visés à l’étape a) ; déterminer les coordonnées ( de chaque point visé à l’étape d) dans le repère (R Caméra ) lié à la caméra qui est dans la deuxième position ; calculer (1008), pour chaque point visé par la caméra aux étapes a) et d), une différence ( entre les coordonnées de ce point dans le repère de base, exprimées en fonction des coordonnées déterminées aux étapes b) et e) en utilisant les première et deuxièmes matrice de passage ( ; calculer (1010, 1018) au moins une valeur globale de différence au moyen d’une fonction (F, G) ayant comme variable une différence calculée à l’étape f) ; et en ce que des variables (X 1 -X 6 ) de coefficients de la première matrice de passage ( et des variables ( de coefficients de la deuxième matrice de passage ( sont déterminées en minimisant (1012, 1018 ; 1021) la valeur globale calculée à l’étape g). Méthode selon la revendication 1, caractérisée en ce qu’une première valeur globale de différence (F) est calculée à l’étape g) et en ce que la méthode comprend des étapes postérieures à l’étape g) et consistant à : déterminer (1012) les variables (X 1 -X 6 ) des coefficients de la première matrice de passage ( en minimisant la première valeur globale calculée à l’étape g) ; mettre à jour (1014) la première matrice de passage ( avec des coefficients qui intègrent les variables (X 1 -X 6 ) déterminées à l’étape h) calculer (1016), pour chaque point (P l ) visé par la caméra aux étapes a) et d), une différence ( entre les coordonnées de ce point dans le repère de base (BF), exprimées en fonction des coordonnées déterminées aux étapes b) et e) en utilisant la première matrice de passage mise à jour à l’étape i) ; calculer (1018) une deuxième valeur globale de différence au moyen d’une fonction (G) ayant comme variable les différences calculées à l’étape i) ; déterminer (1020) les variables ( des coefficients de la deuxième matrice de passage en minimisant la deuxième valeur globale calculée à l’étape j) ; et mettre à jour (1022) la deuxième matrice de passage avec des coefficient qui intègrent les variables ( déterminées à l’étape l). Méthode selon la revendication 1, caractérisée en ce qu’elle comprend une étape postérieure à l’étape g) et consistant à h’) déterminer (1021) les variables (X 1 -X 6 ) des coefficients de la première matrice et les variables ( des coefficients de la deuxième matrice déterminés, en minimisant une valeur globale unique (F+G) calculée à l’étape g). Méthode selon l’une des revendications 2 et 3, caractérisée en ce qu’elle comprend des étapes, postérieures à l’étape h) ou à l’étape h’) et constant à : calculer (1024), pour chaque point (P l ) visé par la caméra aux étapes a) et d), une différence ( entre les coordonnées de ce point dans le repère de base (BF), exprimées en fonction des coordonnées ( déterminées aux étapes b) et e) en utilisant la première matrice de passage ( mise à jour avec les dernières variables (X 1 -X 6 ) de coefficients déterminées pour cette première matrice et la deuxième matrice de passage ( mise à jour avec les dernières variables ( des coefficients déterminées pour cette deuxième matrice. calculer une valeur globale de différence au moyen d’une fonction (F, G) ayant comme variable les différences calculées à l’étape n) ; comparer la valeur globale calculée à l’étape o) à une valeur de seuil (F0, G0); si la comparaison de l’étape p) montre que la valeur globale (F, G) calculée à l’étape o) est supérieure à la valeur de seuil (F0, G0), mettre à nouveau en œuvre les étapes f) et suivantes ; si la comparaison de l’étape p) montre que la valeur globale (F, G) calculée à l’étape o) est inférieure à la valeur de seuil (F0, G0), figer (1026) les première et deuxième matrices de passage avec les dernières variables (X 1 -X 6 , déterminées pour les utiliser dans le modèle géométrique du robot multiaxe. Méthode selon l’une des revendications précédentes, caractérisée en ce que la différence calculée à l’étape f), j) ou n) s’exprime sous la forme où est la différence entre les coordonnées du même point numéroté vu par la caméra à partir de positions k1 et k2 ; est un vecteur représentant la position du point l vu par la caméra à partir de la position k1 ; est un vecteur représentant la position du point l vu par la caméra à partir de la position k2 ; est la première matrice de passage ; ou est la deuxième matrice de passage. Méthode selon l’une des revendications précédentes, caractérisée en ce que la valeur globale de différence calculée à l’étape g), k) ou o) est la somme des carrés des valeurs des différences calculées à l’étape f), i) ou l) et s’exprime sous la forme où est la différence entre les coordonnées du même point numéroté l vu par la caméra à partir d’emplacements k1 et k2 ; l est le numéro d’ordre d’un des points visés par la caméra aux étapes a) et d), compris entre 1 et L k1 est le numéro d’ordre du premier emplacement, compris entre 1 et K ; k2 est le numéro d’ordre du deuxième emplacement, compris entre 1 et K ; K est le nombre d’emplacements que peut prendre la caméra pour viser les points de la cible. Méthode selon l’une des revendications précédentes, caractérisée en ce que la première matrice de passage s’exprime sous la forme où les grandeurs Xi, pour i entier naturel entre 1 et 6, sont les variables déterminées en minimisant la valeur globale calculée à l’étape g) et correspondent pour i entier naturel entre 1 et 3, à des distances entre un centre du repère (R Caméra ) lié à la caméra (30) et le centre du repère (R PG ) lié au poignet (21); pour i entier naturel entre 4 et 6, à des angles représentatifs d’une orientation du repère lié à la caméra dans le repère lié au poignet. Méthode selon l’une des revendications précédentes, caractérisée en ce que la deuxième matrice s’exprime sous la forme où les grandeurs sont des angles mesurés sur l’axe A i du robot multiaxe (20), pour le k-ième point visé par la caméra (30) au cours des étapes a) et d) ; les grandeurs sont les paramètres théoriques de Denavit-Hartenberg modifiés pour l’axe A i du robot multiaxe ; et les grandeurs sont des variables pour les coefficients de la deuxième matrice de passage, déterminées en minimisant la valeur globale calculée à l’étape g). Méthode selon l’une des revendications précédentes, caractérisée en ce que, la minimisation de la valeur globale de différence (F, G) est effectuée par la méthode des moindres carrés avec résolution d’un système de six équations à six inconnues lors de la détermination des variables (X 1 -X 6 ) de coefficient de la première matrice de passage ( et résolution d’un système de vingt-quatre équations à vingt-quatre inconnues lors de la détermination des variables ( de coefficient de la deuxième matrice de passage ( , dans le cas de la méthode de la revendication 2, ou résolution d’un système de trente équations à trente inconnues lors de la détermination des variables (X 1 -X 6 , ) de coefficient des première et deuxième matrices de passage, dans le cas de la méthode de la revendication 3. Application de la méthode selon l’une des revendications précédentes pour l’ajustement du modèle géométrique d’un robot (20) d’application de produit de revêtement équipé d’une tête d’impression (10) ou d’un projecteur de produit de revêtement. Robot (20) équipé d’un outil (10), d’une caméra (30) et d’une unité électronique de contrôle (24, 40) programmée pour mettre en œuvre automatiquement une méthode selon l’une des revendications 1 à 9.