Brevet FR3143814A1 — publié le 20240621

L’invention se rapporte à un procédé de caractérisation de l’environnement d’un porteur mobile, dans lequel, pour chaque itération à un instant , les étapes suivantes sont mises en œuvre : S10) Acquisition d'une pluralité de mesures de distance dans l’environnement par au moins un capteur ; S20) Génération d’une paire de grilles d’occupation à l’instant et à l’instant , chaque grille ( fusionnant les mesures de distance dans une représentation spatiale discrétisée de l’environnement ; S30) Génération d’une grille d’espace statique à l’instant , ou S40) Génération d’une grille d’espace libre à l’instant , Figure pour l’abrégé : Fig. 7 Procédé de caractérisation de l’environnement d’un porteur mobile, produisant une grille d’espace statique et/ou une grille d’espace libre L'invention se rapporte à un procédé permettant de caractériser l’environnement d’un porteur mobile (également appelé robot par la suite) à l'aide de capteurs de distance. Elle concerne également un dispositif pour la mise en œuvre d'un tel procédé. L'invention peut être appliquée à tout type de robots, par exemple des robots terrestres (voitures autonomes), robots aériens (drones), robots sous-marins, etc. Un robot mobile navigue dans l'espace pour effectuer une tâche. Pour naviguer de manière autonome, le robot doit acquérir des connaissances sur son environnement afin de planifier un chemin, éviter les collisions et naviguer en toute sécurité. Le robot s'appuie à cet effet sur un système de perception. Le système de perception récupère des données sur l'environnement grâce à des mesures à partir de capteurs embarqués. Un modèle informatique de l'environnement est ensuite construit et mis à jour en continu par le système de perception basé sur des mesures de capteurs. Le robot raisonnera et prendra des décisions de navigation basées sur le modèle d'environnement. Plusieurs types de modèles d'environnement existent dans l'état de l'art. L'un des plus couramment utilisés en robotique est la grille d’occupation initialement introduite dans l’article de H.P. Moravec et A. Elfes « High resolution maps from wide angle sonar », Proceedings. IEEE ICRA, volume 2, pages 116-121, Mar 1985. La grille d’occupation est une grille de probabilités : chaque cellule a une probabilité d’occupation, donc une valeur comprise entre 0 et 1. Une grille d’occupation est une carte d'une région délimitée de l'environnement entourant le robot. La région est subdivisée en plusieurs sous-régions adjacentes appelées cellules. Une cellule peut avoir trois dimensions (3D), par exemple, un cube de 10 cm × 10 cm × 10 cm. Il peut également être en deux dimensions (2D), par exemple, l'espace délimité par un carré de 10 cm × 10 cm. La collection de toutes les cellules forme une grille. L'objectif du cadre défini par la grille d’occupation est de calculer la probabilité qu'une cellule soit occupée par un obstacle compte tenu des mesures des capteurs embarqués. Une probabilité d'occupation de 1 signifie « la cellule est occupée », 0 signifie « la cellule est libre », et une valeur comprise entre 0 et 1 exprime l'incertitude sur l'occupation de la cellule. En particulier, une probabilité d'occupation de 1/2 signifie que l'occupation de la cellule est inconnue, c'est-à-dire que la cellule peut être libre ou occupée, sans pencher pour une des deux hypothèses. Les grilles d’occupation présentent des propriétés avantageuses pour un système de perception. Ils constituent une méthode robuste pour fusionner plusieurs capteurs de différentes technologies. La demande de brevet FR3041451 fournit un exemple de procédé de fusion multi-capteurs mettant en œuvre des calculs en arithmétique entière. Les technologies de détection couramment utilisées sont les LIDAR, les RADAR, les caméras et les SONAR. Alors que les caméras donnent des couleurs et des textures détaillées, elles sont limitées dans l'obscurité. Les LIDAR peuvent donner une représentation détaillée de nuages de points, mais leurs performances diminuent en cas de mauvais temps (neige, brouillard ou pluie). La fusion multi-capteurs est pour cela essentielle pour surmonter la limitation d'une technologie donnée. Les grilles d'occupation permettent également de gérer les incertitudes de mesure et les bruits par des approches probabilistes. Les grilles de caractéristiques représentent ainsi l'espace libre et l'espace occupé par un ou plusieurs corps matériels. On appelle « corps matériels » (ou obstacle) toute substance ou objet matériel présentant une individualité et pouvant être détecté et identifié par un capteur approprié. Ainsi, sont considérés des corps matériels les objets inanimés, qu'ils soient naturels ou artificiels, les végétaux, les animaux, les êtres humains, mais aussi des particules liquides ou solides en suspension dans l'air, comme les nuages, voire des masses liquides ou gazeuses. Par exemple, pour une voiture autonome, les corps matériels courants sont les piétons, les autres voitures, les bâtiments, la végétation, les cyclistes, etc. Les corps matériels peuvent être mobiles, c'est-à-dire dynamiques ou statiques. L’état d’occupation d’une cellule peut varier dans le temps. En effet, une cellule peut être libre pendant un certain temps, puis subitement occupée par un obstacle mobile, puis à nouveau libre lorsque l'obstacle est parti. Par exemple, un piéton peut se lever sans bouger pendant une courte période, puis marcher dans une direction. Le cadre standard des grilles de caractéristiques calcule les probabilités d'occupation de toutes les cellules d'une grille à chaque itération t, en ne prenant en compte que les mesures instantanées des capteurs à la même itération. Soit la grille de caractéristiques construite à l'itération (ou instant ). Considérons une cellule qui a toujours été estimée libre par toutes les grilles de caractéristiques de à . A l'itération , supposons que la cellule est temporairement occultée des capteurs par un obstacle mobile. L'occlusion fait que l'état d'occupation de la cellule passe soudainement de libre à inconnu à l'itération . Cela signifie que les estimations précédentes ont été brusquement oubliées même s'il est possible que la cellule reste libre à l'itération alors qu'elle est temporairement occluse. Le même effet apparaît pour les cellules statiques, c'est-à-dire les cellules occupées par des obstacles fixes. Ces cellules sont estimées occupées tant qu'elles sont présentes dans le champ de vision des capteurs. Mais, leur état d'occupation devient subitement inconnu dès qu'ils sont obstrués par un obstacle, même en mouvement. Cependant, comme elles sont occupées par un obstacle statique, il est possible que ces cellules restent occupées pendant l'occlusion. Plusieurs approches existent dans l'état de l'art pour tenir compte de l’historique des mesures afin d’estimer les cellules libres et les cellules occupées par un obstacle statique. Les approches les plus prometteuses consistent à utiliser des filtres à particules afin de séparer les cellules occupées par des obstacles dynamiques des cellules occupées par des obstacles statiques. Les filtres à particules sont décrit par exemple dans l’article de Nuss, D., Reuter, S., Thom, M., Yuan, T., Krehl, G., Maile, M., Dietmayer, K. (2018) « A random finite set approach for dynamic occupancy grid maps with real-time application », The International Journal of Robotics Research, 37(8), 841–866. La définition d'une particule dans ces algorithmes peut différer mais en général, les particules sont un ensemble de points répartis dans l'espace. Les particules ont des coordonnées de position et un vecteur de vitesse. Une particule a également un poids utilisé pour estimer les distributions de vitesse. Entre les itérations successives, une particule se déplace en suivant sa vitesse (modulo l'introduction d'une erreur aléatoire). L'idée de base de ces algorithmes est de répandre au hasard des particules dans l'espace de sorte que leur mouvement suive le mouvement des corps matériels dynamiques. Les particules sont distribuées au sein de la grille de sorte qu'un grand nombre de particules sont présentes dans les cellules occupées par des corps matériels dynamiques. Ces approches estiment à la fois où sont les cellules libres, où sont les statiques et où sont les dynamiques, aboutissant à une solution avec une forte demande de calcul. Par conséquent, ces approches nécessitent un matériel informatique à hautes performances de calcul tels que les GPU (« Graphics Processing Unit » ou processeur graphique) pour atteindre une exécution en temps réel, ce qui constitue un frein à leur intégration dans du matériel embarqué ayant des performances de calcul plus limitées et une consommation d'énergie plus faible, tels que les microcontrôleurs. L’invention vise donc à fournir un procédé qui tient compte de l’historique des mesures afin d’estimer les cellules libres et les cellules occupées par un obstacle statique, et qui puisse être intégré dans du matériel embarqué du type microcontrôleur. Un objet de l’invention est donc un procédé de caractérisation de l’environnement d’un porteur mobile, dans lequel, pour chaque itération à un instant , les étapes suivantes sont mises en œuvre : S10) Acquisition d'une pluralité de mesures de distance dans l’environnement par au moins un capteur ; S20) Génération d’une paire de grilles d’occupation à l’instant et à l’instant , chaque grille fusionnant les mesures de distance dans une représentation spatiale discrétisée de l’environnement ; S30) Génération d’une grille d’espace statique à l’instant, chaque cellule de la grille d’espace statique à l’instant ayant une probabilité postérieure que ladite cellule soit occupée par un corps matériel statique à l’instant , ladite probabilité postérieure étant calculée pour chaque cellule par un filtre bayésien binaire à partir : de la probabilité postérieure que ladite cellule soit occupée par un corps matériel statique à l’instant , ladite probabilité postérieure que ladite cellule soit occupée par un obstacle statique à l’instant étant directement injectée en entrée du filtre bayésien binaire ; d’un modèle de prédiction dit statique comprenant une probabilité que la cellule soit occupée par un corps matériel statique à l’instant en fonction de l’état de la cellule à l’instant ; et d’un modèle inverse dit statique à l’instant qui comprend une probabilité de connaissance de l’occupation de la cellule par un corps matériel statique à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant ; ou S40) Génération d’une grille d’espace libre à l’instant , chaque cellule de la grille d’espace libre à l’instant t ayant une probabilité postérieure que ladite cellule soit libre à l’instant , ladite probabilité postérieure étant calculée pour chaque cellule par un filtre bayésien binaire à partir : de la probabilité postérieure que ladite cellule soit libre à l’instant , ladite probabilité postérieure que ladite cellule soit libre à l’instant étant directement injectée en entrée du filtre bayésien binaire ; d’un modèle de prédiction dit libre comprenant une probabilité que la cellule soit libre à l’instant en fonction de l’état de la cellule à l’instant ; et d’un modèle inverse dit libre à l’instant qui comprend une probabilité de connaissance du fait que la cellule soit libre à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant . Avantageusement, l’étape S30) de génération d’une grille d’espace statique et l’étape S40) de génération d’une grille d’espace libre sont mises en œuvre de façon concomitante. Avantageusement, le procédé comprend une étape S50) de génération d’une grille combinée résultant de la combinaison de la grille d’espace statique et de la grille d’espace libre. Avantageusement, chaque cellule de la grille combinée est calculée par la fusion bayésienne suivante : Dans lequel est la fonction de fusion bayésienne, correspond à une cellule de la grille combinée à l’instant , correspond à une cellule de la grille d’espace statique à l’instant , et correspond à une cellule de la grille d’espace libre à l’instant , Et selon la valeur de : si , la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. Avantageusement : - la probabilité postérieure que la cellule soit occupée par un corps matériel statique à l’instant et le modèle inverse statique à l’instant sont approximés par des valeurs appartenant à un ensemble de cardinalité finie, les valeurs étant identifiées respectivement par un indice de probabilité et par un indice de modèle inverse statique ; - la probabilité postérieure que la cellule soit libre à l’instant et le modèle inverse libre à l’instant sont approximés par des valeurs appartenant à un ensemble de cardinalité finie, les valeurs étant identifiées respectivement par un indice de probabilité et par un indice de modèle inverse libre . Avantageusement, l’indice de modèle inverse statique est calculé de la façon suivante : La fonction renvoyant une valeur positive, correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant , et correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant . Avantageusement, l’indice de modèle inverse libre est calculé de la façon suivante : et correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant , la fonction renvoyant une valeur positive. Avantageusement la valeur positive de la fonction et la valeur positive de la fonction sont des valeurs constantes. Avantageusement, la valeur positive renvoyée par la fonction est égale à . Avantageusement la valeur positive renvoyée par la fonction est égale à . Avantageusement, l’indice de modèle inverse libre est calculé de la façon suivante : Avantageusement, l’indice correspondant à la composante de prédiction du filtre bayésien binaire de la grille d’espace statique est obtenu à l’aide d’une table de consultation comprenant un ensemble fini d’indices de probabilités, et l’indice filtré est obtenu par addition d’indices entiers de probabilités, et dans lequel l’indice correspondant à la composante de prédiction du filtre bayésien binaire de la grille d’espace libre est obtenu à l’aide d’une table de consultation comprenant un ensemble fini d’indices de probabilités, et l’indice filtré est obtenu par addition d’indices entiers de probabilités. Avantageusement, chaque cellule de la grille combinée est calculée par la fusion bayésienne suivante : correspond à une cellule de la grille combinée à l’instant , correspond à une cellule de la grille d’espace statique à l’instant , et correspond à une cellule de la grille d’espace libre à l’instant , Et selon la valeur de : si la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. L’invention se rapporte aussi à un procédé d’évitement d’un corps matériel en mouvement autour d’un corps mobile, caractérisé en ce qu’il met en œuvre le procédé précité pour caractériser l’environnement d’un porteur mobile, et qu’il envoie une commande à un actionneur du porteur mobile afin d’éviter ledit corps matériel. L’invention se rapporte aussi à un dispositif de caractérisation de l’environnement d’un porteur mobile, comprenant: - au moins un port d'entrée pour recevoir une pluralité de signaux représentatifs d'une série temporelle de mesures de distance issues d'un ou plusieurs capteurs de distance, et - un processeur de données configuré pour recevoir en entrée lesdits signaux, pour générer une grille d’espace libre ou une grille d’espace statique à partir desdits signaux en appliquant le procédé précité. Description des figures D’autres caractéristiques, détails et avantages de l’invention ressortiront à la lecture de la description faite en référence aux dessins annexés donnés à titre d’exemple. La est un schéma-bloc illustrant les différentes étapes du procédé selon l’invention. La est un schéma illustrant le filtre bayésien binaire. La est un schéma illustrant le filtre bayésien binaire entier. La est un schéma illustrant le filtre bayésien binaire entier dans lequel le modèle de prédiction est constant. La est un schéma illustrant le procédé de caractérisation de l’environnement selon l’invention par génération d’une grille d’espace statique. La est un schéma illustrant le procédé de caractérisation de l’environnement selon l’invention par génération d’une grille d’espace libre. La est un schéma illustrant le procédé de caractérisation de l’environnement selon l’invention par génération d’une grille combinée. La est un schéma illustrant la façon dont le procédé selon l’invention se distingue d’un algorithme d’estimation du déplacement d’un corps matériel connu de l’homme du métier. Les différentes étapes du procédé selon l’invention sont illustrées sur la . Acquisition des mesures de distance Le procédé fonctionne de façon itérative : à chaque itération correspondant à un instant , on effectue une nouvelle série de mesures (étape S10). Le fait d’effectuer les mêmes calculs avec des mesures issues de plusieurs capteurs ne pose aucune difficulté pour l’homme du métier. L’ensemble de mesures issue d'un ou plusieurs capteurs à l’instant permet de déterminer la probabilité d'occupation d'une cellule parmi les cellules d’une grille d’occupation telle que décrite précédemment. Ainsi, chaque cellule est identifiée par un indice allant de 1 à . La référence désigne en fait un ensemble de mesures lorsque le système est composé de plusieurs capteurs ou lorsqu'un capteur produit plusieurs mesures à chaque itération. Par exemple, un balayage LIDAR produit un nuage de points composé d'un nombre fini de points. Génération de la paire de grilles d’occupation La deuxième étape S20 consiste à générer une paire de grilles d’occupation. A chaque instant , une grille d'occupation fusionne les mesures produites par tous les capteurs du système à la même itération. Il est rappelé que pour une grille d'occupation , une cellule a un état binaire : elle est occupée par un obstacle ou non. Dans le cadre des grilles d’occupation, il est communément supposé que « non occupé » signifie « libre ». A chaque instant , on calcule ainsi la probabilité qu'une cellule soit occupée en prenant en compte toutes les mesures de capteurs produites au même instant . représente la probabilité d'occupation d'une cellule compte tenu des mesures . Ainsi, une nouvelle grille d'occupation est générée l’instant , et la grille d'occupation , générée l’instant , est stockée en mémoire afin d’être réutilisée pour la paire de grilles d’occupation. désigne la séquence ( dans la présente demande. Génération de la grille d’espace statique Dans une troisième étape S30, une grille d’espace statique est générée à l’instant . La grille d’espace statique est un réseau de éléments, dans lequel : Chaque cellule de la grille d’espace statique a une probabilité postérieure que, à l’instant , la cellule soit occupée par un corps matériel statique compte tenu de la séquence de paires de grilles d’occupation . Ladite probabilité postérieure est calculée pour chaque cellule de la grille d’espace statique par un filtre bayésien binaire, dont le principe est rappelé en . Le filtrage bayésien binaire consiste à estimer de façon récursive un état caché (ce qui est le cas si certaines cellules sont obstruées par un obstacle, de sorte que la mesure de distance soit inconnue) au vu d'observations. est appelée probabilité postérieure, car elle dépend de la valeur spécifique des observations . Comme l’illustre la figure 2, la probabilité postérieure est réinjectée en entrée du filtre pour l’itération suivante. Le filtre bayésien binaire comprend essentiellement deux étapes, à savoir une étape de prédiction, qui dépend de l’état des paires aux instants antérieurs (probabilité postérieure à l’instant , et modèle de prédiction à l’instant ) et une étape de mise à jour de la mesure, qui prend en compte le modèle inverse à l’instant et la sortie de la prédiction. On peut considérer, de façon avantageuse, que le modèle de prédiction statique (probabilité que la cellule soit occupée par un corps matériel statique à l’instant en fonction de l’état de la cellule à l’instant ) est égal à une valeur constante . Avec cette hypothèse, l’équation de prédiction peut s’écrire : Ainsi, la valeur de influe la prédiction de la façon suivante : Si , le résultat de la prédiction est , ce qui signifie que garde la même valeur que l’estimation précédente de . Si , le résultat de la prédiction est , ce qui signifie que la prédiction de à partir de résulte d’une estimation complètement incertaine de . Une valeur de comprise entre ½ et 1 exprime le degré d’incertitude entre les cas de figures précédents. Dans le cadre de la présente invention, le modèle inverse à l’instant est un modèle inverse dit statique qui comprend une probabilité de connaissance de l’occupation de la cellule par un corps matériel à l’instant et à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant . L’homme du métier peut se référer à la définition détaillée du modèle inverse dans la demande de brevet FR3041451. A chaque instant , le modèle inverse statique est calculé à partir de la paire de grilles d’occupation. Le fait d’utiliser un filtre bayésien binaire avec un modèle inverse comprenant les dernières deux grilles d’occupation (à l’instant et à l’instant ) permet ainsi de déterminer, pour chaque cellule, si celle-ci est occupée par un corps matériel statique. Cette façon d’opérer diffère du procédé divulgué dans la demande de brevet FR3116640, qui estime le mouvement d’un corps matériel en utilisant un modèle inverse qui ne dépend que d’une seule série d’observations temporelles à un instant . Dans la demande de brevet FR3116640, le mouvement est estimé en faisant appel à une grille d’incohérence, générée en réponse à une détection, entre deux itérations, d’un changement d’état d’occupation dans les grilles de caractéristiques. Comme l’illustre la figure 2, la probabilité postérieure que la cellule soit occupée par un obstacle statique à l’instant est directement injectée en entrée du filtre bayésien binaire, ce qui simplifie considérablement le procédé divulgué dans la demande de brevet FR3116640. En effet, dans ce document, le procédé itératif d’estimation du mouvement d’un corps matériel met en œuvre une étape de propagation des indices de probabilités avant que ceux-ci ne soient réinjectés en entrée du filtre. Ainsi, le procédé selon l’invention permet de caractériser les cellules occupées par des obstacles statiques, dans l’environnement d’un porteur mobile, avec une complexité calculatoire inférieure par rapport aux procédés connus. Génération de la grille d’espace libre Dans une étape S40, le même principe peut être appliqué pour générer une grille d’espace libre à l’instant . Ainsi, chaque cellule de la grille d’espace libre à l’instant a une probabilité postérieure que la cellule soit libre à l’instant . La probabilité postérieure est calculée pour chaque cellule par un filtre bayésien binaire à partir : - de la probabilité postérieure que ladite cellule soit libre à l’instant , ladite probabilité postérieure que ladite cellule soit libre à l’instant étant directement injectée en entrée du filtre bayésien binaire ; - d’un modèle de prédiction dit libre comprenant une probabilité que la cellule soit libre à l’instant en fonction de l’état de la cellule à l’instant ; et - d’un modèle inverse dit statique à l’instant qui comprend une probabilité de connaissance du fait que la cellule soit libre à l’instant et à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant . Ainsi, le procédé selon l’invention permet de caractériser les cellules libres, dans l’environnement d’un porteur mobile, également avec une complexité calculatoire inférieure par rapport aux procédés connus. Selon un mode de réalisation, le procédé selon l’invention comprend uniquement la génération d’une grille d’espace statique (c’est-à-dire sans la génération d’une grille d’espace libre), ou, inversement uniquement la génération d’une grille d’espace libre (c’est-à-dire sans la génération d’une grille d’espace statique). En effet, les deux grilles sont indépendantes, et le procédé peut itérer en boucle les étapes S10-S20-S30, ou les étapes S10-S20-S40. Génération de la grille combinée En variante, la grille d’espace statique et la grille d’espace libre sont mises en œuvre de façon concomitante. Ainsi, à chaque itération, une grille d’espace statique et une grille d’espace libre sont générées, à partir des mêmes mesures de distance, avec les mêmes paires de grilles d’occupation. Avantageusement, à chaque itération, la grille d’espace statique et la grille d’espace libre sont combinées entre elles (combinaison cellule par cellule), de façon à obtenir une grille combinée (étape S50, ). Chaque cellule de la grille combinée peut être calculée par la fusion bayésienne suivante : représente la fonction de fusion bayésienne. Ainsi, selon le signe de : si , la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. Ainsi, on regarde l’estimation qui « l’emporte » entre l’estimation de la grille d’espace statique et l’estimation de la grille d’espace libre . Le procédé ainsi décrit utilise le modèle du filtre bayésien binaire. Une façon d’améliorer l’implémentation et l’intégration du procédé dans un microprocesseur serait d’appliquer un filtre bayésien binaire entier, comme l’illustre la figure 3. Le cas où le modèle de prédiction statique est égal à une valeur constante est illustré par la . Pour cela, on identifie la probabilité postérieure que la cellule soit occupée par un corps matériel statique à l’instant par un indice de probabilité . La probabilité postérieure que la cellule soit occupée par un corps matériel statique à l’instant est identifiée par un indice de probabilité . De même, le modèle inverse statique à l’instant est identifié par un indice de modèle inverse statique . De plus, Il est considéré que le modèle de prédiction statique a une valeur constante au cours du temps. Ces indices appartiennent à un ensemble de cardinalité finie de classes de probabilité, et les étapes de prédiction et de mise à jour tirent profit de la fonction de fusion entière pour la fonction de filtrage, en utilisant uniquement une arithmétique entière, telle que décrite dans la demande de brevet FR3116640. En effet, l’ensemble de cardinalité finie de classes de probabilité est formé par l’union d’un ou plusieurs sous-ensembles tels que, lors de l’étape de mise à jour de la mesure, une fusion de deux classes de probabilité appartenant à un même sous-ensemble fournisse un résultat appartenant aussi audit sous-ensemble. L’indice correspondant à la composante de prédiction du filtre bayésien binaire entier de la grille d’espace statique est obtenu à l’aide d’une table de consultation comprenant un ensemble fini d’indices de probabilités, et l’indice filtré est obtenu par addition d’indices entiers de probabilités. Algorithme du filtre bayésien binaire entier de la grille d’espace statique Le filtre bayésien binaire entier de la grille d’espace statique est un algorithme comprenant une fonction d’initialisation du filtre, et une fonction d’application itérative du filtre, avec rebouclage direct depuis la sortie du filtre. L'étape de prédiction, dans la fonction d’initialisation du filtre, est calculée par une table de consultation, et la mise à jour de mesure est calculée par une somme d'indices entiers, comme indiqué ci-dessous. La fonction d’initialisation dépend des paramètres d’entrée suivants : , et est une constante qui correspond au paramètre d’approximation de la séquence de nombres réels entre 0 et 1 utilisés dans la fonction de fusion bayésienne entière, dont les paramètres et sont des valeurs de probabilités, et définie ainsi : L’ensemble d’indices est un sous-ensemble de permettant de définir la séquence . Ainsi, au lieu de travailler sur une séquence infinie , on travaille uniquement dans une séquence finie . Par exemple, pour décoder les indices en 8 bits seulement, on prend . est une constante pour la fonction d’initialisation, avec , et qui correspond à la valeur numérique du modèle de prédiction statique . La première étape de la fonction d’initialisation consiste à déterminer chaque entrée de la table de consultation en partant par un vecteur vide de la taille de , puis, pour chaque dans , l’élément est un indice de probabilité calculé de la manière suivante : où . La fonction approxime le résultat de l’opération par un élément de la séquence . C’est l’indice m de qui est stocké dans l’élément de la table de consultation . La fonction d’application du filtre bayésien binaire entier est définie de la façon suivante : Prédiction: Mise à jour de la mesure: Le procédé selon l’invention calcule, avantageusement l’indice de modèle inverse statique de la façon suivante, illustré par la (bloc « OG vers SIM ) : Si la cellule est probablement occupée à la fois à l’instant et à l'instant , la valeur de l’indice du modèle inverse statique est définie par la fonction qui retourne un indice positif. Sinon, la valeur de l’indice du modèle inverse statique est égale à zéro. correspond à l’indice de probabilité dans la grille d’occupation de la cellule : Selon un mode de réalisation, la fonction renvoie une valeur constante . En variante, on peut définir . Algorithme du filtre bayésien binaire entier de la grille d’espace libre Le filtre bayésien binaire entier de la grille d’espace libre repose sur le même principe que le filtre bayésien binaire entier de la grille d’espace statique . La probabilité postérieure que la cellule soit libre à l’instant est identifiée respectivement par un indice de probabilité et le modèle inverse libre à l’instant est identifié par un indice de modèle inverse libre . Comme l’illustre la figure 6, le filtre bayésien binaire entier de la grille d’espace libre calcule l’indice de probabilité de façon itérative, en fonction de l’indice de probabilité . La fonction d’initialisation du filtre dépend des paramètres d’entrée suivants : , et La constante et l’ensemble d’indices sont identiques à ceux utilisés pour le filtre bayésien binaire entier de la grille d’espace statique . est une constante, avec , et qui correspond à la valeur numérique du modèle de prédiction libre . Le détail des calculs exposé précédemment pour le filtre bayésien binaire entier de la grille d’espace statique s’applique, de la même façon, au filtre bayésien binaire entier de la grille d’espace libre ; il ne sont donc pas repris ici. Les entrées de la fonction d’application du filtre sont l’indice de probabilité qui est calculé récursivement, avec rebouclage direct depuis la sortie du filtre, et l’indice de modèle inverse libre . Avantageusement, l’indice de modèle inverse libre est calculé de la façon suivante (bloc « OG vers FIM », ) : Si la cellule est probablement libre à l'instant à partir de la grille d’occupation à l’instant , la valeur de l’indice de modèle inverse libre est déterminée par la fonction qui renvoie une valeur positive. L’indice du modèle inverse libre est égal à zéro dans les autres cas, c’est-à-dire : correspond à l’indice de probabilité dans la grille d’occupation de la cellule : Selon un mode de réalisation, la fonction h renvoie une valeur constante . En variante, on peut définir Ou encore, selon une autre variante : La génération d’une grille combinée résultant de la combinaison de la grille d’espace statique et de la grille d’espace libre peut être implémentée en utilisant les filtres bayésiens binaires entiers, comme l’illustre la figure 7, la combinaison se faisant, dans ce cas, entre les indices de probabilité et , de façon à obtenir un indice pour chaque cellule de la grille combinée. Dans ce cas, utilisant les filtres bayésiens binaires entiers, chaque cellule de la grille combinée peut être calculée par la fusion bayésienne suivante : Ainsi, selon le signe de : si , la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. Ainsi, on regarde l’estimation qui « l’emporte » entre l’estimation de la grille d’espace statique et l’estimation de la grille d’espace libre . En variante, il est possible d’appliquer un critère de valeur maximale. Ainsi, chaque cellule de la grille combinée peut être calculée de la façon suivante : Grâce à l’utilisation de filtres bayésiens binaires entiers, l’estimation de mouvement peut être effectuée dans des microprocesseurs, par exemple avec des architectures de type ARM (marque déposée), alors que dans l’art antérieur, l’estimation est très gourmande en ressources informatiques, et nécessite habituellement des processeurs graphiques (ou GPU pour « Graphics Processing Unit »). La illustre la complémentarité entre le procédé selon l’invention et celui décrit dans la demande de brevet FR3116640. Une pluralité de capteurs (capteur 1, …, capteur n). Les mesures sont fusionnées en mettant en œuvre le procédé tel que décrit dans la demande de brevet FR3041451, ce qui permet d’obtenir des grilles d’occupations OG (une probabilité d’occupation par cellule). Les grilles d’occupations peuvent être utilisées aussi bien pour générer les grilles d’espace statique et d’espace libre selon la présente invention, ou pour générer la grille de mouvement décrite dans la demande de brevet FR3116640. Néanmoins, les calculs pour aboutir à la grille de mouvement sont fondamentalement différents de ceux utilisés pour aboutir aux grilles d’espace statique et d’espace libre, comme indiqué précédemment. Le procédé selon l’invention peut être mis en œuvre par au moins un capteur de distance, par exemple un LIDAR, un RADAR ou encore un SONAR, ou par une caméra apte à extraire des informations de distance à partir d’une scène acquise par le capteur. Les grilles de caractéristiques sont établies à partir des mesures de distance effectuées par le capteur. Le capteur peut être embarqué à bord d’un porteur mobile, tel qu’un véhicule, en particulier un véhicule autonome, ou bien non autonome, et dans ce cas le procédé selon l’invention permet d’assister le conducteur pour la détection d’obstacles. Le porteur mobile peut également être tout dispositif apte à se mouvoir de façon autonome, par exemple un appareil domestique tel qu’un aspirateur autonome, ou encore un appareil de jardinage tel qu’une tondeuse à gazon autonome. Le porteur mobile comprend par ailleurs des actionneurs, aptes à corriger la trajectoire du porteur mobile autonome en fonction de la caractérisation de l’environnement d’un porteur mobile qui a été effectuée grâce au procédé selon l’invention. L’invention a été essentiellement décrite en discrétisant l’ensemble des probabilités dans un ensemble de cardinalité finie, afin de ne pas avoir à effectuer des calculs utilisant l'arithmétique en virgule flottante. Cette disposition permet en particulier de d’implémenter le procédé selon l’invention dans un système embarqué, ayant de fortes contraintes d’intégration. Le procédé selon l’invention pourrait également utiliser l'arithmétique en virgule flottante, en particulier si les capacités calculatoires du dispositif mettant en œuvre le procédé permettent d’effectuer la caractérisation de l’environnement en temps réel. Procédé de caractérisation de l’environnement d’un porteur mobile, dans lequel, pour chaque itération à un instant , les étapes suivantes sont mises en œuvre : S10) Acquisition d'une pluralité de mesures de distance dans l’environnement par au moins un capteur ; S20) Génération d’une paire de grilles d’occupation à l’instant et à l’instant t , chaque grille ( fusionnant les mesures de distance dans une représentation spatiale discrétisée de l’environnement ; S30) Génération d’une grille d’espace statique à l’instant t , chaque cellule de la grille d’espace statique à l’instant ayant une probabilité postérieure que ladite cellule soit occupée par un corps matériel statique à l’instant ( , ladite probabilité postérieure étant calculée pour chaque cellule par un filtre bayésien binaire à partir : de la probabilité postérieure que ladite cellule soit occupée par un corps matériel statique à l’instant t , ladite probabilité postérieure que ladite cellule soit occupée par un obstacle statique à l’instant étant directement injectée en entrée du filtre bayésien binaire ; d’un modèle de prédiction dit statique comprenant une probabilité que la cellule soit occupée par un corps matériel statique à l’instant en fonction de l’état de la cellule à l’instant ; et d’un modèle inverse dit statique à l’instant qui comprend une probabilité de connaissance de l’occupation de la cellule par un corps matériel statique à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant ; et/ou S40) Génération d’une grille d’espace libre à l’instant t , chaque cellule de la grille d’espace libre à l’instant ayant une probabilité postérieure que ladite cellule soit libre à l’instant , ladite probabilité postérieure étant calculée pour chaque cellule par un filtre bayésien binaire à partir : de la probabilité postérieure que ladite cellule soit libre à l’instant t , ladite probabilité postérieure que ladite cellule soit libre à l’instant étant directement injectée en entrée du filtre bayésien binaire ; d’un modèle de prédiction dit libre comprenant une probabilité que la cellule soit libre à l’instant en fonction de l’état de la cellule à l’instant ; et d’un modèle inverse dit libre à l’instant qui comprend une probabilité de connaissance du fait que la cellule soit libre à l’instant à partir de la paire de grilles d’occupation à l’instant et à l’instant . Procédé selon la revendication 1, dans lequel l’étape S30) de génération d’une grille d’espace statique et l’étape S40) de génération d’une grille d’espace libre sont mises en œuvre de façon concomitante. Procédé selon la revendication 2, comprenant une étape S50) de génération d’une grille combinée résultant de la combinaison de la grille d’espace statique et de la grille d’espace libre . Procédé selon la revendication 3, dans lequel chaque cellule de la grille combinée est calculée par la fusion bayésienne suivante : Dans lequel est la fonction de fusion bayésienne, correspond à une cellule de la grille combinée à l’instant , correspond à une cellule de la grille d’espace statique à l’instant , et correspond à une cellule de la grille d’espace libre à l’instant , Et selon la valeur de : si , la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. Procédé selon l’une des revendications 1 à 3, dans lequel : - la probabilité postérieure que la cellule soit occupée par un corps matériel statique à l’instant et le modèle inverse statique à l’instant sont approximés par des valeurs appartenant à un ensemble de cardinalité finie, les valeurs étant identifiées respectivement par un indice de probabilité et par un indice de modèle inverse statique ; - la probabilité postérieure que la cellule soit libre à l’instant et le modèle inverse libre à l’instant sont approximés par des valeurs appartenant à un ensemble de cardinalité finie, les valeurs étant identifiées respectivement par un indice de probabilité et par un indice de modèle inverse libre . Procédé selon la revendication 5, dans lequel l’indice de modèle inverse statique est calculé de la façon suivante : La fonction renvoyant une valeur positive, correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant , et correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant . Procédé selon l’une des revendications 5 ou 6, dans lequel l’indice de modèle inverse libre est calculé de la façon suivante : et correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant , la fonction renvoyant une valeur positive. Procédé selon les revendications 6 et 7, dans lequel la valeur positive de la fonction et la valeur positive de la fonction sont des valeurs constantes ( . Procédé selon la revendication 6, dans lequel la valeur positive renvoyée par la fonction est égale à . Procédé selon la revendication 7, dans lequel la valeur positive renvoyée par la fonction est égale à . correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant , et correspondant à l’indice d’occupation de la cellule de la grille d’occupation à l’instant . Procédé selon l’une des revendications 5 à 11, dans lequel l’indice correspondant à la composante de prédiction du filtre bayésien binaire de la grille d’espace statique est obtenu à l’aide d’une table de consultation comprenant un ensemble fini d’indices de probabilités, et l’indice filtré est obtenu par addition d’indices entiers de probabilités, et dans lequel l’indice correspondant à la composante de prédiction du filtre bayésien binaire de la grille d’espace libre est obtenu à l’aide d’une table de consultation comprenant un ensemble fini d’indices de probabilités, et l’indice filtré est obtenu par addition d’indices entiers de probabilités. 13. Procédé selon la revendication 5 à 12 en combinaison avec la revendication 3, dans lequel chaque cellule de la grille combinée est calculée par la fusion bayésienne suivante : correspond à une cellule de la grille combinée à l’instant , correspond à une cellule de la grille d’espace statique à l’instant , et correspond à une cellule de la grille d’espace libre à l’instant , Et selon la valeur de : si la cellule est probablement statique compte tenu de la séquence de paires de grilles d’occupation si , la cellule est probablement libre, si , la cellule n'est ni libre ni statique. Procédé d’évitement d’un corps matériel en mouvement autour d’un corps mobile, caractérisé en ce qu’il met en œuvre le procédé selon l'une des revendications 1 à 13 pour caractériser l’environnement d’un porteur mobile, et qu’il envoie une commande à un actionneur du porteur mobile afin d’éviter ledit corps matériel. Dispositif de caractérisation de l’environnement d’un porteur mobile, comprenant: - au moins un port d'entrée pour recevoir une pluralité de signaux représentatifs d'une série temporelle de mesures de distance issues d'un ou plusieurs capteurs de distance, et - un processeur de données configuré pour recevoir en entrée lesdits signaux, pour générer une grille d’espace libre ou une grille d’espace statique à partir desdits signaux en appliquant un procédé selon l'une des revendications 1 à 14.