Source : preprint arXiv:2507.08349v1 [cs.RO], 11 juillet 2025 · pas de DOI dans le PDF · Macau University of Science and Technology ; Institute for AI Industry Research (AIR), Tsinghua University · accès OPEN · PDF ExtrinsicCalibration2025.pdf, sha256 c27fe427c0dc6f0a…, 10 pages Lu le : 2026-11-02. LECTURE COMPLÈTE : texte intégral des pages 1 à 10 (les tableaux I à IV ont été relus sur les images des pages 8 et 9) ; la liste de références (p. 10) n’a été que parcourue. Statut summarized. Thème : fusion · calibration extrinsèque multi-LiDAR / GNSS-INS Métadonnées vérifiées dans le PDF : auteurs Junhui Wang, Yan Qiao, Chao Gao, Naiqi Wu ; année 2025 ; titre conforme (le PDF écrit « of Ground Vehicles » en fin de titre) ; pas de DOI.

Contexte et question

Les camions de mine à ciel ouvert portent plusieurs LiDAR pour couvrir les angles morts, et il faut les exprimer dans le repère du GNSS-INS (GINS) pour fusionner les données (p. 1). Les méthodes existantes exigent des cibles, des champs de vue qui se recouvrent, ou une bonne trajectoire ; et le mouvement quasi plan dégrade l’estimation de la translation verticale (p. 1). Question : calibrer sans cible, sans recouvrement, pour des LiDAR mécaniques comme à état solide.

Méthode

• Chaîne en six étapes (Fig. 1, p. 3) : initialisation, calibration LiDAR-GINS, alignement au sol, calibration multi-LiDAR, analyse du terrain, optimisation conjointe.
• Initialisation (§IV-B) : recherche globale par l’algorithme DIRECT, reprise de LiDAR-align ; seule la rotation est optimisée, la translation est fixée à zéro ; 200 itérations en exemple.
• LiDAR-GINS (§IV-C) : optimisation par lot d’un facteur d’appariement généralisé (distance de Mahalanobis, covariance par point sur k = 20 voisins), noyau de Huber, association par K-D tree.
• Facteur de hauteur d’installation du GINS (§III-B, Eq. 1 à 3) : mesure de la hauteur du GINS au-dessus du sol, qui contraint la direction non observable en mouvement plan. Covariance : cb = 10^8 sur cinq dimensions, cs = 10^-3 sur la hauteur.
• Alignement au sol (§IV-D) : on extrait le plan du sol vu par chaque LiDAR ; la différence d’ordonnée à l’origine donne une hauteur relative Δh qui corrige la translation LiDAR-GINS (Eq. 9 à 11).
• Analyse du terrain (§IV-F) : sol extrait par Patchwork++, rugosité locale (Eq. 16), le tronçon le plus lisse sert à estimer l’extrinsèque entre le LiDAR de base et un « LiDAR virtuel » dont le plan xy coïncide avec le sol (Eq. 17-18).
• Optimisation conjointe (§IV-G) : trois familles de facteurs (hauteur, appariements généralisés, cohérence de mouvement), optimisées par Levenberg-Marquardt (Eq. 23). Les poses GINS sont aussi des variables, car les vibrations en terrain accidenté les faussent (§IV-C).
• Données : images clés tous les 2 m ou 30° de virage ; capteurs synchronisés par horodatage matériel (§III-D).

Données

• Simulation (§V-A, Fig. 6) : trajectoire en huit ; un Livox Avia (état solide), un Velodyne 16 voies (rotatif) servant de LiDAR de base ; champs de vue sans recouvrement ; vérité terrain connue.
• Réel (§V-A, Fig. 7) : camion de mine à ciel ouvert, cinq LiDAR (L0 Livox Horizon, FoV horizontal 81,7° ; L1 à L4 Robosense Bpearl, FoV 360°) et un GINS de précision de position d’environ 3 cm, 10 Hz synchronisés. Pas de vérité terrain pour LiDAR-GINS ; pour le multi-LiDAR, référence obtenue en recalant chaque LiDAR sur un LiDAR 64 voies.
• Calcul : Intel i9-13980HK.

Résultats chiffrés

Métriques : MME et MPV (plus bas = carte plus cohérente), erreurs de rotation (°) et de translation (m) au sens de Lie (§V-B).

Simulation (Tab. II, p. 8)

Pour « Ours* », les auteurs ont ajouté du bruit à l’extrinsèque LiDAR de base / LiDAR virtuel pour simuler l’absence d’analyse du terrain (p. 8) ; ce n’est pas une exécution distincte de l’algorithme.

Réel, cohérence de carte LiDAR-GINS (Tab. I, p. 8) : sur L0, MME -3,24945 et MPV 0,182268 pour la méthode, contre -2,80735 et 0,236921 pour LiDAR-align et -2,81774 et 0,236816 pour HECalib. Sur L1 à L4, MME de -3,73137 à -4,26837 et MPV de 0,11642 à 0,163125 ; LiDAR-align donne MME de -2,92135 à -3,32815 et MPV de 0,189405 à 0,229419. HECalib n’estime que L0 (les autres sont notés en échec).

Réel, multi-LiDAR (Tab. III, p. 9) : RMSE globale de la méthode 0,293° et 0,081 m, contre 0,348° et 0,116 m pour GICP, 0,473° et 0,241 m pour VGICP, 0,778° et 0,246 m pour NDT. GICP fait mieux sur L2 (0,158° ; 0,022 m contre 0,348° ; 0,037 m). Sur L4, la méthode donne 0,422° et 0,149 m ; les auteurs jugent la précision de L4 médiocre pour toutes les méthodes, à cause du FoV limité du LiDAR arrière (p. 9). HECalib échoue sur L1 à L4 (p. 9). Pour GICP, VGICP et NDT, les auteurs ont choisi à la main des nuages favorables au recalage ; cela n’a pas été fait pour leur méthode (p. 8-9).

Temps de calcul (Tab. IV, p. 9) : 12,279 s au total en simulation (29x2 scans), 148,165 s en réel (19x5 scans), dont 103,185 s pour l’optimisation conjointe.

Limites déclarées par les auteurs

• Sans analyse du terrain, la performance baisse (p. 8).
• La charge de calcul croît « exponentiellement » avec le nombre de capteurs et de données ; une fenêtre glissante la réduirait, au prix de la précision (p. 9).
• L’extrinsèque du LiDAR arrière L4 reste peu précis pour toutes les méthodes (p. 9).
• Travaux futurs : ajout de caméras (p. 9).
• Le recalage manuel des nuages favorise les méthodes de comparaison, pas la leur (p. 8-9).

Ce que ça change pour Teleforge (interprétation, hors de l’article)

• L’idée du facteur de hauteur est transposable : toute plateforme quasi plane laisse une translation verticale non observable, et une mesure indépendante de hauteur la contraint. C’est la même logique d’observabilité que pour nos capteurs inertiels au sol.
• Le papier n’utilise que les poses du GINS (Eq. 6 et 22), pas de mesures inertielles brutes : le GINS est traité comme une trajectoire fournie, puis raffinée. Il ne règle donc pas la calibration d’un accéléromètre embarqué, seulement une idée de méthode.
• Chiffres à relativiser : un seul véhicule, un seul site réel, pas de vérité terrain pour LiDAR-GINS, et comparaison à des méthodes parfois mal servies (HECalib en échec, nuages choisis à la main pour GICP, VGICP, NDT).
• Code public annoncé par les auteurs (p. 2) : à tester avant de reprendre quoi que ce soit.