Brevet FR3013867A1 — publié le 20150529

L'invention concerne un procédé de simulation de propagation d'ondes, notamment électromagnétiques ou encore acoustiques. Il est de plus en plus nécessaire de recourir à la simulation, que ce soit pour le développement de nouveaux produits, ou encore pour la qualification ou l'homologation de ceux-ci. En effet dans le cas où un produit interagit avec des ondes, il peut être nécessaire de simuler la propagation de ces ondes et leur interaction avec le produit. Les ondes peuvent dans certains cas être produites par le produit lui-même : c'est le cas par exemple pour des capteurs à ultrasons, des radars ou des lidars. Ce cas survient donc notamment lors du développement de radars ou lidars embarqués à bord de véhicules. De tels capteurs sont de plus en plus souvent prévus à bord de véhicules pour assurer la détection d'obstacles et ainsi contribuer à la réalisation de fonctions de sécurité active, d'assistance à la conduite voire de conduite automatique. Etat de l'art Pour ce qui est de la simulation de radars, les procédés de simulation existants sont fondés habituellement, pour les applications temps réel ou semi temps-réel, sur le lancer d'un petit nombre de rayons, dans un nombre de directions limité et un angle solide restreint, les objets présents dans la scène étant pris en compte de manière simplifiée sous la forme de leurs Surfaces Equivalentes Radar (SER). Par ailleurs, des modèles statistiques issus de campagne de test sont parfois également utilisés pour la simulation de radars. Or, si l'utilisation de tables de SER est souvent acceptable pour les radars aériens, pour lesquels le nombre d'objets à modéliser est très restreint, inversement cette méthode est inapplicable pour des radars embarqués à bord de véhicules terrestres, dont l'environnement proche comporte une très grande variété d'objets possibles, et de surcroit mobiles les uns par rapport aux autres. En effet dans de tels environnements, la simulation de propagation d'ondes radar sur la base de tables de SER, s'avère en pratique insuffisante, car elle ne permet de gérer que de manière relativement grossière le masquage des éléments et le déphasage induit par les rebonds multiples avant que l'onde ne revienne à l'antenne. Le procédé ne permet donc pas de simuler les ondes reçues en retour de manière suffisamment proche de la réalité. Présentation de l'invention Par suite, il existe un besoin d'un procédé de simulation de propagation d'ondes, notamment d'ondes électromagnétiques (en particulier, d'ondes hors du spectre visible, et/ou d'ondes non cohérentes) ou acoustiques, qui soit susceptible de fournir une représentation de la propagation d'ondes simulées de manière plus réaliste que les simulations réalisées en utilisant des SER. Cet objectif est obtenu grâce au fait que le procédé comporte les étapes suivantes : a) on fournit des données de scène représentant une scène tridimensionnelle comportant au moins un objet présentant une pluralité de surfaces, chacune desdites surfaces présentant une normale ; b) on calcule une pluralité de rayons primaires émis dans des directions de propagation respectives ; chacun desdits rayons primaires étant défini au moins par une puissance transportée, une position, et une direction de propagation ; d) en fonction d'un point de réception, on calcule des rayons diffusés primaires, chacun desdits rayons diffusés primaires étant émis par une surface dudit au moins un objet de la scène atteinte par un rayon primaire. Un rayon diffusé par une surface atteinte par un rayon, dit rayon incident, désigne ici un rayon émis par la surface atteinte sous l'effet du rayon incident, ayant une direction de propagation passant par le point de réception, et ayant une puissance transportée au moins fonction d'une orientation relative entre la direction de propagation du rayon incident et la normale à la surface atteinte. Le terme 'calcul d'un rayon' désigne ici le calcul des attributs ou paramètres caractéristiques de ce rayon. La « position » d'un rayon peut être notamment toute information 25 définissant le point d'émission de ce rayon. Elle peut être obtenue par l'intermédiaire des attributs du point d'émission du rayon. Le procédé de simulation permet donc de déterminer les différents rayons reçus par un observateur placé en un point donné par rapport à la scène, dit point de réception. 30 De manière connue en soi, le procédé fait appel à la technique du « lancer de rayons » : l'onde émise est considérée comme étant la réunion d'un ensemble de rayons émis, appelés ici « rayons primaires ». Une caractéristique importante du procédé est le fait que lors de l'étape d), on calcule des rayons diffusés par les surfaces des objets présents dans la scène. 35 Pour ces rayons, la puissance transportée est au moins fonction de l'orientation relative entre la direction de propagation du rayon incident et la normale à la surface atteinte (Cette orientation relative est appelée par la suite « angle d'incidence » du rayon sur la surface atteinte). Ainsi, l'étape d) permet de simuler le rayonnement perçu en retour par l'observateur placé au point de réception en prenant en compte la géométrie des objets présents dans la scène, cette géométrie étant prise en compte par le fait que pour chaque rayon diffusé, la puissance transportée par le rayon diffusé dépend de l'angle d'incidence du rayon primaire sur la surface atteinte. Le rayonnement perçu en retour au point de réception est alors la superposition ou la somme des rayons diffusés par les différentes surfaces des objets de la scène en direction de ce point. Modélisation des objets En général, les objets présents dans la scène sont enregistrés sous forme d'un ensemble de surfaces planes : la normale à la surface de l'objet atteinte par le rayon incident est alors simplement la normale de la surface plane atteinte par le rayon incident. Les zones planes peuvent être des facettes planes, notamment triangulaires. D'autres types de modélisation de la géométrie des objets présents dans la scène peuvent notamment être envisagés, notamment par des surfaces non planes (par exemple définies par des B-splines, des courbes de Bézier, etc.). Modélisation des rayons ; interaction rayons/surfaces La qualité de la simulation réalisée dépend en outre du modèle retenu pour l'interaction entre le rayon incident et la surface atteinte par celui-ci. Le modèle prédéterminé d'interaction entre un rayon incident et une surface d'un objet définit les propriétés du ou des rayons éventuellement émis par la surface, en fonction des propriétés d'une part du rayon incident qui atteint la surface, et d'autre part de la surface. Le ou les rayons émis par la surface sont (exclusivement) soit des rayons diffusés tels que définis précédemment, soit des rayons dits rayons « ré-émis ». Le calcul du ou des rayons émis comporte le calcul de différents paramètres. Ces paramètres peuvent inclure notamment, pour le rayon émis : sa direction de propagation ; son point d'émission (c'est le point de la surface atteint par le rayon incident) ; l'énergie ou la puissance transportée par le rayon émis ; la polarisation du rayon émis ; la phase du signal porté par le rayon émis ; la distance parcourue par le rayon, depuis l'émission du rayon primaire d'origine ; la fréquence du signal porté par le rayon émis. Ils peuvent également comporter, éventuellement, des caractéristiques de forme transverse du rayon, c'est-à-dire des caractéristiques de la forme d'une section du rayon dans un plan transverse à sa direction de propagation. Les valeurs de chacun des paramètres du rayon peuvent varier en fonction d'une part, du trajet entre le point d'émission du rayon incident, et le point atteint sur la surface atteinte (notamment éventuellement en fonction de paramètres de l'atmosphère entre ces deux points) ; et d'autre part en fonction de propriétés de la surface atteinte par le faisceau incident. Notamment, la puissance transportée par un rayon ré-émis peut être au moins fonction d'une orientation relative entre le rayon incident et la normale à la surface atteinte par celui-ci. Par ailleurs, les propriétés des rayons diffusés ou ré-émis peuvent dépendre, non seulement des positions ou orientations relatives du rayon incident par rapport à la surface atteinte de l'objet, mais aussi de propriétés de l'objet lui-même. Dans un mode de mise en oeuvre, on utilise une propriété, notamment physique d'un matériau constituant une surface d'un objet de la scène pour calculer les rayons diffusés. La propriété considérée ici peut être un coefficient de réflexion ou de transmission du matériau de la surface atteinte de l'objet. Par exemple, la propriété de la surface peut notamment être définie par des coefficients de réflexion, de transmission, de diffusion, du rayon incident, soit par les perméabilités électrique et magnétique du matériau considéré, qui peuvent être de simples nombres réels. La propriété de la surface peut également être définie de manière bien plus complexe, en étant par exemple fonction de la fonction de distribution de réflectance bidirectionnelle, ou BDRF (de l'Anglais « Bidirectional reflectance distribution function »), du matériau de la surface considérée. Dans ce mode de mise en oeuvre, il suffit donc d'attribuer certains paramètres, certaines propriétés à des matériaux, et de définir ensuite que des surfaces d'objets présents dans la scène sont constituées par le ou les matériaux ainsi définis. Pour le calcul d'un rayon diffusé, on peut prendre en compte un modèle prédéterminé d'interaction entre un rayon incident et une surface atteinte par celui- ci. Ainsi par exemple, dans un mode de mise en oeuvre, et dans le cas où les ondes sont électromagnétiques, à l'étape d) on calcule la puissance ou l'énergie transportée par le(s) rayon(s) diffusé(s) en utilisant l'équation intégrale du champ électrique (de l'Anglais : Electric field integral equation (EFIE)). Avantageusement, cette équation permet de calculer les propriétés des rayons diffusés de manière relativement simple. Les étapes a), b) et d) du procédé de simulation décrites précédemment permettent de simuler de manière performante la propagation de rayons dans une scène tridimensionnelle. Cependant, le procédé de simulation peut être rendu encore plus réaliste en prenant en compte les réflexions multiples des rayons frappant la scène. Dans ce but, dans un mode de mise en oeuvre, le procédé de simulation comporte en outre une étape complémentaire c) durant laquelle : c) on simule, notamment de manière itérative, la propagation de rayons incidents, de telle manière qu'à chaque itération, pour chaque rayon incident considéré : . on détermine si le rayon incident, par propagation en ligne droite, atteint une surface d'un objet de la scène ; et . si en fonction d'un modèle prédéterminé d'interaction entre un rayon incident et une surface d'un objet, au moins un rayon ré-émis (rayon émis par une surface de la scène, autre qu'un rayon diffusé) est émis après l'interaction dudit rayon incident avec ladite surface atteinte : . on détermine le point atteint sur la surface atteinte ; et, . on calcule ledit au moins un rayon ré-émis depuis la surface atteinte ; lesdits rayons incidents considérés étant : . lors de la première itération, les rayons primaires calculés à l'étape b) ; . lors de chacune des itérations ultérieures, les rayons ré-émis calculés à l'itération précédente. En outre à l'étape d), on calcule de plus les rayons diffusés émis par les surfaces du ou des objet(s) de la scène atteintes par un rayon ré-émis. Dans ce mode de mise en oeuvre, le procédé peut comporter une ou plusieurs itérations à l'étape c). Par ailleurs, comme pour les rayons diffusés on peut utiliser une propriété, notamment physique d'un matériau constituant une surface d'un objet de la scène pour calculer les rayons ré-émis. Les rayons ré-émis tels que définis à l'étape c) peuvent correspondre à une onde transmise ou à une onde réfléchie par l'objet dont la surface est atteinte par le rayon incident. D'autres types de rayons ré-émis peuvent être pris en considération, en fonction de la nature des ondes dont la propagation est simulée. A l'étape c), la détermination du point atteint sur la surface atteinte peut consister à déterminer un point, défini par des coordonnées tridimensionnelles, ou encore par des coordonnées bidimensionnelles sur la surface atteinte. Mais cela peut consister aussi de manière plus large à déterminer la surface atteinte, ou une partie de la surface atteinte, dès lors que cette opération est suffisamment précise pour permettre le calcul des caractéristiques des rayons émis suite à l'interaction du rayon incident avec l'élément atteint (point, portion de surface, ou surface atteint(e)). Ainsi dans ce mode de mise en oeuvre du procédé selon l'invention, la technique du 'lancer de rayons' est mise en oeuvre de manière beaucoup plus réaliste que dans les modes de mise en oeuvre antérieurs, dans lesquels l'interaction entre un objet et un rayon incident était modélisée de manière très simplifiée au moyen de SER. En effet à l'étape c), le calcul des rayons ré-émis permet d'obtenir une représentation beaucoup plus précise de la propagation des ondes. Les rayons ré-émis sont généralement calculés en prenant en compte la géométrie du ou des objets présents dans la scène. Ainsi par exemple, une ou plusieurs propriétés du ou des rayons ré-émis peuvent être fonction de la normale de la surface atteinte. La normale de la surface atteinte de l'objet est calculée naturellement au point de la surface atteint par le rayon incident. Dans le cas où les ondes sont des ondes électromagnétiques, le modèle d'interaction entre le rayon incident et la surface peut être fondé sur les équations de l'optique physique. Dans ce cas, dans un mode de mise en oeuvre, la puissance transportée par un rayon ré-émis est calculée par la loi de Snell-Descartes. Pour calculer la direction de propagation du ou des rayons ré-émis, dans un mode de mise en oeuvre, la direction de propagation d'au moins un rayon émis - et en général, de chacun des rayons ré-émis - dépend uniquement de la direction de propagation du rayon incident et de la normale de la surface atteinte. Par suite à l'étape c), et dans le cas où les rayons sont de nature électromagnétique, pour chaque rayon incident, on calcule seulement zéro ou un rayon réfléchi par la surface (réflexion spéculaire) et éventuellement un rayon transmis à l'intérieur de la surface, si la surface est transparente. Ensuite à l'étape d), on calcule s'il y a lieu le rayon diffusé par la surface sous l'effet du rayon incident. Aussi, le procédé selon l'invention avantageusement est un procédé à complexité sensiblement constante ; c'est-à-dire que le nombre de rayons généralement est relativement constant tout au long des itérations. En effet les rayons incidents qui ne rencontrent aucune surface ne donnent naissance à aucun rayon émis ; inversement, les surfaces transparentes (qui ne forment généralement qu'une très faible proportion des surfaces) peuvent donner naissance à deux rayons ré-émis pour un rayon incident. Pour déterminer la fréquence d'un rayon émis calculé à l'étape c), dans un mode de mise en oeuvre, on prend en compte une vitesse relative de la surface atteinte par le rayon incident, par rapport au point de réception (II s'agit d'une vitesse relative car le point de réception peut éventuellement être lui-même animé d'une vitesse propre). Structure de données ; moyens de calcul : GPU En plus des paramètres physiques qui gouvernent la propagation des rayons étudiés, le choix de certains moyens de calcul spécifiques permet de mettre en oeuvre le procédé de manière particulièrement efficace. En particulier, du fait notamment du choix de l'algorithme de lancer de rayons pour la propagation de l'onde émise, le procédé peut avantageusement être mis en oeuvre sur une carte graphique. Ainsi dans un mode de mise en oeuvre particulièrement intéressant de l'invention, l'étape c) et/ou l'étape d) est réalisée sur une (unique) carte graphique (CGU), et les calculs réalisés pour chacun des rayons sont réalisés en parallèle. Ce mode de mise en oeuvre permet donc de réaliser les calculs de l'étape c) et/ou d) de manière particulièrement rapide. Grâce à cela, avantageusement le procédé peut être mis en oeuvre de manière temps réel ou semi temps-réel. Pour mettre en oeuvre le procédé sur une carte graphique comme indiqué précédemment, les données de scène sont de préférence constituées principalement par au moins une matrice multidimensionnelle, ou bien par un ensemble de données individuelles tridimensionnelles. Dans le premier cas, la géométrie des objets n'est pas enregistrée sous forme de coordonnées tridimensionnelles. La position des objets au contraire est enregistrée sous la forme des coordonnées (numéro de ligne et numéro de colonne) d'une donnée dans une matrice, ainsi que sous la forme d'une valeur de profondeur nommée habituellement 'z'. La géométrie est donc stockée dans une ou plusieurs matrices appelée habituellement « Z-buffer ». Cette ou ces matrices contenant les informations de géométrie des objets de la scène font partie des données de scène, enregistrées dans une matrice multi-dimensionnelle, appelée encore 'G-buffer' (buffer ou tampon graphique). Par matrice « multi-dimensionnelle », on désigne ici un ensemble de matrices ayant le même nombre de lignes et de colonnes. Avantageusement, une matrice multi-dimensionnelle permet de stocker, pour une même position (i,j) dans la matrice multi-dimensionnelle, non pas seulement une information scalaire (comme par exemple une information de profondeur en 'z'), mais un nombre quelconque d'informations relatives à l'élément référencé à l'emplacement (i,j). Dans le second cas au contraire, la géométrie est enregistrée sous forme de données tridimensionnelles. Dans ce cas, les données de scène comportent principalement des entités qui contiennent les coordonnées tridimensionnelles des différentes surfaces des objets contenus dans la scène. Elles peuvent comporter en outre d'autres données, enregistrées par exemple dans des matrices multidimensionnelles, notamment les différents coefficients caractérisant les surfaces des objets de la scène. En outre, de même que la scène peut être enregistrée dans une matrice multi-dimensionnelle, avantageusement les rayons (primaires, et/ou diffusés, et/ou ré-émis) peuvent l'être également. Ainsi dans un mode de mise en oeuvre, les rayons primaire, les rayons ré-émis, et/ou les rayons diffusés, sont enregistrés dans une matrice multidimensionnelle. De préférence à l'étape c) on calcule la direction de propagation des rayons ré-émis au moyen des données de scène et des caractéristiques des rayons incidents, sans faire intervenir d'autres données représentant la géométrie de la 15 scène. Par ailleurs, dans un mode de mise en oeuvre du procédé selon l'invention, les calculs requis à l'étape c) et/ou d) peuvent être programmés sur carte graphique multiprocesseurs au moyen de « shaders ». Un shader est un programme, écrit dans un langage - soit assembleur, soit 20 un langage de plus haut niveau -, directement exécutable par une carte graphique (en Anglais Graphics processing Unit, « GPU ») et qui remplace certaines parties du pipeline d'exécution habituel. Il est possible notamment d'utiliser des shaders programmés en langage GLSL (« openGL Shading Langage »). L'utilisation de « shaders » nécessite impérativement d'enregistrer ou 25 modéliser les données ou informations sous forme de matrices multidimensionnelles, ce qui permet alors de les traiter par des « shaders » de manière massivement parallèle, dans un espace dit « espace image ». Ce mode de programmation est donc particulièrement adapté au cas où la scène et/ou les rayons (primaires, ré-émis, et /ou diffusés) sont représentés par des 30 matrices multidimensionnelles (c'est-à-dire par un G-buffer tel que présenté précédemment). Dans un autre mode de mise en oeuvre, les calculs requis à l'étape c) et/ou d) peuvent être programmés sur carte graphique multiprocesseurs en utilisant un langage de programmation de carte graphique multiprocesseurs autorisant un 35 accès direct aux instructions et à la mémoire des différents processeurs de calcul parallèle de la carte graphique, comme par exemple le langage de programmation Cuda (marque déposée). Ce mode de programmation peut notamment être choisi lorsque la scène est représentée - sous forme de données tridimensionnelles. Simulation de capteurs Une application particulièrement importante du procédé de simulation de propagation d'ondes selon l'invention est la simulation de capteurs. L'invention concerne ainsi également un procédé de simulation de capteur, notamment d'un capteur pour véhicule automobile, le capteur étant prévu pour émettre des ondes et pour produire un signal de sortie en fonction d'ondes reçues en retour suite à ladite émission d'ondes, le procédé comprenant les étapes suivantes : i) on simule l'émission d'ondes par le capteur, et la propagation desdites ondes, de manière à déterminer des rayons reçus par le capteur, en mettant en oeuvre un procédé de simulation de propagation d'ondes tel que défini précédemment ; et ii) on détermine le signal de sortie du capteur en fonction des rayons reçus par le capteur et de caractéristiques prédéterminées du capteur. La simulation de capteur se fait donc d'abord en simulant l'émission des ondes par le capteur, puis leur propagation (étape i). Lors de la simulation de propagation d'ondes à cette étape i), on simule la propagation des ondes de manière à calculer les rayons diffusés que va recevoir le capteur placé au point de réception. On calcule alors, en fonction des caractéristiques du capteur, le signal de sortie du capteur. Ce calcul peut être réalisé de différentes manières. Dans un mode de mise en oeuvre, le signal de sortie du capteur est déterminé à l'étape ii) en réalisant les deux étapes e) et f) suivantes : e) on calcule les signaux respectifs d'au moins une partie desdits rayons diffusés ; et f) on calcule un signal reçu au point de réception en additionnant les signaux respectifs calculés à l'étape e2). Le signal de sortie du capteur est alors déterminé en fonction des signaux des différents rayons diffusés retournés par la scène et reçus par le capteur au point de réception. En particulier, pour limiter les temps de calcul, il est souhaitable de ne calculer les signaux transportés par les rayons diffusés que pour un nombre réduit de rayons diffusés. Dans ce but, de préférence l'étape e) comporte les deux sous-étapes suivantes : el) on sélectionne un sous-ensemble de rayons diffusés en fonction d'un critère prédéterminé, notamment un critère prenant en compte la puissance transportée des rayons diffusé, et/ou au moins une caractéristique de rayons voisins telle qu'une distance totale parcourue et/ou une phase du signal transporté ; e2) on calcule les signaux respectifs des rayons diffusés sélectionnés seulement pour le sous-ensemble sélectionné de rayons diffusés. Typiquement, à l'étape e), la sélection du sous-ensemble de rayons diffusés est faite en triant les rayons diffusés par ordre de puissance ou d'énergie transportée. Cependant, les signaux respectifs des différents rayons diffusés peuvent éventuellement être calculés pour tous les rayons diffusés. L'addition des signaux respectifs des différents rayons diffusés prend en compte les caractéristiques respectives des rayons diffusés de manière à produire un signal cumulé représentatif de la réalité. Différents paramètres peuvent être pris en compte pour augmenter la qualité de la simulation. Dans un mode de mise en oeuvre, pour chaque rayon ré-émis calculé à l'étape c), on détermine une distance parcourue par le rayon depuis le point d'émission du rayon primaire ayant généré le rayon ré-émis considéré. Cela permet notamment à l'étape f), de calculer pour chaque rayon diffusé sélectionné la puissance du signal reçu par le capteur, en prenant en compte la distance totale effectivement parcourue par le rayon diffusé, et par conséquent la perte de puissance (exprimée par exemple en dB) qui s'ensuit. La qualité de la simulation peut en outre être améliorée en ajoutant un bruit aux signaux respectifs des différents rayons diffusés sélectionnés, et/ou au signal reçu ainsi calculé. La qualité de la simulation peut également être améliorée en prenant en compte de manière précise le fonctionnement de l'antenne d'émission d'ondes du capteur (Le terme 'antenne' est utilisé ici pour désigner tout système d'émission d'ondes que comprend le capteur). Dans un mode de mise en oeuvre, à l'étape b), on fournit un diagramme d'émission indiquant des pertes énergétiques en fonction de la direction de propagation ; et la puissance transportée par les rayons primaires est calculée en fonction du diagramme d'émission. Le diagramme d'émission permet ainsi de prendre en compte de manière réaliste les propriétés réelles d'émission de l'antenne du capteur. L'invention concerne également un procédé de simulation d'un véhicule automobile comportant au moins un capteur, dans lequel une propagation d'ondes émises et/ou reçues par le capteur est simulée en mettant en oeuvre un procédé de simulation de propagation d'ondes tel que défini précédemment, en particulier un procédé dans lequel le capteur est simulé par le procédé de simulation de capteur tel que défini précédemment. L'invention concerne également un programme d'ordinateur comportant des instructions pour l'exécution des étapes de l'un des procédés de simulation définis précédemment. L'invention concerne également un support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes de l'un des procédés de simulation définis précédemment. L'invention concerne enfin également un simulateur de propagation d'ondes, comprenant : a) une mémoire, apte à stocker des données de scène représentant une scène tridimensionnelle comportant au moins un objet présentant une pluralité de surfaces, chacune desdites surfaces présentant une normale ; b) des moyens de calcul primaire, aptes à calculer une pluralité de rayons primaires émis dans des directions de propagation respectives ; chacun desdits rayons primaires étant défini au moins par une puissance transportée, une position, et une direction de propagation ; d) des moyens de calcul de rayons diffusés, aptes à calculer en fonction d'un point de réception des rayons diffusés primaires, chacun desdits rayons diffusés primaires étant émis par une surface dudit au moins un objet de la scène atteinte par un rayon primaire ; un rayon diffusé par une surface atteinte par un rayon, dit rayon incident, étant un rayon réfléchi par la surface atteinte sous l'effet du rayon incident, ayant une direction de propagation passant par le point de réception, et ayant une puissance transportée au moins fonction d'une orientation relative entre le rayon incident et la normale à la surface atteinte. Dans un mode de réalisation, ce simulateur comprend en outre : c) des moyens de calcul de rayons ré-émis aptes, notamment de manière itérative, à chaque itération et pour chaque rayon incident considéré : à déterminer si le rayon incident, par propagation en ligne droite, atteint une surface d'un objet de la scène ; et . si en fonction d'un modèle prédéterminé d'interaction entre un rayon incident et une surface d'un objet, au moins un rayon ré-émis est émis après l'interaction dudit rayon incident avec ladite surface atteinte : . à déterminer le point atteint sur la surface atteinte ; et . à calculer ledit au moins un rayon ré-émis depuis la surface atteinte ; lesdits moyens de calcul de rayons ré-émis étant agencés de manière à prendre en considération en tant que rayons incidents : . lors de la première itération, les rayons primaires calculés à l'étape b) ; . lors de chacune des itérations ultérieures, les rayons ré-émis calculés à l'itération 10 précédente. En outre, les moyens de calcul de rayons diffusés sont aptes, à l'étape d), à calculer des rayons diffusés émis par les surfaces dudit au moins un objet de la scène atteintes par un rayon ré-émis. 15 L'invention sera bien comprise et ses avantages apparaîtront mieux à la lecture de la description détaillée qui suit, de modes de réalisation représentés à titre d'exemples non limitatifs. La description se réfère aux dessins annexés, sur lesquels : - la figure 1 est une vue schématique d'une route sur laquelle évoluent trois 20 véhicules ; - la figure 2 est une vue schématique représentant une structure de données permettant d'enregistrer la géométrie d'une scène ; - la figure 3 est une vue schématique représentant l'interaction d'un rayon incident avec une surface d'un objet ; 25 - la figure 4 est une vue schématique du dessus représentant la propagation de rayons primaires projetés par le radar dans la scène de la figure 1 ; - la figure 5 est une vue schématique du dessus représentant la propagation de rayons ré-émis émis par les véhicules présents dans la scène de la figure 1 après réception des rayons primaires ; 30 - la figure 6 est une vue schématique du dessus représentant la propagation de rayons diffusés par les véhicules présents dans la scène de la figure 1 après réception des rayons primaires ; - la figure 7 est une vue schématique d'un signal transporté par un rayon de radar ; - la figure 8 est une vue schématique d'un signal reçu par le radar ; et 35 - la figure 9 est une représentation schématique du procédé de simulation de propagation d'onde selon l'invention. Un mode de réalisation du procédé et du dispositif selon l'invention va maintenant être décrit en faisant référence aux figures. Dans le mode de réalisation présenté, le procédé de propagation d'ondes est un procédé de propagation d'ondes électromagnétiques, utilisé pour simuler le 5 fonctionnement d'un radar embarqué sur un véhicule automobile, par exemple une voiture. Un tel radar est mis en oeuvre typiquement dans une scène telle que celle représentée sur la figure 1. Sur cette figure sont représentées trois voitures A, B, C circulant sur une 10 route 10. La voiture C est équipée d'un radar 12. Celui-ci émet depuis un point P (plus précisément, depuis une zone considérée comme sensiblement ponctuelle) un rayonnement électromagnétique. La partie de ce rayonnement qui est exploitée est celle qui se propage dans l'angle solide S représenté sur la figure 1. Cet angle solide S se divise lui-même en 15 une matrice d'angles solides élémentaires Su où i désigne la ligne et varie de 1 à n, et j désigne la colonne et varie de 1 à p. La partie de l'espace incluse dans l'angle solide S, et limitée du côté avant de la voiture C par un plan L situé à une distance prédéterminée de la voiture C (En pratique, à environ 200 mètres de celle-ci) constitue une scène 14. La fonction du 20 radar 12 est de détecter des obstacles à l'intérieur de cette scène. Dans un premier temps, les grandes étapes du procédé de simulation de propagation de rayons selon l'invention vont être présentées en relation avec la figure 9. Cette figure présente la simulation de radar comme un processus 25 périphérique d'un processus principal qui est la simulation de fonctionnement d'un véhicule automobile. La circulation de la voiture C sur la route 10 en effet est simulée par un simulateur de véhicule automobile 100. Ce simulateur 100 comporte un calculateur central 102 qui exécute un 30 programme informatique 104 de simulation de véhicule, dit 'moteur de simulation'. Ce programme 104 simule la circulation d'un véhicule, et notamment dans le cas présent, la circulation ou le déplacement de la voiture C. Pour cela, le programme 104 simule l'acquisition des différentes informations acquises par les capteurs de la voiture C. Il simule donc notamment l'acquisition des 35 informations fournies par le radar 12. Cette dernière simulation (simulation de radar) est réalisée par un calculateur de radar 110. Celui-ci comprend une unité centrale 112 avec un processeur principal et une mémoire de stockage, ainsi qu'une carte graphique 114. Le programme 104 fonctionne de manière itérative, avec une boucle de traitement en temps réel ou en semi-temps réel simulant l'évolution de la circulation à intervalle régulier, par exemple de 40 ms (en temps réel). A chaque boucle de traitement, le programme 104 transmet au calculateur de radar 110 la description géométrique CG', Fig.9) de la scène qui se trouve en face du radar 12 à l'instant considéré. Cette description inclut les vitesses relatives par rapport au radar 12 des différents objets présents dans la scène. Dans la carte graphique est enregistré un programme de traitement noté également 114. La carte graphique exécute le programme 114 de manière à effectuer les traitements suivants : a) A réception des données relatives à la scène 14, le programme 114 constitue des données de scène représentant la scène placée en face du radar 12 ; b) le programme 114 réalise alors un 'lancer de rayons', et calcule ainsi les caractéristiques des différents rayons primaires qui sont émis par le radar ; c) de manière itérative, le programme 114 simule la propagation des rayons incidents dans la scène 14. Ces rayons incidents sont : . lors de la première itération, les rayons primaires émis par le radar ; et . lors de chacune des itérations ultérieures, les rayons ré-émis calculés à l'itération précédente. Aussi, à chaque itération et pour chacun des rayons incidents, le programme 114 détermine si le rayon incident, par propagation en ligne droite, atteint une surface d'un objet de la scène. Si c'est le cas, le programme 114 détermine alors si un ou plusieurs rayons sont ré-émis après l'interaction du rayon incident considéré avec la surface atteinte. L'émission d'un ou de rayons après cette interaction, et les caractéristiques du ou des rayons ré-émis, sont déterminées en fonction du modèle d'interaction entre un rayon incident et une surface d'un objet, qui est choisi par avance. Le programme 114 détermine donc d'une part, le point atteint sur la surface atteinte ; et, d'autre part, les caractéristiques du ou des rayons ré-émis depuis la surface atteinte sous l'effet du rayon incident considéré. d) le programme 114 calcule les rayons diffusés par les différentes surfaces présentes dans la scène 14 en direction du point de réception, c'est-à-dire ici en direction du radar 12 ; el) le programme 114 sélectionne un sous-ensemble de rayons diffusés en fonction d'un critère prédéterminé prenant en compte la puissance ou l'énergie transportée des rayons diffusé ; il transmet alors cette information à l'unité centrale 112 du calculateur 110. Dans l'unité centrale est enregistré un programme de traitement noté également 112. Le programme 112 effectue les traitements suivants : e2) le programme 112 calcule des signaux respectivement pour chacun des rayons diffusés sélectionnés à l'étape el). f) le programme 112 calcule un signal reçu par le radar en additionnant les signaux respectifs calculés à l'étape e2). Il calcule alors le signal de sortie Sr du radar en fonction du signal reçu par le radar. a) Données de scène Pour simuler le fonctionnement du radar 12, on utilise ou on crée tout d'abord des données de scène représentant la scène 14. Ces données comportent au minimum un ensemble de données géométriques incluant en principe tout objet présent dans la scène, à condition que celui-ci ait une taille suffisante pour être détectable dans un angle solide élémentaire Si. Pour chaque objet, les données de scène comportent donc, pour ses différentes surfaces susceptibles d'être illuminées par un rayon émis par le radar, une description simplifiée de la surface, comportant notamment les coordonnées du vecteur normal à la surface. Les données de scène comportent également de préférence différents coefficients ou paramètres caractérisant la surface (ou du matériau constituant la surface) et servant au calcul des rayons ré-émis suite à l'interaction d'un rayon incident avec la surface. Ces coefficients peuvent être notamment des coefficients caractérisant les propriétés de la surface relatives à la réflexion, la transmission et la diffusion de rayons incidents. Ces coefficients peuvent être par exemple des indices de réfraction, une épaisseur, une rugosité... En outre, les données de scène peuvent inclure, pour chaque surface, l'indication de la vitesse instantanée de la surface par rapport au radar. Ces différents coefficients sont en général enregistrés par une référence au matériau constituant la surface. Dans ce cas, on enregistre pour chaque matériau constituant une surface des objets présents dans la scène 14, des coefficients caractérisant la réponse du matériau lorsqu'une surface constitué par ce matériau reçoit un rayon incident, ce rayon incident ayant la fréquence et la forme d'onde produites par le radar dont le fonctionnement est simulé. Les données de scène sont enregistrées dans la mémoire de la carte graphique 114. Elles peuvent être enregistrées de différentes manières. Dans le mode de réalisation présenté ici, elles sont enregistrées sous forme d'une matrice multidimensionnelle M (Fig.2), encore appelée en infographie « Gbuffer » ou « Graphics-buffer ». La matrice M comporte le même nombre de lignes et colonnes que la matrice des angles solides élémentaires ; mais elle est constituée d'un certain nombre (q) de matrices individuelles, notées M1, M2,... Mg. Ainsi, la localisation de chaque surface d'un objet de la scène est enregistrée seulement d'une part, via les composantes (i,j) de l'élément de matrice dans lequel est enregistrée la surface, et d'autre part par une information de profondeur 'z', qui est enregistrée dans la matrice M, par exemple dans la matrice M1. Il s'agit donc d'un Z- buffer. Les matrices suivantes M2, M3, ... Mu servent à enregistrer les autres propriétés de la surface (ces propriétés sont enregistrées dans cet exemple par u-1 nombres réels). Ces autres propriétés peuvent consister en l'indication du matériau constituant la surface, les propriétés de ce matériau étant enregistrées par ailleurs. Si dans un angle solide élémentaire donné, la scène comporte plusieurs surfaces successives à des distances croissantes du point P, ces différentes surfaces sont enregistrées dans la matrice M, dans des matrices individuelles Mk, où k commence par exemple à la valeur u+1. De manière alternative, les données de scène peuvent être enregistrées comme une scène tridimensionnelle comprenant des objets eux-mêmes composés d'un ensemble de facettes, notamment triangulaires, des propriétés de surface étant associées aux facettes sous forme notamment de textures, au moyen de coordonnées de texture de manière connue en soi. b) Calcul des rayons primaires La simulation du fonctionnement du radar nécessite d'abord de calculer les rayons émis par le radar depuis le point d'émission P (rayons primaires). Pour simplifier le calcul, on calcule un rayon primaire Pu pour chaque angle solide élémentaire S. Les rayons primaires sont donc enregistrés dans une matrice multidimensionnelle similaire à la matrice M décrite précédemment. La figure 4 représente de manière schématique les rayons primaires Pu émis sur une ligne horizontale d'angles solides Su. Pour chaque rayon primaire Pu sont enregistrées les informations suivantes : - sa direction de propagation (Celle-ci est enregistrée implicitement par les coordonnées (i,j) du rayon P1 considéré); - son point d'émission (pour les rayons primaires, c'est le point P) ; - la puissance transportée par le rayon émis. Cette information est enregistrée sous forme d'un gain - négatif - en dB par rapport à la puissance maximale initiale d'un rayon primaire ; - la polarisation du rayon primaire ; - la phase du signal porté par le rayon primaire ; - la distance parcourue par le rayon depuis l'émission du rayon primaire d'origine (distance nulle, dans le cas des rayons primaires) ; - la fréquence du signal porté par le rayon primaire. - Le coefficient doppler associé au rayon primaire (coefficient nul dans le cas des rayons primaires). La valeur de gain n'est pas la même pour tous les rayons primaires. En effet, pour chaque rayon primaire, le gain est fourni par le diagramme d'émission de l'antenne d'émission du radar 12. Ce diagramme d'antenne indique, en fonction de l'angle solide élémentaire Su, l'énergie ou la puissance effective - et donc le gain du rayon que le radar émet dans l'angle solide élémentaire. C'est donc cette valeur qui est enregistrée comme gain initial pour chaque rayon primaire. En pratique, le gain des rayons primaires est élevé au centre du faisceau émis par le radar, et plus faible sur les côtés. c) Propagation des rayons La propagation des rayons est fondée sur un modèle d'interaction avec la surface connu en soi (Fig.3). Conformément à ce modèle, en fonction de paramètres de réflexion, de transmission et de diffusion de la surface, lorsqu'un rayon incident RI frappe une surface en un point E, il peut être absorbé sans émission d'aucun rayonnement ; il peut donner naissance à un rayon réfléchi Rr ; il peut donner naissance à un rayon transmis Rt. Les caractéristiques de ces différents rayons dépendent chacune de la direction d'incidence (et donc de la direction de propagation du rayon incident) et de la normale n à la surface. Par ailleurs, un rayonnement peut être diffusé depuis le point E dans la direction de réception (qui peut être la direction d'émission des rayons primaires) ; il y a alors projection d'un rayon diffusé Rd en direction du point P. En outre, dans ce mode de mise en oeuvre l'interaction du rayon incident avec la surface atteinte repose sur les hypothèses suivantes: On considère que le rayon incident se comporte comme une onde plane monochromatique, et cela même si le signal transporté par le rayon incident (par exemple, dans des radars à modulations FMCW, FSK, à impulsions, ...) est en réalité la superposition de plusieurs ondes monochromatiques. On considère de plus que lorsqu'un front d'onde électromagnétique interagit à la surface d'un objet (ou à l'interface entre deux milieux), on considère localement cette surface comme plane et on l'assimile à un plan infini pour le calcul du rayon émis. (Cependant, le calcul peut éventuellement prendre en compte la courbure de la surface). Un exemple d'une interaction des rayons primaires avec les objets de la scène est illustré par les figures 4 à 6. Dans un des plans d'émission de rayons primaires Pu (j=Cte), illustré par la figure 4, les rayons primaires frappent le véhicule A dans trois angles solides élémentaires et le véhicule B également dans trois angles solides élémentaires (flèches en gras, Fig.4). Pour chaque rayon primaire Pu, lors de la première itération du programme 114 à l'étape c), on évalue l'interaction des différents rayons primaires Pu avec les surfaces des objets présents dans la scène. On notera que l'interaction avec le sol constituant la route 10 doit normalement être prise en compte ; cependant pour simplifier les explications cette interaction n'est pas discutée dans l'exemple présenté. Sur le véhicule A, les surfaces atteintes sont des surfaces opaques situées à l'arrière du véhicule ; aucun rayon n'est transmis à travers la surface. En revanche, ces surfaces donnent naissance à trois rayons réfléchis ré-émis RrA- Sur le véhicule B, les surfaces atteintes sont les surfaces transparentes du pare-brise avant du véhicule. Ces surfaces donnent tout d'abord naissance à trois rayons réfléchis ré-émis Rrg. Etant transparentes, elles donnent également naissance à trois rayons transmis ré-émis RtI3. L'évaluation de l'interaction entre les rayons primaires (ou, aux étapes ultérieures de calcul, les rayons incidents) et les surfaces des objets de la scène conduit donc à créer de nouveaux rayons, les rayons ré-émis RrA, Rr13/ RtB. Ces rayons n'ont pas les mêmes caractéristiques que les rayons primaires. Leurs caractéristiques sont calculées de la manière suivante : - la direction de propagation est calculée suivant le modèle d'interaction avec la surface considérée, en fonction de la direction du rayon incident et de la direction normale de la surface au point frappé par le rayon incident ; - le point d'émission du rayon ré-émis est le point de la surface frappé par le rayon incident ; - le gain, la polarisation du rayon ré-émis sont calculés suivant le modèle d'interaction avec la surface considérée, en fonction des caractéristiques de la surface, et éventuellement en fonction de la direction du rayon incident et de la direction normale de la surface au point frappé par le rayon incident ; - la phase du signal porté par le rayon ré-émis, et la distance parcourue par le rayon depuis l'émission du rayon primaire d'origine, sont calculées en fonction de la position du point frappé sur la surface, par rapport au radar, et éventuellement en fonction de propriétés de la surface ; et - la fréquence du signal porté par le rayon ré-émis est calculée en fonction de la vitesse relative de la surface par rapport au radar. Un paramètre important est le gain du rayon ré-émis, qui correspond à la puissance du signal transporté par le rayon. Le programme 114 calcule le gain des rayons ré-émis en utilisant l'équation de l'optique physique : La puissance transportée par un rayon ré-émis est calculée par la loi de Snell-Descartes. Naturellement, un grand nombre de rayons primaires ne rencontrent aucune surface d'un objet de la scène 14 ; ils ne donnent alors naissance à aucun rayon émis. A l'étape c), après la première interaction de chacun des rayons primaires avec les surfaces des objets de la scène 14 a été évaluée, le programme 114 peut réaliser une ou plusieurs itérations de manière à calculer une ou plusieurs interactions ultérieures entre les rayons ré-émis et les surfaces des objets de la scène. Dans l'exemple présenté cependant, le programme 114 est paramétré de manière à passer à l'étape d) dès que la deuxième itération de l'étape c) a été réalisée, c'est-à-dire dès que tous les rayons primaires ont atteint les surfaces qu'ils peuvent atteindre par propagation en ligne droite depuis le point P. De manière plus générale, le programme 114 peut être paramétré de manière à passer à l'étape d) soit après un nombre fixe d'itérations, soit lorsqu'un autre critère a été atteint. Eventuellement, l'étape c) peut ne pas être réalisée. Le programme 114 passe directement de l'étape b) à l'étape d). Les calculs sont réalisés de manière massivement parallèle sur la carte graphique 114. Le nombre de rayons primaires est choisi de manière à engendrer un nombre de rayons pouvant être traité par la carte graphique à la fréquence de traitement souhaitée. d) Calcul des rayons diffusés A l'étape d), le programme 114 calcule les rayons diffusés par les différentes surfaces présentes dans la scène 14 en direction du radar 12. Ce calcul est fait pour chacune des surfaces de la scène 14. Une même surface peut éventuellement donner naissance à plusieurs rayons émis (diffusés ou ré-émis) si elle est atteinte par plusieurs rayons incidents, notamment lors de plusieurs itérations effectuées à l'étape c) de propagation. Dans l'exemple présenté, l'étape c) est arrêtée à la fin de la première itération. Le programme 114 calcule donc six rayons diffusés : Trois rayons RdA diffusés par la paroi arrière de la voiture A, et trois rayons RdB diffusés par la paroi avant ou la paroi latérale de la voiture B. Le programme 114 ne calcule pas de rayon Rdg pour les autres surfaces de la voiture B, car à la fin de la première itération (Fig.5), seules la paroi avant et une paroi latérale de la voiture B ont été atteintes par des rayons incidents (Puisqu'il s'agit de la première itération, les rayons incidents sont les rayons primaires Pu). Pour chaque surface de la scène 14, le programme 114 calcule le ou les rayon(s) diffusé(s) (éventuel(s)) émis par la surface, en prenant en compte l'ensemble des rayons incidents ayant atteint la surface considérée au cours des différentes itérations de l'étape c). Les rayons diffusés comportent sensiblement les mêmes paramètres que les rayons primaires. el) Tri et sélection des rayons diffusés Ensuite à l'étape el), le programme 114 sélectionne un sous-ensemble de rayons diffusés en fonction d'un critère prédéterminé prenant en compte la puissance transportée par les rayons diffusés. Ce critère est fondé plus précisément sur la puissance retournée à l'antenne, qui est fonction de la puissance transportée par les rayons diffusés. La puissance retournée à l'antenne par un rayon diffusé est calculée sur la base de caractéristiques de l'antenne de réception du radar 12, notamment les dimensions ou la surface de celle-ci. Elle est calculée éventuellement en fonction du gain de l'antenne, ce gain pouvant dépendre de la direction de provenance du rayon diffusé. Plus concrètement, cette puissance est calculée sur la base de l'équation intégrale du champ électrique, en intégrant le flux du vecteur de Poynting du rayon diffusé à travers la surface de l'antenne. Le calcul pour chaque rayon diffusé de la puissance retournée à l'antenne par le rayon diffusé permet de procéder à la sélection des rayons diffusés qui seront pris en compte pour le calcul du signal reçu effectivement par le radar. En effet, cette sélection consiste en pratique à ne prendre en compte dans la suite des calculs que les rayons qui vont retourner une puissance non négligeable au radar. Dans ce but, on commence par trier les rayons diffusés en fonction de la puissance retournée à l'antenne par chacun ; ensuite, on ne conserve que les rayons diffusés de plus grande puissance retournée à l'antenne. Ce travail de sélection permet de réduire la quantité de données à transmettre à l'unité centrale 112 du calculateur de radar 110, et cela sensiblement sans perte de performance pour la simulation, dans la mesure où les rayons diffusés éliminés sont des rayons dont la contribution au signal radar est très faible et ainsi négligeable. Selon une variante, toujours dans le but de limiter la quantité d'information remontée par le simulateur de canal de propagation RADAR accéléré par le GPU, on sélectionne les rayons diffusés de la manière suivante : On forme d'abord des groupes de rayons présentant des similitudes en termes de position, de distance totale parcourue par le rayon, de déphasage, et/ou de fréquence. Un tel regroupement des rayons peut être réalisé par des algorithmes de partitionnement (en Anglais : « clustering ») directement sur la carte graphique et de manière parallèle, par le programme 114. Chaque groupe de rayons diffusés est alors caractérisé par les valeurs statistiques caractéristiques du groupe : moyenne et variance de la distance totale parcourue, moyenne et variance du déphasage, moyenne et variance de la direction de réception, puissance totale. Dans ce cas, pour sélectionner une partie des rayons diffusés, on sélectionne les groupes de rayons de plus grande importance. Les rayons isolés, ou les groupes de faible importance sont assimilés à du bruit radar (en Anglais : clutter) et ainsi éliminés. On peut aussi réaliser la sélection en prenant en compte, éventuellement en plus des critères cités précédemment, la distance totale parcourue. Par exemple, on peut ne sélectionner que les N rayons (N est un entier) ayant la plus faible distance totale parcourue avant le retour au point de réception. Les informations correspondants aux groupes de rayons sélectionnés peuvent alors être transmises au programme 112 du calculateur de radar 110 de manière synthétique, en ne transmettant que des valeurs agrégées correspondants aux groupes de rayons sélectionnés. e2) Calculs des signaux des rayons diffusés Les données relatives aux rayons diffusés ainsi sélectionnés sont alors transmises à l'unité centrale 112 du calculateur de radar 110. A réception de ces informations (étape e2)), le programme 112 calcule un signal Sd pour chacun des rayons diffusés sélectionnés. Ce signal est calculé sur la base du signal initialement émis par le radar, tel que représenté par la figure 7. Les signaux respectifs Sd ainsi générés pour les différents rayons diffusés sélectionnés sont calculés en faisant l'hypothèse que les rayons primaires sont émis pendant un même laps de temps entre des instants 0 et T. La durée T est choisie de manière à être grande par rapport aux déphasages induits par les différences de distance de vol entre les différents rayons (les échelles sur les figures 7 et 8 ne sont pas représentatives). Bien que la figure 7 représente un signal sinusoïdal à fréquence constante, on comprend que le procédé peut être utilisé avec tout type de signal radar, notamment des signaux de fréquences variables, continûment ou par paliers par exemple (FSK, FMCW, ...). Les signaux individuels correspondant à chacun des rayons diffusés ont une forme voisine du signal radar initial, par exemple dans le cas présent, une forme sinusoïdale. Cependant, ils sont modifiés par rapport au signal initial en prenant notamment en compte les paramètres suivants du rayon diffusé : la phase du rayon ; l'énergie ou la puissance transportée par le rayon (ou la puissance retournée à l'antenne), la polarisation du rayon, le cas échéant, la distance parcourue par le rayon depuis l'émission du rayon primaire ; la fréquence du signal transporté par le rayon. f) Calcul du signal reçu par le radar et du signal de sortie du radar Ensuite à l'étape f), le programme 112 calcule le signal Sr reçu par le radar en additionnant les signaux respectifs Sd des différents rayons diffusés sélectionnés, calculés à l'étape e2). Ce calcul prend en compte des paramètres de réception propres au radar lui-même, c'est-à-dire des paramètres d'antenne de réception. Sur la base du signal Sr reçu par le radar, le programme 112 calcule alors le signal de sortie du radar Ss. Le signal de sortie du radar Ss ainsi obtenu est alors transmis au programme 104 de simulation de déplacement de véhicule (Fig.9). Ce signal Ss peut alors faire l'objet de traitements complémentaires, puis être fusionné avec d'autres informations acquises par d'autres capteurs du véhicule pour alimenter, par exemple, un programme déterminant automatiquement les commandes à appliquer au véhicule pour assurer le pilotage ou la conduite de celui-ci. Enfin, bien que l'exemple présenté, un radar ayant une unique antenne pour l'émission et la réception ait été utilisé, l'invention peut être mise en oeuvre avec une pluralité de sources de rayons. Elle peut de plus être mise en oeuvre en évaluant les rayons diffusés (étape d)) non pas pour une seule position de réception, mais pour plusieurs positions de réception. Dans l'exemple présenté, la carte graphique 114 constitue à la fois les moyens de calcul primaire, les moyens de calcul de rayons diffusés, et les moyens de calcul de rayons ré-émis, au sens de l'invention. REVENDICATIONS1. Procédé de simulation de propagation d'ondes, notamment d'ondes électromagnétiques ou acoustiques, dans lequel : a) on fournit des données de scène représentant une scène tridimensionnelle (14) comportant au moins un objet (A,B) présentant une pluralité de surfaces, chacune desdites surfaces présentant une normale ; b) on calcule une pluralité de rayons primaires (P11) émis dans des directions de propagation respectives ; chacun desdits rayons primaires étant défini au moins par une puissance transportée, une position, et une direction de propagation ; le procédé se caractérisant en ce qu'il comporte en outre une étape d) dans laquelle : d) en fonction d'un point de réception (P), on calcule des rayons diffusés primaires (RdA,RdB), chacun desdits rayons diffusés primaires étant émis par une surface dudit au moins un objet de la scène atteinte par un rayon primaire ; un rayon diffusé (Rd) par une surface atteinte par un rayon, dit rayon incident, étant un rayon émis par la surface atteinte sous l'effet du rayon incident, ayant une direction de propagation passant par le point de réception, et ayant une puissance transportée au moins fonction d'une orientation relative entre le rayon incident et la normale (n) à la surface atteinte. 2. Procédé de simulation de propagation selon la revendication 1, dans lequel les données de scène sont constituées principalement par au moins une matrice multidimensionnelle (M), ou bien par un ensemble de données individuelles tridimensionnelles. 3. Procédé de simulation de propagation selon la revendication 1 ou 2, dans le cas où les ondes sont électromagnétiques, dans lequel à l'étape d), on calcule la puissance transportée par ledit au moins un rayon diffusé en utilisant l'équation intégrale du champ électrique. 4. Procédé de simulation de propagation selon l'une quelconque des revendications 1 à 3, comportant en outre une étape complémentaire c) durant laquelle : c) on simule, notamment de manière itérative, la propagation de rayons incidents, de telle manière qu'à chaque itération, pour chaque rayon incident considéré : . on détermine si le rayon incident, par propagation en ligne droite, atteint une surface d'un objet de la scène ; et. si en fonction d'un modèle prédéterminé d'interaction entre un rayon incident et une surface d'un objet, au moins un rayon ré-émis (RrA,RrgiRtg) est émis après l'interaction dudit rayon incident avec ladite surface atteinte : . on détermine le point atteint sur la surface atteinte ; et, . on calcule ledit au moins un rayon ré-émis depuis la surface atteinte ; lesdits rayons incidents considérés étant : . lors de la première itération, les rayons primaires calculés à l'étape b) ; . lors de chacune des itérations ultérieures, les rayons ré-émis calculés à l'itération précédente ; procédé dans lequel : à l'étape d), on calcule en outre des rayons diffusés émis par les surfaces dudit au moins un objet de la scène atteintes par un rayon ré-émis. 5. Procédé de simulation de propagation selon la revendication 4, dans lequel au moins une propriété d'au moins un rayon ré-émis est fonction de la normale (n) de la surface atteinte. 6. Procédé de simulation de propagation selon la revendication 4 ou 5, dans lequel la direction de propagation de chacun des rayons ré-émis dépend uniquement de la direction de propagation du rayon incident et de la normale (n) de la surface atteinte. 7. Procédé de simulation de propagation selon l'une quelconque des revendications 4 à 6, dans le cas où les ondes sont électromagnétiques, dans lequel la puissance transportée par un rayon ré-émis est calculée par la loi de Snell-Descartes. 8. Procédé de simulation de propagation selon l'une quelconque des revendications 1 à 7, dans lequel l'étape c) et/ou l'étape d) est réalisée sur une unique carte graphique (114), et les calculs réalisés pour chacun des rayons sont réalisés en parallèle. 9. Procédé de simulation de propagation selon l'une quelconque des revendications 1 à 8, dans lequel les rayons primaires, les rayons ré-émis, et/ou les rayons diffusés, sont enregistrés dans une matrice multidimensionnelle. 10. Procédé de simulation de propagation selon l'une quelconque des revendications 1 à 9, dans lequel on utilise une propriété d'un matériau constituant une surface d'un objet de la scène pour calculer les rayons ré-émis et/ou les rayons diffusés. 11. Procédé de simulation d'un capteur (12), notamment d'un capteur pour véhicule automobile, le capteur étant prévu pour émettre des ondes et pour produire un signal de sortie en fonction d'ondes reçues en retour suite à l'émission des ondes 5 émises, le procédé comprenant les étapes suivantes : i) on simule l'émission d'ondes par le capteur, la propagation desdites ondes, de manière à déterminer des rayons reçus par le capteur, en mettant en oeuvre un procédé selon l'une quelconque des revendications 1 à 10 ; et ii) on détermine le signal de sortie (S5) du capteur en fonction des rayons (Rd) reçus 10 par le capteur et de caractéristiques prédéterminées du capteur. 12. Procédé de simulation d'un capteur selon la revendication 11, dans lequel, pour déterminer le signal de sortie du capteur à l'étape ii), on réalise des étapes el), e2) et f) durant lesquelles : 15 el) on sélectionne un sous-ensemble de rayons diffusés (Rd) en fonction d'un critère prédéterminé, notamment un critère prenant en compte la puissance transportée des rayons diffusé, et/ou au moins une caractéristique de rayons voisins telle qu'une distance totale parcourue et/ou une phase du signal transporté ; e2) on calcule les signaux respectifs des rayons diffusés sélectionnés seulement 20 pour le sous-ensemble sélectionné de rayons diffusés ; et f) on calcule un signal reçu au point de réception en additionnant les signaux respectifs calculés à l'étape e2). 13. Procédé de simulation d'un capteur selon la revendication 11 ou 12, dans lequel 25 à l'étape b), on fournit un diagramme d'émission indiquant des pertes énergétiques en fonction de la direction de propagation ; et la puissance transportée par les rayons primaires est calculée en fonction du diagramme d'émission. 14. Programme d'ordinateur comportant des instructions de code de programme 30 pour l'exécution des étapes du procédé de simulation selon l'une quelconque des revendications 1 à 13 lorsque ledit programme est exécuté sur un calculateur. 15. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes 35 d'un procédé de simulation selon l'une quelconque des revendications 1 à 13. 16. Simulateur de propagation d'ondes, comprenant :a) une mémoire, apte à stocker des données de scène représentant une scène tridimensionnelle comportant au moins un objet présentant une pluralité de surfaces, chacune desdites surfaces présentant une normale ; b) des moyens de calcul primaire (114), aptes à calculer une pluralité de rayons primaires émis dans des directions de propagation respectives ; chacun desdits rayons primaires étant défini au moins par une puissance transportée, une position, et une direction de propagation ; d) des moyens de calcul de rayons diffusés (114), aptes à calculer en fonction d'un point de réception des rayons diffusés primaires, chacun desdits rayons diffusés primaires étant émis par une surface dudit au moins un objet de la scène atteinte par un rayon primaire ; un rayon diffusé par une surface atteinte par un rayon, dit rayon incident, étant un rayon réfléchi par la surface atteinte sous l'effet du rayon incident, ayant une direction de propagation passant par le point de réception, et ayant une puissance transportée au moins fonction d'une orientation relative entre le rayon incident et la normale (n) à la surface atteinte. 17. Simulateur selon la revendication 16, comprenant en outre c) des moyens de calcul de rayons ré-émis (114) aptes, notamment de manière itérative, à chaque itération et pour chaque rayon incident considéré : à déterminer si le rayon incident, par propagation en ligne droite, atteint une surface d'un objet de la scène ; et . si en fonction d'un modèle prédéterminé d'interaction entre un rayon incident et une surface d'un objet, au moins un rayon ré-émis est émis après l'interaction dudit rayon incident avec ladite surface atteinte : . à déterminer le point atteint sur la surface atteinte ; et, . à calculer ledit au moins un rayon ré-émis depuis la surface atteinte ; lesdits moyens de calcul de rayons ré-émis étant agencés de manière à prendre en considération en tant que rayons incidents : . lors de la première itération, les rayons primaires calculés à l'étape b) ; . lors de chacune d'itérations ultérieures, les rayons ré-émis calculés à l'itération précédente ; les moyens de calcul de rayons diffusés étant en outre aptes, à l'étape d), à calculer des rayons diffusés émis par les surfaces dudit au moins un objet de la scène atteintes par un rayon ré-émis.