FR 3073655 A1 20190517 FR 1701185 A 20171116 I Procédé de détermination d'un vecteur vitesse visuel d'un engin mobile, produit programme d'ordinateur et système de détermination associés La présente invention concerne un procédé de détermination d'un vecteur vitesse visuel d'un engin mobile. La présente invention concerne également un produit programme d'ordinateur et un système de détermination associés. Dans le domaine aéronautique, l'utilisation d'un vecteur vitesse appelé également vecteur de trajectoire de vol ou vecteur FPV (de l'anglais « Flight Path Vector ») pour piloter l'aéronef, est déjà connue. Ce vecteur permet de présenter au pilote de l'aéronef une information sur la trajectoire instantanée de l'appareil et en particulier, de connaitre la dérive (connue sous le terme anglais « drift ») et la pente (connue sous le terme anglais « Flight Path Angle ») de l'aéronef par rapport à un repère géographique. De manière générale, ces informations sont utiles au pilote pour surveiller la pente de l'aéronef ou bien, associées à une représentation de l'environnement extérieur (relief ou piste de destination), pour s'assurer que l'aéronef suit bien la trajectoire prévue. Actuellement, pour calculer le vecteur vitesse, il est nécessaire de calculer des vitesses de l'aéronef dans un repère géographique. Un tel calcul est traditionnellement réalisé à partir des informations fournies par un système de positionnement par satellites tel que le système GPS ou une centrale inertielle embarquée. Il est par ailleurs connu d'utiliser une hybridation des informations fournies par ces deux systèmes afin notamment d'améliorer la qualité du résultat. Toutefois, l'utilisation de ces systèmes se heurte au problème d'intégrité et de précision des informations fournies par ces systèmes. En effet, les vitesses géographiques fournies par une centrale inertielle sont très intègres mais peu précises, alors que celles fournies par le système GPS sont très précises mais ne présentent aucune intégrité. De plus, lors d'un affichage conforme au pilote, notamment via un afficheur tête haute de type HUD (de l'anglais « Head-Up Display »), la précision de l'affichage du vecteur vitesse par rapport à l'environnement extérieur joue un rôle particulièrement important. En effet, une erreur de cap minime (de l'ordre de 1°), peut entraîner une non-conformité de la projection du vecteur vitesse sur l'afficheur tête haute et donc par rapport à l'environnement extérieur, préjudiciable au pilotage. 2 La présente invention a pour but de proposer un procédé de détermination d'un vecteur vitesse visuel présentant une bonne intégrité et une bonne précision, tout en étant particulièrement bien adapté à un affichage conforme. À cet effet, l'invention a pour objet un procédé de détermination d'un vecteur 5 vitesse visuel d'un aéronef se déplaçant dans un milieu environnant, le vecteur vitesse visuel représentant la direction de déplacement de l'engin mobile dans ce milieu. Le procédé comporte les étapes suivantes : - acquisition d'une séquence d'images du milieu environnant à partir d'un capteur d'images orienté selon une direction de prise d'images et définissant un centre optique de 10 projections sensiblement fixe par rapport à l'engin mobile ; - analyse d'au moins deux images successives afin de déterminer sur chacune des deux images un point, dit épipôle, représentant sur cette image la position du centre optique du capteur d'images au moment de la prise de l'autre image ; - pour chaque image analysée, détermination de la position de l'épipôle de cette 15 image sur un afficheur utilisable pour piloter l'engin mobile et affichage sur l'afficheur d'un symbole représentatif du vecteur vitesse visuel dans cette position. Suivant d'autres aspects avantageux de l'invention, le procédé comprend une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : 20 - l'étape d'analyse d'au moins deux images comprend les sous-étapes suivantes : - détermination d'une pluralité de points de repère sur l'une des deux images, dite première image, et d'une pluralité de points de repère sur l'autre image, dite deuxième image ; - association de points de repère entre la première et la deuxième images ; 25 - à partir des points de repère associés, détermination d'une matrice fondamentale définissant pour chaque point de la première image une droite sur la deuxième image, dite droite épipolaire, sur laquelle se trouve le point de la deuxième image correspondant à ce point de la première image ; - détermination de l'épipôle sur la première image en tant que le noyau de la 30 matrice fondamentale ; - détermination de l'épipôle sur la deuxième image en tant que le noyau de la matrice fondamentale transposée ; - le symbole représentatif du vecteur vitesse visuel est affiché sur l'afficheur de manière conforme, ledit afficheur étant de préférence un afficheur tête haute ; 3 - le capteur d'images est monté dans l'engin mobile de sorte que la direction de prise d'images corresponde sensiblement à la direction de déplacement de l'engin mobile ; - le capteur d'images est monté sur la tête d'un pilote de l'engin mobile ; 5 - le procédé comprend en outre une étape de comparaison du vecteur vitesse visuel avec un vecteur vitesse de référence déterminé à partir d'informations fournies par un système de positionnement par satellites et/ou un capteur inertiel ; - le procédé comprend en outre une étape de comparaison du vecteur vitesse visuel avec le vecteur vitesse de référence comprend les sous-étapes suivantes, mise en 10 oeuvre pour chaque couple d'images successives analysées : - détermination d'une première vitesse visuelle correspondant à la vitesse de déplacement du centre optique du capteur d'images entre les images dudit couple ; - détermination d'une deuxième vitesse visuelle correspondant à la première 15 vitesse visuelle calculée dans le repère de l'engin mobile ; - détermination à partir de la deuxième vitesse visuelle, d'une composante de dérive visuelle de l'engin mobile et d'une composante de pente visuelle de l'engin mobile en utilisant le roulis et le tangage l'engin mobile ; - détermination à partir du vecteur vitesse de référence, d'une composante de 20 pente de référence de l'engin mobile et comparaison de cette composante avec la composante de pente visuelle ; - l'étape de comparaison du vecteur vitesse visuel avec le vecteur vitesse de référence qui comprend en outre pour chaque couple d'images successives analysées, la sous-étape de détermination, à partir du vecteur vitesse de référence et du cap de l'engin 25 mobile, d'une composante de dérive de référence et de comparaison de cette composante avec la composante de dérive visuelle ; - l'étape de comparaison du vecteur vitesse visuel avec le vecteur vitesse de référence comprend en outre les sous-étapes suivantes, mise en oeuvre pour chaque couple d'images successives analysées : 30 - détermination d'une troisième vitesse visuelle correspondant à la première vitesse visuelle calculée dans un repère géographique ; - détermination à partir de la troisième vitesse visuelle, d'une composante d'angle de route visuel de l'engin mobile ; - détermination à partir du vecteur vitesse de référence, d'une composante d'angle de route de référence de l'engin mobile et comparaison de cette composante avec la composante d'angle de route visuel ; 4 - l'étape de comparaison du vecteur vitesse visuel avec le vecteur vitesse de référence qui comprend en outre pour chaque couple d'images successives analysées, en cas d'un mouvement en rotation de l'engin mobile, des sous-étapes préliminaires suivantes : 5 - acquisition d'angles de rotation de la direction de prise d'images ; - en fonction des angles de rotation acquises, transformation homographique de l'une des images dudit couple pour aligner la direction de prise d'images correspondant à cette image avec celle correspondant à l'autre image. - l'engin mobile est un aéronef. 10 L'invention a également pour objet un produit programme d'ordinateur comportant des instructions logicielles qui, lorsque mises en oeuvre par un équipement informatique, mettent en oeuvre le procédé tel que défini précédemment. L'invention a également pour objet un système de détermination d'un vecteur vitesse visuel d'un engin mobile, comportant des moyens de mise en oeuvre du procédé 15 tel que décrit plus haut. Ces caractéristiques et avantages de l'invention apparaitront à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif, et faite en référence aux dessins annexés, sur lesquels : - la figure 1 est une vue schématique en perspective d'un engin mobile se 20 déplaçant dans un milieu environnant ; - la figure 2 est une vue schématique de côté de l'engin mobile de la figure 1, l'engin mobile comportant un système de détermination selon l'invention ; - la figure 3 est une vue schématique du système de détermination de la figure 2; - la figure 4 est un organigramme d'un procédé de détermination selon l'invention, 25 le procédé étant mis en oeuvre par le système de détermination de la figure 2; - la figure 5 est vue schématique illustrant la mise en oeuvre de l'une des étapes du procédé de détermination de la figure 4. Dans l'exemple de réalisation de la figure 1, l'engin mobile est un aéronef 10 tel qu'un avion, notamment un avion commercial, piloté par au moins un pilote. 30 Toutefois, dans le cas général, l'engin mobile peut correspondre à tout autre engin se déplaçant sensiblement sur la surface terrestre ou volant au moins dans l'atmosphère terrestre et pilotable par au moins un pilote directement à partir de celui-ci ou à distance. Ainsi, l'engin mobile peut être un véhicule, un avion miliaire, un hélicoptère, un drone, etc. Sur la figure 1, l'aéronef 10 se déplace dans un milieu environnant par rapport à 35 un repère géographique XYZ d'origine O. 5 Ce repère géographique XYZ correspond par exemple au référentiel connu sous le terme anglais « North-East-Down ». Dans ce cas, l'axe OX s'étend suivant la direction nord, l'axe OY suivant la direction est et l'axe OZ suivant la verticale orientée vers la surface terrestre. Le plan formé par les axes OX et OY est appelé plan horizontal. 5 Un repère d'origine 01, dit repère aéronef X1Y1Z1, est par ailleurs associé à l'aéronef 10 de sorte que l'axe 01X1 s'étend dans la direction longitudinale de l'aéronef 10, l'axe 01Y1 se situe dans la plan de l'aéronef 10, est perpendiculaire à l'axe 01X1 et s'étend à gauche de celui-ci, et l'axe 01Z1 complète un repère orthogonal droit direct. Le repère aéronef Xllyi est alors mobile en translation et en rotation par rapport 10 au repère géographique XYZ. L'orientation du repère aéronef Xilyi par rapport au repère géographique XYZ définit l'attitude de l'aéronef 10. Cette attitude est déterminée notamment par l'angle de roulis de l'aéronef 10 correspondant à l'angle entre l'axe 01Y1 et le plan horizontal, et par l'angle de tangage 0 de l'aéronef 10 correspondant à l'angle entre l'axe 01X1 et le 15 plan horizontal. Le vecteur de vitesse V de l'aéronef 10 est défini comme la vitesse de déplacement du point 01 par rapport au point 0. Dans le repère géographique XYZ, le vecteur de vitesse V peut s'écrire sous la forme [VN; VE; VD]t où la composante VN est appelée vitesse nord, la composante VE 20 vitesse est et la composante VD (dénotée également Vz) vitesse verticale. Il est également possible de définir une vitesse horizontale Vh' de l'aéronef 10 correspondant à la projection du vecteur de vitesse V sur le plan horizontal, une vitesse longitudinale 17109 correspondant à la projection du vecteur de vitesse V sur l'axe longitudinal 01X1 et une vitesse latérale Viat correspondant à la projection du vecteur de 25 vitesse V sur l'axe 01Y1. Dans la suite de la description, la pente de l'aéronef 10, appelée également angle FPA (de l'anglais « Flight Path Angle »), est définie comme l'angle entre la vitesse verticale et la vitesse horizontale. Cet angle peut alors être déterminé comme suit : FPA = atan Vz/Vh'. La dérive de l'aéronef 10, connue également sous le terme anglais « drift », est 30 définie comme l'angle entre la vitesse latérale et la vitesse longitudinale, les deux vitesses étant projetées sur le plan horizontal. Cet angle peut alors être déterminé comme suit : Drift = atan Viat/Vi'g. 6 Le cap p de l'aéronef 10, connu également sous le terme anglais « heading », est défini comme l'angle entre l'axe longitudinal 01X1 de l'aéronef projeté sur le plan horizontal et l'axe OX, c'est-à-dire la direction nord. La route de l'aéronef 10, connue également sous le terme anglais « Track », est 5 définie comme la projection du déplacement de l'aéronef sur le plan horizontal. L'angle de route, connue également sous le terme anglais « Track Angle », est défini comme l'angle entre la vitesse est et une vitesse projetée sur un axe de référence, par exemple la vitesse nord lorsqu'elle n'est pas nulle. Cet angle peut être déterminé comme suit : Track = atan VE/VN. 10 Il est à noter que les termes précités sont connus de l'homme du métier et peuvent être définis, le cas échéant, autrement, en fonction notamment du repère géographique choisi. En référence à la figure 2, l'aéronef 10 comporte un système de navigation 12, un afficheur 14, un capteur d'images 16 et un système de détermination 18 d'un vecteur 15 vitesse visuel, selon l'invention. Le système de navigation 12 est connu en soi et est configuré pour fournir au pilote des informations relatives à la position et à la vitesse de l'aéronef 10. Ainsi, par exemple, le système de navigation 12 comprend un récepteur de signaux radioélectriques 21 et un capteur inertiel 22. 20 Le récepteur 21 permet d'acquérir notamment des signaux radioélectriques issus de satellites faisant partie d'un système de positionnement par satellites, tel que le système GPS, et de traiter ces signaux afin de déterminer la vitesse et la position de l'aéronef 10. La vitesse de l'aéronef 10 déterminée par ce récepteur sera dénotée par la suite V Gps. 25 Le capteur inertiel 22 est par exemple une centrale inertielle de basse ou de haute précision apte à fournir les attitudes de l'aéronef 10, c'est-à-dire le roulis 4) et le tangage 0, et le cap IL L'afficheur 14 est utilisable par le pilote pour piloter l'aéronef 10 et présente par exemple un afficheur tête haute de type HUD. 30 Cet afficheur 14 permet notamment de visualiser différents moyens visuels d'aide au pilotage de l'aéronef 10 de manière conforme avec le milieu environnant. À cet effet, l'afficheur 14 comprend par exemple un écran transparent sur lequel sont projetés ces moyens visuels. Selon un exemple de réalisation, l'afficheur 14 est monté dans le casque du pilote. 7 Les moyens visuels d'aide au pilotage affichés sur l'afficheur 14 comprennent notamment un vecteur vitesse connu sous le terme anglais « Flight Path Vector ». De manière connue en soi, ce vecteur permet au pilote de connaitre la dérive et la pente de l'aéronef 10 qui sont déterminées notamment à partir des données fournies par 5 le système de navigation 12. Un tel vecteur vitesse sera dénoté par la suite par le terme « vecteur vitesse de référence » pour le différencier du vecteur vitesse visuel décrit par la suite. Le capteur d'images 16 permet d'acquérir des images du milieu environnant. Pour ce faire, le capteur d'images 16 est monté soit à l'extérieur de l'aéronef 10 soit à l'intérieur 10 de celui-ci dans un endroit donnant une visibilité extérieure, comme par exemple le cockpit. Le capteur d'images 16 est orienté selon une direction de prise d'images qui correspond par exemple sensiblement à la direction de déplacement de l'aéronef 10. Le capteur d'images 16 présente par exemple le même champ de vision que celui 15 du pilote. Avantageusement, le capteur d'images 16 est monté au-dessus de la tête du pilote, par exemple sur son casque. Dans ce dernier cas, le capteur 16 est sensiblement fixe par rapport à l'aéronef 10. Autrement dit, dans ce cas, certaines variations de la position du capteur 16 par rapport à l'aéronef sont possibles. Dans l'exemple décrit ci-dessous, le capteur d'images 16 est une caméra 20 permettant de prendre des images successives du milieu environnant avec une fréquence prédéterminée ou réglable. Ces images sont en couleurs ou monochromes. Toutefois, dans le cas général, le capteur d'images 16 peut être un capteur infrarouge, un radar, un lidar ou tout autre appareil permettant de différencier au moins certains objets du milieu environnant. Il est alors à noter que le terme « image » selon 25 l'invention, correspond à toute représentation du milieu environnant faite par un tel capteur. Selon un exemple de réalisation, le capteur d'images 16 se présente sous la forme de plusieurs capteurs élémentaires disposés éventuellement dans des endroits différents. Dans ce cas, le terme « image » utilisé dans la suite de la description est relatif 30 à une image de synthèse composée des images prises par l'ensemble des capteurs élémentaires. Le capteur d'images 16 définit notamment un centre optique de projection présentant le point d'intersection de toutes les droites de projection, chaque droite de projection présentant la direction suivant laquelle chaque point de l'espace est projeté sur 35 l'image correspondante. 8 Le capteur d'images 16 définit en outre un repère local, dit repère capteur d'origine 0, coïncidant avec le centre optique. La matrice de transformation du repère capteur Xcly, vers le repère aéronef XifiZi sera dénotée par la suite par Rd. 5 Le système de détermination 18 selon l'invention sera désormais expliqué en référence à la figure 3. Ainsi, comme cela est visible sur la figure 3, le système de détermination 18 comprend un module d'acquisition 25, un module de traitement 26 et un module de sortie 27. 10 Le module d'acquisition 25 est raccordé au système de navigation 12 et au capteur d'images 16. Il est ainsi apte à acquérir des informations issues du système de navigation 12 et les images prises par le capteur d'images 16. Le module de traitement 26 permet de traiter des données acquises par le module 15 25 afin de mettre en oeuvre un procédé de détermination selon l'invention expliqué en détail par la suite. En particulier, le module de traitement 26 permet de déterminer un vecteur vitesse visuel à partir d'au moins deux images prises par le capteur 16. Finalement, le module de sortie 27 est raccordé à l'afficheur 14 et permet de traiter 20 des données produites par le module de traitement 26 afin de les afficher sur l'afficheur 14, comme cela sera expliqué par la suite. Le système de détermination 18 se présente par exemple sous la forme d'un calculateur muni d'un processeur apte à exécuter des logiciels stockés dans une mémoire de celui-ci. Dans ce cas, les modules 26 à 27 se présentent au moins partiellement sous 25 la forme des logiciels exécutables par ce processeur. Le système de détermination 18 peut comprendre en outre des parties matérielles, comme par exemple des circuits logiques programmables, pour mettre en oeuvre au moins certaines étapes ou sous-étapes du procédé de détermination. Le procédé de détermination du vecteur vitesse visuel selon l'invention sera 30 désormais expliqué en référence à la figure 4 présentant un organigramme de ses étapes. Ce procédé comprend notamment deux phases qui peuvent être mises en oeuvre de manière indépendante ou bien, l'une après l'autre. Une première phase P1 du procédé permet de déterminer le vecteur vitesse visuel afin qu'il puisse être affiché sur l'afficheur 14 par exemple à la place ou en complément du 35 vecteur vitesse de référence. 9 Une deuxième phase P2 permet de déterminer les composantes de ce vecteur, à savoir la dérive et la pente, afin qu'elles puissent être comparées avec celles du vecteur vitesse de référence. Lors d'une phase initiale 110 de la première phase P1, le module d'acquisition 25 5 acquiert une pluralité d'images prises par le capteur 16 et transmet ces images au module de traitement 26. Lors de la phase 120 suivante, le module de traitement 26 analyse au moins deux images successives afin de déterminer sur chacune d'entre elles un point, dit épipôle. En particulier, l'épipôle sur chacune des deux images représente sur cette image 10 la position du centre optique 0, du capteur d'images 16 au moment de la prise de l'autre image. Les épipôles el et e2 des deux images successives I et /2 prises respectivement aux moments temporels t1 et t2 (t1 étant antérieur à t2), sont illustrés sur la figure 5. Sur cette figure 5, le point 01.- correspond au point optique 0, du capteur d'images 15 16 au moment t1 et le point O correspond au point optique 0, du capteur d'images 16 au moment t2. Ainsi, les points d'intersection de la droite OM avec les plans formant les images correspondent aux épipôles el et e2 de ces images. En outre, de manière connue en soi, l'épipôle de chaque image est le point d'intersection de toutes les droites épipolaires de cette image. 20 En particulier, dans l'exemple de la figure 5, un point Q quelconque de l'espace ayant le point q1 comme la projection sur l'image /1 peut se projeter à n'importe quel point de la droite q2 de l'image /2 . La droite q2 est alors appelée droite épipolaire. Pour déterminer les épipôles sur chacune des deux images, le module de traitement met en oeuvre les sous-étapes 121 à 125 décrites ci-dessous. 25 Lors de la sous-étape 121, le module de traitement 26 détermine une pluralité de points de repère, appelés également amers, sur chacune des deux images. En particulier, par un point de repère, on entend un point remarquable de l'espace qui peut être reconnu sur les deux images. Pour ce faire, il est possible d'associer par exemple à chaque point de repère déterminé sur l'une des images, un descripteur 30 permettant de le reconnaitre sur l'autre image. Puis, lors de la sous-étape 122, le module de traitement 26 associe des points de repère entre les deux images. Cette association est par exemple effectuée en analysant les descripteurs de ces points et/ou en analysant des valeurs de corrélation dans ces points, selon des méthodes connues en soi. 35 Puis, lors de la sous-étape 123, à partir des points de repère associés, le module de traitement 26 détermine une matrice fondamentale F. 10 Cette matrice fondamentale F est telle que pour chaque point de repère x' de l'une des deux images, dite première image, et le point de repère x lui associé sur l'autre image, dite deuxième image, la matrice fondamentale F satisfait la relation suivante : x't. F. x = 0. Autrement dit, la matrice fondamentale F définit pour chaque point de la première 5 image une droite épipolaire sur la deuxième image. Puis, lors de la sous-étape 124, le module de traitement 26 détermine l'épipôle el sur la première image en tant que le noyau de la matrice fondamentale F. Autrement dit, l'épipôle el est déterminé à partir de la relation suivante : F.e1 = 0. Puis, lors de la sous-étape 125, le module de traitement 26 détermine l'épipôle e2 10 sur la deuxième image en tant que le noyau de la matrice fondamentale transposée Ft. Autrement dit, l'épipôle e2 est déterminé à partir de la relation suivante : Ft. e2 = 0. Finalement, le module de traitement 26 transmet les positions des épipôles déterminés au module de sortie 27. Lors de l'étape 130 finale de la première phase P1, le module de sortie 27 15 détermine sur l'afficheur 14, la position conforme de l'épipôle de l'image correspondant au moment courant. L'afficheur 14 affiche alors dans cette position un symbole représentatif du vecteur vitesse visuel. Ce symbole est par exemple analogue à celui du vecteur vitesse de référence et est affiché à la place de celui-ci au moins pendant certaines phases de vol de l'aéronef 10. 20 La deuxième phase P2 comprend notamment une étape 140 lors de laquelle le module de traitement 26 compare le vecteur vitesse visuel avec le vecteur vitesse de référence. Cette étape comprend plusieurs sous-étapes décrites ci-dessous qui sont par exemple mises en oeuvre pour le même couple d'images que celles mentionnées 25 précédemment. Lors de la première sous-étape 141, le module de traitement 26 analyse des données fournies par le capteur inertiel 22 pour détecter une rotation de l'aéronef 10 entre les moments de prise des deux images. Si c'est le cas, le module de traitement 26 met en oeuvre les étapes 142 et 143. 30 Dans le cas contraire, le module de traitement 26 passe directement à l'étape 145. Lors des étapes 142 et 143, le module de traitement 26 analyse les deux images pour supprimer l'effet de rotation sur l'une d'entre elles et pour garder donc seulement un mouvement en translation entre ces deux images. 11 Pour ce faire, lors de la sous-étape 142, le module de traitement 26 acquiert les angles de rotation de la direction de prise d'images entre les deux moments temporels. Ces angles sont fournis soit directement par le capteur inertiel 22 et correspondent aux angles de rotation relative entre les moments de prise des deux images, soit sont 5 prétraités par module d'acquisition 25 en fonction des positions respectives du centre optique du capteur d'images 16 et du capteur inertiel 22. Puis, lors de la sous-étape 143, le module de traitement 26 fait une transformation homographique par exemple de la première image pour aligner la direction de prise d'images correspondant à cette image avec celle correspondant à la deuxième image. 10 Cette transformation est effectuée en fonction des angles de rotation acquis selon par exemple des méthodes connues en soi. Lors de l'étape 145, le module de traitement 26 analyse un couple d'images correspondant aux deux images d'origine lorsqu'aucun mouvement de rotation de l'aéronef 10 n'a été détecté lors de l'étape 141, ou un couple d'images avec l'une des 15 images transformée lors des étapes 142 et 143, dans le cas contraire. En particulier, lors de la sous-étape 145, le module de traitement 26 détermine une première vitesse visuelle V' correspondant à la vitesse de déplacement du centre optique du capteur d'images 16 entre les images du couple. Cette vitesse visuelle V' est calculée dans le repère capteur X,Y,Z, en supposant 20 alors que seul un mouvement de translation du capteur a eu lieu entre les prises des deux images. Il est à noter en outre que, compte tenu de la définition du repère capteur X,Y,Z,, la vitesse visuelle V' est déterminée à un facteur d'échelle près tout en tenant compte de la position de l'épipôle sur au moins l'une des deux images. 25 Ainsi, la première vitesse visuelle V' peut être définie à partir de la relation suivante : telx - )co u V' = el)/ - Yo , v 1 où : eixet el), sont les coordonnées de l'épipôle sur la première image ; u et y sont les distances focales en pixels du capteur 16; et 30 xo et yo sont les coordonnées de la position du centre optique sur la première image. 12 Puis, lors de la sous-étape 146, le module de traitement 26 détermine une deuxième vitesse visuelle Vt,1 correspondant à la première vitesse visuelle V' calculée dans le repère aéronef XiYiZi, en utilisant par exemple la relation suivante : Vv1 = Rcilivc- Puis lors de la sous-étape 147, le module de traitement 26 détermine à partir de la 5 deuxième vitesse visuelle V,1, une composante de dérive visuelle Drift, de l'aéronef 10 et une composante de pente visuelle FPA, de l'aéronef en utilisant le roulis et le tangage 6, de l'aéronef 10 qui sont alors fournies par le capteur inertiel 22. En particulier, pour calculer la pente visuelle FPA,, le module de traitement 26 détermine la vitesse horizontale visuelle Vh'v et la vitesse verticale visuelle 17, de 10 l'aéronef 10 dans le repère géographique XYZ. La pente visuelle FPA, est alors déterminée selon l'expression suivante : FPA, = atan17,,/Viwn, . Pour calculer la dérive visuelle Drift, le module de traitement 26 détermine la vitesse latérale visuelle Via, et la vitesse longitudinale visuelle 17/(mg, de l'aéronef 10 dans le repère géographique XYZ. La dérive visuelle Drift, est alors déterminée selon 15 l'expression suivante : Drift, = atan Viatv/Vi'g, . Puis, lors de la sous-étape 148, le module de traitement 26 détermine une composante de pente de référence FPAref de l'aéronef 10 et compare cette composante avec la composante de pente visuelle FPA,. La pente de référence FPAref est par exemple déterminée à partir du vecteur 20 V Gps, c'est-à-dire à partir du vecteur de vitesse de l'aéronef déterminé via le système GPS. Pour ce faire, le module de traitement 26 détermine la vitesse horizontale de référence V horre f à partir de la vitesse nord VN et la vitesse est VE du vecteur V Gps ainsi que la vitesse verticale V, qui est donnée directement par le vecteur V Gps. Ainsi, la pente 25 de référence FPAref est déterminée selon l'expression suivante : FPAref = atan liz/Vhorref Puis, lors de la même sous-étape 148, le module de traitement 26 compare la pente de référence FPAref avec la pente visuelle FPA, et en fonction du résultat de cette comparaison, associe par exemple un niveau de confiance à l'une ou à l'autre de ces mesures. Ce niveau de confiance est par exemple communiqué au pilote ou à un système 30 avionique pouvant exploiter ce résultat. 13 Lors de la sous-étape 149 suivante, le module de traitement 26 détermine une composante de dérive de référence Drift'f et compare cette composante avec la composante de dérive visuelle Drift. Pour ce faire, le module de traitement 26 acquiert d'abord via le module 25, le cap 5 lp de l'aéronef fourni par exemple par le capteur inertiel 21. En connaissant le cap 1/), le module de traitement 26 détermine la vitesse latérale de référence Viat'f et la vitesse longitudinale de référence Viong'f de l'aéronef 10 à partir par exemple du vecteur V Gps. La dérive de référence Drift'f est ainsi calculée selon l'expression suivante : Drif t'f =- atan Viatref. 10 À la fin de cette sous-étape 149, le module de traitement 26 compare la dérive de référence Driftref avec la dérive visuelle Draft, et comme dans le cas précédent, en fonction du résultat de cette comparaison, associe par exemple un niveau de confiance à l'une ou à l'autre de ces mesures. De manière analogue au cas précédent, ce niveau de confiance est par exemple communiqué au pilote ou est exploité par un système 15 avionique. En variante, lorsque la mesure du cap 1,/, n'est pas disponible ou n'est pas suffisamment précise, au lieu de calculer la dérive de référence Drift'f, le module de traitement 26 calcule lors de cette sous-étape 149, une composante d'angle de route de référence de l'aéronef directement à partir du vecteur V Gps. L'angle de route de référence 20 est déterminé à partir de la relation suivante : Track'f = atan VE/VN. Puis lors de la même sous-étape, le module de traitement 26 compare l'angle de route de référence Track'f avec un angle de route visuel Tracki, calculé à partir du vecteur de vitesse visuel. Pour déterminer l'angle de route visuel Track', le module de traitement 26 25 détermine d'abord une troisième vitesse visuelle V, correspondant à la première vitesse visuelle Vvc calculée dans le repère géographique XYZ. Cette troisième vitesse visuelle V, est par exemple calculée en utilisant la relation suivante : V' = R,V',, où R, est la matrice de passage du repère caméra X,Y,Z, vers le repère 30 géographique XYZ. Puis, l'angle de route visuel Track, est déterminé à partir de la relation suivante : Track'f = atan VEv/VN', 14 où VE, et VNv sont respectivement la vitesse est et la vitesse nord faisant partie du vecteur de la troisième vitesse visuelle V,. La matrice R, peut par exemple être déterminée à l'initialisation du procédé en comparant les projections du vecteur de vitesse V Gps et du premier vecteur de vitesse 5 visuelle À la fin de la sous-étape 149, le module de traitement 26 compare l'angle de route visuel Track, avec l'angle de route de référence Track'f de manière analogue à celle expliquée précédemment. Bien entendu, l'algorithme de détermination de l'angle de route visuel Track, peut 10 être également appliqué pour déterminer la pente visuelle lors de la sous-étape 147. On conçoit alors que l'invention présente un certain nombre d'avantages. Tout d'abord, le procédé selon l'invention permet de déterminer un vecteur vitesse comparable au vecteur vitesse classiquement utilisé dans le domaine aéronautique. Ce vecteur vitesse visuel présente un bon équilibre entre l'intégrité et la précision 15 surtout lors du décollage et de l'atterrissage de l'aéronef. En effet, lors de ces phases, le milieu environnant de l'aéronef présente plus de repères visuels extérieurs que par exemple la phase de croisière. De plus, étant donné que lors de ces phases l'aéronef à une vitesse réduite, l'utilisation du vecteur vitesse visuel est particulièrement avantageuse par rapport au vecteur vitesse classique. 20 En outre, le vecteur vitesse visuel inclut intrinsèquement les rotations instantanées de l'aéronef et peut être ainsi utilisé davantage pour piloter l'aéronef par rapport au vecteur vitesse classique. En effet, il est connu que le vecteur vitesse classique prend en compte les rotations de l'aéronef avec un certain retard ce qui rend généralement impossible d'utiliser 25 uniquement ce vecteur pour piloter l'aéronef. En outre, il est possible de calculer les composantes exactes, à savoir la dérive et la pente, du vecteur de vitesse visuel pour les comparer par exemple avec les composantes analogues du vecteur vitesse classique. La dérive et la pente de l'aéronef calculées à partir du vecteur vitesse visuel sont 30 particulièrement précises même lorsqu'un capteur d'inertiel de basse précision est utilisé. Finalement, l'utilisation du vecteur vitesse visuel est particulièrement avantageuse lorsqu'un afficheur tête haute à affichage conforme est utilisé. En effet, le vecteur vitesse visuel est déterminé directement à partir des images conformes du milieu environnant ce qui permet de le positionner de manière particulièrement précise sur l'afficheur REVENDICATIONS1.- Procédé de détermination d'un vecteur vitesse visuel d'un aéronef (10) se déplaçant dans un milieu environnant, le vecteur vitesse visuel représentant la direction de déplacement de l'engin mobile (10) dans ce milieu ; le procédé comportant les étapes suivantes : - acquisition (110) d'une séquence d'images du milieu environnant à partir d'un capteur d'images (16) orienté selon une direction de prise d'images et définissant un centre optique de projections (Or) sensiblement fixe par rapport à l'engin mobile (10) ; - analyse (120) d'au moins deux images successives afin de déterminer sur chacune des deux images un point (e1, e2), dit épipôle, représentant sur cette image la position du centre optique (Or) du capteur d'images (16) au moment de la prise de l'autre image ; - pour chaque image analysée, détermination (130) de la position de l'épipôle e2) de cette image sur un afficheur (14) utilisable pour piloter l'engin mobile (10) et affichage (130) sur l'afficheur (14) d'un symbole représentatif du vecteur vitesse visuel dans cette position. 2.- Procédé selon la revendication 1, dans lequel l'étape d'analyse (120) d'au moins deux images comprend les sous-étapes suivantes : - détermination (121) d'une pluralité de points de repère sur l'une des deux images, dite première image, et d'une pluralité de points de repère sur l'autre image, dite deuxième image ; - association (122) de points de repère entre la première et la deuxième images ; - à partir des points de repère associés, détermination (123) d'une matrice fondamentale (F) définissant pour chaque point de la première image une droite sur la deuxième image, dite droite épipolaire, sur laquelle se trouve le point de la deuxième image correspondant à ce point de la première image ; - détermination (124) de l'épipôle e2) sur la première image en tant que le noyau de la matrice fondamentale (F) ; - détermination (125) de l'épipôle e2) sur la deuxième image en tant que le noyau de la matrice fondamentale transposée (Ft). 3.- Procédé selon la revendication 1 ou 2, dans lequel le symbole représentatif du vecteur vitesse visuel est affiché sur l'afficheur (14) de manière conforme, ledit afficheur (14) étant de préférence un afficheur tête haute. 16 4.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le capteur d'images (16) est monté dans l'engin mobile (10) de sorte que la direction de prise d'images corresponde sensiblement à la direction de déplacement de l'engin mobile 5 (10). 5.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le capteur d'images (16) est monté sur la tête d'un pilote de l'engin mobile (10). 10 6.- Procédé selon l'une quelconque des revendications précédentes, comprenant en outre une étape de comparaison (140) du vecteur vitesse visuel avec un vecteur vitesse de référence déterminé à partir d'informations fournies par un système (21) de positionnement par satellites et/ou un capteur inertiel (22). 15 7.- Procédé selon la revendication 6, dans lequel l'étape de comparaison (140) du vecteur vitesse visuel avec le vecteur vitesse de référence comprend les sous-étapes suivantes, mise en oeuvre pour chaque couple d'images successives analysées : - détermination (145) d'une première vitesse visuelle (V') correspondant à la vitesse de déplacement du centre optique (Or) du capteur d'images (16) entre les images 20 dudit couple ; - détermination (146) d'une deuxième vitesse visuelle (17,1) correspondant à la première vitesse visuelle (V') calculée dans le repère (XiYiZi) de l'engin mobile (10) ; - détermination (147) à partir de la deuxième vitesse visuelle (V1), d'une composante de dérive visuelle (Drift) de l'engin mobile (10) et d'une composante de 25 pente visuelle (FPA,) de l'engin mobile (10) en utilisant le roulis (0) et le tangage (4)) l'engin mobile (10) ; - détermination (148) à partir du vecteur vitesse de référence, d'une composante de pente de référence (FPA'f) de l'engin mobile (10) et comparaison (148) de cette composante avec la composante de pente visuelle (FPA,). 30 8.- Procédé selon la revendication 7, dans lequel l'étape de comparaison (140) du vecteur vitesse visuel avec le vecteur vitesse de référence comprend en outre pour chaque couple d'images successives analysées, la sous-étape (149) de détermination, à partir du vecteur vitesse de référence et du cap (e) de l'engin mobile (10), d'une 35 composante de dérive de référence (Drift'f) et de comparaison (149) de cette composante avec la composante de dérive visuelle (Drift). 17 9.- Procédé selon la revendication 7 ou 8, dans lequel l'étape de comparaison (140) du vecteur vitesse visuel avec le vecteur vitesse de référence comprend en outre les sous-étapes suivantes, mise en oeuvre pour chaque couple d'images successives 5 analysées : - détermination (149) d'une troisième vitesse visuelle (Vu) correspondant à la première vitesse visuelle (V') calculée dans un repère géographique (XYZ); - détermination (149) à partir de la troisième vitesse visuelle (Vu), d'une composante d'angle de route visuel (Tracky) de l'engin mobile (10) ; 10 - détermination (149) à partir du vecteur vitesse de référence, d'une composante d'angle de route de référence (Track'f) de l'engin mobile (10) et comparaison (149) de cette composante avec la composante d'angle de route visuel (Trackv). 10.- Procédé selon l'une quelconque des revendications précédentes, dans lequel 15 l'étape de comparaison (140) du vecteur vitesse visuel avec le vecteur vitesse de référence comprend en outre pour chaque couple d'images successives analysées, en cas d'un mouvement en rotation de l'engin mobile (10), des sous-étapes préliminaires suivantes : - acquisition (142) d'angles de rotation de la direction de prise d'images ; 20 - en fonction des angles de rotation acquises, transformation (143) homographique de l'une des images dudit couple pour aligner la direction de prise d'images correspondant à cette image avec celle correspondant à l'autre image. 11.- Procédé selon l'une quelconque des revendications précédentes, dans lequel 25 l'engin mobile (10) est un aéronef (10). 12.- Produit programme d'ordinateur comportant des instructions logicielles qui, lorsque mises en oeuvre par un équipement informatique, mettent en oeuvre le procédé selon l'une quelconque des revendications précédentes. 30 13.- Système de détermination (18) d'un vecteur vitesse visuel d'un engin mobile (10), comportant des moyens de mise en oeuvre du procédé selon l'une quelconque des revendications précédentes.