FR 2828623 A1 20030214 FR 0110724 A 20010810 calcule les dits flux (42 - 45) à partir des densités spécifiques validées. 1 2828623 Procédé d'établissement d'une carte de couverture radio La présente invention concerne les réseaux de radiotéléphonie cellulaire, et plus particulièrement, lors de la définition d'un réseau prévu par un opérateur, I'estimation prévisionnelle de la couverture radio de chacune des cellules, c'est-à-dire 1'effficacité de 1'offre de service de communication pour les terminaux mobiles. Un réseau de radiotéléphonie cellulaire est constitué d'une pluralité de stations radio terrestres de base qui assurent la couverture radio de zones définissant des cellules respectives. Pour pouvoir établir les communications o entre les terminaux radio mobiles d'une cellule et la station, les circuits récepteurs de ceux-ci doivent recevoir les signaux au-dessus de leur seuil de sensibilité. Le bruit électronique des composants des circuits récepteurs s'ajoute au bruit de propagation radio et peut rendre inutilisable la liaison, par s exemple, dans une liaison numérique, en provoquant une détection erronée de certains des bits transmis. Quel que soit le type de modulation utilisé, le risque d'erreur de détection des bits cro^'t rapidement dès que le niveau de signal utile reçu, donc le rapport signal à bruit, descend à une valeur de seuil n'offrant qu'une marge de sécurité minimale contre les impuisions de bruit. o Un codage d'émission de paquets de bits de données avec redondance, par code autocorrecteur, permet bien au récepteur de détecter des erreurs et de les corriger si leur nombre est limité. Cela permet de tolérer une baisse supplémentaire limitée du niveau de réception, mais, au-delà, le nombre d'erreurs excessif empêche toute correction. Même si la liaison s log ique avec la station peut i n itialement être établie et ma intenue pa r des signaux de service, pour réserver une voie de transmission utile, le taux excessif de bits utiles erronés rend inutilisable cette dernière, pour la voix numérisée et pour les données autres. D'une façon générale, il faut donc prévoir, pour les stations et les o terminaux, des circuits récepteurs à haute sensibilité, c'est-à-dire à seuil de détection relativement bas, et des circuits émetteurs puissants. 2 2828623 En réception, les composants des circuits utilisés sont sélectionnés pour ne présenter qu'un faible niveau de bruit intrinsèque. On peut donc fortement amplifier le signal utile reçu, jusqu'à ce que le niveau de bruit intrinsèque, amplifié avec le signal utile, commence à devenir un perturbateur. Le seuil de sensibilité du récepteur est alors atteint. Les distorsions de propagation radio constituent aussi une cause d'erreur, limitant les possibilités de détection des signaux faibles. En émission, la puissance des stations doit rester limitée, pour éviter les perturbations entre cellules ou autres systèmes radio. La puissance o d'émission des terminaux doit de même rester limitée, pour préserver l'autonomie de leur batterie et pour garantir la santé des utilisateurs. L'opérateur ne peut donc pas garantir que tous les terminaux d'une cellule vont recevoir un signal d'un niveau au moins égal au seuil de sensibilité pour assurer une liaison de radio, en particulier compte tenu des obstacles radio, tels que murs d'un appartement ou carrosserie d'une voiture transportant l'utilisateur. Il en est de même dans le sens montant, vers la station. En pratique, comme le niveau de champ radio décro'^t avec la distance selon une loi d'atténuation connue, la taille ou "rayon" des cellules est limité à la distance pour laquelle le niveau de champ est descendu sensiblement au seuil de sensibilité des terminaux, avec toutefois une marge déterminée de sécurité pour tenir compte des obstacles temporaires dans la propagation radio, provoquant un évanouissement du signal. L'opérateur du réseau prend aussi en compte les reliefs de la cellule et la nature du sursol, immeubles, forêts et autres, pour corriger les calcuis de niveau de champ effectués d'après la loi d'atténuation. En fonction de ces calcuis, il détermine ainsi des zones masquées dans la cellule, c'est-à-dire des zones dans lesquelles le niveau de réception des signaux de la station est inférieur au seuil de sensibilité des récepteurs. o En d'autres termes, et si chaque cellule est divisoe en un maillage d'une pluralité de pixels représentant chacun une surface de quelques dizaines ou centaines de mètres carrés, la carte de couverture radio est 3 2828623 constituée d'une mosafque de pixels "noirs" ou "blancs" selon que le champ de la station y est respectivement supérieur ou inférieur au seuil de bon fonctionnement des terminaux. Le problème de liaison montante vers la station est moins critique car celle-ci est équilibrée par le paramétrage des s conditions d'émission et la conception du réseau avec la liaison descendante. La demanderesse s'est posé le problème de la validité d'une telle approche binaire pour représenter, comme évoqué ci- dessus, ou autrement, la carte de couverture relative au service radiotéléphonique. o Elle a alors observé que le mode ci-dessus de détermination de l'existence de la couverture radio ne prenait pas en compte les aspects dynamiques de la propagation radio, c'est-à-dire l'évolution de celle-ci au cours du temps, et donc les alternances de possibilité de liaison radio et d'impossibilité dans la zone de chaque pixel, en fonction des fluctuations du s niveau reçu autour du seuil de sensibilité des récepteurs. La présente invention vise ainsi à proposer un procédé de détermination de la carte de couverture radio d'une cellule qui traduise mieux la réalité. A cet effet, I'invention concerne un procédé d'établissement d'une carte de couverture radio d'au moins une cellule d'un réseau de radiotéléphonie cellulaire, au moyen d'un système de calcul disposant de données préalablement mémorisées de carte géographique comportant la position de la cellule et d'une station radio de base associée, disposant de données préalablement mémorisées spécifiant un niveau d'émission de la s station et un niveau de seuil de sensibilité de réception radio de terminaux de la cellule, et disposant de données préalablement mémorisées d'une loi d'atténuation de propagation radio, procédé dans lequel: - le système de calcul effectue, à partir des données de carte, un maillage de la carte de la cellule en une pluralité déterminée de pixels de la carte occupant des positions dont il mémorise des données représentatives, - il calcule, par comparaison entre les données de carte, spécifiant la position de la station, et les données de position de chaque pixel et d'après 4 2828623 la loi d'atténuation, une pluralité d'atténuations de propagation radio du niveau d'émission de la station pour les pixels respectifs d'o il déduit, à partir du niveau d'émission de la station, une pluralité de niveaux moyens de réception, associés aux pixels respectifs, qu'il mémorise, et ii compare la pluralité de niveaux moyens de réception mémorisés au seuil de sensibilité pour déterminer une pluralité de sens d'inégalité entre la valeur du seuil et la valeur du niveau moyen de réception associé à chaque pixel, et en déduire la carte de couverture radio de la cellule, procédé caractérisé par le fait que o - des données représentant une loi de fluctuation temporelle de l'atténuation de propagation radio ayant été préalablement mémorisées et rendues accessibles au système de calcul, ce dernier applique la loi de fluctuation à la pluralité de niveaux moyens de réception, - il en détermine une pluralité correspondante de probabilités pour que s lesdites inégalités respectives aient un sens déterminé, et il établit la carte de couverture radio en associant, à chaque pixel, des données représentant ladite probabilité. Ainsi, et pour reprendre la formulation imagée des pixels noirs et blancs de l'art antérieur, I'invention propose des pixels à différents niveaux o de gris, c'est-à-dire une approche dynamique, qui fait dispara^'tre les distorsions qui existaient entre la réalité et l'estimation binaire, statique, de la couverture. Un pixel dont le niveau "statique" ou "moyen" de réception, calculé d'après la loi de propagation, est inférieur au seuil de sensibilité sera ainsi associé à une probabilité de possibilité de liaison, ou taux local de couverture, compris entre O et 50 %, ce dernier cas correspondant sensiblement au seuil. Un pixel blanc de l'art antérieur devient ainsi gris clair. Inversement, un pixel noir devient gris foncé, puisqu'il existe une probabilité de O à 50 % pour que le niveau de réception correspondant passe so temporairement sous le seuil. La couverture radio de chaque pixel, établie par la probabilité de dépassement ou non du seuil, fonction de la loi d'atténuation et de la loi de s 2828623 fluctuation du niveau de réception, traduit ainsi avec une bonne précision l'influence des phénomènes affectant la propagation radio. L'inégalité considérée, et la probabilité associée, peut indiquer que le niveau reçu dépasse le seuil, et alors elle exprime le fait qu'il y a couverture s radio, ou bien elle traite le cas opposé, en indiquant une non-couverture radio. Les deux cas correspondent donc à des calcuis se référant à des axes de pourcentages de couverture ou non-couverture de directions opposées, la couverture radio étant, dans les deux cas, parfaitement exprimée, directement par sa valeur ou bien son complément. Le cas de la parfaite o égalité des niveaux comparés est un cas particulier sans signification particulière, qui peut être assimilé à l'un ou bien l'autre des cas d'inégalité. On conçoit que le procédé de l'invention peut s'appliquer à plusieurs cellules et, globalement, à toute région géographique ou administrative. De préférence, au moins une autre cellule radio, avec une autre station de base, occupant une position lui permettant de recouvrir une zone de bord de la cellule considérée, le système de calcul détermine aussi la carte de couverture de l'autre cellule et il affecte, aux pixels de la zone de bord, une probabilité composite d'inégalité représentant la composition de la probabilité associée à l'émission de la station de la cellule considérée et une o probabilité homologue associée à l'émission de l'autre station. Ainsi, la modélisation du réseau est améliorée par prise en compte d'un effet d'entraide entre cellules, si bien que le niveau estimé de couverture moyen des cellules est accru. Comme le niveau d'émission de la station est en pratique déterminé par la distance aux zones de bord, le fait d'accro^'tre, s par composition, leur probabilité de couverture évite ainsi la nécessité de remonter le niveau d'émission de la station, qui ne présente par ailleurs aucun intérêt pour les zones proches de la station. On peut donc, le cas échéant, réduire le niveau prévu d'émission de la station, et donc les perturbations vers les autres cellules, si la valeur globale visée de couverture so radio de la cellule est alors dépassée, ou bien augmenter la taille de la cellule. 6 2828623 Avantageusement, la carte comportant des données de densité de population, le système de calcul en détermine la densité de population de chaque pixel et il la pondère par la probabilité associée à chaque pixel. Ainsi, la carte de couverture traduit une réalité fonctionnelle et commerciale, c'est-à-dire l'estimation du pourcentage de clientèle disposant du service radio à un instant déterminé, ou encore le taux de disponibilité fonction nel le des terminaux. On a u ra toutefo is comp ris q u 'il s'ag it de probabilités moyennes s'appliquant à l'ensemble de la population de terminaux, et que la probabilité de disponibilité d'un terminal déterminé o dépend de sa position dans la cellule, et en particulier est a priori plus faible au bord de celle-ci. Dans un but de vérification d'un cahier des charges, par exemple, le système de calcul peut calculer la somme des probabilités de la pluralité de pixels pour déterminer un taux moyen de couverture radio de la cellule. Le calcul peut porter sur la carte de couverture purement géographique ou sur la carte de couverture de population, donc sur les probabilités pondérées. L'invention sera mieux comprise à l'aide de la description suivante d'un mode préféré de mise en _uvre du procédé de l'invention en vue de définir les cellules d'un réseau de radiotéléphonie cellulaire, en référence aux o dessins annexés, dans lesquels: - la figure 1 représente schématiquement une carte géographique portant une partie d'un réseau de radiotéléphonie cellulaire en cours d'élaboration, - la figure 2 représente une gaussienne illustrant une loi de variation du niveau d'un signal radio reçu, - la figure 3 représente des fluctuations de signal reçu illustrant la gaussienne de la figure 2, et - la figure 4 illustre la pondération de la carte de couverture géographique par une carte de densité de population. o La figure 1 représente, sur un dispositif d'affichage d'un système de calcul 30, une carte 31 d'une zone géographique d'un réseau de radiotéléphonie cellulaire en cours de définition, sur laquelle ont été 7 2828623 reportées les positions prévues de quatre cellules radio, respectivement référencées 1, 2, 3 et 4, entourant des stations de base respectives 11, 12, 13 et 14 de positions indiquées. Comme représenté, les cellules 1 à 4 débordent les unes sur les autres, la cellule 1 ayant des zones de bord 120, 130 et 140 respectivement communes avec les cellules 2, 3 et 4. Comme on le comprendra, le recouvrement entre cellules 1 à 4 n'est qu'un cas particulier dans le cadre de l'invention. La carte 31 est une représentation de données d'une carte numérique (32) stockée dans une mémoire 32 d'une base de données 320 contenant o diverses mémoires 32, 33 et 35 à 38 du système de calcul ou calculateur 30, dont les moyens de traitement de calcul, ou unité centrale, de mémorisation, de saisie et de présentation de données, par un clavier eVou un écran d'affichage interactif avec souris ou autre, n'ont pas été représentés en détail. La carte numérique (32) spécifie les reliefs naturels et artificiels, ainsi que leur nature, telle que forêt, immeubles et autres, ce qui permet de calculer une estimation de l'atténuation radio de liaisons affectées par les reliefs. L'opérateur du réseau à implanter dispose d'une plage de fréquences et y a défini, et mémorisé dans le calculateur 30, une pluralité de valeurs de fréquences de porteuses constituant des canaux de transmission, numériques dans cet exemple, tels que ceux du réseau GSM ou UMTS, qui ont respectivement été affectées aux stations 11 à 14. Les valeurs des fréquences affectées ont été mémorisées dans la base de données 320, en association avec les stations concernées. Comme l'illustrent les flèches 41 et 42, un terminal radio mobile 21 situé dans la zone de bord avec recouvrement 120 peut recevoir les émissions de plusieurs stations, ici les deux stations 11 et 12 dans l'exemple du dessin. Parmi les deux signaux reçus, ou plus, le signal de niveau le plus élevé est considéré comme étant le signal utile, par exemple celui de la o station 11, et l'autre, provenant donc de la station 12, est un signal parasite, qui peut toutefois devenir le signal utile si les conditions de propagation radio inversent l'inégalité des niveaux. Accessoirement, il est en outre prévu une 8 2828623 gestion avec hystérésis de niveau pour limiter le nombre des commutations de gestion (hand-over en anglais) de la mobilité des terminaux 21 entre cellules 1 à 4. Ainsi, un signal de plus fort niveau perdant sa primauté de niveau reste toutefois considéré comme étant le signal utile tant que l'écart s avec le signal devenu le plus fort ne dépasse pas une valeur déterminée. Ainsi, d'une façon générale, les limites des cellules 1 à 4 fluctuent avec les conditions de propagation radio. Un bilan de la liaison radio correspondante peut être estimé d'après les positions relatives des stations 11 à 14 par rapport à la position de o terminal mobile 21 qui a été choisie par le calculateur 30, comme position du terminal mobile 21, c'est-à-dire un calcul de l'atténuation de transmission. Cette atténuation dépend en particulier de la distance entre l'émetteur considéré et le terminal récepteur 21 et peut facilement être calculée pour une propagation en espace libre dans l'air, à partir d'une valeur d'atténuation s linéique qui vaut environ 20 dB/km pour la plage de fréquences considérée et pour les 500 premiers mètres, et 30 dB/km audelà, ou même 40 dB/km en milieu urbain. Le cas échéant, en cas d'obstacles radio, les renseignements sur la nature des reliefs fournis par la base de données 32 permettent d'estimer le bilan de liaison sur une trajectoire de propagation déterminée. o La carte 31 est maillée en une pluralité de mailles ou pixels adjacents 34, représentant chacun ici un carré de 5 mètres de côté, la cellule 1 comportant N1 mailles. La carte renseignée 31 permet donc au calculateur de calculer un niveau prévisionnel de réception S1m des signaux des stations 11 à 14 en tout point ou pixel 34, si l'on a préalablement fixé leur s niveau d'émission, et il peut alors mémoriser la pluralité de niveaux de réception prévisionnels, associés à la pluralité de pixels 34, dans une mémoire 38 de la base de données 320. La référence 30 représente, de façon générale, un ensemble de traitement de données dont l'unité centrale a accès aux données nécessaires, contenues dans la base de données 320 (données (32) de carte géographique, lois d'atténuation et autres), ces données ayant été mémorisées dans des mémoires internes au calculateur 30 lui-même ou 9 2828623 étant disponibles dans des mémoires externes, par exemple des bases de données accessibles à travers un réseau de transmission de données. Ayant choisi a priori les positions des stations 11 à 14 et leurs puissances d'émission, le tout mis en mémoire 35, ainsi que la loi d'atténuation de propagation radio, mise en mémoire 33, et les données correspondantes étant accessibles au calculateur 30, ce dernier peut ainsi déterminer le niveau de réception S1m des signaux de la station 11 en tout point ou pixel 34 de la carte 31, dans la cellule 1, par exemple. Le niveau S1 m est en fait un niveau moyen de réception, comme expliqué plus loin. Le o calculateur 30 dispose aussi, dans une mémoire 36, d'une valeur de seuil de sensibilité de réception SO du terminal 21, c'est-à- dire du niveau minimal pour lequel le taux de bits erronés atteint un maximum admissible. Le calculateur 30 peut ainsi déterminer l'écart du niveau de réception S1m (niveau d'émission de la station 11 diminué de la valeur du bilan de liaison) s par rapport au seuil de sensibilité SO (fig. 2). Si l'écart S1m - SO est positif, la réception sera assurée une majorité du temps et, sinon, elle ne sera au mieux qu'occasionnelle. Comme évoqué au début, la propagation radio est tributaire d'aléas qui modulent le bilan de liaison au cours du temps t. Ainsi, un obstacle o mobile peut s'interposer sur le trajet des ondes, et les atténuer, ou au contraire produire une réflexion renforçant le signal reçu par le terminal 21, en " éclairant " par exemple une zone d'ombre radio. De telles fluctuations ont été étudiées et leurs caractéristiques statistiques ont été déterminées. La figure 2 représente une gaussienne illustrant une loi 370 de s fluctuation statistique du niveau de réception radio. La gaussienne représente les diverses probabilités P respectivement associées à divers niveaux de réception S1i, portés en abscisse S1, qui s'écartent du niveau moyen S1m. En pratique, il a été observé qu'il s'agissait d'une loi log- normale. A titre indicatif, la somme des probabilités associées aux niveaux so s'écartant de plus de 1 écart-type, EC, d'un côté vaut environ 16%. En d'autres termes, la probabilité est de 100 - 2 x 16 = 68 % pour qu'un niveau instantané ne soit pas écarté de plus de 1 écart-type EC par rapport à la 2828623 valeur moyenne S1m calculée par le bilan de liaison, sensiblement la plus probable. On notera toutefois que l'invention n'est pas limitée à cette loi particulière, susceptible d'étre modifiée par des conditions particulières, par exemple des passages répétés de mobiles perturbateurs, obstacles ou réflecteurs. D'une façon générale, il suffit que le calculateur 30 dispose en mémoire 37 de la loi de variation ou fluctuation statistique 370 pour déterminer, d'après le niveau de réception moyen S1m, la probabilité pour qu'un niveau instantané S1 i soit d'un côté ou de l'autre du seuil de sensibilité S0. La courbe de la figure 2 représente donc la forme de gaussienne de la loi o 370 considérée, forme qui est normée par le niveau de réception moyen calculé S1 m. Sur la figure 2, on a reporté le seuil de sensibilité S0 du terminal 21, de valeur fixe, et on a (fictivement) fait glisser la gaussienne selon l'axe des abscisses, pour positionner l'abscisse S1 m de son sommet au niveau s déterminé d'après le bilan de liaison du point considéré dans la cellule 1. La gaussienne 370 est ainsi coupée en deux parties, inférieure et supérieure, qui représentent chacune la probabilité pour que le signal instantané soit inférieur, respectivement supérieur, au seuil de sensibilité S0. Ainsi, si le niveau moyen reçu S1 m (sommet de la gaussienne) est. comme ici, o supérieur d'au moins 1 écart-type EC au seuil S0, la probabilité de couverture radio est d'au moins 84 %, c'est-àdire que le risque d'un évanouissement partiel d'amplitude négative dépassant 1 écart-type EC est inférieur à 16 % à tout instant quelconque. Comme le montre la concavité, tournée vers le haut, de la branche inférieure de la gaussienne, la probabilité de non-couverture décro^'t rapidement au-delà de 1 écart-type EC. Il est aussi représenté le cas opposé, pour lequel le niveau moyen S1m est supposé être inférieur au seuil de sensibilité du terminal 21, alors référencé S'0. Pour la commodité du dessin, on a en fait laissé fixe le niveau moyen S1 m et on a remonté le seuil S0 jusqu'au niveau S'0. Ces deux cas so d'inégalité entre le niveau moyen S1m et le seuil S0 ou S'0 sont illustrés par la figure 3. 11 2828623 La figure 3 est un exemple de fluctuations du niveau reçu S1. Comme le suggèrent, très schématiquement, les deux flèches partant respectivement des deux zones de gaussienne situées hors du "c_ur" de la gaussienne, allant de S1m - EC à S1m + EC dans la figure 2, le signal S1 reste statistiquement 16 % du temps t au-dessous du niveau S1m - EC et 16 % du temps au-dessus de S1m + EC. On voit donc que, pour un signal moyen reçu S1m supérieur au seuil S0 d'au moins 1 écart-type EC, la probabilité, pour qu'une impuision parasite négative abaisse le signal reçu instantané sous le seuil S0, est inférieure à 16 %. Symétriquement, si le niveau moyen o reçu S1 m est inférieur au seuil S'0 d'au moins 1 écart-type EC, la probabilité, pour qu'une impuision parasite positive translate le signal reçu instantané au- dessous du seuil S'0, est inférieure à 16 %. Comme expliqué plus haut, on a, pour la facilité du dessin, supposé deux niveaux de seuil S0, S'0 en laissant inchangé le niveau moyen reçu S1m, alors qu'en fait c'est ce dernier qui s varie avec l'atténuation, le seuil de sensibilité S0, S'0 étant fixe. Un exposé détaillé des calcuis de probabilité de couverture radio selon la marge de sécurité du niveau moyen de réception S1m se trouve dans l'article " Microwave Mobile Communication ", de Jakes Jr, John Wiley & Sons, New-York, 1974. o Le calculateur 30 détermine ainsi, pour chacun de la pluralité N1 d'une matrice de points ou pixels 34 de la carte 31 représentant toute la cellule 1, le niveau moyen S1m de réception des signaux de la station 11 et les fluctuations du niveau moyen S1m d'après la loi de variation statistique de la mémoire 37, d'o il détermine, par comparaison avec le seuil de sensibilité SO, la probabilité pour que le niveau instantané S1 i soit d'un côté, déterminé, du seuil S0. Au-dessus de ce dernier, à un instant quelconque, il y a couverture radio à l'instant considéré, et, audessous, le service radio n'est pas assuré. Comme déjà indiqué, le cas particulier d'égalité, très ponctuel, n'a pas de signification spécifique et peut être rattaché à l'un ou l'autre des o deux cas principaux. Ainsi, selon le procédé de l'invention, d'établissement d'une carte de couverture radio d'une cellule 1 du réseau de radiotéléphonie cellulaire, au 12 2828623 moyen du calculateur 30 disposant des données (32) préalablement mémorisées définissant la carte géographique 31 comportant la position de la cellule 1 et d'une station radio de base associée 11, disposant des données (35, 36) préalablement mémorisées spécifiant le niveau d'émission de la station 11 et le niveau de seuil de sensibilité de réception radio SO de terminaux 21 de la cellule 1, et disposant des données préalablement mémorisées d'une loi d'atténuation de propagation radio (33) , - le calculateur 30 effectue, à partir des données de carte (32), un maillage de la carte 31 de la cellule 1 en la pluralité déterminée N1 de pixels o 34 de la carte 31 occupant des positions dont il mémorise des données représentatives, - il calcule, par comparaison entre les données de carte (32), spécifiant la position de la station 11, et les données de position de chaque pixel 34 et d'après la loi d'atténuation (33), une pluralité N1 d'atténuations de propagation radio du niveau d'émission (35) de la station 11 pour les pixels respectifs 34, d'o il déduit, à partir du niveau d'émission (35) de la station 11, une pluralité de niveaux moyens de réception S1m, associés aux pixels 34 respectifs, niveaux S1 m qu'il mémorise en mémoire 38, et - il compare la pluralité de niveaux moyens de réception mémorisés o S1m au seuil de sensibilité SO des terminaux 21 pour déterminer une pluralité N1 de sens d'inégalité entre la valeur du seuil SO et la valeur du niveau moyen de réception S1m associé à chaque pixel 34, et en déduire la carte de couverture radio de la cellule 1. En outre, - d es don n ées (37) représentant u ne loi de fluctuation temporel le de l'atténuation de propagation radio ayant été préalablement mémorisées et rendues accessibles au calculateur 30, ce dernier applique la loi de fluctuation à la pluralité N1 de niveaux moyens de réception S1 m, - il en détermine une pluralité correspondante N1 de probabilités P so (311, fig. 4) pour que lesdites inégalités respectives aient un sens déterminé, et 13 2828623 - il établit la carte de couverture radio en associant, à chaque pixel 34, des donnces représentant ladite probabilité P (311). On notera que i'appellation " carte de couverture radio " désigne fondamentalement des données 311 indiquant la probabilité de couverture radio de chaque pixel 34. Ces données peuvent, si on le souhaite, être représentées sous la forme d'une carte, et/ou être exploitées numériquement pour calculer un taux moyen de couverture dans la cellule 1, par sommation des taux de couvertures ponctuels des pixels respectifs 34,le résultat étant normé par division par N1. On remarquera en outre que les calcuis ci-dessus o de niveau absolu reçu S1 sont équivalents à des calcuis de niveau relatif. En d'autres termes, le niveau de sensibilité S0 pourrait être défini par son écart (atténuation maximale tolérable) par rapport au niveau d'émission de la station 1 1. Des calcuis classiques, à partir d'une formule présentée par Jakes pour relier la probabilité ponctuelle, de couverture en bordure, à la probabilité intégrée globale dans la cellule 1, montrent que, si l'écarttype EC est par exemple de 10 (ou 8) dB, et si l'on considère que la limite de la cellule 1 est définie par un taux de couverture de 50 % (S1m = S0), la probabilité de couverture de la cellule 1 (supposée homogène en propagation radio) est de o 75 (ou 78) %. Inversement, si l'on part d'un taux de couverture radio global moyen visé, par exemple 95 % dans la cellule 1, il faut alors limiter sa taille, ou augmenter la puissance de la station 11, pour que la probabilité de couverture en limite de cellule 1 atteigne 85 (ou 80) %. La comparaison des deux exemples ci-dessus montre qu'un écart type EC réduit est favorable, en ce sens qu'une probabilité de 80 % en bordure suffit pour satisfaire à l'objectif global moyen de 95 %. Dans le cas de l'écart-type EC de 10 dB, la valeur de probabilité de 85 % en bordure conduit à y prévoir une marge (S1 m -S0)de 11,3dB. On notera que le seuil S0 peut varier avec le service radio offert, en o particulier un débit accru peut remonter ce seuil. Par exemple, dans le système de troisième génération UMTS, pour des paquets à un débit moyen 14 2828623 utile de 384 kb/s, un niveau en bordure S1m de -86 dBm a été retenu, incluant une marge de 12 dB. Le seuil S0 est donc de -98 dBm. Le procédé ci-dessus est naturellement applicable aux autres cellules 2 à 4, de façon mutuellement indépendante. Toutefois, en plus, il est ici s prévu une amélioration des résultats du procédé, consistant à prendre en compte un effet de bord, c'est-à-dire le fait qu'il y a une entraide entre stations 11 à 14. En effet, si l'on suppose que le terminal 21 est situé dans la zone de bord 120 de la cellule 1 recouverte par, ou simplement voisine de, une zone de bord d'une cellule adjacente 2, il existe une probabilité pour to qu'un niveau instantané S2 de réception des signaux de la station 12 par le terminal 21 soit supérieur au seuil S0. Le terminal 21 pourra alors se rattacher à la station 1 2 si sa liaison avec la station 11 présente un affaiblissement excessif. On conçoit cependant que la probabilité du secours est limitée, sinon faible, puisqu'il s'agit de liaisons sensiblement en limite de s portée radio, donc avec un sommet de gaussienne dont la position du niveau associé de réception moyen S1m, relativement proche du seuil S0, n'offre pas qu'une faible marge de sécurité. En pareil cas, la probabilité composite Pc de couverture radio du pixel 34 considéré est égale à Pc=A+ B-A.B, c'est-à-dire la somme des deux probabilités de couverture du pixel 34 considéré, A et B (fig. 1), associées aux stations respectives 11 et 12, diminuée de la valeur du produit A.B représentant le fait que l'une des deux probabilités B (ou A) est inutile lorsque l'autre probabilité A (ou B) est s effective, c'est-à-dire lorsque la liaison avec la station 11 (ou 12) fonctionne. A titre d'illustration, sur la figure 1, le pixel 34 repéré précisément par la référence 120 comporte un côté représenté renforcé en noir, intérieurement au pixel, sur toute sa longueur par une barre horizontale dont l'épaisseur est proportionnelle à la probabilité A, et, de même, un côté adjacent o perpendiculaire renforcé en noir, intérieurement au pixel, sur toute sa longueur par une barre verticale dont l'épaisseur est proportionnelle à la probabilité B. La probabilité composite Pc est représentée par la surface 2828623 totale noircie, qui est donc la somme des deux surfaces de trait épais, moins leur surface de recouvrement. On peut prendre en compte encore d'autres stations voisines 13, 14 pour remonter la valeur calculée du taux de couverture du pixel de zone de s bord. Ainsi, au moins une autre cellule radio 2, 3, 4, avec station de base 12, 13, 14, occupant une position lui permettant de recouvrir une zone de bord 120, 130, 140 de la cellule considérée 1, le calculateur 30 détermine aussi la carte de couverture de chaque autre cellule 2, 3, 4 et il affecte, aux opixels 34 de la zone de bord respective 120, 130 et 140, la probabilité composite Pc d'inégalité représentant la composition de la probabilité associée à l'émission de la station 11 de la cellule considérée 1 et une probabilité homologue associée à la station 12, 13, 14 de l'autre cellule 2, 3, 4. La figure 4 illustre un autre complément du procédé de l'invention, permettant de fournir un taux de couverture de la population prévue de terminaux 21 dans la cellule 1 considérée. La référence 311 sur la carte 31 représente les N 1 d iverses probab il ités de couvertu re d es pixels 34 respectifs, c'est-à-dire la carte géographique de couverture radio, ici o restreinte à quatre pixels 34 pour la clarté de l'exposé. La référence 312, sur la même carte 31, à nouveau représentée, désigne les N1 diverses densités prévues de terminaux 21, ou encore la densité de population, éventuellement modulée pour tenir compte du taux d'équipement en terminaux 21 dans la région de la cellule 1 considérée. La rétérence 313 représente, encore sur la s carte 31, les N1 divers produits de la densité statistique ou probabilité de couverture radio locale 311 par la densité 312 des terminaux 21 dans le pixel 34 considéré. Dans le calculateur 30, chaque pixel 34 est associé à une information, représentée ici par des rayures plus ou moins serrées, ou couleur, niveau de gris, rayures ou autres, représentant le produit résultant o de densité composite. Cette densité composite peut être exprimée par un nombre exprimant effectivement une densité relative entre pixels 34 ou par un nombre absolu si les pixels 34 sont de différentes tailles. On aura compris 16 2828623 que les données correspondantes de carte 312 sont accessibles au calculateur 30, y ayant de préférence été préalablement stockées, et que ce dernier effectue les calcuis voulus de multiplication ou pondération pour chaque pixel 34 et ensuite, le cas échéant, d'intégration des N1 résultats de s la carte 313 pour fournir une valeur intégrée de couverture prévisionnelle de la population estimée des utilisateurs, donc des terminaux 21. Ainsi, la carte 31 comportant des données (312) de densité de population, le calculateur 30 en détermine la densité de population de chaque pixel 34 et il pondère le taux de couverture radio des pixels 34 par o leur densité de population respective, pour obtenir la carte 313 de couverture radio de la population. Le lecteur comprendra aisément que le procédé de l'invention s'applique tout aussi bien à un ensemble de cellules d'une zone donnée (cette zone étant définie par un contour fictif [arbitraire ou simplement s commercial, administratif, etc., contour d'un département par exemple]). Il doit être évident, pour les personnes versées dans l'art, que la présente invention permet d es mod es de réal isation sous de nombreuses autres formes spécifiques sans sortir du domaine d'application de l'invention comme revendiqué. Par conséquent, les présents modes de réalisation o doivent être considérés à titre d'illustration mais peuvent être modifiés dans le domaine défini par la portée des revendications jointes. 17 2828623 REVENDICATIONS 1. Procédé d'établissement d'une carte de couverture radio d'au moins une cellule (1) d'un réseau de radiotéléphonie cellulaire, au moyen d'un système de calcul (30) disposant de données (32) préalablement mémorisées de carte géographique (31) comportant la position de la cellule (1) et d'une station radio de base associée (11), disposant de données préalablement mémorisées spécifiant un niveau d'émission (35) de la station (11) et un niveau de seuil de sensibilité (36) de réception radio (S0) de terminaux (21) de la cellule (1), et disposant de données préalablement o mémorisées d'une loi d'atténuation de propagation radio (33), procédé dans lequel: - le système de calcul (30) effectue, à partir des données de carte (32), un maillage de la carte (31) de la cellule (1) en une pluralité déterminée (N1) de pixels (34) de la carte (31) occupant des positions dont il mémorise des données représentatives, - il calcu le. par compa raison entre les don nées de carte (32), spécifiant la position de la station (11), et les données de position de chaque pixel (34) et d'après la loi d'atténuation (33), une pluralité (N1) d'atténuations de propagation radio du niveau d'émission (35) de la station (11) pour les o pixels respectifs (34), d'o il déduit, à partir du niveau d'émission (35) de la station (11), une pluralité de niveaux moyens de réception (S1m), associés aux pixels (34) respectifs, qu'il mémorise (38), et - il compare la pluralité de niveaux moyens de réception (S1 m) mémorisés au seuil de sensibilité (S0) pour déterminer une pluralité de sens d'inégalité entre la valeur du seuil (S0) et la valeur du niveau moyen de réception (S1m) associé à chaque pixel (34), et en déduire la carte de couverture radio de la cellule (1). procédé caractérisé par le fait que 18 2828623 - des données (37) représentant une loi de fluctuation temporelle (370) de l'atténuation de propagation radio ayant été rendues accessibles au système de calcul (30), ce dernier applique la loi de fluctuation (370) à la pluralité (N1) de niveaux moyens de réception (S1 m), - ii en détermine une pluralité correspondante (N1) de probabilités (P. 311) pour que lesdites inégalités respectives aient un sens déterminé, et - il établit la carte de couverture radio en associant, à chaque pixel (34), des données représentant ladite probabilité (P. 311). 2. Procédé selon la revendication 1, dans lequel, au moins une autre o cellule radio (2, 3, 4), avec une autre station de base (12, 13, 14), occupant une position lui permettant de recouvrir une zone de bord (120, 130, 140) de la cellule considérée (1), le système de calcul (30) détermine aussi la carte de couverture de chaque autre cellule (2, 3, 4) et il affecte, aux pixels (34) de la zone de bord respective (120, 130, 140), une probabilité composite d'inégalité (Pc) représentant la composition de la probabilité (A) associée à l'émission de la station (11) de la cellule considérée (1) et une probabilité homologue (B) associée à l'émission de l'autre station (12,13,14). 3. Procédé selon l'une des revendications 1 et 2, dans lequel, la carte (31) comportant des données de densité de population, le système de calcul o (30) en détermine la densité de population (312) de chaque pixel (34) et il la pondère (313) par la probabilité associée à chaque pixel (34). 4. Procédé selon l'une des revendications 1 à 3, dans lequel le système de calcul (30) calcule la somme des probabilités (P. 311, 313) de la pluralité de pixels (34) pour déterminer un taux moyen de couverture radio de la cellule (1). 5. Procédé selon la revendication 4, dans lequel le calcul porte sur ia