FR 2890774 A1 20070316 FR 0552736 A 20050909 L'invention a pour domaine celui de l'assistance à la conduite d'un véhicule. Plus particulièrement, l'invention est relative à une fonction d'assistance permettant d'anticiper les collisions entre le véhicule et des objets présents dans l'environnement du véhicule considéré. Des fonctions d'assistance à la conduite sont disponibles sur les véhicules actuels pour pallier aux éventuelles erreurs de conduite. Par exemple, la fonction ABS (Anti-Braking System) permet de limiter le glissement des roues et donc de réduire la distance de freinage. Le but ultime de ces fonctions d'assistance est de réduire l'insécurité routière. A ce titre, on a cherché à mettre au point des fonctions permettant de maintenir une distance de sécurité entre deux véhicules. Ainsi, le système ACC ( Adaptative Cruise Control ou contrôle adaptatif de la vitesse de croisière ) est disponible sur certains véhicules. Il permet de détecter le véhicule précédent le véhicule équipé d'un laser ou d'une caméra, d'estimer la vitesse et la distance du véhicule précédent, puis de contrôler la vitesse du véhicule équipé pour maintenir la distance de sécurité. Mais, le système ACC ne peut être utilisé que sur des autoroutes lorsque les situations de conduite sont simples. Il y a donc un besoin pour une fonction d'assistance, permettant d'éviter les collisions, qui fonctionne dans un environnement complexe, en particulier dans un environnement urbain. En effet, dans ce dernier cas, le véhicule risque de rentrer en collision avec des objets de nature variés, parfois vulnérables et qui ont un comportement imprévisible. En particulier, un objet peut être temporairement occulté par un autre objet de la scène. C'est par exemple le cas d'un piéton masqué derrière un véhicule en stationnement. Pour répondre à ce problème, certains travaux s'orientent vers un procédé de poursuite multipistes permettant de maintenir une liste de pistes qui représentent chacune un objet supposé être présent dans l'environnement. A chaque nouvelle série d'observations, il s'agit de réaliser une étape d'association des nouvelles mesures avec les pistes existantes; une étape de modification de la liste de pistes en créant, en confirmant ou en supprimant une piste; une étape de prédiction de l'état instantané de chaque piste au moyen d'un filtre du type filtre de Kalman; et une étape de fenêtrage permettant d'associer à une piste les observations à venir autour de l'état estimé de cette piste. Mais cette méthode multipistes présente des problèmes de fiabilité, en particulier lors de l'étape d'association et de l'étape de modification de la liste des pistes. De plus, la représentation de l'environnement au moyen d'une liste de pistes est la cause d'une perte importante de l'information disponible. Par exemple, le procédé multipiste ne permet pas de retenir l'information sur les zones de l'environnement qui sont masquées par un objet. L'invention a pour but de fournir une fonction d'assistance permettant d'anticiper les collisions qui puisse être mise en ceuvre dans un environnement complexe. L'invention a également pour but de répondre aux inconvénients précités et en particulier de tenir compte des occultations d'objets. Pour être effectivement embarqué à bord d'un véhicule, l'invention doit permettre une prise en compte aussi précise que possible des objets de la scène, en limitant au maximum les erreurs sur la présence d'un objet en un lieu donné de la scène. De plus l'invention doit pouvoir être mise en ceuvre quasiment en temps réel de manière à prendre en compte des modifications rapides de la scène. L'invention sera mieux comprise et d'autres buts, détails, caractéristiques et avantages de celle-ci apparaîtront plus clairement au cours de la description donnée uniquement à titre illustratif et non limitatif en référence au dessin annexé. La Figure 1 représente de manière schématique un véhicule embarquant le dispositif d'aide à la conduite selon l'invention. A Description générale d'un mode de réalisation Le procédé d'assistance à la conduite selon l'invention est mis en ceuvre sur un véhicule 1 représenté schématiquement sur la figure 1. De manière connue, le véhicule 1 comporte quatre roues 2a-d. Les roues avant 2a et 2b sont par exemple entraînées par un moteur thermique. Chacune des roues 2a-d est équipée de moyens de freinage 3a-d. Par exemple, les roues avant 2a et 2b sont équipées de freins hydrauliques, alors que les roues arrière 2c et 2d sont équipées de freins électromécaniques 3c et 3d. Le véhicule 1 comporte une pluralité de capteurs. Des capteurs de vitesse 4a-d sont disposés sur chacune des roues 2a-2d. Chaque capteur de vitesse permet de mesurer la vitesse de rotation de la roue qu'il équipe. Le véhicule 1 est également équipé d'un capteur d'accélération 5. Les différents capteurs délivrent des signaux qui sont acheminés via un réseau de communication (par exemple du type CANBus) jusqu'à une unité centrale de calcul 6. L'unité de calcul 6 comporte un processeur apte à exécuter les instructions d'un programme informatique. L'unité de calcul 6 comporte également des moyens de mémorisation permettent de stocker de manière permanente lesdites instructions et des paramètres prédéfinis ou de manière éphémère les résultats des calculs effectués par ledit processeur. L'unité de calcul comporte également une horloge permettant de cadencer et de synchroniser son fonctionnement. De plus l'unité de calcul 6 comporte des interfaces entrée/sortie permettant la réception en entrée des signaux provenant des différents capteurs et l'émission, en sortie, de signaux de contrôle en direction d'actionneurs. L'unité de calcul 6 calcule à partir des différents signaux émis par les capteurs de vitesse 4a-d, la vitesse longitudinale du véhicule 1. L'unité de calcul 6 peut comporter des fonctions de gestion du freinage, par exemple du type ABS ( Anti-lock braking system pour "anti-blocage des roues"). En fonction de la vitesse, l'unité de calcul 6 détermine le glissement de chacune des roues. Lorsque le glissement d'une des roues est supérieur à un seuil prédéterminé, l'unité de calcul 6 émet un signal d'actionnement des freins de manière à réduire le glissement sur la roue considérée. Ce freinage supplémentaire n'est pas du tout géré par le conducteur du véhicule pour qui cette opération est totalement transparente. Selon l'invention, le véhicule 1 est équipé d'au moins un capteur permettant de détecter la présence d'objet. Différents types de capteurs peuvent être utilisés. Par exemple, sur la figure 1, le capteur 10 représenté est un télémètre laser. Il permet de détecter des objets présents dans un champ de vision ou zone de détection correspondant, dans un plan horizontal, à un cône d'angle d'ouverture totale ôo d'environ 90 sur une profondeur po d'environ 20 m en avant du véhicule. En variante de ce mode de réalisation, le capteur embarqué sur le véhicule pourrait être un détecteur du type laser ou sonar, ou encore une camera CCD. Par caméra CCD, il faut comprendre non seulement le détecteur, mais également l'électronique de pré-traitement de l'image qui produit un signal comportant une information exploitable. Par exemple, l'électronique de pré-traitement est apte à générer un cadre autour d'un objet en mouvement en regroupant l'ensemble des points de l'image animés de la même vitesse. Sur la figure 1, le capteur a été représenté en position avant sur le véhicule 1 pour que la zone de détection soit en avant de celui-ci. En variante, par exemple dans le cas d'un bus pour lequel les angles morts sont importants, le capteur peut être disposé de telle sorte qu'il permette une assistance lors d'une marche arrière ou de la fermeture des portes du bus, la zone de détection étant alors en arrière du véhicule ou bien sur le côté de celui-ci. Dans encore une autre variante de réalisation, les capteurs de détection de présence d'objets peuvent être disposés de manière fixe dans le paysage. Sur la figure 1, nous avons représenté un tel capteur 10' qui peut être fixé sur un élément de l'infrastructure routière tel qu'un abri de bus, une glissière de sécurité, un panneau de signalisation ou l'équivalent. Le capteur 10' est équipé de moyens d'émission, qui, soit directement soit au travers de moyens relais, permettent de transmettre un signal de détection au véhicule 1. Ce dernier, équipé de moyens de réception schématisés par l'antenne 8 est apte à replacer l'objet détecté dans son propre référentiel. Pour déterminer sa position et sa vitesse relativement au capteur 10', le véhicule 1 comporte des moyens de positionnement du type GPS avec éventuellement des balises de positionnement locales. Le télémètre laser comporte une source 10a émettant un faisceau directionnel. En une période de temps repérée par un entier k et ayant une période d'échantillonnage de l'ordre de 10 ms, le faisceau balaie le cône d'ouverture 00. Une partie de l'énergie du faisceau émis est réfléchie par un objet en direction du véhicule 1. Un récepteur 10b détecte le faisceau réfléchi. L'électronique d'acquisition de données 10c du capteur 10 détermine, la distance p et la position angulaire c de l'objet détecté. Un télémètre laser ne permet pas de fournir une mesure de vitesse. C'est pourquoi l'électronique d'acquisition peut éventuellement effectuer une étape de prétraitement permettant d'estimer une valeur de la vitesse de l'objet détecté, sur la base par exemple d'un algorithme de plus proche voisin entre deux mesures successives de la profondeur et de la position angulaire de l'objet détecté. Finalement, le capteur 10 émet un signal correspondant à une mesure en position et, éventuellement, en vitesse de l'objet détecté à l'instant k. Les mesures sont faites relativement au référentiel du véhicule 1. Dans le cas d'un capteur statique disposé dans le paysage, tel que le capteur 10', les données de positionnement d'un objet sont transmises au véhicule 1 qui, de manière simple, peut déterminer la position de l'objet détecté par changement de référentiel, compte-tenu de sa position, de sa vitesse et de son accélération par rapport au capteur statique, valeurs dynamiques qui sont obtenues comme indiqué ci-dessus au moyen d'un système de positionnement type GPS. La mesure faite par le capteur 10 est transmise, via des lignes de communication, en entrée de l'unité de calcul 6. L'unité de calcul 6 exécute alors les instructions de la partie logicielle de l'invention avec une période égale à la période d'échantillonnage du capteur, à savoir Êt=10ms. L'exécution de ces instructions permet à chaque instant k, d'abord de calculer la probabilité d'occupation d'une région de l'espace par un objet, puis d'associer à cette région de l'espace occupée un niveau de danger. Pour éviter une collision, l'utilisation de l'information correspondant au niveau de danger, peut être faite sous la forme d'une alarme visuelle et/ou sonore indiquant la présence de l'objet au conducteur ou sous la forme d'une commande d'actionnement des freins pour atteindre une vitesse de consigne permettant d'éviter le danger détecté. II est nécessaire que l'exécution de la partie logicielle selon l'invention se fassent entièrement pendant la période At, c'est-à-dire quasiment en temps réel. Cette spécification impose des contraintes sur l'algorithme de la présente invention compte tenu des puissances de calcul disponibles actuellement. Typiquement, l'unité de calcul 6 est équipée d'un processeur du type Pentium IV fonctionnant à une fréquence de 2 GHz et d'une capacité mémoire de l'ordre Go. B Notion de grille d'occupation Le domaine spatial dans lequel on cherche à déterminer la présence d'un objet correspond à une zone rectangulaire en avant du véhicule 1. Cette zone à deux dimensions est repérée par des coordonnées cartésiennes: X, dans la direction longitudinale 1 et X2 dans la direction transversale. Pour obtenir une grille G, ce domaine spatial est numérisé. Il est découpé en cellules ayant un pas régulier prédéterminé AX1 selon la direction X, et AX2 selon la direction X2. Il s'agit donc de la discrétisation régulière et statique de l'environnement du véhicule sous la forme d'un ensemble de N cellules X (X,, X2). La grille à deux dimensions actuellement utilisée possède de l'ordre de 100 000 cellules. Dans le mode de réalisation de la figurel, elle s'étend typiquement de 20 m selon X, et de 20 m selon X2. Comme cela est représenté sur la figure 1, le domaine spatial ne se superpose pas exactement à la zone de détection couverte par le capteur 10. On obtient alors une grille d'occupation en associant à chacune des cellules de la grille un nombre représentant la probabilité qu'un objet occupe cette cellule. On cherche alors à déterminer la probabilité que la cellule x soit dans l'état occupé EkX 1 à l'instant k, compte tenu de la série d'observations z zk réalisées au cours du temps. Cette probabilité d'occupation est notée P(EkXIz ...zk x). Selon une convention classique du domaine des probabilités, P(AIB) signifie la probabilité d'obtenir A, la valeur de B étant connue. Il est à noter que l'observation zk n'est pas directement la mesure en p et c donnée par le capteur, mais correspond plutôt à cette mesure rapportée au niveau de la cellule: si une mesure indique la présence d'un objet en p et cl), zk vaut 1 pour la cellule xk associée aux coordonnées p et c, et vaut 0 ailleurs. Une hypothèse est faite quant à l'indépendance de l'état EX d'une cellule de la grille par rapport aux autres cellules, et en particulier des cellules voisines. Cette hypothèse a pour but de réduire les corrélations entre cellules et donc, de simplifier le calcul de la probabilité d'occupation sur l'ensemble des N cellules X de la grille G. En effet, cela revient alors à répéter N fois le calcul permettant d'estimer la probabilité d'occupation d'une seule cellule. Avantageusement, l'algorithme selon l'invention est codé de manière à être entièrement parallélisable. Sa mise en ceuvre peut être faite au moyen d'une carte électronique comportant N processeurs fonctionnant en parallèle. Chaque processeur calculant la probabilité d'occupation de la cellule de la grille qu'il représente. L'algorithme selon l'invention se subdivise de manière schématique en plusieurs modules: - le premier module consiste à modéliser la réponse du capteur. Ce modèle capteur se fonde sur une distribution de probabilité P(ZkIEXk X) qui indique la probabilité de l'observation zk connaissant l'état EkX de la cellule X à l'instant k. Il est à noter que les majuscules indiquent une variable pouvant prendre un ensemble de valeurs possibles indiquées en minuscule. Par instanciation d'une variable A, il faut entendre la valeur particulière a que prend cette variable A. -le deuxième module est un module de prédiction. Il permet de prévoir l'état EkX de la cellule X à l'instant k compte tenu de l'état de la grille à l'instant k-1. - le troisième module est un module d'estimation permettant de calculer la probabilité d'occupation EkX de la cellule X à l'instant k compte tenu des observations Zk à l'instant k, et du résultat de l'étape de prédiction obtenu en sortie du deuxième module. - un quatrième module permet enfin d'associer aux cellules de la grille dont la probabilité d'occupation est élevée une grandeur indiquant un risque à collision et d'émettre un signal adapté. C Modélisation d'un capteur télémétrique L'information récupérée par le capteur est limitée en précision. De nombreux facteurs affectent la précision de la mesure: température, éclairage, couleur de l'objet, etc. Pour exploiter au maximum les mesures, il est important de tenir compte du bruitage de la grandeur observée. Selon l'invention, la réponse du capteur est décrite par un modèle probabiliste spécifié sous forme d'un programme bayésien. Les variables pertinentes du problème sont l'observation Z; la cellule X; la variable EX signifiant qu'un objet existe dans la cellule X; et, la variable DX signifiant qu'un objet a été détecté dans la cellule X. La conjonction des variables EX et DX permet d'expliciter les quatre situations possibles suivantes: -[EX =1] A [DX =1] : un objet occupe effectivement la cellule X, et il a été détecté par le capteur; - [EX = 0] A [DX = 0] : la cellule est vide, et le capteur n'y détecte effectivement rien - [EX =1] A [DX = 0] : un objet occupe la cellule X, mais il n'est pas détecté par le capteur. C'est par exemple le cas lorsque la cellule X est en dehors de la zone de détection du capteur, que la cellule X est masquée par un autre objet, ou encore lorsque le capteur est défaillant; -[EX = 0] A [DX =1] : la cellule X est vide et pourtant le capteur y détecte un objet. Ces deux dernières situations correspondent à une défaillance du capteur. Le probabilité conditionnelle P(Z, X, EX,DX) , exprimant la conjonction des différentes variables sélectionnées, se décompose avantageusement sous la forme P(X) et P(EXIX) représentent des connaissances a priori sur l'environnement du 35 véhicule. Dans un premier temps, afin de ne privilégier aucune situation particulière, ces distributions sont choisies uniformes. Mais, comme cela sera explicité ci-dessous, P(EXIX) permettra de tenir compte de l'historique de la scène. La probabilité P(DXIEX X) représente la capacité du capteur à détecter ou non une cible, ou bien à générer des fausses alarmes: - P([DX =1] [EX =1], X) = PD (X) représente la probabilité de détection. - P([DX =0] [EX =1],X)=1-PD(X) représente la probabilité que le capteur manque une cible existante. - P([DX =1] [EX = 0],X) = PFA(X) représente la probabilité de fausse alarme. 10 - P([DX =0] [EX =0],X)=1 PFA(X). Ces distributions de probabilités dépendent du capteur et des conditions dans lesquelles il est utilisé. PD(X) et PFA(X) peuvent être estimées à partir d'une modélisation des caractéristiques physiques des capteurs et des objets détectés. Elles peuvent également être "apprises" expérimentalement dans les conditions particulières d'utilisation du capteur au cours d'une première étape de calibration. Il reste à définir les différentes formes paramétriques associées à la famille de distribution P(Z 1 DX, EX, X) suivant les valeurs de Dx et de E. P(Z [DX =1][EX =0],X) et P(Z [DX =1][EX =1],X) indiquent qu'une détection a eu lieu pour la cellule X. Il s'agit d'une autre formulation de la fonction de vraisemblance du capteur. Le plus souvent une famille de distributions normales est utilisée: P(Z [Dx =1] Ex X) = .nr(Z#(X)#D(X) ). où, pour une valeur de X donnée, la moyenne p(X) de la distribution normale représente la réponse attendue du capteur. La matrice de covariance E(X) permet de représenter les variations possibles de la réponse du capteur autour de la réponse attendue. En ce sens, elle représente la précision du capteur. Les paramètres du modèle capteur, comme la moyenne et la matrice de covariance dans le cas Gaussien, peuvent être acquis au cours d'une deuxième étape de calibration du capteur. P(Z[DX =0][EX = 0],X) indique que la cellule X est vide, et qu'aucune détection n'a eu lieu dans cette cellule. En considérant uniquement ces informations, il est difficile de dire quelle peut être la valeur de l'observation Z dans une autre cellule. On choisit une distribution uniforme. 2890774 9 P(Z [DX = 0][EX =1],X) indique que la cellule X est occupée par un objet, et qu'aucune détection n'a eu lieu. Soit un objet masque la cellule X et a été détecté par le capteur, soit l'objet présent dans la cellule X n'est pas masqué, et n'a pas été détecté. On décide que la probabilité pour la mesure Z d'être dans une zone masquée par la cellule X est différente de celle des zones non masquées. Une fois les paramètres du modèle donnés, la description précédente peut être utilisée pour calculer la probabilité capteur: _ EDx P(ZxeDx) P(e,, P(ZIe,x) P(,x)P(e,. x)EDx (Dx I,x)P(Z Dx e} =1 E P(Dx e, ,x)P(Z Dx nx D Grille d'occupation équivalente et compression Dans le mode de réalisation actuellement envisagé de l'invention, la grille d'occupation est une grille spatiale. En conséquence, pour chaque valeur de la vitesse (deux dimensions) il faut stocker une grille d'occupation à deux dimensions. Il s'ensuit qu'il est nécessaire de stocker un grand nombre de données pour disposer d'une représentation précise de l'environnement, i.e. surveiller de larges zones de l'environnement avec un pas de discrétisation fin pour y distinguer les objets. Outre, le problème de la taille mémoire requise, la vitesse de traitement des données est un paramètre crucial dans une application temps réel. Selon l'invention, on utilise une structure de données représentant la grille d'occupation mais qui possède une taille mémoire réduite. Dans un premier temps, nous allons décrire la manière dont la probabilité d'occupation d'une cellule est mise à jour à l'instant k à partir d'une nouvelle observation. Soit Z une variable aléatoire dont l'ensemble des événements possibles correspond à l'ensemble des mesures possible d'un capteur. Soit P(ZIEX) la distribution de probabilité conditionnelle sur Z sachant l'état d'occupation EX de la cellule considérée x. On a alors la distribution conjointe P(,x)P(e, x) 1 P(Dx e, ,x)P(Z Dx e, z, Dx P(Z,EX) = P(Ex) P(ZIEX), où p(Ex) correspond à un a priori sur le fait que la cellule x est occupée. Pour une nouvelle mesure zk et une valeur d'occupation ex de la cellule x, la règle de Bayes donne: p(ex Zk) p(ex)p(z) ex) p(z avec: p(zk) = p(ex =1)p(zk ex =1) + p(ex = 0)p(zk ex = 0) La mise à jour de la probabilité d'occupation de cellule x correspond à : k ex pk+1(ex) = p(ex zk) e z Lpk(ex)p(z ex) (1) Lors de la construction d'une grille d'occupation, la série d'observation 10 Zk ={z ,z1,...,zk} est perçue par le système, et on note 0 1 k-1 pk(ex)=p(exz,z,...,z). On utilise cette chaîne de mise à jour avec la distribution initiale qui est choisie uniforme puisque aucune information à propos de l'environnement n'est a priori disponible. Comme la variable Ex est binaire, on a la relation: pk (ex = 1) = 1 pk (ex = 0) . Il est alors équivalent de stocker en mémoire lpour chacune des cellules de la grille la variable qk définie par: 1k = pk(ex 1) En injectant la relation (1) dans la définition ci-dessus, il vient la relation 20 recursive: _ p(zk 1 e =1) k-1 p(z ex =1) qk k-1 x qk1 q0 i p(z ex = 0) i=o p(z ex = 0) Avantageusement le logarithme de la variable qk peut être pris en compte de manière à transformer le produit en une somme dans la relation précédente: k-1 p(z e =1) log(qk) =1og(go)+L i=o p(zi ex = 0) La mise à jour de la variable log(qk) pour une cellule donnée x, et par voie de conséquence, la mise à jour d'une grille d'occupation équivalent (chaque Pk (ex 0) cellule de cette nouvelle grille comportant la valeur numérique de log(gk)) est alors une simple opération d'addition à chaque pas temporelle. L'avantage de passer par la variable log(gk), c'est qu'il est alors facile de compresser la grille d'occupation équivalente en effectuant une transformée en ondelettes, par exemple au moyen d'ondelettes de Haar, puis en supprimant les coefficients de faible amplitude. La notion de transformée en ondelettes est largement décrite dans la littérature et a fait l'objet de nombreuses publications aussi bien scientifiques que de vulgarisation. On notera qu'avantageusement l'espace des signaux transformés en ondelettes est un espace vectoriel, qui est en tant que tel stable par les opérations d'addition et de multiplication par un scalaire. De plus, la transformation en ondelettes est une opération linéaire par rapport au nombre de données, ce qui garantie une exécution rapide de l'algorithme de compression. L'algorithme suivant est actuellement envisagé pour mettre en oeuvre cette compression par ondelettes de la grille d'occupation équivalente: acquérir les données provenant du capteur à l'instant k sur tout le champ de vue; - construire dans le champ de vue, la grille des logarithmes t' (zk-1 ex 1) 1og(gk (x)) P(zk l e = 0) x - trouve les carrés c de support de fonction d'échelles de Haar maximaux inclus entièrement dans el champ de vue du capteur; - Faire une transformée en ondelettes de la fonction log(gk(x)) sur chacun des carrés précédemment trouvés: W, - Ajouter chaque Wc à la grille d'occupation globale de Haar -Si le nombre de coefficients est trop important, effectuer une compression non linéaire en enlevant tous les coefficients dont la norme est négligeable par rapport à la norme de toute la représentation. E. Cas multicapteur Il est avantageux d'utiliser plusieurs capteurs pour accroître la robustesse des observations. On peut ainsi s'affranchir des limitations propres à chaque capteur (par exemple un capteur laser fonctionne moins bien par temps de pluie), limiter la dépendance vis-àvis des défaillances d'un capteur (plusieurs capteurs du même type peuvent être utilisés), et surtout augmenter la précision de la détection d'un objet dans l'environnement en utilisant la redondance des informations. Lorsque S capteurs sont utilisés, chacun permettant une mesure zks à l'instant k, l'hypothèse est faite que, connaissant l'état d'une cellule, les observations des différents capteurs sont indépendantes. Alors la probabilité conjointe se décompose selon: Pi où P(Zj EX, X) est le modèle du ième capteur tel que décrit en détail précédemment. Etant donné une série d'observations à l'instant k, Z, ... Zs, la probabilité que la cellule X soit occupée devient: F Grille d'occupation et filtres bayésien Selon l'invention, les modules de prédiction et d'estimation d'une probabilité d'occupation à l'instant k mettent en ceuvre la notion de filtre bayésien appliqué aux grilles d'occupation. Un filtre bayésien permet de prédire l'évolution temporelle d'un système. Le but est ici de conférer à la grille d'occupation une certaine robustesse lui permettant de tenir compte d'une défaillance d'un capteur ou de l'occultation passagère d'un objet dans l'évaluation du risque. D'une manière schématique, nous voulons estimer de manière récursive la probabilité d'occupation d'une cellule à l'instant k. Ceci est fait en deux étapes. Une étape de prédiction dont le but est d'estimer a priori l'état d'occupation de la cellule. Une étape d'estimation dont le but est de calculer la probabilité d'occupation en prenant en compte l'estimation a priori et l'observation instantanée. Les variables considérées sont: - X qui précise sur quelle cellule nous travaillons; - Ekx qui représente l'état de la cellule X à l'instant k; z k:k qui z = z q représentent la série des observations depuis un instant initial Pour l'étape de prédiction, on montre que la probabilité de prédiction s'écrit: JP(E,X (El-,1,X') où P(EX,,X' z Ek 1) est la probabilité d'occupation estimé de la cellule X' à l'instant k-1; où P(EX, X EX,',X') est la probabilité de transition d'un état intermédiaire (EX71, X') quelconque du système à l'instant k-1 vers l'état considéré (EX,X), i. e. correspond au modèle dynamique du système; et où l'on somme sur l'ensemble des états intermédiaires possibles (EX-1,X'). Pour l'étape d'estimation, on montre que la probabilité d'estimation s'écrit: P(EX,X ZO:k) = -P(Zk EX,X)P(EX,X Z0:k-1) a (3) où apparaît la probabilité du modèle capteur P(Zk TEX X) et la probabilité de prédiction à l'instant k: P(EX,X zo:k-1) Ces deux équations permettent d'avoir unestructure du type boucle prédiction/estimation. G - Discrétisation non continue des dimensions de vitesse: Tout ce qui vient d'être décrit est facilement transposable à une grille d'occupation à plus de 2 dimensions. Il est possible de construire une grille à quatre dimensions, ayant les deux dimensions d'espace X, et X2, comme ce qui vient d'être décrit, auxquelles on ajoute deux dimensions de vitesse X,' et X2'. Si le capteur ne donne pas directement une mesure de la vitesse de l'objet détecté, celle-ci peut être calculée à partir de la variation de la position spatiale de l'objet =Ex,/Et et X2'=LIx2/LIt, par exemple, en prenant le temps d'échantillonnage du détecteur entre deux mesures successives de la position de l'objet comme intervalle de temps At. Les dimensions de vitesse sont ensuite discrétisées de manière régulière et continue avec un pas EX,' selon X,' et EX2' selon X2'. Dans cette approche où l'espace des vitesses est discrétisé de manière continue survient le problème suivant. Soit une grille spatiale associée à une valeur de vitesse X'o(X,'o et X2'o), ci après dénommé plan de vitesse. Soit une cellule de ce plan de vitesse occupée à l'instant k-1 par un objet. Si la vitesse X'o ne correspond pas à un multiple entier du pas de discrétisation spatiale EX divisé par l'intervalle de temps At, alors, à l'instant k, l'objet n'occupera pas une unique cellule du plan de vitesse mais mordra sur un ensemble de cellules. Or, une cellule de la grille sur laquelle mord un objet est considérée comme P(EX,X z0:k-1) X E 1, X X')P(E,X' zo:k 1) (2) occupée. A l'intervalle de temps suivant, ces cellules occupées seront propagées en augmentant encore le volume des cellules occupées. Ainsi l'erreur sur l'état de la grille, due essentiellement à la discrétisation, est amplifiée lors de l'étape de prédiction. L'erreur de discrétisation est croissante et n'est seulement réduite que par la phase d'observation. Selon l'invention, dans un premier temps, l'espace cartésien à deux dimensions est discrétisé de manière continue, puis, dans un second temps, le nombre de plans de vitesse possibles est restreint. On ne considère que les plans de vitesse pour lesquels la vitesse peut s'écrire, en fonction des pas de discrétisation Ax, et Ax2 et de l'intervalle d'échantillonnage At, de la manière suivante: X', = p/n (Axi/At) et X'2 =q/n (Ax2/At) avec p et q appartenant à l'ensemble des entiers relatifs. Ainsi, les mouvements autorisés sont ceux qui font exactement passer une cellule à une autre cellule en un nombre n de fois l'intervalle d'échantillonnage At. Un plan de vitesse est une grille spatiale à deux dimensions qui est associé à un triplet (p,q,n) définissant une vitesse. On a certes un espace des vitesses discrétisé, mais qui n'est pas continu. Alors, pour calculer l'intégrale sur l'ensemble des cellules qui intervient dans la détermination de la probabilité d'occupation prédite, une hypothèse de vitesse constante est faite. Selon cette hypothèse, entre deux instants successifs, la vitesse d'un objet est constante et la probabilité de transition est nulle d'un plan de vitesse à un autre. C'est-à-dire que, si la cellule Xk n(X,, X2) du plan de vitesse X'(p/n (Axi/At), q/n (Ax2/At)) est occupée à l'instant k-n, alors la cellule Xk(X,+ p/nAX1; X2+ q/n Ax2) sera occupée à l'instant k. Inversement, lorsqu'on calcule la probabilité d'occupation prédite de la cellule Xk, par plan de vitesse possible (p,q,n), seule la probabilité d'occupation estimée à l'instant k-n de la cellule Xk-n antécédent de la cellule Xk intervient dans le calcul. Avantageusement, à chaque instant d'observation, le facteur 1/n permet de ne considérer le plan de vitesse (p,q,n) qu'à une fréquence 1/n. Ainsi pour un plan de faible vitesse, le calcul de l'intégral ne se fera qu'à une fréquence faible. En revanche, un plan dans lequel se situe un objet à grande vitesse sera observé à une fréquence plus élevée de manière à prendre en considération l'évolution rapide du comportement de cet objet. La mise en ceuvre de l'invention se fait de la manière suivante: Dans un premier temps, une région spatiale (2D) autour d'une cellule centrale X, type est délimitée; un nombre maximum nmax de pas de temps est défini. L'ensemble des 5 cellules de cette région sont celles qui peuvent atteindre la cellule centrale X, en 1, 2, ... nmax pas de temps. L'ensemble des vitesses Vf possibles de la forme (p/n (Axi/At) , q/n (0x2/Ot)) qui permettent à l'une des cellules de cette région de se propager vers la cellule centrale Xc en au plus nmax pas de temps est déterminé ; L'ensemble Vf est ensuite découpé en sous-ensembles de vitesses permettant à une cellule de rejoindre en n pas de temps la cellule centrale: Vf =UVfi avec Vf ={ v=(p,q,n)E Vf n=i} Puis, à chaque élément (p,q,i) appartenant à l'un des Vfi est associée une grille spatiale à deux dimensions GZ(p,q). Alors l'antécédent d'une cellule (I, m) de 15 cette grille est la cellule (I-p, m-q). Un compteur k, donnant l'indice du pas de temps, incrémenté d'une unité à chaque itération évolue entre la valeur nulle et la valeur correspondant au plus petit commun multiple des nombres 2 et nmax. (En effet, il faut considérer le plus petit commun multiple puisqu'il s'agit d'un ensemble de boucles imbriquées les unes dans les autres). En début de processus, il est instancié à la valeur nulle; La probabilité d'occupation de la cellule Puis pour i allant de 1 à nmax, si k est égal à 0 modulo i, alors, pour chacune des grilles d'occupation G1(p,q) d'un sous-ensemble Vfi Déterminer la cellule antécédent de la cellule considérée; -attribuer la probabilité d'occupation estimée donnée par GZ(p,q) pour ledit antécédent en tant que valeur de la probabilité d'occupation prédite pour la cellule considérée (l'équation 2 ci-dessus) ; - mettre à jour la grille en déterminant la probabilité d'occupation estimée à partir de la probabilité d'occupation prédite et de la probabilité capteur tenant compte de l'observation à ce pas de temps; - puis incrémenter k d'une unité et retourner au début de la boucle. H Cas multiobjets Dans ce qui précède, seul a été décrit formellement le cas d'un unique objet présent dans l'environnement du véhicule. Dans le cas d'un environnement complexe, il est nécessaire de pourvoir tenir compte de la présence de plusieurs objets dans cet environnement. Lorsque plusieurs objets sont détectés, le capteur réalise O observations à l'instant k: zk,, zk2... zko. On introduit alors la variable M dite d'association. M est un entier entre 0 et O. Si M=O, alors aucune des observations n'est utilisée pour mettre à jour la probabilité d'occupation de la cellule X. en revanche, si M=i, alors la ième observation est utilisée pour mettre à jour la probabilité d'occupation de la cellule X. La probabilité conjointe se décompose de la manière suivante: 0 X P( ,rx 1 [f = 2890774 15 la valeur de la variable M étant indépendante de l'état EX et de la cellule X considérée, P(MI Ex X) = P(M) ; une simplification a été opérée en faisant une hypothèse d'indépendance conditionnelle des différentes observations z; entre elles connaissant l'état de la cellule considérée. Pour l'observation z;, si M est égal à i, alors la distribution de probabilité à considérer est le modèle capteur dérivé ci-dessus; en revanche, si M est différent de i, alors une distribution uniforme est choisie (U). Dans la présente application, on cherche à déterminer: P(EXIzi...zo x), qui s'écrit, en utilisant la variable M: P(L,. {w. x) NE, ;x Cette relation est une relation statique, qui est rendu dynamique en introduisant la boucle estimation/prédiction décrite précédemment sur un jeu de variables entre un instant initial 0 et un instant présent k. En particulier la variable d'association Mk, indexée par le temps, permet d'associer à l'estimation de l'état, l'une des Ok observations faites à l'instant k. Dans cette modélisation on retrouve l'équation de la prédiction indiquée ci-dessus. En revanche, la relation de l'estimation devient: p(n xk x,:k) (P(mk) n a M.k---O 1'.,. km.kEx xk) x F'(El xk(,:x} Pour éviter un objet, il faut combiner au calcul de la probabilité d'occupation d'une cellule P(EkXXIZk) avec une probabilité de risque à collision P(DkXIEkXX). Ce risque à collision est calculé à partir d'un critère de temps à collision. Il devient alors possible de réguler la vitesse d'un véhicule Bien que l'invention ait été décrite en référence à un mode de réalisation particulier, elle n'est nullement limitée à ce mode de réalisation. Elle comprend tous les équivalents techniques des moyens décrits ainsi que leurs combinaisons qui entrent dans le cadre de l'invention. REVENDICATIONS 1. Procédé d'assistance à la conduite d'un véhicule comportant au moins un capteur apte à émettre un signal de sortie indiquant la présence d'au moins un 5 objet dans une zone de détection, ledit procédé comportant les étapes consistant à : - définir une grille correspondant à la discrétisation en cellules d'un espace des paramètres cinématiques dudit objet par rapport audit véhicule; et, à chaque instant k d'échantillonnage: - déterminer une observation zk à partir dudit signal de sortie du capteur; et, pour chaque cellule x de la grille: -déterminer, à partir de ladite observation zk, la valeur d'une probabilité capteur P(zklEkxX) modélisant le comportement dudit capteur; - calculer une grandeur d'occupation estimée à l'instant k fonction de la probabilité P(EkxXlz .k=Z...zk) que ladite cellule x soit dans l'état Ekx à l'instant k, étant données ladite observation zk; ladite grandeur d'occupation estimée à l'instant k étant fonction de ladite probabilité capteur P(zklEkxX) et d'une grandeur d'occupation prédite P(EkxXIz :k-1=z ...zk-1) à l'instant k; - déterminer ladite grandeur d'occupation prédite à l'instant k P(EkxXIz:k- 1=z zk-1) fonction de la probabilité que la cellule X considérée soit dans l'état occupé Ekx à l'instant k, en prenant en compte les observations depuis un instant initial 0 jusqu'à l'instant k-1; ladite grandeur d'occupation prédite à l'instant k étant déterminée en approchant une relation du type: P(EX, X ZO:k-1) f P(EX, x Exk71, X')P(EX,, X' Zo:k 1) (Ex,i,X') dans laquelle intervient la probabilité d'occupation estimée à l'instant k-1 pour chacune des cellules X' de la grille, en considérant que la probabilité de transition d'une cellule quelconque vers la cellule considérée est non nulle uniquement pour une cellule dite antécédent de ladite cellule considérée; - évaluer, pour chacune des cellules de la grille, une probabilité de risque à collision à partir de ladite grandeur d'occupation estimée à l'instant k; - émettre un signal de commande d'un actionneur en fonction des probabilités de risque à collision des cellules de la grille de manière à éviter ledit objet détecté. 2. Procédé selon la revendication 1, caractérisé en ce que ladite grandeur d'occupation est égale à ladite probabilité d'occupation. 3. Procédé selon la revendication 1, caractérisé en ce que ladite grandeur d'occupation est compatible avec une transformation de compression. 4. Procédé selon la revendication 1, caractérisé en ce que, une grille d'occupation étant définit qui associe à chacune des cellules de la grille une première grandeur fonction de la probabilité d'occupation de ladite cellule, une transformation de compression est appliquée sur ladite grille d'occupation pour obtenir une représentation équivalente de ladite grille d'occupation, ladite grandeur d'occupation étant une deuxième grandeur d'occupation associée à ladite représentation équivalente de ladite grille d'occupation. 5. Procédé selon la revendication 3 ou la revendication 4, caractérisé en ce que ladite première grandeur d'occupation estimée à l'instant k est le logarithme du rapport entre la probabilité d'occupation estimée à l'instant k que la cellule considérée soit occupée et la probabilité d'occupation estimée à l'instant k que la cellule considérée soit vide. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que ladite grille comporte plus de 10 million de cellules. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que, ladite grille étant une grille spatiale obtenue en discrétisant avec un pas prédéfini l'espace des positions dudit objet, ledit procédé comporte une étape supplémentaire de détermination d'un ensemble de vitesses possibles, lesdites vitesses possibles étant celles qui permettent de passer exactement d'une cellule de la grille à une autre cellule de la grille en un nombre fini n de pas d'échantillonnage. 8. Procédé selon la revendication 7, caractérisé en ce qu'une grille spatiale est associée à chacune desdites vitesses possibles, et en ce que lesdites étapes dudit procédé sont réalisées pour chaque grilles prise indépendamment des autres, une hypothèse de vitesse constante dudit objet détecté étant faite. 9. Procédé selon l'une des revendications 7 à 8, caractérisé en ce que ladite grille spatiale comporte plus de 1 million de cellules. 10. Procédé selon l'une des revendications précédentes, caractérisé en ce que ladite période d'échantillonnage est inférieure à 10ms. 11. Dispositif d'assistance à la conduite d'un véhicule, caractérisé en ce qu'il comporte des moyens aptes à la mise en ceuvre du procédé d'assistance à la conduite selon l'une des revendications 1 à 10. 12. Dispositif selon la revendication 11, caractérisé en ce qu'il comporte N processeurs, N étant le nombre de cellules de ladite grille, chacun desdits processeurs fonctionnant en parallèle desdits autres processeurs et calculant à chaque instant k d'échantillonnage la grandeur d'occupation estimée de la cellule X qu'il représente.