Source : The Baltic Journal of Road and Bridge Engineering 6(3): 163–168 · DOI 10.3846/bjrbe.2011.21 · Department of Transport, University of Ljubljana (Slovénie) · accès OPEN · PDF Lipar2011-GPS-Curvature.pdf, sha256 df69b7dc18c2c6fe…, 6 pages Lu le : 2026-11-02. LECTURE COMPLÈTE (texte extrait, pages 4 à 6 relues par vision pour les figures et équations 12 à 14). Les équations 1 à 11 ne sont pas reproduites ici. Statut summarized. Thème : fusion · courbure de la géométrie routière depuis le GPS seul Métadonnées vérifiées dans le PDF : titre, auteurs (4), année 2011 et DOI conformes. L’ancienne fiche ne citait que « Lipar, P., et al. ».

Contexte et question

Les applications avancées (avertissement de survitesse en courbe, alerte de sortie de voie, navigation à la voie) demandent une géométrie routière précise, axe et courbure. Les auteurs écartent les véhicules instrumentés et la photo-interprétation aérienne, et retiennent l’agrégation de traces GPS bruitées de flottes. La question : estimer la courbure d’un axe routier à partir de points GPS bruts, non équidistants, sans définition stricte de la courbure d’une courbe discrète.

Méthode

• Lissage des points bruts par splines cubiques (B-splines, jonction C2), puis rééchantillonnage en sommets équidistants. Les écarts entre points GPS bruts « vary from few to several 10 m » (p. 164).
• Segmentation de la courbe en arcs de cercle : on part de 10 à 20 points, on ajoute un point à chaque étape, on suit l’erreur (distance maximale des points à l’arc) et on ouvre un nouvel arc quand elle dépasse la tolérance prescrite (exemple : 1 m) (§3, p. 164).
• Ajustement d’un arc en deux temps : projection stéréographique des points sur la sphère unité, ajustement par moindres carrés du plan dont l’intersection avec la sphère épouse le mieux les points, puis projection stéréographique inverse du cercle obtenu vers le plan équatorial. On obtient l’arc, son rayon (éq. 12) et son centre (éq. 13) (p. 166).
• Le plan est trouvé par moindres carrés d’une fonction affine, puis par itération dans une base orthonormée dont le dernier vecteur est normal au plan ; la convergence est dite « extremely fast », quelques pas suffisent (p. 166).
• Un segment quasi rectiligne donne un grand rayon sans problème numérique, d’après les auteurs (p. 164). La méthode vaut en dimension supérieure ; le cas 3D est présenté numériquement.
• Implémentation donnée en Mathematica (Fig. 4, p. 166).

Données

• Courbe synthétique : deux arcs de rayons 300 m et 200 m reliés par un segment droit, points espacés de 20 à 35 m, chaque point perturbé uniformément dans un carré de 2 m par 2 m, censé correspondre à l’erreur de position d’un GPS (§3, p. 166, Fig. 5).
• Routes réelles : exemples 2D en Figs. 6 (route urbaine, petits rayons et segments droits courts) et 7 (route avec points bruts, points équidistants à 1 m, extrémités d’arcs, rayons, abscisses de segmentation).
• Modèle 3D réel du réseau routier d’État slovène, « more than 6000 km » (§4, p. 167). Ni le récepteur GPS, ni la fréquence, ni la vérité terrain de ce jeu ne sont décrits.

Résultats chiffrés

Ancres des chiffres principaux :

• Rayons retrouvés : « of R = 299 m and R = 200 m respectively » (p. 166).
• Segment droit : « R = 64 205 m » (p. 166).
• Perturbation : « perturbed inside a square of 2 m by » 2 m (p. 166).
• Tolérance : « exceeds the prescribed error (e.g. 1 m) » (p. 164).
• Réseau : « (more than 6000 km), and the » performance est jugée « very satisfactory » (p. 167).

Le seul essai chiffré est la courbe synthétique. Sur les données réelles, l’article ne donne aucune erreur de rayon ni comparaison à une référence : la performance y est qualifiée, pas mesurée.

Limites déclarées par les auteurs

• « Estimation of clothoide transition lines still has to be explored » (p. 167) : seuls arcs et droites sont ajustés, pas les clothoïdes.
• Les auteurs notent qu’il n’existe pas de méthode numérique définitive d’estimation de la courbure d’une courbe discrète (p. 164).

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

• La méthode traite la géométrie de l’axe à partir de positions seules : elle fournit une référence GNSS-seul pour comparer les estimations de courbure issues d’une fusion GNSS+IMU.
• Elle est conçue pour des traces agrégées de flotte et un tracé fixe de l’axe, pas pour une trace unique à fréquence faible : à ne pas lire comme une validation sur trace individuelle.
• Le seul chiffre d’exactitude est synthétique (bruit uniforme borné de 2 m) ; le rayon de 64 205 m pour une droite montre qu’un segment rectiligne sort avec un très grand rayon, pas avec une courbure nulle exacte : prévoir un seuil de courbure minimale.
• La tolérance de segmentation (1 m dans l’exemple) est un paramètre de réglage qui fixe la finesse des arcs : à rendre explicite si la méthode est reprise.
• L’article ne compare pas à une autre méthode ni à une courbure de référence mesurée : toute reprise exige un banc propre.