FR 2828620 A1 20030214 FR 0110718 A 20010810 13) Valider le plan de fréquences. 1 2828620 Procédé d'estimation prévisionnelle de la couverture radio d'u ne cell u le La présente invention concerne la phase de définition des réseaux de téléphonie cellulaire à implanter sur un territoire et en particulier un procédé d'estimation prévisionnelle de couverture radio du territoire, en vue de déterminer des paramètres de fonctionnement et des positions optimales des stations de base ou répéteurs du réseau, c'est-à-dire les limites des cellules radio correspondantes. On rappellera qu'un réseau radiotéléphonique cellulaire est constitué d'une pluralité de stations terrestres de base qui sont interconnectées à travers le réseau téléphonique filaire et auxquelles les terminaux mobiles peuvent accéder lorsqu'ils se trouvent dans la cellule radio de la station. La propagation radio dans une cellule doit satisfaire aux deux exigences essentielles que sont l'émission d'une puissance non excessive par la station de base et la réception par les terminaux de signaux radio de puissance 1 5 suffisante. En effet, la portée de chaque station de base doit, tout d'abord, être sufffisante pour que la cellule considérée déborde sur les voisines, afin d'éviter tout risque de coupure de communication lorsqu'un terminal mobile change de cellule. Cela nécessite d'augmenter la puissance d'émission des stations au delà du strict nécessaire. Ensuite, comme les liaisons radio ont un trejet sensiblement linéaire au niveau du sol, il faut couvrir les zones d'ombre radio de la station, dues aux reliefs naturels ou aux immeubles. Une zone d'ombre est une zone dans laquelle l'atténuation de propagation radio, entre un terminal mobile qui s'y trouve et la station, dépasse la limite de la spécification de sensibilité des circuits radio, si bien que le niveau de réception est insuffisant pour détecter co rrecteme nt les paq u ets de bits émis, représe nta nt la pa role ou d es do n nées à échanger. Or, il est exclu d'augmenter les niveaux d'émission. En effet, côté station de base, toute augmentation accroA'trait la taille de la cellule, provoquant des interférences excessives entre cellules voisines. Côté 2 2828620 terminaux mobiles, leur puissance maximale est limitée par des impératifs de sécurité des personnes et d'autonomie de leur batterie. Il faut aussi éviter de multiplier inutilement le nombre de stations de base ou répéteurs pour les micro-cellules que représentent les zones d'ombre, afin de limiter les coGts et les interférences. De façon classique, pour une cellule, on effectue 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 l'lnstitut Géographique National (IGN), entre autres, et représentant la carte de la zone géographique considérée avec les bâtiments et autres reliefs du sursol qu'elle comporte. Différents mots de code définissent le type de sursol, tel que forêt, eau, pavillons, spécifiés en coordonnées Lambert et en hauteur au dessus de l'altitude locale du sol par rapport à la mer. Pour estimer l'atténuation prévisionnelle en tout point de la cellule, on effectue une simulation de propagation de signaux radioélectriques. On modélise la propagation en lançant fictivement, par calcul dans un calculateur d'une chane de production de prévisions de couverture radio, un rayon électromagnétique depuis la station de base, dans une direction déterminée, et on calcule ses conditions de propagation dans l'angle solide élémentaire qu'il occupe. Hormis la propagation en espace libre lorsque la station est en vue directe d'un terminal radio fictif, propagation d'atténuation linéaire connue, la trejectoire du rayon rencontre des obstacles qui l'atténuent ou 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, le rayon peut être dévié par réflexion ou réfraction. L'ouverture angulaire de son angle solide peut méme s'en trouver accrue. On répète ces calcuis pour une pluralité d'angles solides élémentaires répartis dans un angle solide global de visée de toute la cellule, tel qu'une couronne globalement horizontale, afin d'ainsi échantillonner les diverses conditions de propagation de 1'espace de la cellule. En chaque point de la trejectoire de chaque rayon, le calculateur consulte la base de données vectorielle pour déterminer s'il comporte un obstacle. Les 3 2828620 calcuis vectoriels correspondants nécessitent une puissance de calcul considérable et s'étendent sur une journée de travail, en pratique ils doivent être lancés le soir, lorsqu'on utilise des ordinateurs classiques. La présente invention vise à réduire la puissance de calcul nécessaire au calcul de la couverture radio de telles cellules, que ce soient des cellules de grande taille, ou macrocellules, ou des microcellules. A cet effet, I'invention concerne un procédé d'estimation prévisionnelle de la couverture radio d'une cellule, d'un réseau de radiotéléphonie cellulaire, par une station radio de gestion de trafic de la cellule, au moyen d'une base de données de carte en relief spécifiant des positions et natures de reliefs de la cellule, procédé dans lequel, un système de calcul dans une phase d'exploitation, - lance fictivement dans la cellule un rayon d'échantillonnage représentatif des conditions de propagation radioélectrique, selon un tronçon initial de trejectoire depuis une position d'origine et selon une direction et des conditions de propagation déterminées, - par lecture de la base de données, le système de calcul compare le tronçon de trajectoire aux données de carte en relief pour identifier la position et la nature de tout point courant d'impact du tronçon de trejectoire avec un relief, - d'après la nature spécifiée dans la base de données du relief considéré, le système de calcul détermine de nouvelles conditions de propagation sur un tronçon aval de trajectoire au-delà du point d'impact, - le système de calcul réitère, le cas échéant, les deux étapes précédentes un nombre déterminé de fois pour d'autres impacts de tronçons aval de trajectoire, - d'après les conditions de propagation sur la totalité de la trajectoire, le système de calcul détermine une atténuation cumulée pour tout point sélectionné de celle-ci, et - le système de calcul répète toutes les étapes ci-dessus une pluralité de fois pour une pluralité de directions initiales afin d'échantillonner toute la cellule et ainsi déterminer une carte d'atténuations des points sélectionnés, procédé caractérisé par le fait que, dans une phase initiale, 4 2828620 - le système élabore la base de données sous la forme d'une matrice de données, à accès direct, de spécification locale et autonome des positions et natures des reliefs respectifs d'une pluralité de mailles prédéterminées d'un maillage correspondant de la carte, et, dans la phase d'exploitation, - compare la position du point courant au maillage pour identifier une maille d'impact, et - calcule les conditions de propagation après impact à partir des données de spécification locale et autonome de relief de la maille d'impact. Les données de propagation ayant ainsi été fragmentées en une pluralité de données autonomes, et donc de volume réduit, il est alors facile d'y accéder rapidement pour les lire sous une forme exploitable. Ainsi, le calcul des conditions de propagation en aval du point d'impact s'effectue à partir uniquement de données locales et de volume limité, à accès direct et non pas séquentiel comme dans le cas d'une base vectorielle, données qui définissent, pour le point d'impact, les conditions d'anisotropie radioélectrique déterminant, par atténuation, réflexion ou diffraction avec éventuellement diffusion, une éventuelle déviation angulaire et une nouvelle atténuation de propagation. En l'absence de relief formant obstacle, les conditions de propagation en espace libre sont parfaitement connues. En particulier dans ce cas, une lecture des données d'une maille traversée montre immédiatement qu'il n'y a pas d'obstacle et on passe donc à la maille suivante, sans nécessiter les calcuis très lourds, de l'art antérieur, de reconstitution vectorielle des donnces locales à partir d'une base de données globale On notera que le terme relief ici considéré désigne globalement tout obstacle à la propagation, y compris le sol, même horizontal, qui est susceptible d'absorber ou de réfléchir au moins partiellement les rayons. Avantageusement, le système mémorise, dans la phase initiale, des données d'orientation géographique des reliefs dans la base de données matricielle et, en exploitation, calcule une direction d'un rayon réfléchi de tronçon aval à partir des données d'orientation. 2828620 De préféren ce, dans ce cas, la mémorisation se l i mite à des do n nées d'azimut et, en exploitation, la direction du rayon réfléchi de tronçon aval est calculée en considérant que les reliefs de réflexion sont verticaux. Les calcuis sont ainsi limités. En particulier, il est préférable d'intégrer, dans la base de données matricielle, des données spécifiant une nature de sursol et, en exploitation, le système calcule les conditions de propagation sur le tronçon aval d'après ladite nature de sursol. La précision des calcuis de propagation est ainsi meilleure. Pour tenir compte de détails des reliefs, dans la phase initiale, on peut intégrer des données relatives à des arétes des reliefs dans la base de données matricielle et, en exploitation, le système de calcul calcule une direction de tronçon aval de rayon réfracté à partir desdites données d'arétes. De préférence encore, dans la phase initiale, des données d'atténuation des reliefs sont mémorisées dans la base de données matricielle et, en exploitation, le système calcule une atténuation de propagation sur le tronçon aval à partir desdites données d'atténuation. En pareil cas, les donnces d'atténuation peuvent concerner la réflexion sur les reliefs et elles sont utilisées pour calculer l'atténuation sur le tronçon aval, eVou les données d'atténuation concernent la propagation à travers les reliefs et elles sont utilisoes pour calculer l'atténuation des rayons qui s'y propagent. Avantageusement, les données d'atténuation de propagation à travers les reliefs comportent en outre des données de transition entre milieux de propagation, spécifiant des atténuations de pénétration dans les reliefs, qui sont utilisées pour déterminer une atténuation locale de pénétration. Pour une modélisation de qualité accrue des conditions de propagation, et disposant d'un algorithme de calcul de dispersion angulaire du rayon après impact, le système calcule une pluralité de directions de tronçons aval, avec des atténuations spécifiques, formant un angle solide de diffusion du rayon au delà de l'impact. On peut en outre prévoir de compter les impacts successifs sur la trajectoire et, en cas de deuxième impact, le système considère, pour calculer 6 2828620 lesdites conditions de propagation au-delà, que le rayon a été polarisé lors du premier impact. Afin de limiter les calcuis requis, le système compte les impacts successifs sur la trejectoire et compare le total à une valeur de seuil haut pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint, ou bien, à chaque impact, le système détermine l'atténuation de propagation cumulée et la compare à une valeur de seuil d'atténuation maximale pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint. Pour exploiter au mieux les résultats, la carte des atténuations est mémorisée en trois dimensions. On peut ainsi estimer un niveau de qualité prévisionnel des liaisons selon les étages des immeubles. En particulier, pour élaborer la base de données, le système représente la carte par un faisceau de chanes de pixels à extension verticale et découpe les chanes de pixels pour constituer des volumes élémentaires empilés de mailles volumiques ayant chacun des donnces particulières. De préférence, le système prend, comme position d'origine, une position prévue pour la station, mais, toutefois, peut prévoir de prendre, comme position d'origine7 une position quelconque dans la cellule et la direction de lancement du rayon est choisie, d'après les positions et natures des reliefs proches, pour que le rayon passe à proximité de la station. Ayant calculé, selon le procédé de l'invention, les conditions de propagation dans une microcellule en contact avec une cellule, donc plus grande et encore appelée macrocellule, le système peut alors calculer des conditions de propagation dans la cellule et effectue ensuite un calcul de lissage des résultats des deux calcuis relatifs à une zone de frontière entre cell u l e et microcel l u le. Le procédé permet ainsi de mieux intégrer fonctionnellement les deux types de cellules. L'invention sera mieux comprise à l'aide de la description suivante d'un mode préféré de mise en _uvre de l'invention, en référence au dessin annexé, sur lequel 7 2828620 - la figure 1 représente une carte géographique en relief formant une base de données, sur laquelle a été reportée la position d'une cellule d'un réseau de radiotélophonie cellulaire en cours de définition, - la figure 2 représente une courbe d'atténuation radio dans la cellule, sur une trejectoire comportant des obstacles et dont la station de base constitue l'une des extrémités, - la figure 3 est une vue en plan du contour d'un immeuble sur une portion de la carte, - la figure 4 est homologue de la figure 3 et illustre une réflexion et une réfraction sur un immeuble, et la figure 5 est une vue en coupe verticale illustrant la propagation entre des immeubles d'une macrocellule. La figure 1 représente la position d'une cellule 1, d'un réseau de radiotélophonie cellulaire prévu pour des terminaux comme celui référencé 21, reportée sur une portion de carte géographique 31 comportant aussi la position prévue d'une station radio de base 11 de gestion de la cellule 1. La carte 31 sert à déterminer des contours prévisionnels d'une pluralité de cellules constituant le réseau radio à implanter, en ajustant leur nombre, tailles et positions afin d'optimiser le volume de matériel tout en assurant la couverture radio voulue avec une qualité de service spécifiée. La référence 31 désigne en fait ici la présentation des éléments du terrain, sol et sursol, de la zone considérée. Les données géographiques correspondantes permettant d'engendrer cette présentation sont stockées dans une base de données de relief 32, 33 d'un calculateur 30. Cette base de données comporte un bloc mémoire 32 spécifiant la forme et le relief du terrain, sol et sursol, associé à un bloc mémoire 33 de morphologie du sol et sursol spécifiant les natures ou caractéristiques de propagation radio des diverses occupations du sol, dans une plage de fréquences déterminée correspondant aux fréquences utilisées par les stations de base, ces données de morphologie du bloc mémoire 33 présentent les divers points du terrain considéré avec les reliefs ou formes du bloc géographique 32. On notera que le terme "relief" doit être pris dans une acception large pour désigner tout obstacle susceptible 8 2828620 d'être atteint par un rayon, direct ou préalablement dévié, de la station 11. Il s'agit donc globalement du sol et du sursol. Outre les bâtiments, collines et équivalents, il peut donc aussi s'agir de flancs de vallées ou encore de plaines ou plans d'eau. Comme expliqué plus loin, les données du bloc 33 permettent de déterminer, directement ou par une table de correspondance entre la nature et des données de caractéristiques de propagation radio, la perturbation apportée à un rayon radio incident, afin de déterminer la direction et l'atténuation ou amplitude d'un rayon aval correspondant. La carte 31 en elle-même n'a donc ici qu'un but didactique, puisque les données la définissant sont contenues dans les blocs 32, 33, qu'un calculateur exploite. En variante, la station de base 11 pourrait être remplacée par une station de même fonction, mais de portée réduite, pour définir une microcellule. Comme indiqué au début, les microcellules sont implantées dans les zones de forts reliefs ou en ville, pour couvrir les zones d'ombre radio des cellules class iq u es. On suppose ici que la trejectoire d'un rayon radio, ou faisceau étroit, reliant la station 11 au mobile 21 dans la cellule 1, rencontre des obstacles 41 et 42, respectivement un bâtiment et les arbres d'une forêt. Pour la simplicité de I'exposé, on suppose ici que les obstacles 41, 42 de la figure 1 ne dévient pas la trajectoire, qui reste donc rectiligne, sans réflexion ou réfraction par ceux-ci. La figure 2 représente en ordonnée, en décibels (dB), le niveau S de signal radio en fonction de la distance X parcourue sur la trejectoire, portée en abscisse. Partant du niveau d'émission, la décroissance du niveau avec la distance de propagation représente donc l'atténuation. Cette atténuation est la somme des atténuations de divers segments ou tronçons de la trajectoire, chaque tronçon correspondant à un milieu de propagation spécifique. On distingue ici cinq tronçons, référencés dans l'ordre 51 à 55, correspondant respectivement pour le premier 51 à l'atténuation due au trajet aérien allant de la station 11 à un bâtiment 41, pour le second 52 à l'atténuation due à la traversée du bâtiment 41, pour le troisième 53 à l'atténuation due au trajet aérien qui suit jusqu'au bord d'une forêt 42, pour le quatrième 54 à 9 2828620 I'atténuation due à la traversée de la forêt 42 et pour le cinquième à l'atténuation due au trajet aérien jusqu'au terminal mobile 21. L'atténuation linéique de propagation radio en espace libre dans l'air, représentée par la pente du niveau sur la figure 2, est une constante physique bien connue, pour une fréquence de porteuse déterminée, et le calculateur 30 peut donc, à partir des distances de propagation relevées dans la base de données représentative de la carte 31, calculer les trois atténuations correspondantes. Par contre, les obstacles 41, 42 correspondent à des atténuations accrues. En outre, comme évoqué, certains types d'obstacles peuvent aussi dévier la trajectoire, comme illustré par la figure 4. Or, pour vérifier la couverture radio prévisionnelle, il faut pouvoir estimer par le calcul le niveau de champ en divers points répartis dans toute la cellule 1, pour les liaisons entre la station de base 11 et le terminal mobile 21. Le bilan de liaison prévisionnel ainsi calculé, ou somme des atténuations associées aux tronçons de trejectoire 51 à 55, ne doit pas excéder l'écart entre un niveau maximal Nm d'émission de la station 11 et un seuil prédéterminé Sm de sensibilité du terminal 21. Il en est de même, pour le sens dit montant de communication du terminal 21 vers la station 11, en ce qui concerne un niveau d'émission maximal du terminal 21, ici 2 watts, et un niveau de sensibilité de la station 11. Ces niveaux de sensibilité tiennent compte d'un code de détection des erreurs de propagation, et de correction d'un nombre limité de bits erronés, dans les paquets de bits échangés à travers des voies temporelles d'une trame radio. Dans un but d'efficacité, les calcuis d'atténuation sur ordinateur ne doivent nécessiter qu'une durée limitée, bien moindre qu'une journce de travail. Pour ce faire, la carte 31 se présente sous la forme d'une carte géographique numérique en relief représentant la cellule 1 et la base de données 32, 33 spécifie les positions des reliefs comme 41 et 42 et autres, et leurs natures déterminant les conditions de propagation radio, telles que bâtiments élevés, forêt, zone pavillonnaire, lac et autres. Dans une phase initiale, avant la phase d'exploitation de calcul des atténuations, un maillage prédéterminé de la carte 31 est effectué et, à chaque maille 34 de la matrice 2828620 géographique ainsi déterminée, on associe les données numériques des reliefs qu'elle comporte, pour ainsi disposer d'une base de données matricielle, ou mosarque, de spécification des reliefs, utilisable dans la phase d'exploitation. En pratique, un opérateur humain ou le calculateur 30 définit le maillage et stocke, dans les blocs 32 et 33, respectivement les données de forme de relief de chaque maille 34 et les données de nature des reliefs. Comme évoqué plus haut, ces données identifient les reliefs, tels que forét, bâtiments, et une table générale de correspondance permet d'en déterminer des valeurs de caractéristiques radio, c'est-à-dire de paramètres déterminant les conditions de propagation des rayons atteignant le relief considéré. En variante, les valeurs des paramètres de propagation sont directement stockées dans le bloc 33, sans nécessité de stocker les types des reliefs. Dans la phase d'exploitation, exécutée en pratique par le calculateur 30, pour l'estimation prévisionnelle de la couverture radio de la cellule 1 par la station de base radio 11 de gestion de trafic de la cellule 1, au moyen de la base de données 32, 33 de carte en relief spécifiant les positions et natures de reliefs comme 41 et 42 de la cellule 1. Le calculateur lance fictivement dans la cellule un rayon d'échantillonnage représentatif des conditions de propagation radioélectrique, selon un tronçon initial de trajectoire par exemple 11, 41 depuis une position d'origine et selon une direction et des conditions de propagation déterminées, c'est-à-dire en pratique à travers l'air. Par lecture de la base de données, le calculateur compare le tronçon 11, 41 de trajectoire aux données de carte en relief, précisément ici celles du bloc 32, pour identifier la position et la nature de tout point courant d'impact du tronçon de trajectoire avec un relief 41 dans l'exemple. D'après la nature spécifiée du relief considéré 41, le calculateur détermine de nouvelles conditions de propagation sur une portion aval de trejectoire 41, 42 (tronçon 52) au-delà du point d'impact en 41. Le calculateur réitère, le cas échéant, les deux étapes précédentes un nombre déterminé de fois pour d'autres impacts de portions aval 53, 54 de trajectoire. 11 2828620 D'après les conditions de propagation sur la totalité de la trajectoire 11 à 21, le calculateur détermine une atténuation cumulée pour tout point sélectionné de celle-ci, et répète toutes les étapes ci-dessus une pluralité de fois pour une pluralité de directions initiales afin d'échantillonner toute la cellule et ainsi déterminer une carte d'atténuations des points sélectionné. En outre, ayant préalablement élaboré, dans la phase initiale, la base de données 32, 33 sous la forme d'une matrice, ou mosaTque, de données de spécification locale et autonome de la pluralité de mailles prédéterminces 34 du maillage correspondant de la carte 31, cela permet, dans la phase d'exploitation, - de comparer la position du point courant au maillage pour identifier la maille d'impact 34, et - de calculer les conditions de propagation après impact (tronçon 41-42 par exemple) à partir des données de spécification locale et autonome de relief de la maille d'impact 34. Ainsi, d'une façon générale, la base de données 32, 33 spécifie, dans chaque maille, notamment dans la mémoire 33 I'anisotropie radio locale, c'est- à-dire les distorsions apportées au rayon, telles que déflexion, atténuation, diffraction polarisation et autres. Il est en quelque sorte défini, pour chaque maille 34, un ellipsode d'anisotropie déterminant les conditions de propagation dans l'espace, ellipsoide rapporté à trois directions propres, par exemple trois vecteurs unitaires orthogonaux d'abscisse, ordonnée, tels que parallèle, méridien et verticale locaux. L'ellipsode ci-dessus est en fait multiple, car il spécifie les valeurs de plusieurs variables de propagation, permettant de calculer par exemple la direction d'émergence d'un rayon incident en fonction de son angle d'arrivée, sur la carte 31, ou encore l'atténuation correspondante de traversée de la maille 34, en fonction des deux directions, d'incidence et d'émergence. Il s'agit donc d'une matrice de transformation des conditions de propagation des diverses mailles, dont les données sont stockées par zones indépendantes dans le bloc 33. 12 2828620 La lecture, dans le bloc 33, de la zone appropriée associée à la maille d'impact 34 permet donc de déterminer rapidement les conditions de propagation de sortie de la maille à partir du tronçon de trajectoire d'entrée par exemple 11, 41. En particulier, un bit par maille 34 peut spécifier si la maille 34 considérée contient ou non un obstacle. Par lecture directe de celui-ci, éventuellement rangé avec ses homologues dans un registre de taille limitée et à accès rapide, une absence d'indication d'obstacle est ainsi immédiatement détectée et le calculateur 30 passe immédiatement à l'examen de la maille suivante, sans calcul de nouvelles conditions de propagation. Ainsi, hormis la lecture du bit indicateur de présence d'obstacle, le calculateur 30 ne consulte alors pas le bloc de données 33 spéciflant la nature des obstacles radio. Un calcul d'atténuation de point courant, de maille en maille, sur la trajectoire n'est pas nécessaire. L'atténuation cumulée peut en effet être calculée uniquement en cas d'obstacle, par calcul de la distance entre les deux mailles 34 d'extrémité du tronçon considéré. L'expérience a montré que les calcuis pour une cellule de 500 mètres de rayon moyen nécessitent environ 1 minute de calculateur de puissance moyenne. Le calculateur 30 ayant fixé ou calculé l'angle du rayon incident, dans l'espace de la carte 31, peut ainsi directement lire, dans le bloc 33, toutes les valeurs correspondantes des paramètres de propagation du rayon sortant de la maille 34, si le rayon peut effectivement en sortir. On notera que, la carte 31 étant en relief, le maillage des donnces est de préférence, comme dans cet exemple, effectué en trois dimensions indépendantes, comme celles évoquées ci-dessus. En d'autres termes, on peut définir un maillage bidirectionnel de mailles 34 dites horizontales, à chaque maille horizontale 34 correspondant un volume à extension verticale découpé à diverses altitudes, éventuellement spécifiques à chaque maille horizontale 34, par des plans ou autres surfaces, afin de déterminer des volumes élémentaires ayant chacun des données de propagation particulières stockées dans une zone du bloc 33. Dans la plupart des cas, il suffit de deux volumes élémentaires par maille horizontale 34, celui de dessous contenant par exemple la totalité d'un immeuble, et celui de dessus correspondant à l'espace libre. Dans le cas 13 2828620 par contre de reliefs avec surplomb, tels qu'immeubles en arche ou ponts, un troisième volume élémentaire, à propagation libre, est à prévoir sous un volume élémentaire contenant l'obstacle considéré. En d'autres termes, on représente la carte 31 par un faisceau de chanes depixels à extension verticale, pour chaque maille 34, et on découpe les chanes de pixels pour constituer des volumes élémentaires empilés de mailles volumiques ayant chacun des données particulières. Les donnces de spécification des caractéristiques ou nature des reliefs de chaque maille 34, intégrées lors de la phase in itiale, peuvent correspond re à I'une ou plusieurs des données suivantes. Les données de spécification des reliefs de la mémoire 32 peuvent consister en des données d'orientation géographique des reliefs, indiquant par exemple un plan de réflexion radio. Connaissant la portion ou segment de trejectoire incidente, on en déduit par calcul la direction de la portion de trajectoire aval réfléchie, qui en est le symétrique par rapport à la normale au plan de réflexion au point d'impact. Dans l'air, en l'absence d'obstacle dans le premier ellipsode de Fresnel (propagation directe), I'atténuation en champ proche est d'environ 20 dB/km sur les 500 premiers mètres de trajectoire; au delà, on prend la valeur 30 dB/km. La figure 3 illustre ainsi le contenu d'une table de donnses en matrice du bloc 33, contenu qui est présenté schématiquement sous la forme graphique d'une partie de la carte 31, afin d'être expliqué plus clairement. La figure 3 représente quelques mailles horizontales 34 sur lesquelles on a reporté la surface occupée par un immeuble. L'immeuble considéré comporte quatre murs rectilignes, référencés 61 à 64, et présente une forme de trapèze avec quatre coins de murs ou arêtes verticales 65 à 68. Dans le bloc 33, les quatre mailles comportant l'une des arêtes de mur comportent des données spécifiant cette caractéristique. Il peut en outre être prévu de spécifier la valeur de l'angle de l'arête et même l'orientation de ses côtés. Les autres mailles 34 traversoes par l'un des murs 61 à 64 comportent des données spécifiant cette particularité. En pratique ici, ces données spécifient l'orientation du mur considéré, c'est-à-dire sa direction azimutale. Il 14 2828620 peut être prévu de spécifier aussi l'inclinaison du relief lorsqu'il s'agit de reliefs autres, par exemple naturels. En pareil cas, ces données d'orientation du plan du relief peuvent être définies par la normale au relief, exprimoe par exemple selon les coordonnées tridimensionnelles évoquées plus haut. Pour toutefois limiter la taille du bloc 33, il peut être prévu de limiter les données d'orientation du relief aux données d'azimut, et, en exploitation, la direction du rayon réfléchi de tronçon aval est calculée en considérant que les reliefs de réflexion du rayon considéré sont verticaux, ce qui est le cas général en ville. On peut aussi intégrer, dans le bloc 33 de la base de données matricielle, des données spécifiant une nature de sursol, tel que forêt, zone pavillonnaire et autres, et, en exploitation, le système 30 calcule les conditions de propagation sur le tronçon aval d'après la nature de sursol telle que spécifiée. Dans la phase initiale, si l'on intègre les données d'arêtes des reliefs dans le bloc 33 de la base de données matricielle, cela permet, en exploitation, de calculer une direction de tronçon aval du rayon réfracté à partir des données d'arêtes. Comme indiqué, la réfraction est provoquée en général par des coins ou arêtes verticales de bâtiments. Toutefois, des bords et arêtes de crête de toitures peuvent de même étre spécifiés dans le bloc 33 pour déterminer de même une direction de réfraction du rayon incident, ainsi défléchi vers le bas par exemple. Les déflexions cidessus accroissent la taille de la zone de couverture de la cellule 1 puisqu'elles dirigent le rayon réfracté vers un volume qui, en propagation linéaire, serait une zone d'ombre. Dans la phase initiale, le système intogre dans la mémoire des données d'atténuation des reliefs dans la base de donnces matricielle 33 et, en exploitation, le système détermine une atténuation du rayon à partir des données d'atténuation ci-dessus. Les données d'atténuation peuvent concerner la réflexion sur les reliefs 60 et elles sont alors utilisées pour calculer l'atténuation de rayons réfléchis, par exemple environ 7 dB, valeur qui dépend de la nature morphologique de la façade, telle que verre, brique ou autre. 2828620 En outre ou à la place, les données d'atténuation peuvent concerner la propagation à travers les reliefs 41, 42, 60 et elles sont alors utilisoes pour calculer l'atténuation des rayons qui s'y propagent, comme illustré sur la figure 2. En particulier, les données d'atténuation de propagation à travers les reliefs 41, 42 peuvent comporter en outre des données de transition de milieux de propagation, spécifiant des atténuations de pénétration dans les reliefs, ou de changement de milieu de propagation, qui sont utilisées pour déterminer une atténuation locale de pénétration, air / bâtiment par exemple. La figure 4 est homologue de la figure 3 mais l'immeuble 70 qu'elle représente a, en vue de dessus, une forme simplifiée, triangulaire dans cet exemple, pour la simplicité de l'exposé. Un rayon incident 81 atteint un point d'une façade 71 de l'immeuble 70, point d'impact qui est situé dans une maille 34 traversée par la façade 71. D'après les données correspondantes de la maille du bloc 33, indiquant la direction azimutale de la façade 71 et indiquant que la maille 34 considérée est (entièrement) traversée par la façade 71, le calculateur 30 détermine que le rayon incident 81 s'y réfléchit selon un rayon 82 dont il calcule la direction d'émergence, direction qui forme, avec la normale locale 72 à la façade 71, un angle égal et opposé à celui du rayon incident 81. Comme représenté, le rayon réfléchi 82 définit en fait la direction principale d'un lobe 84 de rayons annexes 83 qui l'entourent, formant donc, dans l'espace, un angle solide. En effet, le rayon 81 excite la zone qu'il atteint et il engendre ainsi une source secondaire qui diffuse le rayonnement 2s électromagnétique selon un diagramme plus diffus, c'està-dire plus isotrope, qu'une source primaire classique. Afin de mieux modéliser la propagation, dans cet exemple le système détermine, par un algorithme de calcul de dispersion angulaire, la pluralité de directions des rayons 82 et 83, avec des atténuations spécifiques, formant I'angle solide de diffusion du rayon incident 81 audelà du point d'impact. Ceci peut concerner les rayons réfléchis par des surfaces inhomogènes, telles que façades d'immeubles avec fenêtres et balcons, et les rayons diffractés. Pour 16 2828620 illustrer ces derniers, un rayon incident 91 atteint encore la façade 71 mais dans une zone d'arête verticale 73. Les données du bloc 33 de la maille 34 considérée spécifient alors la présence de l'arête, en indiquant une direction de diffraction principale pour un rayon émergent 92 et un angle solide 94 de rayons diffractés annexes 93. L'explication ci-dessus est aussi valable pour des arêtes inclinées sur la verticale. On peut donc en particulier aussi considérer la figure 4 comme étant une coupe verticale d'un relief, à travers des rangées horizontales superposées de mailles, illustrant le fait que le rayon diffracté peut " plonger vers le sol, dans un volume que l'on pourrait a priori supposer être une zone d'ombre radio. On notera que le bloc 33 peut contenir à la fois des données d'orientation pour calculer une réflexion partielle (82) et des données de calcul d'un rayon réfracté 92, si l'on considère que le rayon incident 81, 91 présente une section de l'ordre de grandeur d'une maille 34, qui n'est affectée que partiellement par la présence de l'arête 73. En référence à la figure 5, I'atténuation de diffraction sur un toit Ltàm (toit à mobile), peut se calculer par la formule suivante: Ltàm = -16,9 10 log (W) + 10 log (f) + 20 log (hb- Hm) ou pris = 0, si le calcul cidessus donne Ltàm I'atténuation totale, I'atténuation a une valeur réduite lors des impacts suivants. Le rayon a ainsi en quelque sorte été adapté aux obstacles. L'atténuation de réflexion passe ainsi, en moyenne, de 7 dB pour la première réflexion à 3 dB pour les suivantes. Afin de disposer de données de couverture radio aux divers étages des immeubles et pour aussi traiter les formes en surplomb évoquées plus haut, dans cet exemple on établit en trois dimensions la carte ou base terrain des atténuations. Pour éviter de prolonger inutilement le temps de calcul, le système compte les impacts successifs du rayon et on compare le total à une valeur de seuil haut pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint. Dans le même but, en complément ou à la place, à chaque impact, le système détermine l'atténuation cu mulée et la compare à u ne valeu r de seu il d'atténuation maximale pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint. On remarquera que l'on peut, en variante de cet exemple particulier, appliquer le principe de retou r inverse du rayon, c'est-à-d ire l'émettre depu is la position quelconque du terminal 21 en direction de la station 11, en déplaçant comme précédemment la position du terminal 21 dans toute la cellule 1. Toutefois, en pareil cas, il existe une certaine incertitude pour le succès de chaque essai, puisqu'un rayon émis en direction de la station 11 risque d'être dévié et que, inversement, un rayon émis dans une autre direction subira une déviation par un relief et l'atteindra. Il faudra donc prévoir un nombre accru de lancers de rayons, par exemple dans un grand angle solide contenant la station 11, compte tenu des déviations prévisibles, par exemple un toit d'immeuble à effet de réfraction, proche du terminal 21. Ainsi, on peut prendre, comme position d'origine, une position prévue pour la station 11, ou bien on prend, comme position d'origine, une position quelconque dans la cellule 1 et la direction de lancement du rayon est choisie, d'après les positions et natures des reliefs proches, pour qu'elle passe à proximité de la station 11. Ainsi, une proportion sumsante de calcuis est exploitable. 18 2828620 Le procédé peut en particulier étre utilisé pour des microcellules en contact avec, si non incluses dans, des (macro)cel l u les. Com pte ten u de la hauteur élevée des antennes de ces dernières, la propagation y est moins perturbée et les calcuis peuvent y étre effectués par un procédé classique de S propagation. Ayant aussi calculé, selon le présent procédé, les conditions de propagation dans une microcellule en contact avec une (macro)cellule, on calcule des conditions de propagation dans la cellule et on effectue ensuite un calcul de lissage des résultats des deux calcuis relatifs à une zone de frontière entre cell u le et m icrocell u le La carte 31 a ici un pas de maillage du terrain d'environ 5 mètres selon des parallèles et méridiens. Elle peut en particulier étre obtenue par interpolation linéaire à partir d'une carte altimétrique de l'IGN à plus petite échelle, à mailles sensiblement sous la forme de grands carrés de 50 mètres de côté, ayant une base de données vectorielle définissant, outre l'altimétrie du terrain, les positions des reliefs et leurs natures. Le système de calcul maille celle-ci au pas de 5 mètres, en découpant chaque grand carré en cent petits carrés de 5 mètres de côté. La surface délimitée par chaque petit carré détermine ainsi un sous- ensemble correspondant des données définissant les positions et natures des reliefs. Ensuite, le système de calcul effectue un lissage, ou filtrage spatial passe bas, par un calcul d'interpolation prenant en compte les données d'altimétrie ci dessus des grands carrés voisins de celui considéré, ceci en ligne, en colonne et en diagonale. Le calculateur 30 module ainsi par exemple les données d'altitude du sol moyen du grand carré considéré pour y déterminer une valeur locale la plus probable pour chaque petit carré, pour obtenir ainsi un sous-bloc matriciel de données d'altimétrie, constituant l'une de la pluralité de zones du bloc 32. D'autres données d'ordre global, par exemple de spécification du sursol, tel que forét, zone pavillonnaire, terrain vide ou autres, peuvent étre calculées de cette façon. Les données d'ordre plus ponctuel, spécifiant par exemple une orientation de façade d'immeuble pour le bloc 33, sont par contre, de préférence, établies d'après des relevés effectués sur le terrain, par exemple 19 2828620 des photos adbennes. Le sous-oc mMcl de donndes de sumol em de pn=, ammo id, ample par un sous-bc matMcl de donndes de hour du sumoL opaques par In entry #Jude de demier Fa##ude du soL par mppod 18 men Les donnAes des Boa 32,33 sont nsi plus praises et tour. 2828620 Revendications 1. Procédé d'estimation prévisionnelle de la couverture radio d'une cellule (1), d'un réseau de radiotéléphonie cellulaire, par une station radio (11) de gestion de trafic de la cellule (1), au moyen d'une base de données (32, 33) de carte en relief (31) spécifiant des positions et natures de reliefs (41, 42) de la cellule (1), procédé dans lequel, un système de calcul dans une phase d'exploitation, - lance fictivement dans la cellule (1) un rayon d'échantillonnage représentatif des conditions de propagation radioélectrique, selon un tronçon initial de trajectoire (51) depuis une position d'origine et selon une direction et des conditions de propagation déterminées, - par lecture de la base de données (32, 33), le système de calcul compare le tronçon de trajectoire (51) aux données de carte en relief (31) pour identifier la position et la nature de tout point courant d'impact du tronçon de trajectoire (51) avec un relief (41, 42), - d'après la nature spécifiée dans la base de données du relief considéré (41, 42), le système de calcul détermine de nouvelles conditions de propagation sur u n tronçon aval de trejectoire (52) au- delà du point d' impact, - le système de calcul réitère, le cas échéant, les deux étapes précédentes un nombre déterminé de fois pour d'autres impacts de tronçons aval (52 - 55) de trajectoire, - d'après les conditions de propagation sur la totalité de la trajectoire, le système de calcul détermine une atténuation cumulée pour tout point sélectionné de celle-ci, et - le système de calcul répète toutes les étapes ci-dessus une pluralité de fois pour une pluralité de directions initiales afin d'échantillonner toute la cellule (1) et ainsi déterminer une carte d'atténuations des points sélectionnés, 21 2828620 procédé caractérisé par le fait que, dans une phase initiale, - le système élabore la base de données (32, 33) sous la forme d'une matrice de données, à accès direct, de spécification locale et autonome des positions et natures des reliefs respectifs (41, 42) d'une pluralité de mailles prédéterminées d'un maillage correspondant de la carte (31), et, dans la phase d'exploitation, - compare la position du point courant au maillage pour identifier une maille d'impact, et - calcule les conditions de propagation après impact à partir des données de spécification locale et autonome de relief de la maille d'impact. 2. Procédé selon la revendication 1, dans lequel le système mémorise, dans la phase initiale, des données d'orientation géographique des reliefs dans la base de données matricielle (32, 33) et, en exploitation, calcule une direction d'un rayon réfléchi de tronçon aval (52) à partir des données d'orientation. 3. Procédé selon la revendication 2, dans lequel la mémorisation se limite à des données d'azimut et, en exploitation, la direction du rayon réfléchi de tronçon aval (82) est calculée en considérant que les reliefs de réflexion (71) sont verticaux. 4. Procédé selon l'une des revendications 1 à 3, dans lequel on intègre, dans la base de données matricielle (32, 33), des donnces spécifiant une nature de sursol et, en exploitation, le système calcule les conditions de propagation sur le tronçon aval (52) d'après ladite nature de sursol. 5. Procédé selon l'une des revendications 1 à 4, dans lequel, dans la phase initiale, on intègre des données d'arêtes (73) des reliefs dans la base de données matricielle (32, 33) et, en exploitation, le système de calcul calcule une direction de tronçon aval de rayon réfracté (92) à partir desUites données d'arêtes. 22 2828620 6. Procédé selon l'une des revendications 1 à 5, dans lequel, dans la phase initiale, des données d'atténuation des reliefs (41, 42) sont mémorisoes dans la base de donnces matricielle (32, 33) et, en exploitation, le système calcule une atténuation de propagation sur le tronçon aval à partir desdites données d'atténuation. 7. Procédé selon la revendication 6, dans lequel les données d'atténuation concernent la réflexion sur les reliefs (41, 42) et sont utilisées pour calculer l'atténuation sur le tronçon aval (52). 8. Procédé selon l'une des revendications 6 et 7, dans lequel les données d'atténuation concernent la propagation à travers les reliefs (41, 42) et sont utilisées pour calculer l'atténuation des rayons qui s'y propagent. 9. Procédé selon la revendication 8, dans lequel les données d'atténuation de propagation à travers les reliefs (41, 42) comportent en outre des données de transition entre milieux de propagation, spécifiant des atténuations de pénétration dans les reliefs (41, 42), qui sont utilisées pour déterminer une atténuation locale de pénétration. 10. Procédé selon l'une des revendications 1 à 9, dans lequel, disposant d'un algorithme de calcul de dispersion angulaire du rayon après impact, le système calcule une pluralité de directions de tronçons aval, avec des atténuations spécifiques, formant un angle solide (84) de diffusion du rayon au- delà de l'impact. 11. Procédé selon l'une des revendications 1 à 10, dans lequel le système compte les impacts successifs sur la trejectoire et, en cas de deuxième impact, le système considère, pour calculer lesdites conditions de propagation au-delà, que le rayon a été polarisé lors du premier impact. 12. Procédé selon l'une des revendications 1 à 11, dans lequel le système compte les impacts successifs sur la trejectoire et compare le total à une valeur de seuil haut pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint. 23 2828620 13. Procédé selon l'une des revendications 1 à 12, dans lequel, à chaque impact, le système détermine l'atténuation de propagation cumulée et la compare à une valeur de seuil d'atténuation maximale pour cesser d'exécuter les étapes du procédé lorsque le seuil est atteint. 14. Procédé selon l'une des revendications 1 à 13, dans lequel la carte des atténuations est mémorisée en trois dimensions. 15. Procédé selon l'une des revendications 1 à 14, dans lequel, pour élaborer la base de données (32, 33), le système représente la carte (31) par un faisceau de chanes de pixels à extension verticale et découpe les chanes de pixels pour constituer des volumes élémentaires empilés de mailles volumiques ayant chacun des données particulières. 16. Procédé selon l'une des revendications 1 à 15, dans lequel le système prend, comme position d'origine, une position prévue pour la station (11). 17. Procédé selon l'une des revendications 1 à 15, dans lequel le système prend, comme position d'origine, une position quelconque dans la cellule et la direction de lancement du rayon est choisie, d'après les positions et natures des reliefs proches, pour que le rayon passe à proximité de la station (11). 18. Procédé selon l'une des revendications 1 à 17, dans lequel, ayant ainsi calculé les conditions de propagation dans une microcellule en contact avec une cellule, le système calcule des conditions de propagation dans la cellule et effectue ensuite un calcul de lissage des résultats des deux calcuis