FR 3003425 A1 20140919 FR 1352242 A 20130313 - Domaine technique La présente invention s'applique avantageusement à la phase de caractérisation et d'établissement des réseaux de téléphonie mobile à installer sur un territoire et plus particulièrement elle concerne un procédé de prédiction du niveau de champs radio, électromagnétique ou sonore utilisant un modèle de propagation basé sur la technique de lancer de rayons en exploitant des bases de données terrain à haute résolution dans des environnements tels que, par exemple, les environnements urbains à forte densité. Ainsi, la mise en oeuvre de l'invention permettra de déterminer, aisément et efficacement, des paramètres de fonctionnement et des positions optimales des stations ou répéteurs du réseau de manière à définir précisément les limites des cellules radio correspondantes et ainsi optimiser l'utilisation des ressources radio allouées au service mobile. Technique antérieure Les modèles de propagation à lancer de rayons sont largement connus et très couramment utilisés dans une foultitude d'applications, notamment dans l'évaluation prédictive de la couverture radio d'une cellule en phase de définition de réseaux. Ce sont des modèles déterministes, ces modèles reposant sur une connaissance précise de la réalité environnementale et nécessitant de disposer de bases de données géographiques de type contours de bâtiments. Ils permettent de prédire les différents trajets de propagation dans une configuration donnée. Après ajustement de la bande de fréquence considérée, ces modèles permettent également par simulation de réaliser des études paramétriques pour analyser, par exemple, l'influence des diagrammes d'antennes ou des caractéristiques des SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 2 - matériaux, ce qui est très significativement moins onéreux que le fait de lancer un grand nombre de campagnes de mesures. Pour rappel, la propagation radio dans une cellule doit satisfaire aux deux exigences essentielles que sont l'émission d'une puissance non excessive par la station et la réception par les terminaux de signaux radio de puissance suffisante, chacune de ces exigences répondant à des spécifications reposant sur un certain nombre de valeurs caractéristiques à ne pas dépasser. Par exemple, une augmentation de la puissance d'émission d'une station accroîtrait inévitablement la taille de la cellule en provoquant des interférences entre cellules voisines. Egalement, côté terminaux mobiles, la puissance maximale doit être limitée par des impératifs de sécurité des utilisateurs et d'autonomie de leurs batteries. Divers procédés de prédiction du niveau de champs radio, électromagnétique ou sonore, sont connus et utilisés parmi lesquels celui d'une grande efficacité qui consiste à effectuer un calcul prévisionnel des atténuations en une pluralité de points de la cellule en utilisant une base de données vectorielle fournie, par exemple, par l'Institut Géographique National (IGN), base vectorielle représentant la carte de la zone géographique analysée avec les bâtiments et autres reliefs du sursol qu'elle comporte. Pour estimer l'atténuation prévisionnelle en tout point de la cellule, il est effectué une simulation de propagation de signaux radioélectriques. La propagation est modélisée en lançant fictivement, au moyen d'un calculateur d'une chaîne de production de prévisions de couverture radio, un rayon électromagnétique depuis la station, dans des directions déterminées, puis ses conditions de propagation dans l'angle solide élémentaire qu'il occupe sont calculées. Excepté la propagation en espace libre lorsque SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 3 - la station est en vue directe d'un terminal radio fictif qui est une propagation d'atténuation linéaire connue, la trajectoire des rayons rencontrent des obstacles qui l'atténuent et la dévient de façon supplémentaire, en particulier dans les microcellules dont les stations sont souvent à une altitude inférieure à celle des toits des immeubles. Ainsi, dans une rue, un rayon peut être dévié par réflexion ou diffraction/diffusion et l'ouverture de son angle solide peut même s'en trouver accrue. La demande internationale WO-A-97/4497 divulgue un tel procédé. Cependant, cette technique bien qu'efficace présente des inconvénients majeurs du simple fait qu'elle nécessite la disponibilité de bases de données vectorielles pour procéder à la modélisation de tous les contours des bâtiments. Or, ces bases de données vectorielles sont onéreuses à produire et à maintenir, et sont en outre extrêmement lourdes à gérer et à exploiter impliquant des temps de calcul prohibitifs lors de leur utilisation pour le traitement d'une application de prédiction du niveau de champ radio dans les réseaux de téléphonie mobile. Exposé de l'invention La présente invention a pour but de proposer une solution efficiente aux problèmes et inconvénients 25 présentés par les diverses techniques de l'art antérieur et en particulier la technique exploitant des bases de données vectorielles. La présente invention vise ainsi à réduire très significativement la puissance de calcul nécessaire pour la prédiction de la couverture radio d'une cellule 30 dans un tel réseau mobile tout en présentant la même efficacité qu'avec les techniques exploitant directement des bases de données vectorielles et en permettant, lorsque cela est désiré, d'exploiter le procédé de prédiction du SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 4 - niveau de champ radio pour les réseaux de téléphonie mobile lors d'une analyse en temps réel. Pour cela, selon l'invention, le procédé de prédiction du niveau de champ radio d'un réseau de téléphonie mobile pour la couverture d'une cellule utilisant la technique de lancer de rayons et exploitant des bases de données géographiques est remarquable en ce que la base de données directement exploitée est une base matricielle composée de pixels, chaque pixel portant au moins deux informations, une information propre au type de sursol et une information propre à la hauteur du sursol, à la base matricielle est appliqué un traitement pour linéariser les contours d'obstacles rencontrés alors que des rayons sont fictivement lancés, les rayons selon les obstacles rencontrés étant soit, réfléchis soit, diffractés et diffusés tandis que l'atténuation globale est calculée. Ainsi, le procédé de prédiction du niveau de champ radio selon l'invention qui utilise un modèle de propagation à base de lancer de rayons utilisant avantageusement des données matricielles est extrêmement rapide, tout en étant précis, opérationnel et d'une importante facilité d'exploitation. Ce procédé est applicable à toutes les technologies radio, étant basé sur un concept tout à fait nouveau de modélisation de la propagation dans des environnements urbains denses avec des bases de données terrain à haute résolution. En résumé, les principaux avantages que ce nouveau procédé offre sont les suivants. La rapidité de calcul, caractéristique essentielle du fait qu'il utilise une technique basée sur un échantillonnage matriciel et non vectoriel des angles des façades et des coins des immeubles. Une grande précision car le modèle prend en compte tous les aspects et phénomènes de la propagation radio pour délivrer des résultats avec une extrême SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 5 - précision. Le modèle exploité est précis pour les antennes au-dessus des toits pour lesquelles le signal provient essentiellement de la diffraction verticale, mais également pour les antennes au-dessous des toits pour lesquelles le 5 signal provient majoritairement de la diffraction horizontale. Enfin le procédé est immédiatement opérationnel et facile à utiliser car il est compatible avec toutes les technologies mobiles telles que la 2G, la 3G, le WiMax, la 4G ou LIE, le DVB, etc. et l'ensemble des 10 fréquences utilisées. Le système mettant en oeuvre le procédé peut être fourni avec des paramètres pré-calibrés qui sont intuitifs et aisément accessibles et compréhensibles. Enfin ce système supporte favorablement une interface intuitive pour l'auto-calibrage. L'auto- 15 calibrage est une fonction qui permet de trouver d'une manière automatique les meilleurs paramètres d'un modèle de prédiction et ceci en utilisant des mesures réelles réalisées sur le terrain. Ainsi le modèle de prédiction contient une fonction qui permet d'analyser les points de 20 mesures et trouver les meilleurs paramètres du modèle qui feront que les mesures soient le plus proche possible des prédictions, la comparaison se faisant sur l'erreur moyenne et l'écart-type par exemple. Suivant un mode de réalisation du procédé de 25 prédiction du niveau de champ radio selon l'invention, avant le lancer fictif de rayons, un prétraitement de vectorisation des contours d'obstacles est appliqué aux bases matricielles pour linéariser les contours d'obstacles rencontrés en ne prenant en compte que les objets 30 représentant les obstacles et ainsi retrouver lesdits contours d'obstacles, le prétraitement de vectorisation étant adapté pour entraîner le lancer de rayons une fois que les informations vectorielles propres aux types, c'est-à-dire, intérieur, façade ou coin d'obstacle et propres aux SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 6 - angles, c'est-à-dire, angle de façade ou angles de coin, ce pour chaque pixel, ont été extraites, pour former une matrice de types et une matrice des angles et prédire en fonction de ces matrices la réflexion, les diffraction et diffusion pour chaque rayon lancé et permettre le calcul de l'atténuation globale. Suivant un second mode de réalisation du procédé selon l'invention, un rayon fictif est lancé directement sur la base matricielle dans une direction d'un angle solide qui couvre complètement l'espace, son trajet étant en chaque point scruté et chaque fois que ledit rayon rencontre un obstacle et qu'un point d'impact est détecté, une analyse de la base de données matricielle est effectuée en temps réel pour créer des objets vectoriels adaptés au processus de lancer de rayons, une étude locale de l'environnement du point d'impact étant réalisée pour déterminer d'une part, si le point d'impact est un élément d'une façade ou d'un coin et d'autre part, l'angle de ladite façade ou les deux angles du coin de deux façades qui construisent le point, une fois l'environnement proche déterminé, le rayon continue son trajet, soit réfléchi si le point d'impact est un élément d'une façade, soit diffracté et diffusé si le point d'impact est le coin formé par deux façades. Suivant une variante avantageuse du second mode de réalisation du procédé selon l'invention, lorsque le point d'impact est détecté au moment de lancer l'analyse de la base de données matricielle en temps réel et avant qu'une étude locale de l'environnement du point d'impact ne soit effectuée, il est vérifié que ce point d'impact n'a pas encore été analysé, si ce n'est pas le cas le rayon continue son trajet, si c'est le cas l'analyse du point d'impact est autorisée, puis une fois l'analyse exécutée l'environnement du point d'impact est mémorisé pour être SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 7 - réutilisé dans la suite du processus de lancer de rayons et le rayon continue son trajet. Le procédé de prédiction du niveau de champ radio est également remarquable en en ce que les rayons sont lancés 5 dans toutes les directions horizontales en utilisant un incrément réduit, chaque lancer de rayon s'effectuant de manière optimisée, c'est-à-dire qu'un rayon n'est lancé selon un angle donné que s'il est prévu que ce rayon aura un impact significatif sur la prédiction, une telle 10 prévision étant obtenue par apprentissage lorsque des lancers de rayons précédents selon des angles voisins ont eu un impact sur la prédiction. Ainsi, le fait de lancer des rayons dans toutes les directions horizontales avec un incrément réduit, c'est-à- 15 dire un pas suffisamment petit et de créer cet apprentissage a pour effet avantageux d'accélérer encore et très significativement les calculs. Description sommaire des figures 20 D'autres avantages et caractéristiques ressortiront mieux de la description qui va suivre de modes préférés de mises en oeuvre du procédé selon l'invention en référence aux figures annexées données à titre d'exemple non limitatif. 25 Les figures la et lb donnent de manière schématique un exemple d'une étape d'un premier mode de réalisation du procédé selon l'invention. Les figures 2a et 2b montrent chacune un exemple d'une seconde étape du premier mode de réalisation du 30 procédé selon l'invention. Les figures 3a et 3b représentent, pour un second mode de réalisation du procédé selon l'invention, la description d'une première étape de détection d'un point SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 8 - d'impact entre un rayon lancé et un obstacle rencontré et d'une seconde étape d'analyse dudit point d'impact. Les figures 4a et 4b présentent une troisième étape propre à la continuation du lancer de rayon du second mode 5 de réalisation du procédé selon l'invention. Meilleure manière de réaliser l'invention technique Pour une meilleure appréhension de l'objet de l'invention suivent ci-après quelques rappels relatifs aux 10 diffractions verticale et horizontale et aux atténuations dues auxdites diffractions. L'estimation et le calcul relatifs à l'atténuation Lv due à la diffraction verticale sur les toits sont basés sur la technique de l'extraction de radiale optimisée 15 consistant à extraire un profil vertical le long d'une radiale, ceci en mettant à jour uniquement les données qui ont changé par rapport à une radiale voisine précédemment extraite, puis à calculer les atténuations résultantes des diffractions/diffusions sur les obstacles détectés sur le 20 profil vertical et ainsi mettre à jour les atténuations calculées en chaque point parcouru par ladite radiale.. L'atténuation Lv résulte d'une diffraction de type hybride qui repose sur une combinaison du modèle dit de diffraction couteau ou « knife-edge diffraction » et du modèle dit de 25 Walfish-Ikegami, modèle valide pour des distances relativement courtes généralement à l'intérieur de microcellules. L'estimation et le calcul relatifs à l'atténuation L11 due à la diffraction horizontale sont basés sur une 30 technique, unique et ultrarapide dans la présente demande, de lancer de rayons ou « ray-tracing » qui repose d'une part sur un échantillonnage des angles des obstacles, obstacles qui, en milieu urbain, sont constitués pour la plupart d'immeubles, et d'autre part à lancer des rayons SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 9 - dans toutes les directions horizontales avec un pas suffisamment petit, ce lancer de rayon s'effectuant de manière optimisée c'est-à-dire qu'un rayon n'est lancé selon un angle donné que si on prévoit que ce rayon aura un impact significatif sur la prédiction, une telle prévision pouvant être obtenue par la détection après apprentissage lorsque des lancers de rayons précédents selon des angles voisins ont eu un impact sur la prédiction En effet, le fait de lancer des rayons dans toutes les 10 directions horizontales avec un pas suffisamment petit et de créer cet apprentissage présente l'important avantage d'accélérer très significativement les calculs. Pour cela, dans un exemple préféré de réalisation de prévision qu'un lancer de rayon selon un angle donné aura 15 un impact significatif sur la prédiction, ladite prévision résulte de la comparaison à un seuil du nombre de pixels pour lesquels la prédiction a été mise à jour par le rayon. L'estimation et le calcul propres à l'atténuation Lv due à la diffraction verticale n'intéressent pas la 20 présente demande, étant connus en soi, la présente invention concerne essentiellement la diffraction horizontale et l'estimation et le calcul de sa composante 141 d'atténuation. Le procédé de l'art antérieur habituellement utilisé 25 est ici tout d'abord succinctement évoqué, procédé qui nécessite la disponibilité et l'exploitation de bases de données vectorielles modélisant l'ensemble des contours de bâtiments. Comme précédemment précisé, ces bases de données sont volumineuses, de ce fait elles sont très lourdes à 30 gérer et à exploiter, enfin elles requièrent d'énormes investissements pour les constituer, les produire et les maintenir. Ces bases sont particulièrement lourdes à gérer dans le cadre d'un traitement de type lancer de rayons, car il est nécessaire de retrouver les contours qui se SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 10 - retrouvent sur le chemin des rayons parmi la multitude de contours disponibles dans la base de données, cette contrainte présentant des problèmes d'une grande complexité pour ce qui est de l'indexation des contours et par conséquent impliquant des traitements lourds relativement à la puissance de calcul nécessaire. Dans ce type de bases de données, chaque élément du sursol est défini par un contour qui implique les informations suivantes : - la liste des points constituant les vecteurs qui définissent le contour du bâtiment, - un type de sursol dans le contour : immeuble, forêt, etc., - la hauteur du sursol. De cette manière, la façade d'un immeuble est assimilée à une surface verticale ou respectivement un coin d'un immeuble est assimilé à deux parties de surfaces verticales dont les projections au sol correspondent à une droite ou respectivement à deux droites sécantes formant un angle et son sommet. Dans ces conditions, le rayon lancé est soit réfléchi lorsqu'il rencontre une façade donc une droite, soit diffracté et diffusé lorsqu'il rencontre un coin donc le sommet d'un angle et ceci se produisant et se calculant selon les lois de l'optique géométrique classique. Pour résumer, afin que les prédictions et les calculs soient justes il est donc essentiel que, d'une part, lorsqu'un rayon incident lancé frappe une façade, l'angle exact de la façade soit connu et d'autre part, lorsqu'un rayon incident lancé frappe le sommet de l'angle formé par deux façades, les deux angles des deux façades soient également exactement connus. Contrairement à cela, pour assurer efficacement la couverture d'une cellule au sein d'un réseau de téléphonie mobile ou cellulaire, le procédé de prédiction du niveau de SOFTW6-FR-1 TEXTE DEPOSE HG/ML champ radio selon l'invention utilise également la technique de lancer de rayons en exploitant des bases de données géographiques, mais de manière nouvelle, remarquable et avantageuse, les bases de données exploitées sont directement des bases matricielles et non plus des bases vectorielles. Une base matricielle est composée de pixels et dans cette base chaque pixel porte au moins deux informations, une information propre au type de sursol (bâtiments, arbres, ponts, etc.) et une information propre à la hauteur dudit sursol. De ce fait, les contours d'obstacles ne sont pas linéaires, comme avec l'emploi d'une base vectorielle, car ils sont formés de pixels, un traitement doit alors être appliqué à la base matricielle pour linéariser les contours d'obstacles et former ainsi des droites dont les positions et les différents angles intéressants sont connus avec précision. Deux modes de réalisation principaux du procédé de prédiction du niveau de champ radio selon l'invention seront ci-après décrits, un premier mode de réalisation avec un traitement préalable et un second mode de réalisation exploitant directement les bases de données matricielles en temps réel pendant le lancer de rayons. Les figures la et lb donnent de manière schématique un exemple d'une étape d'un premier mode de réalisation du procédé avec prétraitement selon l'invention pour lequel une configuration d'obstacles simple en milieu urbain est proposée, obstacles vers lesquels les rayons incidents fictifs seront lancés. Dans ce premier mode réalisation, deux étapes sont décrites pour mettre en oeuvre efficacement le procédé de prédiction du niveau de champ radio. Une première étape consiste en un prétraitement spécifique et unique de vectorisation des contours d'obstacles appliqué aux bases matricielles pour linéariser les contours d'obstacles SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 12 - rencontrés en ne prenant en compte que les objets représentant les obstacles et ainsi retrouver lesdits contours d'obstacles. Ce prétraitement est réalisé d'une manière unique car tout à fait différemment des méthodes classiques, puisqu'est acquise une connaissance du type de résultat, il est évident par exemple que les immeubles ont quasiment tous des façades verticales. De cette manière, connaissant la méthode de traitement et la précision de calcul du modèle lui-même, la méthode de prétraitement spécifique ne s'embarrasse pas de précisions inutiles. De même, ce prétraitement est basé sur l'hypothèse très commune que les bâtiments ont en grande majorité des façades planes avec projection en ligne droite, ce qui implique qu'il est aisé de retrouver l'angle d'une façade par un procédé de vectorisation par segments de droites, un tel procédé pouvant par exemple consister à trouver la ligne droite qui minimise la distance moyenne entre cette ligne droite et le centre des pixels détectés comme appartenant à cette façade, cette dernière détection étant effectuée d'une part en ne considérant que les pixels ayant un type de sursol correspondant à un bâtiment et adjacents à un pixel ne faisant pas partie d'un bâtiment, tout en excluant les pixels de façade qui dégradent la distance moyenne avec la meilleure ligne droite possible. Comme cela a été préalablement expliqué, la figure la représente de manière schématique, à titre d'exemple, trois immeubles, Il, 12 et 13 d'un élément d'une base matricielle correspondant à une matrice géographique composée de points ou de pixels. Chaque point ou pixel comporte deux informations, la première relative au type de sursol, bâtiments, arbres, ponts, etc., la seconde relative à la hauteur dudit sursol. Une telle base, si les contours d'obstacles n'étaient pas linéarisés, ne pourrait être efficacement exploitée pour le lancer de rayons puisque les SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 13 - angles des façades ne correspondant pas à la réalité ne présentent que des angles droits relatifs aux pixels. La figure lb montre l'élément de la base matricielle, avec les immeubles II, 12 et 13, de la figure la auquel le prétraitement de vectorisation spécifique a été appliqué afin de linéariser les contours d'obstacles. Ainsi et en outre, le prétraitement de vectorisation exploité est particulièrement bien adapté pour le lancer de rayons car l'algorithme de vectorisation utilisé pendant cette première étape est conçu pour prendre en compte également la future seconde étape de lancer de rayons. Les données de vectorisation ne seront pas conservées pour éviter d'avoir à traiter directement des vecteurs qui, comme cela a été précédemment expliqué, sont excessivement lourds à gérer ce qui a pour effet de considérablement ralentir le calcul. Néanmoins, à partir desdits vecteurs les deux informations vectorielles suivantes sont extraites pour chaque pixel, informations propres au type, type 1 : intérieur d'immeuble, type 2 : façade d'immeuble, type 3 : coin d'immeuble et informations propres aux angles : dans l'exemple figure lb, angles de façades des immeubles II, 12, 13 : 450 ou 135° ou les deux angles de coins d'immeuble : 135° et 45°, 45° et 135°. Ainsi sont obtenues une matrice des types et une matrice des angles qui vont permettre de prédire en fonction des informations qu'elles contiennent la réflexion, les diffraction et diffusion pour chaque rayon lancé et permettre le calcul de l'atténuation globale. Une fois ce traitement appliqué à la base, dans une seconde étape, les rayons sont fictivement lancés, par conséquent non sur des bases vectorielles mais sur des bases matricielles. Deux phénomènes sont alors constatés concernant lesdits rayons qui, selon les obstacles rencontrés, sont soit réfléchis soit diffractés et diffusés SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 14 - alors que leurs trajets pourront alors être formellement prédits tandis que l'atténuation globale sera précisément calculée. Les figures 2a et 2b montrent chacune, conservant la configuration d'environnement des figures la et lb, une seconde étape du premier mode de réalisation du procédé selon l'invention, seconde étape qui correspond donc au lancer de rayons. La figure 2a propose un exemple de réflexion d'un rayon incident Ri lancé à partir d'une antenne ANT en direction de l'immeuble II et qui rencontre ce dernier en un point A. Le fait d'avoir créé la matrice des types et la matrice des angles permet de connaître précisément les informations propres à ce point A. Ainsi, selon la matrice des types, ce point A est de type 2, il est situé sur une façade de l'immeuble II. Dans cet exemple simple, la façade forme un angle de 45° avec la direction du rayon incident, information connue de la matrice des angles. De cette manière, il est facile de calculer l'angle du rayon réfléchi Rr à partir de la connaissance de l'angle du rayon incident Ri, pris comme référence 0° et de l'angle de la façade en appliquant la formule traditionnelle de l'optique géométrique classique : ar = 2 (cf - ai), dans laquelle ar représente l'angle de réflexion, cf représente l'angle de la façade et ai représente l'angle du rayon incident. Comme ici l'angle ai = 0° et l'angle af = 45°, on obtient aisément ar = 90°. Ainsi sont effectués de proche en proche les calculs relatifs au phénomène de réflexion. La figure 2b propose un exemple de diffraction/diffusion d'un rayon incident Ri lancé à partir de l'antenne ANT en direction de l'immeuble 13 et qui rencontre ce dernier en un point B. La matrice des types et la matrice des angles permettent de connaître les informations propres à ce point B. Ainsi, selon la matrice SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 15 - des types, ce point B est de type 3, il correspond à un coin de l'immeuble 13. Dans cet exemple simple, le coin est formé par deux angles, un premier angle 131 = 135° et un second angle 132 = 45° relativement à la direction du rayon 5 incident Ri, informations connues de la matrice des angles. Le calcul des directions des rayons diffractés/diffusés Rd est réalisé en appliquant les lois connues de l'optique géométrique classique, ce calcul dépend donc ici des deux angles 131 et 132 qui définissent au point B le coin de 10 l'immeuble 13, ces deux angles étant précisément connus de la matrice des angles. Suit la description d'un second mode de réalisation du procédé selon l'invention se dispensant également de l'utilisation de lourdes bases vectorielles. Les figures 3a 15 et 3b donnent de manière schématique un exemple du second mode de réalisation exploitant directement des bases de données matricielles pendant la totalité du déroulement du processus de lancer de rayons, ceci pour mettre en oeuvre en temps réel le procédé selon l'invention. La configuration 20 d'obstacles simple en milieu urbain proposée est identique à celle déjà proposée avec le premier mode de réalisation, obstacles vers lesquels les rayons incidents fictifs sont lancés. Ainsi, les rayons fictifs sont lancés directement sur 25 la représentation de la base matricielle, chacun dans une direction d'un angle solide qui couvre complètement l'espace. L'algorithme du lancer de rayons sur la base matricielle ici utilisé consiste à scruter son trajet en chaque point et chaque fois que ledit rayon rencontre un 30 obstacle une analyse de la base de données matricielle est effectuée en temps réel pour créer des objets vectoriels adaptés au processus de lancer de rayons, une étude locale de l'environnement du point d'impact étant pour cela SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 16 - réalisée, puis le processus de lancer de rayons est poursuivi, ceci de proche en proche. De manière remarquable trois étapes sont parcourues durant le processus de lancer de rayons pendant lesquelles 5 les bases matricielles sont directement exploitées. Une première étape est ainsi constituée par la détection du point d'impact entre le rayon incident lancé et un obstacle rencontré, une seconde étape consiste en l'analyse du point d'impact pour créer des objets vectoriels adaptés au 10 processus de lancer de rayons alors qu'une troisième étape représente la continuation du lancer de rayon. Les figures 3a et 3b vont permettre de suivre le déroulement de la première étape de détection d'un point d'impact et de la seconde étape d'analyse dudit point 15 d'impact, la figure 3a concernant l'exemple d'un point d'impact sur la façade d'un immeuble alors que la figure 3b intéresse un point d'impact sur le coin d'un immeuble. Dans une première étape, chaque point du rayon incident lancé Ri à partir de l'antenne ANT est suivi pour 20 détecter si ce rayon rencontre un immeuble II, 12, 13, et lorsque c'est le cas, le point d'impact est ainsi détecté et déterminé. Deux cas se présentent selon que le point est situé sur la façade d'un immeuble II, point A dans l'exemple de la figure 3a ou sur le coin d'un immeuble 13, 25 point B, dans l'exemple de la figure 3b. Une fois le point d'impact détecté et déterminé lors de la première étape, commence la seconde étape dédiée à la création d'objets vectoriels adaptés au processus de lancer de rayons au cours de laquelle une analyse de l'environnement proche du 30 point A ou du point B est réalisée afin d'une part de spécifier si le point A ou respectivement le point B est un élément d'une façade ou respectivement d'un coin formé par deux façades et d'autre part de préciser l'angle de la façade ou bien les deux angles du coin sur lequel est SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 17 - localisé le point A ou bien le point B. Sur l'exemple de la figure 3a, le point A est un élément d'une façade de l'immeuble Tl dont l'angle avec le rayon incident Ri est, dans cet exemple, égal à 45°. Sur la figure 3b, le point B est un élément d'un coin formé par deux façades dont les angles avec le rayon incident Ri sont, dans cet exemple, respectivement de 135° et 45°. Un algorithme de vectorisation ou de linéarisation locale est pour cela mis en oeuvre, par exemple l'algorithme décrit précédemment utilisant l'adaptation de lignes droites aux pixels de façade, lequel algorithme est parfaitement adapté au lancer de rayons étant conçu pour prendre en compte la troisième étape subséquente de continuation de lancer de rayons. Sur les figures 4a et 4b est représentée la troisième étape propre à la continuation du lancer de rayon. Ainsi, l'environnement proche du point A ou du point B étant parfaitement et précisément déterminé lors du déroulement des première et seconde étapes, à la troisième étape le lancer du rayon peut être poursuivi selon les méthodes normales. Sur la figure 4a est proposé un exemple de réflexion d'un rayon incident Ri lancé à partir d'une antenne ANT en direction de l'immeuble Tl et qui rencontre ce dernier en un point A. La seconde étape a permis de connaître précisément les informations propres à ce point A. Ainsi, il a été déterminé que le point A est situé sur une façade de l'immeuble Tl et que, dans cet exemple simple, la façade forme un angle de 45° avec la direction du rayon incident. De la même manière qu'avec le premier mode de réalisation, il est facile de calculer l'angle du rayon réfléchi Rr à partir de la connaissance de l'angle du rayon incident Ri, pris comme référence 0° et de l'angle de la façade en appliquant la formule traditionnelle de l'optique SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 18 - géométrique classique : ar = 2 (cf - ai), dans laquelle ar représente l'angle de réflexion, cf représente l'angle de la façade et ai représente l'angle du rayon incident. Comme ici l'angle ai = 00 et l'angle af = 45°, on obtient aisément ar = 90°. Ainsi, le lancer de rayon peut être poursuivi et de proche en proche sont effectués les calculs d'atténuation relatifs au phénomène de réflexion. La figure 4b propose un exemple de diffraction/diffusion d'un rayon incident Ri lancé à partir de l'antenne ANT en direction de l'immeuble 13 et qui rencontre ce dernier en un point B. La seconde étape a permis de connaître précisément les informations propres à ce point B. Ainsi, il a été déterminé que le point B correspond à un coin de l'immeuble 13. Dans cet exemple simple, le coin est formé par deux angles, un premier angle 131 = 135° et un second angle 132 = 45° relativement à la direction du rayon incident Ri. Le calcul des directions des rayons diffractés/diffusés Rd est réalisé en appliquant les lois connues de l'optique géométrique classique, ce calcul dépend également ici des deux angles 131 et 132 qui définissent au point B le coin de l'immeuble 13, ces deux angles étant précisément connus, le lancer de rayon peut être poursuivi. Une variante avantageuse du second mode de réalisation du procédé selon l'invention est à présent décrite. Cette variante est en fait une extension directe du second mode de réalisation mais qui combine les avantages des deux modes de réalisation précédemment exposés. Cette extension vient s'insérer entre la première et la troisième étape du second mode de réalisation. Ainsi, lorsque le point d'impact est détecté, au moment de lancer l'analyse de la base de données matricielle en temps réel et avant qu'une étude locale de l'environnement du point d'impact ne soit effectuée, il est SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 19 - vérifié que ce point d'impact n'a pas encore été analysé, si ce n'est pas le cas et donc que l'environnement du point d'impact a déjà été analysé le rayon continue son trajet, si par contre c'est le cas, l'analyse du point d'impact est autorisée, une fois l'analyse exécutée l'environnement du point d'impact est mémorisé pour être réutilisé dans la suite du processus de lancer de rayons et le rayon continue son trajet. En conséquence, l'ensemble des étapes du second mode de réalisation reste quasiment identique, la seule modification étant l'utile mise en mémoire des résultats des calculs propres à l'environnement proche d'un point d'impact déjà analysé et donc des informations concernant les types d'obstacles et les angles formés par lesdits obstacles rencontrés. Les données matricielles sont alors linéarisées de proche en proche. De cette façon, grâce à cette variante il est permis de sauvegarder sous forme de données matricielles les données vectorielles calculées dans l'environnement proche du point d'impact, c'est-à-dire autour dudit point d'impact, afin que ces dernières, une fois sauvegardées, puissent être réutilisées subséquemment lorsqu'un autre rayon lancé vient rencontrer l'environnement proche du point d'impact déjà analysé. Ainsi, la base initiale d'objets contenant des informations vectorielles est avantageusement enrichie au fur et à mesure de l'exécution des calculs de lancer de rayons, ce qui a pour principal avantage de présenter un gain significatif en rapidité du procédé puisque ne sont pas ré-exécutés les calculs déjà effectués autour d'un point d'impact lors d'un lancer de rayon précédent. En effet selon cette variante, chaque fois qu'un rayon est lancé vers ou à proximité d'un point d'impact pour lequel des calculs ont déjà été effectués, la seconde SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 20 - étape d'analyse de l'environnement proche d'un point d'impact du second mode de réalisation du procédé n'est pas mise en oeuvre, ce qui revient à appliquer à l'identique le procédé selon le premier mode de réalisation avec prétraitement. Pour rappel, le calcul de l'atténuation globale Lb dans le cas de phénomène de diffraction horizontale basée sur le lancer de rayons est traditionnellement effectué en utilisant la formule suivante : Lb = Gant + los-k r.eiqk formule dans laquelle : - Gant représente le gain d'antenne, n est le nombre des différents segments de propagation en espace libre entre l'émetteur et le récepteur, chaque segment étant créé par une diffraction ou une réflexion, - Llos-k est l'atténuation en espace libre du segment d'indice k, - Diffk est la perte due à la diffraction entre le segment d'indice k et le segment d'indice k+1, - Reflk est la perte due à la réflexion entre le segment d'indice k et le segment d'indice k+1. Egalement et avantageusement, selon une mise en oeuvre préférée du procédé de prédiction du niveau de champ radio, lorsque le point d'impact d'un rayon est détecté, le couplet position du point d'impact et angle d'incidence du rayon sur le point d'impact peut être mémorisé dans une base dite de couplets, le rayon poursuivant alors son trajet par réflexion, ou diffraction/diffusion suite au point d'impact uniquement si le couplet point d'impact et angle d'incidence du rayon est significativement différent des couplets mémorisés lors des lancers de rayons effectués précédemment. SOFTW6-FR-1 TEXTE DEPOSE HG/ML - 21 - Pour mémoire, le procédé selon l'invention s'applique avantageusement à la prédiction du niveau de champ radio d'un réseau de téléphonie mobile, mais il pourrait tout aussi bien s'appliquer, simplement par analogie, à un procédé de simulation du niveau de champ électromagnétique basé sur un processus de lancer de rayons sur des données matricielles ou encore à un procédé de simulation du niveau de champ sonore également basé sur un processus de lancer de rayons sur des données matricielles. Diverses modifications à la portée de l'homme de métier sont bien entendu possibles sans pour autant sortir du cadre de l'invention. SOFTW6-FR-1 TEXTE DEPOSE HG/ML REVENDICATIONS1. Procédé de prédiction du niveau de champ radio d'un réseau de téléphonie mobile pour la couverture d'une cellule utilisant la technique de lancer de rayons et exploitant des bases de données géographiques caractérisé 5 en ce que la base de données directement exploitée est une base matricielle composée de pixels, chaque pixel portant au moins deux informations, une information propre au type de sursol et une information propre à la hauteur du sursol, à la base matricielle est appliqué un traitement pour 10 linéariser les contours d'obstacles rencontrés alors que des rayons sont fictivement lancés, les rayons selon les obstacles rencontrés étant soit réfléchis soit diffractés et diffusés, tandis que l'atténuation globale est calculée. 15 2. Procédé de prédiction du niveau de champ radio selon la revendication 1 caractérisé en ce que, avant le lancer fictif de rayons, un prétraitement de vectorisation des contours d'obstacles est appliqué aux bases matricielles pour linéariser les contours d'obstacles 20 rencontrés en ne prenant en compte que les objets représentant les obstacles et ainsi retrouver lesdits contours d'obstacles, le prétraitement de vectorisation étant adapté pour entraîner le lancer de rayons une fois que les informations vectorielles propres aux types, c'est- 25 à-dire, intérieur, façade ou coin d'obstacle et propres aux angles, c'est-à-dire, angle de façade ou angles de coin, ce pour chaque pixel, ont été extraites, pour former une matrice de types et une matrice des angles et prédire en fonction de ces matrices la réflexion, les diffraction et 30 diffusion pour chaque rayon lancé et permettre le calcul de l'atténuation globale. 3. Procédé de prédiction du niveau de champ radio selon la revendication 1 caractérisé en ce que, un rayon SOFTW6-FR-1 TEXTE DEPOSE HG/ML- 23 - fictif est lancé directement sur la base matricielle dans une direction d'un angle solide qui couvre complètement l'espace, son trajet étant en chaque point scruté et chaque fois que ledit rayon rencontre un obstacle et qu'un point d'impact est détecté, une analyse de la base de données matricielle est effectuée en temps réel pour créer des objets vectoriels adaptés au processus de lancer de rayons, une étude locale de l'environnement du point d'impact étant réalisée pour déterminer d'une part, si le point d'impact est un élément d'une façade ou d'un coin et d'autre part, l'angle de ladite façade ou les deux angles du coin de deux façades qui construisent le point, une fois l'environnement proche déterminé, le rayon continue son trajet, soit réfléchi si le point d'impact est un élément d'une façade, soit diffracté et diffusé si le point d'impact est le coin formé par deux façades. 4. Procédé de prédiction du niveau de champ radio selon la revendication 3 caractérisé en ce que lorsque le point d'impact est détecté au moment de lancer l'analyse de la base de données matricielle en temps réel et avant qu'une étude locale de l'environnement du point d'impact ne soit effectuée, il est vérifié que ce point d'impact n'a pas encore été analysé, si ce n'est pas le cas le rayon continue son trajet, si c'est le cas l'analyse du point d'impact est autorisée, puis une fois l'analyse exécutée l'environnement du point d'impact est mémorisé pour être réutilisé dans la suite du processus de lancer de rayons et le rayon continue son trajet. 5. Procédé de prédiction du niveau de champ radio selon l'une des revendications précédentes caractérisé en en ce que les rayons sont lancés dans toutes les directions horizontales en utilisant un incrément réduit, chaque SOFTW6-FR-1 TEXTE DEPOSE HG/ML- 24 - lancer de rayon s'effectuant de manière optimisée, c'est-à-dire qu'un rayon n'est lancé selon un angle donné que s'il est prévu que ce rayon aura un impact significatif sur la prédiction, une telle prévision étant obtenue par apprentissage lorsque des lancers de rayons précédents selon des angles voisins ont eu un impact significatif sur la prédiction. 6. Procédé de prédiction du niveau de champ radio selon la revendication 5 caractérisé en ce que la prévision qu'un lancer de rayon selon un angle donné aura un impact significatif sur la prédiction résulte de la comparaison à un seuil du nombre de pixels pour lesquels la prédiction a été mise à jour par le rayon. 7. Procédé de prédiction du niveau de champ radio selon l'une des revendications précédentes caractérisé en ce que lorsque le point d'impact d'un rayon est détecté, le couplet position du point d'impact et angle d'incidence du rayon sur le point d'impact est mémorisé dans une base de couplets et le rayon poursuit son trajet par réflexion, ou diffraction/diffusion suite au point d'impact uniquement si le couplet point d'impact et angle d'incidence du rayon est significativement différent des couplets mémorisés lors des lancers de rayons effectués précédemment. SOFTW6-FR-1 TEXTE DEPOSE HG/ML