Brevet FR3022049A1 — publié le 20151211

Procédé d'analyse d'une scène dynamique, module d'analyse et programme d'ordinateur associés La présente invention concerne un procédé d'analyse d'une scène dynamique observée à l'aide d'un bloc de capteur(s) comprenant une étape de définition d'une grille découpée en cellules et correspondant à la scène observée ; et la mise en oeuvre itérative par ordinateur des étapes : collecter au moins une nouvelle observation du bloc de capteur(s) à la kème itération et déterminer, en fonction de la nouvelle observation collectée, une première probabilité d'occupation de chaque cellule de la grille modélisant le fonctionnement du bloc de capteur(s) ; déterminer, pour chaque cellule, à la kème itération, une deuxième probabilité d'occupation de la cellule et un ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération. Selon les algorithmes de perception Bayésienne BOF (en anglais « Bayesian Occupancy Filter »), la scène observée est représentée par une grille d'occupation subdivisée en cellules et est analysée en termes de cellules occupées et de leur évolution dans le temps. Pour chaque cellule de cette grille d'occupation et à chaque itération de l'algorithme, une probabilité d'occupation de la cellule est calculée, ainsi qu'une distribution de probabilités de mouvement du contenu de la cellule. Cette distribution est représentée sous forme d'une grille de distribution de mouvement, encore nommé histogramme de transition de voisinage. Sur la figure 1, une vue d'une grille d'occupation de l'environnement G,' est représentée, de 16 cellules sur 16 cellules. Dans cette représentation, plus les hachures couvrant une cellule sont rapprochées, plus la valeur de la probabilité d'occupation de la cellule est grande. Un histogramme de transition Hist(C) déterminé pour une cellule, référencée C, de la grille Oc, est aussi représenté en figure 1. Chaque case de l'histogramme déterminé pour la cellule C représente, pour le pas de temps considéré, la transition de la cellule C vers ladite case ; à chaque case de l'histogramme est ainsi associé un vecteur vitesse de valeur distincte ; dans le cas présent, la cellule au centre de l'histogramme correspond à la transition de la cellule C vers elle-même. Ces algorithmes ont l'avantage de fournir une analyse fiable et synthétique d'une scène dynamique, mais présentent un certain nombre de limitations. Dans le cas d'une grille en mouvement, par exemple correspondant à l'environnement analysé depuis un véhicule en mouvement, des problèmes de crénelages (en anglais « aliasing ») spatio-temporels surgissent. En outre, il est nécessaire d'utiliser des grilles d'occupation et des histogrammes de mouvement à haute résolution pour détecter les mouvements lents et les petits objets, et aussi d'évaluer la vitesse d'objets très dynamiques (par exemple les autres véhicules en mouvement). Typiquement, si la norme de la vitesse maximale du véhicule est de Vmax, la taille de la grille de distribution de mouvement maintenue au cours de l'algorithme doit permettre de représenter des vitesses de -2V'x à + 2Vmax : chaque itération de l'algorithme implique donc d'estimer la probabilité de distribution des vitesses dans chaque case de cette grille pour chaque cellule. Le volume de calcul nécessaire est donc très important. Ce problème est crucial, notamment dans les cas d'une analyse embarquée dans un véhicule. A cet effet, suivant un premier aspect, l'invention propose un procédé d'analyse d'une scène dynamique observée à l'aide d'un bloc de capteur(s) du type précité, caractérisé en ce que la deuxième probabilité d'occupation de la cellule et de l'ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération est déterminée par résolution de l'équation 0 -iv-iP (CAO 0-11(V -1Z) P (011 IZC) EA00-ivy-i P (CAO -11717 -1Z) où C : identifiant de la cellule considérée ; A : identifiant de de la cellule qui contenait, à la (k-1)ème itération ce que contient la cellule considérée à la kème itération ; O : état d'occupation de la cellule considérée, parmi les états « vide » et « occupé » ; 0-1 : état d'occupation de la cellule A à la (k-1)ème itération ; V : vitesse de la cellule C considérée ; V-1 : vitesse de la cellule A à la (k-1)ème itération ; Z : observations des capteurs depuis la première itération jusqu'à la kème itération ; - la vitesse et la position respectives d'un ensemble de particules fictives dans la grille sont déterminées à la kème itération en fonction des vitesses, positions des particules déterminées à la (k-1)ème itération et de la probabilité P(V1V-1) ; - le procédé comprend une étape de détermination des particules localisées en chaque cellule en fonction des positions déterminées et en ce que la résolution de l'équation, pour une cellule, est scindée en la résolution d'une partie statique correspondant à P(0=vide, V=01ZC) et P(0=occupé, V=01ZC) et en la résolution d'une partie dynamique correspondant aux P(0=occ, V= i = 1 à nk, où nk est le nombre de particules déterminées dans la cellule C pour la kème itération ; la partie statique de la cellule C à la kème itération étant déterminée en fonction de la partie statique de la cellule C déterminée à la (k-1)ème itération et de P(010-1) ; la probabilité de P(O=occupé, V= vk ,1ZC) de la partie dynamique de la cellule C étant déterminée à la kème itération en fonction de la probabilité P(O=occupé, V= -1)-1,1ZA) calculée à la (k-1)ème itération pour la partie dynamique de la cellule A et de P(010-1), où la particule p, déterminée dans la cellule C à la kème itération avec une vitesse Vk se trouvait dans la cellule A à la (k-1)ème itération avec une vitesse -1)-1. Le procédé selon l'invention permet de diminuer nettement les ressources de calcul et mémoire nécessaires pour la détermination de l'occupation des cellules et des distributions de probabilités de mouvement lors des itérations. Dans des modes de réalisation, le procédé d'analyse suivant l'invention comporte en outre une ou plusieurs des caractéristiques suivantes : - à l'issue de la kème itération, un ou des couples (p,, (-1) )), où xk est la position de la particule p, à la kème itération, sont dupliqués au sein d'une cellule (C) ou supprimé de sorte que le nombre de couples par cellule soit fonction de la partie dynamique déterminée pour la cellule ; - le procédé comporte une étape de sélection du couple à dupliquer/supprimer, la sélection dudit couple étant réalisée en fonction de la probabilité P(O=occupé, V= -1)1ZC) déterminée à la kème itération, où vk est la composante de vitesse du couple ; - le nombre total de particules dans la grille est constant au fil des itérations ; P(O=vide, V=01ZC) déterminé à la kème itération et noté Pk(0=vide, V=01ZC), est déterminé en fonction du produit d'un premier terme fonction de la première probabilité d'occupation et d'un deuxième terme, ledit deuxième terme étant fonction de : - Pk-1(0=OCC, V=01ZO).(1 -c) Pk-1(0=vide, V=01ZC) .£, et/ou - coeff ( k ). pk-1(0=occupé, V= ik ,1ZA) .(1 - où coeff (-1) est une fonction décroissante de Fkll, et/ou - une probabilité d'apparition pa d'un nouvel objet dans la scène observée. Suivant un deuxième aspect, la présente invention propose un programme d'ordinateur à installer dans un module d'analyse d'une scène dynamique observée à l'aide d'un bloc de capteur(s) comprenant une étape de définition d'une grille découpée en cellules et correspondant à la scène observée, ledit programme comportant des instructions pour mettre en oeuvre les étapes d'un procédé suivant le premier aspect de l'invention lors d'une exécution du programme par des moyens de traitement du module d'analyse. Suivant un troisième aspect, la présente invention propose un module d'analyse d'une scène dynamique observée à l'aide d'un bloc de capteur(s) comprenant une étape de définition d'une grille découpée en cellules et correspondant à la scène observée ; ledit dispositif étant adapté pour, sur la base d'une grille découpée en cellules et correspondant à la scène observée, mettre en oeuvre itérativement le traitement de : - collecter au moins une nouvelle observation du bloc de capteur(s) à une kème itération et déterminer, en fonction de la nouvelle observation collectée, une première probabilité d'occupation de chaque cellule de la grille et modélisant le fonctionnement du bloc de capteur(s) ; - déterminer, pour chaque cellule, à la kème itération, une deuxième probabilité d'occupation de la cellule et un ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération, ledit module d'analyse étant caractérisé en ce qu'il est adapté pour déterminer ladite deuxième probabilité d'occupation de la cellule et l'ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération par résolution de l'équation 13(07IZG) = Li10 0-1VV-1 P (CAO 0-11717 -1Z) où C : identifiant de la cellule considérée ; A : identifiant de de la cellule qui contenait, à la (k-1)ème itération ce que contient la cellule considérée à la kème itération ; O : état d'occupation de la cellule considérée, parmi les états « vide » et « occupé » ; 0-1 : état d'occupation de la cellule A à la (k-1)ème itération ; V : vitesse de la cellule C considérée ; V-1 : vitesse de la cellule A à la (k-1)ème itération ; Z : observations des capteurs depuis la première itération jusqu'à la kème itération ; ledit module étant adapté pour déterminer, à la kème itération, la vitesse et la position respectives d'un ensemble de particules fictives dans la grille en fonction des vitesses, des positions des particules déterminées à la (k-1)ème itération et de la probabilité P(V1V1) ; ledit module étant adapté pour déterminer des particules localisées en chaque cellule en fonction des positions déterminées et pour scinder la résolution de l'équation, pour une EAU-1V-1 P (C'il0 0 -1171/7- Z) cellule, en la résolution d'une partie statique correspondant à P(O=vide, V=01ZC) et P(O=occupé, V=01ZC) et en la résolution d'une partie dynamique correspondant aux P(O=occ, V= i = 1 à nk, où nk est le nombre de particules déterminées dans la cellule C pour la kème itération ; ledit module étant adapté pour déterminer la partie statique de la cellule C à la kème itération en fonction de la partie statique de la cellule C déterminée à la (k-1)ème itération et de P(010-1) ; et pour déterminer la probabilité de P(O=occupé, V= vk ,1ZC) de la partie dynamique de la cellule C à la kème itération en fonction de la probabilité P(O=occupé, V= -1)-1,1ZA) calculée à la (k-1)ème itération pour la partie dynamique de la cellule A et de P(010-1), où la particule pi déterminée dans la cellule C à la kème itération avec une vitesse se trouvait dans la cellule A à la (k-1)ème itération avec une vitesse -1)-1. Ces caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple, et faite en référence aux dessins annexés, sur lesquels : la figure 1 représente une grille d'occupation et l'histogramme de transition d'une cellule de la grille ; la figure 2 représente une vue d'un véhicule et d'une grille couvrant une scène autour du véhicule dans un mode de réalisation de l'invention ; la figure 3 représente un dispositif d'analyse dans un mode de réalisation de l'invention ; la figure 4 représente un ensemble d'étapes mises en oeuvre dans un mode de réalisation de l'invention ; la figure 5 représente une grille d'occupation et la représentation des probabilités de mouvement associées à une cellule de la grille dans un mode de mise en oeuvre de l'invention. Conformément à une convention usuelle du domaine des probabilités, si U,V,W, variables quelconques, P(U) désigne la probabilité d'obtenir U, et P(U I V) désigne la probabilité d'obtenir U, la valeur de V étant connue. Par ailleurs, l'expression de la probabilité conjointe, qui permet l'estimation des fonctions de probabilité de combinaisons de variables, est la suivante : P(UVW) P (W1 U) w P(UVW) Un algorithme BOF estime, au cours d'itérations successives, et à l'aide d'une grille spatiale décomposée en cellules et représentant la scène analysée, l'occupation des cellules de la grille et le mouvement du contenu des cellules, à l'aide de probabilités fonction de variables aléatoires représentant les grandeurs à étudier, qui sont associées à chaque cellule et dont les évolutions au fil des itérations décrivent le déroulement de la scène observée (cf. notamment Pierre Bessière, Emmanuel Mazer, Juan Manuel Ahuactzin-Larios, et Kamel Mekhnacha, « Bayesian Programming », CRC Press, December 2013, ou MK Tay, Kamel Mekhnacha, Cheng Chen, and Manuel Yguel, « An efficient formulation of the bayesian occupation filter for target tracking in dynamic environments », International Journal of Vehicle Autonomous Systems), 6(1):155{171, 2008). Les variables considérées dans le cas présent sont décrites ci-après : - C : identifiant de la cellule considérée ; - A : identifiant de l' « antécédent » : l'antécédent est la cellule qui contenait, à l'itération précédente, ce que contient la cellule considérée à l'itération courante ; 0 : occupation de la cellule considérée ; l'occupation prend un état parmi « vide » et « occupé » ; 0-1 : occupation de l'antécédent à l'itération précédente ; V : vitesse (vecteur de dimension 2) de la cellule considérée ; V-1 : vitesse de l'antécédent à l'itération précédente Z : observations des capteurs depuis la première itération jusqu'à l'itération courante. La probabilité conjointe est exprimée à l'aide de l'équation (1) ci-dessous : P(C A 0 0-1 V V-1 Z) = P(A) P(0-1 V-11A) P(0 V10-1 V-1) P(C 1A V) P(Z O V C), où : les probabilités P(A) représentent la distribution sur tous les antécédents possibles A; elle est, dans le mode de réalisation considéré, choisie uniforme, une cellule étant considérée comme atteignable depuis toutes les autres cellules de la grille avec une probabilité égale ; F,(01 v ", 1 A) est la distribution jointe conditionnelle sur l'occupation et le mouvement des antécédents ; les probabilités correspondant à cette distribution sont calculées à l'itération précédant l'itération considérée ; P(0 V 1 0-1 V-1) est le modèle dynamique de prédiction ; le modèle choisi peut être décomposé comme le produit de deux termes : P(01 0-1)P (V1 V-10), où : P(01 0-1) est la distribution conditionnelle sur l'occupation de la cellule considérée, l'occupation de l'antécédent étant connue ; cette 1-e e distribution est définie par la matrice de transition T, T = e 1-e qui autorise le changement d'état d'occupation d'une cellule avec une probabilité basse pour prendre en compte des erreurs d'évaluation et d'observation ; P (V 1 V-10) est la distribution conditionnelle sur la vitesse de l'itération courante, la vitesse à l'itération précédente étant connue ; usuellement elle suit une loi normale, centrée sur la vitesse à l'itération précédente pour représenter un modèle d'accélération gaussien ; cette distribution peut être adaptée pour correspondre à différents objets observés ; et comme il est considéré qu'une cellule vide ne peut pas bouger, une fonction de Dirac est ajoutée pour interdire un mouvement dans ce cas ; - P(C 1 A V) est la distribution qui indique si la cellule C considérée est accessible depuis l'antécédent connu A = a et avec la vitesse connue V= v ; cette distribution est une fonction de Dirac avec une valeur égale à 1 si et seulement si a+v.dt C C, où dt est le pas de temps séparant l'itération courante de l'itération précédente ; - P(Z 1 O V C) est la distribution sur les valeurs d'observation du capteur, encore appelé « modèle capteur ». Un algorithme BOF estime, au cours d'itérations successive, l'occupation des cellules et la vitesse du contenu des cellules en fonction d'observations courantes fournies par le(s) capteur(s) pour chaque cellule, par la détermination de P(OV1ZC). Par discrétisation des cellules-antécédents et des vitesses, P(OV1ZC) peut s'écrire à l'aide de l'équation (2) suivante : P(OVIZC = ZAD-1L-1 P (CAO 0-1171/ -1Z) ) LAO G-117V- P (CAO 0 -117V-1Z) Sur la figure 2, un véhicule 10 est représenté en vue de dessus. Le véhicule 10, outre les moyens usuels de locomotion propres aux véhicules, comporte un bloc capteur 11 et un dispositif d'analyse 20 dans un mode de mise en oeuvre de l'invention. Comme représenté schématiquement en figure 3, le dispositif d'analyse 20, appelé aussi module HSBOF (pour « Hybrid Sampling BOF ») comporte un microprocesseur 21 et une mémoire 22. La mémoire 22 comporte des instructions logicielles, qui lorsqu'elles sont exécutés par le microprocesseur 21, mettent en oeuvre les étapes indiquées ci-dessous et réalisées par le dispositif d'analyse 20. Dans le cas considéré, le dispositif d'analyse 20 comporte en outre un accélérateur de calcul parallèle adapté pour traiter en parallèle des opérations mises en oeuvre dans le procédé détaillé ci-dessous. Le bloc capteur 11 est adapté d'une part pour estimer les déplacements du véhicule 10, sous la forme de translations et rotations. Ces estimations de déplacement sont fournies au dispositif d'analyse 20. Le bloc capteur 11 comprend en outre un capteur embarqué (ou plusieurs capteurs embarqués) permettant d'effectuer des mesures et, en fonction de ces mesures, de détecter la présence de zones occupées par exemple assimilées selon certains critères à des obstacles dans l'environnement du véhicule 10 et leurs positions par rapport à un référentiel (X1, X2) lié au véhicule (on notera que dans d'autres modes de réalisation, une carte instantanée de profondeur est déterminée en fonction de ces mesures). Un tel capteur comprend par exemple un télémètre laser dotée d'une source émettant un faisceau laser et qui effectue un balayage d'une surface cp dans le plan (X1, X2), et/ou une caméra stéréo, etc. et détecte la présence d'obstacle(s) dans la scène balayée, et leur position en fonction de l'analyse de faisceaux réfléchis reçus par le télémètre. Ces données d'observation de la scène comportant la position d'obstacles détectés sont fournies par le capteur 11 au dispositif d'analyse 20. Dans des modes de réalisation, des données d'observations de la scène indiquant la présence d'obstacle(s) dans la scène balayée, et leur position, sont fournies au dispositif d'analyse 20 par un ou des capteurs extérieurs 13 au véhicule 10, par exemple un capteur fixé sur le bord de la route sur laquelle le véhicule 10 circule. On considère par ailleurs une grille spatiale 14 correspondant à une portion de l'environnement du véhicule 10, par exemple une vue de dessus d'une zone de largeur et longueur fixées dans le plan (X1, X2), par exemple de 20 mètres chacune, qui s'étend par exemple depuis le véhicule vers l'avant en direction de l'axe longitudinal X1 du véhicule. Cette grille 14 est découpée en cellules C, par exemple 100 000 cellules, de manière à constituer une grille spatiale représentant une vue de dessus de l'environnement du véhicule. Le dispositif d'analyse 20 est adapté pour déduire, à chaque instant tk=to+k.dt, de ces données d'observation de l'environnement du véhicule reçues du ou des capteurs 11, la distribution sur les valeurs d'observations reçues, encore appelé « modèle probabiliste capteur » : P(Z O V C). Ainsi soit zk les observations de l'environnement du bloc capteur 11 à considérer pour l'instant tk. La probabilité P(Z O V C) déterminée à l'instant tk pour chaque cellule C de la grille 14 par le module d'analyse 20 en fonction des observations z1 à zk reçues aux instants d'itération t1 à tk modélise le comportement du capteur 11 et prend en compte le bruitage de la mesure par le bloc capteur 11, les facteurs d'imprécision tels que la température, l'éclairage, les couleurs des objets, les masquages d'objets, la capacité du capteur à détecter ou non les objets et à détecter des objets quand il n'y en a pas (« fausses détections »). On notera que l'observation Z ne prend pas dans le cas considéré la forme des coordonnées des objets détectés dans le référentiel lié au véhicule et d'axes X1, X2, mais prend la valeur 1 dans une cellule C ou un obstacle a été détecté et 0 ailleurs. Le modèle probabiliste capteur est selon les cas estimé à partir d'une modélisation des caractéristiques physiques du capteur et des objets à détecter, ou « appris » expérimentalement. Des exemples de modèles de capteurs sont par exemple décrits dans FR 0552736 ou encore à l'adresse http://ieeexplore.ieee.orq/stamp/stamp.isp?tp=&arnumber= 30720&tab=1 . Dans le cas présent, en considérant que le module d'analyse 20 collecte à un instant tk des données d'observation actualisées fournies par m capteurs, chacun fournissant des observations ,j = 1 à m, le vecteur observation zk = (z1, z2, zm). Typiquement, un laser donne une observation correspondant à une liste de valeurs de distance d'impact pour différentes directions. Le module d'analyse 20 obtient à chaque instant tk un vecteur d'observations actualisé zk. A partir de ces vecteurs instantanés, le dispositif d'analyse 20 déduit de zk et du "modèle de capteur" (en l'occurence ici multicapteur) la probabilité P( Z = zk O V C), ce qui donne une grille instantanée (indiquant pour chaque C, si la case était occupée, et avec telle vitesse, quelle serait la probabilité d'avoir ces observations) qui est fournie ensuite comme donnée d'entrée de l'étape 103 détaillée plus bas et dans laquelle les probabilités P( O V z' ...zk C) sont calculées récursivement en utilisant les évaluations P( O VI z' ...zk-' C) calculées à tk-i. L'invention repose sur les principes suivants : décomposition de chaque cellule en une partie statique correspondant à la partie non occupée de la cellule et à la partie de la cellule occupée par des objets de vitesse nulle, et en une partie dynamique correspondant à la partie de la cellule occupée par des objets en mouvement ; représenter la distribution de mouvement d'une cellule de la grille d'occupation par un ensemble de particules représentant la partie dynamique de la cellule, à la place d'une grille régulière de distribution de mouvement, ce qui permet une représentation plus compacte. Ainsi la détermination de P(OV1ZC) revient à calculer la probabilité associée à chaque particule P(O=occ, V=vk I ZC) et les probabilités associées au domaine statique P(O=occ, V=0 I ZC) et P(O=vide, V=0 ZC). En utilisant la décomposition (1) et en considérant le dénominateur de l'équation (2) comme un facteur d'échelle, cette équation (2) peut être réécrite de la manière suivante : P(OVIZC) « P(ZIOC) P(A)P(0-'V1A)P(00-1)P(VOV-1)P(CIAV) où le sigle « « » signifie « est proportionnel à ». (3) On note a(o,v) = P(A)P(0-117-1k1)P(o10-1)P(+17-1)P(CIAv) (4) A0-117-1 où o a pour valeur l'état d'occupation « occupé » (nommé « occ ») ou « vide », et v est une vitesse. Procédé 100 Un ensemble 100 d'étapes 101-106 mises en oeuvre dans une kème itération à l'instant tk=to+k.dt, par le dispositif d'analyse 20 est décrit ci-après en référence à la figure 4 dans un mode de réalisation de l'invention (kk1). Il est considéré qu'à l'issue de l'itération précédente, i.e. (k-1)ème itération, pour chaque cellule C de la grille, il a été obtenu les probabilités associées au domaine statique : P(O=occ, V=0 ZC), nommé Pk-1(0=occ, V=0 ZC), P(O=vide, V=0 I ZC), nommé Pk-1(0=vide, V=0 ZC), et qu'un ensemble de particules a été attribué à la cellule C, tel que chaque particule dans cet ensemble a une vitesse respective -1)-1 et une position respective .x-1 et qu'il a été obtenu les probabilités associées au domaine dynamique P(0=occ, V= k-1 I I ZC), nommé pk1(0=OCC, V= 1)-1 ZC). Dans le mode de réalisation considéré, la position .x-1 permet de positionner exactement les particules au sein de la cellule. L'introduction de positions précises des particules au sein d'une cellule permet de résoudre des problèmes de crénelage spatio- temporel. Dans une étape 101 de détermination de la propagation des particules, pour chaque particule, le modèle de mouvement est appliqué, conformément à la distribution P(V I V-1) : k' k-1 V, =V, +G k' k-1 k' = di. V, où o est un bruit gaussien centré en 0 de covariance E. Dans une étape 102, le déplacement de la grille correspondant à un déplacement du véhicule 10 le cas échéant est pris en compte, en fonction de la rotation Rk et de la translation Tk effectuées entre l'itération (k-1) et l'itération courante k et déterminées ar le bloc capteur 11 :: k k' V, = Rk . xk Rk xk' + Tk Dans une étape 103, les observations zk de la scène à considérer pour l'instant tk sont collectées (sous-étape 103_0) et les probabilités P(Z I O , V, C) issues de ces observations zk et modélisant le capteur, comme décrit plus haut, sont déduites de ces observations. Puis les probabilités d'occupation et de mouvement prenant en compte ces observations collectées sont mises à jour (sous-étapes 103_1 à 103_3 détaillées ci-après). La partie statique de chaque cellule C est mise à jour, dans la sous-étape 103 1, en utilisant la part statique de la cellule considérée C déterminée à l'issue de la (k-1)ème itération, i.e. Pk-1(occ, 0 Z, C) et Pk-1(vide,01Z,C). En effet, P(0-1 =occ,17-1 =01Z,A=C)= Pk-1(occ, 0 IZ,C), et P(0-1 = vide,V-1 = 0 Z, A= C) = Pk-1 (vide, 0 Z, C). Il découle que les coefficients a(occ,0) et a(vide,0) d'après la définition donnée à l'équation (4) sont obtenus pour la cellule C considérée par le calcul suivant : a(occ,0) = Pk-1(occ,0 I Z, C) .(1 -e) + Pk-1 (vide, 0 IZ, C) .c a(vide,0) = Pk-1 (occ, 01Z , C) . c + Pk-1(vide, 0 1Z, C) .(1-c), OÙ c = P(O = occ 0-1 = vide) = P(0 =vide0-1 = occ). La partie dynamique des cellules est à son tour mise à jour, en considérant les particules. Ainsi, pour chaque cellule C de la grille, dans une étape 103_2 : - on détermine quelles sont les particules, parmi les ensembles de particules associées aux cellules, qui se retrouvent dans la cellule C après la mise en oeuvre des étapes 101 et 102. Considérons que nck particules se retrouvent dans la cellule C considérée. Chacune de ces particules a une vitesse respective vk et une position courante xk au sein de la cellule C et était dans une cellule, nommée a, à l'itération k-1, telle que déterminées à l'issue de l'étape 102 ; et on définit le coefficient de pondération associé 142-1 à cette particule : 142-1 = P(0-1 = occ, V-1 = vk 1 1 A =a,, z1,..., zk-1) qui est égal à Pk-1(0=occ, V= -1)-1 I Z, a,), déterminé à l'itération k-1 (sous-étape 103_3) relativement à la cellule a, : - le coefficient a(occ, ) d'après la définition donnée à l'équation (4) est alors déterminé pour la cellule C considérée comme égal à : a(occ, ) = w,k- 1 . (1 - c) ; une cellule vide est supposée de vitesse nulle : a(vide, vik ) =o. Dans une sous-étape 103_3 de mise à jour des probabilités, pour chaque cellule C, pour terminer la résolution de l'équation (3), on calcule [3(occ,0), [3(vide,0) et [3(occ, vk ) de la façon suivante, à l'aide des probabilités d'observation actualisées P(Z I O , V, C) issues de la sous-étape 103_0 et des coefficients a(occ,0), a(vide,0) et a(occ, )) ) calculées lors des sous-étapes 103_1 et 103_2 : [3(occ,0) = P(Z I occ , 0, C).a(occ,0) [3(vide,0) = P(Z I vide , 0, C).a(vide,0) [3(occ, vik ) = P(Z I occ , vik , C).a(occ, vik ). Le dénominateur de l'équation (2), nommé facteur d'échelle « d », est calculé : d = [3(occ,0) + [3(vide,0) + E, 3(occ, vk ), où E, [3(occ, ) est la somme des coefficients [3(occ, ) calculés pour les particules se trouvant dans la cellule C. Finalement, on obtient pour la cellule C les probabilités suivantes qui représentent I"occupation de filtrage bayésien et les probabilités de mouvement associés : c Pk(occ, 0 I Z,C) - fl(odc,0) pour la vitesse nulle, Pk(occ, )6(occ,v1:) I Z,C) - pour chaque vitesse vik , les vitesses vk étant les d vitesses des particules déterminées dans la cellule C, Pk(vide, 0 I Z,C) - 13(vide,0)pour la vitesse nulle. d Une étape 104 (optionnelle) de réorganisation des particules est ensuite mise en oeuvre pour adapter le nombre de particules par cellule et par exemple dupliquer les particules à haute probabilité, ceci afin de permettre une meilleure efficacité de l'ensemble d'étapes 100, en mettant en oeuvre le processus Proc_reorg indiqué ci-dessous. On notera que la précision de la représentation de la distribution de vitesse dans une cellule croît avec le nombre de particules dans la cellule. Dans le mode de réalisation considéré, le nombre total de particules, nommé nbPart, est maintenu constant au fil des itérations (la charge de calcul est ainsi la même pour chaque itération, la mémoire à allouer est constante et connue). Ce nombre total peut évoluer dans d'autres modes de réalisation (par exemple, le nombre total pourrait être augmenté pour une scène très dynamique et être diminué pour une scène plus statique). On considère pour chaque cellule C l'ensemble pondéré des nck particules défini par les nck triplets { xk , vk , wl associés chacun à une particule, où wk = Pk(O=occ, V= vk Z, C). Selon le processus Proc reorg, on réitère un nombre égal à nbPart l'ensemble d'étapes successives suivantes i, ii et iii, après avoir initialisé à un ensemble vide, par cellule de la grille, un nouvel ensemble de particules attribué à la cellule : i/ tirer un identifiant de cellule selon une loi d'occurrence proportionnelle au poids de la cellule (le poids d'une cellule C est égal à la somme des nck poids wk des particules dans la cellule et représente donc Pk(0=occ, V#01Z, C)) ; ii/ sélectionner aléatoirement une particule parmi les nck particules de cette cellule, selon une loi d'occurrence proportionnelle au poids wk de la particule ; iii/ ajouter la particule sélectionnée en ii/ au nouvel ensemble de particules attribué à la cellule dont l'identifiant a été tiré en i/. Une fois ces étapes réitérées, le poids des particules est normalisé par cellule, i.e. le poids de chaque particule dans le nouvel ensemble de particules attribué à la cellule est divisé par la somme des poids des particules dans ce nouvel ensemble ; ce sont ces poids normalisés qui sont considérés comme les poids wk pour l'itération suivante. Dans une étape 106, la valeur de k est augmentée de 1. Lorsqu'un nouvel objet apparaît sur la grille, il est nécessaire d'initialiser une ou des particules correspondante(s) avec une distribution de vitesse. Pour cela, dans un mode de réalisation, une probabilité d'apparition constante avec une distribution uniforme entre les états occupé/vide et entre statique/dynamique est ajoutée. Par exemple, dans l'étape 103 1, on calcule a(occ,0), respectivement a(vide,0),en ajoutant pa/4, respectivement pa/2, soit a(occ,0) = Pk-1(occ,0 I Z, C) .(1-c) + Pk-1 (vide, 01Z ,C) .0 + pa/4 a(vide,0) = Pk-1 (OCC, 01Z , C) . + Pk-1 (vide, 0 Z, C) .(1-c) + pa/2 ; et dans l'étape 103 2, on considère pour la cellule un terme dynamique supplémentaire a(occ, V = « inconnue »)= pa/4. Ainsi dans l'étape 104 de réorganisation des particules, quand la vitesse « inconnue » est sélectionnée (étape ii/), la vitesse est déterminée à partir d'une distribution uniforme, en sélectionnant uniformément une vitesse sur [-V', +V'x]x[-V'x, +V'x] et en l'associant à la cellule. Dans un mode de réalisation qui peut se combiner ou non avec le précédent, pour empêcher des particules de suivre un objet statique, chaque particule contribue à la part statique représentée par a(occ,0) calculé à l'étape 103 1, par un coefficient de valeur décroissant avec la norme de vitesse : a(occ,0) = Pk-1(occ, 01Z, C) .(1-c) + 13k-1 (vide, 01Z ,C) .c + à e 2cr>. 142-1 *(1- e) ; -0v=` et à l'étape 103_2 : a(occ, ) = (1- e 2°1 ). "142-1 . (1 - £), où as est un coefficent constant représentant la contribution de vitesse statique (pour 111)k > 3 Dans une étape d'initialisation préalable à la première itération (k=1) de l'ensemble d'étape 100, il est procédé, dans un mode de réalisation, comme pour les antécédents de type "inconnu" : tirage uniforme de la vitesse sur [ - Vmax , Vmax ] x [ - Vmax , Vmax ]. La position est prise en tirant uniformément sur la grille (tirage d'une position flottante sur [- gridWidth/2, gridWidth/2] x [-gridHeight/2, gridHeight/2], où gridWidth est la largeur de la grille et gridHeight est la hauteur de la grille), l'occupation est fixée à 0,5 (inconnue). Les w° sont égaux à 1 / nc° , ce nombre ng de particules dans la cellule C étant un résultat du tirage aléatoire précédent. La probabilité statique initiale est un paramètre sans grand impact, car la réévaluation a lieu très rapidement. Ainsi selon l'invention, la représentation classique BOF comprenant une grille d'occupation et associant à chaque cellule de la grille un histogramme de vitesse est remplacée par une grille d'occupation avec, pour chaque cellule, une probabilité d'occupation, un coefficient statique correspondant à une probabilité de vitesse nulle, un ensemble de particules dans la cellule, ayant chacune une vitesse respective v, et un ensemble de probabilités P(V = y) . L'invention permet de construire à chaque itération, via la grille, une carte probabiliste discrète de la scène, en distinguant les composantes statiques et dynamiques. Le mouvement des composantes dynamiques est estimé à l'aide des particules. La qualité des estimations d'occupation et de vitesse est améliorée, par une meilleure précision des vitesses calculées et une diminution des effets de crénelage spatio-temporel. A titre d'illustration, la grille d'occupation 14 obtenue dans un mode de réalisation de l'invention est représentée. La cellule 15 de la grille d'occupation 14 comporte par exemple 3 particules pu, à des positions respectives dans la cellule, ayant des vitesses respectives vu, j = 1 à 3. Les observations et estimations de déplacement du véhicule sont effectuées par les capteurs et fournies au dispositif d'analyse 20, puis traitées en temps réel par le dispositif d'analyse 20. Par exemple le pas d'itération dt est égal à 0,06 seconde., les observations de rafraîchissement de la grille et estimations de déplacement étant fournies au dispositif d'analyse 20 selon des pas de temps de durée inférieurs ou égaux à dt. Les valeurs d'occupations et de vitesses associées aux particules déterminés pour chaque cellule à chaque itération peuvent, dans des modes de réalisation de l'invention, être exploitées de différentes manières dans une étape 105 : action/décision courante en fonction des valeurs d'occupations et de champs de vitesse déterminés pour chaque cellule (par prédiction d'état(s) futur(s) de l'environnement (par exemple évaluation de risques de collisions entre le véhicule 10 et des objets correspondant à des cellules dont la probabilité d'occupation est supérieure à un seuil d'occupation donné et la distance observée est au-dessous d'un seuil de distance donné), puis par exemple action/décision courante en fonction d'un ou des états futurs déterminés (freinage automatique du véhicule 10 par exemple) ; on notera que les distributions de vitesses inférées par l'invention sont directement utilisables pour prédire l'évolution de la scène, et donc déduire l'intersection entre le véhicule et des obstacles, y compris les objets dynamiques. Les étapes du procédé 100 sont aisément parallélisables. La détermination de la grille d'occupation et des champs de vitesse est réalisable pour un environnement analysé comportant toute sorte d'objets préalablement connus ou non à l'aide de capteurs même hétérogènes, grâce à l'utilisation de modèles de capteur probabilistes. La présente invention permet de diminuer nettement les ressources en calcul et mémoire nécessaires pour la détermination de l'occupation des cellules et des distributions de probabilités de mouvement. Elle permet d'atteindre gain d'un facteur 100. Elle permet de compenser le mouvement propre des capteurs (calculer les opérations de rotation/translation des grilles et des particules sur les champs de vitesses lorsqu'ils étaient considérés sous forme d'histogrammes BOF était en pratique impossible du fait d'un crénelage trop important et d'une interpolation 4D trop coûteuse en calculs).35 REVENDICATIONS1.- Procédé d'analyse d'une scène dynamique observée à l'aide d'un bloc (11) de capteur(s) comprenant une étape de définition d'une grille (14) découpée en cellules (C) et correspondant à la scène observée ; et la mise en oeuvre itérative par ordinateur des étapes : collecter au moins une nouvelle observation du bloc de capteur(s) à la kème itération et déterminer, en fonction de la nouvelle observation collectée, une première probabilité d'occupation de chaque cellule de la grille modélisant le fonctionnement du bloc de capteur(s) ; déterminer, pour chaque cellule, à la kème itération, une deuxième probabilité d'occupation de la cellule et un ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération, ledit procédé étant caractérisé en ce que la deuxième probabilité d'occupation de la cellule et de l'ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération est déterminée par résolution de l'équation p(o 1v-1 P (CAO 0-1171/7-1-Z) - v izc) = A° la vitesse et la position respectives d'un ensemble de particules fictives dans la grille 30 sont déterminées à la kème itération en fonction des vitesses, positions des particules déterminées à la (k-1)ème itération et de la probabilité P(V1V-1) ; le procédé comprend une étape de détermination des particules localisées en chaque cellule en fonction des positions déterminées et en ce que la résolution de l'équation, pour une cellule, est scindée en la résolution d'une partie statique correspondant à 35 P(0=vide, V=01ZC) et P(0=occupé, V=01ZC) et en la résolution d'une partie AO 0 - ivv-i P (CAO 0 -1171/7-1Z) où C : identifiant de la cellule considérée ; A : identifiant de de la cellule qui contenait, à la (k-1)ème itération ce que contient la cellule considérée à la kème itération ; O : état d'occupation de la cellule considérée, parmi les états « vide » et « occupé » ; 0-1 : état d'occupation de la cellule A à la (k-1)ème itération ; V : vitesse de la cellule C considérée ; V-1 : vitesse de la cellule A à la (k-1)ème itération ; Z : observations des capteurs depuis la première itération jusqu'à la kème itération ;dynamique correspondant aux P(O=occ, V= i = 1 à nk, où nk est le nombre de particules déterminées dans la cellule C pour la kème itération ; la partie statique de la cellule C à la kème itération étant déterminée en fonction de la partie statique de la cellule C déterminée à la (k-1)ème itération et de P(010-1) ; la probabilité de P(O=occupé, V= vk ,1ZC) de la partie dynamique de la cellule C étant déterminée à la kème itération en fonction de la probabilité P(O=occupé, V= -v-1,1ZA) calculée à la (k-1)ème itération pour la partie dynamique de la cellule A et de P(010-1), où la particule p, déterminée dans la cellule C à la kème itération avec une vitesse se trouvait dans la cellule A à la (k-1)ème itération avec une vitesse 2.- Procédé selon la revendication 1, selon lequel à l'issue de la kème itération, un ou des couples (p,, )), où xk est la position de la particule p, à la kème itération, sont dupliqués au sein d'une cellule (C) ou supprimé de sorte que le nombre de couples par cellule soit fonction de la partie dynamique déterminée pour la cellule. 3.- Procédé selon la revendication 2, comportant une étape de sélection du couple à dupliquer/supprimer, la sélection dudit couple étant réalisée en fonction de la probabilité P(O=occupé, V= -vIZC) déterminée à la kème itération, où vk est la composante de vitesse du couple. 4.- Procédé selon la revendication 2 ou 3, selon lequel le nombre total de particules dans la grille (14) est constant au fil des itérations. 5.- Procédé selon l'une des revendications précédentes, selon lequel P(O=vide, V=01ZC) déterminé à la kème itération et noté Pk(0=vide, V=01ZC), est déterminé en fonction du produit d'un premier terme fonction de la première probabilité d'occupation et d'un deuxième terme, ledit deuxième terme étant fonction de : - Pk-1(0=occ, V=01ZC).(1-e) + Pk-1(0=vide, V=01ZC) .c, et/ou - coeff ( ). Pk-1(0=occupé, V= 1 ) ,1ZA) .(1 - où coeff ( v est une fonction décroissante de 11-v11 , et/ou - une probabilité d'apparition pa d'un nouvel objet dans la scène observée. 6.- Programme d'ordinateur à installer dans un module d'analyse (20) d'une scène dynamique observée à l'aide d'un bloc (11) de capteur(s) comprenant une étape de définition d'une grille (14) découpée en cellules (C) et correspondant à la scène observée, ledit programme comportant des instructions pour mettre en oeuvre les étapes d'un procédé selon l'une des revendications 1 à 5 lors d'une exécution du programme par des moyens de traitement du module d'analyse. 7.- Module d'analyse (20) d'une scène dynamique observée à l'aide d'un bloc (11) de capteur(s) comprenant une étape de définition d'une grille (14) découpée en cellules (C) et correspondant à la scène observée ; ledit dispositif étant adapté pour, sur la base d'une grille (14) découpée en cellules (C) et correspondant à la scène observée, mettre en oeuvre itérativement le traitement de : - collecter au moins une nouvelle observation du bloc de capteur(s) à une kème itération et déterminer, en fonction de la nouvelle observation collectée, une première probabilité d'occupation de chaque cellule de la grille et modélisant le fonctionnement du bloc de capteur(s) ; - déterminer, pour chaque cellule, à la kème itération, une deuxième probabilité d'occupation de la cellule et un ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération, ledit module d'analyse étant caractérisé en ce qu'il est adapté pour déterminer ladite deuxième probabilité d'occupation de la cellule et l'ensemble de probabilité(s) de mouvement du contenu de la cellule en fonction de la première probabilité d'occupation de la cellule déterminée à la kème itération par résolution de l'équation P-(Q EA0-1v- P (CAC 0-11717 Z) I/ IZ C) = v, LAO 0-117V-1 P (CA00-11717-12) où C : identifiant de la cellule considérée ; A : identifiant de de la cellule qui contenait, à la (k-1)ème itération ce que contient la cellule considérée à la kème itération ; O : état d'occupation de la cellule considérée, parmi les états « vide » et « occupé » ; 0-1 : état d'occupation de la cellule A à la (k-1)ème itération ; V : vitesse de la cellule C considérée ; V-1 : vitesse de la cellule A à la (k-1)ème itération ; Z : observations des capteurs depuis la première itération jusqu'à la kème itération ;ledit module étant adapté pour déterminer, à la kème itération, la vitesse et la position respectives d'un ensemble de particules fictives dans la grille en fonction des vitesses, des positions des particules déterminées à la (k-1)eme itération et de la probabilité P(VIV- 1 )' ledit module étant adapté pour déterminer des particules localisées en chaque cellule en fonction des positions déterminées et pour scinder la résolution de l'équation, pour une cellule, en la résolution d'une partie statique correspondant à P(0=vide, V=ORC) et P(O=occupé, V=OIZC) et en la résolution d'une partie dynamique correspondant aux P(0=occ, V= v, ,IZC), i = 1 à nk, où nk est le nombre de particules déterminées dans la cellule C pour la kème itération ; ledit module étant adapté pour déterminer la partie statique de la cellule C à la kème itération en fonction de la partie statique de la cellule C déterminée à la (k-1)ème itération et de P(0104) ; et pour déterminer la probabilité de P(O=occupé, V= VIk ,IZC) de la partie dynamique de la cellule C à la kème itération en fonction de la probabilité P(O=occupé, V= ,IZA) calculée à la (k-Dème itération pour la partie dynamique de la cellule A et de P(010-1), où la particule pi déterminée dans la cellule C à la kème itération avec une vitesse viic se trouvait dans la cellule A à la (k-1)ème itération avec une vitesse . 8.- Module d'analyse selon la revendication 7, adapté pour, à l'issue de la kème itération, dupliquer ou supprimer un ou des couples (pi, ( xik )), où xi: est la position de la particule pi à la kème itération, au sein d'une cellule (C) de sorte que le nombre de couples par cellule soit fonction de la partie dynamique déterminée pour la cellule. 9.- Module d'analyse selon la revendication 8, adapté pour sélectionner ledit couple à dupliquer/supprimer, en fonction de la probabilité P(0=occupé, V= vik IZC) déterminée à la kème itération, où 1,/,' est la composante de vitesse du couple. 10.- Module d'analyse selon la revendication 8 ou 9, selon lequel le nombre total de particules dans la grille (14) est constant au fil des itérations. 11.- Module d'analyse selon l'une des revendications 7 à 10, adapté pour déterminer P(O=vide, V=01ZO) à la kème itération, noté Pk(0=vide, V=01ZO), en fonction du produit d'un premier terme fonction de la première probabilité d'occupation et d'un deuxième terme, ledit deuxième terme étant fonction de : - Pk-1(0=occ,V=01ZO).(1-e) + Pk-1(0=vide, V=01ZO) .c ; et/ou - coeff (vk ). Pk-1(0=occupé, V= 1 ) ,1ZA) .(1 - où coeff (vk) est une fonction décroissante de Ovk , et/ou - une probabilité d'apparition pa d'un nouvel objet dans la scène observée. 15