FR 3114174 A1 20220318 FR 2009291 A 20200914 [Procédé de détection coopérative de corps matériels(O 1 , O 2 ) dans un environnement mettant en œuvre une pluralité de capteurs de distance équipant des porteurs respectifs (V 1 , V 2 , V 3 ), comprenant: A) pour chaque capteur, l’acquisition d’une série de mesures de distance d’un plus proche corps matériel selon une ligne de vue ; l’application d’un modèle inverse du capteur sur une grille d’occupation locale (GO 1 , GO 2 , GO 3 ) ; et la construction d’une grille d’occupation locale consolidée par fusion bayésienne des probabilités d’occupation ainsi déterminées ; B) le calcul des probabilités d’occupation des cellules d’une grille d’occupation globale (GO G ) par fusion bayésienne des probabilités d’occupation de cellules correspondantes d’au moins certaines desdites grilles d’occupation locales consolidées. Système pour la mise en œuvre d’un tel procédé. Figure pour l’abrégé : 12 Procédé et système de détection coopérative de corps matériels dans un environnement L’invention porte sur un procédé de détection coopérative de corps matériel dans un environnement, au moyen de capteurs équipant une pluralité de porteurs, dont au moins certains sont mobiles, qui évoluent dans un environnement. La détection est « coopérative » au sens que des données issues de plusieurs capteurs équipant des porteurs distincts sont fusionnées pour créer un modèle de l’environnement. Ce dernier est donc issu de la « coopération » des différents porteurs. L’invention s’applique d’une manière générale au domaine de la perception et préférentiellement à celui de la navigation de robots, drones, véhicules autonomes, etc. Il existe deux principales familles de techniques de perception : les méthodes géométriques, qui visent à identifier la géométrie des objets de l’espace environnant, et celles à base de grille d’occupation, qui visent à déterminer si un certain emplacement est occupé par un obstacle (plus généralement, par un corps matériel). L’invention relève des techniques à base de grille d’occupation. Les fondements théoriques des méthodes de perception basées sur des grilles d’occupation sont décrits dans l’article de (Elfes 1990). L’approche décrite dans cet article permet en particulier la fusion d’informations issues de plusieurs capteurs, qui peuvent ne pas être co-localisés. L’article (Rakotovao 2016) et la demande de brevet FR 3 041 451 divulguent une méthode permettant d’effectuer de façon efficace mais néanmoins précise, par utilisation d’une arithmétique entière, les calculs de fusion des probabilités d’occupation requis par la méthode de (Elfes 1990). Les articles (Coué 2006), (Pathak 2007), (Adarve 2012) et la thèse (Fernandez 2015) divulguent des perfectionnements de la méthode de (Elfes 1990) et son application à la fusion d’informations issues de plusieurs capteurs, homogènes ou hétérogènes, portés par un même véhicule. Ces méthodes permettent seulement d’exploiter des mesures issues de capteurs équipant un même porteur, par exemple un véhicule autonome. Or, dans de nombreuses applications, plusieurs porteurs équipés de capteurs évoluent dans une même région de l’espace ; il suffit de penser, par exemple à plusieurs véhicules automobiles à un carrefour. Conformément à l’art antérieur, ces véhicules calculent chacun une grille d’occupation couvrant une partie limitée de l’environnement. En outre, si l’environnement est dynamique et évolue à une cadence élevée, ou lorsque le porteur se déplace, chaque grille d’occupation devra être recalculée fréquemment pour représenter au mieux l’environnement autour du porteur. Si ce n’est pas le cas, on observera des incertitudes élevées sur l’état d’occupation des cellules de la grille. L’invention vise à surmonter, au moins en partie, ces inconvénients de l’art antérieur. Plus particulièrement elle vise à améliorer la qualité – en termes d’étendue spatiale et d’incertitude – de la perception de l’environnement dans le cas où plusieurs porteurs équipés de capteurs évoluent dans un environnement, et dont la position est connue au moins de manière approximative (dans le cas de véhicules automobiles, une estimation de la position peut être fournie, par exemple, par un système de positionnement par satellite). Conformément à l’invention, ce but est atteint par une approche coopérative, dans laquelle les grilles d’occupation « locales » correspondant aux mesures effectuées par les différents porteurs sont fusionnées en une grille d’occupation « globale ». Un aspect important de l’invention est que, lorsque des cellules de différentes grilles locales recouvrent une même zone de l’environnement, les probabilités d’occupation des cellules correspondantes de la grille globale sont calculées par fusion bayésienne des probabilités d’occupation des cellules des grilles locales. Cela permet de réduire l’incertitude sur l’état d’occupation des cellules de la grille locale, en utilisant de manière optimale toute l’information disponible. Le document (Birk 2006) décrit un procédé coopératif dans lequel plusieurs robots mobiles autonomes acquièrent conjointement une carte de l’environnement dans lequel ils évoluent. Dans ce procédé, les robots évoluent de manière indépendante, explorant différentes régions de l’environnement, avec un certain recouvrement, et en génèrent des cartes locales. Les robots ignorent leurs positions relatives, mais l’identification de régions identiques dans les cartes locales de différents robots permet de les recaler pour parvenir à une carte globale. Le procédé ne permet pas d’affiner de manière coopérative la détermination des probabilités d’occupation des cellules correspondant aux régions de recouvrement des cartes locales. En effet, la notion de probabilité d’occupation n’est pas exploitée : les cellules sont considérées soit occupées avec certitude, soit inoccupées avec certitude, soit dans un état totalement inconnu, auquel cas elles sont ignorées. L’environnement est considéré statique ou quasi-statique. En outre, le procédé ne permet pas d’exploiter une connaissance approximative de la position des différents robots à chaque instant. Le document EP 3 121 762 décrit un procédé dans lequel un véhicule fusionne les informations fournies par ses propres capteurs avec des données transmises par d’autres véhicules proches, également équipés de capteurs. Ce procédé coopératif n’utilise cependant pas de grilles d’occupation, mais une approche basée sur des caractéristiques (« feature-based » en anglais) qui ne permet pas d’aboutir à une véritable cartographie et nécessite un environnement structuré. Contrairement aux procédés utilisant des grilles d’occupation, cette approche ne permet donc pas de cartographier des environnements inconnus, et d’adresser donc la situation de porteurs évoluant dans un environnement inconnu. Un objet de l’invention est donc un procédé de détection coopérative de corps matériels dans un environnement mettant en œuvre une pluralité de capteurs de distance équipant des porteurs respectifs situés dans ledit environnement, au moins un desdits porteurs étant mobile, chaque capteur ayant une position et une orientation connues au moins de manière approchée avec un niveau d'incertitude connu, le procédé comprenant: A) pour chaque capteur, l'acquisition d'une série de mesures de distance d'un plus proche corps matériel selon une ligne de vue ; l'application, à chaque dite mesure de distance, d'un modèle inverse du capteur sur une grille d'occupation locale fournissant une représentation spatiale discrétisée de son champ de vue, pour déterminer une probabilité d'occupation par un corps matériel d'un ensemble de cellules de ladite grille d'occupation ; et la construction d'une grille d'occupation locale consolidée dont chaque cellule présente une probabilité d'occupation calculée par fusion bayésienne des probabilités d'occupation ainsi déterminées ; B) le calcul des probabilités d'occupation des cellules d'une grille d'occupation globale par fusion bayésienne des probabilités d'occupation de cellules correspondantes d'au moins certaines desdites grilles d'occupation locales consolidées. Selon des modes de réalisation particuliers d’un tel procédé : - L'étape B) peut comprendre, pour chaque cellule de la grille de probabilité globale: b1) l'identification d'un sous-ensemble de cellules de chaque grille d'occupation locale consolidée satisfaisant une condition de proximité avec ladite cellule de la grille d'occupation globale ; b2) le calcul d'un niveau de proximité de chaque cellule ainsi identifiée avec ladite cellule de la grille d'occupation globale ; b3) pour chaque grille d'occupation locale, le calcul d'une contribution à la probabilité d'occupation par moyenne pondérée des probabilités d'occupation des cellules identifiées lors de la sous-étape b1), les coefficients de pondération étant fonction des niveaux de proximité correspondants ; et b4) la fusion bayésienne des contributions aux probabilités d'occupation calculés lors de la sous-étape b3). - La sous-étape b1) peut comprendre l'identification des cellules de chaque grille d'occupation locale consolidée qui se trouvent en tout ou en partie à l'intérieur d'une région d'un rayon donné, centrée sur la cellule de la grille de probabilité globale. En particulier, ledit rayon peut être d'autant plus grand que l'incertitude sur la position du capteur correspondant à la grille d'occupation locale consolidée est élevée. - Le niveau de proximité calculé lors de la sous-étape b2) peut être représentatif de la densité de probabilité que le centre de la cellule identifiée coïncide avec le centre de la cellule de la grille d'occupation globale. - L'étape A) peut être mise en œuvre par une pluralité d'unités de traitement locales, associées chacune à un porteur et l'étape B) peut être mise en œuvre par une unité de traitement centrale, le procédé comprenant également la transmission des grilles d'occupation locales consolidées desdits capteurs à ladite unité de traitement centrale. - Le procédé peut comprendre une étape préalable de détermination d'une position de référence et d'une orientation définissant un repère global de référence, l'étape A) peut comprendre, pour chaque capteur, la construction d'une dite grille d'occupation locale consolidée alignée et localisée dans le repère global de référence, les valeurs de probabilité d'occupation de ladite grille d'occupation locale consolidée étant calculées en prenant en compte l'incertitude sur la position et l'orientation du capteur. Plus particulièrement, un capteur peut présenter une position et une orientation fixes et connues sans incertitude, prises comme position de référence et comme orientation du repère global de référence. En outre : - une unité de traitement centrale peut déterminer la position de référence et l'orientation définissant un repère global de référence et les transmettre à une pluralité d'unité de traitement locales, associée chacune à un porteur ; - l'étape A) peut être mise en œuvre par lesdites unités de traitement locales ; - les grilles d'occupation locales consolidées peuvent être transmises par les unités de traitement locales à l'unité de traitement centrale ; et - l'étape B) peut être mise en œuvre par ladite unité de traitement centrale. Chaque dit modèle inverse de capteur peut être un modèle discret, associant à chaque cellule de la grille d'occupation correspondante, et pour chaque mesure de distance, une classe de probabilité choisie à l'intérieur d'un même ensemble de cardinalité finie, chaque dite classe de probabilité étant identifiée par un indice entier ; les opérations de fusion bayésienne de probabilité étant effectuées au moyen de calculs entiers effectués sur les indices des classes de probabilités déterminées lors de ladite étape b). Un autre objet de l’invention est un système de détection coopérative de corps matériels comprenant : - une pluralité de capteurs de distance équipant des porteurs respectifs, au moins un desdits porteurs étant mobile ; - une pluralité de dispositifs de localisation équipant au moins ledit ou lesdits porteurs mobiles, lesdits dispositifs de localisation étant adaptés pour déterminer une position et une orientation du porteur, ainsi que leurs niveaux d'incertitude ; et - un système de traitement configuré ou programmé pour : - pour chaque capteur : acquérir une série de mesures de distance d'un plus proche corps matériel selon une ligne de vue ; appliquer, à chaque dite mesure de distance, un modèle inverse du capteur sur une grille d'occupation locale fournissant une représentation spatiale discrétisée de son champ de vue, pour déterminer une probabilité d'occupation par un corps matériel d'un ensemble de cellules de ladite grille d'occupation ; construire une grille d'occupation locale consolidée dont chaque cellule présente une probabilité d'occupation calculée par fusion bayésienne des probabilités d'occupation ainsi déterminées ; et - calculer des probabilités d'occupation des cellules d'une grille d'occupation globale par fusion bayésienne des probabilités d'occupation de cellules correspondantes d'au moins certaines desdites grilles d'occupation locales consolidées. Selon des modes de réalisation particuliers d’un tel système : - Le système de traitement peut comprendre une pluralité d'unités de traitement locales, associées chacune à un porteur, une unité de traitement centrale et un réseau de communication reliant l'unité de traitement centrale à chaque unité de traitement locale, et dans lequel : - chaque unité de traitement locale est configurée ou programmée pour construire une dite grille d'occupation locale consolidée et la transmettre à l'unité de traitement centrale ; et - l'unité de traitement centrale est configurée ou programmée pour calculer les probabilités d'occupation des cellules de ladite grille d'occupation globale. - Plus particulièrement, l'unité de traitement centrale peut être configurée ou programmée pour, pour chaque cellule de la grille d'occupation globale: - identifier un sous-ensemble de cellules de chaque grille d'occupation locale consolidée satisfaisant une condition de proximité avec ladite cellule de la grille d'occupation globale ; - calculer un niveau de proximité de chaque cellule ainsi identifiée avec ladite cellule de la grille d'occupation globale ; - pour chaque grille d'occupation locale, calculer une probabilité d'occupation par moyenne pondérée des probabilités d'occupation des cellules identifiées, les coefficients de pondération étant fonction des niveaux de proximité correspondants ; et - effectuer la fusion bayésienne des probabilités ainsi calculées. - L'unité de traitement centrale peut être configurée ou programmée pour déterminer une position de référence et une orientation définissant un repère global de référence, et les transmettre aux unités de traitement locales ; et chaque unité de traitement locale peut être configurée ou programmée pour construire une dite grille d'occupation locale consolidée alignée et localisée dans le repère global de référence, et calculer des valeurs de probabilité d'occupation de ladite grille d'occupation locale consolidée en prenant en compte l'incertitude sur la position et l'orientation du capteur respectif. Les dessins annexés illustrent l’invention : représente une grille d’occupation générée par un véhicule équipé d’un capteur de distance ; illustre la fusion de deux grilles d’occupation non alignées (définies dans des repères différents et avec des tailles de cellule différentes) conformément à un premier mode de réalisation de l’invention ; , illustrent la sélection de cellules de la grille d’occupation d’un véhicule « proches » d’une cellule d’une grille d’occupation globale, et la détermination de leur niveau de proximité, selon un deuxième mode de réalisation de l’invention ; est un organigramme d’un procédé selon ledit deuxième mode de réalisation de l’invention ; , illustrent les étapes d’un procédé selon un troisième mode de réalisation de l’invention ; et est un schéma de principe d’un système selon un quatrième mode de réalisation de l’invention. La montre une grille d’occupation GO associée à un porteur (ici un véhicule automobile) V. La grille GO est définie dans un repère solidaire du véhicule V et occupe une région de l’environnement qui correspond au moins en partie avec le champ de vue d’un ou plusieurs capteurs de distance embarqués sur le véhicule. La grille constitue une représentation spatiale discrétisée de l’environnement du véhicule, et plus précisément du champ de vue de son ou de ses capteurs de distance. Elle est constituée d’un ensemble de cellules adjacentes, chacune caractérisée par une probabilité d’occupation, c’est-à-dire par la probabilité que la cellule soit occupée par un obstacle. Ces probabilités (représentées sur la figure par des nuances de gris) sont calculées à partir des mesures issues du ou des capteurs, comme expliqué en détail dans l’article (Rakotovao 2016) et la demande de brevet FR 3 041 451 précités. Plus précisément : - Le (ou chaque) capteur embarqué sur le véhicule V est considéré « incertain ». Cela signifie que si un obstacle se trouve dans une région de l’environnement correspondant à une cellule c i de la grille d’occupation, alors la mesure de distance z issue du capteur suivra une distribution de probabilité p(z|o i ), où p(a|b) est la densité de probabilité de la variable aléatoire a conditionnée à l’évènement b et o i représente l’évènement « la région de l’environnement correspondant à la cellule c i de la grille d’occupation contient un obstacle » (par abus de langage on dira « la cellule c i est occupée »). La distribution de probabilité p(z|o i ) est appelée le « modèle direct » du capteur. - A partir du modèle direct du capteur et des probabilités d’occupation a priori des cellules P(o i ) on calcule le « modèle inverse » du capteur P(o i |z), qui exprime la probabilité que chaque cellule c i soit occupée, conditionnée à l’acquisition d’une mesure de distance z : Par exemple, on peut prendre des probabilités a priori P(o i ) identiquement égales à 1/2 , ce qui exprime une incertitude totale. Concrètement, seules les cellules situées sur la ligne de visée du capteur voient leur probabilité d’occupation modifiée par rapport à sa valeur a priori. Le calcul direct du modèle inverse d’un capteur pour une grille d’occupation déterminée présente une complexité exponentielle par rapport aux dimensions de la grille. L’article (Rakotovao 2016) et la demande de brevet FR 3 041 451 précités enseignent une méthode approchée beaucoup plus efficiente, présentant une complexité linéaire. - Le modèle inverse, qui est stocké dans une mémoire, est appliqué à chaque mesure issue du (ou de chaque) capteur. En général, plusieurs mesures sont disponibles, issues simultanément ou successivement d’un ou plusieurs capteurs ; il est alors possible de calculer la probabilité d’occupation de chaque cellule par fusion bayesienne de ces mesures. Si l’on considère deux mesures z 1 et z 2 on a par exemple : La généralisation à plus de deux capteurs est immédiate – il suffit de considérer comme le modèle inverse d’un capteur « virtuel » et de le fusionner avec la mesure fournie par un troisième capteur, et ainsi de suite. L’équation ci-dessus n’est valable que si P(o i ) = P(v i ) = 1/2, où v i correspond à « l’état de la cellule c i est vide », mais sa généralisation à d’autres hypothèses est sans difficulté, voir par exemple EP 3 364 213. D’une manière générale, la fusion de deux probabilités P 1 , P 2 est donnée par . On peut considérer que chaque mesure de distance z 1 , z 2, … permet de construire une grille d’occupation « provisoire », souvent en n’actualisant que les probabilités d’occupation des cellules situées sur une ligne de visée du capteur de distance, et que la fusion bayesienne de ces grilles conduit à une grille d’occupation « consolidée ». Dans la suite, par « grille d’occupation locale » on entendra, sauf mention contraire, la grille d’occupation consolidée pour un porteur. La fusion bayesienne nécessite des calculs en virgule flottante, couteux en ressources informatiques ; cependant l’article (Rakotovao 2016) et la demande de brevet FR 3 041 451 précités enseignent une méthode approchée beaucoup plus efficiente, nécessitant uniquement des calculs sur des variables entières. Une idée à la base de l’invention consiste à construire, de la manière décrite ci-dessus, des grilles d’occupation consolidées pour une pluralité de porteurs, chacun équipé d’un capteur de distance respectif (la généralisation au cas où tous ou certains des porteurs sont équipés de plusieurs capteurs ne pose pas de difficulté), puis de fusionner ces grilles d’occupation – que l’on peut qualifier de « locales » – pour construire une grille d’occupation « globale ». Par exemple, les porteurs peuvent être des véhicules automobiles évoluant sur un réseau routier ainsi que, éventuellement, des éléments d’infrastructure (fixes) tels des feux de signalisation équipés également de capteurs de distance. La illustre la fusion de deux grilles d’occupation locales GO1, GO2 pour engendrer une grille d’occupation globale GOG, donc seule une cellule cGk est représentée. Typiquement, les grilles d’occupation locales sont calculées par des ordinateurs embarqués sur des porteurs respectifs et transmises par des moyens de communication – par exemple radio – à un ordinateur « central », avec des informations sur la position et l’orientation du porteur (obtenues, par exemple, au moyen d’un système de navigation satellitaire). L’ordinateur central – qui peut être dans une station fixe ou être embarqué sur un porteur équipé de capteurs, qui engendre sa propre grille d’occupation locale – calcule la grille d’occupation globale. On remarque que, en général, les grilles d’occupation locales ne sont pas alignées car elles présentent des orientations différentes, car elles sont définies dans des référentiels attachés aux porteurs qui se déplacent de manière indépendante. Par ailleurs, comme dans le cas de la figure, elles peuvent présenter des cellules de taille différente. La grille d’occupation globale est définie dans un repère dit de référence. Selon un premier mode de réalisation de l’invention, illustré sur la , on identifie à l’intérieur de la cellule cGk de la grille d’occupation globale plusieurs régions correspondant à des superpositions de différentes cellules des grilles locales GO1, GO2. Ces régions sont numérotées de 1 à 21 sur la figure. Pour chacune de ces régions on détermine une probabilité d’occupation par fusion bayesienne des probabilités d’occupation des cellules locales considérées (pour mémoire, une cellule est considérée « occupée » si un obstacle a été détecté dans cette cellule, sans nécessairement que l’obstacle occupe l’intégralité de la cellule). Puis, on calcule une probabilité d’occupation de la cellule cGk à partir des probabilités d’occupation de ces régions. Par exemple, on peut considérer que la cellule toute entière prend la probabilité d’occupation la plus élevée parmi celles calculées pour ses régions. L’inconvénient de cette approche est sa complexité computationnelle : dans l’exemple de la , une opération de fusion bayesienne est nécessaire pour calculer la probabilité d’occupation de chaque région de la cellule cGk de la grille globale, soit 21 opérations de fusion bayesienne à calculer pour cGk. Il faut ensuite utiliser un mécanisme qui calcule une valeur unique de probabilité d’occupation pour la cellule cGk, par exemple avec une fonction maximum. Et le nombre d’opérations augmente de manière exponentielle avec celui des grilles d’occupation locales que l’on souhaite fusionner. Un deuxième et un troisième mode de réalisation de l’invention permettent de surmonter cette difficulté. Le calcul d’une grille d’occupation globale selon un deuxième mode de réalisation est illustré par les à [Fig. 6]. Il comprend les étapes suivantes, mises en œuvre pour chaque cellule cGk de la grille d’occupation globale : Premièrement, on identifie des cellules de chaque grille d’occupation locale (consolidée) qui peuvent être considérées « proches » - par application d’un critère de proximité prédéfini – avec ladite cellule de la grille d’occupation globale. Seules ces cellules contribuent au calcul de la probabilité d’occupation de cGk. Par exemple, comme illustré sur les et [Fig. 4], on peut considérer « proches » de la cellule cGk les cellules de la grille d’occupation locale GO1 associée au porteur (véhicule) V1 qui se trouvent au moins en partie à l’intérieur d’un cercle CP1 de rayon R1 dont le centre coïncide avec celui de la cellule cGk. Avantageusement, le rayon R1 est d’autant plus grand que l’incertitude sur la position de V1 est élevée (les systèmes de navigation satellitaire, par exemple, permettent non seulement d’estimer une position, mais également de connaître l’incertitude affectant cette estimation). La [Fig. 3] montre le disque de frontière CP1 et la [Fig. 4] montre la région RPi1 de la grille d’occupation GO1 formée par les cellules qui intersectent le disque de frontière CP1 et sont donc considérées proches de cGk. Ensuite, pour chaque cellule de RP i 1 , on calcule un niveau de proximité avec la cellule c G k . Il convient de noter que le critère de proximité utilisé pour identifier les cellules faisant partie de RP i 1 est binaire – il est satisfait ou ne l’est pas – tandis que le niveau de proximité peut prendre une pluralité de valeurs discrètes ou continues. Par exemple, le niveau de proximité peut être la densité de probabilité normalisée, calculée au centre de la cellule identifiée, ce qui dénote une valeur de probabilité pour que cette cellule coïncide avec le centre de la cellule de la grille d’occupation globale, ou une fonction bornée entre 0 et 1 (typiquement monotone) de cette probabilité. En effet, étant donné que la position et l’orientation du porteur V 1 ne sont connues qu’avec une certaine imprécision, la position par rapport à c G k du centre de chaque cellule de GO 1 est décrite par une distribution de probabilité. Puis, une contribution de la grille d’occupation locale GO 1 à la probabilité d’occupation de la cellule c G k est calculée par moyenne pondérée des probabilités d’occupation des cellules de RP i 1 , les coefficients de pondération étant fonction des niveaux de proximité correspondants (un coefficient est d’autant plus grand que la cellule de RP 1i correspondante est proche de c G k ; par exemple, il peut dépendre de l’inverse du niveau de proximité). Enfin, les différentes contributions provenant de différentes grilles d’occupation locale sont fusionnées (fusion bayesienne) pour trouver la probabilité d’occupation de la cellule c G k . Pour cela, on peut utiliser par exemple les méthodes décrites dans (Coué 2006), (Pathak 2007), (Adarve 2012), (Fernandez 2015) ou encore, en appliquant avantageusement la technique de « fusion entière » divulguée par (Rakotovao et al. 2016) et WO2017/050890. Au lieu de procéder cellule par cellule de la grille d’occupation globale, il est possible d’associer à chaque grille d’occupation locale GO j une grille d’occupation locale intermédiaire, dont la géométrie (position, orientation, nombre, forme et agencement des cellules) est identique à celle de la grille d’occupation globale. La probabilité d’occupation de chaque cellule de la grille d’occupation locale intermédiaire est calculée à partir des probabilités d’occupation des cellules de GO j , comme décrit plus haut. En d’autres termes, on projette les grilles d’occupation locales dans le repère de référence de la grille d’occupation globale. Enfin, toutes les grilles d’occupation locale intermédiaires sont fusionnées pour aboutir à la grille d’occupation globale, ce qui ne pose pas de difficulté car ces grilles sont définies dans le même repère. Les deux approches sont similaires et diffèrent seulement par l’ordre de réalisation des opérations. Dans un cas, on calcule la probabilité d’occupation d’une cellule c g k de la grille d’occupation, en prenant en compte toutes les contributions à la fois. Dans l’autre cas, d’abord on calcule les contributions de chaque grille d’occupation locale à toutes les cellules de la grille d’occupation globale, puis on fusionne ces contributions. Les calculs nécessités par la mise en œuvre d’un procédé selon ce deuxième mode de réalisation de l’invention seront maintenant décrits plus en détails à l’aide de la . Soit un porteur V j positionné en sur la grille globale et une cellule c i (sélectionnée au préalable selon la procédure décrite précédemment) de sa grille (locale) de repère j associé au porteur V j , positionnée en par rapport au repère de référence RR . Les deux repères diffèrent par une rotation d’un angle j . On considère que la position d’un porteur est définie par la position d’un point à l’intérieur dudit porteur, par exemple son barycentre, et celle d’une cellule est définie par son centre (en variante, on pourrait aussi choisir un de ses vertex). La relation entre les positions du porteur V j et de la cellule c i dans le repère de référence est donnée par les équations où est l’estimation de la position du porteur qui peut être modélisée par des variables aléatoires (par exemple des variables gaussiennes de moyenne et variances , et est l’estimation de l’orientation de la droite reliant la position du porteur avec le centre de la cellule c i dans le repère de référence qui peut également être modélisée par une variable aléatoire (par exemple une variable gaussienne tronquée de moyenne de variance ). La distance entre le porteur et la cellule c i est considérée constante, car la grille locale qui contient la cellule c i est attachée au véhicule. A partir des équations ci-dessus et en connaissant les distributions de probabilité de et on peut déterminer la distribution de probabilité conjointe des variables à savoir, Pour déterminer on peut utiliser différentes approches : histogramme (numérique), fonction caractéristique, matrice jacobienne, formes différentielles. Si l’on considère par exemple que l’estimation de la position du porteur est modélisée par une variable aléatoire gaussienne de densité de probabilité avec et indépendantes et , l’orientation du véhicule suit une loi de probabilité gaussienne tronquée donnée par En utilisant les formes différentielles, soit l’élément de probabilité où et est le produit extérieur. On peut déduire que en prenant en compte la propriété (i.e. le produit extérieur est alterné). Ainsi, l’élément de probabilité peut être réécrit comme : avec indépendant de et Par conséquence, où L’intégrale peut être résolue numériquement. La fonction permet de quantifier le niveau de proximité entre la cellule c i de la grille GO i et la cellule de la grille globale GO G (en noir). Plus précisément on définit le niveau de proximité par la valeur où est la position de la cellule cGk dans repère de référence. La mesure dénote une probabilité pour laquelle le centre de la cellule ci de la grille GOi coïncide avec celui de la cellule cGk de la grille globale GOG. On détermine le niveau de proximité pour chaque cellule sélectionnée ci aux alentours de chaque cellule cGk du repère de référence de porteur Vj , comme cela a été expliqué plus haut en référence à la et à la [Fig. 4]. Dans le cas d’un porteur fixe (élément d’infrastructure) dont la position est connue avec certitude, on définit directement pour toutes les cellules sélectionnées. Dans le cas d’un porteur mobile, la normalisation de par le est nécessaire car peut être supérieur à 1 ; sans normalisation on arriverait à la situation paradoxale d’un niveau de proximité supérieur à celui atteignable dans le cas d’un porteur fixe. Une fois les niveaux de proximité calculés, on procède avec la projection de la probabilité d’occupation dans le repère d’un porteur V j sur le repère de référence. Soit la probabilité d’occupation de la cellule c i du repère , déterminée à partir des capteurs de distances installés sur le porteur ; on définit alors la probabilité d’occupation de la cellule du repère de référence pour le porteur V j par : La probabilité est donc une moyenne pondérée des valeurs de probabilité d’occupation des cellules sélectionnées préalablement dans la grille GO i dont les poids sont les mesures de proximité de position de ces cellules à l’égard de la cellule c Gk du repère de référence. Avec cette procédure, on arrive à construire les grilles d’occupation de tous les porteurs dans un même repère avec une complexité de calcul maîtrisée. En particulier, pour réduire le coût de calcul, le résultat de certains calculs intermédiaires pourrait être sauvegardé dans une table. A ce point, la fusion de ces grilles d’occupation ne pose pas de difficulté particulière et peut être réalisée en appliquant avantageusement la technique de « fusion entière » divulguée par (Rakotovao et al. 2016) et WO2017/050890. La représente un organigramme d’un procédé selon le deuxième mode de réalisation de l’invention. Le procédé de la figure se rapporte à une situation où M porteurs Vj (identifiés par l’indice j, j prenant les valeurs de 1 à M) fournissent chacun une grille d’occupation locale. La grille d’occupation globale comprend N cellules c. Pour chaque porteur (boucle d’itération sur j), on projette la grille d’occupation locale correspondante sur une grille d’occupation locale intermédiaire ayant la même géométrie que la grille d’occupation globale (boucle d’itération imbriquée sur chaque cellule c) ; puis on procède à la fusion des grilles d’occupation locale intermédiaires. Comme expliqué plus haut, en variante, on pourrait aussi procéder différemment en considérant les cellules de la grille d’occupation globale une par une et, pour chacune d’entre elles, en calculant et fusionnant les contributions à la probabilité d’occupation provenant de toutes les grilles locales d’occupation. Dans le deuxième mode de réalisation de l’invention qu’on vient de décrire, chaque porteur (ou plus précisément l’ordinateur embarqué dans chaque porteur) génère une grille d’occupation locale dans son propre repère et la transmet à l’ordinateur central, qui projette toutes ces grilles locales dans un repère de référence pour les fusionner et obtenir la grille d’occupation globale. Dans le troisième mode de réalisation, l’orientation et la position du repère de référence sont transmises par l’ordinateur central aux porteurs, qui génèrent alors directement leurs grilles d’occupation locales dans ce repère. Les grilles d’occupation locales ainsi générées peuvent ensuite être envoyées à l’ordinateur central qui les fusionne. La montre trois porteurs mobiles (voitures) V1, V2, V3 et un porteur fixe (élément d’infrastructure) EI avec des grilles d’occupation locales respectives GO1, GO2, GO3 GOEI « vides », c’est-à-dire avant que les probabilités d’occupation de leurs cellules n’aient été calculées. Deux obstacles O1, O2 sont également présents. Un ordinateur central, par exemple co-localisé avec l’élément d’infrastructure EI, transmet aux ordinateurs embarqués sur les porteurs mobiles (et, le cas échéant, sur l’élément d’infrastructure) les données qui définissent un repère global : une origine et les orientations des axes xG, yG. Une grille d’occupation globale GOG – elle aussi « vide » - est représenté sur la figure par un quadrillage. Chaque porteur mobile détermine, par exemple au moyen d’un système de navigation satellitaire, sa position par rapport à l’origine du repère global, ainsi que l’incertitude affectant sa position (pour un porteur statique tel EI, la position est supposée connue sans incertitude). Une distribution de probabilité de la position du porteur V i numéro « i » p(x i , y i ) peut donc être calculée, par exemple en faisant l’hypothèse que cette distribution suit une loi normale multidimensionnelle. L’ordinateur embarqué dans chaque porteur reçoit une succession de mesures de distance et s’en sert pour calculer la grille d’occupation locale, en utilisant le modèle inverse du capteur de distance et la probabilité de position p(x i , y i ). On tient donc compte de l’incertitude sur la position. Soit z t-1 la suite de mesures de distance {z 1 , …z t-1 } acquises avant l’instant t et z la mesure acquise au temps t. On désigne par P (o j |z, z t-1 ) la probabilité que la cellule « j » de la grille d’occupation soit occupée, conditionnée par les mesures z et z t-1 . On démontre que : avec où est le modèle inverse du capteur, calculé en supposant que la position du porteur soit (x i , y i ), P(x i , y i ) est la probabilité que la position du porteur soit effectivement (x i , y i ) et la somme est calculée sur toutes les positions possibles (on considère ici que x i et y i prennent seulement des valeurs discrètes, ce qui justifie le passage d’une densité de probabilité – représentée par un « p » minuscule – à une probabilité représentée par un « P » majuscule). On a considéré ici uniquement l’incertitude sur la position du porteur, mais l’incertitude sur son orientation peut être prise en considération de la même manière (on peut parler de « pose » pour désigner à la fois la position et l’orientation d’un porteur). Les , [Fig. 9], [Fig. 10] et [Fig. 11] montrent les grilles d’occupation locale GO1, GO2, GO3 et GOEI des porteurs V1, V2, V3 et EI, respectivement. On remarque, dans le cas de V3, que les cellules de la grille d’occupation locale sont alignées sur le repère de référence, et non sur celui du porteur. La [Fig. 12] montre la grille d’occupation globale GOG obtenue par fusion. Quel que soit le mode de réalisation considéré, la grille d’occupation globale peut être transmise par l’ordinateur central aux ordinateurs embarqués. L’ordinateur central, par ailleurs, n’est pas nécessairement « fixe » : il est possible d’utiliser en tant qu’ordinateur central un ordinateur embarqué sur un porteur. L’invention a été décrite en considérant des grilles d’occupation bidimensionnelles, mais la généralisation au cas tridimensionnel ne pose pas de difficulté. L’invention a été décrite en considérant le cas où chaque porteur est équipé d’un seul capteur de distance, mais la généralisation au cas où au moins certains porteurs sont équipés de plusieurs capteurs ne pose pas de difficulté. Bien que cela ne soit pas essentiel, il est avantageux d’appliquer la technique de « fusion entière » divulguée par (Rakotovao et al. 2016) et WO2017/050890, qui est brièvement rappelée ci-après. L’idée à la base de cette technique est que, bien que les probabilités à fusionner puissent en principe prendre n’importe quelle valeur dans l’intervalle [0, 1], en pratique les calculs sont nécessairement effectués sur des valeurs discrètes. Par conséquent, une probabilité dans l’intervalle [0, 1] doit nécessairement être représentée de manière discrétisée, par des « classes de probabilité » qui peuvent être identifiées par des indices entiers. Un sous-ensemble dénombrable de l’intervalle [0, 1] dont les éléments p n sont identifiés par un indice entier relatif « n » est appelé un « système de classes de probabilités » . Certains systèmes de classes de probabilités ont la propriété que la fusion de deux classes de probabilité quelconques du système appartient, elle aussi, au système. Formellement , . On parle dans ce cas de « systèmes de classes sans erreurs » car l’opération de fusion n’introduit aucune erreur ou approximation (la discrétisation, cependant, est une approximation). L’intérêt de ces classes est que le résultat de la fusion est également identifié par un indice entier. Le problème de la fusion de probabilité se réduit donc à déterminer une fonction F d qui associe à une paire d’indices entiers un autre indice entier. Des opérations en virgule flottante sur des valeurs de probabilité sont donc remplacées par des opérations entières sur des indices. Un schéma de discrétisation avantageux est basé sur un ensemble dénombrable de classes de probabilité S k donné par l’union des deux sous-ensembles suivants : où k est un entier positif ( ). S k discrétise l’intervalle [0, 1] de manière non uniforme, avec un intervalle de discrétisation qui est relativement grand autour de la valeur 0,5 et de plus en plus petit lorsqu’on s’approche des valeurs extrêmes 0 et 1. Cela est doublement avantageux ; en effet : - Une probabilité de 0,5 correspond à une incertitude totale ; il n’est donc pas très utile d’être précis autour d’une telle valeur ; par contre, il est utile d’être précis autour des valeurs extrêmes qui correspondent à des quasi-certitudes de non-occupation (0) ou d’occupation (1). - au-delà d’une certaine valeur de |n|, les valeurs discrètes de probabilité deviennent extrêmement proches les unes des autres ; l’erreur introduite en tronquant l’ensemble est donc négligeable. Bien que S k ne soit pas un système sans erreurs, Sk+ et Sk- le sont. Un autre système de quantification utile, sans erreurs, est défini par récursion : p 0 =0,5 ; p 1 =p ; p 2 =F(p,p); p 3 =F(p 2 ,p); … p n+1 =F(p n ,p) p -1 =1-p ; p -2 =F(p -1 , p -1 ); p -3 =F(p -2 , p -1 ); … p n-1 =F(p n , p -1 ) où “p” est un paramètre tel que 0,5 . La représente très schématiquement un système pour la mise en œuvre d’un procédé selon l’invention. Le système comprend des capteurs de distance CD1, CD2, des dispositifs de localisation L1, L2 et des ordinateurs (plus généralement des unités de traitement locales) UT1, UT2, équipant des porteurs respectifs V1, V2. Les capteurs de distance et les dispositifs de localisation fournissent à l’unité de traitement du porteur les mesures de distance et informations de position et orientation nécessaires pour construire les grilles d’occupation locales. Les unités de traitement locales échangent des données avec une unité de traitement centrale UTC par l’intermédiaire d’un réseau de communications RCM, typiquement sans fil. L’unité de traitement centrale UTC construit une grille d’occupation globale à partir des grilles d’occupation locale et de la position des porteurs V1 et V2. Les unités de traitements locales et globale peuvent être des ordinateurs, des microprocesseurs ou microcontrôleurs, des circuits numériques dédiés ou programmables etc., programmés ou configurés pour mettre en œuvre les traitements décrits plus haut. Les dispositifs de localisation peuvent être ou comprendre des récepteurs de signaux de navigation satellitaire, des centrales inertielles etc. Références bibliographiques Elfes 1990 : A. Elfes, “Occupancy grids: a stochastic spatial representation for active robot perception” (Sixth Conference on Uncertainty in AI, 1990). Rakotovao 2016 : T Rakotovao, J Mottin, D Puschini, C Laugier, “Multi-sensor fusion of occupancy grids based on integer arithmetic”, Robotics and Automation (ICRA), 2016 IEEE International Conference on, 1854-1859. Coué 2006 : Coué, C., Pradalier, C., Laugier, C., Fraichard, T. and Bessière, P. “Bayesian Occupancy Filtering for Multitarget Tracking: An Automotive Application”. The International Journal of Robotics Research, Janvier 2006. Pathak 2007 : Pathak, K., Birk, A., Poppinga, J., and Schwertfeger, S., “3D forward sensor modeling and application to occupancy grid based sensor fusion”. IEEE/RSJ International Conference on Intelligent Robots and Systems, Décembre 2007. Adarve 2012 : Adarve, J. D., Perrollaz, M., Makris, A. and Laugier, C. “Computing Occupancy Grids from Multiple Sensors using Linear Opinion Pools”. IEEE International Conference on Robotics and Automation, Mai 2012. Fernandez 2015 : Fernandez, C. G. D. P., ”Grid-based multi-sensor fusion for on-road obstacle detection: application to autonomous driving”. Thèse de Master, KTH-CVAP et Volvo Car Corporation, Août 2015. Birk 2006 : Birk, A, and Carpin, S.; “Merging occupancy grid maps from multiple robots”; Proceedings of the IEEE, s vol. 94, no 7, pp. 1384-1397, Juillet 2006. Procédé de détection coopérative de corps matériels (O 1 , O 2 ) dans un environnement mettant en œuvre une pluralité de capteurs de distance (CD 1 , CD 2 ) équipant des porteurs respectifs (V 1 , V 2 ) situés dans ledit environnement, au moins un desdits porteurs étant mobile, chaque capteur ayant une position et une orientation connues au moins de manière approchée avec un niveau d’incertitude connu, le procédé comprenant: A) pour chaque capteur, l’acquisition d’une série de mesures de distance d’un plus proche corps matériel (O 1 , O 2 ) selon une ligne de vue ; l’application, à chaque dite mesure de distance, d’un modèle inverse du capteur sur une grille d’occupation locale fournissant une représentation spatiale discrétisée de son champ de vue, pour déterminer une probabilité d’occupation par un corps matériel d’un ensemble de cellules de ladite grille d’occupation ; et la construction d’une grille d’occupation locale consolidée (GO 1 , GO 2 ) dont chaque cellule présente une probabilité d’occupation calculée par fusion bayésienne des probabilités d’occupation ainsi déterminées ; B) le calcul des probabilités d’occupation des cellules d’une grille d’occupation globale (GO G ) par fusion bayésienne des probabilités d’occupation de cellules correspondantes d’au moins certaines desdites grilles d’occupation locales consolidées. Procédé selon la revendication 1 dans lequel l’étape B) comprend, pour chaque cellule (c Gi ) de la grille de probabilité globale: b1) l’identification d’un sous-ensemble de cellules (RP i1 ) de chaque grille d’occupation locale consolidée (GO 1 ) satisfaisant une condition de proximité avec ladite cellule de la grille d’occupation globale ; b2) le calcul d’un niveau de proximité de chaque cellule ainsi identifiée avec ladite cellule de la grille d’occupation globale ; b3) pour chaque grille d’occupation locale, le calcul d’une contribution à la probabilité d’occupation par moyenne pondérée des probabilités d’occupation des cellules identifiées lors de la sous-étape b1), les coefficients de pondération étant fonction des niveaux de proximité correspondants ; et b4) la fusion bayésienne des contributions aux probabilités d’occupation calculés lors de la sous-étape b3). Procédé selon la revendication 2 dans lequel la sous-étape b1) comprend l’identification des cellules de chaque grille d’occupation locale consolidée qui se trouvent en tout ou en partie à l’intérieur d’une région d’un rayon donné (R 1 ), centrée sur la cellule de la grille de probabilité globale. Procédé selon la revendication 3 dans lequel ledit rayon est d’autant plus grand que l’incertitude sur la position du capteur correspondant à la grille d’occupation locale consolidée est élevée. Procédé selon l’une quelconque des revendications 2 à 4 dans lequel le niveau de proximité calculé lors de la sous-étape b2) est représentatif de la densité de probabilité que le centre de la cellule identifiée coïncide avec le centre de la cellule de la grille d’occupation globale. Procédé selon l’une quelconque des revendications précédentes dans lequel l’étape A) est mise en œuvre par une pluralité d’unités de traitement locales (UT 1 , UT 2 ), associée chacune à un porteur (V 1 , V 2 ) et l’étape B) est mise en œuvre par une unité de traitement centrale (UTC), le procédé comprenant également la transmission des grilles d’occupation locales consolidées desdits capteurs à ladite unité de traitement centrale. Procédé selon la revendication 1 comprenant une étape préalable de détermination d’une position de référence et d’une orientation définissant un repère global de référence (RR), dans lequel l’étape A) comprend, pour chaque capteur, la construction d’une dite grille d’occupation locale consolidée alignée et localisée dans le repère global de référence, les valeurs de probabilité d’occupation de ladite grille d’occupation locale consolidée étant calculées en prenant en compte l’incertitude sur la position et l’orientation du capteur. Procédé selon la revendication 7 dans lequel un capteur présente une position et une orientation fixes et connues sans incertitude, prises comme position de référence et comme orientation du repère global de référence. Procédé selon l’une quelconque des revendications 7 ou 8 dans lequel : - une unité de traitement centrale (UTC) détermine la position de référence et l’orientation définissant un repère global de référence (RR) et les transmets à une pluralité d’unité de traitement locales (UT 1 , UT 2 ), associées chacune à un porteur ; - l’étape A) est mise en œuvre par lesdites unités de traitement locales ; - les grilles d’occupation locales consolidées sont transmises par les unités de traitement locales à l’unité de traitement centrale ; et - l’étape B) est mise en œuvre par ladite unité de traitement centrale. Procédé selon l’une des revendications précédentes, dans lequel chaque dit modèle inverse de capteur est un modèle discret, associant à chaque cellule de la grille d’occupation correspondante, et pour chaque mesure de distance, une classe de probabilité choisie à l’intérieur d’un même ensemble de cardinalité finie, chaque dite classe de probabilité étant identifiée par un indice entier ; les opérations de fusion bayésienne de probabilité étant effectuées au moyen de calculs entiers effectués sur les indices des classes de probabilités déterminées lors de ladite étape b). Système de détection coopérative de corps matériels comprenant : - une pluralité de capteurs de distance (CD 1 , CD 2 ) équipant des porteurs respectifs (V 1 , V 2 ), au moins un desdits porteurs étant mobile ; - une pluralité de dispositifs de localisation (L 1 , L 2 ) équipant au moins ledit ou lesdits porteurs mobiles, lesdits dispositifs de localisation étant adaptés pour déterminer une position et une orientation du porteur, ainsi que leurs niveaux d’incertitude ; et - un système de traitement (UT 1 , UT 2 , UTC) configuré ou programmé pour : - pour chaque capteur : acquérir une série de mesures de distance d’un plus proche corps matériel (O 1 , O 2 ) selon une ligne de vue ; appliquer, à chaque dite mesure de distance, un modèle inverse du capteur sur une grille d’occupation locale fournissant une représentation spatiale discrétisée de son champ de vue, pour déterminer une probabilité d’occupation par un corps matériel d’un ensemble de cellules de ladite grille d’occupation ; construire une grille d’occupation locale consolidée (OG 1 , OG 2 ) dont chaque cellule présente une probabilité d’occupation calculée par fusion bayésienne des probabilités d’occupation ainsi déterminées ; et - calculer des probabilités d’occupation des cellules d’une grille d’occupation globale (OG G ) par fusion bayésienne des probabilités d’occupation de cellules correspondantes d’au moins certaines desdites grilles d’occupation locales consolidées. Système selon la revendication 11 dans lequel le système de traitement comprend une pluralité d’unités de traitement locales (UT 1 , UT 2 ), associées chacune à un porteur, une unité de traitement centrale (UTC) et un réseau de communication (RCM) reliant l’unité de traitement centrale à chaque unité de traitement locale, et dans lequel : - chaque unité de traitement locale est configurée ou programmée pour construire une dite grille d’occupation locale consolidée et la transmettre à l’unité de traitement centrale ; et - l’unité de traitement centrale est configurée ou programmée pour calculer les probabilités d’occupation des cellules de ladite grille d’occupation globale. Système selon la revendication 12 dans lequel l’unité de traitement centrale (UTC) est configurée ou programmée pour, pour chaque cellule (c Gi ) de la grille d’occupation globale: - identifier un sous-ensemble (RP i1 ) de cellules de chaque grille d’occupation locale consolidée satisfaisant une condition de proximité avec ladite cellule de la grille d’occupation globale ; et - calculer un niveau de proximité de chaque cellule ainsi identifiée avec ladite cellule de la grille d’occupation globale ; - pour chaque grille d’occupation locale, calculer une probabilité d’occupation par moyenne pondérée des probabilités d’occupation des cellules identifiées, les coefficients de pondération étant fonction des niveaux de proximité correspondants ; et - effectuer la fusion bayésienne des probabilités ainsi calculées. Système selon la revendication 12 dans lequel l’unité de traitement centrale (UTC) est configurée ou programmée pour déterminer une position de référence et une orientation définissant un repère global de référence (RR), et les transmettre aux unités de traitement locales (UT 1 , UT 2 ) ; et chaque unité de traitement locale est configurée ou programmée pour construire une dite grille d’occupation locale consolidée alignée et localisée dans le repère global de référence, et calculer des valeurs de probabilité d’occupation de ladite grille d’occupation locale consolidée en prenant en compte l’incertitude sur la position et l’orientation du capteur respectif.