FR 2929394 A1 20091002 FR 0801789 A 20080401 Planification de chemins en présence de courants forts L'invention concerne un procédé de planification de chemins pour un mobile entre un point de départ et un point d'arrivée situés dans un espace donné en présence de courants, notamment en présence de courants forts, c'est-à-dire des courants dont la vitesse est supérieure à celle du mobile se déplaçant dans cet espace. L'invention s'applique à la recherche de chemins possibles ou optimales pour des mobiles se déplaçant dans des conditions de courant non uniforme, les engins mobiles étant par exemple des avions sans pilote, également appelés drones, ou des robots autoguidés sous-marins. Pour les premiers, les courants impliqués sont des courants aériens alors que pour les seconds, il s'agit de courants marins. Dans les deux cas, les engins mobiles se déplacent généralement à des vitesses relativement faibles, ces vitesses pouvant être inférieures aux vitesses des courants du milieu dans lequel ils évoluent. Par conséquent, les courants doivent impérativement être pris en compte afin de fournir un chemin physiquement réalisable compte tenu de la vitesse de l'engin mobile. Actuellement, les procédés de planification de chemins en présence de courants utilisent principalement des algorithmes génétiques et des algorithmes de propagation d'onde. Les algorithmes génétiques miment le processus de sélection naturelle proposé par Darwin. Dans le présent contexte, les chemins jouent le rôle d'individus d'une population. Les meilleurs d'entre eux, c'est-à-dire ceux dont le temps de parcours est le plus court, sont sélectionnés pour survivre et se reproduire au cours des générations. Cependant, les algorithmes génétiques souffrent d'un problème majeur dit de convergence prématurée. Cela signifie que, de façon difficilement contrôlable, l'ensemble de la population peut s'appauvrir el: converger vers un individu unique, appelé super-individu. Toutefois, ce super-individu peut être de piètre qualité, le chemin présentant par exemple de nombreux détours, ou même s'avérer physiquement irréalisable. De nombreuses solutions ont été proposées, telles le front de Pareto, la formation de sous-populations stables dans le voisinage des solutions optimales, appelée méthode des niches, et la méthode d'élitisme ; mais ces solutions ne font que limiter le problème sans le résoudre. Par ailleurs, les algorithmes génétiques nécessitent le réglage de nombreux paramètres, notamment la taille de la population, le nombre de générations, le type d'opérateurs et les probabilités associées. Ce réglage, généralement effectué de façon empirique, est fastidieux et ne fournit aucune garantie sur le comportement de l'algorithme. Les algorithmes de propagation d'onde sont utilisés depuis longtemps pour la planification de chemins en l'absence de courants. Ils furent initialement introduits par Raymond JARVIS et sont notamment présentés dans JARVIS, R. Collision-Free Trajectory Planning Using Distance Transforms, Mechanical Engineering Transactions, Journal of Institution of Engineers, Vol. ME10, No. 3, 1985, p. 187-191. Selon ces algorithmes, l'espace est discrétisé en un nombre fini d'éléments appelés cellules. Chaque cellule comprend un point de passage possible pour le chemin. Le point de passage correspond généralement au centre de la cellule. Lors d'une première étape de l'algorithme de propagation, tous les points de passage des cellules adjacentes à la cellule contenant le point de départ sont évalués en termes de temps de parcours, également appelés coûts. Lors des étapes ultérieures, on considère le point de passage pour lequel le coût est minimal et on évalue le temps de parcours entre ce point de passage et les points de passage des cellules adjacentes non encore explorées, une cellule étant considérée comme explorée lorsqu'un coût lui a déjà été attribué. Le temps de parcours entre le point de départ et les points de passage des cellules adjacentes non explorées, autrement dit le coût de ces cellules, peut ainsi être déterminé par une simple addition. Dans le cas où plusieurs points de passage ont un coût égal, on considère successivement ces points de passage. A chaque étape, il est possible de distinguer les cellules explorées des cellules non explorées. La frontière entre les cellules explorées et les cellules non explorées forme un front d'onde à partir duquel les points de passage suivants peuvent être évalués. En d'autres termes, le front d'onde se propage progressivement dans les cellules, jusqu'à atteindre la cellule contenant le point d'arrivée. Le chemin optimal est ensuite construit en partant du point d'arrivée, en utilisant un algorithme couramment appelé "hill climbing" dans la littérature anglo-saxonne. Cet algorithme est l'équivalent discret de la méthode du gradient et consiste à se déplacer de proche en proche au sein des cellules jusqu'à atteindre le point de départ, chaque déplacement étant déterminé afin de :3 choisir la cellule de coût minimal dans le voisinage direct de la cellule courante. La propagation d'onde a ensuite été adaptée pour pouvoir prendre en compte la présence de courants. Une première adaptation est proposée dans PETRES, C., Y. PAILHAS, J. EVANS, Y. PETILLOT et D. LANE. Underwater path planning using fast marching algorithms, Proceedings of Oceans 2005 ù Europe, Vol. 2, Brest, France, juin 2005, p. 814-819. Une deuxième adaptation est décrite dans GARAU, B., A. ALVAREZ et G. OLIVER, Path planning of autonomous underwater vehicles in current fields with complex spatial variability: an A* approach, Proceedings of the IEEE Conference on Robotics and Automation (ICRA 2005), Barcelone, Espagne, avril 2005, p. 546-556. Une troisième adaptation est proposée dans SOULIGNAC, M. Planification de trajectoires pour engins volants autonomes, Rencontre Nationale des Jeunes Chercheurs en Intelligence Artificielle, p. 195-212. Ces adaptations donnent des résultats satisfaisants en présence de courants modérés, mais elles ne sont pas fiables en présence de courants forts. En particulier, les première et deuxième adaptations présentent des problèmes d'incorrection, l'incorrection étant définie comme la fourniture d'une solution physiquement irréalisable compte tenu de la vitesse du mobile et de celle du courant. La troisième adaptation a permis de corriger ce problème mais présente un problème d'incomplétude. L'incomplétude est définie comme l'échec de fourniture d'un chemin alors qu'au moins un chemin valide existe. En résumé, aucun procédé de planification de chemins en présence de courants ne permet de trouver, à coup sûr lorsqu'il existe, un chemin à la fois optimal et physiquement réalisable. Dans le contexte de l'invention, on entend par chemin optimal, le chemin parcouru en un minimum de temps. Un but de l'invention est notamment de pallier tout ou partie des inconvénients précités en proposant un procédé de planification de chemins en présence de courants exempte des problèmes d'incorrection et d'incomplétude. A cet effet, l'invention a pour objet un procédé de planification de chemins pour un mobile entre un point de départ et un point d'arrivé situés dans un espace donné en présence de courants, ces courants pouvant varier en vitesse et en direction dans l'espace, le mobile étant soumis aux courants et se déplaçant à une vitesse donnée par rapport aux courants. Selon l'invention, le procédé comprend les étapes suivantes : - discrétiser l'espace en zones de courant élémentaire, chaque zone étant délimitée par des bords, - associer à chaque bord un point de passage possible pour le mobile, chaque point de passage pouvant translater sur son bord, - déterminer une séquence de points de passage permettant au mobile de relier le point de départ au point d'arrivée en un minimum de temps. L'invention a notamment pour avantage qu'elle permet de trouver, même en présence de courants forts, un chemin optimal entre un point de départ et un point d'arrivée, dès lors qu'au moins un chemin physiquement réalisable existe. Selon une forme particulière de réalisation, une zone de courant élémentaire est une zone de forme polygonale, indivisible et dans laquelle le courant est considéré comme homogène sur la totalité de la surface de la zone. Les zones de courant élémentaire peuvent être construites en calculant le diagramme de Voronoï autour de points pour lesquels on connaît une vitesse et une direction de courant. Selon une forme particulière de réalisation, la détermination de la séquence de points de passage permettant au mobile de relier le point de départ au point d'arrivée en un minimum de temps utilise un algorithme de propagation d'onde adapté aux points de passage. Selon une forme particulière de réalisation, un chemin comprend une succession de segments de droite entre le point de départ et le point d'arrivée, une fonction de coût étant associée à chaque segment de droite entre deux points de passage voisins et définie par la relation : J t3 f+l = û cXl.dJ j 2 cy .d yJ + où d' (dx ; d ) représente le (v ) ù (ci )2 déplacement du mobile entre les points de passage voisins Vi et ccx;cÿ) représente le vecteur vitesse du courant dans la zone de courant élémentaire dans laquelle se déplace le mobile, vj(vx;vÿ) représente le vecteur vitesse du mobile par rapport au courant dans cette même zone de courant élémentaire et où A=(vj)2.((dx)2 +(dÿ 2)ù(cx.dÿ ùcÿ .dÿ Y. Cette forme particulière de réalisation a notamment pour avantage qu'elle permet de ne trouver que les seuls chemins physiquement réalisables, la fonction de coût tj'j+' prenant en compte la vitesse effective de déplacement du mobile par rapport à l'espace. Selon cette forme particulière de réalisation, la détermination de la séquence de points de passage permettant au mobile de relier le point de départ au point d'arrivée en un minimum de temps comprend la détermination de coûts Ti pour des points de passage V;, ces coûts T,. étant définis par la relation : T. 4,(4 j=0 Les coûts T, peuvent être déterminés en utilisant une méthode du gradient projeté paramétrée par un paramètre de tolérance E . L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description détaillée d'un mode de réalisation donné à titre d'exemple, description faite en regard de dessins annexés qui représentent : - la figure la, un premier exemple de discrétisation d'un espace donné ; - la figure lb, un deuxième exemple de discrétisation de 25 l'espace ; - la figure 2, deux zones de courant élémentaire comprenant des bords et des points de passage ; - la figure 3, l'ensemble des directions de propagation effectivement accessibles à partir d'un point de passage ; 30 - la figure 4, différentes étapes de l'algorithme de propagation d'onde selon l'invention. Le procédé de planification de chemins selon l'invention comprend plusieurs étapes. Les figures 1 a et 1 b représentent deux exemples différents résultant de la première étape du procédé. Pour la suite de la description, on considère un espace plan de dimensions finies. Cependant, le procédé selon l'invention peut être étendu à un espace à trois dimensions moyennant des adaptations exposées plus loin. Dans ces deux exemples, on considère un espace 1 carré délimité par une frontière 2 à l'intérieur de laquelle un mobile peut se déplacer. Le procédé selon l'invention s'applique néanmoins tout aussi bien à un espace polygonal. Dans tous les cas, un repère cartésien R(O,x,y), par exemple orthogonal, est associé au plan de l'espace 1. La répartition des courants dans cet espace 1 est connue. Pour différents points 3 de l'espace 1, on connait des données relatives aux courants, notamment une vitesse et une orientation du courant. Sur les figures la et 1 b, la vitesse et l'orientation du courant en chaque point 3 sont représentées par un vecteur 4 de courant centré sur le point 3 considéré. La vitesse et l'orientation du courant en chaque point 3 peuvent être déterminées par prévision ou par mesure. La distribution des points 3 dans l'espace 1 peut être régulière, comme représenté à la figure la ou irrégulière, comme représenté à la figure 1 b. Dans la première étape du procédé, l'espace 1 est discrétisé en zones 5 de courant élémentaire. Chaque zone 5 de courant élémentaire est délimitée par des bords 6i. L'indice j varie de 1 à N, N étant le nombre de bords 6j des zones 5 de courant élémentaire, et permet de repérer et différencier les différents bords 6j. Avantageusement, la discrétisation est fonction de la distribution des points 3, une zone 5 de courant élémentaire de forme polygonale étant définie pour chaque point 3. Dans chaque zone 5 de courant élémentaire, la vitesse et l'orientation du courant est supposée uniforme. Dans une forme particulière de réalisation, les zones 5 de courant élémentaire sont déterminées en calculant le diagramme de Voronoï autour des points 3. Le diagramme de Voronoï comprend un ensemble de segments de droites 7, chaque droite étant équidistante de deux points 3. Dans le cas général, c'est-à-dire dans le cas d'une distribution irrégulière des points 3, les segments de droite 7 forment des zones 5 de courant élémentaire de forme polygonale. Dans le cas particulier d'une distribution régulière des points 3, les zones 5 de courant élémentaire sont des carrés égaux. La longueur de côté des carrés correspond à la distance séparant deux points 3 consécutifs. Chaque carré est centré sur un point 3. La figure 2 illustre la deuxième étape du procédé selon l'invention. Selon cette étape, un point de passage Vi est associé à chaque bord 6i des zones 5 de courant élémentaire. Selon la forme particulière de réalisation dans laquelle les zones 5 de courant élémentaire sont déterminées en calculant le diagramme de Voronoï, les bords 6i des zones 5 de courant élémentaire correspondent aux segments de droite 7 de ce diagramme. Avantageusement, les bords 6i des zones 5 de courant élémentaire correspondant à une frontière 2 ne comportent pas de point de passage V. La position du point de passage Vi sur un bord 6i est repérée par une abscisse curviligne I~. Un chemin peut être déterminé par une séquence de Q points de passage Vi entre un point de départ Vo et un point d'arrivée Vo_1. En d'autres termes, un chemin peut être déterminé par une séquence de bords 6i et d'abscisses curvilignes I~ associées à chaque point de passage V. Les bords 6i peuvent être modélisées de la manière suivante : - pour chaque bord 6i, on détermine une équation de droite du type a1.x+bj.y+ci =0, deux bornes BJ =(xi ;y;) et B; =(x;;y,) correspondant aux points d'intersection des segments de droite 7 et un domaine Dj _ [0;11 pour les valeurs possibles de lb entre les bornes B; et B - pour chaque point de passage on peut repérer ses 25 coordonnées V. _ (x1 ; y1) dans le repère cartésien R(O, x, y) par les relations : x =x~ +,Jq(1+(a1/b.)) (1) yj = y; + (ù .x~ ù c~ )/b Par ailleurs, pour chaque bord 6i, on note A. = {ECAk;ECAI} la liste des zones 5 de courant élémentaire adjacentes au bord 6i, où ECAk et 30 ECA, sont des zones 5 de courant élémentaire, les indices k et I pouvant varier de 1 à m, m étant le nombre de zones 5 de courant élémentaire. La notion de voisinage entre deux points de passage peut être définie ainsi : deux points de passage V; et Vi sont voisins si et seulement si A. (1 Ai ~ { }, où { } désigne l'ensemble vide. Dans une troisième étape, on détermine une séquence de points de passage Vi permettant au mobile de relier le point de départ au point d'arrivée en un minimum de temps. En d'autres termes, on détermine simultanément une séquence de bords 6j et les abscisses curvilignes Ij associées à chaque point de passage V. Selon une forme particulière de réalisation, cette troisième étape utilise un algorithme de propagation d'onde, cet algorithme étant adapté aux points de passage Vi selon l'invention. Un point de passage associé à une cellule d'un algorithme de propagation d'onde classique est défini à l'avance, au moment du processus de discrétisation. Généralement, ce point de passage est fixé au centre de la cellule ou au milieu de l'un de ses côtés. Cet aspect statique est à l'origine des problèmes d'incomplétude évoqués précédemment. Si un point de passage est situé dans une zone non atteignable, par exemple à cause d'un courant opposé au déplacement du mobile, l'ensemble de la cellule est considéré comme inaccessible même si certaines régions de la cellule sont en réalité accessibles. Un point de passage Vi selon l'invention introduis: un degré de liberté supplémentaire, permettant de rendre accessible l'ensemble des directions de propagation effectivement accessibles à l'intérieur d'une zone 5 de courant, comme représenté à la figure 3. Sur cette figure, on considère un point de passage V; situé sur un bord 6;. L'ensemble des directions de propagation effectivement accessibles à partir du point de passage V1 est représenté par une zone située entre deux droites 31. Par conséquent, l'ensemble des points de passage possibles suivants sont situés sur les bords 6, et 6y appartenant à cette zone, les indices x et y appartenant à l'ensemble de nombres entiers compris entre 1 et N. Selon l'invention, si un chemin physiquement réalisable existe, il peut être calculé avec une précision arbitraire, cette précision étant directement liée à la finesse d'évaluation des points de passage V. L'algorithme de propagation d'onde est adapté aux points de 35 passage Vi selon l'invention, par exemple de la façon décrite ci-après. Cette adaptation est nécessaire dans la rnesure où les points de passage Vi possèdent un degré de liberté. Ce degré de liberté doit être fixé pour chaque point de passage Vi sur le chemin d'un point de passage V; avant de pouvoir déterminer le point de passage V;min de coût minimal. Le bord associé à ce point de passage V;mi" est noté 6ri". Lors d'une première étape, tous les points de passage V; associés aux bords 6; de la zone 5 de courant élémentaire contenant le point de départ sont évalués en termes de coût. Contrairement à un algorithme de propagation d'onde classique, l'évaluation du coût d'un point de passage V; comprend la détermination de sa position sur son bord 6; permettant de minimiser le coût. Lors d'une deuxième étape, on considère le point de passage V;'i" pour lequel le coût est minimal et on évalue le coût entre le point de départ et les points de passage V; associés aux bords 6; de la zone 5 de courant élémentaire adjacente au bord 6;min. De même, l'évaluation d'un coût pour un point de passage V; comprend la détermination de la position de ce point de passage V; ainsi que de tous les points de passage Vi sur le chemin entre le point de départ et le point de passage ces positions minimisant le coût du point de passage V. Le même processus est répété jusqu'à ce que le chemin atteigne un bord 6; de la zone 5 de courant élémentaire contenant le point d'arrivée : on considère le point de passage V;m'" pour lequel le coût est minimal et on évalue le coût entre le point de départ et les points cle passage V; associés aux bords 6; de la zone 5 de courant élémentaire adjacente au bord 6;m'" Dans le cas où plusieurs points de passage ont un coût égal, on considère successivement ces points de passage. La séquence de points de passage Vr" est enregistrée au fur et à mesure de la propagation de l'algorithme afin d'obtenir le chemin optimal entre le point de départ et le point d'arrivée. La suite de la description est consacrée à la détermination des points de passage Vi par un algorithme de propagation d'onde adapté. On 30 considère les notations suivantes : d(dx;dy) représente un déplacement du mobile dans le repère R(O, x, y), c(cX;cy) represente le vecteur vitesse du courant dans une zone 5 de courant élémentaire, v(vx;vy) représente le vecteur vitesse du mobile par rapport au courant. Etant donné que les zones 5 de courant élémentaire sont une décomposition exacte de l'espace en cellules, un chemin P entre le point de départ et le point d'arrivée passe nécessairement par une suite de zones 5 de courant élémentaire adjacentes. Par ailleurs, le courant étant uniforme dans chaque zone 5 de courant élémentaire, le mobile doit se déplacer en ligne droite entre deux points de passage Vi pour minimiser le temps de parcours. Par conséquent, un chemin P peut être défini comme une succession de segments de droite entre le point de départ et le point d'arrivée. On associe une fonction de coût t'•'' à chaque segment de droite entre deux points de passage Vi et voisins. Cette fonction de coût t'''' correspond au temps de parcours effectif du mobile entre les points de passage Vi et Vi+1. De manière générale, le temps de parcours t pour effectuer un mouvement d par un mobile se déplaçant à la vitesse v dans un courant c" vérifie la relation : d = (-5 +-d ).t . Projetée sur les axes x et y , cette relation devient : dx = (vx + cx ).t dy =(vy+cy).t La relation v2 =vx +vÿ permet d'éliminer vx et vy des équations précédentes et d'aboutir à la relation : (dx -cx.t)2 +(dy -cy.t)2 =v2.t2 (3) La résolution de cette équation du second degré donne les solutions : ù (cx.dx + cy.dy) t = 2 2 v ùc ' (~) 25 où A=v2.(dx +dy)ù(cx.dy ùcy.dx)2. Le temps de parcours étant forcément positif, la solution est donnée par la relation : - (c.d + c.d )+~ t= x x y y (5) 2 2 v ùc Par conséquent, la fonction de coût entre les points de passage Vj et 30 V~+1 peut être définie par la relation : (2) ù(cxj dj +cjYdj)+-FA (i)2 ù(Cj)2 où d j (dx ; d') représente le déplacement du mobile entre les points de passage Vi et Vi+0 (d,, = ùx,; d = y)+, ù y . ), cj k; ci,) représente le vecteur vitesse du courant dans la zone 5 de courant élémentaire dans 5 laquelle se déplace le mobile, vj (vx ; vÿ) représente le vecteur vitesse du mobile par rapport au vecteur vitesse cj du courant dans cette même zone 2 2 2 de courant élémentaire et où A=wj) .((di)2 +(dÿ) ù(cX.dÿ ùcy.d') . Le temps de parcours optimal pour un chemin entre un point de passage Vo, par exemple le point de départ, et un point de passage V;, par 10 exemple le point d'arrivée, peut être déterminé par un coût T. défini par la relation : 7 ;=. (7) d;,dv j=0 Cette relation signifie que le chemin entre Vo et V; est déterminé en minimisant de façon globale les fonctions de coût tj'j+' du coût T,, , la position 15 des points de passage Vi sur les bords 6i étant déterminée afin de minimiser le temps de parcours entre Vo et V. Dans une forme particulière de réalisation, les variables xj et yj pour chaque point de passage Vi sont substituées par une variable unique, l'abscisse curviligne I. Cette forme de réalisation permet de ne faire 20 dépendre le coût T. que des abscisses curvilignes Ii et donc de diviser par deux le nombre de variables du coût T. Le coût T, peut être déterminé par une technique d'optimisation numérique appelée méthode du gradient projeté. Cette technique est similaire à une descente du gradient, mais permet en plus de gérer des contraintes de 25 bornes sur les variables. En l'occurrence, les abscisses curvilignes li sont limitées à leur domaine Dj 4;1;1 La technique du gradient projeté converge généralement vers la solution optimale, toutefois sans l'atteindre. Pour stopper cette convergence, il est possible d'appliquer un paramètre de tolérance e sur l'annulation du gradient de T,.. Ce paramètre de tolérance e 30 est prédéterminé ; il peut être choisi par un utilisateur. Il a un impact direct sur la qualité du chemin planifié. Plus ce paramètre de tolérance s est (6) 12 proche de zéro, plus le temps de parcours du chemin planifié se rapproche de l'optimal. L'utilisation du paramètre de tolérance s peut être mathématisée par la relation suivante : I grad (T. l La troisième étape du procédé selon l'invention est illustrée par la figure 4. Sur cette figure, un espace 1 est représenté à différentes étapes, non successives, de l'algorithme de propagation d'onde selon l'invention. Les différentes étapes sont présentées chronologiquement de gauche à droite et de haut en bas. La répartition du courant dans l'espace 1 est identique à celle de la figure 1 a avec une vitesse de courant c =150 km/h . La vitesse du mobile est fixée à v =100 km/h . Les segments de droite 41 reliant les différents points de passage Vi représentent différents chemins P possibles entre le point de départ Vo et un point de passage V. A chaque étape de l'algorithme de propagation d'onde, un front d'onde 42 se propage avec l'algorithme de propagation d'onde selon l'invention. Ce front d'onde 42 délimite les zones 5 de courant élémentaire dont au moins un bord 6j a été exploré des autres zones 5 de courant élémentaire. Comme représenté à la figure 4, la position des points de passage Vi sur leur bord 6j, de même que la séquence de points de passage Vb autrement dit le chemin, peut évoluer 25 au cours de la détermination du coût T,. L'évolution de la séquence de points de passage V; et de leur position illustre la recherche d'un minimum global pour le coût T,,. La dernière étape représentée à la figure 4 correspond à l'étape où un chemin P est trouvé entre le point de départ Vo et le point d'arrivée VQ_1. 30 Le procédé selon l'invention ne s'applique pas seulement à un espace 1 bidimensionnel ; il s'applique également à un espace tridimensionnel. Dans un espace tridimensionnel, les zones 5 de courant élémentaire sont remplacées par des polyèdres. Chaque face du polyèdre 35 correspond à un bord 6i tel qu'il est défini dans un espace à deux REVENDICATIONS1. Procédé de planification de chemins (P) pour un mobile entre un point de départ (Vo) et un point d'arrivé (VQ_1) situés dans un espace (1) donné en présence de courants, ces courants pouvant varier en vitesse et en direction dans l'espace (1), le mobile étant soumis aux courants et se déplaçant à une vitesse (v) donnée par rapport aux courants, le procédé étant caractérisé en ce qu'il comprend les étapes suivantes : - discrétiser l'espace (1) en zones (5) de courant élémentaire, chaque zone (5) étant délimitée par des bords (6i), - associer à chaque bord (6i) un point de passage (Vi) possible 10 pour le mobile, chaque point de passage (Vi) pouvant translater sur son bord (6i), - déterminer une séquence de points de passage (Vi) permettant au mobile de relier le point de départ (Vo) au point d'arrivée (VQ_1) en un minimum de temps. 2. Procédé de planification de chemins selon la revendication 1, caractérisé en ce qu'une zone (5) de courant élémentaire est une zone de forme polygonale, indivisible et dans laquelle le courant est considéré comme homogène sur la totalité de la surface de la zone. 20 3. Procédé de planification de chemins selon la revendication 2, caractérisé en ce que les zones (5) de courant élémentaire sont construites en calculant le diagramme de Voronoï autour de points (3) pour lesquels on connait une vitesse et une direction de courant. 25 4. Procédé de planification de chemins selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination de la séquence de points de passage (Vi) permettant au mobile de relier le point de départ (Vo) au point d'arrivée (VQ_1) en un minimum de temps utilise un 30 algorithme de propagation d'onde adapté aux points de passage (Vi). 5. Procédé de planification de chemins selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un chemin (P ) comprend une succession de segments de droite (41) entre le point de 15départ (Vo) et le point d'arrivée (Vo..1), une fonction de coût étant associée à chaque segment de droite (41) entre deux points de passage Vi et voisins et définie par la relation : l ! ' ' t>>+1 =ù cX.dj 2cy.dy2+,1L~ où d'(dx;dÿ) représente le (v J ù (c ) déplacement du mobile entre les points de passage voisins V~ et Vi+1, c' (cx ; cÿ) représente le vecteur vitesse du courant dans la zone (5) de courant élémentaire dans laquelle se déplace le mobile, v' (vX ; vÿ) représente le vecteur vitesse du mobile par rapport au courant dans cette même zone (5) de courant élémentaire et où 0 = (v ` )2 .�d' )2 + (cl ÿ )2 )ù (ci .d ÿ cÿ .d y )2 10 6. Procédé de planification de chemins selon la revendication 5, caractérisé en ce que la détermination de la séquence de points de passage (Vi) permettant au mobile de relier le point de départ (Vo) au point d'arrivée (Vol) en un minimum de temps comprend la détermination de coûts T, pour 15 des points de passage V;, ces coûts T, étant définis par la relation : T. = min E 7. Procédé de planification de chemins selon la revendication 6, caractérisé en ce que les coûts T. sont déterminés en utilisant une méthode 20 du gradient projeté paramétrée par un paramètre de tolérance c.