FR 3058533 A1 20180511 FR 1750507 A 20170123 1 Procédés de détermination de la position d'une balise de géolocalisation, produit programme d'ordinateur, balise, station de base et dispositif de calcul de position correspondants 1 DOMAINE TECHNIQUE La présente description se rapporte au domaine des balises de géolocalisation (également appelées dispositifs de géolocalisation ou terminaux de géolocalisation), et plus particulièrement aux méthodes adressant la réduction de consommation de telles balises. 2 ARRIERE-PLAN TECHNOLOGIQUE Présentés comme la « troisième révolution de l'Internet », les objets connectés sont en train de s'imposer dans tous les domaines de la vie quotidienne et de l'entreprise. La plupart de ces objets sont destinés à produire des données grâce à leurs capteurs intégrés afin de fournir des services à valeur ajoutée pour leur propriétaire. De par les applications visées, la plupart de ces objets connectés sont nomades. En particulier, ils doivent pouvoir transmettre les données produites, régulièrement ou à la demande, à un utilisateur déporté. Pour ce faire, la transmission radio longue portée du type radio mobile cellulaire (2G/3G/4G...) a été une technologie de choix. Cette technologie permettait en effet de bénéficier d'une bonne couverture réseau dans la plupart des pays. Cependant, l'aspect nomade de ces objets s'accompagne souvent d'un besoin en autonomie d'énergie. Or, même basés sur une des technologies radio mobile cellulaire les plus économes en énergie, les objets connectés actuels continuent de présenter une consommation rédhibitoire pour permettre un déploiement à grande échelle à un coût raisonnable. Face à la problématique de la consommation du lien radio pour de telles applications nomades, de nouvelles technologies radio basse consommation et bas débit dédiées spécifiquement aux réseaux « Internet des Objets », c'est-à-dire des technologies radio pour des réseaux dits LPWAN (pour « Low-Power Wide-Area Networks » en anglais), sont en train de voir le jour. A titre d'exemple, nous pouvons citer la technologie de la société Sigfox (www.sigfox.com), ou bien la technologie LoRa 3058533 2 (www.lora-alliance.org), ou encore la technologie de la société Qowisio (www.qowisio.com). Dans la suite de la description, on appelle « réseau LPWAN » n'importe quel type de réseau du type basse consommation et bas débit dédié aux objets connectés (réseau 5 Sigfox, réseau loRa, etc.). Par ailleurs, on appelle « lien radio LPWAN » un lien radio vers un réseau LPWAN. Les technologies radio LPWAN ont la particularité d'être très peu consommatrices d'énergie tout en fournissant une couverture nationale voire internationale grâce à une portée radio significative permettant un maillage régulier du 10 territoire à couvrir. Pour arriver à ce résultat, ces technologies proposent des capacités très faibles dans les sens montant et descendant, en adéquation avec les besoins des objets connectés. A titre d'exemple, le réseau Sigfox offre la possibilité : - de transmettre au maximum 144 trames radio par jour dans le sens montant, 15 i.e. de l'objet connecté vers le réseau Sigfox, chacune de ces trames radio comprenant une partie utile de 12 octets destinée aux données à transmettre ; et - de transmettre de 1 à 4 trames radio chacune dans le sens descendant, i.e. du réseau Sigfox vers l'objet connecté, chacune de ces trames radio comprenant 20 une partie utile de 8 octets. Pour la technologie LoRa, la longueur de la partie utile des trames radio est variable dans le sens montant, et peut théoriquement aller jusqu'à 59 octets. Cependant, la consommation en courant pour l'émission de la trame est quasiment proportionnelle à la longueur de la trame. De ce fait, la longueur des trames radio est 25 alors limitée en pratique à 12 ou 16 octets de données utiles. De même, le débit dans le sens descendant est fortement réduit. Par exemple en Europe, où ce type de standard suit la réglementation de la bande de fréquence « Industriel, Scientifique et Médical », dite ISM, les stations de base ne peuvent émettre plus de 1% du temps. La capacité du réseau LoRa dans le sens descendant est alors environ cinquante fois plus faible que 30 dans le sens montant. 3058533 3 Il apparaît que l'utilisation d'une telle technologie devrait permettre de diminuer drastiquement la consommation de la partie radio longue portée d'un objet connecté tout en répondant aux besoins de transmission faible débit des données générée par la plupart de ces objets. 5 Parmi de tels objets connectés nomades, les balises de géolocalisation occupent une position de premier plan. En effet, les besoins en géolocalisation restent nombreux, aussi bien dans le domaine industriel (e.g. suivi de véhicules, de marchandises, de personnes) que dans le domaine grand public (e.g. suivi des animaux domestiques). Classiquement, de telles balises déterminent leur position en embarquant un 10 dispositif de géolocalisation basé sur un système de positionnement par satellites (ou plus généralement « véhicules spatiaux »), du type GNSS (pour « Global Navigation Satellite System » en anglais). Il peut s'agir, par exemple, d'un récepteur GPS, Glonass, Beidou, bientôt Galiléo, etc. Comme pour la plupart des objets connectés nomades, une réduction de 15 consommation significative est attendue pour une telle balise, en adoptant une technologie de transmission des données via un lien radio LPWAN plutôt qu'un lien radio du type mobile cellulaire (lien radio GSM par exemple). Cependant, ce que consomme la balise pour déterminer sa position avec une technologie de géolocalisation du type GPS représente un tiers de ce qu'elle consomme 20 pour la transmission de cette position sur le lien radio GSM. Ainsi, le passage à une technologie radio du type LPWAN pour la transmission longue portée des informations de géolocalisation ne permet qu'un gain de facteur quatre sur la consommation de la balise complète, et donc sur son autonomie au final. En pratique, il apparaît qu'un tel gain reste insuffisant pour pouvoir proposer 25 une balise de géolocalisation ayant une autonomie de plusieurs années à un coût permettant une diffusion à plus grande échelle Il existe donc un besoin pour une solution permettant d'améliorer l'autonomie des balises de géolocalisation, quand bien même ces balises transmettent leurs données via un réseau LPWAN. 30 3 RESUME 3058533 4 Dans un mode de réalisation de l'invention, il est proposé un procédé de détermination de la position d'une balise de géolocalisation connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de 5 géolocalisation provenant d'une constellation de M véhicules spatiaux. Le procédé, mis en oeuvre par la balise de géolocalisation, comprend les étapes suivantes : réception de signaux de géolocalisation provenant de N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' 2 ; détermination d'une mesure de phase de code pour chacun des N' véhicules 10 spatiaux ; obtention de N informations de phases de codes sur la base des N' mesures de phases de codes, avec 2 2, lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la 30 contrainte que le nombre total de bits alloués aux N informations de phases de 3058533 7 codes reste inférieur ou égal à un nombre prédéterminé de bits ; horodatage de ladite trame radio reçue, ou une desdites trames radio reçues, formant une trame radio de référence, avec une information d'horodatage indiquant un instant de réception de ladite trame radio de référence par la 5 station de base ; transmission de la ou les trame(s) reçue(s), y compris ladite trame radio de référence horodatée, vers un dispositif de calcul de position configuré pour déterminer la position de la balise de géolocalisation en fonction des N informations de phases de codes, de l'information d'horodatage et 10 d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux. Ainsi, la station de base horodate, à la date correspondant à l'instant ti3Ts de sa réception, au moins une trame radio de référence parmi la ou les trames radio reçues qui encapsule(n)t une même trame de données. La station de base, transmet ensuite la 15 ou les trames radio reçues, dont cette trame radio de référence horodatée, au dispositif de calcul (qui va calculer la position de la balise). Le dispositif de calcul est alors capable de déterminer, sur la base de cette information d'horodatage et du délai D1 précité, l'instant de réception tR par la balise des signaux de géolocalisation correspondant aux données transmises dans la trame de données correspondante. 20 Dans un autre mode de réalisation de l'invention, il est proposé un procédé de détermination de la position d'une balise de géolocalisation connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux, le procédé étant 25 mis en oeuvre par le dispositif de calcul de position et comprenant les étapes suivantes au moins une fois : réception, via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' 30 véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N 3058533 8 > 2, lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phase de 5 codes reste inférieur ou égal à un nombre prédéterminé de bits ; décodage desdites N informations de phases de codes encodées ; obtention d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux ; et détermination de la position de la balise de géolocalisation en fonction des N 10 informations de phases de codes et des informations d'éphéméride et/ou d'almanach. Ainsi, le dispositif de calcul de position détermine la position de la balise sur la base des mesures de phases de codes transmises par la balise et d'informations d'éphéméride et/ou almanach obtenues préalablement, par exemple auprès d'un 15 serveur dédié. Ainsi, comme déjà mentionné plus haut, la balise ne détermine pas elle- même les éphémérides des véhicules spatiaux dont elle reçoit les signaux. Ceci minimise le temps d'acquisition des signaux reçus en provenance des véhicules spatiaux par la balise, et donc la consommation en énergie de cette dernière. De même, les calculs les plus conséquents pour la détermination de la position de la balise sont effectués par le 20 dispositif de calcul, ce qui permet de réduire la consommation d'énergie de la balise. Selon un mode de réalisation particulier de l'invention, ladite trame radio reçue, ou une desdites trames radio reçues, est une trame radio de référence horodatée et associée avec une information d'horodatage indiquant un instant de réception de ladite trame radio de référence par une station de base dudit réseau de communication, et le 25 dispositif de calcul de position : détermine en fonction de l'information d'horodatage un instant de réception tR, par la balise de géolocalisation, des signaux de géolocalisation provenant des N' véhicules spatiaux ; et détermine la position de la balise de géolocalisation également en fonction dudit 30 instant de réception tR. 3058533 9 Ainsi, dans un mode de réalisation, le dispositif de calcul est apte à déterminer l'instant de réception tR par la balise, des signaux émis par les véhicules spatiaux correspondant aux données qu'il reçoit, données correspondant donc à une mesure de position de la balise à un instant donné. Cette détermination peut se faire sur la base de 5 la connaissance du délai D1, et de l'information d'horodatage de la trame radio de référence. Le dispositif de calcul peut alors utiliser la connaissance de cet instant tR pour déterminer la position de la balise sur la base d'un nombre réduit de mesures de phases de codes. Le nombre N de mesures de phases de codes que la balise doit transmettre est ainsi minimisé, permettant de réduire encore la consommation de la balise et/ou 10 d'allouer plus de bits de codage à chacune des mesures de phases de codes transmises (dans le but d'améliorer la précision obtenue sur la position au final). Selon une variante, le dispositif de calcul de position : connaît une information prédéterminée de délai indiquant un délai (D1) entre un instant de réception, par la balise de géolocalisation, des signaux de géolocalisation 15 provenant des N' véhicules spatiaux et un instant de transmission, par la balise de géolocalisation, de ladite trame radio de référence ; et détermine ledit instant de réception tR en fonction de l'information d'horodatage et ladite information de délai. Ainsi, dans une variante, le dispositif de calcul de position connaît a priori le 20 délai D1, minimisant par là-même la quantité d'information devant transiter sur le lien radio LPWAN. Ceci permet de réduire encore la consommation de la balise et/ou d'allouer plus de bits de codage à chacune des mesures de phases de codes transmises (dans le but d'améliorer la précision obtenue sur la position au final). Selon une autre variante, le dispositif de calcul de position : 25 lit dans ladite au moins une trame de données reçue une information de délai indiquant un délai (D1) entre un instant de réception des signaux de géolocalisation provenant des N' véhicules spatiaux et un instant de transmission de ladite trame radio de référence par la balise de géolocalisation ; et - détermine ledit instant de réception tR en fonction de l'information 30 d'horodatage et ladite information de délai. 3058533 10 Ainsi, dans une autre variante, le dispositif de calcul de position ne connaît pas a priori le délai D1, mais reçoit une information provenant de la balise indiquant ce délai Dl. Ce délai n'étant pas prédéterminé, la balise est alors apte à supporter des délais de transmission variables, ce qui permet par exemple de supporter des protocoles sans 5 accusé de réception systématique nécessitant de transmettre plusieurs fois le même message afin d'assurer la réception d'au moins un message. Selon un aspect particulier de l'invention, le procédé selon l'invention est caractérisé en ce que : N 8. Ainsi, selon un mode de réalisation de l'invention, le procédé est dimensionné 10 pour ne tenir compte que d'un nombre volontairement restreint de mesures de phases de codes, réduisant ainsi la taille des trames de données transmises par la balise au dispositif de calcul via le réseau LPWAN. Par exemple, il est prévu de ne traiter que les signaux des véhicules spatiaux visibles depuis un point donné à la surface de la Terre à un instant donné. Pour mémoire, dans le système GPS, huit véhicules spatiaux sont 15 visibles simultanément en moyenne. Selon une caractéristique particulière de l'invention, ladite au moins une trame de données comprend au maximum kmax octets, avec kmax inférieur ou égal à 20 et préférentiellement égal à 12. Ainsi, selon un mode de réalisation de l'invention, la longueur des trames de 20 données telles qu'encodées par la balise est adaptée à la longueur du champ de données utiles des trames radio de longueur préférentielle de la plupart des réseaux LPWAN, comme par exemple les réseaux Sigfox ou LoRa. Selon un mode de réalisation particulier de l'invention, ladite au moins une trame de données comprend, pour chaque information de phases de codes parmi 25 lesdites N informations de phases de codes : un champ d'en-tête contenant un identifiant pour un véhicule spatial donné parmi lesdits N' véhicules spatiaux ; et un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule spatial donné. 3058533 11 Ainsi, selon un mode de réalisation de l'invention, seuls les champs d'identifiant correspondant à une mesure de phase de code effectivement transmise sont inclus dans la trame de données. L'utilisation globale de la trame de données se trouve optimisée, les bits de données étant redistribués aux mesures de phases de codes lorsque peu 5 d'entre elles, et donc peu de champs d'identifiants, sont transmis au réseau LPWAN. Selon une caractéristique particulière de l'invention, ladite au moins une trame de données comprend : un champ d'en-tête, commun et indiquant un identifiant d'un véhicule spatial donné, parmi lesdits N' véhicules spatiaux, pour chacune des N informations de 10 phases de codes ; et pour chacun des N identifiants de véhicules spatiaux, un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule spatial. Ainsi, selon un mode de réalisation de l'invention, le dispositif de calcul a la 15 connaissance des véhicules spatiaux auxquels correspondent chacune des mesures de phases de codes qu'il a en sa possession. La détermination de la position de la balise s'en trouve simplifiée. Dans une variante, le champ d'en-tête commun est un masque comprenant une séquence de M bits, chaque bit possédant un rang au sein de ladite séquence qui 20 correspond à un identifiant d'un des M véhicules spatiaux de la constellation, et qui est mis à 1 pour indiquer que ladite au moins une trame de données comprend un champ de phase de code pour ledit véhicule spatial et à 0 sinon. L'ordre des champs de phase de code dans ladite au moins une trame de données correspond à un ordre prédéterminé, préférentiellement un ordre croissant ou décroissant, des N identifiants 25 de véhicules spatiaux. Ainsi, dans une variante, un simple masque binaire suffit à identifier les véhicules spatiaux dont les mesures de phases de codes sont transmises, simplifiant par là-même le décodage des identifiants au niveau du dispositif de calcul. Selon une caractéristique particulière, une desdites N' mesures de phases de 30 codes, associée à un des N' véhicules spatiaux, est sélectionnée comme phase de code 3058533 12 de référence. Les N informations de phases de codes sont N'-1 différences entre les N'-1 mesures de phase de code non sélectionnées comme phase de code de référence, et ladite phase de code de référence. Ladite au moins une trame de données comprend lesdites N'-1 différences. La balise et le dispositif de calcul de position partagent une 5 information représentative d'un identifiant dudit véhicule spatial associé à ladite phase de code de référence. Ainsi, la quantité d'information transmise par la balise au dispositif de calcul via le réseau LPWAN se trouve être encore réduite, permettant un gain additionnel en consommation. 10 Dans une autre variante, ladite au moins une trame de données ne comprend pas de champ d'en-tête contenant un ou plusieurs identifiants de véhicule spatial. Pour chacun des N véhicules spatiaux, ladite au moins une trame de données comprend un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule spatial, et l'ordre des champs de phases de codes dans ladite au moins 15 une trame de données correspond à un ordre prédéterminé, préférentiellement un ordre croissant ou décroissant, des identifiants des N véhicules spatiaux. Ainsi, dans cette autre variante, la trame de données ne comprend pas de champ d'entête contenant les identifiants des véhicules spatiaux dont les mesures de phases de codes sont transmises, permettant par là-même de réduire la taille de la 20 trame de données et/ou augmenter le nombre de bits disponibles pour transmettre ces mesures de phases de codes. Par ailleurs, les mesures de phases de codes transmises sont rangées dans la trame de données selon un ordre prédéterminé des identifiants des véhicules spatiaux correspondants. Ainsi, le nombre de combinaisons possibles à tester au niveau du 25 dispositif de calcul pour déterminer à quel véhicule spatial correspond chaque mesure de phase de code se trouve minimisé. Selon un mode de réalisation particulier de l'invention, le nombre N est fixe et prédéterminé. Ceci permet de simplifier le décodage au niveau du dispositif de calcul. 3058533 13 Dans une variante, ladite au moins une trame de données comprend un champ contenant le nombre N. Ainsi, dans une variante, le dispositif de calcul n'a pas à connaître à l'avance le nombre N de mesures de phases de codes transmises dans une trame de données, 5 permettant par là-même d'adapter ce nombre en fonction des conditions opérationnelles de la balise, par exemple en fonction du nombre de véhicules spatiaux dont elle reçoit effectivement les signaux. Ceci permet également de réduire les calculs effectués par le dispositif de calcul (pas nécessaire de tester des combinaisons non existantes). 10 Selon une caractéristique particulière de l'invention, chacun des champs de phase de code est codé de manière adaptative en fonction de la valeur de N, en maximisant le nombre de bits utilisés, sous la contrainte que le nombre de bits alloués aux champs de phase de code reste inférieur ou égal à (8*kmax - H), avec H le nombre de bits alloués au(x) champ(s) d'en-tête, ou à (8*kmax) en l'absence de champ d'en-tête. 15 Ainsi, selon un mode de réalisation de l'invention, le nombre de bits disponibles dans la partie utile de la trame radio est optimisé en fonction du nombre de mesures de phases de codes transmises. L'information transmise est ainsi optimisée pour un nombre total de bits utiles disponibles dans la trame radio LPWAN utilisée. Dans une variante, chaque information de phase de code est composée de E bits 20 pour la partie entière, préférentiellement E=10, et F bits pour la partie fractionnaire, et : si N est inférieur à un seuil S, préférentiellement S=5, la partie entière de l'information de phase de code est codée sur E bits, pour chacun des champs de phase de code; et si N est supérieur ou égal audit seuil S, la partie entière de l'information de 25 phase de code est codée sur E' bits, préférentiellement E'=7, résultant d'une suppression d'au moins un des bits de poids fort parmi lesdits E bits, pour chacun des champs de phase de code. Ainsi, dans une variante, la partie entière des mots binaires codant les mesures de phases de codes est transmise avec troncature lorsque, préférentiellement, plus de 30 quatre d'entre elles sont transmises. En effet, la position de la balise pouvant être 3058533 14 connue approximativement de part son appartenance à une cellule du réseau LPWAN, les bits de poids fort de cette partie entière peuvent être tronqués dans une certaine mesure, tant que l'ambiguïté spatiale correspondant à la dynamique codée par cette partie entière reste inférieure à la taille d'une de ces cellules réseau. La quantité 5 d'information transmise sur le lien radio LPWAN par la balise se trouve ainsi optimisée et la consommation en énergie de la balise également par là-même. Dans une autre variante, chaque information de phase de code est composée de E bits pour la partie entière, préférentiellement E=10, et F bits pour la partie fractionnaire, et : 10 si N est inférieur à un seuil S, préférentiellement S=5, la partie entière de l' information de phase de code est codée sur E bits, pour chacun des champs de phase de code; et si N est supérieur ou égal audit seuil S : pour S-1 champs de phase de code, la partie entière de l'information de 15 phase de code est codée sur E bits ; et pour les champs de phase de code restants, la partie entière de l'information de phase de code est codée sur 0 bit. Ainsi, dans une autre variante, la partie entière des mots binaires codant les mesures de phases de codes est transmise sans troncature pour quatre d'entre elles 20 préférentiellement. Ainsi, la position déterminée par le dispositif de calcul sur la base de ces quatre mesures de phases de codes reste dénuée d'ambiguïté qui serait liée à une dégradation de la dynamique totale du codage de la partie entière. A contrario, une fois cette ambiguïté levée grâce aux quatre premières mesures de phases de codes transmises, il est possible de dégrader la partie entière des mesures de phases de codes 25 suivantes de manière à préserver le nombre de bits utilisés pour coder leur partie fractionnaire. Cette partie fractionnaire est en effet la partie d'intérêt pour améliorer la résolution spatiale et donc la précision sur la position de la balise au final. La quantité d'information transmise sur le lien radio LPWAN par la balise se trouve ainsi optimisée pour une résolution spatiale donnée et la consommation en énergie de la balise 30 également par là-même. 3058533 15 Selon différents modes de réalisation, chacune desdites N informations de phases de codes appartient au groupe comprenant : une desdites N' mesures de phase de code ; et une différence entre une desdites N' mesures de phase de code et une phase de 5 code de référence. Dans un autre mode de réalisation de l'invention, il est proposé un programme d'ordinateur, comprenant des instructions de code de programme pour la mise en oeuvre d'un procédé tel que décrit précédemment, selon l'un quelconque de ses différents modes de réalisation, lorsqu'il est exécuté sur un ordinateur. 10 Dans un autre mode de réalisation de l'invention, il est proposé un médium de stockage lisible par ordinateur et non transitoire, stockant un produit programme d'ordinateur pour la mise en oeuvre d'un procédé tel que décrit précédemment, selon l'un quelconque de ses différents modes de réalisation. Dans un autre mode de réalisation de l'invention, il est proposé une balise de 15 géolocalisation connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux. La balise de géolocalisation comprend : des moyens de réception de signaux de géolocalisation provenant de N' 20 véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' 2; des moyens de détermination d'une mesure de phase de code pour chacun des N' véhicules spatiaux ; des moyens d'obtention de N informations de phases de codes sur la base d'au 25 moins N mesures de phases de codes parmi lesdites N' mesures de phases de codes, N' N 2 ; des moyens d'encodage des N informations de phases de codes dans au moins une trame de données sous la contrainte que le nombre total de bits alloués aux N informations de phases de codes reste inférieur ou égal à un nombre 30 prédéterminé de bits, ledit encodage étant adaptatif en fonction de N ; et 3058533 16 des moyens de transmission de ladite au moins une trame de données, dans au moins une trame radio via le réseau de communication, à destination d'un dispositif de calcul de position configuré pour déterminer la position de la balise de géolocalisation en fonction des N informations de phases de codes et d'informations d'éphéméride 5 et/ou d'almanach décrivant la constellation de véhicules spatiaux. Ainsi, l'invention concerne également une balise de géolocalisation apte à mettre en oeuvre le procédé selon l'invention (selon l'un quelconque des différents modes de réalisation précités). Dans un autre mode de réalisation de l'invention, il est proposé une station de 10 base d'un réseau de communication du type basse consommation et bas débit dédié aux objets connectés auquel est connectée une balise de géolocalisation, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux. La station de base comprend : des moyens de réception, via ledit réseau de communication, d'une ou plusieurs 15 trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N > 2, lesdites N informations de phases de codes étant encodées dans ladite au moins 20 une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phase de codes reste inférieur ou égal à un nombre prédéterminé de bits ; des moyens d'horodatage de ladite trame radio reçue, ou une desdites trames radio reçues, formant une trame radio de référence, avec une information 25 d'horodatage indiquant un instant de réception de ladite trame radio de référence par la station de base ; des moyens de transmission de la ou les trame(s) reçue(s), y compris ladite trame radio de référence horodatée, vers un dispositif de calcul de position configuré pour déterminer la position de la balise de géolocalisation en fonction 30 des N informations de phases de codes, de l'information d'horodatage et 3058533 17 d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux. Ainsi, selon un autre aspect, l'invention concerne également une station de base d'un réseau de communication du type basse consommation et bas débit dédié aux 5 objets connectés apte à mettre en oeuvre le procédé selon l'invention (selon l'un quelconque des différents modes de réalisation précités). Dans un autre mode de réalisation de l'invention, il est proposé un dispositif de calcul de position d'une balise de géolocalisation connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, 10 ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux. Le dispositif de calcul de position comprend : des moyens de réception, via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N 15 informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N > 2, lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la 20 contrainte que le nombre total de bits alloués aux N informations de phase de codes reste inférieur ou égal à un nombre prédéterminé de bits ; des moyens de décodage desdites N informations de phases de codes encodées ; des moyens d'obtention d'informations d'éphéméride et/ou d'almanach 25 décrivant la constellation de véhicules spatiaux ; et des moyens de détermination de la position de la balise de géolocalisation en fonction des N informations de phases de codes et des informations d'éphéméride et/ou d'almanach. 3058533 18 Ainsi, selon encore un autre aspect, l'invention concerne également un système de base de données apte à mettre en oeuvre le procédé selon l'invention (selon l'un quelconque des différents modes de réalisation précités). 4 LISTE DES FIGURES 5 D'autres caractéristiques et avantages de la solution proposée apparaîtront plus clairement à la lecture de la description suivante de modes de réalisation particuliers, donnés à titre de simples exemples illustratifs et non limitatifs, et des dessins annexés, parmi lesquels : la figure la représente un synoptique d'un système selon un mode de réalisation 10 particulier de l'invention, comprenant notamment une balise de géolocalisation en liaison avec un dispositif de calcul de position ; la figure lb illustre des délais de traitement au sein de la balise de géolocalisation ; les figures 2a (étapes côté balise) et 2b (étapes côté réseau) illustrent un 15 procédé de détermination de position selon un mode de réalisation particulier de l'invention (avec des variantes); les figures 3a et 3b illustrent les interactions entre différentes entités du système de la figure la, lors de la mise en oeuvre du procédé des figures 2a et 2b; 20 les figures 4a à 4f illustrent des formats de codage adaptatif des données transmises par la balise selon différents modes de réalisation de l'invention ; et les figures 5a, 5b et 5c présentent des exemples de structures du dispositif encodeur de la balise, de la station de base et du dispositif de calcul de position respectivement, selon un mode de réalisation particulier de l'invention. 25 5 DESCRIPTION DETAILLEE Sur toutes les figures du présent document, les éléments et étapes identiques sont désignés par une même référence. Le principe général de la technique décrite consiste en la transmission, par une balise de géolocalisation, d'informations nécessaires à sa géolocalisation et reçues de 30 véhicules spatiaux d'une constellation d'un système de géolocalisation. Cette 3058533 19 transmission s'effectue via un réseau LPWAN à destination d'un dispositif de calcul de position qui détermine alors la position de la balise à l'aide des informations transmises. Par ailleurs, les données sont encodées par la balise de manière adaptative en fonction du nombre de données à transmettre et sous la contrainte que le nombre total de bits 5 alloués aux données à transmettre reste inférieur ou égal à un nombre prédéterminé de bits. Ainsi, la consommation d'énergie de la balise se trouve améliorée. On décrit maintenant, en relation avec les figures la et lb, la configuration d'une balise de géolocalisation en liaison avec un dispositif de calcul de position ainsi que les délais de traitement associés, selon un mode de réalisation de l'invention. 10 La balise de géolocalisation 100 comprend un module 101 apte à recevoir et à traiter les signaux reçus de véhicules spatiaux 110 appartenant à un système de géolocalisation (e.g. GPS, Galiléo, Glonass, Beidou, etc.). Le module 101 est apte à fournir certaines informations sur les signaux reçus, entre autres l'identifiant du véhicule spatial dont le signal est reçu et traité par le 15 module, ainsi que la mesure de phase de code (« code-phase measurement » en anglais) du code associé à ce véhicule. En effet, dans tous les systèmes de géolocalisation existants, un code à étalement de spectre du type CDMA (pour « Code Division Multiple Access » en anglais) est utilisé par les véhicules spatiaux de la constellation. Plus précisément, chaque véhicule spatial de la constellation utilise un tel code pour moduler 20 les signaux qu'il émet. La mesure de phase de code correspond alors à l'information de décalage temporel entre le signal reçu par le module en provenance du véhicule spatial, et une référence temporelle interne à ce module. Selon une technique connue de l'homme du métier, le signal reçu étant modulé par le code utilisé par le véhicule spatial en question, cette information de décalage temporel peut être obtenue par 25 l'observation du pic résultant de la corrélation du signal reçu avec le code en question. Ce décalage temporel représente alors la « phase » par rapport à cette référence temporelle interne au module, du code utilisé par le véhicule spatial ayant émis les signaux à destination du module. Par ailleurs, dans la plupart des systèmes de géolocalisation, un code CDMA différent est associé à chacun des véhicules spatiaux de 30 la constellation. Dans ce cas, l'identifiant du véhicule spatial peut se déduire du code 3058533 20 ayant conduit à une corrélation non nulle, i.e. présentant un tel pic. C'est en particulier le cas dans le système GPS, comme décrit par exemple dans l'ouvrage de Mohinder S. Grewal, Lawrence R. Weill et Angus P. Andrews, « Global Positioning Systems, Inertial Navigation, and Integration », second edition, publié par Wiley (IBSN 978-0-470-04190- 5 1). Cependant d'autres techniques peuvent se rencontrer pour la détermination de cet identifiant. Par exemple, le système Glonass associe une fréquence porteuse différente à chaque véhicule spatial. Dans un mode de réalisation particulier où la balise de géolocalisation utilise le standard GPS, le module 101 peut par exemple être le module « BD970 » proposé par la 10 société Trimble (http://www.navtechgps. com/departments/equipment/receivers/oem/trimble_bd970_ gnss_receiver/). Par ailleurs, il apparaît que les phases de codes sont représentatives de la distance entre la balise réceptionnant les signaux émis par des véhicules spatiaux, et les 15 véhicules spatiaux en question, mais à une constante prêt (cette constante étant liée notamment à l'instant d'échantillonnage du signal reçu au niveau de la balise). De la sorte, les traitements classiques de détermination de la position de la balise en question ne font appel qu'à la différence entre les phases de codes disponibles. Ainsi, dans un mode de réalisation, il est tiré parti de ce traitement particulier 20 des phases de code pour la détermination de la position de la balise. Plus particulièrement, ce sont non pas directement les phases de codes qui sont encodées pour être transmises via le lien radio, mais leurs différences avec une phase de code de référence (e.g. une mesure de phase de code choisie comme référence parmi les mesures de phase de code préalablement déterminées) de manière à minimiser encore 25 plus la quantité d'information que la balise doit transmettre via le lien radio. La balise de géolocalisation comprend en outre un dispositif encodeur 102 apte à mettre en forme les données (c-à-d les informations de phases de codes du type comprenant les mesures de phases de codes ou leurs différences avec une phase de code de référence, et éventuellement les identifiants des véhicules spatiaux) générées 30 par le module 101 (après traitement des signaux reçus des véhicules spatiaux). Plus 3058533 21 précisément, Le dispositif encodeur 102 encode ces données sous la forme de trames de données de manière à pouvoir les transmettre de manière efficace, en termes de consommation d'énergie, via une transmission radio LPWAN jusqu'à la station de base 104 d'un réseau LPWAN. Ainsi, les trames de données sont fournies à un module radio 5 103 du type LPWAN, apte à communiquer avec la station de base 104 via des trames radio. Le module radio encapsule alors les trames de données dans des trames radio, i.e. les trames de données deviennent les données utiles des trames radio LPWAN, dans le but de les transmettre à la station de base. Le module radio 103 du type LPWAN peut être par exemple un module Sigfox proposé par la société Adeunis RF 10 (http://www.adeunis-rf.com/fr/produits/module-sigfox) ou un module LoRa proposé par la même société (http://www.adeunis-rf.com/fr/produits/module-rf-lora). Dans un premier temps, on décrit plus particulièrement le cas où une trame de données délivrée par le dispositif encodeur 102 et correspondant à une mesure de position de la balise (i.e. les informations de phases de codes du type comprenant les 15 mesures de phases de codes ou leurs différences avec une phase de code de référence, et éventuellement les identifiants des véhicules spatiaux, obtenus via le traitement des signaux reçu de ces véhicules spatiaux à un instant donné) constitue la partie utile d'une trame radio donnée. Dans ce cadre, selon un mode de réalisation, les trames radio reçues par la 20 station de base sont horodatées par cette dernière à l'instant ti3Ts auquel elle les reçoit. Ainsi, un serveur de réseau 105 (dédié à la gestion du réseau LPWAN dont fait partie la station de base 104) peut délivrer à un dispositif de calcul de position 120 des trames radio horodatées, c'est-à-dire associées chacune à une information d'horodatage indiquant l'instant ti3Ts auquel cette trame radio (transmise par la balise) a été reçue par 25 la station de base. Dans une variante, le serveur réseau 105 est un mini-serveur inclus dans la station de base 104. Le dispositif de calcul de position 120, par exemple embarqué dans un serveur applicatif, est apte à obtenir les éphémérides et/ou almanachs des véhicules spatiaux du système de géolocalisation, auprès d'un serveur d'éphémérides et/ou almanachs 106. 3058533 22 Le dispositif de calcul de position est également apte, pour chaque trame radio qu'il reçoit (provenant de la balise via la station de base), à déterminer l'instant de réception tR auquel le module 101 de la balise a effectivement reçu les signaux des véhicules spatiaux correspondant aux données contenues la trame de donnée qui 5 constitue la partie utile de cette trame radio. L'instant de réception tR est référencé à l'antenne de réception du module 101. Pour déterminer l'instant de réception tR, le dispositif de calcul de position combine l'instant tin-s et un délai de transit Dl dans la balise (par exemple en soustrayant l'un à l'autre, i.e. tR = tin-s - Dl). Le délai Dl est défini comme suit : D1= tT - tR, avec 10 - tR l'instant de réception, par le module 101, de signaux de véhicules spatiaux ; et - tT un instant de transmission, par le module radio 103, des données générées par le module 101 (à partir des signaux reçus) et encodées par le dispositif encodeur 102). Dans un premier exemple de réalisation, le délai de transit Dl est prédéterminé 15 (imposé) et identique pour toutes les balises (ou au moins celles transmettant leurs données à destination du même dispositif de calcul de position). Les balises doivent alors implémenter des mécanismes de garantie de ce délai, par exemple via la gestion de files d'attente des informations à transmettre, ou bien via l'utilisation d'horloges et de compteurs de temps internes à un microcontrôleur implémentant ce mécanisme. 20 Dans un deuxième exemple de réalisation, le délai de transit Dl est prédéterminé (imposé) mais sans être identique pour toutes les balises de géolocalisation. Dans ces premier et deuxième exemples de réalisation, le dispositif de calcul de position 120 connaît le délai de transit Dl de chaque balise. 25 Dans un troisième exemple de réalisation, le délai de transit Dl n'est pas préalablement connu du dispositif de calcul de position 120. Cette situation se produit par exemple si le délai de transit Dl n'est pas prédéterminé (i.e. pas fixe). C'est notamment le cas dans certains réseaux LPWAN qui demandent de transmettre plusieurs fois le même message afin d'assurer la réception d'au moins un message (c'est 30 en particulier le cas lorsqu'il s'agit de protocoles sans accusé de réception 3058533 23 systématique). À titre d'exemple, les réseaux Sigfox demandent de transmettre trois fois le même message. Plus précisément, dans un réseau Sigfox, chaque message (trame radio) dans le sens montant dure environ deux secondes et la triple transmission du même message (i.e. la transmission de trois messages identiques) prend alors six 5 secondes selon le protocole de ce réseau. Dans le troisième exemple de réalisation, afin que le dispositif de calcul de position 120 prenne connaissance du délai de transit Dl associé à une trame radio donnée, il est nécessaire que la balise transmette (en plus des données précitées : informations de phases de codes du type comprenant les mesures de phases de codes 10 ou leurs différences avec une phase de code de référence, et éventuellement identifiants des véhicules spatiaux) une information indiquant le délai de transit Dl. L'information indiquant le délai de transit Dl est par exemple : - le numéro de la retransmission (i.e. un, deux ou trois dans l'exemple ci-dessus concernant le cas d'un réseau Sigfox ou LoRa) ; ou 15 - le délai de transit effectif, i.e. le délai Dl effectivement écoulé entre l'instant tR de réception des signaux en provenance des véhicules spatiaux, et l'instant tT de transmission des données correspondantes via une trame radio sur le lien radio LPWAN. Ainsi, dans le cas d'un réseau LPWAN qui demande de transmettre trois fois le 20 même message, la balise selon le troisième exemple de réalisation transmet trois messages qui ne sont pas identiques : ils contiennent les mêmes données mais chaque message contient une information indiquant un délai de transit D1 différent (par exemple, 2s, 4s ou 6s selon qu'il s'agit du premier, second ou troisième message). Selon un mode de réalisation particulier, l'instant de réception tR peut être 25 estimé par le dispositif de calcul de manière suffisamment précise pour l'application visée en soustrayant le délai de transit Dl à l'information d'horodatage tBTs- En effet, le délai de propagation D2 de la trame radio sur le lien radio LPWAN (D2 = tBTs - tT), entre la balise de géolocalisation et la station de base, reste négligeable sur la précision finale obtenue sur la localisation de la balise. À titre d'exemple, pour une distance de 30km 30 entre la balise de géolocalisation et la station de base, le délai de propagation D2 reste 3058533 24 inférieur à 1001.1.s. Or, dans un système GPS par exemple, les véhicules spatiaux se déplacent à la vitesse par rapport à la surface de la Terre de 3,7km/s. Une erreur de 1001.1.s sur le temps de propagation induit alors une erreur dans le positionnement des véhicules spatiaux de 37cm au maximum. Cette erreur, déjà faible en soi, se trouve 5 encore réduite lors du calcul du positionnement de la balise par des effets de moyennage entre les erreurs entachant la position des différents véhicules spatiaux. Selon un mode de réalisation de l'invention, le dispositif de calcul de position a à sa disposition d'une part l'instant de réception tR (qu'il a pu calculer selon la technique décrite ci-dessus), et d'autre part les éphémérides et/ou almanachs des véhicules 10 spatiaux. Ainsi, si les données transmises par la balise de géolocalisation contiennent les informations de phases de codes (e.g. les mesures de phases de codes, ou leurs différences avec une phase de code de référence), ainsi que les identifiants d'au moins quatre véhicules spatiaux différents, le dispositif de calcul de position peut déterminer la position de la balise à l'aide d'une technique connue de l'homme du métier, comme par 15 exemple celle décrite dans l'ouvrage de Mohinder S. Grewal, Lawrence R. Weill et Angus P. Andrews précité. Dans d'autres variantes de réalisation, le dispositif de calcul de position a à sa disposition des informations additionnelles relatives à la localisation de la balise, ces informations lui permettant de réduire le nombre d'inconnues dans le système 20 d'équations qu'il a à résoudre pour déterminer la position de cette balise. Par exemple, si le dispositif de calcul de position connaît l'altitude de la balise, la connaissance des mesures de phases de codes et des identifiants de seulement trois véhicules spatiaux différents suffit pour le calcul de position. Dans d'autres exemples, de telles informations additionnelles permettent de réduire à deux ce nombre minimal de 25 véhicules spatiaux différents dont il faut connaître les mesures de phases de codes ainsi que les identifiants. De manière non limitative, on considère dans la suite de la description le cas classique où les informations relatives à au moins quatre véhicules spatiaux différents sont transmises par la balise. Alternativement, et comme détaillé par la suite, la balise de géolocalisation peut 30 ne pas fournir les identifiants des véhicules spatiaux dont les informations de phases de 3058533 25 codes (e.g. les mesures de phases de codes, ou leurs différences avec une phase de code de référence), sont transmises. Le dispositif de calcul de position détermine alors la position de la balise en testant les différentes hypothèses possibles sur les identifiants de ces véhicules spatiaux. 5 De même, il existe des variantes de réalisation dans lesquelles le dispositif de calcul de position n'a pas accès à l'information d'horodatage ti3Ts et/ou la balise ne fournit pas son délai de transit Dl. Dans ce cas, le dispositif de calcul de position ne peut reconstruire directement l'instant de réception tR. Cependant, il peut tout de même déterminer la position de la balise si celle-ci lui transmet les informations de phases de 10 codes (e.g. les mesures de phases de codes, ou leurs différences avec une phase de code de référence), d'au moins cinq véhicules spatiaux en mettant en oeuvre une technique de récupération de temps connue de l'homme du métier, par exemple celle décrite dans l'article de Syrjarinne, J., « Time recovery through fusion of inaccurate network timing assistance with GPS measurements, » publié dans Proceedings of the Third International 15 Conference on Information Fusion, vol.2, pp.WED5/3-WED5/10 vol.2, 10-13 July 2000. On décrit maintenant, en relation avec les figures 2a et 2b, les étapes d'un procédé de détermination de la position d'une balise de géolocalisation, selon des modes de réalisation particuliers de l'invention. Plus particulièrement, on décrit en relation avec la figure 2a, les étapes de ce 20 procédé qui sont mises en oeuvre dans la balise 100 dont on veut déterminer la position. Lors d'une étape 200, des signaux émis par les véhicules spatiaux appartenant à la constellation d'un système de géolocalisation sont reçus, à l'instant tR, par le module 101 de la balise. Lors d'une étape 201, pour chacun d'au moins quatre véhicules spatiaux dont 25 elle reçoit un signal, le module 101 de la balise détermine les données suivantes : une mesure de phase de code ainsi qu'un identifiant du véhicule spatial. Sur la base des mesures de phases de codes ainsi déterminées, la balise obtient des informations de phases de codes qu'elle va transmettre sur le lien radio LPWAN à destination de la station de base 104 du réseau LPWAN lors de l'étape 204 décrite ci-dessous. Dans un 30 mode de réalisation, une information de phase de code correspond directement à une 3058533 26 mesure de phase de code préalablement déterminée. Dans un autre mode de réalisation, une information de phase de code correspond à une différence de mesure de phase de code avec une phase de code de référence. Dans cet autre mode de réalisation, la balise 100 et le dispositif de calcul de position 120 partagent une 5 information représentative de la phase de code de référence afin de permettre au dispositif de calcul de position 120 d'effectuer les calculs de détermination de la position de la balise 100. Dans une implémentation particulière de cet autre mode de réalisation, la phase de code de référence est sélectionnée parmi les mesures de phase de code 10 préalablement déterminées. De la sorte, une information de phase de code se trouve être nulle et n'est pas transmise, ce qui minimise encore la quantité d'information transmise par la balise 100 sur le lien radio LPWAN à destination de la station de base 104 du réseau LPWAN. Dans cette implémentation particulière, l'information partagée entre la balise 100 et le dispositif de calcul de position 120 est par exemple 15 représentative de l'identifiant du véhicule spatial auquel correspond l'information de phase de code nulle non transmise comme décrit ci-dessous en relation avec les figures 4a à 4f. Lors d'une étape 202, le dispositif encodeur 102 de la balise encode dans une trame de données les données déterminées à l'étape 201, pour permettre leur 20 transmission sur le lien radio LPWAN via une trame radio. Comme décrit en relation avec la figure la, les données sont encodées de manière adaptative, en fonction du nombre d'informations de phases de codes à transmettre, et sous la contrainte que le nombre total de bits alloués aux informations de phases de codes reste inférieur ou égal à un nombre prédéterminé de bits. Les informations de phases de codes peuvent ainsi être 25 transmises de la manière la plus efficace, en termes de consommation d'énergie, sur le lien radio LPWAN. Des modes de réalisation particuliers de cet encodage adaptatif sont décrits ci-après en relation avec les figures 4a à 4f. Lors d'une étape 203, le dispositif encodeur 102 de la balise égalise le délai de transit Dl existant entre l'instant tR de réception par la balise des signaux émis par les 30 véhicules spatiaux, et l'instant tT de transmission de la trame radio encapsulant la trame 3058533 27 de données (dans laquelle sont encodées les données correspondantes) sur le lien radio LPWAN. En d'autres termes, la balise de géolocalisation transmet la trame radio après le délai Dl, avec : Dl = tT - tR. Comme décrit en relation avec la figure 1, ce délai Dl est : 5 - soit prédéterminé et connu du dispositif de calcul de position (dans ce cas, toutes les balises partagent le même délai Dl ou bien chaque balise possède un délai Dl différent); - soit non prédéterminé et donc non connu du dispositif de calcul de position. Dans ce cas, le dispositif encodeur 102 de la balise insère dans la trame de 10 données une information indiquant le délai de transit Dl associé à ces données, permettant ainsi au dispositif de calcul de position de connaître ce délai Dl. Dans une variante, l'information indiquant le délai de transit Dl est incluse dans la trame radio (par exemple dans son en-tête) directement par le module radio LPWAN. 15 Lors d'une dernière étape 204, la balise transmet la trame radio sur le lien radio LPWAN à destination de la station de base 104 du réseau LPWAN. On décrit maintenant, en relation avec la figure 2b, les étapes du procédé qui sont mises en oeuvre du côté réseau. Lors d'une étape 205, la station de base du réseau LPWAN reçoit la trame radio 20 transmise par la balise. Lors d'une étape 206, la trame radio reçue par la station de base est horodatée avec une information d'horodatage indiquant l'instant de sa réception tin-s. Lors d'une étape 207, le dispositif de calcul de position 120 reçoit la trame radio horodatée et la décode, de manière à retrouver (reconstruire) les informations de 25 phases de codes préalablement encodées de manière adaptative avant d'être transmises par la balise comme décrit ci-dessus en relation avec l'étape 202. De telles techniques de reconstruction des informations de phases de codes sont décrites ci-après en relation avec les figures 4a à 4f. Dans certaines variantes (détaillées ci-dessus et/ou ci-après), le dispositif de calcul de position retrouve également les identifiants des 3058533 28 véhicules spatiaux correspondant aux informations de phases de codes transmises et/ou l'information indiquant le délai de transit Dl. Lors d'une étape 208, l'instant tR de réception par la balise des signaux émis par les véhicules spatiaux est calculé par le dispositif de calcul, à partir de l'information 5 d'horodatage ti3Ts et du délai de transit Dl. Lors d'une étape 209, le dispositif de calcul obtient les éphémérides et/ou almanachs des véhicules spatiaux, par exemple auprès du serveur d'éphémérides et/ou almanachs 106. Lors d'une étape 210, le dispositif de calcul 120 calcule la position de la balise 10 sur la base de l'instant de réception calculé tR, des éphémérides et/ou almanachs des véhicules spatiaux, ainsi que les données contenues dans la trame de données reçue via la trame radio (i.e. les informations de phases de codes, ainsi qu'éventuellement les identifiants des véhicules spatiaux). On décrit maintenant, en relation avec les figures 3a et 3b, les interactions entre 15 les différentes entités mises en jeu lors de la mise en oeuvre du procédé selon les organigrammes des figures 2a et 2b. Plus particulièrement, on décrit en relation avec la figure 3a, les interactions impliquant la balise de géolocalisation. Lors d'une étape 300, des signaux émis par les véhicules spatiaux 110 20 appartenant à la constellation d'un système de géolocalisation sont reçus par le module 101 de la balise 100. Lors d'une étape 301, le module 101 de la balise détermine la mesure de phase de code ainsi que l'identifiant de chaque véhicule spatial dont il reçoit un signal. Sur la base des mesures de phases de codes ainsi déterminées, la balise obtient des 25 informations de phases de codes qu'elle va transmettre sur le lien radio LPWAN à destination de la station de base 104 du réseau LPWAN lors de l'étape 306 décrite ci-dessous. Dans un mode de réalisation, une information de phase de code correspond directement à une mesure de phase de code préalablement déterminée. Dans un autre mode de réalisation, une information de phase de code correspond à une différence de 30 mesure de phase de code avec une phase de code de référence. Dans cet autre mode de 3058533 29 réalisation, la balise 100 et le dispositif de calcul de position 120 partagent une information représentative de la phase de code de référence afin de permettre au dispositif de calcul de position 120 d'effectuer les calculs de détermination de la position de la balise 100. 5 Dans une implémentation particulière de cet autre mode de réalisation, la phase de code de référence est sélectionnée parmi les mesures de phase de code préalablement déterminées. De la sorte, une information de phase de code se trouve être nulle et n'est pas transmise, ce qui minimise encore la quantité d'information transmise par la balise 100 sur le lien radio LPWAN à destination de la station de base 10 104 du réseau LPWAN. Dans cette implémentation particulière, l'information partagée entre la balise 100 et le dispositif de calcul de position 120 est par exemple représentative de l'identifiant du véhicule spatial auquel correspond l'information de phase de code nulle non transmise comme décrit ci-dessous en relation avec les figures 4a à 4f. 15 Lors d'une étape 302, ce module 101 transmet les données brutes (données obtenues dans l'étape 301) au dispositif encodeur 102 de la balise. Lors d'une étape 303, le dispositif encodeur 102 encode les données dans une trame de données à transmettre via une trame radio sur le lien radio LPWAN. Comme décrit en relation avec les figures 1 et 2a, ces données sont encodées de manière 20 adaptative, en fonction du nombre d'informations de phases de codes à transmettre, et sous la contrainte que le nombre total de bits alloués aux informations de phases de codes reste inférieur ou égal à un nombre prédéterminé de bits. Les informations de phases de codes peuvent ainsi être transmises de la manière la plus efficace, en termes de consommation d'énergie, sur le lien radio LPWAN. Des modes de réalisation 25 particuliers de cet encodage sont décrits ci-après en relation avec les figures 4a à 4e. Lors d'une étape 304, le dispositif encodeur 102 égalise le délai de transit Dl (voir définition plus haut) ou, dans une variante, insère dans la trame de données une information indiquant ce délai de transit Dl. 3058533 30 Lors d'une étape 305, le dispositif encodeur 102 fournit la trame de données au module radio 103 du type LPWAN qui encapsule alors la trame de données dans une trame radio. Lors d'une dernière étape 306, le module radio 103 transmet la trame radio sur 5 le lien radio LPWAN à destination de la station de base 104. On décrit maintenant, en relation avec la figure 3b, les interactions impliquant les équipements du côté réseau (station de base 104 et dispositif de calcul 120 notamment). Lors d'une étape 310, la station de base 104 reçoit la trame radio transmise par 10 le module radio 103 de la balise. Lors d'une étape 311, la trame radio reçue est horodatée par la station de base 104, avec une information indiquant l'instant de réception tBTs. Lors d'une étape 312, la trame radio horodatée est transmise au serveur réseau 105 dédié à la gestion du réseau LPWAN, qui la fait suivre au dispositif de calcul de 15 position 120 lors d'une étape 313. Lors d'une étape 314, les données horodatées sont décodées par le dispositif de calcul de position 120 de manière à reconstruire les informations de phases de codes préalablement encodées de manière adaptative avant d'être transmises par la balise dont on veut déterminer la position et/ou les identifiants des véhicules spatiaux 20 correspondants. De telles techniques de reconstruction des données transmises sont décrites ci-après en relation avec les figures 4a à 4f. Dans des variantes, le dispositif de calcul reconstruit également les identifiants des véhicules spatiaux dont les informations de phases de codes ont été transmises et/ou l'information indiquant le délai de transit Dl. 25 Lors d'une étape 315, l'instant tR de réception par la balise des signaux émis par les véhicules spatiaux est calculé par le dispositif de calcul, à partir de l'information d'horodatage ti3Ts et du délai de transit Dl. Lors d'une étape 316, le dispositif de calcul de position 120 obtient les éphémérides et/ou almanachs des véhicules spatiaux, par exemple auprès du serveur 30 d'éphémérides et/ou almanachs 106. 3058533 31 Lors d'une étape 317, la position de la balise est calculée par le dispositif de calcul de position 120 sur la base de l'instant de réception calculé tR, des éphémérides et/ou almanachs des véhicules spatiaux, ainsi que des informations de phases de codes transmis ainsi qu'éventuellement des identifiants des véhicules spatiaux. 5 Dans les modes de réalisation décrits précédemment en relation avec les figures 1 à 3b, une trame de données correspond à une trame radio, i.e. une trame de données délivrée par le dispositif encodeur 102 constitue la partie utile d'une trame radio. Cependant, dans des variantes il peut être envisagé, par exemple, que les données correspondant à une mesure de position de la balise (i.e. les informations de phases de 10 codes, et éventuellement les identifiants des véhicules spatiaux, obtenus via le traitement des signaux reçu de ces véhicules spatiaux à un instant donné) sont encapsulées dans au moins deux trames radio. Ce peut être le cas, par exemple, si ces données sont encodées dans deux trames de données distinctes par le dispositif encodeur 102, mais également si les données sont encodées de manière adaptative par 15 le dispositif encodeur dans une seule trame de données qui se trouve être de longueur supérieure à la partie utile de la trame radio prévue par l'implémentation faite du standard LPWAN dans le module radio 103. Dans ce dernier cas, le module radio va répartir la trame de données qu'il reçoit du dispositif encodeur 102 de manière à la transmettre via au moins deux trames radio. 20 Dans le cas où les données correspondant à une mesure de position de la balise sont transmises via au moins deux trames radio, le dispositif encodeur 102 de la balise égalise le délai de transit Dl existant entre l'instant tR de réception par la balise des signaux émis par les véhicules spatiaux, et l'instant tT de transmission d'une trame radio de référence parmi les au moins deux trames radio encapsulant la trame de données. 25 Comme décrit précédemment, le délai Dl peut être prédéterminé et connu du dispositif de calcul de position. Dans ce cas, le dispositif de calcul de position peut déterminer l'instant de réception tR en fonction d'une part de l'information d'horodatage de la trame de référence telle que délivrée par la station de base du réseau LPWAN, i.e. de la connaissance de l'instant tBTs auquel la station de base reçoit la 30 trame de référence, et d'autre part de l'information de délai Dl. 3058533 32 Dans d'autres cas, le délai Dl peut être non prédéterminé et donc non connu du dispositif de calcul de position. La balise de géolocalisation insère alors dans la trame de données une information indiquant le délai Dl entre l'instant tR de réception par la balise des signaux émis par les véhicules spatiaux, et l'instant tT de transmission d'une 5 trame radio de référence parmi les au moins deux trames radio encapsulant la trame de données. Dans une variante, l'information indiquant le délai de transit Dl est inclus dans une des trames radio, par exemple dans la trame radio de référence, directement par le module radio LPWAN. Le dispositif de calcul de position détermine alors l'instant de réception tR en 10 fonction d'une part de l'information d'horodatage de la trame radio de référence telle que délivrée par la station de base du réseau LPWAN, et d'autre part de l'information de délai Dl incluse soit dans la trame de données, soit directement dans une des trames radio (par exemple dans un en-tête). Dans d'autres modes de réalisation, il peut être envisagé de transmettre 15 plusieurs trames de données, correspondant par exemple aux données correspondant à différentes mesures de position de la balise (i.e. les informations de phases de codes, et éventuellement les identifiants des véhicules spatiaux, obtenus via le traitement des signaux reçu des véhicules spatiaux à différents instants), dans la même trame radio. Ce peut être par exemple le cas dans un réseau LoRa dont la partie utile des trames peut 20 atteindre 59 octets comme décrit précédemment. Cette approche permet par exemple de minimiser les pertes de capacité dans le réseau LPWAN en réduisant le nombre d'en-têtes de trames radio utilisés pour une quantité d'information utiles transmises donnée. Dans ce cas, le dispositif de calcul de position détermine l'instant de réception tR en fonction d'une part de l'information d'horodatage de la trame radio telle que 25 délivrée par la station de base du réseau LPWAN, et d'autre part d'une information de délai Dl, connue du dispositif de calcul de position (lorsque cette information est prédéterminée) et spécifique à chacune des trames de données présente dans la partie utile de la trame radio. Lorsque l'information de délai Dl est non prédéterminée, et donc non connu du dispositif de calcul de position, elle est incluse, pour chacune des 3058533 33 trames de données présente dans la partie utile de la trame radio, dans la trame de données correspondante. Ces différents modes de réalisation peuvent être combinés. Par exemple, certaines trames radio émises par la balise peuvent encapsuler plusieurs trames de 5 données correspondant à différentes mesures de position de la balise à différents instants, alors que d'autres trames de données correspondant chacune à une mesure de position de la balise à un instant donné peuvent être transmises chacune via plusieurs trames radio. Dans chaque cas, les principes détaillés ci-dessus s'appliquent à chaque configuration correspondante rencontrée dans l'ensemble des trames (trames de 10 données et trames radio) destinées à transmettre l'information vers la station de base du réseau LPWAN. On décrit maintenant en relation avec les figures 4a à 4f différents modes de réalisation de l'encodage adaptatif des données devant être transmises sur le lien radio LPWAN (l'encodage étant adaptatif en fonction du nombre de données à transmettre, 15 sous la contrainte que le nombre total de bits alloués aux données reste inférieur ou égal à un nombre prédéterminé de bits), ainsi que les techniques de décodage et de levée des ambiguïtés correspondantes par le dispositif de calcul de position. Plus particulièrement, la figure 4a présente un format d'encodage binaire des mesures de phases de codes devant être transmises via le lien radio selon un mode de 20 réalisation de l'invention. Le mot binaire 402 représentant une mesure de phase de code peut être décomposé en une partie entière 4021 et une partie fractionnaire 4022, qui définissent la position d'une virgule dans le mot binaire 402. Cette virgule permet de définir le poids associé aux bits du mot binaire. Selon un mode de réalisation de l'invention, le bit de poids le plus faible, ou LSB (pour « Least 25 Significant Bit » en anglais), de la partie entière 4021 est défini comme correspondant à la durée Tc d'un chip des codes CDMA utilisés dans le système de géolocalisation en question. À titre d'exemple, dans le système GPS, cette durée chip est de Tc = 1/(1,023MHZ) = 0,9971.1.s. Alternativement, le poids du LSB peut également être interprété en termes de 30 résolution spatiale, i.e. comme la distance équivalente Dref parcourue par le signal radio 3058533 34 émis par les véhicules spatiaux à destination d'une balise pendant cette durée Tc. Reprenant l'exemple du système GPS, il apparaît alors que le LSB de la partie entière du mot binaire codant la mesure de phase de code peut être interprété comme représentant une résolution spatiale sur la distance entre la balise et le véhicule spatial 5 ayant émis le signal (et donc sur la position de ce véhicule spatial) de Dref = C*Tc = 3e9 * Tc = 293m, avec C représentant la vitesse de la lumière dans le vide. Le poids de chaque bit du mot binaire 402 peut alors se déduire de cette valeur de LSB de la partie entière, correspondant à la position de la virgule dans ce mot. Le Tableau 1 présente une résolution spatiale correspondant à un bit de rang 10 donné dans le mot binaire codant la mesure de phase de code du véhicule spatial considéré. Tableau 1 Rang du bit considéré dans le Résolution spatiale Résolution spatiale mot binaire codant une équivalente (Dref) équivalente dans le cas GPS mesure de phase de code (m) 2 bits avant la virgule Drer2A2 1172m 1 bit avant la virgule Drer2A1 586m 0 bit avant la virgule Dref 293m 1 bit après la virgule Dref/2"1 146,5m 2 bits après la virgule Dref/2"2 73,22m 3 bits après la virgule Dref/2"3 36,63m 4 bits après la virgule Dref/2"4 18,3m 5 bits après la virgule Dref/2"5 9,15m 6 bits après la virgule Dref/2"6 4,58m 7 bits après la virgule Dref/2"7 2,29m 8 bits après la virgule Dref/2"8 1,14m Selon la convention de ce mode de réalisation, seule la valeur du poids des bits change d'un système de géolocalisation à l'autre. Par exemple, le système Galiléo ayant 15 une fréquence chip de 10,23MHz, le poids du LSB de la partie entière 4021 du mot binaire 402 codant la mesure de phase de code d'un signal émis par un véhicule spatial appartenant à cette constellation est de Tc = 1/(10,23MHZ) = 0,099711s, ou Dref = 29,3m. Toutes les valeurs présentées dans le Tableau 1 peuvent alors être adaptées à ce système en appliquant la même opération. Il en va de même pour n'importe quel autre 20 système de géolocalisation (e.g. Glonass, Beidou, etc.). 3058533 Par ailleurs, dans certains modes de réalisation, ce sont non pas directement les phases de codes qui sont encodées pour être transmises via le lien radio, mais leurs différences avec une phase de code de référence. Dans une variante, la phase de code de référence peut être choisie comme étant 5 une des mesures de phases de codes déterminées par la balise. Ainsi, la différence avec elle-même conduit à un résultat nul et seulement L-1 différences de phases de codes doivent être encodées pour transmission via le lien radio quand L phases de codes devaient l'être. Par exemple, si 5 mesures de phases de codes sont disponibles au niveau de la 10 balise, chacune d'elle étant codée sur 20 bits (e.g. 10 bits pour la partie entière et 10 bits pour la partie fractionnaire), il apparaît que 5*20 = 100 bits sont nécessaires pour transmettre l'information contenue dans les 5 mesures de phases de codes via le lien radio. A contrario, les 4 différences de phases de codes obtenues par différence avec l'une d'entre elles peuvent être codées sur 21 bits chacune (si l'on considère un résultat 15 signé pour la soustraction), ce qui donne 4*21 = 84 bits au total pour transmettre l'information via le lien radio au lieu des 100 bits dans le cas de la transmission directe des 5 mesures de phases de codes. Par ailleurs, il peut même être envisagé de ne coder les 4 différences de phases de codes que sur 20 bits chacune, i.e. sans tenir compte du bit de signe. Dans ce cas, le dispositif de calcul de position 120 traite les différences de 20 phases de codes modulo 20 bits (dont 10 bits de partie entière dans l'exemple considéré), ce qui correspond à une ambiguïté de 210 * Tc * C = 297km dans le système GPS. Une telle ambiguïté peut être gérée par le dispositif de calcul de position 120 (e.g. sur la base d'informations concernant la position précédente de la balise, ou d'un emplacement grossier de cette dernière, donnée par exemple par la localisation de la 25 station de base 104 à laquelle elle transmet la trame de donnée, etc.). Par la suite, les modes de réalisation décrits ci-dessous en rapport avec les données présentées dans les tableaux 2 à 7 sont détaillés pour le cas de l'encodage adaptatif de mesures de phases de codes. Cependant, les principes divulgués d'encodage adaptatif sur un nombre prédéterminé de bits s'appliquent de manière 30 générale à l'encodage adaptatif des informations de phases de codes telles que définies 3058533 36 ci-dessus (e.g. correspondant directement à des mesures de phases de codes, ou à leur différence avec une phase de code de référence). De même, les principes divulgués pour le décodage de ces données et la levée des ambiguïtés correspondantes par le dispositif de calcul de position lors de la détermination de la position de la balise de 5 géolocalisation s'appliquent également dans le cas de la transmission de différences de mesures de phases de codes avec une phase de code de référence. On décrit maintenant en relation avec la figure 4b un premier mode de réalisation particulier de l'encodage des données devant être transmises sur le lien radio LPWAN. 10 Selon ce mode de réalisation, les données sont encodées dans une trame de données 400b destinée à constituer la partie utile d'une trame radio de niveau supérieure, i.e. la partie utile d'une trame radio au format du standard LPWAN utilisé pour transmettre les données au réseau. Plus précisément, cette trame de données 400b comprend un entête 401b 15 permettant d'identifier les véhicules spatiaux dont les mesures de phases de codes 402 sont transmises dans cette trame de données. Dans ce mode de réalisation, cet entête est formé d'un masque de bits numérotés de un à M, avec M le nombre de véhicules spatiaux existant dans la constellation du système de géolocalisation considéré. Par exemple : 20 - le rang i de chaque bit mis à un dans ce masque indique que la mesure de phase de code du véhicule spatial de la constellation dont i est l'identifiant est inclus dans la trame de données ; - un bit à zéro indique que la mesure de phase de code du véhicule spatial correspondant n'est pas inclus. 25 Dans une variante, le poids des bits peut être inversé par rapport à cette convention. À titre d'exemple, il y a 32 véhicules spatiaux dans le système GPS. Ainsi, dans un mode de réalisation particulier de l'invention, l'entête 401b est constitué de 4 octets formant un masque de 32 bits numérotés de 1 à 32. L'indice de chaque bit de ce masque 30 correspond alors à un véhicule spatial de la constellation GPS. 3058533 37 Selon d'autres modes de réalisation, il est aussi possible d'utiliser des codages plus évolués pour l'entête 401b. En effet, pour une position donnée d'une balise de géolocalisation sur Terre, seule une fraction des véhicules spatiaux appartenant à la constellation d'un système de géolocalisation donné est visible en même temps par 5 cette balise. Par exemple, dans le système GPS, huit véhicules spatiaux seulement sont visibles simultanément en moyenne. Ainsi, n'ayant qu'un nombre restreint de bits codés par exemple à la valeur « 1 » dans l'entête 401b, il est envisageable d'utiliser un codage additionnel afin de compresser également cette information d'entête. A titre d'exemple, si seulement 5 mesures de phases de codes correspondant à 5 10 véhicules spatiaux différents de la constellation GPS sont transmis dans la trame de données, le nombre d'entêtes 401b possibles est de C(32,5) = 201376 (correspondant au nombre de tirages possibles de 5 éléments parmi 32, sans répétition). Comme 218 = 262144, il apparaît que cette l'information à transmettre dans l'entête 401b pour une configuration de véhicules spatiaux donnée peut alors se coder avec un minimum de 18 15 bits à l'aide d'un codage entropique (ou codage statistique à longueur variable), par exemple un codage de Huffman. Dans une variante, le nombre de mesures de phases de codes transmises dans la trame de données est codé à l'aide d'un premier mot binaire, par exemple sur 4 bits. Les identifiants des véhicules spatiaux correspondants sont alors codés de manière différentielle, i.e. on transmet la suite des entiers représentant la 20 différence entre 2 identifiants de véhicules spatiaux consécutifs, cette suite étant elle- même codée à l'aide d'un codage entropique. Selon un autre exemple, l'entête 401b (comprenant les 5 bits indiquant les identifiants des 5 véhicules spatiaux dont on transmet les mesures de phases de codes dans la trame de données 400b) est encodé à l'aide d'un code du type CRC (pour 25 « Contrôle de Redondance Cyclique », ou « Cyclic Redundancy Check » en anglais) ou code de hachage. Le mot binaire résultant de cet encodage est alors transmis dans l'entête 401b. Une telle information n'étant que parcellaire, il apparaît que le dispositif de calcul de position 120 doit mettre en oeuvre un traitement particulier afin de déterminer les identifiants des 5 véhicules spatiaux en question. Par exemple, le 30 dispositif de calcul de position 120 détermine tout d'abord les combinaisons possibles 3058533 38 de véhicules spatiaux dont la balise de géolocalisation 100 qui a généré la trame de données a pu acquérir les signaux (généralement entre 3 et 7 véhicules spatiaux parmi les 8 à 10 qui sont visibles simultanément par la balise de géolocalisation 100). Le dispositif de calcul de position 120 calcule ensuite le mot binaire résultant de l'encodage 5 des identifiants des véhicules spatiaux pour toutes les combinaisons précédemment déterminées, en utilisant le mode de codage considéré (e.g. CRC ou code de hachage). Enfin, le dispositif de calcul de position 120 compare les mots binaires de toutes les combinaisons considérées avec le mot binaire reçu dans l'entête 401b afin de déterminer quelle est la bonne combinaison de véhicules spatiaux dont les mesures de 10 phases de codes sont transmises dans la trame de données 400b. Dans certains cas, au moins deux combinaisons d'identifiants de véhicules spatiaux pourraient conduire à un même mot binaire après encodage selon le mode de codage considéré (e.g. CRC ou code de hachage). Il est cependant notable que cette possibilité est rare, notamment lorsque le mot binaire en question est codé sur un nombre suffisant de bits (e.g. de 16 à 18 bits). 15 Par ailleurs, le dispositif de calcul de position 120 peut dans ce cas déterminer les au moins deux combinaisons d'identifiants de véhicules correspondantes et choisir celle qui est la plus pertinente sur la base d'un critère de vraisemblance (e.g. sur la base d'informations concernant la position précédente de la balise de géolocalisation 100, ou d'un emplacement grossier de cette dernière, donnée par exemple par la localisation de 20 la station de base 104 à laquelle elle transmet la trame de donnée, etc.). Dans le premier mode de réalisation particulier de l'encodage représenté sur la figure 4b, les mesures de phases de codes sont placées de manière incrémentale à la suite de l'entête. Cela signifie que le premier champ après cet entête est occupé par le mot binaire codant la mesure de phase de code du véhicule spatial correspondant à 25 l'indice du premier bit non nul de cet entête. Le deuxième champ après cet entête est occupé par le mot binaire codant la mesure de phase de code du véhicule spatial correspondant à l'indice du deuxième bit non nul de cet entête, et ainsi de suite. Dans l'exemple de la figure 4b, les mesures de phases de codes des véhicules spatiaux dont les identifiants sont 1, 3, 9, 20 et 27 sont donc présentes dans cet ordre à la suite de 30 l'entête de 32 bits. 3058533 39 Dans une variante, les mesures de phases de codes devant être transmises sont placées dans la trame de données selon un autre ordre prédéterminé, par exemple décroissant, des identifiants des véhicules spatiaux correspondants. Par ailleurs, dans le mode de réalisation où ce sont non pas directement les 5 phases de codes qui sont encodées pour être transmises via le lien radio, mais leurs différences avec une phase de code de référence, et où la phase de code de référence est sélectionnée parmi les mesures de phases de codes déterminées par la balise de géolocalisation 100, seulement L-1 différences de phases de codes sont encodées pour transmission via le lien radio. Cependant, le dispositif de calcul de position 120 doit 10 connaître l'identifiant du véhicule spatial auquel correspond la différence nulle non transmise afin de pouvoir implémenter les algorithmes connus de triangulation pour la détermination de la position de la balise de géolocalisation 100. Pour ce faire, la balise de géolocalisation 100 et le dispositif de calcul de position 120 partagent une information représentative de l'identifiant du véhicule spatial auquel correspond la 15 différence nulle non transmise (i.e. l'identifiant du véhicule spatial auquel correspond la mesure de phase de code choisie comme référence en pratique). Par exemple, dans le cas où l'entête 401b contient L bits mis à un alors que seulement L-1 différences de phases de codes sont transmises dans la trame de données correspondante, la balise de géolocalisation 100 et le dispositif de calcul de position 120 partagent une information 20 représentative du rang du bit identifiant le véhicule spatial auquel correspond la différence nulle non transmise dans la série de bits mis à un dans l'entête 401b. Dans une variante, l'information représentative du rang du bit identifiant le véhicule spatial auquel correspond la différence nulle non transmise est prédéterminée. Dans une autre variante, cette information est transmise par la balise de géolocalisation 100 au 25 dispositif de calcul de position 120 via le lien radio. Dans un mode de réalisation de l'invention, il est proposé de réduire de manière adaptative le nombre de bits sur lequel sont codées les mesures de phases de codes lorsque le nombre total de mesures de phases de codes transmises augmente. Ainsi, il est proposé de conserver un nombre total de bits permettant de transmettre toutes les 30 informations voulues qui reste inférieur, bien que le plus proche possible, de la taille de 3058533 la trame de données 400b, et ce, quel que soit le nombre de mesures de phases de codes transmises. Selon le mode de réalisation présenté dans le Tableau 2 ci-dessous, le nombre de bits alloués est identique pour chaque mesure de phase de code présente dans la 5 trame de données. Par ailleurs, ce nombre alloué par mesure de phase de code est tel que le nombre total de bits de la trame de données reste le plus proche possible tout en étant inférieur à un nombre prédéterminé de bits, ici 64 bits (8 octets) lorsque l'on tient compte d'un entête de 32 bits comme décrit ci-dessus dans le cas d'une constellation comprenant 32 véhicules spatiaux comme le GPS. Ainsi, cette trame de données peut 10 par exemple être transmise sur une seule trame radio selon le protocole Sigfox. On note par ailleurs que dans ce mode de réalisation, le bit de poids fort du mot codant chaque mesure de phase de code correspond à dix bits avant la virgule dans le format décrit ci-dessus en rapport avec la figure 4a. En effet, ce mode de réalisation étant particulièrement adapté à la transmission sur un réseau Sigfox de mesures de 15 phases de codes provenant d'un système GPS, il est notable que la longueur des codes CDMA utilisés dans un tel système est de 1023 Chips. Ainsi, la dynamique totale de la partie entière des mesures de phases de codes ne nécessite que dix bits pour être codée. Dans le mode de réalisation proposé en relation avec le Tableau 2, il apparaît alors que cette dynamique totale est préservée quel que soit le nombre de mesures de 20 phases de codes transmises. Seule la résolution est dégradée lorsque ce nombre de mesures de phases de codes transmises augmente. Le Tableau 2 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre de mesures de phases de codes transmises, et une résolution spatiale obtenue, dans un premier mode de réalisation du codage des 25 mesures de phases de codes. Tableau 2 Nombre de véhicules Nombre de bits alloués Résolution obtenue en gardant le bit de poids fort du mot codant chaque mesure de phase de code spatiaux dont les pour coder chaque correspondant à 10 bits avant la mesures de phases de mesure de phase de virgule codes sont transmises code 4 16 Dref/2"6 5 12 Dref/2"2 6 10 Dref 3058533 41 7 9 Drer2A1 8 8 Drer2A2 Cependant, il apparaît dans le codage proposé dans le Tableau 2 que la résolution spatiale sur la position des véhicules spatiaux diminue proportionnellement au nombre de mesures de phases de codes transmises au réseau. Même si la résolution obtenue sur la position de la balise restera meilleure lors de la détermination de celle-ci 5 du fait du moyennage des erreurs entachant la position des véhicules spatiaux, il peut être intéressant de garder une bonne résolution sur les mesures de phases de codes, et plutôt sacrifier les bits de poids fort de la partie entière du mot binaire codant ces mesures de phases de codes. Comme décrit ci-dessus, ces bits de poids forts de la partie entière du mot 10 binaire codant ces mesures de phases de codes donnent la dynamique maximum accessible pour les mesures de phases de codes, et donc au final sur la position de la balise. Dans le cas du GPS, nous avons vu ci-dessus que cette dynamique maximum correspond à dix bits avant la virgule, i.e. à une position modulo Dref*1024 = Dref*2^10 = 300km. 15 Il apparaît ainsi qu'un tel encodage adaptatif des mesures de phases de codes basé sur une réduction des bits de poids fort de la partie entière nécessite des traitements particuliers au niveau du dispositif de calcul de position afin de reconstruire l'information manquante et lever les ambiguïtés correspondantes sur la position de la balise ainsi déterminée par le dispositif de calcul de position. 20 Pour ce faire, il est proposé dans un mode de réalisation d'utiliser la position des stations de base du réseau LPWAN en question. En effet, il apparaît par exemple que les réseaux LPWAN fournissent dans la plupart des cas, via la position des antennes de leur réseau, une position approximative des balises présentes dans leur réseau. Plus particulièrement, le réseau Sigfox fournit les latitudes et longitudes au degré près des 25 stations de base de son réseau, ce qui permet, connaissant dans quelle cellule du réseau une balise émet, de déterminer la position de cette balise dans un rayon maximum de cent kilomètres. De même, dans le réseau LoRa la position est fournie par la position GPS des stations de base. La portée des stations de base étant de quelques kilomètres à environ dix kilomètres, la position des balises est donc aussi connue de manière 3058533 42 grossière. Cette information de positionnement grossière liée à la connaissance de la cellule réseau dans laquelle la balise de géolocalisation se trouve peut alors être transmise au dispositif de calcul de position 120. La transmission de cette information de positionnement se fait par exemple via le serveur réseau 105 qui transmet déjà la 5 trame radio horodatée. Dans d'autres variantes, la position de la balise de géolocalisation est également connue approximativement de par l'application visée. C'est par exemple le cas lors d'un déploiement de balises sur une zone localisée comme sur un site industriel. Dans ce cas, l'information de positionnement grossière peut être préalablement fournie au dispositif 10 de calcul de position 120. Ainsi, il apparaît que les bits de poids fort des dix bits avant la virgule utilisés pour représenter la dynamique maximum de la partie entière du mot codant une mesure de phase de code (correspondant à une détermination de la position modulo 300 kilomètres dans le cas GPS comme décrit ci-dessus) peuvent ne pas être utiles en 15 pratique. De même, lorsque le nombre de véhicules spatiaux est supérieur aux quatre requis, il sera possible de déterminer a posteriori les bits de poids forts de la partie entière de la mesure de phase de code à partir de l'information de position liée à la cellule du réseau dans laquelle se trouve la balise. Il peut également être envisagé de coupler cette information avec un test effectué à partir de différentes hypothèses pour 20 les bits de poids fort de cette partie entière. Dans un mode de réalisation de l'invention, il est ainsi proposé, pour un nombre de véhicules spatiaux supérieur ou égal à cinq, de limiter les bits de poids fort pour cette partie entière comme présenté dans le Tableau 3 ci-dessous. Comme dans le mode de réalisation décrit en relation avec le Tableau 2, il apparaît ici que le nombre de bits 25 alloués par mesure de phase de code est tel que le nombre total de bits de la trame de données reste le plus proche possible tout en étant inférieur à un nombre prédéterminé de bits, ici 64 bits (8 octets) lorsque l'on tient compte d'un entête de 32 bits comme décrit ci-dessus dans le cas d'une constellation comprenant 32 véhicules spatiaux comme le GPS. Ainsi, cette trame de données peut par exemple être transmis sur une 30 seule trame radio selon le protocole Sigfox. 3058533 43 Le Tableau 3 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre de mesures de phases de codes transmises, dynamique, et une résolution spatiale obtenue, dans un second mode de réalisation du codage des mesures de phases de codes. 5 Tableau 3 Nombre de Nombre de Nombre de Dynamique Nombre de Résolution véhicules bits alloués bits alloués obtenue en bits alloués obtenue spatiaux dont pour coder pour coder gardant le bit de pour coder la les mesures chaque la partie poids faible de partie de phases de mesure de entière de la partie entière fractionnaire codes sont phase de chaque correspondant à de chaque transmises code mesure de la virgule mesure de phase de phase de code code 4 16 bits 10 bits Drer2"10 6 bits Dref/2"6 5 12 bits 7 bits Drer2A7 5 bits Dref/2"5 6 10 bits 7 bits Drer2A7 3 bits Dref/2"3 7 9 bits 7 bits Drer2A7 2 bits Dref/2"2 8 8 bits 7 bits Drer2A7 1 bit Dref/2"1 Selon un autre mode de réalisation décrit en relation avec le Tableau 4 ci-dessous, il peut être intéressant de garder une résolution complète pour quatre des mesures de phases de codes transmises, ce qui permettra de déduire la position de la balise avec une bonne résolution (e.g. quelque dizaine de mètres près) tout en gardant 10 la dynamique totale pour lever toute ambiguïté de localisation. Seules les parties fractionnaires des mesures de phases de codes suivantes seront alors transmises, et sur un nombre réduit de bits, dans le but d'améliorer encore la résolution obtenue sur la position de la balise au final. Comme dans les modes de réalisation décrits ci-dessus en relation avec les 15 Tableaux 2 et 3, il apparaît ici que le nombre de bits alloués par mesure de phase de code est tel que le nombre total de bits de la trame de données reste le plus proche possible tout en étant inférieur à un nombre prédéterminé de bits, ici 64 bits (8 octets) lorsque l'on tient compte d'un entête de 32 bits comme décrit ci-dessus dans le cas d'une constellation comprenant 32 véhicules spatiaux comme le GPS. Ainsi, cette trame 20 de données peut par exemple être transmis sur une seule trame radio selon le protocole Sigfox. 3058533 44 Le Tableau 4 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre L de mesures de phases de codes transmises, dynamique, et une résolution spatiale obtenue, dans un troisième mode de réalisation du codage des mesures de phases de codes. 5 Tableau 4 A B C D E F G H I J K 4 16 bits 10 bits Drer2"10 6 bits Dref/2"6 N/A N/A N/A N/A N/A 5 15 bits 10 bits Drer2"10 5 bits Dref/2"5 4 bits 0 bits Dref 4 bits Dref/2"4 6 14 bits 10 bits Drer2"10 4 bits Dref/2"4 4 bits 0 bits Dref 4 bits Dref/2"4 7 13 bits 10 bits Drer2"10 3 bits Dref/2"3 4 bits 0 bits Dref 4 bits Dref/2"4 8 10 bits 10 bits Drer2"10 0 bits Dref 3 bits 0 bits Dref 3 bits Dref/2"3 Notation : A : Nombre de véhicules spatiaux dont les mesures de phases de codes sont transmises B : Nombre de bits alloués pour coder chacune des 4 premières mesures de 10 phases de codes transmises C : Nombre de bits alloués pour coder la partie entière de chacune des 4 premières mesures de phases de codes transmises D : Dynamique obtenue en gardant le bit de poids faible de la partie entière correspondant à la virgule 15 E : Nombre de bits alloués pour coder la partie fractionnaire de chacune des 4 premières mesures de phases de codes transmises F : Résolution obtenue pour chacune des 4 premières mesures de phases de codes transmises G : Nombre de bits alloués pour coder chacune des (L-4) mesures de phases de 20 codes restantes transmises H : Nombre de bits alloués pour coder la partie entière de chacune des (L-4) mesures de phases de codes restantes transmises I : Dynamique obtenue en gardant le bit de poids faible de la partie entière correspondant à la virgule 25 J : Nombre de bits alloués pour coder la partie fractionnaire de chacune des (L-4) mesures de phases de codes restantes transmises 3058533 K : Résolution obtenue pour chacune des (L-4) mesures de phases de codes restantes transmises Il apparaît qu'un tel encodage adaptatif des mesures de phases de codes basé sur une réduction du nombre de bits utilisés pour coder la partie entière des mesures de 5 phases de codes au-delà d'un certain nombre de mesures de phases de codes transmises nécessite des traitements particuliers au niveau du dispositif de calcul de position afin de reconstruire l'information manquante et lever les ambiguïtés correspondantes sur la position de la balise ainsi déterminée par le dispositif de calcul de position. 10 Pour ce faire, il est proposé dans un mode de réalisation d'utiliser les 4 premières mesures de phases de codes transmises pour déterminer une première estimation de la position de la balise sans problème d'ambiguïté (les 4 premières mesures de phases de codes étant transmises avec une partie entière significative). Cette estimation de la position de la balise permet ainsi au dispositif de calcul de 15 position de reconstruire la partie entière des mesures de phases de codes des (L-4) mesures de phases de codes restantes (qui sont quant à elle transmises sans partie entière). Sur la base de l'ensemble des mesures de phases de codes disponibles à l'issue de cette reconstruction (i.e. les 4 premières mesures de phases de codes et les (L-4) 20 mesures de phases de codes dont les parties entières ont été reconstruites), le dispositif de calcul de position détermine la position de la balise sur la base de L mesures de phases de codes quand bien même le nombre total de bits utilisés pour transmettre l'information sur le lien radio était réduit. On décrit maintenant en relation avec la figure 4c un deuxième mode de 25 réalisation de l'encodage des données. Seules les mesures de phases de codes sont transmises, sans aucun entête indiquant à quels véhicules spatiaux de la constellation elles appartiennent. Ainsi, la trame de données 400c est utilisée de manière optimale pour transmettre les mots binaires 402 codant les mesures de phases de codes avec le maximum de résolution et/ou de dynamique. 3058533 46 Ainsi, par exemple, dans le cas d'une trame de données pouvant être transmise via une trame radio comprenant 64 bits (8 octets) de données utiles, et pour des mesures de phases de codes dont la partie entière est codée de manière optimale sur dix bits, comme dans un système GPS, l'on obtient le codage donné dans le Tableau 5 ci- 5 dessous. Le Tableau 5 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre de mesures de phases de codes transmises, dynamique, et une résolution spatiale obtenue, dans un quatrième mode de réalisation du codage des mesures de phases de codes. 10 Tableau 5 Nombre de véhicules spatiaux dont les mesures de phases de codes sont transmises Nombre de Nombre de bits Dynamique Nombre de Résolution bits alloués alloués pour obtenue en bits alloués obtenue pour coder coder la partie gardant le bit de pour coder la chaque entière de poids faible de partie mesure de chaque mesure la partie entière fractionnaire phase de de phase de correspondant à de chaque code code la virgule mesure de phase de code 4 24 bits 10 bits Drer2"10 14 bits Dref/2"14 5 19 bits 10 bits Drer2"10 9 bits Dref/2"9 6 16 bits 10 bits Drer2"10 6 bits Dref/2"6 7 13 bits 10 bits Drer2"10 3 bits Dref/2"3 8 12 bits 10 bits Drer2"10 2 bits Dref/2"2 Il est à noter que cette approche est compatible avec celle présentée ci-dessus où l'on garde une résolution et dynamique optimale pour les quatre premières mesures de phases de codes et où l'on dégrade cette dynamique et/ou résolution pour les mesures de phases de codes suivantes. 15 Cependant, au moment de la détermination de la position de la balise, il reste nécessaire de pouvoir retrouver à quel véhicule spatial appartiennent les mesures de phases de codes obtenues afin de pouvoir effectuer le calcul de la position de la balise. Pour cela, il est possible d'envisager effectuer toutes les combinaisons possibles en sélectionnant les solutions non aberrantes. Dans des variantes, il est également 20 envisageable de réduire le nombre de combinaisons possibles, par exemple en : - disposant les mesures de phases de codes dans l'ordre croissant des identifiants de véhicules spatiaux correspondants dans la trame de donnée ; 3058533 47 - utilisant la connaissance a priori des véhicules spatiaux qui peuvent être vus à une date donnée par la balise. Par exemple, dans un système de géolocalisation donné, seuls les véhicules spatiaux étant du bon côté de la Terre à un instant donné peuvent être détectés par la balise. À titre d'exemple, 5 sous des latitudes européennes, ce nombre maximum de véhicules spatiaux détectables est de douze pour le système GPS dont la constellation comprend 32 véhicules spatiaux. On décrit maintenant en relation avec la figure 4d un troisième mode de réalisation de l'encodage des données. L'entête 401d indique maintenant simplement le 10 nombre de mots binaires 402 codant des mesures de phases de codes et qui sont transmis dans la trame de données 400d correspondante. À titre d'exemple, si le système de géolocalisation en question est le GPS, cet entête comprend 4 bits permettant d'indiquer le nombre de mesures de phases de codes transmises. En effet, même si la constellation de ce système de géolocalisation 15 comprend 32 véhicules spatiaux, seule la moitié au plus peut être détectée à un point donné de la surface de la Terre à un instant donné. Ainsi, cet entête de 4 bits permet de coder cette fraction utile de véhicules spatiaux visibles. De la sorte, pour maintenir l'ensemble des mots binaires codant les mesures de phases de codes transmises sur un format de 64 bits (8 octets) par exemple, un codage tel que présenté dans le Tableau 6 20 ci-dessous peut être envisagé. Le Tableau 6 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre de mesures de phases de codes transmises, dynamique, et une résolution spatiale obtenue, dans un cinquième mode de réalisation du codage des mesures de phases de codes. 25 Tableau 6 Nombre de Nombre de Nombre de Dynamique Nombre de Résolution véhicules spatiaux dont les mesures de phases de codes sont transmises bits alloués bits alloués obtenue en bits alloués obtenue pour coder pour coder la gardant le bit pour coder chaque partie de poids faible la partie mesure de entière de de la partie fractionnaire de chaque mesure de phase de phase de code chaque entière code mesure de correspondant phase de à la virgule code 4 23 10 bits Dref*2"10 13 Dref/2"13 3058533 48 5 18 10 bits Dref*2"10 8 Dref/2"8 6 15 10 bits Dref*2"10 5 Dref/2"5 7 13 10 bits Dref*2"10 3 Dref/2"3 8 12 10 bits Dref*2"10 2 Dref/2"2 Un tel encodage adaptatif des mesures de phases de codes sur un nombre prédéterminé de bits nécessite des traitements particuliers au niveau du dispositif de calcul de position, notamment afin de lui permettre d'associer un identifiant de véhicule spatial à une mesure de phase de code reçues quand bien même cet identifiant ne lui a 5 pas été transmis. Par exemple, si l'on sait que 5 mesures de phases de codes sont transmises par la balise au dispositif de calcul de position (sans que les identifiants des véhicules spatiaux correspondants n'aient été transmis), sachant par ailleurs qu'au maximum 12 véhicules spatiaux peuvent être vus en un point de la Terre par la balise de 10 géolocalisation, il existe C(5,12)= 792 combinaisons différentes possibles de véhicules spatiaux de la constellation correspondant aux mesures de phases de codes reçues par le dispositif de calcul de position. Dans un mode de réalisation il est proposé de déterminer une estimation de la position de la balise de géolocalisation pour chacune des combinaisons de véhicules 15 spatiaux possible. Puis, pour chacune des positions estimées, il est proposé de corréler la position estimée avec au moins une indication additionnelle de position. Dans une variante, l'indication additionnelle de position correspond à une altitude possible à la surface de la Terre pour une balise de géolocalisation (e.g. -100 mètres à +10000 mètres). Ainsi les combinaisons de véhicules spatiaux conduisant à une 20 altitude de la balise de géolocalisation qui n'est pas dans cet intervalle prédéterminé d'altitudes admissibles sont écartées. Dans une autre variante, l'indication additionnelle de position correspond à l'erreur résiduelle obtenue lorsque les mesures de phases de code d'au moins 6 véhicules spatiaux (ou au moins 5 véhicules spatiaux lorsque le temps absolu du système 25 de géolocalisation considéré est déjà connu) sont disponibles au niveau du dispositif de calcul de position. En effet, si la position estimée de la balise à partir de 5 mesures de phases de code parmi les au moins 6 (ou à partir de 5 mesures de phases de code parmi les au moins 5 lorsque le temps absolu du système de géolocalisation considéré est déjà 3058533 49 connu) est fausse, l'erreur résiduelle obtenue en prenant en compte la ou les mesures de phases de codes supplémentaire(s) est très forte. Un critère de sélection des combinaisons de véhicules spatiaux admissibles sur la base d'un seuil sur cette erreur résiduelle est alors possible. 5 Dans encore une autre variante, l'indication additionnelle de position correspond à un critère de divergence d'un algorithme de détermination du temps absolu du système de géolocalisation considéré. En effet, si la combinaison de véhicules spatiaux considérée ne correspond pas aux véhicules spatiaux dont les signaux ont été effectivement reçus par la balise, l'algorithme de détermination du temps absolu du 10 système de géolocalisation considéré va diverger. Un critère de sélection des combinaisons de véhicules spatiaux admissibles sur la base d'un seuil sur cette erreur est alors possible. Il est bien sur envisageable de combiner ces différents critère afin d'obtenir une sélection plus fiable des combinaisons de véhicules spatiaux admissibles afin de 15 déterminer la combinaison de véhicules spatiaux qui correspond effectivement aux mesures de phases de codes reçues par le dispositif de calcul de position. On décrit maintenant en relation avec la figure 4e un quatrième mode de réalisation de l'encodage des données. Un entête 401e indique maintenant un identifiant du véhicule spatial dont la mesure de phase de code est transmise dans la 20 même trame de données. Ainsi, chaque entête de ce type présent dans la trame de données est associé à un mot binaire 402 codant une mesure de phase de code. Dans un mode de réalisation, chaque entête est suivi par le mot binaire codant la mesure de phase de code correspondante. Dans d'autres modes de réalisation, le champ de phase peut être placé avant l'entête correspondant. Dans d'autres modes de 25 réalisation, les champs de phase de code sont disposés dans la trame de données dans un ordre prédéterminé, préférentiellement un ordre croissant ou décroissant, des identifiants des véhicules spatiaux dont les identifiants sont fournis via les entêtes 401e. Dans une variante, l'entête associé à chaque mot binaire codant une mesure de phase de code est un mot de 4 bits, i.e. a le même format que l'entête présenté dans un 30 mode de réalisation décrit ci-dessus en relation avec la figure 4d. Ainsi, un tel entête est 3058533 adapté à coder la fraction utile de véhicules spatiaux visibles depuis un point donné à la surface de la Terre dans un système de géolocalisation du type GPS. Par exemple, pour des mesures de phases de codes dont la partie entière est codée de manière optimale sur dix bits, comme dans le système GPS, le codage donné 5 dans le Tableau 7 ci-dessous permet d'envisager une transmission optimale sur un format de trame de données 400e de 128 bits, adapté par exemple à une trame radio du réseau LoRa. Le Tableau 7 présente un nombre de bits alloués pour coder chaque mesure de phase de code en fonction du nombre de mesures de phases de codes transmises, 10 dynamique, et une résolution spatiale, obtenue dans un sixième mode de réalisation du codage des mesures de phases de codes Tableau 7 Nombre de Nombre de Nombre de Dynamique Nombre de Résolution véhicules bits alloués bits alloués obtenue en bits alloués obtenue spatiaux dont pour coder pour coder la gardant le bit pour coder les mesures de chaque partie de poids faible la partie phases de codes mesure de entière de de la partie fractionnaire sont transmises phase de code chaque entière de chaque mesure de correspondant mesure de phase de à la virgule phase de code code 4 27 10 bits Dref 17 Dref/2"17 5 20 10 bits Dref 10 Dref/2"10 6 16 10 bits Dref 6 Dref/2"6 7 13 10 bits Dref 3 Dref/2"3 8 11 10 bits Dref 1 Dref/2"1 On décrit maintenant en relation avec la figure 4f un cinquième mode de réalisation de l'encodage des données. Un entête additionnel 401f contient une 15 information indiquant le délai de transit Dl dans le but de permettre au dispositif de calcul de position de connaître ce délai Dl, par exemple lorsque celui-ci n'est pas égalisé dans la balise. Comme décrit en relation avec la figure 1, cette information indiquant le délai de transit peut être par exemple le numéro de la retransmission (i.e. un, deux ou trois dans le cas d'un réseau Sigfox ou LoRa), ou bien le délai de transit effectif écoulé 20 entre l'instant tR de réception des signaux en provenance des véhicules spatiaux, et l'instant tT de transmission des données correspondante via une trame radio sur le lien radio LPWAN. 3058533 51 La trame de données 400f obtenue selon ce mode de réalisation comprend l'entête additionnel 401f auquel vient s'ajouter une trame de donnée 400b/c/d/e selon l'un quelconque des modes de réalisation décrits précédemment en relation avec les figures 4a, 4b, 4c et 4d. 5 Dans des variantes, l'entête additionnel 401f peut être inséré dans une trame de donnée 400b/c/d/e selon l'un quelconque des modes de réalisation décrit précédemment en relation avec les figures 4a, 4b, 4c ou 4d, plutôt que d'y être accolé. Les figures 5a, 5b et 5c présentent des exemples de structures du dispositif encodeur 102 (compris dans la balise 100), de la station de base 104 et du dispositif de 10 calcul de position 120, permettant la mise en oeuvre d'au moins un des procédés des figures 2a, 2b, 3a et 3b. Le dispositif encodeur 102, la station de base 104 ainsi que le dispositif de calcul de position 120 comprennent une mémoire vive 513, 523, 533 (par exemple une mémoire RAM), une unité de traitement 512, 522, 532, équipée par exemple d'un 15 processeur, et pilotée par un programme d'ordinateur stocké dans une mémoire morte 511, 521, 531 (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 513, 523, 533 avant d'être exécutées par le processeur de l'unité de traitement 512, 522, 532. 20 Ces figures 5a, 5b et 5c illustrent seulement une manière particulière, parmi plusieurs possibles, de réaliser le dispositif encodeur 102, la station de base 104 ainsi que le dispositif de calcul de position 120, afin qu'ils effectuent certaines étapes du procédé détaillé ci-dessus, en relation avec les figures 2a, 2b, 3a et 3b (dans l'un quelconque des différents modes de réalisation). En effet, ces étapes peuvent être 25 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 encodeur 102 et/ou la station de base 104 et/ou le 30 dispositif de calcul de position 120 est(/sont) réalisé(s) avec une machine de calcul 3058533 52 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 une disquette, un CD-ROM ou un DVD-ROM) ou non, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. 5 REVENDICATIONS1. Procédé de détermination de la position d'une balise de géolocalisation (100) connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux (110), caractérisé en ce que la balise de géolocalisation effectue les étapes suivantes : réception (200, 300) de signaux de géolocalisation provenant de N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' 2 ; détermination (201, 301) d'une mesure de phase de code pour chacun des N' véhicules spatiaux ; obtention (201, 301) de N informations de phases de codes sur la base desdites N' mesures de phases de codes, avec 2 2. Procédé selon la revendication 1, caractérisé en ce que - si ladite au moins une trame de données est transmise dans une unique trame radio, la balise de géolocalisation effectue la transmission de ladite unique trame radio après un délai (D1) prédéterminé depuis un instant de réception des signaux de géolocalisation provenant des N' véhicules spatiaux ; - si ladite au moins une trame de données est transmise dans au moins deux trames 3058533 54 radio, la balise de géolocalisation effectue la transmission d'une trame radio de référence parmi lesdites au moins deux trames radio après un délai (D1) prédéterminé depuis un instant de réception des signaux de géolocalisation provenant des N' véhicules spatiaux ; 5 ledit délai (D1) étant connu dudit dispositif de calcul de position. 3. Procédé selon la revendication 1, caractérisé en ce que : - si ladite au moins une trame de données est transmise dans une unique trame radio, la balise de géolocalisation insère dans ladite au moins une trame de données une 10 information indiquant un délai (D1) entre un instant de réception des signaux de géolocalisation provenant des N' véhicules spatiaux et un instant de transmission de ladite unique trame radio ; - si ladite au moins une trame de données est transmise dans au moins deux trames radio, la balise de géolocalisation insère dans ladite au moins une trame de données une 15 information indiquant un délai (D1) entre un instant de réception des signaux de géolocalisation provenant des N' véhicules spatiaux et un instant de transmission d'une trame radio de référence parmi lesdites au moins deux trames radio. 4. Procédé de détermination de la position d'une balise de géolocalisation (100) 20 connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux (110), caractérisé en ce qu'une station de base (104) dudit réseau de communication effectue l'étape suivante au moins une fois : 25 réception (205, 310), via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N > 2, 30 lesdites N informations de phases de codes étant encodées dans ladite au moins 3058533 55 une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phases de codes reste inférieur ou égal à un nombre prédéterminé de bits ; horodatage (206, 311) de ladite trame radio reçue, ou une desdites trames radio 5 reçues, formant une trame radio de référence, avec une information d'horodatage indiquant un instant de réception de ladite trame radio de référence par la station de base ; transmission (207, 312) de la ou les trame(s) reçue(s), y compris ladite trame radio de référence horodatée, vers un dispositif de calcul de position configuré 10 pour déterminer la position de la balise de géolocalisation en fonction des N informations de phases de codes, de l'information d'horodatage et d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux. 15 5. Procédé de détermination de la position d'une balise de géolocalisation (100) connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux (110), caractérisé en ce qu'un dispositif de calcul de position (120) effectue les 20 étapes suivantes au moins une fois : réception (207, 313), via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la 25 constellation, avec N' N > 2, lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phase de codes reste inférieur ou égal à un nombre prédéterminé de bits ; 30 décodage (207, 314) desdites N informations de phases de codes encodées ; 3058533 56 obtention (209, 316) d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux ; et détermination (210, 317) de la position de la balise de géolocalisation en fonction des N informations de phases de codes et des informations 5 d'éphéméride et/ou d'almanach. 6. Procédé selon la revendication 5, caractérisé en ce que ladite trame radio reçue, ou une desdites trames radio reçues, est une trame radio de référence horodatée et associée avec une information d'horodatage indiquant un instant de réception de ladite 10 trame radio de référence par une station de base dudit réseau de communication, et en ce que le dispositif de calcul de position : détermine en fonction de l'information d'horodatage un instant de réception tR, par la balise de géolocalisation, des signaux de géolocalisation provenant des N' véhicules spatiaux ; et 15 détermine la position de la balise de géolocalisation également en fonction dudit instant de réception tR. 7. Procédé selon la revendication 6, caractérisé en ce que le dispositif de calcul de position : 20 connaît une information prédéterminée de délai indiquant un délai (D1) entre un instant de réception, par la balise de géolocalisation, des signaux de géolocalisation provenant des N' véhicules spatiaux et un instant de transmission, par la balise de géolocalisation, de ladite trame radio de référence ; et détermine ledit instant de réception tR en fonction de l'information 25 d'horodatage et ladite information de délai. 8. Procédé selon la revendication 6, caractérisé en ce que le dispositif de calcul de position : lit dans ladite au moins une trame de données reçue une information de délai 30 indiquant un délai (D1) entre un instant de réception des signaux de géolocalisation 3058533 57 provenant des N' véhicules spatiaux et un instant de transmission de ladite trame radio de référence par la balise de géolocalisation ; et détermine ledit instant de réception tR en fonction de l'information d'horodatage et ladite information de délai. 5 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que : N 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que 10 ladite au moins une trame de données comprend au maximum kmax octets, avec kmax inférieur ou égal à 20 et préférentiellement égal à 12. 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que ladite au moins une trame de données comprend, pour chaque information de phases 15 de codes parmi lesdites N informations de phases de codes : un champ d'en-tête contenant un identifiant pour un véhicule spatial donné parmi lesdits N' véhicules spatiaux ; et un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule spatial donné. 20 12. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que ladite au moins une trame de données comprend : un champ d'en-tête, commun et indiquant un identifiant d'un véhicule spatial donné, parmi lesdits N' véhicules spatiaux, pour chacune des N informations de 25 phases de codes ; et pour chacun des N identifiants de véhicules spatiaux, un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule spatial. 30 13. Procédé selon la revendication 12, caractérisé en ce que le champ d'en-tête 3058533 58 commun est un masque comprenant une séquence de M bits, chaque bit possédant un rang au sein de ladite séquence qui correspond à un identifiant d'un des M véhicules spatiaux de la constellation, et qui est mis à 1 pour indiquer que ladite au moins une trame de données comprend un champ de phase de code pour ledit véhicule spatial et à 5 0 sinon, et en ce que l'ordre des champs de phase de code dans ladite au moins une trame de données correspond à un ordre prédéterminé, préférentiellement un ordre croissant ou décroissant, des N identifiants de véhicules spatiaux. 14. Procédé selon l'une quelconque des revendications 1 à 13, 10 caractérisé en ce qu'une desdites N' mesures de phases de codes, associée à un des N' véhicules spatiaux, est sélectionnée comme phase de code de référence, en ce que lesdites N informations de phases de codes sont N'-1 différences entre les N'1 mesures de phase de code non sélectionnées comme phase de code de référence, et ladite phase de code de référence, 15 en ce que ladite au moins une trame de données comprend lesdites N'-1 différences, et en ce que ladite balise et ledit dispositif de calcul de position partagent une information représentative d'un identifiant dudit véhicule spatial associé à ladite phase de code de référence. 20 15. Procédé selon l'une quelconque des revendications 11 à 14, caractérisé en ce que ladite au moins une trame de données ne comprend pas de champ d'en-tête contenant un ou plusieurs identifiants de véhicule spatial, en ce que, pour chacun des N véhicules spatiaux, ladite au moins une trame de données comprend un champ de phase de code contenant l'information de phase de code déterminée pour ledit véhicule 25 spatial, et en ce que l'ordre des champs de phases de codes dans ladite au moins une trame de données correspond à un ordre prédéterminé, préférentiellement un ordre croissant ou décroissant, des identifiants des N véhicules spatiaux. 16. Procédé selon la revendication 15, caractérisé en ce que le nombre N est fixe et 30 prédéterminé. 3058533 59 17. Procédé selon la revendication 15, caractérisé en ce que ladite au moins une trame de données comprend un champ contenant le nombre N. 5 18. Procédé selon l'une quelconque des revendications 11 à 17 quand elles dépendent de la revendication 10, caractérisé en ce que chacun des champs de phase de code est codé de manière adaptative en fonction de la valeur de N, en maximisant le nombre de bits utilisés, sous la contrainte que le nombre de bits alloués aux champs de phase de code reste inférieur ou égal à (8*kmax - H), avec H le nombre de bits alloués 10 au(x) champ(s) d'en-tête, ou à (8*kmax) en l'absence de champ d'en-tête. 19. Procédé selon la revendication 18, caractérisé en ce que chaque information de phase de code est composée de E bits pour la partie entière, préférentiellement E=10, et F bits pour la partie fractionnaire, et en ce que : 15 si N est inférieur à un seuil S, préférentiellement S=5, la partie entière de l'information de phase de code est codée sur E bits, pour chacun des champs de phase de code; et si N est supérieur ou égal audit seuil S, la partie entière de l'information de phase de code est codée sur E' bits, préférentiellement E'=7, résultant d'une 20 suppression d'au moins un des bits de poids fort parmi lesdits E bits, pour chacun des champs de phase de code. 20. Procédé selon la revendication 18, caractérisé en ce que chaque information de phase de code est composée de E bits pour la partie entière, préférentiellement E=10, et 25 F bits pour la partie fractionnaire, et en ce que : si N est inférieur à un seuil S, préférentiellement S=5, la partie entière de l'information de phase de code est codée sur E bits, pour chacun des champs de phase de code; et si N est supérieur ou égal audit seuil S : 30 * pour S-1 champs de phase de code, la partie entière de l'information de 3058533 60 phase de code est codée sur E bits ; et pour les champs de phase de code restants, la partie entière de l'information de phase de code est codée sur 0 bit. 5 21. Procédé selon l'une quelconque des revendications 1 à 20, caractérisé en ce que chacune desdites N informations de phases de codes appartient au groupe comprenant : une desdites N' mesures de phase de code ; et une différence entre une desdites N' mesures de phase de code et une phase de code de référence. 10 22. Produit programme d'ordinateur, comprenant des instructions de code de programme pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 21, lorsque ledit programme est exécuté sur un ordinateur. 15 23. Balise de géolocalisation (100) connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux (110), et ladite balise de géolocalisation étant caractérisée en ce qu'elle comprend : 20 des moyens de réception (101) de signaux de géolocalisation provenant de N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' 2; des moyens de détermination (101) d'une mesure de phase de code pour chacun des N' véhicules spatiaux ; 25 des moyens d'obtention (102) de N informations de phases de codes sur la base d'au moins N mesures de phases de codes parmi lesdites N' mesures de phases de codes, N' N > 2 ; des moyens d'encodage (102) des N informations de phases de codes dans au moins une trame de données sous la contrainte que le nombre total de bits 30 alloués aux N informations de phases de codes reste inférieur ou égal à un 3058533 61 nombre prédéterminé de bits, ledit encodage étant adaptatif en fonction de N ; et des moyens de transmission (103) de ladite au moins une trame de données, dans au moins une trame radio via le réseau de communication, à destination 5 d'un dispositif de calcul de position configuré pour déterminer la position de la balise de géolocalisation en fonction des N informations de phases de codes et d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux. 10 24. Station de base (104) d'un réseau de communication du type basse consommation et bas débit dédié aux objets connectés auquel est connectée une balise de géolocalisation (100), ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux (110), ladite station de base étant caractérisée en ce qu'elle comprend : 15 des moyens de réception (522), via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N > 2, 20 lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phase de codes reste inférieur ou égal à un nombre prédéterminé de bits ; des moyens d'horodatage (522) de ladite trame radio reçue, ou une desdites 25 trames radio reçues, formant une trame radio de référence, avec une information d'horodatage indiquant un instant de réception de ladite trame radio de référence par la station de base ; des moyens de transmission (522) de la ou les trame(s) reçue(s), y compris ladite trame radio de référence horodatée, vers un dispositif de calcul de position 30 configuré pour déterminer la position de la balise de géolocalisation en fonction 3058533 62 des N informations de phases de codes, de l'information d'horodatage et d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux. 5 25. Dispositif de calcul de position d'une balise de géolocalisation connectée à un réseau de communication du type basse consommation et bas débit dédié aux objets connectés, ladite balise de géolocalisation étant configurée pour recevoir des signaux de géolocalisation provenant d'une constellation de M véhicules spatiaux, ledit dispositif de calcul de position étant caractérisé en ce qu'il comprend : 10 des moyens de réception (532), via ledit réseau de communication, d'une ou plusieurs trame(s) radio comprenant au moins une trame de données comprenant N informations de phases de codes provenant de la balise de géolocalisation et relatives à N' véhicules spatiaux parmi les M véhicules spatiaux de la constellation, avec N' N > 2, 15 lesdites N informations de phases de codes étant encodées dans ladite au moins une trame de données de manière adaptative en fonction de N, sous la contrainte que le nombre total de bits alloués aux N informations de phase de codes reste inférieur ou égal à un nombre prédéterminé de bits ; des moyens de décodage (532) desdites N informations de phases de codes 20 encodées ; des moyens d'obtention (532) d'informations d'éphéméride et/ou d'almanach décrivant la constellation de véhicules spatiaux ; et des moyens de détermination (532) de la position de la balise de géolocalisation en fonction des N informations de phases de codes et des informations 25 d'éphéméride et/ou d'almanach.