Brevet FR3109827A1 — publié le 20211105

L’invention concerne un procédé de géolocalisation d’un équipement utilisateur situé à une altitude donnée, recevant une pluralité de faisceaux radiofréquences émis selon au moins une première et une deuxième direction de propagation respectivement par au moins une antenne d’une première station de base et par au moins une antenne d’une deuxième station de base selon au moins une troisième et une quatrième direction de propagation. Ledit procédé met en œuvre : - une obtention (E300) d’une puissance mesurée par ledit équipement terminal pour des faisceaux radiofréquences émis par la première station de base au moins selon une première et une deuxième direction de propagation et des faisceaux radiofréquences émis par la deuxième station de base au moins selon une troisième et une quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; - un calcul (E310) d’au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative pour la première station de base et deuxième information de puissance relative pour la deuxième station de base, la première puissance étant associée pour première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de rayonnement ; et - une détermination (E320) d’une position dudit équipement utilisateur à l’altitude donnée à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative et, d’autre part d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction de la direction d’observation de l’équipement utilisateur par ladite antenne. FIGURE 4 Procédé de géolocalisation d’un équipement utilisateur, dispositif, équipement utilisateur, station de base, système et programme d’ordinateur correspondants. Domaine technique de l'invention La présente invention appartient au domaine général des télécommunications. Elle concerne plus particulièrement un réseau de communication sans fil et un procédé de localisation d’un utilisateur d’un équipement terminal dans un tel réseau de télécommunications. Art antérieur De nombreux services ou applications proposé(e)s aujourd’hui par les opérateurs de réseaux de télécommunications s’appuient sur une localisation géographique des terminaux des utilisateurs. Ces services sont connus sous le nom de services LSB (pour « Location Based Services », en anglais) et suscitent un grand intérêt aussi bien pour des applications commerciales (jeu, publicité, etc) que pour des applications liées à la sécurité (interventions d’urgence, médical, services d’alerte, etc). On connaît en effet différentes techniques pour déterminer la position géographique d’un utilisateur, c’est-à-dire pour le géolocaliser. Une technique très répandue s’appuie sur des signaux diffusés par des satellites en orbite autour de la Terre. Selon cette technique, le système le plus connu est le système GPS (pour « Global Positioning System », en anglais). Chaque satellite envoie des signaux indiquant sa position dans l'espace ainsi que l'heure et la date d'émission desdits signaux. Un récepteur GPS embarqué dans un terminal de l’utilisateur réceptionne ces signaux, calcule ensuite la durée de parcours de chaque signal entre le satellite émetteur et le récepteur, et en déduit enfin, par trilatération, sa position, en latitude, longitude et altitude, sur le globe. Pour que le système fonctionne, il faut compiler les données reçues de quatre satellites : trois pour la position et un pour la synchronisation. Avec un tel système, le terminal d’un utilisateur, qu'il soit sur terre, sur mer ou dans les airs, peut connaître sa position à toute heure et en tout lieu sur la surface ou au voisinage de la surface de la Terre, dès lors qu'il est équipé d'un récepteur et du logiciel nécessaire au traitement des informations reçues. La précision de l’estimation fournie par le GPS par rapport à sa localisation réelle dépend du nombre de satellites visibles par le système GPS, qui peut fortement varier en fonction des conditions climatiques, mais pour un terminal mobile de type téléphone intelligent (pour « smartphone », en anglais), elle peut atteindre 5 à 10 m. En revanche, les systèmes GPS sont gourmands en termes de ressources ce qui peut s’avérer problématique lorsqu’ils équipent des terminaux sur batterie, comme les téléphones mobiles. C’est précisément le cas des équipements terminaux d’utilisateurs situés sur des bateaux naviguant en mer. Il existe donc un besoin d’une technique de localisation d’un équipement utilisateur qui reste suffisamment précise pour constituer une bonne alternative au GPS, tout en étant plus économe en ressources. Présentation de l'invention L'invention répond à ce besoin en proposant un procédé de géolocalisation d’un équipement utilisateur situé à une altitude donnée, recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne d’une première station de base configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne d’une deuxième station de base configurée pour émettre selon au moins une troisième et une quatrième direction de propagation. Le procédé met en œuvre : - une obtention d’une puissance mesurée par ledit équipement terminal pour des faisceaux radiofréquences émis par la première station de base selon au moins la première et la deuxième direction de propagation et des faisceaux radiofréquences émis par la deuxième station de base selon au moins la troisième et la quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; - un calcul d’au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative pour la première station de base et deuxième information de puissance relative pour la deuxième station de base, la première puissance étant associée pour première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de propagation; et - une détermination d’une position dudit équipement utilisateur à l’altitude donnée à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction de la direction d’observation de l’équipement utilisateur par ladite antenne. Ainsi l’invention repose sur une approche tout à fait nouvelle et inventive pour géolocaliser un équipement utilisateur situé à une altitude donnée et à portée d’une première et d’une deuxième stations de base (e.g. connecté à l’une des stations de base) d’un réseau de radiocommunications (e.g. un réseau de radiocommunications conformes aux normes 3GPP). Plus particulièrement, la position de l’équipement utilisateur est déterminée à partir de mesures de puissance de différents signaux émis depuis deux stations de base capables d’émettre chacune selon deux directions de propagation distinctes (e.g. des signaux émis via deux antennes sectorielles de chaque station de base, ou des signaux correspondant à différents faisceaux d’une antenne de type matrice d’éléments rayonnants équipant chaque station de base). Par ailleurs, cette détermination s’appuie sur le calcul d’informations de puissances relatives à partir des puissances mesurées. De telles informations sont représentatives d’un rapport (lorsque les puissances sont exprimées en unité naturelle) ou d’une différence (lorsque les puissances sont exprimées en unité logarithmique) entre les deux puissances de deux faisceaux radiofréquences donnés. De la sorte, les effets de l’atténuation de propagation (ou « Path Loss » en anglais), des évanouissements rapides du canal de propagation (ou « Fast Fading » en anglais), ainsi que des effets d’occultation ou de masque (ou « Shadowing » en anglais) sont annulés. En effet, ce type de phénomènes impacte de la même façon les deux faisceaux radiofréquences considérés, ces faisceaux étant émis par des éléments rayonnants localisés sensiblement en un même point géographique de la même station de base. Le calcul d’une puissance relative permet ainsi d’annuler les effets de ces phénomènes sur les signaux reçus par l’équipement utilisateur. On note que selon une technologie d’antennes de type SISO (pour « Single Input Single Output » en anglais), dans le cas où les antennes ne sont pas exactement localisées au même point, une moyenne des puissances mesurées sur une période donnée permet d’éliminer les effets de fading. Selon une configuration d’antennes de type MIMO (pour « Multiple-Input Multiple-Output » en anglais), un tel moyennage n’est pas nécessaire. Enfin, l’utilisation de modèles de diagrammes de rayonnement permet de remonter de manière simple et robuste depuis la ou les informations de puissances relatives jusqu’à la position de l’équipement utilisateur, sans qu’il soit nécessaire de connaître les caractéristiques des antennes de la deuxième station de base. L’invention s’applique à la géolocalisation de tout type d’équipement utilisateur pourvu qu’il soit simultanément à portée de faisceaux émis par une première station de base respectivement une deuxième station de base selon deux directions de propagation distinctes et que son altitude soit connue. Elle est donc particulièrement bien adaptée à la localisation de bateaux en mer qui peuvent, du fait que les antennes d’une station de base terrestre sont plutôt tournées vers les terres, ne recevoir qu’une partie des signaux émis par cette station de base, comme par exemple les faisceaux émis par deux des trois antennes de cette station de base. Néanmoins, elle s’applique aussi à la géolocalisation d’équipements utilisateurs situés dans tout autre type d’environnement géographique dans lequel l’installation de stations de base à deux antennes ou à une seule antenne comprenant une matrice d’éléments rayonnants, est préconisée ou rencontrée. Selon un aspect de l’invention, ladite détermination met en œuvre, pour au moins une information de puissance relative donnée parmi la ou les informations de puissances relatives, la résolution d’une équation dont les membres sont fonction, d’une part, de ladite information de puissance relative calculée et, d’autre part, d’une valeur attendue de ladite information de puissance relative donnée, fonction, pour chaque faisceau radiofréquences associé à ladite information de puissance relative donnée, du modèle de diagramme de rayonnement caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative donnée, en fonction de la direction d’observation de l’équipement utilisateur par ladite antenne. Ainsi, la position de l’équipement utilisateur est déterminée de manière simple et robuste en mettant en adéquation la ou les informations de puissances relatives calculées à partir des puissances mesurées avec la ou les valeurs attendues telles que prédites par les modèles de diagramme de rayonnement des antennes utilisées. Selon un autre aspect de l’invention, la première station de base comprend une première antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la première direction de propagation et une deuxième antenne sectorielle configurée pour émettre des faisceaux radiofréquences selon la deuxième direction de propagation, et la deuxième station de base comprend une première antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la troisième direction de propagation et une deuxième antenne sectorielle configurée pour émettre des faisceaux radiofréquences selon la quatrième direction de propagation, ladite obtention comprend l’obtention d’une première et d’une deuxième puissances mesurées correspondant respectivement à un premier faisceau radiofréquences émis par la première, respectivement la deuxième antenne sectorielle de ladite première, respectivement deuxième station de base et à un deuxième faisceau radiofréquences émis par la première, respectivement la deuxième antenne sectorielle de ladite première , respectivement deuxième station de base. Le calcul délivre une première, respectivement une deuxième information de puissance relative associée auxdits premier et deuxième faisceaux radiofréquences pour chaque station de base, et ladite détermination comprend la détermination d’un angle représentatif d’une longitude dudit équipement utilisateur dans un repère de la première, respectivement de la deuxième, station de base à partir, d’une part, de ladite première, respectivement de la deuxième information de puissance relative et d’autre part, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction desdits angles de longitude et de latitude de l’équipement utilisateur, desdits premier et deuxième faisceaux radiofréquences. Un avantage de ce mode de réalisation est d’obtenir assez simplement un angle de longitude de l’équipement utilisateur dans le repère de chaque station de base. Selon encore un autre aspect de l’invention, ladite résolution comprend, pour chaque station de base, la mise en œuvre de l’équation : θ1i = - (M1i dB /24 ) θ 2 3dB / θb i + θb i /2, où : i est un entier qui vaut 1 lorsqu’il désigne la première station de base et 2 lorsqu’il désigne la deuxième station de base ; θ1i représente ledit angle représentatif d’une longitude dudit équipement utilisateur dans lerepère de la i-ème station de base, M1i dB représente ladite i-ème information de puissance relative exprimée en décibels, θ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement de chacune desdites première et deuxième antennes dans un plan de définition dudit angle θ1i, et θb i représente l’angle de dépointage entre lesdites première et deuxième antennes de la i-ème station de base dans ledit plan de définition dudit angle θ1i. Ainsi, l’angle représentatif de la longitude de l’équipement utilisateur est déterminé de manière simple et robuste dans le repère de chacune des deux stations de base, lorsqu’elles implémentent une technologie du type SISO sur différents secteurs de couverture. Selon encore un autre aspect de l’invention, lesdits premier et deuxième faisceaux radiofréquences sont rayonnés par lesdites première et deuxième antennes de la première station de base selon un premier angle d’inclinaison, ladite obtention comprend l’obtention d’une troisième puissance mesurée correspondant à un troisième faisceau radiofréquences émis par ladite première antenne de la première station de base selon un troisième angle d’inclinaison, distinct du premier angle d’inclinaison, ledit calcul délivre une troisième information de puissance relative, représentative d’un rapport ou d’une différence entre la troisième puissance mesurée et la première puissance du premier jeu de puissances mesurées ; et ladite détermination comprend une première détermination d’un angle représentatif d’une latitude dudit équipement utilisateur dans lerepère de la première station de base à partir, d’une part de ladite troisième information de puissance relative et, d’autre part d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction desdits angles de longitude et de latitude de l’équipement utilisateur, desdits premier, deuxième, et troisième faisceaux radiofréquences de la première station de base. Ainsi, l’angle noté ci-après φ1 représentatif d’une latitude de l’équipement utilisateur est déterminé de manière simple et précise à partir des faisceaux émis par les deux antennes de la première station de base lorsqu’elles émettent des faisceaux radiofréquences selon deux angles d’inclinaison différents. On note que l’angle d’inclinaison des antennes de la deuxième station de base peut être identique à celui de la première. Selon encore un autre aspect de l’invention, ladite détermination comprend la mise en œuvre de l’équation : φ11 = M12 dB /24 φ 2 3dB / (φt1– φt3) + (φt1+ φt3)/2, où : φ11 représente ledit angle représentatif d’une latitude dudit équipement utilisateur dans le repère de la première station de base, M12 dB représente ladite deuxième information de puissance relative exprimée en décibels, φ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement de chacune desdites première et deuxième antennes dans un plan de définition dudit angle φ11, φt1 représente ledit premier angle d’inclinaison, et φt3 représente ledit troisième angle d’inclinaison. Selon encore un autre aspect de l’invention, ladite détermination d’une position comprend une deuxième détermination d’une valeur de l’angle de latitude de l’équipement utilisateur dans le repère de la première station de base, en fonction d’une hauteur de la première antenne de la première station de base, et des angles de longitude de l’équipement utilisateur dans les repères respectifs de la première et de la deuxième station de base et une sélection de la valeur des angles de longitude de l’équipement utilisateur minimisant une erreur entre la première et la deuxième détermination de l’angle de latitude. Ainsi, on utilise la connaissance de la position relative des deux stations de base et des relations trigonométriques entre les angles pour lever une ambiguïté sur la valeur d’angle de longitude de l’équipement utilisateur et déterminer l’angle de latitude de l’équipement utilisateur à l’imprécision de mesure de puissance près . Selon encore un autre aspect de l’invention, ladite détermination comprend le calcul des coordonnées (x1,y1,z1) de l’équipement utilisateur dans le repère de la première station de base à partir des angles de longitude déterminés. Ainsi, les coordonnées (x1, y1, z1) de l’équipement utilisateur par rapport à la première station de base sont déterminées de manière simple et robuste, sans recourir à des caractéristiques de la deuxième station de base, comme par exemple la hauteur de ses antennes. Selon un autre aspect de l’invention, la première station de base comprend une antenne comprenant une matrice d’éléments rayonnants, configurée pour émettre des faisceaux radiofréquences selon au moins la première direction de propagation et la deuxième direction de propagation et la deuxième station de base comprend une antenne comprenant une matrice d’éléments rayonnants configurée pour émettre des faisceaux radiofréquences selon au moins la troisième direction de propagation et la quatrième direction de propagation et chaque faisceau radiofréquences de ladite pluralité de faisceaux radiofréquences est rayonné par la matrice d’éléments rayonnants de ladite première, respectivement deuxième station de base. Ainsi, la méthode selon l’invention s’applique aussi dans le cas d’une station de base implémentant une technologie du type MIMO. Selon un autre aspect de l’invention, ladite résolution d’équation comprend la mise en œuvre d’une méthode de résolution numérique de ladite au moins une équation. Selon encore un autre aspect, ladite résolution met en œuvre pour ladite au moins une information de puissance relative donnée : une obtention de ladite valeur attendue de ladite information de puissance relative donnée pour un jeu de directions d’observation de l’équipement utilisateur délivrant un jeu de valeurs attendues correspondant chacune à une direction d’observation de l’équipement utilisateur parmi une pluralité de directions d’observations; une comparaison entre, d’une part, ladite information de puissance relative calculée et, d’autre part, chaque valeur attendue dudit jeu de valeurs attendues délivrant une direction d’observation pour laquelle la valeur attendue de ladite puissance relative donnée est la plus proche de l’information de puissance relative calculée ; et ladite position dudit équipement utilisateur étant fonction de la direction d’observation délivrée. Ainsi, la résolution de la ou des équations mises en jeu se fait de manière simple et robuste par exploration de l’espace des solutions possibles. Dans des modes de réalisation, la direction d’observation possible correspond à la valeur attendue ou théorique d’information de puissance relative parmi les valeurs attendues du jeu de valeurs attendues qui est la plus proche de l’information de puissance relative calculée à partir des puissances mesurées, à l’imprécision de mesure de puissance près. Selon un autre aspect de l’invention, pour la première, respectivement la deuxième, station de base, ledit modèle de diagramme de rayonnement d’antenne caractérisant la puissance d’un i-ème faisceau radiofréquences en fonction d’une direction d’observation de l’équipement utilisateur (θ, φ) s’exprime en fonction de : avec : où : θ représente un angle représentatif d’une longitude de l’équipement utilisateur dans un repère de la première, respectivement deuxième, station de base, φ représente un angle représentatif d’une latitude de l’équipement utilisateur dans le repère de la première, respectivement deuxième station de base, (θ i,escan , φ i,etilt ) représente un couple d’angles caractéristique de la direction de propagation principale voulue pour ledit i-ème faisceau radiofréquences dans un repère de l’antenne, N H représente le nombre d’éléments rayonnants de ladite matrice d’éléments rayonnants dans une direction horizontale, N V représente le nombre d’éléments rayonnants de ladite matrice d’éléments rayonnants dans une direction verticale, d V représente l’espacement vertical entre deux éléments rayonnants, d H représente l’espacement horizontal entre deux éléments rayonnants, λ représente la longueur d’onde du faisceau radiofréquences, et A E (θ,φ) représente le diagramme de rayonnement de chaque élément rayonnant de ladite matrice d’éléments rayonnants ; et ladite détermination comprend une première détermination d’au moins deux couples de valeurs de l’angle de longitude et de l’angle de latitude (θ1, φ1) de l’équipement utilisateur dans le repère de la première station de base à partir dudit modèle de rayonnement, une deuxième détermination d’au moins deux couples de valeurs de l’angle de longitude et de l’angle de latitude de l’équipement utilisateur dans le repère de la deuxième station de base à partir dudit modèle de rayonnement, une troisième détermination d’une deuxième valeur de l’angle de latitude de l’équipement utilisateur (UE) dans le repère de la première station de base, en fonction d’une hauteur de la première station de base et des angles de longitude de l’équipement utilisateur dans les repères respectifs de la première et de la deuxième station de base et une sélection des valeurs des angles de longitude de l’équipement utilisateur minimisant une erreur entre la valeur de l’angle de longitude associé à l’angle de latitude issue de la première détermination et la deuxième valeur de l’angle de latitude issue de la troisième détermination. Ainsi, on détermine des couples de valeurs d’angles de longitude et de latitude possibles de l’équipement utilisateur en appliquant le modèle de diagramme de rayonnement de chaque élément rayonnant de chaque station de base et on utilise ensuite la connaissance de la position relative des deux stations de base et des relations trigonométriques entre les angles pour lever l’ambiguïté sur la valeur d’angle de longitude de l’équipement utilisateur et l’angle de latitude φ1 de l’équipement utilisateur, à l’imprécision de mesure de puissance près. De cette manière, on peut obtenir assez simplement les angles de longitude et de latitude de l’équipement utilisateur à partir d’un système d’équations complexe. L’invention concerne également un produit programme d'ordinateur comprenant des instructions de code de programme pour la mise en œuvre d’un procédé de géolocalisation d’un équipement utilisateur selon l’invention, tel que décrit précédemment, lorsqu’il est exécuté par un processeur. L’invention vise également un support d’enregistrement lisible par un ordinateur sur lequel sont enregistrés les programmes d’ordinateur tels que décrits ci-dessus. Un tel support d'enregistrement peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une clé USB ou un disque dur. D'autre part, un tel support d'enregistrement peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens, de sorte que le programme d’ordinateur qu’il contient est exécutable à distance. Le programme selon l'invention peut être en particulier téléchargés sur un réseau par exemple le réseau Internet. Alternativement, le support d'enregistrement peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé de géolocalisation précité. L’invention concerne aussi un dispositif de géolocalisation d’un équipement utilisateur recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne d’une première station de base configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne d’une deuxième station de base configurée pour émettre selon au moins une troisième et une quatrième direction de propagation. Ledit dispositif comprend une machine de calcul reprogrammable ou une machine de calcul dédiée, configurée pour : - obtenir une puissance mesurée par ledit équipement utilisateur pour des faisceaux radiofréquences émis par la première station de base selon au moins la première et la deuxième direction de propagation et des faisceaux radiofréquences émis par la deuxième station de base selon au moins la troisième et quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; - calculer au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative pour la première station de base et deuxième information de puissance relative pour la deuxième station de base, la première puissance étant associée pour la première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de propagation; et - déterminer une position dudit équipement utilisateur à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction d’une direction d’observation de l’équipement utilisateur par ladite antenne. Avantageusement, ledit dispositif est configuré pour mettre en œuvre le procédé de géolocalisation d’un équipement utilisateur précité, selon ses différents modes de réalisation. Avantageusement, ledit dispositif peut être intégré dans un équipement du réseau de communication. Il s’agit par exemple d’une station de base dotée d’aux moins deux antennes selon une technologie de type SISO émettant selon des directions de propagation distinctes ou d’une antenne comprenant une matrice d’éléments rayonnants selon une technologie de type MIMO. En variante, le dispositif de géolocalisation selon l’invention est intégré dans un équipement utilisateur apte à se connecter audit réseau de communication. Corrélativement, l’invention concerne enfin un système de géolocalisation d’un équipement utilisateur recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne d’une première station de base configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne d’une deuxième station de base configurée pour émettre selon au moins une troisième et une quatrième direction de propagation. Un tel système comprend au moins la première et la deuxième stations de base et un dispositif de géolocalisation précités Le système de géolocalisation, l’équipement utilisateur, la station de base, le dispositif de géolocalisation et le programme d'ordinateur correspondants précités présentent au moins les mêmes avantages que ceux conférés par le procédé précité selon les différents modes de réalisation de la présente invention. Brève description des figures D'autres buts, caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante, donnée à titre de simple exemple illustratif, et non limitatif, en relation avec les figures, parmi lesquelles : représente un équipement utilisateur connecté à une station de base d’un réseau de radiocommunications; représente les secteurs couverts par deux antennes sectorielles de la station de base de la fig. 1a ; représente les secteurs couverts par trois antennes sectorielles de la station de base de la fig. 1a ; représente une antenne comprenant une matrice d’éléments rayonnants pouvant équiper la station de base des figures 1a et 1b selon un mode de réalisation de l’invention ; représente une station de base équipée de la matrice d’éléments rayonnants de la figure 2a ; illustre un exemple d’architecture d’un système de géolocalisation d’un équipement utilisateur selon un mode de réalisation de l’invention ; représente les étapes du procédé de géolocalisation selon un mode de réalisation de l’invention ; illustre de façon géométrique le procédé de géolocalisation d’un équipement utilisateur selon un mode de réalisation de l’invention dans le système de la figure 3 ; détaille la détermination d’une position d’un équipement utilisateur selon un mode de réalisation de l’invention mettant en œuvre une technologie de type SISO ; illustre de façon géométrique la détermination d’une position d’un équipement utilisateur selon ce mode de réalisation de l’invention selon ce mode de réalisation de l’invention ; détaille la détermination d’une position d’un équipement utilisateur selon un autre mode de réalisation de l’invention mettant en œuvre une technologie de type MIMO ; et présente un exemple de structure matérielle d’un dispositif permettant la mise en œuvre des étapes du procédé de géolocalisation selon un mode de réalisation de l’invention. Description détaillée de l’invention Principe général de l’invention Le principe général de l'invention repose sur l’obtention de puissances mesurées par un équipement utilisateur situé à une altitude donnée, pour des faisceaux radiofréquences qu’il reçoit d’une première station de base et d’une deuxième station de base, chaque station de base comprenant au moins une antenne configurée pour émettre des faisceaux radiofréquences selon au moins deux directions de propagation, sur le calcul d’au moins une information de puissance relative par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances mesurées et sur une détermination d’une position dudit équipement utilisateur à l’altitude donnée à partir, d’une part, desdites informations de puissance relative et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction d’une direction d’observation de l’équipement utilisateur par l’antenne émettrice, dudit faisceau radiofréquences associé à ladite information de puissance relative. L’invention fonctionne aussi bien avec une technologie SISO que MIMO et trouve de nombreuses applications, en particulier dans la géolocalisation des bateaux en mer. Elle permet plus généralement de géolocaliser tout équipement utilisateur situé à la portée de deux stations de base, à partir de faisceaux radiofréquences émis selon des directions de propagation distinctes par chacune de ces stations de base. Dans la suite de la description, on s’attache à décrire en détails un exemple de mise en œuvre de l’invention. On présente désormais, en relation avec les figures 1a -1c un équipement utilisateur UE recevant les faisceaux radiofréquences émis par une première station de base BS1 d’un réseau de radiocommunications selon un mode de réalisation de l’invention basé sur une technologie de type SISO. Dans ce mode de réalisation, l’équipement utilisateur UE est connecté à la station de base BS1. Dans d’autres modes de réalisation, l’équipement utilisateur UE n’est pas connecté à la première station de base BS1, par exemple lorsqu’il met en œuvre lui-même le procédé de géolocalisation selon l’invention. Dans ce cas, l’équipement utilisateur UE n’a pas besoin de remonter à la première station de base BS1 les mesures de puissances qu’il effectue. De retour aux figures 1a et 1b , le réseau de radiocommunications est un réseau cellulaire, tel que par exemple un réseau 2G, 3G, 4G ou 5G défini par le standard 3GPP ou un autre standard. Dans la présente demande, une station de base est définie comme étant dédiée à la gestion d’un site géographique donné (par exemple un site géographique correspond à une cellule du réseau). Dans le cas des fig. 1a et fig. 1b, la première station de base BS1 gère le site géographique correspondant de manière multisectorielle (ou multi-secteurs). Plus particulièrement, la première station de base BS1 de la figure 1b couvre le site via deux secteurs distincts, chaque secteur étant couvert par une antenne dite sectorielle correspondante A11, A21. Dans la présente demande, on entend par antenne sectorielle une antenne émettant principalement dans une direction de propagation donnée. Par exemple, une cellule du réseau de radiocommunications comprend deux secteurs. Les 2 secteurs sont supposés ici de dimensions identiques. Par souci de simplification, on suppose que chaque secteur est défini par un angle, θb 1 , de dépointage dans le plan horizontal qui vaut dans cet exemple θb1 = 120°. C’est cette valeur de θb1 que l’on utilisera dans l’exemple de réalisation de l’invention détaillé ci-après. On notera que d’autres valeurs d’angles de dépointage peuvent être choisies. Chaque secteur est couvert au moyen d’une unique antenne A11, A21 apte à émettre selon un unique faisceau dans une direction de propagation donnée (ou tout du moins un unique faisceau principal concentrant la majeure partie de la puissance rayonnée par l’antenne) sur une bande de fréquences donnée. Autrement dit, les directions de propagation (ou de rayonnement) des faisceaux émis par les antennes A11, A21 couvrant deux secteurs adjacents de ce site présentent entre elles un angle égal à θb1. Chaque antenne A11, A21 est caractérisée, de façon connue en soi, par un diagramme de rayonnement. Bien sûr, comme déjà évoqué, l’invention s’applique aussi dans le cas où la station de base est dotée de trois antennes A11, A21, A31, comme illustré par la figure 1c. Dans ce cas, il suffit que l’équipement utilisateur UE reçoive les faisceaux radiofréquences émis par deux d’entre elles seulement pour que l’invention puisse s’appliquer. On considère par exemple des antennes telles que décrites dans le document Report ITU-R M.2135-1 de l’ITU-R, intitulé « Guidelines for évaluation of radio technologies for IMT Advanced » de décembre 2009. Le diagramme de rayonnement de chaque antenne A11, A21 présente un angle d’ouverture à trois décibels dans le plan horizontal noté θ 3dB . Les antennes A11 et A21 sont colocalisées au niveau de la station de base BS1, celle-ci se trouvant en un point de la cellule couverte par la station de base BS1. On note que par « colocalisées », on entend que les antennes A11, A21 se trouvent au niveau d’un même site (c’est-à-dire ici d’une même station de base). Elles ne sont toutefois pas nécessairement positionnées en un même point géographique (correspondant à une distance idéalement nulle entre les antennes) et peuvent être séparées de quelques centimètres ou de quelques dizaines de centimètres, voire de quelques mètres. Par exemple, les antennes peuvent être espacées d’une distance inferieure à λ/2 où λ désigne la longueur d’onde des signaux émis par les antennes A11, A21 pour communiquer sur le réseau. En variante, elles peuvent être espacées d’une distance supérieure à λ/2. On note que dans un environnement urbain, on se limitera préférentiellement à un espacement inférieur à une distance allant de 3 à 5 mètres ; dans un environnement rural, un espacement supérieur peut être envisagé, les cellules couvrant des zones plus grandes. Dans l’exemple envisagé ici, par souci de simplification, on suppose que les antennes A11 et A21 sont localisées en un même point au sommet d’un pylône de la station de base BS1. Comme illustré par la figure 1a, l’équipement utilisateur UE est repéré par rapport à la première station de base BS1 dans un repère orthonormé O xyz . L’origine O du repère est ici localisée au pied du pylône de la station de base BS1 supportant les antennes A11, A21, au niveau du sol. L’axe O z est vertical (le long du pylône ici, parallèle à celui-ci) et les axes O x et O y définissent un plan horizontal parallèle au sol et perpendiculaire au pylône. Dans l’exemple envisagé ici, le plan (Oxy) se trouve au niveau du sol et est tangent à la surface de la Terre au point O situé au pied du pylône. L’axe Ox coïncide avec la projection de la direction de propagation principale de l’antenne A11 sur ce plan horizontal parallèle au sol. La projection de la direction de propagation principale de l’antenne A21 dans le plan Oxy est obtenue à partir de l’angle de dépointage θb1, comme illustré sur la figure 1b . On désigne ici par direction d’observation de l’équipement utilisateur UE par une antenne de la station de base sa direction vue par cette antenne. On considère aussi le repère O’xyz dont l’origine O’ est situé au sommet du pylône, et plus particulièrement au point où se trouve l’antenne A11 (et dans l’exemple envisagé ici, les deux antennes sont supposées colocalisées en un même point par souci de simplification, comme mentionné ci-avant). Par exemple, l’équipement utilisateur UE est repéré via des angles d’un système de coordonnées sphériques (θ11, φ11) dans le repère O’xyz en question, et par la distance r11 représentant la projection dans le plan horizontal Oxy de la distance de l’équipement utilisateur UE par rapport à l’origine O’ du repère O’xyz. On suppose ici par souci de simplification que l’équipement utilisateur UE se trouve dans le plan horizontal Oxy, au niveau du sol (on néglige la hauteur à laquelle se trouve l’équipement utilisateur par rapport au sol). En d’autres termes, sur la figure 1a , r11 représente la distance de l’équipement utilisateur UE par rapport au pied du pylône situé en O, supportant les antennes A11, A21, et θ11 et φ11 représentant respectivement la longitude et la latitude de l’équipement utilisateur UE dans le repère O’ xyz . L’angle θ11 est ainsi défini via la projection du vecteur joignant l’origine du repère 0’ xyz à l’équipement utilisateur UE dans le plan horizontal O xy , et mesuré par rapport à la projection de la direction principale de propagation de l’antenne A11 dans ce plan (qui coïncide avec Ox). Les coordonnées (r11, θ11, φ11) définissent la position relative de l’équipement utilisateur UE par rapport à la première station de base BS1, c’est-à-dire en prenant la position de la première station de base BS1 comme référence. Pour obtenir la position absolue de l’équipement utilisateur UE dans un autre référentiel, par exemple le référentiel terrestre, le méridien de Greenwich, etc., il convient de tenir compte de la position de la première station de base BS1 dans ce référentiel. Cette position est connue du réseau de télécommunications auquel appartient la première station de base BS1. On note par ailleurs que l’équipement utilisateur UE peut être situé à une autre altitude que celle du niveau du sol ; par exemple, l’équipement utilisateur UE peut se trouver sur un bateau en mer et être situé au niveau de la mer, le niveau de la mer et le niveau du sol où se trouve le pied du pylône de la station de base BS1 n’étant pas nécessairement les mêmes. Dans ce cas, on pourra soit considérer un repère Oxyz situé à l’altitude de l’équipement utilisateur UE (typiquement au niveau de la mer dans l’exemple envisagé ci-avant) et on aura dans ce repère une altitude z1 de l’équipement utilisateur égale à 0, ou bien conserver la même définition pour le repère (Oxyz) (i.e. origine O au sol, au bas du pylône) et considérer que l’altitude z1 de l’équipement utilisateur dans ce repère est celle du niveau de la mer dans ce repère. Les antennes A11, A21 émettent des faisceaux radiofréquences selon un angle d’inclinaison φt1 (ou angle de « tilt » en anglais) correspondant à un angle de dépointage (une latitude de dépointage ici) de leur diagramme de rayonnement par rapport au plan horizontal O’xy (et donc 0xy ) . De telles antennes sectorielles A11, A21 sont par exemple adaptées à une mise en œuvre dite SISO (pour « Single Input Single Output » en anglais) du réseau de radiocommunications considéré. Sur la figure 1a , les antennes A11 et A21 se trouvent à une hauteur H qui correspond à la distance entre les antennes et le sol. Autrement dit, le plan O’xy se trouve à une hauteur H par rapport au plan Oxy. On présente désormais, en relation avec la figure 2a une antenne comprenant une matrice 200 d’éléments rayonnants 200er (désignée parfois également par réseau d’éléments rayonnants), par exemple du type à balayage électronique, pouvant équiper la station de base BS1 de la figure 2b selon un autre mode de réalisation de l’invention. Plus particulièrement, dans cet exemple, la matrice 200 comprend N H =8 éléments rayonnants 200er dans la direction horizontale (i.e. selon la direction O’ y du repère dans le cas présent) et N V =12 éléments rayonnants 200er dans la direction verticale (i.e. selon la direction 0 z du repère dans le cas présent). On note que dans le présent mode de réalisation, l’origine O’ du repère O’ xyz est localisée au niveau de l’extrémité basse de la matrice 200, comme illustré par la figure 2b. La direction O’x est perpendiculaire au plan de la matrice d’éléments rayonnants et coïncide avec une direction principale de propagation (ou de façon équivalente de rayonnement) de cette matrice 200. Les centres respectifs de deux éléments rayonnants 200er consécutifs sont espacés d’une distance d V dans la direction verticale, et d’une distance d H dans la direction horizontale. En d’autres termes, le pas de la matrice 200 est d V dans la direction verticale et d H dans la direction horizontale. Dans d’autres modes de réalisation, d’autres nombres N H et N V d’éléments rayonnants 200er sont considérés. En relation avec les figures 2a et 2b , une telle matrice 200 d’éléments rayonnants est capable d’émettre différents faisceaux radiofréquences pointant chacun dans une direction de propagation voulue autour de la direction de propagation principale de la matrice 200. Plus particulièrement, des lois de pondération (en amplitude et/ou en phase) de chaque élément rayonnant 200er doivent être mises en œuvre. Des exemples de telles lois sont donnés ci-dessous en relation avec la description de la figure 4 . On considère par exemple une matrice 200 d’éléments rayonnants 200er tel que spécifié dans le document 3GPP TR 37.842 V.13.2.0. Une telle matrice 200 est par exemple adaptée à une mise en œuvre dite MIMO (pour « Multiple-Input Multiple-Output » en anglais) du réseau de radiocommunications considéré. En relation avec la figure 2B, la matrice d’éléments rayonnants de la figure 2a est située au sommet du pylône de la station de base qu’elle équipe. Dans la suite, on désigne par repère de la station de base, comme dans le cas SISO, le repère Oxyz placé au niveau du sol, au pied de ce pylône. On présente maintenant, en relation avec la figure 3 un exemple d’architecture d’un système 10 de localisation d’un utilisateur d’un équipement terminal UE dans un réseau de télécommunications RT selon un mode de réalisation de l’invention. Un tel réseau comprend la station de base BS1, qui vient d’être décrite en relation avec les figures 1a et 1b et une deuxième station de base BS2, présentant des caractéristiques identiques ou similaires à celles décrites pour la première station de base BS1. Néanmoins, on note que les antennes de la deuxième station de base BS2 peuvent avoir des hauteurs différentes de celles des antennes de la première station de base BS1. Le système 10 comprend en outre un dispositif 100 configuré pour géolocaliser l’utilisateur de l’équipement terminal UE dans le réseau RT. Dans un mode de réalisation, un tel dispositif comprend un module OBT. Pij d’obtention d’une puissance mesurée par ledit équipement utilisateur pour chacun des faisceaux radiofréquences émis par la première station de base selon la première et la deuxième direction de propagation et chacun des faisceaux radiofréquences émis par la deuxième station de base selon la troisième et quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; un module OBT. IPR de calcul d’au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative pour la première station de base et deuxième information de puissance relative pour la deuxième station de base, la première puissance étant associée pour la première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de propagation; et un module DET. POS de détermination d’une position dudit équipement utilisateur à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction d’une direction d’observation de l’équipement utilisateur par ladite antenne. Le dispositif 100 met ainsi en œuvre le procédé de géolocalisation selon l’invention qui sera décrit plus en détails en relation avec les figures 4 et 6 . Dans l’exemple de la figure 3 , le dispositif 100 est logé dans le réseau lui-même et intégré à la première station de base BS1, laquelle comprend, classiquement, des mémoires MEM associées à un processeur CPU. Les mémoires peuvent être de type ROM (de l’anglais « Read Only Memory ») ou RAM (de l’anglais « Random Access Memory ») ou encore Flash. La première station de base BS1 comprend en outre un module TX configuré pour piloter l’émission/réception de faisceaux radiofréquences par sa ou ses antennes (non représentées) selon le type de technologie utilisé. En variante, le dispositif 100 pourrait être intégré à la deuxième station de base BS2 ou à un équipement nœud du réseau RT. Selon un autre mode de réalisation, il est intégré à l’équipement terminal UE. On présente désormais, en relation avec la figure 4 les étapes du procédé de géolocalisation selon un mode de réalisation de l’invention dans l’environnement illustré par la figure 5 . Lors d’une étape E300 , la puissance d’une pluralité de faisceaux radiofréquences émis par la première station de base BS1 selon différentes directions de propagation est obtenue. Plus particulièrement, une telle puissance est mesurée par l’équipement utilisateur UE. De la sorte, un jeu correspondant de puissances mesurées est obtenu. On note que les faisceaux radiofréquences dont l’équipement utilisateur UE mesure la puissance ne sont pas nécessairement des faisceaux radiofréquences qui ont été émis par la station de base BS1 à destination de l’équipement utilisateur. En effet, à titre illustratif, dans la configuration SISO représentée aux fig. 1a -1c, l’équipement utilisateur UE, du fait de sa position le cas échéant dans l’un des secteurs couverts par la station de base BS1, reçoit les faisceaux radiofréquences émis par l’antenne de la station de base BS1 couvrant ce secteur, par exemple la première antenne A11 selon une première direction de propagation. Toutefois, il est également en mesure de recevoir des faisceaux radiofréquences émis par une autre antenne de la station de base BS1 couvrant un secteur adjacent à celui dans lequel il se trouve, par exemple de la deuxième antenne A21 selon une deuxième direction de propagation. De même, dans une configuration MIMO s’appuyant par exemple sur une antenne comprenant une matrice d’éléments rayonnants tel que représentée sur les figures 2a et 2b, plusieurs faisceaux radiofréquences peuvent être émis simultanément dans plusieurs directions de propagation distinctes, qui ne coïncident pas nécessairement avec la direction dans laquelle se trouve l’équipement utilisateur UE par rapport à l’antenne en question. Celui-ci peut toutefois être en mesure de recevoir ces faisceaux radiofréquences même s’ils ne lui sont pas directement destinés. Pour mettre en œuvre l’invention, on suppose que les faisceaux radiofréquences émis par la première station de base peuvent être mesurés par l’équipement utilisateur et qu’il est capable de distinguer cette puissance d’une puissance de bruit. Par exemple, l’équipement utilisateur UE identifie quelle antenne A11, A21 a émis le faisceau pour lequel il mesure la puissance à partir d’informations véhiculées par le faisceau en question dans des canaux de signalisation, dits canaux communs. De telles informations sont par exemple des symboles pilotes, connus en soi. Ces informations lui sont transmises par la station de base de façon spontanée ou à sa demande expresse. Selon l’invention, l’équipement utilisateur UE reçoit aussi de la deuxième station de base BS2 un premier faisceau radiofréquences émis selon une troisième direction de propagation et un deuxième faisceau radiofréquences émis selon une quatrième direction de propagation. Ces faisceaux émis par la deuxième station de base peuvent être reçus simultanément ou consécutivement de ceux émis par la première station de base BS1. Par ailleurs, comme illustré par la figure 1a , on suppose que les première et deuxième antennes de la première station de base BS1 sont situées à une même hauteur H et émettent des faisceaux radiofréquences selon un même angle d’inclinaison, ou premier tilt φt1 (i.e. une même latitude, non représentée sur la fig. 5 ) par rapport au plan horizontal O’xy (et donc par rapport au plan horizontal O xy) . De même, on suppose que les première et deuxième antennes de la deuxième station de base BS2, A12 et A22, émettent des faisceaux radiofréquence selon un même angle d’inclinaison, ou deuxième tilt φt2 (i.e. une même latitude, non représentée sur la fig. 5) par rapport à un plan horizontal O 2 ’x 2 y 2 défini par rapport à la deuxième station de base BS2, de manière identique à la façon dont a été défini précédemment le plan O’xy par rapport à la première station de base BS1. Autrement dit, O 2 ’ est situé au sommet du pylône de la station de base BS2, au point où se trouvent les antennes A12 et A22 de la station de base BS2 (qu’on suppose colocalisées en un même point par souci de simplification). L’axe (O 2 ’x 2 ) coïncide avec la projection de la direction de propagation principale de l’antenne A12 sur un plan parallèle au sol, situé au niveau du sol. Par souci de simplification ici, on fait l’hypothèse que les stations de base BS1 et BS2 sont suffisamment proches l’une de l’autre pour négliger le rayon de courbure de la Terre entre ces deux stations de base de sorte que les plans (O’xy) et (O 2 ’x 2 y 2 ) sont supposés confondus. L’inventeur a constaté que cette hypothèse restait valable y compris lorsque les stations de base considérées sont éloignées de quelques dizaines de kilomètres. φt1 et φt2 peuvent avoir des valeurs identiques ou distinctes. On note que, selon l’invention, il n’est pas nécessaire de connaître la hauteur de la ou des antennes de la deuxième station de base BS2. Dans certains modes de réalisation, les puissances mesurées par l’équipement utilisateur UE sont remontées par ce dernier au réseau de radiocommunications (e.g. via une transmission à la première station de base BS1 ou alternativement à la deuxième station de base BS2 ou encore à un équipement nœud EN du réseau RT). C’est le cas par exemple lorsque le dispositif 100 mettant en œuvre le présent procédé de géolocalisation est logé dans le réseau lui-même (e.g. dans un nœud du réseau ou dans l’une des deux stations de base BS1, BS2). Dans ces modes de réalisation, la première station de base BS1 reçoit de l’équipement utilisateur la puissance mesurée par l’équipement utilisateur UE pour chaque faisceau radiofréquences de la pluralité de faisceaux radiofréquences émis par la première station de base BS1 selon la première et la deuxième direction de propagation (formant un « premier jeu (P11,P21) de puissances mesurées » au sens de l’invention pour la première station de base BS1) et la transmet au dispositif 100. En parallèle, la deuxième station de base BS2 reçoit la puissance mesurée par l’équipement utilisateur UE pour chaque faisceau radiofréquences de la pluralité de faisceaux radiofréquences émis par la deuxième station de base BS1 selon la troisième et la quatrième direction de propagation (formant un « deuxième jeu (P12,P22) de puissances mesurées » au sens de l’invention pour la deuxième station de base BS2) et la transmet au dispositif 100. Par exemple, si le dispositif 100 est intégré dans la première station de base BS1, les puissances mesurées (P12, P22) sont transmises par la deuxième station de base BS2 à la première station de base BS1 via le réseau ou par l’intermédiaire de l’équipement utilisateur UE. Dans d’autres modes de réalisation, le présent procédé de géolocalisation est mis en œuvre directement dans l’équipement utilisateur UE. Dans ce dernier cas, le dispositif 100 est intégré à l’équipement utilisateur UE et les puissances mesurées sont stockées localement. Lors d’une étape E310 , au moins une information de puissance relative, représentative d’un rapport (lorsque les puissances sont exprimées en unité naturelle) ou d’une différence (lorsque les puissances sont exprimées en unité logarithmique) entre deux puissances ((P11, P21) pour la station de base BS1) et (P12, P22) pour la station de base BS2) du jeu de puissances mesurées associées à deux faisceaux radiofréquences correspondants, émis par la même station de base, est calculée. Ainsi, le dispositif 100 obtient une première information de puissance relative M11 à partir du jeu de puissances mesurées pour la première station de base BS1 et une deuxième information M12 de puissance relative à partir du jeu de puissances mesurées pour la deuxième station de base BS2. Dans le cas d’une configuration SISO, lorsque les antennes d’une station de base se trouvent en un même point géographique, les effets du fast fading s’éliminent quand on fait le rapport ou la différence des puissances reçues par l’utilisateur. On note que dans certains modes de réalisation, les antennes peuvent ne pas être exactement colocalisées ; dans ce cas, les puissances mesurées par l’équipement utilisateur UE peuvent être moyennées sur une durée déterminée pour éliminer les effets du fast fading. Par exemple, elles sont collectées à une fréquence d’acquisition déterminée (par exemple toutes les millisecondes) et sont moyennées par l’équipement utilisateur UE sur une période déterminée. La durée de cette période peut être déterminée en fonction de différents paramètres, comme par exemple la mobilité éventuelle de l’équipement utilisateur UE et, le cas échéant, sa vitesse, etc. Cette moyenne peut être réalisée à l’aide d’une fenêtre glissante de longueur égale à la période déterminée envisagée. Par exemple, l’inventeur a déterminé que pour une fréquence de 1 GHz, une moyenne effectuée sur une durée de 50 ms de mesures acquises toutes les millisecondes est suffisante pour de nombreuses antennes classiquement utilisées pour obtenir une estimation précise de la position de l’équipement utilisateur. Cette moyenne permet de s’affranchir des phénomènes de variations rapides (ou « fast-fading » en anglais) des canaux de propagation qui peuvent différer légèrement d’une antenne à l’autre lorsque celles-ci sont éloignées de quelques centimètres ou quelques dizaines de centimètres notamment. Néanmoins, la distance entre antennes d’une station de base doit être relativement faible pour que ce moyennage élimine l’effet de la distance r avec le récepteur avec une précision suffisante. A cet égard, on note que la mobilité éventuelle de l’équipement utilisateur UE impose des contraintes plus importantes sur la rapidité des mesures, que dans le cas statique. Dans le cas d’une configuration MIMO, on peut considérer une seule mesure de la puissance, l’ensemble des éléments rayonnants étant utilisé pour émettre la pluralité de faisceaux radiofréquences. En E320 , la position de l’équipement utilisateur UE est déterminée à partir, d’une part, des informations de puissance relative et, d’autre part, pour chacune de ces informations de puissance relative et pour chaque faisceau radiofréquences associé à l’information de puissance relative (i.e. pour chacun des deux faisceaux radiofréquences de directions de propagation différentes dont la puissance est à la base de l’information de puissance relative en question), d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction d’une direction d’observation de l’utilisateur, du faisceau radiofréquences associé à l’information de puissance relative. En relation avec les figures 4 et 6 , des exemples de mise en œuvre d’un tel procédé de géolocalisation sont maintenant détaillés en reconsidérant les configurations des première et deuxième stations de base BS1, BS2 discutées ci-dessus en relation, d’une part, avec les figures 1a et 1b et, d’autre part, avec les figures 2a et 2b . On se réfère en outre aux schémas des figures 5 et 7 qui illustrent géométriquement les principes de l’invention. On détaille d’abord le mode de réalisation de la figure 6, pour lequel une technologie de type SISO est utilisée. Lors de l’étape E300 déjà décrite, le dispositif 100 obtient par exemple une première P11(UE) et une deuxième P21(UE) puissances mesurées correspondant respectivement à un premier faisceau radiofréquence émis par la première antenne A11 selon une première direction de propagation et à un deuxième faisceau radiofréquence émis par la deuxième antenne A21 selon une deuxième direction de propagation. En effet, quand bien même l’équipement utilisateur UE se situe dans un secteur donné, il reçoit des puissances des différentes antennes de la station de base couvrant le site. Comme déjà évoqué, on suppose que les première A11 et deuxième A21 antennes émettent des faisceaux radiofréquences selon des directions de propagation différentes mais selon un même angle d’inclinaison ou premier tilt φt1 par rapport au plan horizontal O’ xy . De la sorte, la première puissance P11(UE) mesurée au niveau de l’équipement utilisateur UE s’exprime de manière générale selon l’expression : P11(UE) = K.P0.r -eta .G11(θ11, φ11).X BS (UE).Y BS (UE) (Eq. 1) où K est une constante, r désigne la distance séparant l’équipement utilisateur UE de l’antenne A11, eta est un facteur d’évanouissement (aussi appelé facteur de « pathloss ») modélisant l’atténuation de propagation, P0 est la puissance émise par l’antenne A11 sur le faisceau considéré, G11(θ11, φ11) est le gain de l’antenne A11 rayonné sur le faisceau considéré dans la direction (θ11, φ11) (qui définit la direction d’observation de l’équipement utilisateur UE par l’antenne A11 au sens de l’invention), X BS (UE) est un paramètre représentant les évanouissements rapides du canal de propagation entre l’antenne A11 et l’équipement utilisateur UE, et Y BS (UE) est un paramètre représentant les effets de masque (ou « shadowing » en anglais) du canal de propagation entre l’antenne A11 et l’équipement utilisateur UE. Une telle modélisation est connue en soi et n’est pas décrite plus en détails ici. On note que le modèle utilisé en équation (1) reliant la puissance reçue par l’équipement utilisateur à la direction d’observation (θ11, φ11) est obtenu en considérant que l’antenne A11 et l’équipement utilisateur UE sont en visibilité directe (ou LOS pour « line-of-sight » en anglais) : dans cette configuration, le faisceau émis par l’antenne A11 selon la direction (θ11, φ11) est reçu selon la même direction par l’équipement utilisateur UE. Cette configuration en visibilité directe se présente avec une forte probabilité lorsque l’équipement utilisateur se trouve sur un bateau en mer ou dans un environnement rural. Lorsque l’équipement utilisateur se trouve dans un environnement urbain, il peut également être considéré, avec une probabilité non nulle, comme étant dans une configuration en visibilité directe. Toutefois, ce modèle peut encore être utilisé lorsque l’équipement utilisateur UE n’est pas considéré comme étant en visibilité directe avec l’antenne A11 (on parle aussi de configuration NLOS pour « Non Line of Sight » en anglais), par exemple en raison de la présence d’obstacles entre l’antenne A11 et l’équipement utilisateur UE. Dans ce cas en effet, il se peut que le faisceau émis par l’antenne A11 selon la direction (θ11, φ11) soit reçu par l’équipement utilisateur UE selon une direction légèrement différente, par exemple (θ11+δθ11,φ11+δφ11). Le modèle décrit ci-dessus et la localisation de l’équipement utilisateur UE qui en résulte conformément à l’invention peuvent alors s’avérer moins précis dans cette configuration. De manière alternative, on peut utiliser en configuration NLOS, un modèle plus précis décrivant la puissance reçue au niveau de l’utilisateur en fonction de la direction d’observation de l’utilisateur par l’antenne, estimé par exemple via des simulations numériques ou expérimentalement.De même, la deuxième puissance P21(UE) mesurée au niveau de l’équipement utilisateur UE s’exprime de manière générale selon l’expression : P21(UE) = K.P0.r -eta .G21(θ21, φ21).X BS (UE).Y BS (UE) (Eq. 2) où G21(θ21, φ21) est le gain de l’antenne A21 rayonné sur le faisceau considéré dans la direction (θ21, φ21) de l’équipement utilisateur UE (qui définit la direction d’observation de l’équipement utilisateur UE en visibilité directe par l’antenne A21 au sens de l’invention). En pratique, si l’on suppose que les antennes A11 et A21 ont un même diagramme de rayonnement et émettent des faisceaux radiofréquences selon un même angle d’inclinaison ou tilt φt1, on a ainsi θ21 = θb 1 - θ11 et φ21 = φ11. On suppose également la même puissance émise P0 pour l’ensemble des antennes de la station de base BS1. Par ailleurs, comme mentionné précédemment par souci de simplification, on suppose ici que les antennes A11 et A21 sont colocalisées en un même point au sommet du pylône de la station de base BS1 : il en résulte que les termes représentant l’atténuation de propagation, les évanouissements rapides (ou « fast fading ») et les effets de masquage sont identiques dans les équations (Eq. 1) et (Eq. 2). Les autres configurations possibles de disposition des antennes A11 et A21, et leur impact sur les évanouissements rapides et les effets de masquage, ont déjà été évoquées précédemment. De la sorte, lors de la mise en œuvre de l’étape E310, la première information de puissance relative M11 associée aux premier et deuxième faisceaux radiofréquences s’exprime uniquement en fonction des gains G11(θ11, φ11) et G21(θb 1 -θ11, φ11). En effet, à partir des équations (Eq. 1) et (Eq. 2), on peut écrire : M11 = P11(UE)/P21(UE) = G11(θ11, φ11)/G21(θb 1 -θ11, φ11) (Eq. 3lin) ou, alternativement, en décibels : M11 dB = (G11(θ11, φ11) dB - G21(θb 1 -θ11, φ11)) dB (Eq. 3dB) On note que lorsque les puissances d’émission des antennes A11 et A21 sont différentes, les mêmes équations peuvent être obtenues à une constante (multiplicative ou additive) près. Cette constante caractérise le rapport (en linéaire) ou la différence (en dB) entre les deux puissances d’émission des antennes considérées. De la sorte, ayant la connaissance d’un modèle (e.g. analytique ou obtenu par mesure de l’antenne en question) de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction d’une direction d’observation de l’équipement utilisateur par l’antenne considérée, du faisceau radiofréquences considéré (ici un modèle pour G11(θ1, φ1) et pour G21(θb 1 -θ11, φ11) en visibilité directe), il est possible de remonter tout ou partie des coordonnées caractérisant la position de l’équipement utilisateur UE en résolvant l’équation (Eq. 3lin) ou l’équation (Eq. 3dB) lors de la mise en œuvre de l’étape E320. En d’autres termes, lors de l’étape E320, la détermination de la position de l’équipement utilisateur UE met en œuvre, pour au moins une information de puissance relative donnée, la résolution d’une équation (équation (Eq. 3lin) ou équation (Eq. 3dB)) dont les membres sont fonction, d’une part, de l’information de puissance relative donnée et, d’autre part, d’une valeur attendue de l’information de puissance relative donnée. Plus particulièrement, la valeur attendue en question est fonction, pour chaque faisceau radiofréquences associé à l’information de puissance relative donnée, du modèle de diagramme de rayonnement caractérisant la puissance, en fonction d’une direction d’observation de l’équipement utilisateur par l’antenne émettrice, du faisceau radiofréquences associé à l’information de puissance relative donnée. Par exemple, reconsidérant un modèle d’antennes telles que décrites dans le document Report ITU-R M.2135-1 de l’ITU-R, intitulé « Guidelines for évaluation of radio technologies for IMTAdvanced » de décembre 2009, on peut écrire : (G11(θ11, φ11)) dB = -12 (θ11/ θ 3dB ) 2 - 12 ((φ11 – φt1) / φ 3dB ) 2 (Eq. 4) où θ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement dans le plan de définition de l’angle θ11 (i.e. dans le plan horizontal Oxy ou de façon équivalente O’xy), φ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement dans le plan de définition de l’angle φ11. On considère un plan vertical défini par l’axe vertical 0 z et par la direction d’observation de l’équipement utilisateur UE considérée par l’antenne A11, c’est-à-dire la direction de l’équipement utilisateur UE dans l’équation (Eq. 4) qui est donnée par les angles θ11 et φ11. Ce plan est parfois appelé plan méridien. φt1 est représentatif de l’angle d’inclinaison des faisceaux radiofréquences ou premier tilt émis par les antennes A11 et A21 par rapport au plan horizontal 0’ xy . On peut ainsi écrire de la même façon : (G21(θ21, φ21)) dB = -12 (θ21/ θ 3dB ) 2 - 12 ((φ21 – φt1) / φ 3dB ) 2 (Eq. 5) Sachant que θ21 = θb 1 - θ11 et φ21 = φ11 comme détaillé ci-dessus, on obtient en injectant les équations (Eq. 4) et (Eq. 5) dans l’équation (Eq. 3dB) : M11 dB = 12 θb 1 ²/ θ 2 3dB - 24 θb 1 .θ11/θ 2 3dB (Eq. 6) De la sorte, l’angle θ11 repérant l’équipement utilisateur UE dans le repère O’xyz considéré centré sur la première antenne A11 de la station de base BS1 est donné par : θ11 = - (M11 dB /24 ) θ 2 3dB / θb 1 + θb 1 /2 (Eq. 7) En d’autres termes, comme illustré par la figure 6 , au cours d’une première sous-étape E321 de l’étape E320, l’angle θ11 est déterminé à partir, d’une part, de la première information de puissance relative M11 (exprimée en unité logarithmique dans l’équation (Eq. 7)) et, d’autre part, du modèle de diagramme de rayonnement caractérisant la puissance, en fonction d’une direction d’observation, en l’espèce la direction de l’équipement utilisateur UE par rapport aux antennes A11 et A21 respectivement, des premier et deuxième faisceaux radiofréquences (ici G11(θ11, φ11) et G21(θ21, φ21)). Cependant une incertitude existe sur le signe de θ11 selon la valeur de M11 dB (ou, alternativement, M11) mesurée. Deux valeurs de θ11 sont en effet possibles, le signe de l’angle θ11 dans le repère Oxyz dépendant du signe de l’information de puissance relative exprimée en décibels M11 dB . Elles définissent donc deux plans méridiens PM11, PM21 dans le repère Oxyz de la première station de base. La position de l’utilisateur de l’équipement utilisateur UE se trouve dans l’un de ces deux plans, situés d’un côté ou de l’autre de la direction principale de la première antenne A11 de la première station de base BS1. Une solution pour lever cette indétermination sera détaillée ci-après. Avantageusement, selon ce mode de réalisation de l’invention, on obtient en E310 une troisième puissance mesurée pour la première station de base BS1 lors d’une sous-étape E301 de l’ étape E300 . Cette troisième puissance a été mesurée pour un faisceau radiofréquence émis par la première antenne A11 de la première station de base BS1 avec un troisième tilt φt3, distinct du premier tilt φt1. Avantageusement, ce changement de tilt et cette nouvelle émission de faisceaux sont déclenchés par la première station de base BS1. En d’autres termes la première station de base BS1 est équipée d’antennes qui sont configurées électroniquement (par exemple à distance) pour émettre selon deux valeurs distinctes de tilt, par exemple alternativement. Dans d’autres modes de réalisation, ce changement de tilt et cette émission supplémentaire de faisceaux par la première station de base avec le troisième tilt, sont déclenchés par l’émission d’une requête en provenance de l’équipement terminal UE. C’est en particulier le cas lorsque le dispositif de géolocalisation 100 est embarqué dans l’équipement terminal UE. Dans tous les cas, les mesures des puissances reçues selon un premier tilt et un troisième tilt sont effectuées à des instants temporels très proches l’un de l’autre. Selon ces modes de réalisation, au moins une troisième information de puissance relative M’11, représentative d’un rapport (lorsque les puissances sont exprimées en unité naturelle) ou d’une différence (lorsque les puissances sont exprimées en unité logarithmique) entre la première puissance P11 du jeu de puissances mesurées en E300 et associées au faisceau radiofréquence émis par la première antenne A11 de la première station de base BS1 avec le premier tilt φt1 et la puissance P11’ du jeu de puissances mesurées et associées au faisceau radiofréquence émis par la première antenne A11 de la première station de base BS1 avec le troisième tilt φt3, est calculée en E310. Au cours d’une deuxième sous-étape E322 de l’étape E320 , on procède de la même manière qu’en E321 pour obtenir, à partir des informations de puissance relative obtenues en E310 à partir des puissances des faisceaux radiofréquences mesurées pour la deuxième station de base BS2. On fait d’abord l’hypothèse que les première et deuxième antennes de la deuxième station de base BS2 émettent avec un deuxième tilt φt2, égal au premier φt1 : φt2= φt1. On peut donc écrire, dans un repère orthonormé O 2 ’x 2 y 2 z 2 de la première antenne de la deuxième station de base (défini de façon identique au repère O’xyz pour la première antenne A11 de la première station de base BS1, en prenant, comme indiqué précédemment, O 2 ’ au sommet du pylône de la deuxième station de base BS2, là où sont colocalisées les antennes A12 et A22 de la deuxième station de base BS2, et en définissant l’axe O 2 ’x 2 par rapport à la direction principale de propagation de l’antenne A12) : M12 = P12(UE)/P22(UE) = G12(θ12, φ12)/G22(θb 2 -θ12, φ12) (Eq. 3lin bis) ou, alternativement, en décibels : M12 dB = (G12(θ12, φ12) - G22(θb 2 -θ12, φ12)) dB (Eq. 3dB bis) où G12(θ12, φ12) est le gain de la première antenne A12 rayonné sur le faisceau considéré dans la direction (θ12, φ12) d’observation de l’équipement utilisateur UE par l’antenne A12 et G22(θb 2 -θ12, φ12) le gain de la deuxième antenne A22 rayonné sur le faisceau considéré dans la direction (θ22, φ22) d’observation de l’équipement utilisateur UE par l’antenne A22. Comme pour M11, on considère ici des gains reflétant une visibilité directe de l’équipement utilisateur UE par les antennes A12 et A22. Conformément à l’équation 4 précédente, on a : (G12(θ12, φ12)) dB = -12 (θ12/ θ 3dB ) 2 - 12 ((φ12 – φt1) / φ 3dB ) 2 (Eq. 8) où θ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement dans le plan de définition de l’angle θ12 (i.e. dans le plan horizontal 0 2 x 2 y 2 , qui est, comme souligné précédemment, confondu avec le plan O xy, comme illustré par les figures 5 et 7), φ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement dans le plan de définition de l’angle φ12. On considère le plan vertical (ou plan méridien) défini par l’axe vertical O 2 z et par la direction d’observation considérée, c’est-à-dire la direction de l’équipement utilisateur UE dans l’équation (Eq. 8). φt1 est représentatif de l’angle d’inclinaison des faisceaux radiofréquences ou premier tilt émis par les antennes A12 et A22 par rapport au plan horizontal 0’ 2 xy et donc par rapport au plan horizontal O 2 xy (pour rappel, φt2= φt1 dans l’exemple envisagé ici). On peut ainsi écrire de la même façon : (G22(θ22, φ22)) dB = -12 (θ22/ θ 3dB ) 2 - 12 ((φ22 – φt1) / φ 3dB ) 2 (Eq. 9) Sachant que θ22 = θb 2 – θ21 et φ22 = φ12 comme détaillé ci-dessus, on obtient en injectant les équations (Eq. 8) et (Eq. 9) dans l’équation (Eq. 3dB bis) : M12 dB = 12 θb 2 ²/ θ 2 3dB - 24 θb 2 .θ12/θ 2 3dB (Eq. 9) De la sorte, l’angle θ12 repérant l’équipement utilisateur UE dans le repère O’ 2 xyz considéré centré sur la première antenne A12 de la deuxième station de base BS2 est donné par : θ12 = - (M12 dB /24 ) θ 2 3dB / θb 2 + θb 2 /2 (Eq. 10) Comme pour la première station de base BS1, la longitude θ12 peut prendre deux valeurs distinctes et définit ainsi deux plans méridiens PM12, PM22, comme illustré par la figure 5. A l’issue des sous-étapes E341 et E342, on a donc 4 plans verticaux PM11, PM21, PM12, PM22 qui se coupent en 4 droites verticales D11, D21, D’11, D’21. La position de l’équipement utilisateur UE se trouve à l’intersection de l’une de ces droites avec le plan horizontal correspondant à son altitude, autrement dit, dans l’exemple envisagé ici d’un équipement utilisateur au niveau du sol, avec le plan (Oxy). Il convient de noter que pour que les équations (6) et (9) qui viennent d’être décrites, s’appliquent, il faut que conditions suivantes soient vérifiées en visibilité directe, comme spécifié dans le document Report ITU-R M.2135-1 de l’ITU-R, intitulé « Guidelines for évaluation of radio technologies for IMTAdvanced » de décembre 2009, déjà cité : Gij = -min(12 (θij/θ3dB)2 + 12 ((φij – φt1) / φ3dB) 2 , Am) (Eq. 11) avec i l’indice de l’antenne, entier égal à 1 ou 2, j l’indice de la station de base, entier égal à 1 ou 2 et Am une constante qui caractérise la valeur minimale de gain de l’antenne. Les angles de longitude et de latitude doivent satisfaire les conditions suivantes : θij φij 12 (θij/θ3dB)2 + 12 ((φij – φtk) / φ3dB)2 Si ces conditions ne sont pas vérifiées, l’équation (11) de diagrammes de rayonnement ne s’applique pas. Il conviendrait alors de choisir une troisième antenne de la station de base, si elle en comporte une troisième, ou une autre station de base, pour laquelle cette équation est valable. Cela ne pose pas en soi de difficultés pour l’équipement utilisateur UE de mesurer les puissances reçues d’une autre station de base, les stations de base d’un réseau de télécommunications étant généralement suffisamment proches les unes des autres pour que l’équipement utilisateur UE reçoive des signaux de plusieurs stations de base. Dans la suite, on considère que, comme c’est généralement le cas, ces conditions sont satisfaites pour les stations de base BS1 et BS2 et leurs antennes respectives A11, A21, et A12, A22. Comme illustré par la figure 7 , la première antenne A12 de la deuxième station de base BS2 peut être située à une hauteur H2 différente de la hauteur H des antennes A11 et A21, et l’angle de latitude φ12 de l’équipement utilisateur UE dans le repère de la deuxième station de base BS2 est, dans le cas général, différent de l’angle de latitude φ11 dans le repère de la première station de base BS1. On note toutefois que leur connaissance n’est pas requise pour mettre en œuvre l’invention. On détaille maintenant une première détermination E323 de l’angle de latitude φ11 de l’équipement utilisateur UE dans le repère Oxyz de la première station de base, selon un mode de réalisation de l’invention. Pour ce faire, lors de la mise en œuvre de l’étape E300, le dispositif 100 obtient une troisième puissance P’11(UE) mesurée correspondant à un troisième faisceau radiofréquences émis par l’antenne sectorielle A11 selon un deuxième angle d’inclinaison ou tilt φt3, par rapport au plan horizontal 0’ xy (et donc par rapport au plan Oxy), différent du premier angle d’inclinaison ou tilt φt1. De la sorte, lors de la mise en œuvre de l’étape E310, une troisième information de puissance relative M’11 associée à ce troisième faisceau radiofréquences et au premier faisceau radiofréquences précité est calculée. Ainsi, lors d’une sous-étape E323 de l’étape E320, l’angle φ11 repérant l’équipement utilisateur UE dans le repère O’xyz centré sur l’antenne A11 de la première station de base BS1 est déterminé à partir, d’une part, de la deuxième information de puissance relative M’11 et, d’autre part, du modèle de diagramme de rayonnement caractérisant la puissance, en fonction d’une direction d’observation, des premier et troisième faisceaux radiofréquences. Par exemple, à partir du modèle de diagramme de rayonnement d’antenne considéré ci-dessus, on peut écrire le gain G’11(θ11, φ11) de l’antenne A11 dans la direction (θ11, φ11) de l’équipement utilisateur UE lorsque l’antenne A11 émet selon un angle d’inclinaison égal à φt3. On obtient : (G’11(θ11, φ11)) dB = -12 (θ11/ θ 3dB ) 2 - 12 ((φ11 – φt3) / φ 3dB ) 2 (Eq. 13) Selon l’équation 4, on a toujours : (G11(θ11, φ11)) dB = -12 (θ11/ θ 3dB ) 2 - 12 ((φ11 – φt1) / φ 3dB ) 2 (Eq. 4) De la sorte, à partir des équations (Eq. 4) et (Eq. 11) on peut écrire que la deuxième information de puissance relative en décibels, M’11 dB = (G11(θ1, φ1) – G’11(θ1, φ1)) dB , s’exprime théoriquement comme : M’11 dB = - 12 ((φ11 – φt1) / φ 3dB ) 2 + 12 ((φ11 – φt3) / φ 3dB ) 2 (Eq. 14) L’angle φ11 est ainsi obtenu comme étant : φ11 = M’11 dB /24 φ 2 3dB / (φt1– φt3) + (φt1+ φt3)/2 (Eq. 15) Ainsi, l’angle de longitude θ11 et l’angle de latitude φ11 repérant l’équipement utilisateur UE dans le repère Oxyz de la première station de base BS1, sont déterminés de manière simple et précise à partir des faisceaux émis par les antennes A11, A12 de la première station de base BS1 selon au moins deux angles d’inclinaison différents et de ceux émis par les antennes A12, A22 de la deuxième station de base BS2 avec au moins un angle d’inclinaison qui peut être égal à l’un des deux précédents. En E324, on réalise une deuxième détermination de l’angle de latitude φ11 dans le repère O’xyz de la première station de base. Comme illustré par les figures 5 et 7 , cet angle de latitude φ11 définit un cône C centré sur l’origine O’ du repère O’xyz . La position de l’équipement utilisateur UE se trouve à l’intersection de ce cône C avec l’une des quatre droites D11, D21, D’11, D’21 et du plan horizontal correspondant à l’altitude z1 de l’équipement utilisateur UE dans le repère Oxyz (dans l’exemple envisagé ici, l’équipement utilisateur est supposé dans le plan Oxy , autrement dit, z1=0). On note (x1,y1,z1) les coordonnées de l’équipement utilisateur UE dans le plan Oxyz de la première station de base BS1. En relation avec la figure 5 , la première antenne A11 de la première station de base BS1 est située à l’origine O’ du repère O’xyz , donc en (0,0,0). La première antenne A12 de la deuxième station de base BS2 est située en O’ 2 (x0, y0, 0) du repère O’xyz. On note que dans cet exemple elle est située à la même hauteur H2 = H que la première antenne A11 de la première station de base BS1, mais que, comme précédemment évoqué, ce n’est pas une condition de mise en œuvre de l’invention. Pour la première station de base BS1, on peut écrire l’équation des plans verticaux PM11, PM21 passant par l’origine O du repère Oxyz et dirigés selon l’angle θ11 comme suit : y = x tg (θ11) (Eq. 16) De même, pour la deuxième station de base BS2, les plans verticaux PM12, PM22 sont définis dans le repère Oxyz par l’équation suivante : y-y0 = (x-x0) tg (θ12) (Eq.17) Comme illustré par la figure 5 , ces plans se coupent deux à deux en une droite (D11, D21, D’11, D’21) répondant au système d’équations suivant : y = x tg (θ11) et y – y0 = (x-x0) tg (θ12) La troisième coordonnée z est quelconque. L’équipement utilisateur UE étant au niveau du sol, on a z1 = 0 et la droite coupe le plan z = 0. Les coordonnées de l’équipement utilisateur UE dans le repère Oxyz , s’expriment donc comme suit : x1= (y0 - x0. tg (θ12) ) /( tg (θ11) - tg (θ12) ) (Eq. 18) y1 = x1 tg (θ11) z1=0 Les deux valeurs possibles de chaque angle de longitude θ11, θ12, conduisent aux coordonnées de 4 points possibles (x1, y1) sur les 4 droites D11, D21, D11’ et D21’. Or, l’équipement utilisateur UE se situe à une distance r11 de l’origine du repère Oxyz telle que r11 = H / tg(φ11). Comme on a :x1 = r11.cos(θ11), on en déduit donc que : x1 = H / tg(φ 11). cos(θ11) et on peut écrire que : tg(φ 11) = H. cos(θ11). [ tg(θ11) - tg(θ12) ] / ( y0 - x0. tg(θ12) ) (19) On a donc : φ 11 = Arctan ((H / ( y0 - x0. tg(θ12) ) ). cos(θ11). [ tg(θ11) - tg(θ12) ]) (20) Ainsi, on dispose à ce stade de deux expressions différentes de φ 11 données respectivement par les équations 15 et 20. On peut ainsi déterminer en E325 la valeur des angles de longitude θ11, θ12 qui conduisent à la valeur de latitude φ11 la plus proche de celle calculée à partir de l’équation (15). On note que la différence existant le cas échéant entre la valeur de latitude obtenue à partir de l’équation (15) et celle issue de l’équation (20) est due à une imprécision de la mesure des puissances des faisceaux radiofréquence reçus par l’équipement utilisateur UE. Une fois l’incertitude levée sur la valeur des angles de longitude θ11, θ12, on peut déduire en E326 les cordonnées x1 et y1 de l’équipement utilisateur UE à partir des équations (18). Dans l’exemple de la figure 7 , la position de l’équipement utilisateur UE correspond géométriquement à l’intersection entre le cône C, le plan horizontal correspondant à l’altitude de l’équipement utilisateur (ici le plan z=0 dans le repère Oxyz ). Le cône C détermine un arc de cercle AC sur le plan Oxy qui ne passe que par une seule des 4 droites verticales déterminées précédemment. On note qu’il existe une très faible probabilité que cet arc de cercle AC passe par plus d’une des 4 droites. En effet, il faudrait une configuration particulière de positionnement de l’équipement utilisateur par rapport aux stations de base BS1 et BS2, par exemple que la direction principale d’une antenne de la première station de base soit alignée avec celle d’une antenne de la deuxième station de base considérée pour déterminer la position de l’équipement utilisateur UE. Dans ce cas, l’arc de cercle pourrait couper deux droites verticales. On note cependant, même dans cette configuration, les contraintes exprimées par les équations (12) limitent fortement les possibilités d’intersection de l’arc de cercle AC à une seule des quatre droites verticales. Si ce n’était pas le cas, il faudrait alors prendre en compte des mesures de puissance de faisceaux radiofréquences émis par une troisième station de base distincte de BS1 et BS2. Dans l’exemple illustratif envisagé précédemment, on a supposé que l’équipement utilisateur UE se trouvait au niveau du sol, autrement dit z1=0. Bien entendu cette hypothèse n’est pas limitative en soi, et l’homme du métier saurait adapter sans difficulté les équations qui précèdent pour prendre en compte l’altitude de l’équipement utilisateur UE. A titre illustratif, on considère maintenant le cas où l’équipement utilisateur UE est sur un bateau qui navigue en mer. On suppose par ailleurs qu’il se trouve à une hauteur h par rapport au niveau de la mer, qui se trouve à une hauteur H M par rapport au niveau du sol. Autrement dit, l’équipement utilisateur UE se trouve à une hauteur h+H M par rapport au niveau du sol. Si le niveau du sol et celui de la mer H M coïncident, la hauteur relative des antennes des stations de base BS1 et BS2 dans le repère Oxyz par rapport à la hauteur de l’équipement utilisateur UE est H-h, où H désigne la hauteur des antennes dans le repère Oxyz. Dans ce cas, il convient de considérer dans le calcul de la tangente de l’angle φ11 la hauteur relative H-h. On a donc : φ11 = Arctan (((H-h) / ( y0 - x0. tg(θ12) ) ).cos(θ11). [ tg(θ11) - tg(θ12) ]) (20). On note que si le niveau de la mer se trouve en-dessous du niveau du sol, sa hauteur H M est alors négative dans le repère Oxyz. De même, l’équipement utilisateur UE peut également se situer à une hauteur négative par rapport au niveau du sol. Dans tous les cas, il convient de considérer dans l’équation 20, la hauteur relative des antennes des stations de base BS1 et BS2 par rapport à celle de l’équipement utilisateur UE, celle-ci étant égale, dans le cas envisagé ici, à H-H M -h. On remplace donc la valeur de H par H-H M -h dans l’équation 20. Dans l’exemple envisagé ici, on a considéré un repère Oxyz situé au niveau du sol et du bas du pylône de la station de base. Toutefois, suivant le contexte d’application de l’invention, on peut envisager une autre référence pour le repère, par exemple un repère situé au niveau de la mer. Alternativement, on décrit maintenant, en relation avec la figure 7 , un exemple de réalisation de l’invention à partir de deux stations de base ayant une configuration d’antenne de type MIMO. Selon cette technologie, chaque station de base BS1, BS2 est équipée d’une antenne comprenant une matrice d’éléments rayonnants 200er comme illustré par les figures 2a et 2b . Cette antenne est configurée pour émettre des faisceaux radiofréquences selon au moins une première direction de propagation et une deuxième direction de propagation qui vont être définies ci-après en fonction d’angles θ i,escan et φ i,etilt . A cet égard, on considère, pour chaque antenne (i.e. chaque matrice d’éléments rayonnants de chaque station de base), le modèle de diagramme de rayonnement d’antenne A A,Beami (θ, φ) caractérisant la puissance d’un i-ème faisceau radiofréquences (« beam » en anglais) émis en fonction d’une direction d’observation (θ, φ) par exemple spécifié dans le document 3GPP TR 37.842 V.13.2.0. Il est précisé toutefois ici que les notations utilisées, i.e. θ et φ, dans la présente demande pour les angles représentatifs de la longitude et de la latitude de l’équipement utilisateur UE dans le cas MIMO s’appuient sur une convention différente de celle utilisée dans le document 3GPP TR 37.842 V.13.2.0 : l’angle θ dans la présente demande s’apparente à l’angle φ décrit dans le document 3GPP désignant la longitude, et l’angle φ dans la présente demande s’apparente à l’angle décrit dans le document 3GPP où θ désigne la colatitude. Il est rappelé que la colatitude est égale à la latitude plus π/2. Ce modèle de diagramme de rayonnement s’exprime comme suit : (Eq. 21) avec : où θ i,escan et φ i,etilt représentent la longitude et la latitude définissant la direction de propagation voulue dans le repère O’xyz de l’antenne considérée (cf. définition donnée précédemment en référence à la figure 2a) pour le i-ème faisceau radiofréquences respectivement dans les plans de définition des angles de latitude θ et de longitude φ de l’équipement utilisateur UE dans ce repère (qui correspondent aux angles θ et φ du cas SISO, et définissent la direction d’observation de l’utilisateur vue par l’antenne, notée (θ 1 , φ 1 ) pour l’antenne de la première station de base BS1, (θ 2 , φ 2 ) pour l’antenne de la deuxième station de base BS2 et (θ, φ) quand on désigne indifféremment l’antenne de la première ou de la deuxième station de base), et A E (θ,φ) représente le diagramme de rayonnement, dans le repère 0’ xyz , de chacun des éléments rayonnants 200er de la matrice 200 (supposés identiques ici). De la sorte, une puissance P200i(UE) mesurée par l’équipement utilisateur UE et correspondant au i-ème faisceau radiofréquences émis par la matrice 200 s’exprime de manière générale selon l’expression : P200i(UE) = K.P0.r -eta . A A,Beami (θ, φ).X BS (UE).Y BS (UE) (21bis) Une telle expression présente la même structure que l’équation (Eq. 1) et les raisonnements détaillés ci-dessus s’appliquent de la même façon. Ainsi, dans le cas présent également, la connaissance du modèle de diagramme de rayonnement de la matrice 200 d’éléments rayonnants 200er, caractérisant la puissance, en fonction d’une direction d’observation (celle de l’équipement utilisateur UE par l’antenne considérée), du faisceau radiofréquences considéré (e.g. l’expression de A A,Beami (θ, φ) pour un i-ème faisceau radiofréquences dont la direction de propagation voulue est θ i,escan et φ i,etilt et un j-ème faisceau radiofréquences dont la direction de propagation voulue est θ i,escan et φ j,estilt ) , permet de déterminer en E320 au moins une équation (similaire aux équations (Eq. 3lin) et équations (Eq. 3dB) détaillées ci-dessus dans le cas SISO) dont les membres sont fonction, d’une part, de l’information de puissance relative considérée calculée par mise en œuvre des étapes E300 et E310 et, d’autre part, d’une valeur attendue de l’information de puissance relative fonction, pour chaque faisceau radiofréquences associé à l’information de puissance relative considérée, du modèle de diagramme de rayonnement caractérisant la puissance de ce faisceau radiofréquences. A titre d’exemple, on considère ici pour chaque station de base BS1 et BS2 une matrice d’éléments rayonnants de dimensions Nh = 2 ; Nv = 2 et on considère des faisceaux radiofréquences émis par ses éléments selon des directions de propagation différentes. Par exemple, la première station de base BS1 émet, via sa matrice d’éléments rayonnants, des faisceaux f1 et f2 caractérisés respectivement par leur direction de propagation (θ 1,iescan φ 1,ietilt ) et (θ 2,esan et φ 2,etilt ). On peut ainsi écrire : P200_1(UE) = K.P0.r -eta . A A,Beam1 (θ, φ).X BS (UE).Y BS (UE) pour le premier faisceau, et P200_2(UE) = K.P0.r -eta . A A,Beam2 (θ, φ).X BS (UE).Y BS (UE) pour le deuxième faisceau. Ainsi, on a, pour la première station de base BS1, l’information de puissance relative M12 suivante : M12 = P200_1(UE)/P200_2(UE) = A A,Beam1 (θ, φ)/ A A,Beam2 (θ, φ) (Eq. 3MIMO) ou, alternativement, en décibels : M12 dB = A A,Beam1 (θ, φ) dB - A A,Beam2 (θ, φ) dB ; avec : A A,Beam1 (θ, φ). = A E (θ, φ) + 10log 10 ( (Module1)²); A A,Beam2 (θ, φ). = A E (θ, φ) + 10log 10 ( (Module2)²) où (Module1)² s’exprime à l’aide de l’équation 21 de la façon suivante. On considère : dh/λ=a et dv/λ=b où λ désigne la longueur d’onde 1/2 = c ; a, b et c désignant des nombres réels. Par ailleurs, on pose : a1 = ; b1 = b . et : v11 = 1 v12 = v21 = v22 = . w111 = c w112 = c. = c w121 = c. = c w122 = c. . donc w122 = c. . Il en résulte l’expression suivante de Module1 2 : Module1² = | w111.v11 +w112.v12 +w121.v21 +w122.v22 |² Module1² = | . + c. . + . . . |² = c²|1 + + + |² On pose Module1² = F1(θ,φ, θ 1,escan ,φ 1,etilt ) donc on obtient : 10log 10 ( Module1²) = F1(θ,φ, θ 1,escan ,φ 1,etilt ) dB De même pour le deuxième faisceau, on obtient 10log 10 ( Module2²) = F2(θ,φ, θ 2,escan ,φ 2,etilt ) dB Dans l’exemple envisagé ici, on considère en outre deux autres faisceaux f3, f4 émis par la première station de base BS1 selon des directions de propagation différentes. On peut alors obtenir, en appliquant ce qui précède aux faisceaux f3 et f4, deux expressions analogues à (Eq. 3dB) d’informations de puissance relative : M12 dB = F1(θ,φ, θ 1,escan ,φ 1,etilt ) dB – F2(θ,φ, θ 2,escan ,φ 2,etilt ) dB M34 dB = F3(θ,φ, θ 3,escan ,φ 3,etilt ) dB – F4(θ,φ, θ 4,escan ,φ 4,etilt ) dB On obtient donc 2 équations à 2 inconnues. On note que les inconnues θ et φ sont des fonctions de sinus et de cosinus. Il y a donc potentiellement plusieurs solutions. En considérant quatre faisceaux f1, f2, f3, f4 (ou plus) au lieu de deux comme décrit dans le cas SISO, il est possible de réduire le nombre de solutions potentielles. Toutefois l’invention s’applique également en considérant uniquement deux faisceaux pour chaque station de base. On comprend à la lumière de cet exemple que la complexité de l’équation (Eq. 21) fait qu’une résolution analytique est difficilement envisageable dans le cas général afin de déterminer la direction (θ, φ) de l’équipement utilisateur UE dans le repère Oxyz de la première station de base BS1 (notée (θ1,φ1)). Ainsi, dans certaines variantes de réalisation, la résolution d’une telle équation comprend la mise en œuvre d’une méthode de résolution numérique. Par exemple, en E721, on considère I faisceaux, avec I entier supérieur ou égal à deux, émis par la première station de base BS1, un i-ème faisceau étant caractérisé par un couple (θ i,etilt , φ i,escan) . Considérons par exemple que l’antenne de la première station de base BS1 possède une matrice de 4 éléments rayonnants (Nh=4, Nv=4) et le cas de I=4 faisceaux émis f1, f2, f3, f4. Pour chaque faisceau, la puissance reçue par l’équipement utilisateur UE répond aux équations 21 et 21bis précédentes. Le rapport ou la différence des puissances, reçues et mesurées par l’équipement utilisateur UE, correspondant aux faisceaux f1 et f2, permet de calculer une valeur d’information de puissance relative M12, tandis que le rapport ou la différence des puissances, reçues par l’équipement utilisateur UE, correspondant aux faisceaux f3 et f4, donne une valeur d’information de puissance relative telle que le rapport M34. Ces deux valeurs satisfont chacune l‘équation 21 (dans laquelle le facteur Ae(θ i , φ i ) a été éliminé). On obtient donc deux équations fonction chacune de (θ i,escan , φ i,etilt ) dont les inconnues sont θ 1 et φ 1 , qui correspondent aux directions de l’équipement utilisateur UE dans le repère Oxyz de la première station de base BS1. La résolution de ces deux équations permet dans cet exemple d’obtenir deux couples de valeurs possibles pour (θ 1 , φ 1 ). On note toutefois que si la matrice 200er comprend un plus grand nombre d’éléments rayonnants, les équations 21 sont plus complexes, ce qui fait que le nombre de couple de valeurs possibles pour (θ 1 , φ 1 ) est supérieur à deux (par rapport au cas SISO décrit précédemment). On note que la détermination d’une position de l’équipement utilisateur UE selon ce mode de réalisation de l’invention est applicable à un nombre plus élevé de couples de valeurs possibles, de préférence inférieur à 10. Toutefois, pour réduire encore davantage le nombre de couples de valeurs possibles pour (θ 1 , φ 1 ), on peut utiliser d’autres faisceaux émis par la matrice d’éléments rayonnants de la première station de base, par exemple f5, f6, pour calculer une information de puissance relative supplémentaire, telle que le rapport M56, de la même façon qu’on a calculé les informations de puissance relative M12 et M34. On établit ainsi une 3ème équation fonction de (θ 5,escan , φ 5,etiltet (θ 6,escan , φ 6,etilt) dont les inconnues sont θ et φ, ce qui permet caractériser de manière plus précise le couple (θ, φ), l’ensemble des valeurs possibles prises par ce couple (θ , φ) dans le système d’équations élargi ainsi obtenu étant plus restreint. Si nécessaire, on répète l’opération avec encore d’autres faisceaux jusqu’à obtenir un système d’équations dont la résolution conduit à un nombre de valeurs possibles du couple (θ 1 , φ 1 ) le plus petit possible. Selon une variante de l’étape E721, la résolution d’un tel système d’équations comprend l’exploration de l’espace des solutions (θ 1 , φ 1 ) afin de déterminer la direction (θ 1 , φ 1 ) de l’équipement utilisateur UE qui est solution ou présente des valeurs proches des solutions de toutes les équations établies avec tous les couples de faisceaux considérés. Par exemple, la résolution met en œuvre, pour une information de puissance relative donnée calculée lors de la mise en œuvre de l’étape E310 : - une obtention de la valeur attendue ou théorique de cette information de puissance relative donnée pour un jeu de différentes directions d’observation (θ1, φ1), telle que dérivée du modèle de rayonnement considéré ci-dessus. Par exemple, pour chaque paire de faisceaux émis par la matrice d’éléments rayonnants de la première station de base, on fait varier (θ1, φ1) degré par degré et on associe à chaque nouvelle valeur du couple, une valeur d’information de puissance relative attendue, c’est-à-dire telle que donnée par l’équation 21. Un jeu de valeurs attendues correspondant chacune à une direction d’observation (θ1, φ1) est ainsi obtenu; et - une comparaison entre, d’une part, l’information de puissance relative calculée par l’étape E310, c’est-à-dire issues des mesures de puissance et, d’autre part, chaque valeur attendue du jeu de valeurs attendues associé au couple de valeurs (θ1, φ1). On sélectionne les valeurs de (θ1, φ1) dont la valeur d’information de puissance relative attendue associée est la plus proche (e.g. au sens d’une norme donnée du type valeur absolue) de la valeur de l’information de puissance relative calculée, issue des mesures de puissances (E300, E310), avec une imprécision de l’ordre de l’erreur de mesure. A l’issue de cette résolution numérique ou exploratoire, on obtient donc un nombre restreint de valeurs possibles de (θ1, φ1), idéalement compris entre 2 et 9. On répète ces opérations en E722 pour la deuxième station de base BS2 et on obtient de façon similaire un nombre restreint de valeurs possibles pour un couple (θ2, φ2) définissant la position de l’équipement utilisateur UE dans le repère O 2 xyz de la deuxième station de base BS2. On note que ce nombre peut différer d’une station de base à l’autre. On note que le fait d’utiliser deux stations de base permet a priori de caractériser la direction de l’équipement utilisateur UE à partir d’un nombre plus restreint de faisceaux émis par chacune d’elles. Un avantage est que le système d’équations établi comprend un nombre d’équations plus réduit, et peut être résolu avec une méthode moins complexe. Ainsi, dans l’exemple illustratif envisagé précédemment, on a considéré 4 faisceaux par station de base, mais on aurait pu se limiter à 2 faisceaux par station de base. D’un point de vue géométrique, les deux valeurs possibles de l’angle de longitude θ1 dans le repère de la première station de base BS1 et de l’angle de longitude θ2 dans le repère de la deuxième station de base BS2 définissent autant de plans verticaux qui se coupent en quatre droites verticales, comme précédemment décrit dans le cas SISO. Chacune de ces droites verticales coupe le plan horizontal correspondant à l’altitude donnée de l’équipement utilisateur UE en un point. A ce stade, l’ensemble de ces quatre points constitue donc l’ensemble des solutions, c’est-à-dire des positions possibles de l’équipement utilisateur UE. Pour déterminer la position de l’équipement utilisateur UE dans le repère Oxyz de la première station de base BS1, on cherche en E723 à déterminer son angle de latitude φ 1 . Pour ce faire, on utilise l’équation (20), déjà décrite pour le cas SISO et que l’on rappelle ci-dessous : φ1’ = Arctan ((H / ( y0 - x0.tg(θ2) )).cos(θ1) .[ tg(θ1) - tg(θ2)] ) (20) où (x0,y0) correspond à la position de la deuxième station de base BS2 dans le repère de la première. H correspond ici à la -hauteur de la matrice d’antennes de la première station de base BS1 dans le repère Oxyz, comme illustré par la figure 2b , θ1 est l’angle de longitude de l’équipement utilisateur UE dans le repère Oxyz de la première station de base BS1 et θ2 l’angle de longitude de l’équipement utilisateur UE dans le repère de la deuxième station de base BS2. Plus précisément, on calcule une valeur de φ1’ qui satisfait cette équation (20) pour chacun des couples de valeurs possibles d’angles de latitude (θ1, θ2) obtenus au cours des étapes précédentes E721, E722 et on compare en E724 les valeurs de φ1’ obtenues aux valeurs de φ1 associées à θ1 dans les couples (θ1, φ1) issus de l’étape E721. Si une seule valeur de φ1’correspond, à l’imprécision de mesure de puissance près, à la valeur de φ1 d’un seul des couples (θ1, φ1), on peut en déduire que les angles de longitude et de latitude de l’équipement utilisateur UE sont déterminés par ce couple (θ1, φ1). Si plus d’une valeur de φ1’correspondait, alors il conviendrait d’utiliser des mesures de puissances issus de faisceaux supplémentaires reçus des deux stations de base comme déjà évoqué ci-dessus. Enfin, en E725, on déduit les coordonnées x1, y1 de l’équipement utilisateur UE de l’équation 18 précédemment décrite : x1= (y0 - x0.tg (θ2) ) /( tg(θ1) - tg(θ2) ) (Eq. 18) y1 = x1.tg (θ1) z1 vaut zéro si l’équipement utilisateur UE se trouve au niveau du sol (sur lequel se trouve l’origine O du repère de la première station de base). Les différents cas considérés ci-dessus dans le cas SISO, notamment pour le cas d’un équipement utilisateur UE situé sur un bateau en mer ou de changement de référence pour les repères associés aux stations de base, s’appliquent de la même façon au cas MIMO présenté ici. Les formules obtenues ci-dessus sont liées à la convention adoptée, notamment en ce qui concerne la définition du repère 0 xyz et des angles associés (i.e. longitude et latitude), mais d’autres formules équivalentes peuvent être envisagées si une autre convention est adoptée. On présente désormais, en relation avec la figure 9 un exemple de structure matérielle du dispositif 100 permettant de mettre en œuvre les étapes du procédé de géolocalisation selon l’invention. Le dispositif 100 comprend une mémoire vive 103 (par exemple une mémoire RAM), une unité de traitement 102 équipée par exemple d'un processeur, et pilotée par un programme d'ordinateur stocké dans une mémoire morte 101 (par exemple une mémoire ROM ou un disque dur). A l'initialisation, les instructions de code du programme d’ordinateur sont par exemple chargées dans la mémoire vive 103 avant d'être exécutées par le processeur de l’unité de traitement 102. Cette fig. 8 illustre seulement une manière particulière, parmi plusieurs possibles, de réaliser le dispositif 100 afin qu’il effectue les étapes du procédé de géolocalisation (selon l’un quelconque des modes de réalisation et/ou variantes décrit(e)s ci-dessus en relation avec les figures 4, 6 et 7 ). En effet, ces étapes peuvent être réalisées indifféremment sur une machine de calcul reprogrammable (un ordinateur PC, un processeur DSP ou un microcontrôleur) exécutant un programme comprenant une séquence d’instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre module matériel). Dans le cas où le dispositif 100 est réalisé avec une machine de calcul reprogrammable, le programme correspondant (c'est-à-dire la séquence d’instructions) pourra être stocké dans un médium de stockage amovible (tel que par exemple un CD-ROM, un DVD-ROM, une clé USB) ou non, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. Dans certains modes de réalisation, le dispositif 100 est inclus dans l’équipement utilisateur UE. Dans certains modes de réalisation, le dispositif 100 est inclus dans un dispositif du réseau de radiocommunications, e.g. dans un nœud du réseau de communications ou dans une des stations de base BS1 ou BS2, comme illustré par la figure 3. Procédé de géolocalisation d’un équipement utilisateur (UE) situé à une altitude donnée, recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne (A11, A21, 200er) d’une première station de base (BS1) configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne (A12, A22, 200er) d’une deuxième station de base (BS2) configurée pour émettre selon au moins une troisième et une quatrième direction de propagation, ledit procédé mettant en œuvre : - une obtention (E300) d’une puissance mesurée par ledit équipement terminal pour des faisceaux radiofréquences émis par la première station de base au moins selon la première et la deuxième direction de propagation et des faisceaux radiofréquences émis par la deuxième station de base selon au moins la troisième et la quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; - un calcul (E310) d’au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative (M11) pour la première station de base et deuxième information de puissance relative (M12) pour la deuxième station de base, la première puissance étant associée pour première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de propagation ; et - une détermination (E320) d’une position dudit équipement utilisateur à l’altitude donnée à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative (M11, M12) et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction de la direction d’observation de l’équipement utilisateur par ladite antenne. Procédé selon la revendication 1, caractérisé en ce que ladite détermination (E330) met en œuvre, pour au moins une information de puissance relative donnée parmi la ou les informations de puissances relatives, la résolution d’une équation dont les membres sont fonction, d’une part, de ladite information de puissance relative calculée et, d’autre part, d’une valeur attendue de ladite information de puissance relative donnée, fonction, pour chaque faisceau radiofréquences associé à ladite information de puissance relative donnée, du modèle de diagramme de rayonnement caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative donnée, en fonction de la direction d’observation de l’équipement utilisateur (UE) par ladite antenne. Procédé selon la revendication 2, caractérisé en ce que, la première station de base comprenant une première antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la première direction de propagation et une deuxième antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la deuxième direction de propagation et la deuxième station de base comprenant une première antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la troisième direction de propagation et une deuxième antenne sectorielle configurée pour émettre des faisceaux radio fréquences selon la quatrième direction de propagation, ladite obtention (E300) comprend l’obtention d’une première et d’une deuxième puissances mesurées correspondant respectivement à un premier faisceau radiofréquences émis par la première, respectivement la deuxième antenne (A11, A12) sectorielle de ladite première, respectivement deuxième, station de base et à un deuxième faisceau radiofréquences émis par la première, respectivement la deuxième antenne sectorielle (A21, A22) de ladite première , respectivement deuxième, station de base, dans lequel ledit calcul (E310) délivre une première, respectivement une deuxième, information de puissance (M11, M12) relative associée auxdits premier et deuxième faisceaux radiofréquences pour chaque station de base (BS1, BS2), et dans lequel ladite détermination (E320) comprend la détermination (E321 ; E322) d’un angle (θ11, θ12) représentatif d’une longitude dudit équipement utilisateur (UE) dans un repère de la première, respectivement de la deuxième, station de base (BS1, BS2) à partir, d’une part, de ladite première, respectivement de la deuxième, information de puissance relative (M11, M12) et, d’autre part, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction desdits angles de longitude et de latitude de l’équipement utilisateur, desdits premier et deuxième faisceaux radiofréquences. Procédé selon la revendication 3, caractérisé en ce que ladite résolution comprend, pour chaque station de base, la mise en œuvre de l’équation : θ1i = - (M1i dB /24 ) θ 2 3dB / θb i + θb i /2, où : i est un entier qui vaut 1 lorsqu’il désigne la première station de base et 2 lorsqu’il désigne la deuxième station de base ; θ1i représente ledit angle représentatif d’une longitude dudit équipement utilisateur dans lerepère de la i-ème station de base, M1i dB représente ladite i-ème information de puissance relative exprimée en décibels, θ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement de chacune desdites première et deuxième antennes dans un plan de définition dudit angle θ1i, et θb i représente l’angle de dépointage entre lesdites première et deuxième antennes de la i-ème station de base dans ledit plan de définition dudit angle θ11. Procédé selon l’une quelconque des revendications 3 et 4, caractérisé en ce que lesdits premier et deuxième faisceaux radiofréquences sont rayonnés par lesdites première et deuxième antennes de la première station de base (BS1) selon un premier angle d’inclinaison (φt1), ladite obtention (E300) comprend l’obtention (E301) d’une troisième puissance mesurée correspondant à un troisième faisceau radiofréquences émis par ladite première antenne (A11) de la première station de base (BS1) selon un troisième angle d’inclinaison (φt3), distinct du premier angle d’inclinaison (φt1), ledit calcul (E310) délivre une troisième information de puissance relative (M’11), représentative d’un rapport ou d’une différence entre la troisième puissance mesurée et la première puissance du premier jeu de puissances mesurées ; et ladite détermination (E320) comprend une première détermination (E323) d’un angle (φ11) représentatif d’une latitude dudit équipement utilisateur dans le repère de la première station de base à partir, d’une part, de ladite troisième information de puissance relative et, d’autre part, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance, en fonction desdits angles de longitude et de latitude de l’équipement utilisateur , desdits premier, deuxième, et troisième faisceaux radiofréquences de la première station de base (BS1). Procédé selon la revendication 5, caractérisé en ce que ladite détermination (E320) comprend la mise en œuvre de l’équation : φ11 = M12 dB /24 φ 2 3dB / (φt1– φt3) + (φt1+ φt3)/2, où : φ11 représente ledit angle représentatif d’une latitude dudit équipement utilisateur dans le repère de la première station de base, M12 dB représente ladite deuxième information de puissance relative exprimée en décibels, φ 3dB représente l’angle d’ouverture à trois décibels du diagramme de rayonnement de chacune desdites première et deuxième antennes dans un plan de définition dudit angle φ11, φt1 représente ledit premier angle d’inclinaison, et φt3 représente ledit troisième angle d’inclinaison. Procédé selon la revendication 6, caractérisé en ce que ladite détermination (E320) comprend une deuxième détermination (E324) d’une valeur de l’angle de latitude (φ11’) de l’équipement utilisateur (UE) dans le repère ( Oxyz ) de la première station de base (BS1), en fonction d’une hauteur (H) de la première antenne de la première station de base (BS1) et des angles de longitude (θ11, θ21) de l’équipement utilisateur dans les repères respectifs de la première et de la deuxième station de base et une sélection (E325) de la valeur des angles de longitude (θ11, θ12) de l’équipement utilisateur (UE) minimisant une erreur entre la première et la deuxième détermination de l’angle de latitude (φ11, φ11’). Procédé selon l’une quelconque des revendications 6 et 7, caractérisé en ce que ladite détermination (E320) comprend le calcul (E326) des coordonnées (x1,y1,z1) de l’équipement utilisateur (UE) dans le repère ( Oxyz ) de la première station de base à partir des angles de longitude (θ11, θ12) déterminés. Procédé selon la revendication 1, caractérisé en ce que, la première station de base comprenant une antenne comprenant une matrice d’éléments rayonnants (200er), configurée pour émettre des faisceaux radiofréquences selon au moins la première direction de propagation et la deuxième direction de propagation et la deuxième station de base comprenant une antenne comprenant une matrice d’éléments rayonnants (200er) configurée pour émettre des faisceaux radiofréquences selon au moins la troisième direction de propagation et la quatrième direction de propagation, chaque faisceau radiofréquences de ladite pluralité de faisceaux radiofréquences est rayonné par la matrice d’éléments rayonnants de ladite première, respectivement deuxième, station de base. Procédé selon la revendication 9, caractérisé en ce qu ’il dépend de la revendication 2, dans lequel ladite résolution d’équation comprend la mise en œuvre d’une méthode de résolution numérique de ladite au moins une équation. Procédé selon la revendication 9, caractérisé en ce qu ’il dépend de la revendication 2, dans lequel ladite résolution met en œuvre pour ladite au moins une information de puissance relative donnée : une obtention de ladite valeur attendue de ladite information de puissance relative donnée pour un jeu de directions d’observation de l’équipement utilisateur délivrant un jeu de valeurs attendues correspondant chacune à une direction d’observation de l’équipement utilisateur parmi une pluralité de directions d’observations; une comparaison entre, d’une part, ladite information de puissance relative calculée et, d’autre part, chaque valeur attendue dudit jeu de valeurs attendues délivrant une direction d’observation pour laquelle la valeur attendue de ladite puissance relative donnée est la plus proche de l’information de puissance relative calculée ; et ladite position dudit équipement utilisateur étant fonction de la direction d’observation délivrée. Procédé selon l’une quelconque des revendications 9 à 11, caractérisé en ce que , pour la première, respectivement la deuxième station de base, ledit modèle de diagramme de rayonnement d’antenne caractérisant la puissance d’un i-ème faisceau radiofréquences en fonction d’une direction d’observation de l’équipement utilisateur (θ, φ) s’exprime en fonction de : avec : où θ est un angle représentatif d’une longitude de l’équipement utilisateur dans un repère de la première, respectivement deuxième station de base, φ représente un angle représentatif d’une latitude de l’équipement utilisateur dans le repère de la première, respectivement deuxième station de base, (θ i,escan, φ i,etilt) représente un couple d’angles caractéristique de la direction de propagation principale voulue pour ledit i-ème faisceau radiofréquences dans un repère ( O’xyz )de l’antenne, N H représente le nombre d’éléments rayonnants de ladite matrice d’éléments rayonnants dans une direction horizontale, N V représente le nombre d’éléments rayonnants de ladite matrice d’éléments rayonnants dans une direction verticale, d V représente l’espacement vertical entre deux éléments rayonnants, d H représente l’espacement horizontal entre deux éléments rayonnants, λ représente la longueur d’onde du faisceau radiofréquences, et A E (θ,φ) représente la diagramme de rayonnement de chaque élément rayonnant de ladite matrice d’éléments rayonnants ; en ce que ladite détermination (E320) comprend une première détermination (E721) d’au moins deux couples de valeurs de l’angle de longitude et de l’angle de latitude (θ1, φ1) de l’équipement utilisateur (UE) dans le repère ( Oxyz ) de la première station de base (BS1) à partir dudit modèle de rayonnement, une deuxième détermination (E722) d’au moins deux couples de valeurs de l’angle de longitude et de l’angle de latitude (θ2, φ2) de l’équipement utilisateur (UE) dans le repère de la deuxième station de base (BS2) à partir dudit modèle de rayonnement, une troisième détermination (E723) d’une deuxième valeur de l’angle de latitude (φ1’) de l’équipement utilisateur (UE) dans le repère de la première station de base (BS1) ), en fonction d’une hauteur (H) de la première station de base (BS1) et des angles de longitude (θ1, θ2) de l’équipement utilisateur dans les repères respectifs de la première et de la deuxième station de base et une sélection (E324) des valeurs des angles de longitude (θ1, θ2) de l’équipement utilisateur (UE) minimisant une erreur entre la valeur de l’angle de latitude (φ1) associé à l’angle de latitude (θ 1 ) issue de la première détermination et la deuxième valeur de l’angle de latitude (φ1’) issue de la troisième détermination. Produit programme d’ordinateur comprenant des instructions de code de programme pour la mise en œuvre du procédé selon l’une quelconque des revendications 1 à 12, lorsque ledit programme est exécuté sur un ordinateur. Dispositif (100) de géolocalisation d’un équipement utilisateur recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne (A11, A21, 200er) d’une première station de base (BS1) configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne (A12, A22, 200er) d’une deuxième station de base (BS2) configurée pour émettre selon au moins une troisième et une quatrième direction de propagation, caractérisé en ce qu’ il comprend une machine de calcul reprogrammable (102) ou une machine de calcul dédiée configurée pour : - obtenir (OBT. Pij) une puissance mesurée par ledit équipement utilisateur pour des faisceaux radiofréquences émis par la première station de base au moins selon la première et la deuxième direction de propagation et des faisceaux radiofréquences émis par la deuxième station de base au moins selon la troisième et quatrième direction de propagation, délivrant respectivement un premier jeu de puissances mesurées pour la première station de base et un deuxième jeux de puissances mesurées pour la deuxième station de base; - calculer (OBT. IPR) au moins une information par station de base, représentative d’un rapport ou d’une différence entre une première et une deuxième puissances de chaque jeu de puissances mesurées, dite première information de puissance relative (M11) pour la première station de base et deuxième information de puissance relative (M12) pour la deuxième station de base, la première puissance étant associée pour la première station de base à un faisceau radiofréquences émis selon la première direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la deuxième direction de propagation, la première puissance étant associée pour la deuxième station de base à un faisceau radiofréquences émis selon la troisième direction de propagation et la deuxième puissance étant associée à un faisceau émis selon la quatrième direction de propagation; et - déterminer (DET. POS) une position dudit équipement utilisateur à partir, d’une part, de ladite au moins une première et deuxième informations de puissance relative (M11, M12) et, d’autre part, pour chaque information de puissance relative et pour chaque faisceau radiofréquences associé à ladite information de puissance relative, d’un modèle de diagramme de rayonnement d’antenne caractérisant la puissance dudit faisceau radiofréquences associé à ladite information de puissance relative, en fonction d’une direction d’observation de l’équipement utilisateur par ladite antenne. Equipement utilisateur (UE) comprenant un dispositif (100) de géolocalisation selon la revendication 14. Station de base (BS1, BS2) comprenant un dispositif (100) de géolocalisation selon la revendication 14. Système (10) de géolocalisation d’un équipement utilisateur recevant une pluralité de faisceaux radiofréquences émis respectivement par au moins une antenne (A11, A21, 200er) d’une première station de base (BS1) configurée pour émettre selon au moins une première et une deuxième direction de propagation, et par au moins une antenne (A12, A22, 200er) d’une deuxième station de base (BS2) configurée pour émettre selon au moins une troisième et une quatrième direction de propagation, caractérisé en ce qu’il comprend au moins la première et la deuxième stations de base et un dispositif de géolocalisation selon la revendication 14.