FR 2788132 A1 20000707 FR 9816730 A 19981231 1 2788132 Détermination de couverture radioélectrique d'émetteur La présente invention concerne la détermination de la couverture radioélectrique autour d'un émetteur à radiofréquence dans le cadre de la planification et la simulation de réseaux de radiocommunication, de radiodiffusion et de télédiffusion. Une carte géographique sur laquelle est dessiné au moins le contour de la couverture radioélectrique est à visualiser sur un écran ou est à éditer. La détermination de la couverture radioélectrique autour d'un émetteur fait appel aux valeurs de champ radioélectrique en des points de mesure réels sur le terrain autour de l'émetteur et aux caractéristiques du site et du signal émis de l'émetteur (Cf. Article "CECILIA WEB, A WWW INTERFACE FOR THE CECILIA APPLICATION" de G. BOURLIER et al., IEEE 47th Vehicular Technology Conference Proceedings, 1997). Dans une région de traitement prédéterminée contenant les points de mesure, le champ radioélectrique est déterminé en un point géographique donné en fonction d'un algorithme basé sur des valeurs du champ mesuré sur le terrain, de la distance de ce point avec les points de mesure réels et d'un modèle de propagation du signal émis depuis l'émetteur en dépendance du relief autour de celui- ci. L'algorithme est d'autant plus efficace que la densité des points de mesure dans la région de traitement est élevée, et plus précisément que la densité des points de mesure autour du point géographique donné est élevée. 2 2788132 Lorsque les points de mesure sont peu nombreux dans une zone relativement étendue ou accidentée dans la région de traitement, des aberrations numériques sur les valeurs de champ radioélectrique dans la zone apparaissent. La présente invention vise à remédier aux aberrations précitées en créant des points de mesure "fictifs" à l'intérieur des zones relativement désertes de points de mesure réels. Par exemple dans des zones accidentées, montagneuses, les points de mesure fictifs augmentent la densité des points de mesure sur lesquels le champ radioélectrique en des points géographiques de la zone est déterminé selon l'invention. A cette fin, un procédé pour déterminer la couverture radioélectrique d'un émetteur sur une carte géographique, comprenant des mesures de champ radioélectrique émis par l'émetteur en des points de mesure réels, est caractérisé par les étapes suivantes pour chaque point géographique: déterminer une variation maximale d'une grandeur physique dans un rayon prédéterminé autour du point géographique, - rechercher un point de mesure le plus proche du point géographique et évaluer la distance entre ces deux points, - rechercher un intervalle de grandeur physique contenant la variation de grandeur physique maximale parmi des intervalles de grandeur physique associés à des distances prédéterminées pour refuser le point géographique comme point de mesure fictif lorsque ladite distance est inférieure à la distance associée audit un intervalle de grandeur physique, et 3 2788132 - admettre le point géographique comme point de mesure fictif lorsque la distance est supérieure à la distance associée audit un intervalle de grandeur physique; puis - par une étape de déduire le champ radioélectrique en chaque point géographique en fonction des champs radioélectriques aux points de mesure réels et fictifs et de distances entre ces derniers et le point géographique. La grandeur physique peut être l'altitude ou le champ radioélectrique. La variation de grandeur physique autour d'un point géographique donné simule le relief autour de ce point. Le point de mesure le plus proche, qui peut être réel ou fictif au fur et à mesure du traitement des points géographiques prédéterminés, simule la densité de points de mesure autour du point géographique donné. Plus le point de mesure le plus proche sera près du point géographique, et plus le relief autour du point géographique donné est accidenté pour que le point géographique donné ait des chances de devenir un point de mesure fictif. Ainsi les champs radioélectriques initialement prédits ne sont véritablement utilisés, selon le procédé de l'invention, que pour les points géographiques réputés ultérieurement points de mesure fictifs, lorsque la grandeur physique n'est pas le champ radioélectrique. Le champ radioélectrique en chaque point géographique est ensuite déduit, c'est-à-dire interpolé, non seulement en fonction des champs et distances relatifs aux points de mesure sur le terrain, mais également en fonction des champs et 4 2788132 distances relatifs à des points de mesure "fictifs" selon l'invention. Afin de ne pas produire des points de mesure fictifs présentant un champ radioélectrique faible qui n'est plus assez élevé pour des radiotéléphones mobiles d'un réseau de radiocommunication, ou bien pour des récepteurs d'un réseau de radiodiffusion ou de télédiffusion, l'étape d'admettre le point géographique comme point de mesure fictif lorsque la distance est supérieure à la distance associée audit un intervalle de grandeur physique peut comprendre: - lire le type réel ou fictif du point de mesure le plus proche, - prédire un champ radioélectrique d'émetteur au point géographique en fonction de caractéristiques de l'émetteur et sur la base d'un modèle de propagation prédéterminé - admettre le point géographique comme point de mesure fictif lorsque le point de mesure le plus proche est un point de mesure réel, - comparer les champs radioélectriques au point géographique et au point de mesure le plus proche à un champ radioélectrique minimal lorsque le point de mesure le plus proche est un point de mesure fictif, et - admettre le point géographique comme point de mesure fictif lorsque les champs radioélectriques au point géographique et au point de mesure fictif le plus proche sont supérieurs au champ radioélectrique minimal. Selon une autre variante, lorsque la grandeur physique est le champ radioélectrique, l'étape de prédire un champ radioélectrique d'émetteur au point géographique peut être réalisée avec l'étape 2788132 préalable de déterminer une variation maximale de champ radioélectrique. De préférence, les limites des intervalles de grandeur physique, altitude ou champ radioélectrique, consécutifs croissent lorsque les distances prédéterminées qui leur sont associées décroissent afin d'améliorer l'efficacité du procédé de l'invention. Le procédé de l'invention peut comprendre préalablement une détermination d'un pas en longitude et d'un pas en latitude entre les points géographiques o le champ radioélectrique est prédit. Le rayon prédéterminé autour du point géographique peut être sensiblement égal au demi-pas, afin de faciliter la mise en euvre partielle du procédé sous forme de logiciels et réduire la durée du procédé. Le procédé de l'invention peut comprendre préalablement une sélection du nombre de zones de couverture radioélectrique à visualiser et une sélection d'intervalles de variation de champ radioélectrique respectivement associés aux zones de couverture afin de mieux appréhender les variations de champ autour du site de l'émetteur. La durée de traitement des points géographiques prédéterminés o le champ radioélectrique est prédit, est diminuée lorsque les points géographiques prédéterminés sont successivement et périodiquement situés le long d'un chemin sensiblement en spirale depuis l'émetteur. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante de plusieurs réalisations préférées de l'invention en référence aux dessins annexés correspondants dans lesquels: 6 2788132 - la figure 1 est un graphe sur un écran d'ordinateur, montrant trois zones de couverture dans une région de traitement autour d'un émetteur; la figure 2 est un algorithme des étapes du procédé selon l'invention; la figure 3 est un graphe illustrant schématiquement la succession de points géographiques de la région de traitement montrée à la figure 1, auxquels le procédé selon l'invention est appliqué; et - la figure 4 est un algorithme détaillant des étapes incluses dans le procédé selon la figure 2 pour générer un point de mesure fictif, en fonction de variables associées à un point géographique courant. Le procédé de détermination de couverture radioélectrique d'un émetteur à radiofréquence S selon l'invention est principalement mis en oeuvre sous forme de logiciel. La configuration matérielle minimale requise est un ordinateur PC INTEL PENTIUM fonctionnant sous le système d'exploitation MICROSOFT WINDOWS NT muni de 32 Mo de mémoire RAM et un disque dur minimal de 1 Go. Le langage de programmation est C++. L'ordinateur PC est configuré dans un mode d'affichage possédant 256 couleurs afin que l'effet de relief d'une carte géographique soit visualisé. Selon une autre variante, le matériel est réparti entre un serveur avec des bases de données et un terminal d'usager reliés par un réseau de transmission numérique à débit élevé (INTERNET) et fonctionne avec un logiciel en langage de programmation JAVA, comme décrit dans l'article précité de G. BOURLIER et al. et illustré à la figure 1 de cet article. 7 2788132 L'ordinateur PC contient également des fichiers de données géographiques d'une carte. Par exemple, les données géographiques sont relatives à une carte de France comprenant les principaux axes AR et RI du réseau routier et du réseau hydrographique et les noms de ville, ainsi que le relief. Chaque point géographique élémentaire de la carte correspond à un pixel sur l'écran de l'ordinateur et est repéré par trois coordonnées X, Y et Z correspondant respectivement à la longitude, la latitude et l'altitude du point géographique. L'altitude Z du point géographique est déduite des deux autres coordonnées X et Y dans le plan de l'écran montré à la figure 1, par lecture d'un tableau mémorisé à deux entrées. Le procédé de détermination de couverture radioélectrique selon l'invention comprend essentiellement onze étapes E1 à E12, comme montré à la figure 2. A la première étape El, un émetteur à radiofréquence S dont la couverture est à déterminer est défini par des caractéristiques géographiques et techniques. Les caractéristiques géographiques de l'émetteur sont la longitude XS, la latitude YS et l'altitude ZS du site de l'émetteur S associées au nom de l'émetteur. Les caractéristiques techniques de l'émetteur S sont la fréquence du signal à radiofréquence émis par l'émetteur, par exemple la fréquence porteuse d'un canal de télévision, la polarisation et la puissance apparente rayonnée du signal émis par l'émetteur S, la hauteur de l'antenne de l'émetteur au-dessus du site de l'émetteur, ainsi que le diagramme azimutal d'émission de l'antenne. 8 2788132 Une campagne de mesures est effectuée sur le terrain à l'étape El pour mesurer le champ radioélectrique du signal reçu en différents points de mesure réels PMRj qui sont géographiquement éparpillés autour de l'émetteur S. L'indice entier j varie de 1 à un nombre élevé J supérieur à 30 par exemple. Chaque point de mesure réel PMRj est défini notamment par une longitude XMj et une latitude YM desquelles est déduite l'altitude ZMj du point PMRj. Le point de mesure PMRj est associé à la date de la mesure et au nom de la commune voisine, ainsi qu'à la valeur de champ radioélectrique EMj mesuré en ce point, exprimée en [dB(iV/m)]. Les caractéristiques de l'émetteur S et les caractéristiques des points de mesure PMR1 à PMRJ sont chargées dans l'ordinateur PC. Une région de traitement RT formant une partie de la carte préalablement enregistrée dans l'ordinateur est déterminée automatiquement en fonction de la puissance d'émission de l'émetteur S et est sensiblement centrée sur l'émetteur. La région de traitement RT englobe tous les points de mesure réels PMR1 à PMRJ. En variante, l'utilisateur peut définir la région de traitement à l'aide de la souris de l'ordinateur et d'un rectangle élastique déplaçable sur l'écran de l'ordinateur de manière que la région de traitement remplisse complètement l'écran. A titre d'exemple, la figure 1 montre une région de traitement rectangulaire RT dans laquelle sont indiqués le site de l'émetteur S et les points de mesure réels PMR1 à PMRJ représentés sous forme de carré avec une croix à l'intérieur, sur un fond géographique de région de traitement en relief contenant des axes routiers AR en traits pointillés 9 2788132 et des rivières RI en traits mixtes. En cliquant sur un pixel, c'est-àdire en sélectionnant un point géographique élémentaire du fond de carte, une première barre d'état inférieure BE1 fournit des informations sur le point de mesure le plus proche de ce point géographique pointé par la souris tels que commune la plus proche, coordonnées du point de mesure et champ radioélectrique du point de mesure le plus proche, et une deuxième barre d'état BE2 fournit des informations sur le point géographique élémentaire pointé par la souris, tels que la longitude X, la latitude Y et l'altitude Z de celui- ci. A une deuxième étape préalable E2, des paramètres de la zone de couverture à déterminer autour de l'émetteur S dans la région de traitement RT illustrée à la figure 1 sont sélectionnés par l'utilisateur au moyen d'une boîte de dialogue. Les paramètres de la zone de couverture sont, outre les caractéristiques précitées de l'émetteur S, un seuil minimal de visualisation de champ radioélectrique SEm exprimé en [dB(pV/m)] en deçà duquel la couverture ne sera pas visualisée, le nombre NZ de zones à visualiser par exemple au plus égal à 6, éventuellement la variation de champ radioélectrique pour chaque zone exprimée en [dB(tV/m)] ou bien le pas de champ radioélectrique entre deux zones consécutives, et le pas p de I points géographiques PGi d'une grille de modèle numérique du terrain contenu dans la région de traitement RT, et optionnellement un champ radioélectrique mininal prédéterminé Emin propre à l'invention, avec Emin 2788132 Le pas p qui définit la distance horizontale en longitude ou verticale en latitude entre deux points géographiques voisins PGi et PGi+1 o le champ radioélectrique EGi, EGi+1 est à déterminer. Le pas p des points géographiques à traiter aux étapes suivantes E3 à Ell est supérieur ou égal au pas des points géographiques élémentaires. Selon un premier exemple, deux zones de couverture sur l'écran de l'ordinateur PC sont à visualiser, la première couvrant des champs radioélectriques supérieurs à 60 [dB(pV/m)], et la deuxième couvrant des champs radioélectriques compris entre Emin = 54 et 60 [dB(ptV/m)]. Selon un deuxième exemple illustré à la figure 1, NZ = 3 zones de couverture Z1, Z2 et ZNZ = Z3 à visualiser sur l'écran de l'ordinateur PC de préférence avec des couleurs ou des nuances différentes sont sélectionnées avec un pas de visualisation de 10 [dB(iV/m)], soit une première zone Z1 autour de l'émetteur S pour des champs radioélectriques supérieurs à 55 [dB(pV/m)], une deuxième zone Z2 pour des champs radioélectriques compris entre 45 et 55 [dB(kV/m)], et une troisième zone Z3 souvent éloignée de l'émetteur S pour des champs compris entre Emin = 35 et 45 [dB(pV/m)]. Le pas p entre deux points géographiques voisins ayant une même longitude ou une même latitude correspond par exemple à 250 m, ou à 500 m ou 1000 m. Plus le pas p est grand, plus le temps pour déterminer la couverture radioélectrique de l'émetteur est faible, mais également la précision est plus faible. Après les étapes préliminaires El et E2 de chargement de caractéristiques de l'émetteur S et des ii 2788132 points de mesure PM1 à PMj et des paramètres de couverture, trois étapes E3, E4, et E5 sont d'abord effectuées suivant cet ordre, ou suivant un autre ordre, pour chaque point géographique PGi de la grille sélectionnée constituant un modèle numérique de la région de traitement RT autour de l'émetteur S, avant de décider à des étapes suivantes E6 à E10 si le point géographique PGi doit constituer un point de mesure fictif PMF selon l'invention. Les points géographiques PG1 à PGI de la grille sélectionnée sont traités successivement depuis le site de l'émetteur S en suivant un chemin en "spirale" vers les bords du fond de carte géographique de la zone de traitement RT, comme montré schématiquement par des flèches dans la figure 3. L'étape E3 prédit le champ radioélectrique EGi au point géographique courant PGi selon un modèle de propagation préalablement enregistré dans l'ordinateur PC, en fonction du relief indiqué par le fond de carte dans une zone sensiblement entre l'émetteur S, dont les caractéristiques de fonctionnement sont connues, et le point géographique courant PGi. Le champ prédit EGi est exprimé en [dB(tV/m)]. Le modèle de propagation peut être déterminé selon la Recommandation 370-5 de l'UIT, Genève. Une variation maximale VGPi d'une grandeur physique prédéterminée GP dans une zone circulaire ayant un rayon prédéterminé typiquement égal à p/2 autour du point PGi est déterminée à l'étape E4 pour caractériser le type de relief autour du point PGi. Selon une première réalisation, la grandeur physique prédéterminée GP est l'altitude A; la variation VGPi est alors une variation d'altitude 12 2788132 maximale VAi, c'est-à-dire la différence d'altitude entre le point le plus haut et le point le plus bas dans la zone circulaire précitée. Selon une deuxième réalisation, la grandeur physique prédéterminée GP est le champ radioélectrique E; la variation VGPi est alors une variation de champ électrique maximale VEi, c'est-à-dire la différence de champ radioélectrique prédit entre le point qui possède la meilleure réception théorique et le point qui possède la plus mauvaise réception théorique dans la zone circulaire précitée. Les champs radioélectriques aux points dans la zone circulaire sont prédits selon le modèle de propagation utilisé à l'étape E3. L'étape suivante E5 recherche le point de mesure le plus proche du point géographique courant PGi. Le type du point de mesure le plus proche, dénoté ci- après PMi, est réel si ce point est un point de mesure PMR1 à PMRj, ou est fictif si ce point est un point géographique d'indice inférieur à i et sélectionné comme point de mesure fictif PMF lors du déroulement des étapes E3 à E10 précédant celles du point géographique courant PGi. Le point de mesure le plus proche PMi est associé à une valeur de champ radioélectrique EMi qui est chargée à l'étape El si le point PMi est un point de mesure réel, ou qui est prédite à une étape précédente E3 si le point PMi est un point de mesure fictif. La distance Di = IIPGi, PMill entre le point géographique PGi et le point de mesure le plus proche PMi est calculée. En référence à la figure 4, les grandeurs physiques EGi, VGPi, EMi et Di associées au point géographique courant PGi et mémorisées aux étapes E3, E4 et E5 sont comparées respectivement à des seuils prédéterminés pour décider si le point géographique 13 2788132 PGi constitue ou ne constitue pas un point de mesure fictif selon l'invention. A l'étape E6, K intervalles de variation de grandeur physique, c'est-àdire d'altitude ou de champ radioélectrique prédit, consécutifs [GPo, GP1] à [GPK_1, GPK] sont associés respectivement à des distances D1 à DK. L'étape E6 recherche l'un [GPk, GPk+1] des intervalles précités, avec 0 radioélectrique, VGPi = VAi ou VGPi = VEi appartient. Sachant que les distances D1 à DK décroissent lorsque les limites GP0 à GPK des intervalles de variation de grandeur physique augmentent, plus le point de mesure le plus proche PMi est éloigné du point géographique PGi, plus la probabilité qu'il soit sélectionné est élevée lorsque la variation de grandeur physique maximale VGPi est faible, et a contrario, plus la variation de grandeur physique maximale VGPi est élevée, plus le point PMi peut être proche du point géographique PGi pour que ce dernier soit sélectionné. Les intervalles de variation de grandeur physique, altitude ou champ radioélectrique, associés respectivement à des distances Dk simulent le relief environnant le point géographique PGi: lorsque le relief est plat, en plaine, peu de points de mesure sont nécessaires, tandis qu'une densité de points de 14 2788132 mesure fictifs et/ou réels sont nécessaires lorsque le relief est accidenté, tel qu'en montagne. Pour la première réalisation avec l'altitude A comme grandeur physique GP, par exemple K = 6 intervalles de variation d'altitude sont sélectionnés [Ao, A1] = [ 0 m, 10 m] avec D1 = 20000 m, [A1, A2] = [ 10 m, 50 m] avec D2 = 12000 m, [A2, A3] = [ 50 m, 100 m] avec D3 = 7000 m, [A3, A4] = [100 m, 200 m] avec D4 = 4000 m, [A4, A5] = [200 m, 400 m] avec D5 = 3000 m, et [A5, A6] = [400 m, 800 m] avec D6 = 2000 m. Selon cet exemple, si VAi = 42 m, et si Di = 6450 m D2, le procédé passe à l'étape suivante E8. Pour la seconde réalisation avec le champ radioélectrique E comme grandeur physique GP, par exemple K = 4 intervalles de variation de champ radioélectrique sont sélectionnés: [E0, E1] = [ 0 dB(tV/m), 2 dB(gV/m)] avec D1 = 20000 m, [E1, E2] = [ 2 dB(gV/m), 5 dB(gV/m)] avec D2 = 10000 m, [E2, E3] = [ 5 dB(gV/m), 12 dB(jV/m)] avec D3= 5000 m, [E3, E4] = [12 dB(.V/m), 30 dB(HV/m)] avec D4= 2000 m, Selon cet exemple, si VCRi = 4 dB(tV/m), et si Di = 6450 m 2788132 le point géographique suivant PGi+1. Par contre si Di est égale par exemple à 10500 m > D2, le procédé passe à l'étape suivante E8. Lorsque la variation de grandeur physique maximale VGPi, VAi ou VEi et la distance Di pour le point géographique PGi satisfont les conditions des étapes E6 et E7, le type, réel ou fictif, du point de mesure le plus proche PMi est lu à l'étape E8. Si le point de mesure PMi est un point de mesure réel, c'est-à-dire ayant été défini lors de la campagne de mesures précédente, le point géographique PGi est réputé point de mesure fictif PMF selon l'invention à l'étape E10, et le procédé retourne à l'étape E3 pour le point géographique suivant PGi+1 si l'indice i est inférieur à I à l'étape Ell. Lorsqu'à l'étape E8 le point de mesure PMi est un point de mesure fictif, le point géographique courant PGi est réputé point de mesure fictif à l'étape E10 si à la fois, à l'étape E9, le champ radioélectrique prédit EGi au point géographique PGi est supérieur au champ radioélectrique minimal prédéterminé Emin et le champ radioélectrique EMi au point de mesure le plus proche PMi est supérieur au champ radioélectrique minimal Emin. Dans le cas contraire lorsque le champ radioélectrique EGi ou EMi est inférieur à Emin, le point géographique PGi n'est pas un point de mesure fictif, et le procédé retourne à l'étape E3 après incrémentation de l'indice i. Puis après l'étape E10, le procédé est rebouclé, comme déjà dit, si l'indice i est inférieur au nombre I de points géographiques de la grille de modèle numérique de terrain, de manière à incrémenter d'une unité l'indice i et à effectuer les étapes E3 à E11 pour chaque point géographique suivant de la grille dans la région de traitement RT balayée suivant la 16 2788132 spirale montrée à la figure 3. Chaque point géographique réputé point de mesure fictif PMF présente ainsi un champ radioélectrique prédit qui peut être supérieur au champ radioélectrique minimal Emin lorsque le point de mesure le plus proche PMi est fictif, ou qui peut être inférieur ou supérieur à Emin lorsque le point PMi est réel. Ainsi, pour un intervalle [GPk, GPk+1] contenant VGPi, tout point géographique PGi avec EGi > Emin est sélectionné comme point de mesure fictif PMF lorsque Di est supérieur à Dk, à l'exception du cas o le point de mesure le plus proche PMi est un point fictif et le champ EMi du point de mesure le plus proche ou le champ prédit EGi est déjà inférieur au champ radioélectrique minimal prédéterminé Emin. Ceci signifie que le point considéré PGi est déjà trop éloigné de la zone de couverture et qu'il est inutile de le sélectionner, ce qui limite le nombre de points de mesure fictifs générés et réduit le temps de déroulement du procédé. En variante, lorsque la grandeur physique GP est différente du champ radioélectrique prédit E, c'est- à-dire est l'altitude A, l'étape E3 pour prédire le champ radioélectrique au point géographique PGi montrée à la figure 2 est supprimée et remplacée par une étape E3' identique à l'étape E3 et introduite entre les étapes E8 et E9 lorsque le point PMi est fictif, et entre les étapes E8 et E10 lorsque le point PMi est réel, comme montré en traits pointillés à la figure 4. Les champs ne sont prédits que pour les points géographiques PGi ayant satisfait la condition Di > Dk à l'étape précédente E7. Que le point de mesure le plus proche PMi soit de type fictif ou de type réel à l'étape E8, la prédiction du 17 2788132 champ EGi est quand même effectuée puisque l'étape suivante E12 a besoin des champs prédits de tous les points de mesure fictifs. Cette variante réduit la durée des étapes E3 à E10 et par conséquent la durée du procédé. Lorsque la grandeur physique GP est égale au champ radioélectrique E, l'étape E3' n'existe pas et l'étape E3 est maintenue puisqu'à l'étape E4, le champ radioélectrique est obligatoirement prédit au point géographique EGi et aux points environnants. En se reportant à la figure 2, après que les étapes E3 à Ell aient été effectuées pour chaque point géographique PG1 à PGI de la grille dans la région de traitement RT, l'étape E12 suit un algorithme d'interpolation basé sur les points de mesure réels PMR1 à PMRj et les points de mesure fictifs PMF générés à l'étape précédente E10 selon l'invention, désignés indifféremment par PMm ci- après. L'interpolation consiste à définir une fonction s de R vers R qui satisfait aux conditions {s(XMm, YMm) = Em; avec 1 18 2788132 entre les points PG1 et PMm, et pm représente une pondération en fonction du point de mesure PMm; les points de mesure réels PMR1 à PMRJ ont un poids par exemple 10 fois supérieur à celui des points de mesure fictifs PMF. L'interpolation par fonction radiale et le type d'algorithme traité sont décrits dans l'article "RADIAL BASIS FUNCTIONS FOR MULTIVARIABLE INTERPOLATION: A REVIEW" de M.J.D. POWELL, pages 143 à 167, dans l'ouvrage "Algorithms for Approximation", édité par Mason & Cox, Clarendon Press Oxford, 1987. La fonction d'interpolation "passe" par les points de mesure réels, "passe" quasiment par les points de mesure fictifs et est calculée sur les autres points géographiques en fonction des distances qui existent avec les points de mesure et des valeurs de champ radioélectrique qui leur sont associées. En pratique, plusieurs campagnes de mesures peuvent compléter des premières déterminations de couverture radioélectrique d'un émetteur donné afin d'affiner la couverture radioélectrique. 19 2788132 REVENDICATIONS 1 - Procédé pour déterminer la couverture radioélectrique d'un émetteur (S) sur une carte géographique (RT), comprenant des mesures (El) de champ radioélectrique (EMj) émis par l'émetteur en des points de mesure réels (PMRj), caractérisé par les étapes suivantes pour chaque point géographique (PGi): - déterminer (E4) une variation maximale (VGPi) d'une grandeur physique (GP) dans un rayon prédéterminé autour du point géographique (PGi), - rechercher (E5) un point de mesure (PMi) le plus proche du point géographique (PGi) et évaluer la distance (Di) entre ces deux points (PMi, PGi), - rechercher (E6, E7) un intervalle de grandeur physique contenant la variation de grandeur physique maximale (VGPi) parmi des intervalles de grandeur physique (GPk, GPk+1) associés à des distances prédéterminées (Dk) pour refuser le point géographique (PGi) comme point de mesure fictif lorsque ladite distance (Di) est inférieure à la distance associée audit un intervalle de grandeur physique, et - admettre (E8, El0) le point géographique (PGi) comme point de mesure fictif (PMF) lorsque la distance (Di) est supérieure à la distance associée audit un intervalle de grandeur physique; puis - par une étape de déduire (El2) le champ radioélectrique en chaque point géographique (PGi) en fonction des champs radioélectriques aux points de mesure réels et fictifs (PMRj, PMF) et de distances entre ces derniers et le point géographique. 2788132 2 - Procédé conforme à la revendication 1, selon lequel la grandeur physique (GP) est l'altitude (A). 3 - Procédé conforme à la revendication 1, selon lequel la grandeur physique (GP) est le champ radioélectrique (E). 4 - Procédé conforme à l'une quelconque des revendications 1 à 4, selon lequel l'étape d'admettre le point géographique comme point de mesure fictif comprend - lire (E8) le type réel ou fictif du point de mesure le plus proche (PMi), - prédire (E3') un champ radioélectrique d'émetteur (EGi) au point géographique (PGi) en fonction de caractéristiques de l'émetteur (S) et sur la base d'un modèle de propagation prédéterminé, - admettre (E10) le point géographique (PGi) comme point de mesure fictif (PMF) lorsque le point de mesure le plus proche (PMi) est un point de mesure réel, - comparer (E9) les champs radioélectriques (EGi, EMi) au point géographique (PGi) et au point de mesure le plus proche (PMi) à un champ radioélectrique minimal (Emin) lorsque le point de mesure le plus proche (PMi) est un point de mesure fictif, - admettre (ElO) le point géographique (PGi) comme point de mesure fictif (PMF) lorsque les champs radioélectriques (EGi, EMi) au point géographique (PGi) et au point de mesure fictif le plus proche (PMi) sont supérieurs au champ radioélectrique minimal (Emin) 21 2788132 - Procédé conforme à la revendication 4, lorsqu'elle dépend de la revendication 3, selon lequel l'étape (E3) de prédire un champ radioélectrique d'émetteur (EGi) au point géographique (PGi) est réalisée avec l'étape (E4) de déterminer une variation maximale de champ radioélectrique (VEi). 6 - Procédé conforme à l'une quelconque des revendications 1 à 5, selon lequel les limites des intervalles de grandeur physique consécutifs (GPk, GPk+l) croissent lorsque les distances prédéterminées (Dk) qui leur sont associées décroissent. 7 - Procédé conforme à l'une quelconque des revendications 1 à 6, comprenant préalablement (E2) une détermination d'un pas (p) en longitude (X) et d'un pas (p) latitude (Y) entre les points géographiques (PGi) o le champ radioélectrique est prédit. 8 - Procédé conforme à la revendication 7, selon lequel le rayon prédéterminé autour du point géographique (PGi) est sensiblement égal au demi-pas (p/2). 9 - Procédé conforme à l'une quelconque des revendications 1 à 8, comprenant préalablement une sélection (E2) du nombre (NZ) de zones de couverture radioélectrique (Z1 à ZNZ) à visualiser et une sélection d'intervalles de variation de champ radioélectrique respectivement associés aux zones de couverture. 22 2788132 - Procédé conforme à l'une quelconque des revendications 1 à 9, selon lequel les points géographiques (PG1 à PGI) o le champ radioélectrique est prédit sont successivement et périodiquement situés le long d'un chemin sensiblement en spirale depuis l'émetteur (S).