FR 3078398 A1 20190830 FR 1851705 A 20180227 DOMAINE TECHNIQUE AUQUEL SE RAPPORTE L'INVENTION La présente invention concerne de manière générale le domaine de la cartographie. Elle concerne plus particulièrement un procédé d'estimation d'une position précise d'un véhicule sur une carte, comportant : - une étape d'acquisition d'au moins une position géolocalisée du véhicule au moyen d'un système de géolocalisation, - une étape de pré-positionnement du véhicule sur la carte, à sa position géolocalisée, et - une étape de filtrage particulaire au cours de laquelle des positions possibles du véhicule, appelées particules, sont traitées pour déterminer la position précise du véhicule sur la carte. Elle concerne également un véhicule comportant : - des moyens de mémorisation d'une carte, - un système de géolocalisation, et - un calculateur adapté à pré-positionner le véhicule sur la carte. ARRIERE-PLAN TECHNOLOGIQUE Pour assurer la sécurité des véhicules autonomes et des véhicules partiellement automatisés, il est nécessaire d'avoir une connaissance poussée de l'environnement dans lequel ces véhicules évoluent. En pratique, la perception par un véhicule de son environnement se fait de deux façons différentes, à savoir : - en utilisant une carte et un système de géolocalisation du véhicule, et - en utilisant des capteurs extéroceptifs (caméras, capteur RADAR ou LIDAR...). Les sociétés qui élaborent les cartes travaillent actuellement sur des cartes dites « à haute définition », permettant d'obtenir des informations très détaillées sur les caractéristiques du réseau routier (largeur des voies, marquages au sol, panneaux de signalisation...). Ces cartes sont embarquées dans des véhicules équipés de systèmes de géolocalisation, ce qui permet à ces véhicules de se situer sur la carte en une position estimée par une longitude et une latitude. 2 On constate malheureusement que cette position n'est pas toujours très précise et fiable, ce qui se traduit alors par la localisation du véhicule hors de l'itinéraire réellement emprunté. Ce problème peut s'avérer particulièrement dangereux dans le cas d'un véhicule autonome qui utilise ces informations pour se 5 diriger. Pour remédier à ce problème, une solution connue du document DE102013217060 consiste à prendre en compte les informations provenant des capteurs extéroceptifs du véhicule pour ajuster la position du véhicule sur la carte. Cette solution utilise une méthode appelée « filtre de particules ». Selon 10 cette méthode, des particules correspondant aux positions probables du véhicule sont envisagées et traitées de façon à tenter de trouver la position précise du véhicule. Dans le document précité, les particules sont progressivement sélectionnées par ce filtre en fonction des positions des marquages de voies vus par les caméras embarqués dans le véhicule, jusqu'à ce qu'une seule et unique 15 particule correspondant probablement à la position réelle du véhicule soit trouvée. Cette solution présente deux inconvénients. Son premier inconvénient est que sa fiabilité dépend très fortement de la visibilité des marquages de route. On comprend en effet qu'en l'absence de visibilité ou en l'absence de marquage de route, il donne des résultats peu fiables. 20 Son second inconvénient est qu'il ne permet de sélectionner in fine qu'une seule et unique particule, si bien que s'il commet une erreur, le véhicule automobile n'est pas en mesure de le savoir, ce qui peut s'avérer très dangereux. OBJET DE L'INVENTION Afin de remédier aux inconvénients précités de l'état de la technique, la 25 présente invention propose un nouveau procédé permettant de déterminer s'il est possible ou non de connaître la position précise d'un véhicule automobile sur une carte, et, si c'est le cas, quelle est cette position précise. Plus particulièrement, on propose selon l'invention un procédé tel que défini dans l'introduction, dans lequel l'étape de filtrage particulaire comprend : 30 - une étape préalable de répartition de particules sur la carte, en fonction de la position géolocalisée du véhicule sur la carte, puis une étape de mise à jour des particules sur la carte, - une étape de calcul de la vraisemblance de chaque particule, à partir au moins de données issues de la carte, 3 - une étape de sélection d'un ensemble restreint de particules, et - si un indicateur, relatif à la vraisemblance des particules sélectionnées et au nombre de particules sélectionnées, passe en dessous d'un seuil, une étape de ré-échantillonnage de particules sur la carte. 5 Ainsi, grâce à l'invention, le calcul de la vraisemblance de chaque particule est réalisé à partir de données issues de la carte. Ce calcul peut donc être indépendant des données relevées par les capteurs extéroceptifs du véhicule. Il peut donc être indépendant de la visibilité des marquages au sol et de la visibilité de l'environnement du véhicule. Il fournit ainsi des résultats dont la 10 fiabilité n'est pas dépendante des conditions extérieures. Egalement grâce à l'invention, l'étape de ré-échantillonnage n'est pas effectuée systématiquement, ce qui évite un appauvrissement des particules, lequel pourrait potentiellement conduire à éliminer une particule correcte. D'autres caractéristiques avantageuses et non limitatives du procédé 15 conforme à l'invention sont les suivantes : - la vraisemblance de chaque particule est calculée en fonction uniquement de données issues de la carte ; - la vraisemblance de chaque particule est calculée en fonction également de données issues de capteurs permettant au véhicule de percevoir 20 son environnement, à condition que ces données soient jugées fiables ; - à l'étape de filtrage particulaire, la position précise est choisie parmi l'ensemble restreint de particules sélectionnées, en fonction de la vraisemblable de chaque particule ; - la carte stockant des données relatives à des tronçons de routes, la 25 vraisemblance de chaque particule est calculée en fonction de la position du plus proche tronçon de route par rapport à la particule ; - à l'étape préalable, les particules sont réparties dans un disque centré sur la position géolocalisée du véhicule ; - le rayon du disque est déterminé en fonction du niveau de protection 30 horizontal affecté à la position géolocalisée du véhicule ; - à l'étape de sélection, les particules sont sélectionnées en fonction de la distance les séparant de la position géolocalisée du véhicule ; - à l'étape de mise à jour, les particules sont déplacées sur la carte en fonction uniquement d'informations relatives à la dynamique du véhicule ; 4 - à l'étape de ré-échantillonnage, les particules sont ré-échantillonnées à l'aide d'une technique à faible variance. L'invention concerne aussi un véhicule comportant : - des moyens de mémorisation d'une carte, 5 - un système de géolocalisation, et - un calculateur adapté à mettre en oeuvre un procédé d'estimation d'une position précise du véhicule sur la carte tel que précité. DESCRIPTION DETAILLEE D'UN EXEMPLE DE REALISATION La description qui va suivre en regard des dessins annexés, donnés à 10 titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. Sur les dessins annexés : - la figure 1 est un schéma illustrant les différentes étapes d'un procédé conforme à l'invention, 15 - la figure 2 est une vue de dessus d'un véhicule roulant sur une route, - la figure 3 est une vue schématique de particules réparties sur une carte, - la figure 4 est une vue schématique de deux particules situées à côté de deux tronçons de route successifs, et 20 - la figure 5 est une vue schématique de quatre particules situées à côté de quatre tronçons de route. Sur la figure 2, on a représenté un véhicule automobile 10 qui se présente sous la forme d'un voiture et qui roule sur une portion de route à quatre voies de circulation V1, V2, V3, V4. 25 Dans la suite de la description, on s'intéressera plus particulièrement à la localisation de ce véhicule automobile 10 sur une carte mais l'invention ne se limitera pas à un tel exemple. Elle s'appliquera ainsi notamment à la localisation de tout véhicule terrestre, maritime, aérien ou spatial sur une carte. Le véhicule automobile 10 ici considéré comporte classiquement un 30 châssis, un groupe motopropulseur, un système de direction, un système de freinage, et une unité de calcul électronique et/ou informatique, appelée ci-après calculateur. Le calculateur est connecté à des capteurs dits « proprioceptifs », qui permettent de mesurer précisément la vitesse du véhicule et la vitesse angulaire 5 en lacet du véhicule. Le calculateur est de préférence aussi connecté à des capteurs dits « extéroceptifs », qui permettent de percevoir l'environnement immédiat du véhicule automobile 10 (il peut s'agir de caméras, de capteurs RADAR, de 5 capteurs L I DA R ...). Le calculateur est par ailleurs connecté à un système de géolocalisation qui permet d'évaluer la position géolocalisée Po du véhicule 10, ici définie par une latitude et une longitude. Il peut par exemple s'agir d'un système GPS. On considérera ici que ce système de géolocalisation est par ailleurs 10 adapté à transmettre au calculateur une donnée appelée « niveau de protection horizontal HPL » (de l'anglais « horizontal protection level »). Cette donnée, bien connue de l'homme de l'art, correspond à l'incertitude de mesure de la position géolocalisée Po. Sa valeur varie par exemple en fonction du nombre de satellites desquels le système de géolocalisation reçoit des données, de la qualité de la 15 réception des signaux, de la qualité du système de géolocalisation utilisé... Dans le même ordre d'idées, on considérera ici aussi que ce système de géolocalisation est adapté à transmettre au calculateur une matrice de covariance relative à cette même incertitude. Le véhicule automobile 10 considéré pourrait être semi automatisé, de 20 façon à ce que son calculateur puisse par exemple déclencher un freinage d'urgence lorsque le conducteur n'a pas perçu un danger et n'a pas lui-même exercé d'action adaptée. Toutefois, dans la suite de cet exposé, on considérera que le véhicule automobile 10 est du type autonome, et que le calculateur est adapté à 25 commander le groupe motopropulseur, le système de direction, et le système de freinage du véhicule. Le calculateur comporte alors une mémoire informatique qui enregistre des données utilisées dans le cadre du pilotage autonome du véhicule, et notamment dans le cadre du procédé décrit ci-dessous. 30 Il mémorise notamment une application informatique, constituée de programmes d'ordinateur comprenant des instructions dont l'exécution par un processeur permet la mise en oeuvre par le calculateur du procédé décrit ci-après. Il mémorise également une carte topographique dite « à haute-définition ». 6 Cette carte stocke de nombreuses données. Elle comprend tout d'abord des informations relatives à la topographie des routes. Cette topographie est ici mémorisée sous forme de tronçons de routes (ou « mailles » ou, en anglais, « link »). Chaque tronçon de route est ici défini 5 comme une portion d'une seule voie de circulation d'une route, dont les caractéristiques sont constantes sur toute sa longueur (forme des marquages au sol identique le long du tronçon de route, largeur constante de ce tronçon de route...). La carte stocke également d'autres données caractérisant chaque 10 tronçon de route, parmi lesquelles la largeur de la voie de circulation, la forme des marquages au sol situés de part et d'autre de la voie de circulation, la position et la forme de chaque panneau bordant la route au niveau du tronçon de route, les identifiants des tronçons de route précédent et suivant... Le procédé mis en oeuvre par le calculateur pour estimer la position 15 précise Pp du véhicule automobile 10 sur la carte comprend deux grandes opérations, dont une opération de filtrage de particules 100, et une opération de sélection d'hypothèses 200 (voir figure 1). L'opération de sélection d'hypothèses 200 utilise les résultats de l'opération de filtrage de particules 100, si bien qu'elle est mise en oeuvre après 20 celle-ci. On peut alors commencer par décrire la première opération de filtrage de particules 100. Cette opération est mise en oeuvre de manière récursive, c'est-à-dire en boucle et à pas de temps réguliers. 25 Elle comprend trois étapes principales. La première étape 101 consiste pour le calculateur à acquérir différentes données via les capteurs auxquels il est connecté. Le calculateur acquiert ainsi la position géolocalisée Po du véhicule automobile 10 et le niveau de protection horizontal HPL qui lui est associé. Ces 30 données sont acquises grâce au système de géolocalisation qui fournit une latitude, une longitude et un niveau de protection horizontal HPL. Le calculateur acquiert également des données relatives à la dynamique du véhicule automobile 10. Il acquiert ainsi la vitesse V du véhicule et sa vitesse angulaire en lacet Y. 7 La seconde étape 102 est une étape de pré-positionnement du véhicule 10 sur la carte, à la position géolocalisée Po acquise. La troisième étape 103 est une étape de filtrage particulaire au cours de laquelle des positions possibles du véhicule (ou plus exactement des postures 5 possibles du véhicule), appelées particules Ph sont traitées en vue de déterminer la position précise Pp du véhicule 10 sur la carte (ou plus exactement la posture précise du véhicule sur la carte). Chaque particule Pi peut être définie par : - deux coordonnées xi, yi permettant de définir la position de la particule 10 dans un repère cartésien (ces coordonnées sont liées aux latitude et longitude acquises), - un angle de lacet permettant de définir l'angle que fait la particule par rapport à une direction donnée telle que le Nord, et - un identifiant de la carte à laquelle la particule Pi est associée. 15 Sur les figures 3 à 5, on a représenté des particules Pi sous la forme de triangles isocèles, chaque triangle présentant un centre Mi qui correspond à la position de la particule sur la carte, et une orientation qui correspond à l'angle de lacet de la particule sur la carte. Comme le montre la figure 1, la troisième étape 103 de filtrage 20 particulaire est plus précisément constituée de plusieurs sous-étapes que l'on peut maintenant décrire plus en détail. La première sous-étape 110 consiste à déterminer si l'on se trouve ou on dans une phase d'initialisation du filtre particulaire, ce qui est par exemple le cas au démarrage du véhicule automobile 10. 25 On peut tout d'abord se placer dans cette situation, auquel cas aucune particule n'a encore été générée. La sous-étape suivante 112 consiste alors à créer et répartir des particules Pi sur la carte, compte tenu de la position géolocalisée Po du véhicule 10. 30 Pour cela, les particules Pi sont réparties dans un disque centré sur la position géolocalisée Po du véhicule 10, dont le rayon est ici égal au niveau de protection horizontal HPL. Elles sont plus précisément réparties selon une spirale, avec un écart angulaire constant. Les caractéristiques de la spirale et l'écart angulaire entre les 8 particules Pi sont choisis en fonction du nombre de particules P, que l'on souhaite générer. Ce nombre est supérieur à 100, et préférentiellement de l'ordre de 1000. Il est déterminé de façon à obtenir une précision suffisante, sans pour autant 5 surcharger le calculateur. A ce stade, les particules Pi ne sont pas encore orientées. Chaque particule Pi correspond ainsi à une position possible que le véhicule pourrait présenter, compte tenu de l'erreur affectant le système de géolocalisation. 10 Certaines particules, comme on peut le constater sur la figure 3, se situent en dehors de la route. Cela illustre le fait que les particules ne sont pas contraintes sur la carte et qu'elles pourront évoluer dans un espace à deux dimensions. Le filtre est donc très souple et permet d'envisager initialement un très grand nombre de solutions différentes, dont les plus absurdes seront ensuite 15 éliminées par le filtre particulaire. Au cours d'une sous-étape 113 suivante, le calculateur associe chaque particule Pi à son plus proche tronçon de route. La méthode choisie est ici du type « point-to-curve ». Elle consiste à associer chaque particule Pi au tronçon de route qui est le plus proche au sens de 20 la distance Euclidienne. A titre d'exemple illustratif, sur la figure 4, on observe ainsi que la particule P1 est associée au tronçon de route AB. A ce stade, le calculateur peut orienter les particules Pi en fonction notamment de l'orientation du tronçon de route auquel chaque particule est 25 associée (et possiblement aussi en fonction de la dynamique du véhicule). Le procédé se poursuit ensuite en une sous-étape 116 qui sera décrite plus loin. Comme cela a été décrit supra, la première sous-étape 110 consistait à déterminer si l'on se trouvait ou non dans une phase d'initialisation du filtre 30 particulaire. On peut maintenant envisager que tel n'est pas le cas et que le processus a déjà été initialisé auparavant. Dans ce cas, au cours d'une sous-étape 114, le calculateur met à jour les particules Pi sur la carte. 9 Pour cela, les particules Pi sont toutes déplacées sur la carte en fonction d'informations relatives à la dynamique du véhicule. Les deux données que sont la vitesse V du véhicule et la vitesse angulaire en lacet Y du véhicule automobile 10 sont en effet employées pour 5 déplacer toutes les particules Pi d'une distance donnée et pour réorienter les particules d'un angle donné. On notera que cette sous-étape n'utilise cette fois pas la position géolocalisée Po du véhicule automobile. Au cours d'une sous-étape 115 suivante, le calculateur ré-associe 10 chaque particule Pi à un tronçon de route. Il détermine plus précisément quelles sont les particules Pi qui doivent être associées à un nouveau tronçon de route, et il identifie ce nouveau tronçon de route. Pour comprendre comment le calculateur opère, on peut se référer à la 15 figure 4 sur laquelle deux particules P1, P2, centrées sur des points Ml, M2, sont représentées et sur laquelle est également représenté un tronçon de route AB. On considère ici qu'au pas de temps précédent, les deux particules P1, P2, étaient associées à un seul et même tronçon de route AB, puis qu'elles ont été déplacées au cours de la sous-étape 114. 20 Le calculateur détermine alors pour chaque particule Pi un ratio r, afin de savoir si chaque particule doit ou non être associée à un nouveau tronçon de route. Ce ratio r est calculé selon la formule suivante : (AB - AM) r= Dans le cas où ce ratio r est compris entre 0 et 1, l'association de la 25 particule Pi avec son tronçon de route d'origine ne doit pas être changée. C'est le cas ici pour la particule P1. Dans le cas où ce rapport est négatif, l'association de la particule Pi avec son tronçon de route doit être changée. Cette particule doit plus précisément être associée avec le tronçon de route précédent ou avec l'un des tronçons de route 30 précédents. Dans le cas où ce rapport est strictement supérieur à 1, l'association de 10 la particule Pi avec son tronçon de route doit être changée. Cette particule doit plus précisément être associée avec le tronçon de route suivant ou avec l'un des tronçons de route suivants. Plusieurs situations peuvent ainsi être rencontrées. 5 Dans la situation de la figure 4 où le tronçon de route AB ne comprend qu'un unique successeur BB', la particule P2 est associée à ce successeur (pour autant que le ratio r soit compris entre 0 et 1 avec ce nouveau tronçon de route, sinon un autre successeur est envisagé). Dans la situation de la figure 5 où le tronçon de route AB comprend de 10 multiples successeurs BC, BD, BE, la particule P2 considérée au pas de temps précédent est clonée en autant de particules P21, P22, P23 qu'il y a de successeurs BC, BD, BE. On peut prévoir aussi de cloner la particule moins de fois si certains des successeurs ne sont pas envisageables, compte tenu de la dynamique du 15 véhicule. Dans une autre situation non représentée sur les figures, il est possible que la particule doive être associée à un autre tronçon de route parallèle au tronçon de route auquel il était associé au pas de temps précédent (ce qui arrivera notamment lorsque le véhicule change latéralement de voie de circulation, par 20 exemple pour doubler un autre véhicule). Cela est rendu possible car les particules ne sont pas contraintes à évoluer uniquement sur le même tronçon de route. Cette situation peut être détectée compte tenu de la nouvelle position de la particule et des données stockées dans la carte (informations de marquage au sol, largeurs des voies de circulation...). Dans une variante, on peut envisager que 25 cette situation soit détectée à l'aide en outre de caméras embarquées dans le véhicule. Au cours d'une sous-étape 116 qui suit aussi bien la sous-étape 115 que la sous-étape 113, le calculateur calcule la vraisemblance wi de chaque particule Pi. 30 La vraisemblance d'une particule est ici exprimée par son poids wi. Plus le poids d'une particule est grand, plus il est probable que la particule considérée corresponde à la position exacte du véhicule automobile 10. Ce poids peut être calculé de différentes façons. Dans un premier mode de réalisation, le poids wi de chaque particule Pi 11 est calculé en fonction uniquement de données issues de la carte. Il est plus précisément déterminé en fonction de la distance Euclidienne qui sépare la particule considérée du tronçon de route auquel elle est associée (ce poids est par exemple inversement égal à cette distance). 5 Dans un second mode de réalisation, le poids wi de chaque particule Pi est calculé en fonction également de données issues des capteurs extéroceptifs, à condition que ces données soient jugées fiables. On peut en effet imaginer d'augmenter ou de réduire le poids de la particule considérée en fonction des informations latérales en provenance des 10 caméras CAM du véhicule. Ces caméras sont effectivement capables de détecter les lignes de marquage au sol et de les retourner au calculateur sous la forme d'un modèle polynomial. Le calculateur peut alors vérifier si la forme de ces lignes correspond à celle des marquages au sol enregistrés dans la carte, et ajuster le poids de la particule en correspondance. 15 On peut noter que les marquages au sol ne sont pas toujours détectés par les caméras. Ceci peut être dû à des conditions difficiles pour les capteurs comme par exemple un mauvais éclairage, une route mouillée, des marquages effacés, etc... Dans ces cas particuliers, la caméra indique au calculateur un niveau de confiance faible et le calcul du poids se base alors uniquement sur les 20 données fournies par la carte comme décrit dans le premier mode de réalisation. On observera ici que le critère de fiabilité d'une donnée est généralement fourni par le capteur qui mesure la donnée, par exemple sous la forme d'n pourcentage (lequel pourcentage sera alors utilisé pour déterminer si la donnée doit ou non être prise en compte). 25 Quelle que soit la méthode utilisée, le procédé se poursuit en une sous- étape 117 de sélection d'un ensemble restreint de particules Ph de façon à éliminer celles qui se seraient trop écartée de la position géolocalisée Po instantanée du véhicule automobile 10. Pour mettre en oeuvre cette sous-étape, le calculateur acquiert la 30 nouvelle position géolocalisée Po du véhicule automobile 10, puis il calcule la distance séparant chaque particule Pi de cette position géolocalisée Po instantanée. Si cette distance est supérieure au niveau de protection horizontal HPL, le poids wi de la particule Pi correspondante est fixé à zéro, ce qui permettra à 12 cette particule d'être automatiquement éliminée par la suite. Dans le cas contraire, le poids wi de la particule Pi correspondante n'est pas modifié. Au cours d'une sous-étape 118 suivante, le calculateur détermine s'il est 5 nécessaire ou non de ré-échantillonner des particules Pi sur la carte. Il utilise pour cela un indicateur Nen, qui est calculé en fonction du poids wi des particules Pi et du nombre de particules P. Si cet indicateur Neff passe en dessous d'un seuil prédéterminé (stocké dans la mémoire morte du calculateur), alors le calculateur ré-échantillonne les 10 particules Pi sur la carte. Sinon, les particules Pi sont conservées dans leur état. Comme on le sait, un ré-échantillonnage consiste à considérer les particules (ci-après appelées particules originales) dans leur ensemble, et à tirer de nouvelles particules dans cet ensemble original. Pour ré-échantillonner les particules, le calculateur pourrait utiliser une 15 méthode classique au cours de laquelle il tirerait aléatoirement un nombre prédéfini de nouvelles particules dans l'ensemble original de particules Ph la probabilité de tirer chaque particule Pi étant proportionnelle au poids wi de cette particule P. Toutefois, cette méthode cause généralement un appauvrissement des particules, puisque ce sont toujours celles qui ont un poids très important qui 20 sont tirées. De façon préférentielle, le calculateur utilise plutôt ici une méthode de ré-échantillonnage dit « à faible variance » (de l'anglais « Low Variance Resampling »). Cette méthode favorise en effet le maintien d'une bonne répartition des particules sur la carte. Cette méthode consiste à tirer aléatoirement un 25 nombre prédéfini de nouvelles particules dans l'ensemble original de particules Ph la probabilité de tirer chaque particule Pi étant fonction du poids wi de cette particule Pi mais n'étant cette fois pas proportionnelle à ce poids. A ce stade, le calculateur pourrait se contenter de recommencer en boucle les sous-étapes 114 à 118 jusqu'à obtenir des particules situées toutes 30 autour d'un seul et même point, lequel serait considéré comme correspondant à la position précise Pp du véhicule automobile 10 sur la carte. Ce n'est toutefois pas ici l'option qui est choisie. Ainsi, comme cela a été exposé précédemment, une fois la sous-étape 118 achevée, il est prévu une opération de sélection d'hypothèses 200. 13 Cette opération de sélection d'hypothèses 200 est mise en oeuvre une fois que l'opération de filtrage de particules 100 a convergé et a donné un nombre restreint de solutions (les particules étant par exemple groupées autour d'un nombre de points inférieur à un seuil prédéterminé). 5 Cette opération de sélection d'hypothèses 200 est mise en oeuvre de manière récursive, c'est-à-dire en boucle et à pas de temps réguliers. Elle comporte plusieurs étapes successives. Au cours d'une première étape 201, le calculateur sélectionne des « hypothèses ». 10 Il considère pour cela les particules Pi en différents ensembles au sein de chacun desquels les particules sont toutes associées à une même voie de circulation (ou, en variante, à un même tronçon de route). L'intérêt de travailler sur des hypothèses est qu'il va alors être possible de sélectionner toutes les hypothèses les plus vraisemblables, ce qui permettra, 15 d'une part, de garder la bonne hypothèse parmi celles sélectionnées, et, d'autre part, de vérifier la validité de chaque hypothèse sélectionnée. Les hypothèses peuvent être formulées sous la forme d'assertions telles que « le véhicule se situe dans la voie de circulation dont la référence est ... ». Pour bien comprendre à quoi correspond une hypothèse au sens du 20 présent exposé, on a regroupé sur la figure 3 les particules en huit ensembles Z1, Z2, Z3, Z4, Z5, Z6, Z7, 43 correspondant chacun à une hypothèse. A titre d'exemple, les particules de l'ensemble Z1 correspondent à l'hypothèse « le véhicule se situe dans la voie de circulation de droite de la route R1 ». 25 Les particules de l'ensemble Z2 correspondent à l'hypothèse « le véhicule se situe dans la voie de circulation de gauche de la route R1 ». Les particules de l'ensemble Z3 correspondent à l'hypothèse « le véhicule se situe dans la voie de circulation de gauche de la route R2 >>. Les particules de l'ensemble Z4 correspondent à l'hypothèse « le 30 véhicule se situe sur le rond-point, entre ses jonctions avec les routes R1 et R2 >>... En considérant qu'il est trouvé un nombre « J » d'hypothèses (sur la figure 3, J=8), chaque hypothèse peut également être exprimé sous la forme d'un vecteur Xj dont les composantes correspondent à la somme des coordonnées des particules Pi de cette hypothèse, pondérées par le poids wi de ces particules. 14 Le calculateur peut affecter à chaque hypothèse un « indice de confiance » égal à la somme des poids wi des particules Pi de cette hypothèse. Au cours d'une seconde étape 202, le calculateur va déterminer la matrice de covariance Y..»-(X) de chaque hypothèse et la matrice de covariance 5 E (XGNss) de la position géolocalisée Po du véhicule 10. Manipuler de telles matrices de covariance permet en effet de caractériser l'incertitude liée à chaque hypothèse et celle liée à la position géolocalisée Po fournie par le système de géolocalisation. Comme cela a été exposé supra, la matrice de covariance E-,---GNSS (X 1 1 , 10 liée à la position géolocalisée Po du véhicule 10 est ici directement transmise au calculateur par le système de géolocalisation. Il s'agit ici d'une matrice 2x2. S'agissant de la matrice de covariance E-1() liée à chaque hypothèse, elle est calculée en fonction des poids wi de l'ensemble des particules Pi associées à cette hypothèse. Il s'agit ici aussi d'une matrice 2x2. 15 Il est ensuite nécessaire de déterminer dans quelle mesure chaque hypothèse est « cohérente », compte tenu de la position géolocalisée Po fournie par le système de géolocalisation et en tenant compte de l'erreur liée à la mesure de cette position géolocalisée. Pour cela, au cours d'une étape 203, on utilise un objet mathématique 20 appelé distance de Mahalanobis DM], dont l'expression est la suivante : Dmj() = - xGNss)T. [I-1() - 1-1(xGNss)] OÙ XGNss correspond au vecteur « position géolocalisée Po ». La distance de Mahalanobis est en effet un objet qui permet d'évaluer la cohérence entre deux situations incertaines, en prenant en compte les covariances des variables (c'est-à-dire le doute lié à chaque variable). 25 Alors, au cours d'une étape 204, il est prévu de sélectionner un premier ensemble restreint (voire vide) d'hypothèses parmi les hypothèses acquises à l'étape 201. On réalise pour cela un test de Khi-deux (X2) pour chaque distance de Mahalanobis DM]. 30 En pratique, chaque distance de Mahalanobis DM] est ici comparée à un seuil critique pour déterminer si l'hypothèse considérée est ou non cohérente avec 15 la position géolocalisée Po. Si l'hypothèse considérée et la position géolocalisée Po sont cohérentes au sens du test de KHI-deux, l'hypothèse est conservée. Au contraire, si l'hypothèse considérée et la position géolocalisée Po ne 5 sont pas cohérentes au sens du test du KHI-deux, l'hypothèse est rejetée. On notera ici qui si une hypothèse est conservée, cela ne signifie pas nécessairement que cette hypothèse est vraie. En effet, à ce stade, plusieurs hypothèses peuvent être conservées. A contrario, si une hypothèse est rejetée, cela ne signifie pas 10 nécessairement que cette hypothèse était fausse. Il peut en effet arriver qu'une grosse erreur affecte la mesure de la position géolocalisée Po. Dans ce cas, une hypothèse vraie peut être rejetée. Comme cela apparaîtra clairement dans la suite de cet exposé, cela n'affectera pas pour autant la fiabilité de la méthode ici proposée. 15 Au cours d'une étape 205 subséquente, il est prévu de sélectionner un second ensemble restreint (voire vide) d'hypothèses parmi les hypothèses sélectionnées à l'étape 204. On notera ici que cette seconde sélection aurait pu être opérée avant la première sélection sans pour autant affecter le déroulement du procédé. 20 Cette seconde sélection consiste à ne conserver que les hypothèses « vraisemblables », pour lesquelles un indicateur, lié aux poids \At; des particules Pi qui composent cette hypothèse, est supérieur à un seuil déterminé. L'objectif est en effet d'éliminer les hypothèses qui ont satisfait le test de cohérence de KHI-deux, mais qui sont peu probables. 25 Pour cela, le calculateur élimine les hypothèses pour lesquelles l'indice de confiance (dont on rappelle qu'il est égal à la somme des poids wi des particules Pi de l'hypothèse considérée) est inférieur à un seuil déterminé. Ce seuil est ici invariable et mémorisé dans la mémoire morte du calculateur. A l'issue de ces deux étapes de sélection d'hypothèses, le calculateur a 30 retenu un nombre N d'hypothèses non seulement cohérentes mais aussi vraisemblables. Au cours d'une étape 206, il est alors prévu de déterminer le caractère utilisable ou non de chaque hypothèse sélectionnée, en fonction de ce nombre N. Trois cas sont alors envisageables. 16 Le premier cas est celui où le nombre N est égal à 1. Dans ce cas, puisqu'une seule hypothèse a été retenue, cette hypothèse est considérée comme juste et utilisable pour générer une consigne de conduite du véhicule autonome. Il s'ensuit que le calculateur peut s'y fier. Dans ce cas, le calculateur peut alors 5 considérer que la particule de plus gros poids de cette hypothèse correspond à la position précise Pp du véhicule 10. Le second cas est celui où le nombre N est strictement supérieur à 1. Dans ce cas, puisque plusieurs hypothèses ont été retenues, aucune n'est considérée comme utilisable pour générer une consigne de conduite du véhicule 10 autonome. Le procédé est donc redémarré en boucle jusqu'à obtenir une seule et unique hypothèse. Le dernier cas est celui où le nombre N est égal à 0. Dans ce cas, puisqu'aucune hypothèse n'a été retenue, aucune particule n'est considérée comme utilisable pour générer une consigne de conduite du véhicule autonome. 15 Par ailleurs, le calculateur peut avantageusement déduire de cette situation qu'il existe une incohérence entre les mesures effectuées par le système de géolocalisation et les hypothèses acquises, ce qui est probablement dû à un problème affectant le système de géolocalisation. Dans cette éventualité, il est prévu une étape 207 d'émission d'une alerte à destination du conducteur et/ou de 20 l'unité de pilotage du véhicule en mode autonome, de façon à ce que celui-ci puisse prendre les mesures qui s'imposent (arrêt d'urgence, pilotage en mode dégradé, REVENDICATIONS1. Procédé d'estimation d'une position précise (Pp) d'un véhicule (10) sur une carte, comportant : - une étape d'acquisition d'au moins une position géolocalisée (Po) du véhicule (10) au moyen d'un système de géolocalisation, - une étape de pré-positionnement du véhicule (10) sur la carte, à sa position géolocalisée (Po), et - une étape de filtrage particulaire au cours de laquelle des positions possibles du véhicule, appelées particules (Pi), sont traitées pour déterminer la position précise (Pp) du véhicule (10) sur la carte, caractérisé en ce que l'étape de filtrage particulaire comprend : - une étape préalable de répartition de particules (Pi) sur la carte, en fonction de la position géolocalisée (Po) du véhicule (10) sur la carte, puis une étape de mise à jour des particules (Pi) sur la carte, - une étape de calcul de la vraisemblance (wi) de chaque particule (Pi), à partir au moins de données issues de la carte, - une étape de sélection d'un ensemble restreint de particules (Pi), et - si un indicateur (Neff), relatif à la vraisemblance (wi) des particules (Pi) sélectionnées et au nombre de particules (Pi) sélectionnées, passe en dessous d'un seuil, une étape de ré-échantillonnage de particules (Pi) sur la carte. 2. Procédé d'estimation selon la revendication 1, dans lequel la vraisemblance (wi) de chaque particule (Pi) est calculée en fonction uniquement de données issues de la carte. 3. Procédé d'estimation selon la revendication 1, dans lequel la vraisemblance (wi) de chaque particule (Pi) est calculée en fonction également de données issues de capteurs permettant au véhicule (10) de percevoir son environnement, à condition que ces données soient jugées fiables. 4. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape de filtrage particulaire, la position précise (Pp) est choisie parmi l'ensemble restreint de particules (Pi) sélectionnées, en fonction de la vraisemblable (wi) de chaque particule (Pi). 5. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, la carte stockant des données relatives à des tronçons de routes, la 18 vraisemblance (wi) de chaque particule (Pi) est calculée en fonction de la position du plus proche tronçon de route par rapport à la particule (Pi). 6. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape préalable, les particules (Pi) sont réparties dans un disque 5 centré sur la position géolocalisée (Po) du véhicule (10). 7. Procédé d'estimation selon la revendication 6, dans lequel le rayon du disque est déterminé en fonction du niveau de protection horizontal (HPL) affecté à la position géolocalisée (Po) du véhicule (10). 8. Procédé d'estimation selon l'une des revendications précédentes, 10 dans lequel, à l'étape de sélection, les particules (Pi) sont sélectionnées en fonction de la distance les séparant de la position géolocalisée (Po) du véhicule (10). 9. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape de mise à jour, les particules (Pi) sont déplacées sur la carte 15 en fonction uniquement d'informations relatives à la dynamique du véhicule. 10. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape de ré-échantillonnage, les particules (Pi) sont ré-échantillonnées à l'aide d'une technique à faible variance. 11. Véhicule (10) comportant : 20 - des moyens de mémorisation d'une carte, - un système de géolocalisation, et - un calculateur adapté à pré-positionner le véhicule (10) sur la carte, caractérisé en ce que le calculateur est adapté à mettre en oeuvre un procédé d'estimation d'une position précise (Pp) du véhicule (10) sur la carte, 25 conforme à l'une des revendications précédentes.