L’application des compensations de phase dans chaque direction visée est synchronisée avec l’application de l’intégration cohérente (11-6). L’étape (11-13) effectue une estimation au temps GNSS courant (synchronisé sur un pulse physique fourni par la référence de temps) des valeurs prédites des phases de code attendues correspondant aux retards de propagation t q des signaux à la position géoréférencée de la prochaine balise virtuelle, au temps GNSS courant t i . La synchronisation du code local (pour chacune des directions de voie) est réalisée par commande du retard de code (11-13) à appliquer sur la génération d'un signal du code local GNSS (modulation type BPSK) par un oscillateur à commande numérique NCO code (11-14), et calculée à partir de l’étape (11-1) de détermination des valeurs attendues en retard Le NCO code (11-14) réalise la synchronisation du code local GNSS sur la valeur du temps GNSS reçu (phase de code) prédit au temps courant pour la position de la balise virtuelle. Cette estimation est réalisée à chaque instant (cadencé par le 1PPS) pour chacun des satellites GNSS visible depuis la position de la balise virtuelle, par calcul de la distance apparente précises vue de la balise virtuelle ainsi que des retards supplémentaires dus à la propagation ionosphérique, troposphériques, des retards de groupe liés à l’émetteur, au récepteur et aux défauts de synchronisation résiduels entre les satellites émetteurs et le récepteur. L’étape (11-15) réalise le filtrage adapté en phase de code selon un des principes de l’invention décrit en , qui utilise le retard prédit des signaux des satellites vu de la balise virtuelle, incluant le retard de propagation ainsi que les biais connus du système (retards de groupe, bras de levier, erreur de synchronisation). Le signal reçu S(t) est multiplié par un multiplicateur (11-15) par le signal généré (11-14) par l'oscillateur à commande numérique de code. L’étape (11-6) réalise l’intégration cohérente après corrélation du signal reçu et du signal local attendu selon la fréquence, le gisement et le retard. La durée de l’intégration cohérente doit être compatible de la stationnarité de phase des signaux reçus et spécifiquement du rythme de transmission des données du message de navigation des signaux GNSS (50 bit/s dans le cas GPS). Elle détermine aussi la longueur apparente de l’antenne synthétique, aussi fonction de la vitesse de déplacement du train. Par exemple, la durée de corrélation élémentaire des codes PRN GPS étant de 1ms, l’équivalence d’un échantillonnage spatial à Shannon ( A l’étape (11-12), l’intégration non cohérente est réalisée sur chaque voie en gisement, après intégration cohérente et détection quadratique (11-11). La puissance reçue peut être alors estimée par sommation des sorties de voie en gisement, sous l’hypothèse d’un déplacement rectiligne pendant la durée d’intégration non cohérente (11-12). En sortie de l’intégration non cohérente (11-12), on dispose de la puissance du signal reçu dans chacune des directions de voies en gisement, adaptée à différentes durées d’intégration non cohérentes. L’étape (11-21) effectue le cumul des puissances intégrées, estimées pour chaque direction gisement, réalisant ainsi un filtrage adapté multivoies adaptées aux directions d’arrivées en gisement des signaux satellite au passage à la position de la balise virtuelle. Lorsque le train arrive à proximité de la position de la balise virtuelle, les directions du filtrage spatial en gisement coïncident avec celles des signaux reçus, ainsi que les retards et fréquences de code et de porteuse. La fonction de corrélation ainsi constituée mesure donc la vraisemblance de la position qui passe par un maximum à la position exacte de la balise virtuelle. Le filtrage spatial réalisé par le traitement d’antenne synthétique apporte un gain supplémentaire sur le rapport signal à bruit en réception des signaux GNSS qui participe à augmenter la précision et la disponibilité de la détection. Le filtrage spatial contribue aussi à améliorer la robustesse du filtrage adapté classique des signaux GNSS en retard et fréquence, en réduisant l’impact des interférences et trajets multiples qui arrivent en dehors du lobe de directivité des voies formées. Dans les zones plus éloignées de la position de la balise virtuelle, l’existence de coïncidences simultanées et stables d’artefacts susceptibles d’apparaitre à la fois en gisement, retard et fréquence est extrêmement rare, et conduit à réduire le risque de non intégrité (risque de détecter une position erronée ou ambiguë sans avertissement par le calculateur électronique. L’étape (11-22) réalise la détection de passage embarquée par un maximum de vraisemblance et la datation précise de ce passage selon différentes approches possibles (par exemple mettant en œuvre un enregistrement des sorties, une interpolation temporelle, un couplage avec un capteur de vitesse externe) non détaillées ici. Applications industrielles possibles de l’invention