FR 2735240 A1 19961213 FR 9506646 A 19950606 La présente invention est relative aux procédés et systèmes pour la localisation radio-satellitaire de points à la surface de la terre. L'utilisation des signaux radioélectriques transmis par les satellites a trouvé, depuis une dizaine d'années, un vaste champ d'applications dans les domaines de la navigation, de la géodésie et de la topographie. Particulièrement, le système GPS (Global Navigation System selon la terminologie anglosaxonne), mondialement utilisé, permet désormais d'accéder, en temps réel, à des précisions centimétriques pour la détermination de points à la surface du globe. Ces grandes précisions sont indispensables pour les topographes dans un grand nombre d'applications de positionnement terrestre. Les techniques mises en oeuvre pour parvenir à ce résultat sont dites de "cinématique temps réel". Elles sont utilisées depuis les années 1993-94 et illustrées en particulier dans le brevet français déposée par la demanderesse sous le n 94-00523. On pourra également, pour une présentation générale de ces techniques, se référer à la publication " GPS Vermessung mit Real-Time-Kinematic Technische Untersicht " - VERWENDUNGEN & TECHNISCHE ANGABEN, B 1-7, TRIMBLE-1993. A l'aide de ces techniques, il est possible, sous la seule contrainte de recevoir simultanément les signaux en provenance d'au moins cinq satellites, d'assurer un travail de qualité et d'atteindre des précisions de détermination de l'ordre du centimètre pour le lever ou l'implantation de points. Par conséquent, les techniques de radiolocalisation satellitaire sont, dans de nombreuses applications, substituées aux techniques de localisation antérieures qui utilisaient des appareillages optiques pour la mesure d'angles ou de distances (distancemètre et tachéomètre). Toutefois, certains points au sol ne sont pas toujours suffisamment dégagés ou convenablement situés pour permettre une réception correcte des signaux en provenance des satellites. Ceci constitue un frein au développement des techniques de radiolocalisation satellitaire et notamment à leur utilisation par les géomètres pour leurs opérations les plus courantes. cette limitation apparaît clairement lorsqu'il s'agit de déterminer la position d'un point sous couvert végétal, à l'aplomb d'un mur, sous un pont ou une trémie, et de façon plus générale, dans les environnements peu propices à la réception correcte des signaux satellitaires. La solution apportée à ce jour à ce type de situation consiste à déterminer, avec toute la précision requise, les coordonnées d'un ou plusieurs points dits "intermédiaires" situés dans un environnement dégagé, à proximité du point objectif masqué, puis à utiliser des techniques classiques de mesure d'angles et de distances (chaîne d'arpenteur, tachéomètre, distancemètre, etc) pour situer le point objectif par rapport à ces points intermédiaires. Cette association de techniques disparates rend les matériels lourds, les procédures complexes et les risques d'erreurs importants; en pratique, elle limite considérablement l'intérêt que présente l'utilisation des techniques satellitaires, puisque de toute façon les techniques classiques et les matériels associés restent indispensables dans de très nombreuses applications. Un but de l'invention est d'accroître considérablement le champ d'application des techniques radio-satellitaires. Elle propose à cet effet une méthode et un système de radiolocalisation basées sur l'exploitation des signaux satellitaires, qui permettent d'accéder à la position précise des points masqués sur lesquels la réception des signaux radio-satellitaires est impossible ou insuffisante. Selon le procédé proposé par l'invention: - on positionne l'antenne d'un mobile de radiolocalisation évoluant librement à proximité du point masqué sur une pluralité de points quelconques o la réception des signaux émis par les satellites est suffisante, - on détermine par radiolocalisation les coordonnées de ces points et - on détermine en fonction de ces coordonnées les coordonnées du point masqué, caractérisé en ce que, à cet effet, - on dispose un émetteur pseudolite sur le point masqué, - on mesure en chacun des points o l'antenne est successivement positionnée la phase du signal émis par l'émetteur pseudolite et reçu en ces points par ladite antenne, - on détermine en fonction de ces mesures et des coordonnées de ces points les coordonnées du point masqué. Ainsi selon l'invention, on exploite les mesures de phase du signal émis par l'émetteur pseudolite et les mesures issues des satellites permettant de calculer les positions successives du mobile localisateur, pour calculer les coordonnées du point masqué. On rappelle ici qu'un "pseudolite" ou "pseudo- satellite" est une unité d'émission au sol qui émet des signaux radioélectriques à codage pseudo-aléatoire qui sont analogues à ceux émis par les satellites du système de radiolocalisation utilisé (satellites de la constellation NAVSTAR pour les systèmes GPS). Les pseudolites sont classiquement connus de l'Homme du Métier. Ils sont principalement utilisés dans le domaine de l'aéronautique pour le guidage d'avions par rapport à des aires d'atterrissage ou pour la levée d'ambiguïtés d'entiers de phase. On pourra avantageusement à cet égard se référer aux brevets US-5.177. 489 et US-5.311.194 et aux publications suivantes: - Real time flight testing using integrity beacons for GPS, category III Precison Landing, CLARK COHEN and al. - Revue Navigation Vol. 41 N 2 - Summer 94 - Institute of Navigation, - Integrity Monitoring for precision approach using Kinematic GPS and a ground-based pseudolite - Boris PERVAN and al. - Revue Navigation - Vol. 41 n 2 - Summer 94. Notamment, dans une variante avantageuse du procédé conforme à l'invention: - pour chaque point o l'antenne est successivement positionnée, on détermine les paramètres d'une équation linéaire reliant les coordonnées du point masqué et le déphasage entre l'émetteur pseudolite et le mobile de radiolocalisation en un instant de référence, ces paramètres étant fonction de la mesure de la phase du signal pseudolite et des coordonnées de ces points, - on détermine les coordonnées du point masqué par résolution des équations linéaires ainsi obtenues. Selon l'invention, il est également proposé un dispositif qui comporte un émetteur pseudolite destiné à être placé sur un point à localiser o la réception des signaux émis par les satellites n'est pas suffisante, l'un des canaux des moyens de radiolocalisation, habituellement dédié à la réception des signaux radio-satellitaires, étant affecté à la réception du signal émis par l'émetteur peusdolite, ce dispositif comportant des moyens de traitement permettant de déterminer les coordonnées du point à localiser en fonction de la mesure de la phase du signal émis par l'émetteur pseudolite en différents points o l'antenne est successivement positionnée et des coordonnées de ces points déterminées par l'unité. D'autres caractéristiques et avantages ressorti- ront encore de la description qui suit. Cette description est purement illustrative et non limitative. Elle doit être lue en regard du dessin unique annexé sur lequel on a illustré schématiquement le principe de la méthode et du système proposé par l'invention. Sur cette figure, on a représenté un point masqué P dont on cherche à relever les coordonnées, qui est situé dans une zone couverte Zc dans laquelle les signaux S radio-satellitaires ne sont pas convenablement reçus. Le dispositif utilisé pour relever les coordonnées XYZ de ce point P comporte un émetteur pseudolite 1 et un mobile 2 de radiolocalisation G.P.S.. Le pseudolite 1 rayonne à une très faible puissance (quelques microwatts), à la fréquence GPS de 1,57542 GHz (L1), ou à une fréquence légèrement décalée, un signal Sp modulé en phase par un code pseudoaléatoire dont le format est celui des codes GPS. Le mobile 2 comporte une antenne de réception 3 montée sur une canne tenue par l'opérateur et une unité de traitement 4 portée par l'opérateur et reliée à ladite antenne 3. L'unité 4 traite les signaux radioélectriques reçus par l'antenne 3 en provenance des satellites GPS, du pseudolite 1, ainsi que d'une balise 5 de référence qui reçoit également les signaux en provenance des satellites GPS et du pseudolite 1 et qui retransmet des signaux de correction différentielle de phase (signaux SB sur la figure). Comme on l'aura compris, l'un des canaux du mobile de radiolocalisation 2, habituellement dédié à la réception des signaux radio-satellitaires, est affecté à la réception du signal Sp émis par le pseudolite 1, de sorte que l'unité 4 dispose d'une mesure de phase qui est fonction de la distance entre le récepteur mobile que constitue le mobile 2 et l'émetteur pseudolite 1. Les traitements que l'unité 4 réalise sur les signaux GPS sont des traitements cinématiques en temps réels, de sorte que la position de l'antenne est connue avec une précision centimétrique. Ces traitements sont par exemple des traitements du type de ceux décrits dans la demande de brevet français déposée par la demanderesse sous le numéro 94- 00.523, à laquelle on pourra avantageusement se reporter. L'opérateur déplace l'antenne 3 du récepteur mobile 2 sur une succession de points Pi o les signaux radio-satellitaires S sont convenablement reçus. En chacun de ces différents points Pi, le récepteur 2 détermine les coordonnées Xi, Yi, Zi du point Pi sur lequel il se trouve et réalise une mesure de phase Ri sur le signal émis par le pseudolite 1 et reçu par le récepteur GPS 2. La connaissance des coordonnées d'au moins quatre points Pi et des mesures de phase Ri correspondantes permet à l'unité 4 de déterminer les coordonnées X, Y et Z du point P. Les traitements mis en oeuvre à cet effet par l'unité 4 vont maintenant être explicités. Au point P0, de coordonnées X0, Y0, Z0, la mesure de phase est: R0 = Do + H0, o Do est une valeur de différence de phase proportionnelle à la distance entre le pseudolite 1 et le récepteur 2 et o H0 est le déphasage d'horloge entre le pseudolite 1 et le récepteur 2. En un temps T1 ultérieur, l'opérateur s'est déplacé en un point P1 de coordonnées X1, Y1, Zl et la mesure de phase que réalise le récepteur 2 est: R1 = D1 + H0 + AH1, o D1 est la nouvelle valeur de différence de phase correspondant à la distance entre le pseudolite 1 et le récepteur 2 et o AHlest l'évolution relative des horloges du récepteur 2 et du pseudolite 1 entre To et T1. En un temps T2, l'opérateur s'est déplacé en un point P2 de coordonnées X2, Y2, Z2. La mesure de phase qu'il réalise est: R2 = D2 + H0 + AH1 + AH2, o D2 est la différence de phase qui correspond à la distance entre le récepteur 2 et le pseudolite 1 en ce point P2, et o AH2est l'évolution relative des horloges du récepteur 2 et du pseudolite 1 entre T1 et T2. Il poursuit ses mesures sur plusieurs points Pi, la mesure de phase en un point Pi donné étant: i Ri = Di + Ho + X AHi j=i o Di est la différence de phase qui correspond à la distance entre le récepteur 2 et le pseudolite 1 en ce point Pi, et o AHiest l'évolution relative des horloges du récepteur 2 et du pseudolite 1 entre Ti_1 et Ti. Les évolutions AHi sont calculées par l'unité 4 par soustraction des évolutions d'horloge AHi(2_5) entre le récepteur 2 et la balise 5 et AHi(5_1) entre la balise 5 et le pseudolite 1. La valeur de l'évolution AHi(2_5> est calculée par l'unité 4 et est un sous-produit de la détermination des coordonnées du point Pi- La deuxième de ces évolutions ( AHi(5-1_)) est mesurée par la balise de correction 5 implantée en vue du pseudolite 1 et à laquelle est affectée de la même façon que pour le récepteur 2 un canal spécifique pour la réception des signaux émis par le pseudolite 1. La distance entre la balise 5 et le pseudolite 1 est en effet constante, de sorte que l'évolution de phase mesurée par la balise 5 sur le signal émis par le pseudolite 1 correspond à l'évolution relative des horloges de la balise 5 et du pseudolite 1. Cette information est transmise par la balise 5 au récepteur 2. L'unité 4 mémorise les valeurs successives des évolutions AHj et calcule en chaque point Pi (ième point) une valeur de phase Ri donnée par i R'i = Ri - AHj j=l c'est à dire: R'i = Di + H0. Cette valeur R'i vérifie l'équation linéaire suivante: (1) R'i2 - R02 + Ki = aiX + biY + ciZ + diH0 ou ai, bi, ci, di et Ki sont des paramètres que l'unité 4 calcule et qui sont fonction des coordonnées des points Pi et Po- Entre Ti et To on a en effet: Di2 _ Do2 + Xo2 - Xi 2 + y02 -Yi 2 + Zo2 - 2 = 2 (Xi - X0) X + 2 (Yi-Yo) Y+ 2 (Zi + Zo) Z, et on a également: Di 2 Do2 = R'i 2 - 2HoR'i + Ho2 - (Ro2 - 2HoRo + Ho2) = R'i 2 - Ro2 - 2 (R'i Ro) Ho On retrouve donc l'équation (1) avec ai, bi, ci, di respectivement égaux à 2 (Xi - Xo), 2 (Yi-Yo), 2 (zi - Z0), 2 (R'i - Ro), Ki étant égal à Xo2 - Xi 2 + y02 -Yi 2 + Zo2 - Zi 2. Comme on l'aura compris, l'équation (1) est une expression linéaire en X, Y, Z, Ho qu'il est possible de résoudre dès que l'on dispose des mesures en quatre points Pi distincts du point P0. En pratique, on utilise la méthode des moindres carrés et on réalise sur les valeurs X, Y, Z successivement déterminées pour un même point P un test de stabilité et l'on arrête les mesures lorsque la stabilité des valeurs de ces coordonnées X, Y, Z est satisfaisante. Par exemple, la solution obtenue est validée lorsque le bruit des trois inconnues X, Y, Z est inférieur à un seuil prédéterminé fonction de la précision requise pour le travail à effectuer (en pratique 1 à quelques cm). Le temps de convergence de la solution est fonction de la géomètrie d'observation du pseudolite et varie de 1 à quelques minutes. Comme on l'aura compris, le procédé et le dispositif qui viennent d'être décrits permettent de déterminer la position d'un point sur lequel la réception des signaux satellitaires est impossible, en ne mettant en oeuvre que des techniques de radiolocalisation. Toutes les mesures sont faites en un seul point (à savoir le centre de phase de l'antenne 3), ce qui élimine les problèmes de raccordement entre mesures issues de différents systèmes ou capteurs. On notera que les traitements qui viennent d'être décrits s'effectuent à partir d'observations sur les variations de phase des signaux reçus, et non sur la phase de ces signaux, ce qui élimine les problèmes liés à la détermination des entiers de phase. Sur le plan opérationnel, le procédé et le système conforme à l'invention présentent de nombreux avantages. La détermination du point objectif peut se faire au cours d'une libre évolution du récepteur mobile, sans aucune contrainte de positionnement initial ou en des points relais. La détermination des coordonnées du point objectif peut se faire simultanément à un levé ou une implantation classique en d'autres points. Le procédé peut être activé ou interrompu sans perturber le fonctionnement normal de la cinématique temps réel en ces points. La liaison entre l'émetteur pseudolite et le récepteur mobile peut être momentanément interrompue, par exemple par masquage, sans qu'il n'en résulte aucune conséquence sur la qualité du résultat. Le système accumule et traite les mesures de temps nécessaire à l'obtention d'un résultat selon des critères de précision qui ont été préalablement déterminés par l'opérateur. Lorsque la qualité requise est obtenue, le calcul se termine automatiquement et l'opérateur en est averti par l'unité 4. On notera également que le système proposé par l'invention est d'une grande simplicité de manipulation: hormi la mise en place de l'émetteur pseudolite, aucune autre manipulation n'est nécessaire. REVENDICATIONS 1. Procédé de localisation radio-satellitaire de points à la surface de la terre, dans lequel, pour localiser un point masqué o la réception des signaux émis par les satellites n'est pas suffisante, - on positionne l'antenne (3) d'un mobile de radiolocalisation évoluant librement à proximité du point masqué sur une pluralité de points quelconques (Pi) o la réception des signaux émis par les satellites est suffisante, - on détermine par radiolocalisation les coordonnées de ces points (Pi) et - on détermine en fonction de ces coordonnées les coordonnées du point masqué (P), caractérisé en ce que, à cet effet, - on dispose un émetteur pseudolite (1) sur le point masqué (P), - on mesure en chacun des points (Pi) o l'antenne (3) est successivement positionnée la phase du signal émis par l'émetteur pseudolite (1) et reçu en ces points par ladite antenne (3), - on détermine en fonction de ces mesures et des coordonnées de ces points (Pi) les coordonnées du point masqué (P). 2. Procédé selon la revendication 1, caractérisé en ce que - pour chaque point (Pi) o l'antenne (3) est successivement positionnée, on détermine les paramètres d'une équation linéaire reliant les coordonnées du point masqué (P) et le déphasage entre l'émetteur pseudolite (1) et le mobile de radiolocalisation en un instant (To) de référence, ces paramètres étant fonction de la mesure de la phase du signal pseudolite (Sp) et des coordonnées de ces points (Pi), - on détermine les coordonnées du point masqué (P) par résolution des équations linéaires ainsi obtenues. 3. Procédé selon la revendication 2, caractérisé en ce que pour le ième point (Pi) o l'antenne (3) est successivement positionnée - on détermine l'évolution AHi relative des horloges de l'émetteur pseudolite (1) et d'une unité (4) de traitement du mobile de radiolocalisation, entre ce point (Pi) et le point précédent (Pi-1), - on calcule la valeur de phase i R'i = Ri - AHj j=l o Ri est la mesure de phase du signal pseudolite (Ps) en ce point (Pi), les valeurs AHj des évolutions relatives aux points (Pj) précédents ayant à cet effet été mémorisées, - on calcule les paramètres Ki, ai, bi, ci, di de l'équation linéaire suivante: (1) R'i2 - R02 + Ki = aiX + biY + ciZ + diH0, o X, Y, Z sont les coordonnées du point masqué (P), H0 est le déphasage entre l'émetteur pseudolite (1) et le mobile de radiolocalisation à l'instant de référence (To), ai, bi, ci, di étant respectivement égaux à 2 (Xi - X0), 2 (Yi-Yo), 2 (Zi - Z0), 2 (R'i - R0o), o Xi, Yi, Zi sont les coordonnées du point (Pi), X0, Y0, Z0 les coordonnées du point P0 o se trouve l'antenne (3) à l'instant de référence, R0 étant la mesure de phase du signal pseudolite à l'instant de référence, Ki étant égal à X02 - Xi 2 + Y02 _yi 2 + Z02 - Zi 2 4. Procédé selon la revendication 3, caractérisé en ce que l'on détermine l'évolution relative des horloges de l'émetteur pseudolite (1) et de l'unité (4) de traitement du mobile de radiolocalisation par soustraction: 10. - de l'évolution d'horloge entre ladite unité (4) et une balise de référence (5), - et de l'évolution d'horloge entre cette balise de référence (5) et l'émetteur pseudolite (1). 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'on réalise sur les valeurs (X, Y, Z) successivement déterminées pour les coordonnées du point à localiser (P) un test de stabilité, la radiolocalisation étant terminée lorsque la stabilité de ces valeurs est suffisante. 6. Dispositif de localisation radio-satellitaire de points à la surface de la terre comportant un mobile de radiolocalisation (2) qui comprend une antenne (3) pour la réception des signaux radio-électriques (S) émis en direction de la terre par une constellation de satellites (GPS) et une unité (4) pour le traitement de ces signaux, ce dispositif étant caractérisé en ce qu'il comporte en outre un émetteur pseudolite (1) destiné à être placé sur un point masqué (P) o la réception des signaux émis par les satellites n'est pas suffisante, l'un des canaux du mobile de radiolocalisation (2), habituellement dédié à la réception des signaux radio-satellitaires (S), étant affecté à la réception du signal (Sp) émis par l'émetteur peusdolite (1), ce dispositif comportant des moyens de traitement (4) permettant de déterminer les coordonnées du point masqué (P) en fonction d'une part de la mesure de la phase du signal émis par l'émetteur pseudolite (1) en différents points (Pi) o l'antenne (3) est successivement positionnée et d'autre part des coordonnées de ces points (Pi) déterminées par l'unité (4). 7. Dispositif selon la revendication 6, caractérisé en ce que les moyens de traitement (4) comportent: - des moyens de calcul pour déterminer, pour chaque point (Pi) o l'antenne (3) est successivement positionnée, les paramètres d'une équation linéaire reliant les coordonnées du point masqué (P) et le déphasage entre l'émetteur pseudolite (1) et le mobile de radio-localisation (2) en un instant (To) de référence, ces paramètres étant fonction de la mesure de la phase du signal pseudolite (Ps) et des coordonnées de points (Pi), et - des moyens de calcul pour résoudre les équations linéaires ainsi obtenues. 8. Dispositif selon la revendication 7, caractérisé en ce qu'il comporte: - des moyens pour déterminer l'évolution relative AHi des horloges de l'émetteur pseudolite (1) et de l'unité (4) de traitement des moyens de radiolocalistion entre le ième point (Pi) o l'antenne (3) est successivement positionnée et le point (Pi-1) précédent, - des moyens de calcul de la valeur de phase i R'i = Ri - XAHj j=l o Ri est la mesure de phase du signal pseudolite en ce point (Pi), les valeurs AHj des évolutions relatives aux points (Pi Pj) précédents ayant à cet effet été mémorisées, - des moyens de calcul des paramètres Ki, ai, bi, ci, di de l'équation linéaire (1) R'i2 - R02 + Ki = aiX + biY + ciZ + diH0, o X, Y, Z sont les coordonnées du point masqué (P), H0 est le déphasage entre l'émetteur pseudolite (1) et le mobile de radiolocalisation en un temps To de référence, ai, bi', ci, di étant respectivement égux à 2 (Xi - X0), 2 (Yi-Yo), 2 (Zi - Z0), 2 (R'i - R0), o Xi, Yi, Zi sont les coordonnées du point Pi, X0, Y0, Z0 les coordonnées du point P0 o se trouve l'antenne (3) à l'instant de référence, Ro étant la mesure de phase du signal pseudolite à l'instant de référence, Ki étant égal à X02 - Xi 2 + y02 -Yi 2 + Z02 - Zi 2. 9. Dispositif selon la revendication 8, caractérisé en ce que les moyens pour déterminer l'évolution relative des horloges de l'émetteur pseudolite (1) et de l'unité (4) de traitement comportent des moyens de soustraction de l'évolution d'horloge entre l'unité (4) de traitement et une balise de référence (5) et de l'évolution d'horloge entre cette balise de référence (5) et l'émetteur pseudolite (1). 10. Dispositif selon l'une des revendications 6 à 9, caractérisé en ce qu'il comporte des moyens pour réaliser un test de stabilité sur les valeurs (X, Y, Z) successivement déterminées pour les coordonnées du point à localiser (P).