FR 3014553 A1 20150612 FR 1362390 A 20131211 L'invention concerne la calibration angulaire de la position d'une caméra embarquée dans un véhicule, notamment d'un véhicule automobile. Elle s'applique en particulier à la calibration initiale d'un système d'aide à la conduite installé aussi bien en usine, "en première monte", sans banc de calibration spécifique qu'en "deuxième monte" par un utilisateur sur son véhicule. Un autre exemple de situation où une calibration de l'orientation de la caméra est nécessaire est l'utilisation d'une application pour smartphone qui utilise l'image captée par la caméra de ce smartphone monté sur un support fixé au tableau de bord ou sur le pare-brise du véhi- 1 0 cule. On connait de nombreux systèmes d'assistance à la conduite mettant en oeuvre l'image captée par une caméra fixée dans le véhicule et qui "regarde" la route. Des algorithmes permettent de réaliser diverses fonctions telles que : détection des sorties de route, estimation de la distance avec 15 les autres véhicules, alerte anticollision, détection d'obstacle, etc. Pour toutes ces applications, il est nécessaire de connaitre précisément la position et l'orientation de la caméra par rapport au véhicule. Dans toutes ces situations, le montage de la caméra peut être très approximatif, ce qui nécessite des corrections importantes à appliquer au si- 20 gnal délivré par la caméra avant de pouvoir analyser convenablement ce- lui-ci. Il peut être également nécessaire de procéder à un réajustement ultérieur des paramètres de calibration, par exemple pour compenser les déformations que l'installation peut subir. L'installation d'une caméra à un endroit précis est difficile du fait de la di- 25 versité des modèles de véhicules et des contraintes techniques. S'il est relativement simple de déterminer la position de la caméra dans le véhicule en termes de degrés de liberté en translation (mesure de la distance entre le point focal de la caméra et le centre du véhicule ou le niveau du sol) il est en revanche beaucoup plus difficile de mesurer précisément 30 l'orientation de la caméra par rapport à celle du véhicule en termes de degrés de liberté en rotation, c'est-à-dire de déterminer l'écart angulaire entre l'axe optique de la caméra et l'axe longitudinal du véhicule. Le but de l'invention est de proposer une méthode de calibration automatique de l'orientation d'une caméra embarquée dans un véhicule, avec 35 une caméra dont les paramètres intrinsèques (longueur focale, résolution, distorsion) sont connus, mais dont l'orientation relative par rapport au véhicule n'est pas connue. Comme on le verra, cette méthode ne nécessite aucune intervention du conducteur et elle est mise en oeuvre directement à partir de l'image de la scène captée par la caméra lorsque le véhicule roule normalement. On verra également que cette méthode est applicable aussi bien à une caméra située à l'avant du véhicule qu'à l'arrière de celui-ci, par exemple pour une caméra de recul ou d'estimation de la distance avec le véhicule qui suit le véhicule équipé. La calibration angulaire consiste à évaluer une rotation (pouvant être ex- primée sous forme matricielle) entre un repère de coordonnées lié à la caméra et un repère de coordonnées lié au véhicule. Cette rotation peut être décomposée en trois rotations élémentaires définies respectivement par des angles de tangage, de roulis et de lacet. Une estimation complète de la rotation entre la caméra et le véhicule nécessite donc de déterminer ces trois angles. De nombreuses méthodes de calibration angulaire ont déjà été proposées. Parmi ces méthodes, un certain nombre d'entre elles permettent d'estimer complètement la rotation entre la caméra et le véhicule, mais nécessitent une intervention humaine et/ou une étape préalable, par exemple du fait du recours à une grille spécifique fixée au système, à des cibles en damier posées au sol ou contre un mur, ou à des marqueurs fixés sur le capot du véhicule et visibles par la caméra. Les autres méthodes, qui n'utilisent pas d'élément extérieur et ne nécessi- tent aucune intervention humaine, n'ont jusqu'à présent permis qu'une calibration partielle, tous les degrés de liberté de la rotation (en tangage, en roulis et en lacet) ne pouvant être estimés ensemble. En particulier, de nombreuses méthodes ont été proposées pour estimer les rotations en tangage et en lacet, par exemple en analysant la position de l'horizon et/ou du point de fuite dans l'image de la scène captée par la caméra, ou bien à partir de points d'intérêt situés sur le sol et dont le mouvement est suivi entre les images successives. Mais ces techniques ne permettent pas d'estimer automatiquement l'angle de roulis entre la caméra et le véhicule. On mentionnera seulement en passant les techniques impliquant des paires de caméras, qui produisent une représentation stéréoscopique de la scène. Il est beaucoup plus simple de calibrer la position d'une paire de caméras, mais les techniques mises en oeuvre ne sont pas transposables à une caméra unique, comme dans le cas présent. La présente invention a pour but de résoudre les problèmes précités, en proposant une technique de calibration angulaire complète - c'est-à-dire selon les trois degrés de liberté en tangage, roulis et lacet - de l'orientation d'une caméra (unique) embarquée dans un véhicule, ceci de façon entièrement automatique pendant une phase de conduite sur route et sans intervention de l'utilisateur ni recours à un quelconque marqueur ou cible lors d'une phase préalable. L'idée de base de l'invention consiste, lorsque le véhicule roule, à capter les images des marquages au sol délimitant les voies de circulation sur la route, et par un processus itératif estimer complètement l'orientation de la caméra par rapport au véhicule à partir de la position de deux voies situées côte-à-côte dans l'image, par exemple la voie centrale (celle sur laquelle roule le véhicule) et la voie à gauche et/ou à droite de celle-ci. En particulier, comme on le verra plus loin, l'invention ne nécessite pas de connaitre précisément la hauteur du véhicule ni la largeur des voies. Seules des quantités approchées sont utilisées pour initialiser ces valeurs, qui seront estimées par l'algorithme. Plus précisément, l'invention propose un procédé comprenant, pendant une phase de déplacement du véhicule sur une route, des étapes de : a) capture d'une séquence d'images video d'une scène de la route, chaque image de la séquence comprenant une pluralité de traces d'intérêt représentatives de marquages au sol de voies de circulation parallèles ; et b) détection de bords dans les images video, pour l'obtention d'une liste de bords rectilignes correspondant, dans chaque image, auxdites traces d'intérêt. De façon caractéristique de l'invention, ce procédé comprend en outre les étapes successives suivantes : c) correction de la liste de bords par application à chaque bord d'une ro- tation selon l'estimée de la matrice de calibration ; d) calcul d'angles de tangage et de lacet résiduels par une méthode de détection du point de fuite appliqué à la liste de bords corrigée obtenue à l'étape c) ; e) mise à jour de l'estimée de la matrice de calibration par application d'une rotation selon les angles de tangage et de lacet résiduels calcu- lés à l'étape d) ; f) calcul d'un angle de roulis résiduel à partir de la liste de bords corrigée obtenue à l'étape c) et de l'estimée de la matrice de calibration obtenue à l'étape e), par : - estimation des largeurs respectives d'au moins deux voies de cir- culation à partir de la liste de bords ; - comparaison des largeurs respectives ainsi estimées et calcul d'au moins une différence de largeur entre voies ; et calcul de l'angle de roulis en fonction de ladite différence de lar- geur ainsi calculée, les largeurs réelles des voies étant supposées identiques ; g) mise à jour de l'estimée de la matrice de calibration par application d'une rotation selon l'angle de roulis résiduel calculé à l'étape f) ; et h) itération, sur au moins une séquence d'images, des étapes c) à g) jusqu'à ce que les angles de tangage, de lacet et de roulis résiduels soient inférieurs à un seuil prédéterminé. Très avantageusement, l'étape f) d'estimation de l'angle de roulis comprend les étapes successives suivantes : f1) correction de la liste de bords par application à chaque bord d'une ro- tation selon l'estimée de la matrice de calibration et d'une rotation se- lon l'estimée d'une matrice de roulis ; f2) pour chaque image, à partir de la liste de bords obtenue à l'étape f1) et en fonction d'une référence de largeur de voie : - calcul de la largeur de chacune des voies de circulation détectée dans l'image, et - calcul d'un angle de roulis résiduel pour chaque couple de voies de circulation détecté dans l'image ; f3) pour la séquence d'images, calcul d'une largeur moyenne en fonction des largeurs de voie estimées à l'étape f2) pour les différentes images de la séquence ; f4) pour la séquence d'images, calcul d'un angle de roulis résiduel moyen à partir des angles de roulis résiduels calculés à l'étape f2) pour les différentes images de la séquence ; f5) mise à jour : - de la matrice de roulis par application d'une rotation selon l'angle de roulis résiduel moyen calculé à l'étape f4), et - de la largeur de référence moyenne par la largeur moyenne calculée à l'étape f3) ; f6) itération des étapes f1) à f5) jusqu'à ce que l'angle de roulis résiduel soit inférieur à un seuil prédéterminé ; et f7) extraction d'un angle de roulis estimé à partir de la matrice de roulis. Dans une forme de mise en oeuvre préférentielle, l'étape f2) de calcul de la largeur de chacune des voies de circulation et de l'angle de roulis pour chaque couple de voies de circulation comprend : - en cas de détection dans l'image courante de plus de deux voies de circulation, la délivrance en sortie de trois valeurs de largeur de voie et de deux valeurs d'angle de roulis résiduel ; - en cas de détection dans l'image courante de deux voies de circulation, la délivrance en sortie seulement de deux valeurs de largeur de voie et d'une valeur d'angle de roulis ; - en cas de détection dans l'image courante d'une seule voie de circulation, la délivrance en sortie seulement d'une valeur de largeur de voie. L'estimation de la largeur de chacune des voies de circulation peut notamment comprendre l'application de la relation : i7v=iiv COS (U)( a,-a, W étant l'estimée de la largeur de la voie, H, étant la hauteur de la caméra par rapport au sol, a étant l'angle du véhicule par rapport à l'axe de la voie sur lequel il circule, et al et a2 étant les coefficients directeurs des bords délimitant la voie de circulation dans l'image. À l'étape f2) l'estimation de l'angle de roulis peut comprendre l'application de la relation : X cos (c0 H, W étant la référence de largeur de voie ou son estimation, vs7y et VÇ étant les estimées de largeur de voies du couple de voies considéré, a étant l'angle du véhicule par rapport à l'axe de la voie sur lequel il cir- cule, et H, étant la hauteur de la caméra par rapport au sol. 0 On va maintenant décrire un exemple de mise en oeuvre de l'invention, en référence aux dessins annexés où les mêmes références désignent d'une figure à l'autre des éléments identiques ou fonctionnellement semblables. La Figure 1 est une représentation schématique d'un véhicule équipé d'une caméra frontale, montrant les repères de coordonnées respectifs impliqués dans le procédé de calibration de l'invention. La Figure 2 est un exemple de scène captée par la caméra, montrant les marquages au sol délimitant les voies de circulation, ainsi que la géométrie des droites correspondant à cette représentation. La Figure 3 est un schéma général par blocs fonctionnels, explicitant le mode opératoire du procédé de calibration de l'invention. La Figure 4 est un schéma par blocs fonctionnels détaillant l'algorithme d'estimation du roulis du schéma de la Figure 3. La Figure 5 est un schéma par blocs fonctionnels détaillant l'algorithme d'analyse de la largeur des voies du schéma de la Figure 4. 0 La Figure 1 illustre les différents repères utilisés pour définir la position et l'orientation de la caméra C par rapport au véhicule. II s'agit des repères suivants : - le repère caméra (X,, Y, 4), qui correspond aux axes de la caméra, l'axe Z, étant l'axe optique et les axes X, et Y correspondant aux axes du plan focal pour un modèle de sténopé ; - le repère véhicule (X,, Y', 4), où l'axe OZ, correspond à l'axe longitudi- nal du véhicule et le plan (0X,, 04) correspond au sol ; et - le repère global (Xg, Yg, Zg) basé sur la route, l'axe Zg correspondant à l'axe de la route. L'objectif de la calibration est de déterminer la rotation qui existe entre le repère véhicule et le repère caméra. Soit un point quelconque. Si l'on appelle P, ses coordonnées dans le re- père véhicule et P, celles dans le repère caméra, on a : X,, X' = RX yr z, Pc= +T =Rx P,,+T) où T est la translation entre les centres O et C des deux repères. Par définition du repère véhicule, on a : T= 0 -H , 0 où H, est la distance de la caméra au point O dans le repère véhicule (hauteur). On peut également écrire : R= Rot xx Rot yx Rot z avec : 1 0 0 coss(ay) 0 sin ( a, ) cos (az) sin (az) 0 Roty= 0 coscrs) sin 0 t 0 Ro z= -sin (az) cos (az 0 0 -sin (a x') cos (a ) -sin (ay) 0 cos (ay ) 0 0 1 ax, ay, et a, correspondant respectivement aux angles de tangage, lacet et roulis. La Figure 2 illustre un exemple de scène captée par la caméra, montrant des marquages au sol délimitant des voies de circulation, avec dans cet exemple une voie centrale VC, une voie de droite VD et une voie de gauche VG. Ces voies sont définies dans l'image par des droites D1 ... D4 correspondant à la trace dans l'image recueillie par la caméra des marquages au sol délimitant ces voies de circulation, marquages qui seront par la suite appelés "bords". Plus généralement, on notera Dk les bords de voies obtenus sur une image, et I(Dal la liste des bords obtenus sur une séquence d'images consécutives. On notera Re l'estimation de la rotation. À partir de cette estimation, on peut facilement calculer la position des bords dans l'image d'une caméra virtuelle en rotation R e par rapport à la caméra. On notera /5k ces bords, ou bords corrigés. On notera )1 la liste des bords corrigés pour une séquence d'images. L'extraction des bords, c'est-à-dire la détermination des paramètres défi- nissant les droites Dk dans l'image, est effectuée par une technique d'analyse d'image conventionnelle, par exemple par détection de contours (qui transforment l'image en un tableau de valeurs binaires auquel est appliqué une technique de calcul des gradients), puis exécution d'une transformation de Hough, qui est une technique bien connue (dont le principe est notamment exposé dans le US 3 069 654 A) permettant de détecter et de caractériser très rapidement des droites dans une image numérisée. Dans le cas d'images présentant une distorsion (c'est-à-dire où les traces des bords ne seraient pas rectilignes dans l'image brute captée par la caméra), cette distorsion sera corrigée avant de procéder à la détection des bords, ces derniers étant supposés rectilignes. Algorithme de calibration automatique selon l'invention On va décrire en référence au schéma général de la Figure 3 (dont cer- taines parties seront explicitées plus en détail par les Figures 4 et 5) le mode opératoire du procédé de calibration de l'invention, permettant de caractériser la position angulaire de la caméra par rapport au véhicule, c'est-à-dire de déterminer les angles respectifs ax, ay, et a, de tangage, de lacet et de roulis définis plus haut. Les schémas des Figures 3 à 5 sont des représentations par blocs fonctionnels illustrant les différentes étapes du processus de calibration. On notera que, bien que ces schémas soient présentés sous forme de blocs interconnectés, la mise en oeuvre des différentes fonctions est logicielle, cette représentation n'étant qu'illustrative. Du fait de son caractère non linéaire, le procédé de l'invention met en oeuvre un algorithme itératif, cet algorithme pouvant être itéré soit sur une séquence d'images, soit sur plusieurs séquences d'images en changeant de séquence à chaque itération. Cette deuxième possibilité permet de li- miter l'impact des erreurs de mesures sur une partie de la vidéo ; on ob- tiendra alors une première estimation des angles de tangage, roulis et lacet, estimation qui sera ensuite affinée au fur et à mesure du temps et des différentes séquences ultérieures. On notera également que le procédé de l'invention peut permettre aussi bien une calibration initiale qu'une recalibration ultérieure par correction fine des déformations mécaniques susceptibles d'apparaitre pendant l'utilisation de la caméra, par exemple un déplacement de cette dernière sur son support en raison des vibrations subies, des variations de température, etc. On notera enfin que le procédé de l'invention est applicable à tout type de caméra installée dans un véhicule et observant la route, qu'il s'agisse aussi bien d'une caméra située à l'avant de la voiture qu'une caméra de recul, dès lors que la scène captée comprend au moins deux voies de circulation parallèles de largeur constante. Essentiellement, la calibration comprend les étapes suivantes : - correction de la position des bords avec l'estimation de la rotation Re (bloc 12) ; - estimation du tangage et du lacet résiduels (bloc 16) ; - mise à jour de la matrice de rotation (bloc 18) ; - estimation du roulis résiduel (bloc 20) ; - mise à jour de la matrice de rotation (bloc 24) ; Ces étapes sont itérées jusqu'à ce que les angles correctifs estimés par chaque module soient négligeables (test du bloc 22). Plus précisément, un traitement de détection des bords (bloc 10) est ap- pliqué à une séquence d'images en entrée, ce qui permet d'obtenir une liste de bords I(Dal . Cette liste de bords est corrigée (bloc 12) par application d'une rotation R qui est l'estimée, pour l'itération courante, de la rotation recherchée. Initialement, cette estimée R e est initialisée par la matrice identité Id, comme cela est schématisé par le commutateur 14. À partir de la liste de bords corrigés r(k)~}, l'algorithme estime (bloc 16) les angles de tangage et de lacet résiduels ax,ay correspondant au signal d'erreur de l'algorithme itératif selon ces deux composantes. Cette estimation opérée par le bloc 16 est effectuée avantageusement par la méthode du point de fuite, qui est une méthode en elle-même bien connue : essentiellement, on considère le point de fuite dans l'image, qui est l'intersection de tous les bords Di, de la route dans l'image. Si l'on désigne (uTp, vi-p) les coordonnées de ce point de fuite dans l'image et f la longueur focale de la caméra, si le véhicule est parfaitement aligné avec la route le repère global (Xg, Yg, Zg) lié au sol et le repère (Xv, 17', 4) lié au véhicule sont identiques, et l'on obtient les résultats suivants : tan (ay) uvp-= cos ,ax)X f v, =tan (ax)x f quel que soit l'angle de roulis a7 . On peut ainsi obtenir les angles av de tangage et ay de lacet. En pratique, le véhicule n'est jamais parfaitement aligné avec l'axe de la route, de sorte que les équations précédentes ne sont plus exactes. Pour compenser ce fait, on utilise une séquence d'images suffisamment longue pour que le véhicule soit considéré comme globalement aligné sur la route sur la durée de cette séquence et on calcule le point d'intersection des bords pour chaque image, puis on réalise une moyenne de tous les points d'intersection ainsi calculés, de manière à définir le point de fuite dont les coordonnées permettront ensuite d'estimer le tangage et le lacet. L'étape suivante (bloc 18) consiste à mettre à jour la matrice de rotation Re en appliquant à la valeur Re de la précédente itération une rotation de compensation par les valeurs de tangage et de lacet résiduelles détermi- nées par le bloc 16 : Re Les valeurs de tangage, lacet et roulis résiduelles à,,dy c et ci , qui cons- tituent le signal d'erreur de la boucle, sont testées (bloc 22) et l'algorithme est itéré si nécessaire par une nouvelle mise à jour Re4-Rot XR, (bloc 24) pour opérer une correction en roulis (le tangage et le lacet ayant déjà été corrigés au niveau du bloc 18). La matrice de rotation résultante Re constituera la nouvelle valeur d'entrée de l'algorithme au bloc 12 pour l'itération suivante. Le processus est ainsi poursuivi jusqu'à ce que les angles de tangage, lacet et roulis résiduels soient considérés comme négligeables (test du bloc 22). Si tel est le cas, cela signifie que l'algorithme a convergé, et la valeur estimée finale Re peut être alors délivrée en sortie, après application d'une dernière mise à jour R,4-Rot (c7,)XR, (bloc 26, identique au bloc 24) pour une correction en roulis par la dernière valeur ci, précédemment déterminée. Module d'estimation de l'angle de roulis La Figure 4 illustre de façon plus détaillée la manière dont est opérée la détermination du roulis résiduel à par le bloc 20 de la Figure 3. Le procédé est basé sur le principe suivant : pour une image donnée, on estime la largeur des voies visibles sur l'image (voie centrale, et voies gauche et/ou droite qui ont été détectées) : on pourra estimer le roulis en utilisant la différence entre les largeurs estimées et la largeur réelle des voies. On notera R7 la rotation due au roulis et W, l'estimation de la largeur de référence des voies, qui n'est pas connue. Afin de converger vers une estimation correcte du roulis et de la largeur réelle, l'algorithme réalise plusieurs itérations. À chaque itération, on utilise l'estimation Rz précédente pour corriger les bords, puis on estime l'angle de roulis résiduel ô a7 . L'algorithme a pour entrées : - une séquence d'images, sur lesquelles le module de détection de bords (bloc 10) est appliqué. On obtient la liste des bords 1( ; - une matrice de correction Re qui provient du module de correction du tangage et du lacet (bloc 16 de la Figure 3), et éventuellement des ité- rations précédentes de l'ensemble de l'algorithme ; - une initialisation (schématisée par le commutateur 30) de la largeur de référence de la route, par exemple We = 3 m et de l'estimation de la rotation due au roulis : Rz=Id . Les étapes suivantes sont itérées jusqu'à ce que le roulis résiduel ôa7 soit suffisamment faible (test du bloc 38) : - correction des bords avec Rz XR, pour corriger le tangage, le lacet et une partie du roulis (bloc 28) ; - Pour chaque image : estimation des largeurs et du roulis pour chaque paire de voies (bloc 32, dont le détail sera décrit plus bas en référence à la Figure 5) ; - moyenne des angles sur toute la séquence, donnant le roulis résiduel ô az (bloc 34) ; - estimation de la nouvelle largeur W à partir des largeurs estimées pendant toute la séquence (bloc 36) ; - mise à jour de Rz (bloc 40). En ce qui concerne plus précisément l'étape d'estimation de la largeur de référence (bloc 36), pour estimer la largeur réelle des voies on utilise toute la séquence d'images, sur laquelle on calcule les largeurs (par image) des voies de droite, de gauche et centrale. On calcule ensuite la moyenne des largeurs par type de voie, et on utilise la largeur minimale obtenue : l'utilisation de cette valeur minimale permet de limiter les risques que le facteur sous la racine soit négatif. Dans la séquence d'images, le module ne détecte pas toujours plusieurs voies. Lorsque sur une image seule la voie centrale est détectée, on utili- sera tout de même l'estimation de sa largeur pour estimer la largeur de référence, mais on ne pourra pas calculer d'estimation du roulis pour cette image. Au moment du calcul de la moyenne des roulis par image, on vérifie qu'il y a eu suffisamment d'images pour lesquelles plusieurs voies ont été dé- tectées. Au fur et à mesure des itérations, l'estimation de la largeur réelle des voies converge (test du bloc 38). L'algorithme s'achève alors par extraction (bloc 42) de la valeur de l'angle a7 tel que R7=Rot (az) , permettant de mesurer l'importance de la correc- tion en roulis. Module d'estimation de la largeur des voies et du roulis résiduel La Figure 5 illustre plus en détail les opérations exécutées par le bloc 32 de la Figure 4. Pour une image donnée, l'estimation est opérée de la manière suivante : - détection des bords de l'image Dk. (bloc 10), le module détectant éventuellement une voie à droite ou à gauche en plus de la voie centrale ; - correction de la rotation estimée, en particulier pour corriger le lacet et le tangage (bloc 28, de la Figure 4) ; - estimation de la largeur de chaque voie visible : voie centrale (bloc 44) et voies de droite et/ou de gauche (blocs 46) si celles-ci sont visibles sur l'image de la caméra ; - pour chaque couple de voies : estimation du roulis résiduel (bloc 48). Si trois voies sont détectées, on estime la largeur de chaque voie meg, ,147 d et deux angles de correction du roulis a et d! . Si seule- ment deux voies sont détectées, alors on estime un seul angle de correction du roulis. Si seule la voie centrale est détectée, on estime sa largeur, mais pas de roulis. Plus précisément, l'estimation de la largeur d'une voie (blocs 44, 46) uti- lise les équations des droites des bords des voies, équations données pour chaque image par le module de détection des bords 10. Plus précisément, une voie est définie par deux bords, à gauche et à droite, avec les équations suivantes : D gaud : u= al y +hl droi : u= a2 +1)2 À partir de ces bords, la largeur peut être estimée par : = H ,,cos( )(a où a représente l'angle entre le véhicule et la route et Hv la hauteur de la caméra. Si la hauteur de la caméra par rapport au sol, est connue, cette valeur peut être utilisée utilement pour calculer l'angle de roulis ; dans le cas contraire, il est possible de prendre une valeur arbitraire quelconque, par exemple = 1 m qui sera ensuite éliminée lors des itérations de l'algo- rithme. En effet, les estimations de la largeur de chaque voie 147, ) sont proportionnelles à IL. Mais comme We est estimé à partir de ces mesures, Tre est aussi proportionnel à et a, ne dépend donc pas de IL. Le seul impact de TI, est lors de l'initialisation de We : si TI, est inconnu, il suffit de lui donner une valeur probable et cohérente avec l'initialisation de We , de sorte qu'il n'est pas nécessaire de connaitre la hauteur véritable de la caméra. Comme bl= b2= f tan (oc) , on peut estimer a pour chaque image par : tan (o.)= b+ b2 2 f La formule précédente pour estimer la largeur d'une voie est valable seulement si les bords sont corrigés parfaitement. Lorsqu'il reste une rotation partielle, ce n'est qu'une estimation de la largeur, mais grâce aux itérations cette estimation convergera vers le bon résultat. Ce calcul est valable pour la voie centrale comme pour les voies de droite et de gauche. En ce qui concerne l'estimation du roulis résiduel pour chaque couple de voies, on dispose de la détection de deux voies parallèles, par exemple voie centrale et voie de gauche. Les largeurs des voies leg et 147m ayant été estimées sur chaque image et W représentant la largeur de référence des voies, on démontre que l'on peut calculer le roulis a, par application de l'expression suivante : cos Cette relation est une approximation au premier ordre, valable lorsque a, est petit. est la hauteur de la caméra par rapport au sol, et a qui repré- sente l'angle entre le véhicule et la route est estimé à partir des équations des bords des voies. La même formule est appliquée pour le couple de voies droite et centrale. REVENDICATIONS1. Un procédé de calibration angulaire de la position d'une caméra video embarquée dans un véhicule automobile, par estimation d'une matrice de calibration (Re) définissant une rotation entre un repère lié à la caméra (X,, Y,, Z,) et un repère lié au véhicule (Xv, Yv, 4,), cette rotation étant dé- finie en termes d'angles de tangage (aï), de lacet (ay) et de roulis (ai), procédé comprenant, pendant une phase de déplacement du véhicule sur une route, des étapes de : a) capture d'une séquence d'images video d'une scène de la route, chaque image de la séquence comprenant une pluralité de traces d'intérêt représentatives de marquages au sol de voies de circulation parallèles ; et b) détection de bords (10) dans les images video, pour l'obtention d'une liste de bords (Dk) rectilignes correspondant, dans chaque image, auxdites traces d'intérêt, procédé caractérisé en ce qu'il comprend en outre les étapes successives suivantes : c) correction (12) de la liste de bords par application à chaque bord d'une rotation selon l'estimée de la matrice de calibration ; d) calcul (16) d'angles de tangage et de lacet résiduels par une méthode de détection du point de fuite appliqué à la liste de bords corrigée ob- tenue à l'étape c) ; e) mise à jour (18) de l'estimée de la matrice de calibration par application d'une rotation selon les angles de tangage et de lacet résiduels calculés à l'étape d) ; f) calcul (20 ; 28, 32, 34, 36) d'un angle de roulis résiduel à partir de la liste de bords corrigée obtenue à l'étape c) et de l'estimée de la matrice de calibration obtenue à l'étape e), par : estimation des largeurs respectives d'au moins deux voies de circulation à partir de la liste de bords ; comparaison des largeurs respectives ainsi estimées et calcul d'au moins une différence de largeur entre voies ; et- calcul de l'angle de roulis en fonction de ladite différence de largeur ainsi calculée, les largeurs réelles des voies étant supposées identiques ; g) mise à jour (40) de l'estimée de la matrice de calibration par applica- tion d'une rotation selon l'angle de roulis résiduel calculé à l'étape f) ; et h) itération, sur au moins une séquence d'images, des étapes c) à g) jusqu'à ce que les angles de tangage, de lacet et de roulis résiduels soient inférieurs à un seuil prédéterminé. 2. Le procédé de la revendication 1, dans lequel l'étape f) d'estimation de l'angle de roulis comprend les étapes successives suivantes : f1) correction (28) de la liste de bords par application à chaque bord d'une rotation selon l'estimée de la matrice de calibration et d'une ro- tation selon l'estimée d'une matrice de roulis ; f2) pour chaque image, à partir de la liste de bords obtenue à l'étape f1) et en fonction d'une référence de largeur de voie (We) : - calcul (32) de la largeur de chacune des voies de circulation détectée dans l'image, et - calcul (32) d'un angle de roulis résiduel pour chaque couple de voies de circulation détecté dans l'image ; f3) pour la séquence d'images, calcul (36) d'une largeur moyenne en fonction des largeurs de voie estimées à l'étape f2) pour les différentes images de la séquence ; f4) pour la séquence d'images, calcul (34) d'un angle de roulis résiduel moyen à partir des angles de roulis résiduels calculés à l'étape f2) pour les différentes images de la séquence ; f5) mise à jour (40) : - de la matrice de roulis par application d'une rotation selon l'angle de roulis résiduel moyen calculé à l'étape f4), et - de la largeur de référence moyenne par la largeur moyenne calculée à l'étape f3) ; f6) itération des étapes f1) à f5) jusqu'à ce que l'angle de roulis résiduel soit inférieur à un seuil prédéterminé ; et f7) extraction d'un angle de roulis estimé à partir de la matrice de roulis. 3. Le procédé de la revendication 2, dans lequel l'étape f2) de calcul de la largeur de chacune des voies de circulation et de l'angle de roulis pour chaque couple de voies de circulation comprend : - en cas de détection dans l'image courante de plus de deux voies de circulation, la délivrance en sortie de trois valeurs de largeur de voie et de deux valeurs d'angle de roulis résiduel ; - en cas de détection dans l'image courante de deux voies de circulation, la délivrance en sortie seulement de deux valeurs de largeur de voie et d'une valeur d'angle de roulis ; - en cas de détection dans l'image courante d'une seule voie de circulation, la délivrance en sortie seulement d'une valeur de largeur de voie. 4. Le procédé de la revendication 2, dans lequel à l'étape f2) l'estimation de la largeur de chacune des voies de circulation comprend l'application de la relation : w= cos W étant l'estimée de la largeur de la voie, H, étant la hauteur de la caméra par rapport au sol, a étant l'angle du véhicule par rapport à l'axe de la voie sur lequel il cir- cule, et al et a2 étant les coefficients directeurs des bords délimitant la voie de circulation dans l'image. 5. Le procédé de la revendication 2, dans lequel à l'étape f2) l'estimation de l'angle de roulis comprend l'application de la relation : cos(oc) H, W étant la référence de largeur de voie ou son estimation, 1,:i7g W et étant les estimées de largeur de voies du couple de voies con- sidéré,a étant l'angle du véhicule par rapport à l'axe de la voie sur lequel il circule, et H, étant la hauteur de la caméra par rapport au sol.