Brevet FR3007832A1 — publié le 20150102

La présente invention se rapporte à un système de navigation et de pointage dans lequel une centrale de navigation est couplée de manière solidaire à un dispositif de pointage. La centrale de navigation fournit des mesures précises de position et d'orientation du système par rapport à un repère global et le dispositif de pointage fournit une mesure d'un paramètre, par exemple de distance, dans une direction de pointage. Plus particulièrement, l'invention se rapporte à un système qui couple sur une même plateforme une centrale de navigation inertielle à un lidar (acronyme pour Light Detection And Ranging) à balayage en deux dimensions (ou lidar 2D) ou en trois dimensions (lidar 3D). Plus particulièrement, l'invention concerne une méthode de calibration de l'erreur angulaire de pointage dans un système de navigation et de pointage, par exemple un lidar terrestre à 2D ou à 3D couplé à une centrale de navigation inertielle. Dans sa conception la plus simple, un lidar est un appareil optique comportant une source laser, qui émet un faisceau laser pulsé dans un axe de visée, appelé aussi direction de pointage, et un système de détection qui mesure le temps de propagation aller et retour du faisceau laser pulsé entre la source laser et la cible, le faisceau laser étant réfléchi ou rétrodiffusé par la cible dans la direction de pointage et en sens inverse du faisceau laser incident. Un lidar permet ainsi de mesurer la distance à une cible. Il existe des systèmes de navigation et de pointage intégrant à la fois un Lidar et une centrale de navigation, qui indique la position et l'orientation géographique du système, pour permettre de localiser une cible dans l'espace. On connaît en particulier des systèmes de navigation et de pointage dans lesquels un lidar est combiné à une centrale navigation inertielle (SNI ou INS pour lnertial Navigation System). Dans ces systèmes de navigation et de pointage, il existe différents types de lidar selon les moyens utilisés pour orienter la direction de pointage du faisceau laser dans l'espace : - un lidar 2D combine un lidar avec un système opto-mécanique à un degré de liberté, par exemple de type miroir rotatif selon un axe orthogonal à l'axe de visée du laser, adapté pour balayer l'orientation du faisceau laser dans un plan de l'espace (i.e. en 2D) et acquérir en retour des mesures suivant une ligne (par exemple par projection sur un plan). Dans le présent document, chaque ligne ainsi acquise s'appelle une ligne de scan ou ligne de balayage ; - un Lidar 3D combine un Lidar avec un système opto-mécanique à deux degrés de liberté, adapté pour balayer l'orientation du faisceau laser dans toutes les directions de l'espace (i.e. en 2D) et acquérir en retour des mesures provenant de n'importe quel point de l'espace. Par exemple un Lidar 3D peut être constitué d'un Lidar 2D monté sur une table tournante permettant des mouvements rotatifs selon un axe orthogonal à celui du miroir de balayage à 2D. On appelle « scanner laser » un système intégrant un lidar-2D ou lidar-3D avec un système de positionnement et un système d'orientation. La figure 1 représente un exemple de système de navigation et de pointage 10 comprenant une centrale de navigation inertielle 11 et un lidar 12 reliés mécaniquement par une pièce rigide 13 de manière à ce que le SNI et le lidar soient solidaires. Le lidar 12 comporte une source laser et un système de détection d'un faisceau rétrodiffusé. Le lidar 12 comporte aussi un miroir plan 15 monté mobile en rotation autour d'un axe de rotation 6. Par exemple, la normale au plan du miroir 15 est inclinée à 45 degrés par rapport à la direction d'émission du faisceau laser qui est de préférence coaxiale à l'axe de rotation 16. Par rotation 26 du miroir 15 autour de l'axe de rotation 16, la direction de pointage 17 du faisceau laser balaye un plan 27, ici appelé plan de balayage. On obtient ainsi un lidar à deux dimensions (lidar-2D) qui permet de mesurer la distance d'une cible situé dans le plan de balayage 27. Connaissant la direction de pointage 17 dans le repère de l'INS, il est ensuite possible de fournir les indications de position et d'orientation de la cible par rapport à un repère global. Selon les dispositifs, la rotation du miroir 15 autour de l'axe de rotation 16 peut être soit continue, et toujours dans le même sens, soit oscillante le miroir 15 effectuant des allers-retours entre deux butées. Le système de navigation et de pointage 10 est par exemple monté sur une plateforme 14 mobile en rotation autour d'un axe transverse à l'axe de rotation 16 du miroir 15. Une rotation de la plateforme 14 entraîne la rotation de l'ensemble du système 10 et donc la rotation du plan de balayage 27. On obtient ainsi un lidar 3D qui permet de mesurer la distance d'une cible dans une direction de pointage 17 quelconque dans l'espace en trois dimensions (lidar-3D). Afin de pouvoir géo-référencer les points LiDAR, il est nécessaire de coupler le lidar avec d'autres systèmes permettant de connaitre la position et l'orientation du LiDAR. La position et l'orientation sont en général fournies par un système de positionnement (noté SP) tel qu'un système combinant un système de localisation GPS et une centrale navigation inertielle {GPS+SNI} ou un système de localisation GPS et 30 une centrale de cap et de verticale {GPS+AHRS pour Attitude and Heading Reference System} ou autre. Dans la suite du document, on entend par système de positionnement tout système adapté pour fournir sa position et son orientation dans un repère global sans faire de distinction entre les différents sous-systèmes possibles. Une centrale AHRS (Attitude Heading Reference Système) ou une centrale SNI 35 (Système de Navigation Inertiel) sont des systèmes fournissant leur orientation dans un repère global tel que le repère ECEF. Reste alors à connaitre l'orientation relative du LiDAR par rapport au système de positionnement et d'orientation qui lui est mécaniquement attaché. Dans la suite du présent document, on entend par « système LiDAR » un 40 système qui combine un instrument de visée (le LiDAR), un système de positionnement et un système d'orientation. Un tel système lidar fournit ainsi des mesures de position, d'orientation et de visée. A titre d'exemple illustratif et non limitatif, dans la suite du présent document, un système fournissant l'orientation est illustré par un SNI. Un repère 21 lié au SNI, appelé ici repère SNI est défini par un centre OsNi et trois axes XSNI, YSNI et ZsNi. Par convention, XSNI pointe vers l'avant, YsNi pointe vers la droite et ZsNi pointe vers le bas. La centrale de navigation 11 permet à chaque instant de déterminer la position et l'orientation du repère SNI par rapport à un repère géographique, ou repère global, par exemple le repère ecef (Earth centered, Earth fixed) défini par un centre géographique OEcEF et trois axes géographiques XECEF, YECEF et ZECEF. Le centre du repère ecef est le centre de la terre. Les axes XECEF et YECEF sont dans le plan équatorial, XECEF pointe dans la direction du méridien de Greenwich et YECEF pointe vers 90deg. Est. Enfin ZECEF pointe vers le pôle. Le repère ecef est un repère global, son centre est fixe. Un repère 22 propre au dispositif de pointage, appelé ici repère lidar est défini par un centre Olidar et trois axes Xlidar, Ylidar et Ziidar. Le LiDAR considéré ici est un LiDAR 2D. Par convention, l'axe Y - lidar est l'axe de rotation 16 du miroir 15 et pointe vers l'interface du LiDAR que l'on considère comme étant l'avant du LiDAR. Ainsi, le LiDAR balaie le plan (Xlidar, Ziidar ) où Xlidar pointe vers la droite du LiDAR et Ziidar pointe vers le haut. Le plan (Xlidar, Ziidar) est aussi repéré plan de balayage 27 de la direction de pointage 17 du faisceau laser. La figure 2 représente schématiquement le repère SNI 21 et le repère lidar 22 du système de la figure 1, par rapport au repère géographique 20. La distance entre le centre OsNi du repère SNI et le centre Olidar du repère lidar est communément appelée bras de levier. Par construction, dans le présent document, l'axe XsNi du repère 21 du SNI est pratiquement parallèle à l'axe Y - lidar du repère 22 du lidar, l'axe YsNi du repère 21 du SNI est pratiquement parallèle à l'axe Xlidar du repère 22 du lidar et l'axe ZsNi du repère 21 du SNI est pratiquement parallèle et opposé à l'axe Ziidar du repère 22 du lidar. Dans le présent document, on entend par pratiquement parallèles, alignés à ± 5 deg., voire à mieux que 1 degré. On définit une matrice de rotation comprenant trois angles d'Euler pour passer du repère SNI au repère lidar. Par construction, on connaît en général ces angles de rotation à mieux que quelques degrés, voire à quelques dixièmes de degrés près. Toutefois, étant donné que le lidar mesure des distances allant de quelques mètres à plusieurs centaines de mètres, une erreur sur l'orientation angulaire de la direction de pointage peut se traduire par une erreur de positionnement, qui n'est pas une erreur de biais constante, mais une erreur variable selon la distance de la cible. Il est donc nécessaire de définir une méthode de calibration pour déterminer les angles de désalignement entre le repère SNI et le repère lidar. On appelle angles de désalignement (boresight error), les angles correspondant aux erreurs angulaires entre le repère estimé du SNI et le repère 21 attaché au SNI. La figure 3 représente schématiquement le changement de repère pour passer du repère 21 propre au SNI à un repère estimé 31. La figure 3 illustre les trois angles d'Euler de désalignement (d), 0, Y) que l'on cherche à calibrer. L'angle de désalignement représente un angle de rotation autour de l'axe ZSNI pour passer du repère (XsNi, YsNi) à un repère (X'sNi, Y'sNi). L'angle de désalignement 0 représente un angle de rotation autour de l'axe Y'sNi pour passer du repère (X'sNi, ZsNi) à un repère (X"sNi, l'sNi), l'axe Y'sNi étant confondu avec l'axe Y"SNI. L'angle de désalignement c1) représente un angle de rotation autour de l'axe X"sNi pour passer du repère (Y"sNi, Z"sNi) au repère estimé 31. Lors de l'installation du système lidar, on connaît l'orientation relative des deux dispositifs à quelques degrés près. Cette orientation relative est exprimée par la matrice de changement de repère C LiSiç/EI,AR L'objet de la calibration des angles de désalignement est de déterminer précisément les trois angles d'Euler de désalignement (d), 0, Y) représentés sur la figure 3. Ce problème de calibration des angles de désalignement a déjà été abordé en particulier dans des applications de lidar aérien (airborne lidar) dans lesquelles un système lidar est embarqué sur un aéronef de manière à scanner ou balayer la surface du sol et des toits des immeubles pendant un déplacement de l'aéronef. Dans ces applications de lidar aérien, l'erreur sur les angles de désalignement est particulièrement critique. En effet, les lidars aériens balayent des distances de plusieurs centaines de mètres. Or plus la distance de pointage est grande, plus l'erreur sur la position des points mesurés par le système lidar est importante, donc facilement observable. La publication K. Kris Morin, Naser EI-Sheimy, 2002, "Post-mission adjustment methods of airborne laser scanning data" In FIG XXII International Congress, Washington, DC USA présente l'impact des erreurs d'angle de désalignement sur le positionnement des points mesurés par un lidar aérien puis expose différentes méthodes de calibration. On distingue deux catégories de méthodes de calibration : d'une part les méthodes dirigées par les données (Data driven), d'autre part les méthodes dirigées par le système (System driven). Une première méthode de calibration dirigée par les données consiste à ajuster manuellement différents acquisitions de lignes ou d'images par balayage du faisceau laser jusqu'à ce que le résultat soit satisfaisant visuellement. Cette méthode n'est absolument pas rigoureuse mathématiquement et ne fournit pas de mesure statistique sur la qualité de la calibration. Une autre méthode de calibration dirigée par les données est décrite dans la publication Chao Gao, Spletzer, J.R., « On-line calibration of multiple lidars on a mobile vehicle platform », Robotics and Automation (ICRA), 2010 IEEE International Conference, vol., no., pp.279, 284, 3-7 May 2010. Selon cette publication, des points remarquables, tels que des lignes ou des contours, sont identifiés dans différents images produites par balayage (ou scan) d'une même zone. Ces points remarquables sont ensuite utilisés afin de calculer la matrice de rotation R et le vecteur translation T minimisant l'écart entre ces points ,AR . Après correction, on obtient les points x7,e,;' : x7,e,;' = R.k7,e,;' +T Les méthodes dirigées par les données ne reposent sur aucun modèle et ne nécessitent pas de récupérer les données position du GPS du véhicule sur lequel est embarqué le système ou d'attitude du SNI mais seulement les points lidar. D'un premier abord simples, ces méthodes ne le sont en réalité pas tant que ça. En effet, la sélection de points remarquables n'est pas triviale car même avec une acquisition à haute densité de points, il est difficile de sélectionner précisément les mêmes points sur plusieurs scans, d'autant plus que ces points sont bruités. Cette sélection de points n'est pas automatique, elle prend beaucoup de temps et le résultat de la calibration dépend grandement de l'utilisateur. Par ailleurs, les méthodes dirigées par les données ne se basant pas sur un modèle d'erreur rigoureux, elles ne fournissent pas d'outils de mesure statistiques permettant de juger de la qualité de l'ajustement. Les méthodes de calibration dirigées par le système reposent sur un modèle d'erreur, tel que présenté par exemple dans les publications Skaloud, J., Litchi, D., 2006, « Rigorous approach to boresight self-calibration in airborne laser scanning », ISPRS Journal of Photogrammetry & remote Sensing 61, 47-59 et Skaloud, J., Schaer, P., « Towards automated lidar boresight self-calibration », in Proc. 5th International Symposium on Mobile Mapping Technology, Padua, Italy, 6 p, May 2007. Ce modèle d'erreur permet d'identifier des erreurs systématiques, dont en particulier l'erreur de désalignement entre le lidar et la centrale de navigation. Ainsi, une même opération de sélection de points remarquables que précédemment est effectuée mais l'estimation des angles de désalignement passe par l'utilisation d'un modèle d'erreur. On utilise ensuite un modèle d'erreur du type : xngAR = xef + cefSN (b+ (Cral CAR) xuee)I où ecef (earth centered earth fixed) désigne un répère géocentré et géostationnaire, C AB représente de manière générale une matrice de passage d'un repère A à un repère B ; X'8, : représente de manière générale une position d'un dispositif a dans le repère b : représente le bras de levier entre le lidar et le système de positionnement. Le modèle d'erreur peut aussi prendre en compte d'autres paramètres d'erreurs tels que les erreurs de bras de levier, le biais sur les angles du SNI, le biais sur la mesure de distance par le lidar... Toutefois, ces méthodes de calibration nécessitent l'intervention d'un utilisateur pour sélectionner des points remarquables dans les différents scans. La publication Skaloud, J., Schaer, P., « Towards automated lidar boresight selfcalibration », in Proc. 5th International Symposium on Mobile Mapping Technology, Padua, Italy, 6 p, May 2007 propose d'automatiser la procédure de calibration d'un système lidar aérien, en sélectionnant un ensemble de points appartenant à une même surface, par exemple de toiture, au lieu de comparer des points remarquables deux à deux. Cette approche permet l'automatisation de l'algorithme de calibration. Il est en effet plus facile de mettre en place des traitements de détection/extraction automatique de plans et d'identification de plans communs entre différents scans. Cette méthode donne des résultats satisfaisants en lidar aérien car la distance aux cibles est importante et rend donc l'erreur facilement observable. Dans les systèmes dits de lidar terrestre, un système lidar est embarqué sur un véhicule terrestre, par exemple sur le toit d'une voiture, et le faisceau laser balaye les façades des immeubles au cours de déplacements du véhicule le long d'une rue par 20 exemple. Les méthodes classiques de calibration de système de navigation et pointage en lidar terrestre s'inspirent très largement l'approche de calibration lidar aérien qui semble donner de bons résultats. Ainsi, la publication Rieger, P., Studnicka, N., Pfennigbauer, M., Zach, G., 25 Boresight alignment method for mobile laser scanning systems, Journal of Applied Geodesy. Volume 4, Issue 1, Pages 13-21 décrit la calibration d'un système lidar terrestre en scannant un plan de manière dynamique, lors d'un déplacement du système dans un sens puis dans l'autre et en faisant varier l'angle de cap du lidar suivant quatre positions. Cette calibration nécessite au minimum 8 scans au total, 30 chaque orientation devant être répétée au moins deux fois. La procédure se déroulant en dynamique, les points doivent être exprimés dans un repère global (ex : ECEF). Une fois les points exprimés dans le repère global, un algorithme détecte les plans communs formés par les points lidar entre les différents scans et estime les angles de désalignement en ajustant ces plans des différents scans 35 les uns par rapport aux autres. Toutefois, les performances de cette calibration de lidar terrestre sont moindres qu'en lidar aérien. Premièrement, les distances de pointage sont beaucoup plus faibles (de l'ordre de 10m) rendant l'erreur sur la position des points lidar beaucoup moins facilement observable. De plus, il n'est pas facile de faire varier la distance de pointage 40 de manière suffisamment importante pour mettre en évidence l'erreur de désalignement. Ensuite, la procédure de calibration étant effectuée de manière dynamique, il est nécessaire d'utiliser un système de positionnement fournissant la position du véhicule en temps réel. Or, au cours des déplacements terrestres, le signal du GPS peut être masqué ou affecté par des trajets multiples, ce qui est moins fréquent dans le cas aérien. D'autre part, en lidar terrestre, la diversité des orientations des plans est moindre : on trouve surtout des plans horizontaux et verticaux alors qu'en lidar aérien les toits des immeubles offrent généralement une gamme d'orientations très variées permettant de diminuer la corrélation entre les différents angles de désalignement. Enfin, cette méthode se révèle sensible au biais sur les angles donnés lo par la centrale de navigation, si bien que les valeurs des angles de désalignement découlant de cette calibration sont en général biaisées. En résumé, les méthodes de calibration des angles de désalignement sont sensibles à différentes sources d'erreurs : les erreurs humaines, lors du choix des points remarquables, les erreurs de GPS, lors du déplacement du système au cours de 15 la calibration, et les erreurs de biais du SNI. De plus, ces méthodes de calibration sont chronophages. Il existe donc un besoin d'une méthode de calibration des angles de désalignement à la fois facile et rapide à mettre en oeuvre et de préférence automatique. Un autre but de l'invention est de fournir une méthode de calibration qui 20 soit immune aux erreurs GPS. Un autre but de l'invention est de fournir une méthode de calibration qui soit robuste aux bras de levier. Encore un autre but de l'invention est de fournir une méthode de calibration qui soit insensible aux biais du SNI afin de déterminer des angles de désalignement non biaisés. La présente invention a pour but de remédier aux inconvénients des méthodes et 25 dispositifs antérieurs. L'invention concerne en particulier un procédé de calibration d'un système de navigation et de pointage comprenant une centrale de navigation reliée solidairement à un dispositif de pointage, la centrale de navigation étant configurée pour mesurer la position et l'orientation d'un repère de ladite centrale de navigation, le dispositif de 30 pointage comprenant des moyens d'émission d'un faisceau dans une direction de pointage et des moyens de balayage de la direction de pointage, le dispositif de pointage étant configuré pour mesurer un signal de pointage, de préférence par réflexion ou rétro-diffusion dudit faisceau sur une cible, en fonction de la direction de pointage dans un repère lié au dispositif de pointage, le procédé de calibration 35 comprenant les étapes suivantes de calibration des angles entre le repère de ladite centrale de navigation et le repère lié au dispositif de pointage : a. Estimation d'un repère approché de navigation, le repère approché de navigation étant incliné par rapport au repère de navigation avec des angles d'Euler de désalignement (d), 0, Y) ; 40 b. Sélection d'une première surface plane, de préférence verticale ; c. Détermination d'un vecteur normal à ladite première surface plane dans le repère de navigation ; d. Orientation du système de navigation et de pointage dans une pluralité d'orientations angulaires distinctes par rapport à la première surface plane ladite pluralité d'orientations angulaires étant déterminées dans le repère de la centrale de navigation ; e. Balayage de la direction de pointage de manière à acquérir, au moyen du dispositif de pointage, pour chacune des orientations de ladite pluralité d'orientations par rapport à la première surface plane, une série de mesures de pointage suivant une ligne de balayage sur la première surface plane ; f. Traitement de chacune desdites séries de mesures de pointage suivant une ligne de balayage sur la première surface plane de manière à déterminer, dans le repère lié au dispositif de pointage, un vecteur directeur approché de chaque ligne de balayage sur la première surface plane ; g. Traitement d'une pluralité de vecteurs directeurs approchés d'une pluralité de lignes de balayage sur la première surface plane pour en déduire un vecteur normal approché de la première surface plane ; h. Calcul d'au moins un angle de désalignement (d), 0, Y) entre le vecteur normal approché déterminé à l'étape précédente dans le repère de navigation approché et le vecteur normal déterminé à l'étape c) dans le repère de la centrale de navigation. Selon un mode de réalisation préféré, le procédé de calibration comporte en outre les étapes suivantes : i. Sélection d'une deuxième surface plane, de préférence horizontale ; j. Détermination d'un vecteur normal à ladite deuxième surface plane dans le repère de navigation ; k. Orientation du système de navigation et de pointage dans une pluralité d'orientations angulaires distinctes par rapport à la deuxième surface plane, ladite pluralité d'orientations angulaires étant déterminées dans le repère de la centrale de navigation ; I. Balayage de la direction de pointage de manière à acquérir, au moyen du dispositif de pointage, pour chacune des orientations de ladite pluralité d'orientations par rapport à la deuxième surface plane, une série de mesures de pointage suivant une ligne de balayage sur la deuxième surface plane ; m. Traitement de chacune desdites séries de mesures de pointage suivant une ligne de balayage sur la deuxième surface plane de manière à déterminer, dans le repère lié au dispositif de pointage, un vecteur directeur approché de chaque ligne de balayage sur la deuxième surface plane ; n. Traitement d'une pluralité de vecteurs directeurs approchés d'une pluralité de lignes de balayage sur la deuxième surface plane pour en déduire un vecteur normal 40 approché de la deuxième surface plane ; o. Calcul d'au moins un autre angle de désalignement (c1), 0, Y) entre le vecteur normal approché déterminé à l'étape précédente dans le repère de navigation approché et le vecteur normal déterminé à l'étape j) dans le repère de la centrale de navigation. Selon un aspect particulier du procédé de l'invention, ladite pluralité d'orientations angulaires comporte les quatre orientations angulaires du système de navigation et de pointage par rapport à la première surface plane suivantes, où r désigne l'angle de roulis, t l'angle de tangage et c l'angle de cap du système de navigation et de pointage par rapport à un repère local tel que les axes X et Y de ce repère local soient dans le plan horizontal, et l'axe X soit parallèle au plan la première surface plane : - Orientation 1 : (r, t, c) = (0, 0, 0± 5 deg.) - Orientation 2 : (r, t, c) = (0, 0, 180± 5 deg.) - Orientation 3 : (r, t, c) = (0, 0, 55 ± 25 deg.) - Orientation 4 : (r, t, c) = (0, 0, -55 ± 25 deg.) et/ou les deux orientations angulaires suivantes du système de navigation et de pointage par rapport à la deuxième surface plane : Orientation 5 : (r, t, c) = (0, 55 ± 25 deg., 0) Orientation 6 : (r, t, c) = (0, -55 ± 25 deg., 0) Selon un aspect particulier du procédé de l'invention : - l'étape e) et/ou l'étape I) comporte une pluralité de balayage d'une même ligne de balayage et en ce que ladite série de mesures de pointage suivant une ligne de balayage est obtenue par un calcul de moyenne de la pluralité de balayage d'une même ligne de scan ; - l'étape g) et/ou l'étape n) de détermination d'un vecteur normal approché comporte une étape de détermination d'un vecteur propre associé à la plus petite valeur propre d'une analyse en composantes principales des séries de mesures de pointage suivant ladite pluralité de lignes de balayage ou une étape de calcul d'une moyenne d'une pluralité de vecteurs normaux mesurés par produit vectoriel des vecteurs directeurs approchés de deux lignes de scan différentes déterminés à l'étape f) et/ou à l'étape m). Selon un mode de réalisation particulier, l'étape h) et/ou l'étape o) de calcul des angles de désalignement (d), 0, Y) comporte une étape de maximisation d'un produit scalaire pour en déduire une estimation d'un angle de désalignement, ladite étape de maximisation comportant une étape de traitement par une méthode de régression des moindres carrés non linéaires ou par une méthode de régression des moindres carrés linéaires ou une méthode de Quasi-Newton. Selon un mode de réalisation particulier et avantageux, le procédé de calibration comporte une étape de calcul d'un indicateur de précision des angles de désalignement (d), 0, Y). L'invention concerne aussi un système de navigation et de pointage comprenant une centrale de navigation reliée solidairement à un dispositif de pointage, la centrale de navigation étant configurée pour mesurer la position et l'orientation d'un repère de ladite centrale de navigation, le dispositif de pointage comprenant des moyens d'émission d'un faisceau dans une direction de pointage et des moyens de balayage de la direction de pointage, le dispositif de pointage étant configuré pour mesurer un signal de pointage, de préférence par réflexion ou rétrodiffusion dudit faisceau sur une cible, en fonction de la direction de pointage dans un repère lié au dispositif de pointage, le système de navigation comprenant des moyens de traitement des mesures de position et d'orientation de la centrale de navigation et des moyens de traitement du signal de pointage, et des moyens de calibration adaptés pour implémenter une calibration angulaire entre le repère de ladite centrale de navigation et le repère de pointage suivant un procédé de calibration selon l'un des modes de réalisation décrits. Selon des aspects particuliers et avantageux du système de navigation et de pointage : - ladite centrale de navigation comporte un système d'orientation parmi une centrale de navigation inertielle ou une centrale d'attitude de type AHRS et/ou un système de positionnement de type GPS ; - le dispositif de pointage comporte un lidar à balayage suivant un plan en deux dimensions, un lidar à balayage en trois dimensions ou un sonar ; De façon avantageuse, le système de navigation et de pointage comporte en outre des moyens d'orientations angulaire dudit système de navigation et de pointage, adaptés pour orienter ledit système de navigation et de pointage dans une pluralité d'orientations angulaires distinctes par rapport à ladite première surface plane et/ou, respectivement, à ladite deuxième surface plane. L'invention trouvera une application particulièrement avantageuse dans les systèmes de navigation et de pointage comprenant une centrale inertielle et un lidar. L'invention s'applique aussi à d'autres systèmes intégrant à la fois un système de positionnement, un système d'orientation et un système de visée, et en particulier un système de visée tel qu'un sonar basé sur l'émission et la réception de signaux acoustiques. La présente invention concerne également les caractéristiques qui ressortiront au cours de la description qui va suivre et qui devront être considérées isolément ou selon toutes leurs combinaisons techniquement possibles. Cette description donnée à titre d'exemple non limitatif fera mieux comprendre comment l'invention peut être réalisée en référence aux dessins annexés sur lesquels : - la figure 1 représente un système lidar à balayage en 2D ; - la figure 2 représente respectivement les repères du SNI et du lidar dans un système lidar ; - la figure 3 représente schématiquement les trois angles d'Euler d'erreur d'alignement que l'on cherche à calibrer pour passer du repère SNI (le vrai repère du dispositif SNI) à un repère sir (repère approximé) ; - la figure 4 illustre l'angle de roulis r d'un système lidar dans une première orientation de calibration face à un mur de référence ; - la figure 5 illustre l'angle de tangage t d'un système lidar dans une deuxième orientation de calibration face à un mur de référence ; - la figure 6 illustre l'angle de cap c d'un système lidar dans une autre orientation face à un mur de référence ; - la figure 7 illustre l'orientation d'un système lidar dans une orientation de calibration par rapport à un sol de référence ; - la figure 8 illustre l'orientation d'un système lidar dans une autre orientation de calibration par rapport à un sol de référence ; - la figure 9 est une vue de dessus pour illustrer l'orientation d'un système lidar pour la calibration par rapport à un mur de référence ; - la figure 10 est une vue de côté pour illustrer l'orientation d'un système lidar pour la calibration par rapport à un sol de référence. Dans la description du mode de réalisation qui va suivre, on considère la calibration d'un système lidar comprenant un lidar 2D. Néanmoins, l'invention n'est nullement limitée à ce mode de réalisation. L'homme du métier transposera aisément la méthode décrite à un système lidar à 3D, en limitant le balayage à 2D. On considère un système lidar 10 placé à une position fixe dans un repère géographique. Le système lidar 10 est monté solidairement comme illustré par exemple sur la figure 2. Avant de débuter l'acquisition des données, la centrale de navigation 11 doit être alignée selon les instructions du constructeur, afin de déterminer précisément l'orientation angulaire du repère 21 (SNI) de la centrale de navigation 11 par rapport à un repère géographique de référence 20, par exemple le repère ecef (earth centered earth fixed). Le procédé de calibration comprend l'orientation angulaire du système lidar dans une pluralité d'orientations par rapport à une surface plane de référence, le système lidar pouvant rester statique, c'est à dire à une même position de l'espace par rapport au repère géographique, car la position du système n'a aucun impact sur la calibration, seule l'orientation angulaire du système lidar par rapport au mur compte. On précise bien qu'il s'agit ici de l'orientation du système lidar, le SNI étant solidaire du lidar. Les différentes orientations du système lidar sont décrites de manière détaillée ci-dessous. Le procédé de calibration comprend l'acquisition de mesures lidar suivant une ligne de scan ou ligne de balayage, pour chacune des orientations du système lidar par rapport à la surface plane de référence. Le procédé de calibration comprend ensuite le traitement des différents scans acquis pour les différentes des orientations du système lidar par rapport à la surface plane de référence, pour en extraire une mesure d'un ou de plusieurs angles d'Euler de désalignement. On définit l'orientation du système lidar 10 par rapport à une première surface plane de référence à l'aide de trois angles d'Euler : (r, t, c). Ces angles sont respectivement les angles de roulis, tangage et cap du SNI définis par rapport à un repère local au site de calibration. Plus précisément, à titre d'exemple illustratif, on choisit comme première surface lo plane de référence un mur 8 supposé plan et vertical dans le repère géographique. On a représenté sur les figures 4-8 le repère du lidar (Xlidar, Ylidar, Zlidar). Dans la suite de l'exposé, l'axe de rotation de balayage du faisceau laser est supposé être l'axe Ylidar. On a aussi représenté sur les figures 4-6 un repère (Xmur, Ymur, Zmur) lié au mur. Le repère mur est propre à chaque scan car il dépend de l'orientation du scan. L'axe Xmur 15 correspond à la normale du mur, les axes Ymur et Zmur sont donc dans le plan du mur avec Zmur colinéaire à la direction du scan (voir figure 4-6). Le repère mur s'applique aux scans de n'importe quelle surface plane, que ce soit un mur vertical (Figures 4-6 et 8) ou un sol horizontal (figure 9). On a représenté sur les figures 4-6 un repère local au site de calibration. Ce 20 repère local est défini pour des raisons pratiques afin d'expliquer clairement les angles utilisés pour la procédure de calibration. Le repère local est seulement utilisé pour exprimer des orientations, l'emplacement de son origine n'a pas d'importance, seule la direction de ses trois axes importe. Ce repère local correspond à un repère de navigation à une rotation près autour de Z. Ainsi )(local et Ylocal sont dans le plan tangent 25 local avec Ylocal pointant vers le mur et )(local vers la gauche du mur. Zocal est confondu avec la verticale locale et pointe vers le bas. La figure 4 illustre schématiquement l'angle de roulis r du système lidar par rapport à un mur 8 plan et vertical. La direction de pointage 17 du faisceau laser du lidar est représentée schématiquement par plusieurs lignes trait-tiret correspondant à 30 différentes orientations angulaires du miroir 5 autour de l'axe Ylidar pour produire un balayage du faisceau laser. Par projection sur le mur 8, le faisceau laser suit une ligne de scan ou ligne de balayage 41. Lorsque le lidar présente une orientation telle qu'illustrée sur la figure 4, avec uniquement un angle de roulis r non nul par rapport à une verticale locale (par exemple, une droite perpendiculaire au sol), les angles de 35 tangage t et de cap c étant nuls, la ligne de scan 41 est une ligne verticale qui est située dans le plan du mur 8. Dans ce cas, l'angle de cap c étant nul, le plan d'incidence du faisceau laser, formé par les directions de pointage 17, est normal au plan du mur 8. La figure 5 illustre schématiquement l'angle de tangage t du système lidar par rapport au mur 8 plan et vertical. De manière analogue à la figure 4, la direction de 40 pointage 17 du faisceau laser du lidar est représentée schématiquement par plusieurs lignes trait-tiret correspondant à différentes orientations angulaires du miroir 5 autour de l'axe Ylidar pour produire un balayage du faisceau laser. Par projection sur le mur 8, le faisceau laser suit une ligne de scan ou ligne de balayage 42. Lorsque le lidar présente une orientation telle qu'illustrée sur la figure 5, avec uniquement un angle de tangage t non nul, les angles de roulis r et de cap c étant nuls, la ligne de scan 42 est une ligne qui est située dans le plan vertical du mur 8 mais qui est inclinée de l'angle t par rapport à la verticale. Dans ce cas, l'angle de cap c étant nul, le plan d'incidence du faisceau laser, formé par les directions de pointage 17, est normal au plan du mur vertical 8. La figure 6 illustre schématiquement l'angle de cap c d'un système lidar par lo rapport au mur 8 plan et vertical. De manière analogue aux figures 4-5, la direction de pointage 17 du faisceau laser du lidar est représentée schématiquement par une ligne trait-tiret correspondant à différentes orientations angulaires du miroir 5 autour de l'axe Ylidar pour produire un balayage du faisceau laser. Par projection sur le mur 8, le faisceau laser suit une ligne de scan ou ligne de balayage 43. Lorsque le lidar présente 15 une orientation telle qu'illustrée sur la figure 6, avec uniquement un angle de cap c non nul, les angles de roulis r et de tangage t étant nuls, la ligne de scan 43 est une ligne verticale qui est située dans le plan du mur 8. Dans ce cas, le plan d'incidence du faisceau laser, formé par les directions de pointage 17, est incliné de l'angle c par rapport à la normale N8 au plan du mur vertical 8. 20 La figure 7 illustre schématiquement l'angle de tangage t du système lidar lors d'un scan sur un plan horizontal 9. La direction de pointage 17 du faisceau laser du lidar est représentée schématiquement par une ligne trait-tiret correspondant à différentes orientations angulaires du miroir 5 autour de l'axe Ylidar pour produire un balayage du faisceau laser. Par projection sur le sol 9, le faisceau laser suit une ligne de scan ou 25 ligne de balayage 45. Lorsque le lidar présente une orientation telle qu'illustrée sur la figure 7, avec uniquement un angle de tangage t non nul, les angles de roulis r et de cap c étant nuls, la ligne de scan 45 est une ligne qui est située dans le plan horizontal du sol 9 et qui est donc horizontale. Dans ce cas, le plan d'incidence du faisceau laser, formé par les directions de pointage 17, est incliné de l'angle (n/2-t) par rapport à la 30 normale au plan du sol horizontal 9. De manière analogue, la figure 8 illustre schématiquement un autre angle de tangage t' du système lidar lors d'un scan sur le plan horizontal 9. Par projection sur le sol 9, le faisceau laser suit une ligne de scan ou ligne de balayage 46. Lorsque le lidar présente une orientation telle qu'illustrée sur la figure 7, avec uniquement un angle de 35 tangage t' non nul, les angles de roulis r et de cap c étant nuls, la ligne de scan 46 est une ligne qui est située dans le plan horizontal du sol 9 et qui est donc horizontale. Dans ce cas, le plan d'incidence du faisceau laser, formé par les directions de pointage 17, est incliné de l'angle (7/2-f) par rapport à la normale au plan du sol horizontal. Le procédé de l'invention exploite les différentes orientations angulaires du 40 système lidar représentées sur les figures 4 à 8, la position du système Lidar étant par ailleurs statique, pour en déduire une mesure précise des angles entre le repère 22 du lidar et le repère 21 du SNI ou plus précisément une mesure précise des angles de désalignement entre le repère estimé 31 du SNI et le repère réel 21 du SNI. Dans un exemple simple de réalisation, le procédé comprend les étapes suivantes : - sélectionner une première surface plane de préférence verticale, par exemple un mur de référence ; - acquérir des mesures au moyen du système lidar en orientant la direction de pointage vers le mur de référence selon une série d'orientations du système lidar prédéterminée ; - appliquer un algorithme de traitement des mesures acquises à l'étape précédente pour en extraire des valeurs de deux angles de désalignement. Avantageusement, le procédé est complété par les étapes suivantes : - sélectionner une deuxième surface plane de préférence hrizontale, par exemple un sol de référence ; - acquérir des mesures au moyen du système lidar- SNI en orientant la direction de pointage vers le sol de référence selon une série d'orientations du système lidar prédéterminée ; - appliquer un algorithme de traitement des mesures acquises à l'étape précédente pour en extraire des valeurs d'un troisième angle de désalignement. Plus précisément, dans un mode de réalisation particulier, le procédé comporte une première étape et une deuxième étape principales. La première étape consiste à acquérir plusieurs lignes de scan successives, de telle manière que le lidar soit orienté par rapport à un mur vertical 8 selon les orientations définies par les angles suivants dans un repère local fixe défini par rapport au mur, où les axes )(local et Ylocal sont dans le plan du sol, l'axe )(local étant parallèle au mur, l'axe Ylocal perpendiculaire au mur et l'axe Ziocai perpendiculaire au sol : Orientation 1 : (r, t, c) = (0, 0, 0) Orientation 2 : (r, t, c) = (0, 0, 180 deg.) Orientation 3 : (r, t, c) = (0, 0, 70 deg.) Orientation 4 : (r, t, c) = (0, 0, -70 deg.) Le passage de l'orientation 1 à l'orientation 2 consiste par exemple à utiliser un miroir 15 double face, à faces planes et parallèles et à appliquer un demi-tour au système lidar. La deuxième étape consiste à scanner un plan horizontal 9, typiquement le sol, en penchant le système lidar dans le plan formé par les vecteurs Y - lidar et Zlidar selon les orientations définies par les angles suivants : Orientation 5 : (r, t, c) = (0, 70 deg., 0) Orientation 6 : (r, t, c) = (0, -70 deg., 0). Dans chacune des orientations 1, 2, 3, 4, 5, 6 définies ci-dessus, le système lidar acquiert au moins une ligne de scan. Avantageusement, l'acquisition de deux scans ou plus par orientation permet de calculer un scan moyen pour chaque orientation. Pour chaque point d'une ligne de scan, le système lidar acquiert une mesure de position et d'orientation angulaire du SNI associée à une mesure de distance entre le lidar et, respectivement, le mur 8 ou le sol 9. On observe que les orientations 1 et 2 n'ont pas besoin d'être précises. Une précision angulaire de ± 5 deg. environ est suffisante pour chacun des angles (r, t, c). D'autre part, dans les orientations 3 à 6, il est conseillé d'avoir des angles de lo tangage t et de cap c supérieurs à 60 degrés en valeur absolue. Toutefois, on observe qu'une valeur absolue d'angle de tangage supérieure à 30 degrés pour les orientations 5 et 6 suffit pour la détermination des angles de désalignement détaillée ci-après. On verra qu'un minimum de 9 observations est nécessaire pour estimer les angles de désalignement (d), 0, Y). Par conséquent, des mesures selon ces 6 15 orientations sont faites au moins deux fois. Les mouvements pour orienter le système de navigation et de pointage SNI-lidar par rapport au mur (ou par rapport au sol) peuvent être effectués manuellement, par exemple en plaçant le système lidar dans une cage qu'il est possible d'incliner ou en plaçant le système lidar sur une table permettant des mouvements de rotations. Dans 20 une alternative, les mouvements d'inclinaison ou de rotation sont motorisés. Eventuellement, on peut aussi embarquer le système lidar sur un véhicule (sur le toit d'une auto par exemple), auquel cas, le véhicule se déplace sur des plans inclinés pour faire les essais dans les différentes orientations du système lidar par rapport à un plan vertical et un plan horizontal. 25 Il n'est pas nécessaire d'utiliser un seul mur 8 de référence pour faire les acquisitions nécessaires à la calibration. Cependant, plus les acquisitions sont longues, plus les erreurs d'orientation de la centrale inertielle peuvent avoir un impact, du fait des dérives d'une centrale inertielle dans le temps. Si la centrale de navigation est éteinte puis rallumée, les biais de la centrale de navigation peuvent aussi affecter les résultats 30 des mesures car les biais auront changé. La série des orientations 1, 2, 3, 4, 5 et 6 définie ci-dessus n'est qu'un exemple illustratif et nullement limitatif. L'homme du métier détermine d'autres séries d'orientations également utiles pour déterminer les angles de désalignement. La sensibilité et la précision des mesures des trois angles de désalignement dépendent de 35 la sélection des différentes orientations. A cet égard, un point important dans la sélection des différentes orientations est d'effectuer des inclinaisons selon les différents axes pour les différentes orientations choisies, de façon à être sensible aux différents angles de désalignement. Et plus les orientations sont différentes l'une de l'autre, autrement dit plus l'amplitude des mouvements est grande pour passer d'une 40 orientation à une autre, et plus le résultat de la calibration sera précis. Algorithme de traitement des mesures L'algorithme de traitement de mesures doit calculer les angles de désalignement entre le SNI et le LiDAR (d), 0, Y). Un algorithme supportant les non linéarités est conseillé car on ne dispose pas de relation linéaire liant les observations aux angles de désalignement. Equation d'observation Il faut d'abord définir l'équation d'observation. Pour chaque scan statique d'un plan, on dispose d'un ensemble de points formant une ligne de scan. Chacune de ces lignes, caractérisée par son vecteur directeuruin , appartient à ce même plan et donc vérifie l'équation : fi = (t'in .Nn = 0 rn SNI (0 0 IA r- lidar )T Nn f = (`- SNI` Y' I lidar "i lidar N normale au plan 8 ou au plan 9 selon le plan auquel appartientui dans le repère [n] lidar : direction de la ligne de scan i dans le repère [lidar] SNI lidar : Matrice de passage approximée du repère [lidar]au repère [SNI] c SNI sNr : Matrice corrigeant le désalignement (d), 0, Y) entre le lidar et le SNI Cn SNI : Matrice de passage du repère [SNI] au repère de [n]. Comme mentionné précédemment, l'équation d'observation ne fait pas apparaitre d'information de position, ni de bras de levier, mais seulement des informations d'orientation. Elle est donc immune aux erreurs GPS. Par ailleurs, le fait de travailler dans le repère de navigation et non dans un repère global rend l'algorithme insensible aux biais du SNI car les erreurs de biais ne sont alors pas projetées lors du changement de repère. Dans cette équation, les paramètres à estimer sont les angles de désalignement (d), 0, Y) et la normale Nn qui n'est a priori pas connue. L'équation d'observation n'étant pas une fonction linéaire de (d), 0, Y), il est conseillé de choisir un algorithme adapté. Un algorithme de moindres carrés non linéaires est ici présenté. Mise en forme des données Les directions de ligne de scan ne sont pas directement fournies par le système LiDAR. Il faut les calculer. Pour chaque ligne de scan, calcul de la direction de la ligne de scan Pour chaque ligne de scan respectivement 41, 42, 43, 45, 46 acquise dans les différentes orientations 1 à 6 indiquées ci-dessus, on cherche à déterminer la direction respective de la ligne de scan 41, 42, 43, 45, 46. Ce calcul peut se faire à l'aide de différentes méthodes, par exemple la détermination du vecteur propre associé à la plus grande valeur propre issue d'une analyse en composantes principales sur l'ensemble des points de chaque ligne de scan, ou par régression linéaire, etc. On note Wide le vecteur directeur d'une ligne de scan i exprimé dans le repère du LiDAR, c'est-à-dire le vecteur directeur que l'on cherche et ûilidar l'approximation qui en est faite, c'est-à-dire le résultat du calcul effectué à partir des mesures de position des points d'une ligne de scan. A partir des quatre orientations 1, 2, 3 et 4 indiquées ci-dessus, on mesure donc quatre lignes de scan respectivement 41, 42, 43, 44. A partir des quatre lignes de scan 41, 42, 43, 44, on calcule donc quatre vecteurs directeurs approchés, respectivement eidar û2lidar û3lidar et û 4lidar qui sont en principe tous dans le plan du mur vertical 8. A partir des deux orientations 5 et 6 indiquées ci-dessus, on mesure deux lignes de scan respectivement 45 et 46. A partir des deux lignes de scan 45 et 46, on calcule deux vecteurs directeurs approchés, respectivement û5lidar et û6lidar qui sont en principe tous les deux dans le plan du sol horizontal 9. Calcul des matrices de passage Cis//,, et C sn : L'orientation relative entre le lidar et le SNI est connue approximativement à quelques degrés voir quelques dixièmes de degré. On utilise cette approximation pour calculer C1i//,. . La matrice de passage du repère SNI au repère de navigation est calculée à partir des angles fournis par le SNI. Application de l'algorithme des Moindres Carrés Non Linéaires L'algorithme des moindres carrés non linéaires est un algorithme itératif. C'est-à-dire que le paramètre contenant le vecteur d'état contenant les paramètres à estimer est mis à jour suivant la formule : Xk+1= Xk+ AX T Le vecteur d'état à l'itération (k) étant X (k) = [ço(k) 9(k) t(k) N8n (k) N9n (k)1 \ -1 Avec AX =(k)T J(k) ) J(k)', -F(k)) Où F(k) est le vecteur contenant les N observations à l'itération (k) : F(k) = fl ( k ) Nxl - - fN (k ) Et J(k) est la jacobienne de F calculée en Xk L'algorithme s'arrête lorsque le critère de convergence est atteint I IAXII D'autres algorithmes à la disposition de l'homme du métier peuvent être utilisés pour estimer les angles de désalignement (d), 0, Y) tel que des méthodes basées sur des méthodes de Quasi-Newton. Description détaillée d'un mode de réalisation du procédé de calibration On cherche ici à déterminer des orientations favorables du système lidar lors de l'acquisition des différentes lignes de scan utilisées pour la calibration. Pour un scan i, on note : ûin = uin fiin OÙ : - uin: représente le « vrai » vecteur directeur - ûin : représente un vecteur directeur estimé flirt: représente l'erreur sur le vecteur directeur, qui est en pratique majoritairement due aux angles de désalignement entre le repère approché du SNI et le repère propre au SNI, car les erreurs spécifiques du lidar peuvent être considérées à moyenne nulle. Pour la calibration des angles de désalignement, on cherche à maximiser l'erreur (uin, Nn : (fiin , Nn = (Ûin Uin , Nn = (Ûin , Nn où l'opérateur (.,.) représente le produit scalaire des deux termes. Cette erreur dépend donc de la normale au plan N et du vecteur directeur u, c'est-à-dire de l'orientation du lidar par rapport à ce plan. Le lidar 12 et la centrale de navigation 11 étant solidaires, on peut aussi dire que l'erreur dépend de l'orientation de la centrale de navigation 11 par rapport au plan. Définissons alors un repère lié au mur 8 : (Xmur, Ymur, Zmur) et représenté sur les figures 4, 5, 6 et 9. L'axe Xmur correspond avec la normale N8 du plan du mur et pointe dans la direction du lidar 12, l'axe Zmur est confondu avec le vecteur directeur u d'une ligne de scan et enfin l'axe Ymur est défini de manière à créer un repère orthonormé direct (Xmur, Ymur, Zmur). L'axe Xmur est donc fixe, tandis que les axes Ymur, Zmur dépendent de l'orientation du scan. Ainsi : LiDAR 0- mur et LiDAR 0 1 Définissons ensuite une matrice de passage du repère (Xmur, Ymur, Zmur) au repère approché [SNI] de la centrale de navigation 11, la matrice de passage étant définie par trois angles d'Euler (ct, p, y), telle que : cosG6)cos(y) - cos(a)sin(y)+ sin(a)sinG6)cos(y) sin(a)sin(y)+ cos(a)sinG6)cos(y) cr, cos(fl)sin(y) cos(a)cos(y)+ sin(a)sinG6)sin(y) - sin(a)cos(y)+ cos(a)sinG6)sin(y) - sinG6) sin(a)cos(ô3) cos(a)cos(ô3) Réécrivons les termes de l'équation d'erreur en fonction de cette matrice de changement de repère : Notons : Csnm, = Où les angles )sont les angles mesurés par la centrale. Ces angles ne 30 sont pas à confondre avec les angles de désalignement (cp, 0, Y). 25 1 mur o o 0 mur û in = Cn R (0 0 0),-,e1 Jidar SNI, boresight SNI' mur O 1 Nn = CSnNI Rboresight (ç , 6 1 ,Vf) )C SNI O Remplaçons ces expressions dans l'équation (iiin,Nn) afin de trouver les orientations pour lesquellesûin ne reste pas perpendiculaire au vecteur normal, sinon l'erreur n'est pas visible : ru mur Csmzr mur ÇCP 9 SNI 1 mur (flin,Nn) SNI 1 CSnNI, Rboresight ,yi.) c MUT 0 mur (flin,Nn) = [0 O c sn NI ilboresight 9 ) mur 0 1 0 0 L'expression peut être alors simplifiée : 1 O mur (t'in , Nn) = ([0 0 CsIr ) Rboresight (V, 8, ,Vf) 1 0 0 mur "\ ( C SNI mur 1 = C1 représente la première colonne de la matrice Cmsl''.` Où CZ et [0 0 1]C;uir = C3T représente la transposée de la troisième colonne de CmsC Dans le cas où les angles d'Euler de désalignement (0,0,0 sont petits, on peut faire l'approximation : ( 0 _y, e 0 e R boresight (0, I + o - O = I + (s-2 x) , avec n = 0 0 -e L'équation devient alors : (t'in , Nn c3T (1 +(nX))Ci (flin,Nn) = C3TC1 + C3T (n x) Les colonnes des matrices de changement de repère sont orthogonales, par conséquent leur produit scalaire est nul. Nn = T 12 (C1XC3) (iiin,Nn) = -52T C2 En développant Cz, on obtient une formule représentant l'erreur en fonction des 3 angles de désalignement : N") = 0(cos(a)sin(y)- sin (a)sin(fi)cos(y)) + 0(- cos(a)cos(y) - sin (a)sin (fi)sin (y)) + W- sin(a)cos(/3)) Appliquons cette formule dans les deux cas concrets évoqués plus haut : - balayage ou scanning d'un plan vertical 8 - balayage ou scanning d'un plan horizontal 9. Sur les figures 9 et 10 on a représenté le repère 21 de la centrale de navigation SNI, et le repère d'un mur pris comme plan vertical sur la figure 9 et du sol pris comme plan horizontal sur la figure 10. Dans la suite, l'orientation relative entre le SNI et le lidar correspond à celle représentée sur la figure 2. Plus précisément, la figure 9 représente 15 une vue de dessus du repère SNI et du repère du mur pendant un balayage du mur 8. La figure 10 représente une vue de côté du repère SNI et du repère du sol pendant un balayage du sol 9. Pour le scan sur un plan vertical 8 (cf. figure 9), on se limite à des variations de l'angle de cap. L'angle est un angle orienté défini par rapport à la parallèle au mur 20 dans un plan horizontal. Il correspond à l'angle de cap c défini plus haut. SNL Dans ce cas là, les valeurs des angles de la matrice de passage `"- a '172- fl=71- y=0 (iiin,Nn = (sin (71 - + cos(71 - = ÇO(COS + 9 ( -sin 25 On observe que l'erreur de cap 1P n'apparait pas dans l'expression ci-dessus. L'erreur de cap n'est donc pas mesurable en utilisant le produit scalaire avec la normale. D'autres orientations sont alors nécessaires pour estimer l'erreur de cap. On en déduit aussi que l'on peut estimer les erreurs de roulis et tangage indépendamment de l'erreur de cap. 30 Les orientations favorables à l'estimation de l'erreur de roulis sont celles correspondant à un SNI parallèle au mur, le plan LiDAR 27 étant perpendiculaire au mur. En effet dans ce cas =0deg. et (iiin,Nn) ço(cos Avec des angles proches de 0 ou 180 deg., on peut même considérer que l'erreur de roulis peut être estimée indépendamment des erreurs de cap et de tangage. Les orientations favorables à l'estimation de l'erreur de tangage sont celles pour lesquelles le SNI est presque perpendiculaire au mur, le plan Lidar 27 et le mur étant presque parallèles, cela veut dire que le mur est scanné de manière très rasante, à des angles proches de 90 ou -90 deg. En pratique, il est difficile d'avoir des angles supérieurs à 70 deg. Cet angle ne peut donc être estimé indépendamment de l'erreur de roulis. En conclusion, dans le cas du balayage d'un mur vertical 8 : - l'erreur de cap ne peut être estimée sur un plan vertical si le système lidar n'a pas de roulis/tangage ; - l'erreur de roulis est estimée en scannant le mur lorsque le système lidar est parallèle à ce dernier ; - l'erreur de tangage est estimée en scannant le mur lorsque le système lidar a une orientation la plus rasante possible par rapport au mur. Pour obtenir l'erreur de cap, on scanne une autre surface, par exemple une surface horizontale 9 au sol, comme illustré sur la figure 10. Dans le scan du plan horizontal 9, on se limite à des variations de l'angle de tangage. En effet, faire varier les angles de roulis et cap ne sert à rien dans cette configuration SNI-lidar, puisque cela ne ferait pas varier les angles (11,13,7). Sur la figure 10, l'angle est un angle orienté défini par rapport au sol. Il correspond à l'angle t défini plus haut (voir figure 7). Dans ce cas, les valeurs des angles de la matrice de passage C ',s1'/'.1 (a ,fl, y) : a - o 7 - -2Y2 (fiin,Nn)= yo(sin(71A -d +yi(cos(71A -d (iiin,Nn)= ço(cos)+ yi(sin) En conclusion, pendant le scan d'un plan horizontal : - l'erreur de tangage n'apparait pas dans l'expression, elle n'est donc pas visible en scannant un plan horizontal ; - les orientations favorables à l'estimation de l'erreur de roulis sont celles correspondant à un système lidar à plat sur le sol. En effet dans ce cas =Odeg. et Kiiin,Nn) :z-- ço(cos). - les orientations favorables à l'estimation de l'erreur de cap sont celles pour lesquelles le système lidar, présente un angle proche de -±7c/2 par rapport au sol. Par conséquent, l'erreur de tangage ne peut pas être estimée sur un plan horizontal, mais l'erreur de cap peut être estimée en scannant le sol avec un angle II suffisamment important. En conclusion, le procédé de calibration est basé sur le scan de deux surfaces planes, de préférence un mur vertical et un sol horizontal, dans différentes orientations du système Lidar. La position du système lidar est de préférence statique entre les différentes orientations, seule l'inclinaison du système peut varier. Ce procédé de calibration permet de déterminer les erreurs d'alignement par une procédure relativement simple consistant par exemple à scanner le mur, le système lidar étant lo parallèle au mur, c'est-à-dire le plan LiDAR 27 étant perpendiculaire au mur et l'angle de roulis nul (cf figure 9 pour = 0), puis à scanner le mur avec une incidence importante pour estimer les angles de désalignement de roulis et de tangage O. Enfin l'estimation de l'angle de désalignement en cap Y est faite grâce à des mesures acquises sur un plan horizontal avec un angle de tangage important du 15 système lidar (cf figures 7-8) afin de décorréler au maximum l'erreur de roulis et l'erreur de cap et de faire ressortir au maximum l'erreur de cap. Par rapport aux solutions généralement utilisées pour la calibration SNI-lidar, l'invention présente notamment les avantages suivants. Le procédé de calibration est insensible aux erreurs humaines. L'utilisateur n'a 20 pas à repérer des points remarquables sur des images, il n'y a donc pas de risques d'erreurs dans le choix de ces points. Le procédé de calibration est insensible aux erreurs de GPS. Le procédé est appliqué en une position statique, il n'est donc pas sensible aux erreurs de position et n'est pas dégradé par un fonctionnement anormal d'un GPS. 25 Le procédé de calibration est insensible aux biais du SNI. Les angles de désalignement (boresight) sont calculés dans un repère local de navigation, le calcul n'est pas affecté par des erreurs de cap ou d'attitude, car celles-ci sont constantes au cours des différentes mesures de scans. Le procédé de calibration supprime ainsi les trois principales sources d'erreurs 30 lors d'une calibration SNI-Lidar traditionnelle et permet donc d'obtenir une plus grande précision. Les méthodes classiques de calibration des angles de désalignement ne permettent pas de calculer un facteur de qualité fiable. Au contraire, le procédé de calibration proposé permet de calculer un écart-type représentatif des incertitudes liées 35 aux mesures de la calibration. Enfin, la mise en oeuvre du procédé de calibration proposé est simple. Les méthodes classiques proposées jusqu'à présent imposent généralement de mettre en place le SNI et le Lidar sur un véhicule ainsi que de trouver un site dégagé pour avoir un signal GPS RTK très précis. La méthode proposée permet de calibrer un système 40 Lidar en intérieur : il suffit de disposer d'un mur et d'une plateforme orientable pour générer plusieurs orientations du système. Ceci est particulièrement bien adapté à la fabrication de systèmes couplés SNI-lidar, qui peut ainsi réaliser la calibration en usine. L'utilisation d'une table à deux axes automatique peut permettre d'automatiser entièrement la calibration. Enfin le procédé permet de réduire considérablement la durée du procédé de calibration. Les méthodes classiques imposent de nombreuses manipulations lors des acquisitions de mesure puis lors du traitement des mesures de calibration. La méthode proposée permet d'effectuer une calibration en quelques minutes une fois le système de navigation correctement aligné. Enfin, il n'est pas nécessaire de disposer d'un lidar ayant une vitesse de rotation du miroir importante, puisque la calibration se fait de manière statique. REVENDICATIONS1. Procédé de calibration d'un système de navigation et de pointage (10) comprenant une centrale de navigation (11) reliée solidairement à un dispositif de pointage (12), la centrale de navigation étant configurée pour mesurer la position et l'orientation d'un repère (21) de ladite centrale de navigation (11), le dispositif de pointage (12) comprenant des moyens d'émission d'un faisceau dans une direction de pointage (17) et des moyens de balayage de la direction de pointage (17), le dispositif de pointage (12) étant configuré pour mesurer un signal de pointage, de préférence par réflexion ou rétro-diffusion dudit faisceau sur une cible, en fonction de la direction de pointage (17) dans un repère (22) lié au dispositif de pointage, le procédé de calibration comprenant les étapes suivantes de calibration des angles entre le repère (21) de ladite centrale de navigation (11) et le repère (22) lié au dispositif de pointage (12) : a. Estimation d'un repère approché (31) de navigation (11), le repère approché (31) de navigation étant incliné par rapport au repère (21) de navigation avec des angles d'Euler de désalignement (d), 0, Y) ; b. Sélection d'une première surface plane (8), de préférence verticale ; c. Détermination d'un vecteur normal (N8) à ladite première surface plane dans le repère (21) de navigation ; d. Orientation du système de navigation et de pointage (10) dans une pluralité d'orientations angulaires distinctes par rapport à la première surface plane (8) ladite pluralité d'orientations angulaires étant déterminées dans le repère (21) de la centrale de navigation (11) ; e. Balayage de la direction de pointage (17) de manière à acquérir, au moyen du dispositif de pointage (12), pour chacune des orientations de ladite pluralité d'orientations par rapport à la première surface plane (8), une série de mesures de pointage suivant une ligne de balayage (41, 42, 43) sur la première surface plane (8) ; f. Traitement de chacune desdites séries de mesures de pointage suivant une ligne de balayage (41, 42, 43) sur la première surface plane (8) de manière à déterminer, dans le repère repère (22) lié au dispositif de pointage, un vecteur directeur approché (fil , û2, û3) de chaque ligne de balayage (41, 42, 43) sur la première surface plane (8) ; g. Traitement d'une pluralité de vecteurs directeurs approchés (û,, û2, û3) d'une pluralité de lignes de balayage (41, 42, 43) sur la première surface plane pour en déduire un vecteur normal approché (1S'i 8) de la première surface plane (8) ; h. Calcul d'au moins un angle de désalignement (d), 0, Y) entre le vecteur normal approché (1S'18) déterminé à l'étape précédente dans le repère de navigation approché et le vecteur normal (N8) déterminé à l'étape c) dans le repère (21) de la centrale de navigation (11). 2. Procédé de calibration d'un système de navigation et de pointage (10) selon la revendication 1 comportant en outre les étapes suivantes : i. Sélection d'une deuxième surface plane (9), de préférence horizontale ; j. Détermination d'un vecteur normal (N9) à ladite deuxième surface plane (9) dans le repère (21) de navigation ; k. Orientation du système de navigation et de pointage (10) dans une pluralité d'orientations angulaires distinctes par rapport à la deuxième surface plane (9), ladite pluralité d'orientations angulaires étant déterminées dans le repère (21) de la centrale de navigation (11) ; I. Balayage de la direction de pointage (17) de manière à acquérir, au moyen du dispositif de pointage (12), pour chacune des orientations de ladite pluralité d'orientations par rapport à la deuxième surface plane (9), une série de mesures de pointage suivant une ligne de balayage (45, 46) sur la deuxième surface plane (9) ; m. Traitement de chacune desdites séries de mesures de pointage suivant une ligne de balayage (45, 46) sur la deuxième surface plane (9) de manière à déterminer, dans le repère repère (22) lié au dispositif de pointage, un vecteur directeur approché (t15, û6) de chaque ligne de balayage (45, 46) sur la deuxième surface plane (9) ; n. Traitement d'une pluralité de vecteurs directeurs approchés ( û5, û6) d'une pluralité de lignes de balayage (45, 46) sur la deuxième surface plane (9) pour en déduire un vecteur normal approché (1S'i 9) de la deuxième surface plane (9); o. Calcul d'au moins un autre angle de désalignement (d), 0, Y) entre le vecteur normal approché (1S'19) déterminé à l'étape précédente dans le repère de navigation approché et le vecteur normal (N9) déterminé à l'étape j) dans le repère (21) de la centrale de navigation (11). 3. Procédé de calibration selon l'une des revendications 1 à 2 dans lequel ladite pluralité d'orientations angulaires comporte les quatre orientations angulaires du système de navigation et de pointage (10) par rapport à la première surface plane (8) suivantes où r désigne l'angle de roulis, t l'angle de tangage et c l'angle de cap du 30 système de navigation et de pointage par rapport à un repère local ayant des axes X et Y dans le plan horizontal, et l'axe X parallèle au plan de la première surface plane : Orientation 1 : (r, t, c) = (0, 0, 0± 5 deg.) Orientation 2 : (r, t, c) = (0, 0, 180± 5 deg.) Orientation 3 : (r, t, c) = (0, 0, 55 ± 25 deg.) 35 Orientation 4 : (r, t, c) = (0, 0, -55 ± 25 deg.) et/ou les deux orientations angulaires suivantes du système de navigation et de pointage (10) par rapport à la deuxième surface plane (9) : - Orientation 5 : (r, t, c) = (0, 55 ± 25 deg., 0) - Orientation 6 : (r, t, c) = (0, -55 ± 25 deg., 0) 4. Procédé de calibration selon l'une des revendications précédentes, dans lequel l'étape e) et/ou l'étape I) comporte une pluralité de balayage d'une même ligne de balayage et en ce que ladite série de mesures de pointage suivant une ligne de balayage (41, 42, 43, 45, 46) est obtenue par un calcul de moyenne de la pluralité de balayage d'une même ligne de scan. 5. Procédé de calibration selon l'une des revendications précédentes, dans lequel l'étape g) et/ou l'étape n) de détermination d'un vecteur normal approché comporte une étape de détermination d'un vecteur propre associé à la plus petite valeur propre d'une analyse en composantes principales des séries de mesures de pointage lo suivant ladite pluralité de lignes de balayage ou une étape de calcul d'une moyenne d'une pluralité de vecteurs normaux mesurés par produit vectoriel des vecteurs directeurs approchés de deux lignes de scan différentes déterminés à l'étape f) et/ou à l'étape m). 6. Procédé de calibration selon l'une des revendications précédentes, dans 15 lequel l'étape h) et/ou l'étape o) de calcul des angles de désalignement (d), 0, Y) comporte une étape de maximisation d'un produit scalaire pour en déduire une estimation d'un angle de désalignement, ladite étape de maximisation comportant une étape de traitement par une méthode de régression des moindres carrés non linéaires ou par une méthode de régression des moindres carrés linéaires ou une 20 méthodes de Quasi-Newton. 7. Procédé de calibration selon l'une des revendications précédentes, comportant une étape de calcul d'un indicateur de précision des angles de désalignement (d), 0, Y). 8. Système de navigation et de pointage (10) comprenant une centrale de 25 navigation (11) reliée solidairement à un dispositif de pointage (12), la centrale de navigation étant configurée pour mesurer la position et l'orientation d'un repère (21) de ladite centrale de navigation (11), le dispositif de pointage (12) comprenant des moyens d'émission d'un faisceau dans une direction de pointage (17) et des moyens de balayage de la direction de pointage (17), le dispositif de pointage (12) étant configuré 30 pour mesurer un signal de pointage, de préférence par réflexion ou rétro-diffusion dudit faisceau sur une cible, en fonction de la direction de pointage (17) dans un repère (22) lié au dispositif de pointage, le système de navigation comprenant des moyens de traitement des mesures de position et d'orientation de la centrale de navigation et des moyens de traitement du signal de pointage, et des moyens de calibration adaptés pour 35 implémenter une calibration angulaire entre le repère (21) de ladite centrale de navigation (11) et le repère (22) de pointage suivant un procédé de calibration selon l'une des revendications 1 à 7. 9. Système de navigation et de pointage (10) selon la revendication 8 dans lequel ladite centrale de navigation (11) comporte un système d'orientation parmi unecentrale de navigation inertielle ou une centrale d'attitude de type AHRS et/ou un système de positionnement de type GPS. 10. Système de navigation et de pointage (10) selon l'une des revendications 8 ou 9 dans lequel le dispositif de pointage (12) comporte un lidar à balayage suivant un plan en deux dimensions, un lidar à balayage en trois dimensions ou un sonar. 11. Système de navigation et de pointage (10) selon l'une des revendications 8 à 10 comportant en outre des moyens d'orientations angulaire dudit système de navigation et de pointage (10), adaptés pour orienter ledit système de navigation et de pointage (10) dans une pluralité d'orientations angulaires distinctes par rapport à ladite première surface plane (8) et/ou, respectivement, à ladite deuxième surface plane (9).