FR 3141519 A1 20240503 FR 2211132 A 20221026 Procédé de recalage sur une pluralité d’amers, produit programme d'ordinateur et dispositif de recalage associés La présente invention concerne un procédé de recalage sur une pluralité d’amers lors de la géolocalisation d’un dispositif mobile, comprenant les étapes suivantes : - acquisition (110) de mesures de gisement et de site d’au moins quatre amers dans un repère du dispositif mobile ; - détermination (120) d’un écart angulaire pour chacun d’au moins trois couples desdits amers ; - détermination (130) d’une approximation initiale de la position du dispositif mobile à partir desdits écarts angulaires ; - détermination (140) d’une approximation améliorée de la position du dispositif mobile à partir de l’approximation initiale et desdits écarts angulaires ; - détermination (150) des distances entre le dispositif mobile et chacun desdits amers ; - détermination (160) d’une première approximation de l’attitude du dispositif mobile à partir desdites distances. Figure pour l'abrégé : Figure 2 Procédé de recalage sur une pluralité d’amers, produit programme d'ordinateur et dispositif de recalage associés La présente invention concerne un procédé de recalage sur une pluralité d’amers. La présente invention concerne également un produit programme d'ordinateur et un dispositif de recalage associés à ce procédé de recalage. Plus particulièrement, l’invention est applicable lors de la géolocalisation d’un dispositif mobile, notamment d’un engin mobile tel qu’un véhicule ou un aéronef, à l’aide des amers. Par « amer », on entend tout objet naturel ou artificiel repérable par le dispositif mobile et dont la position est connue ou déterminable. Selon des méthodes de relèvement déjà connues dans l’état de la technique, il est possible de géolocaliser un dispositif mobile ou de recaler sa position déterminée par un autre moyen, en utilisant une pluralité d’amers. Pour ce faire, le dispositif mobile est généralement apte à mesurer des directions (sites ou gisements) par rapport aux amers et/ou des distances jusqu’à ceux-ci. Ainsi, en analysant ces mesures et les positions connues des amers, il est possible généralement de déterminer ou de préciser la position du dispositif mobile. Toutefois, ces méthodes de relèvement, lorsqu’elles sont utilisées sans dispositif additionnel tel qu’une centrale inertielle, ne permettent généralement de résoudre le problème de géolocalisation du dispositif mobile qu’en cas bidimensionnel. En revanche, il n’existe pas de méthode connue de géolocalisation tridimensionnelle permettant à partir de mesures de relevés angulaires seuls (gisement et site) de réaliser une géolocalisation en position tridimensionnelle et en attitude tridimensionnelle (cap, roulis, tangage). La présente invention a pour but de proposer une solution permettant de géolocaliser un dispositif mobile dans l’espace tridimensionnel en utilisant seulement des mesures de relevés angulaires (gisement et site). À cet effet, l’invention a pour objet un procédé de recalage sur une pluralité d’amers lors de la géolocalisation d’un dispositif mobile comprenant les étapes suivantes : - acquisition de mesures de gisement et de site d’au moins quatre amers dans un repère du dispositif mobile ; - détermination d’un écart angulaire pour chacun d’au moins trois couples desdits amers ; - détermination d’une approximation initiale de la position du dispositif mobile à partir desdits écarts angulaires ; - détermination d’une approximation améliorée de la position du dispositif mobile à partir de l’approximation initiale et desdits écarts angulaires ; - détermination des distances entre le dispositif mobile et chacun desdits amers ; - détermination d’une première approximation de l’attitude du dispositif mobile à partir desdites distances. Suivant d’autres aspects avantageux de l’invention, le procédé comprend des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : l’écart angulaire de chaque couple d’amers est déterminé selon l’orientation respective de ces amers par rapport au dispositif mobile ; l’approximation initiale de la position du dispositif mobile est déterminée par un algorithme de type de recuit simulé ; l’approximation améliorée de la position du dispositif mobile est déterminée par un estimateur du maximum de vraisemblance, de préférence en utilisant un algorithme de type Gauss-Newton ; la première approximation de l’attitude du dispositif mobile est déterminée par un algorithme de type Kabsch en déterminant une matrice de rotation optimale entre un repère tangent et un repère du dispositif mobile ; le procédé comprend en outre une étape de détermination d’une deuxième approximation de l’attitude du dispositif mobile à partir de la première approximation de cette attitude et d’une approximation de la position du dispositif mobile, avantageusement de l’approximation améliorée de la position du dispositif mobile ; la deuxième approximation de l’attitude du dispositif mobile est déterminée par un algorithme de type Levenberg-Marquardt ; et le procédé comprend en outre une étape de détermination d’une précision de la première ou deuxième approximation de l’attitude du dispositif mobile et/ou de l’approximation améliorée de la position de ce celui-ci. La présente invention a également pour objet un produit programme d'ordinateur comportant des instructions logicielles qui, lorsque mises en œuvre par un équipement informatique, mettent en œuvre le procédé tel que défini précédemment. La présente invention a également pour objet un dispositif de recalage sur une pluralité d’amers lors de la géolocalisation d’un dispositif mobile, comprenant des moyens techniques adaptés pour mettre en œuvre le procédé tel que défini précédemment. Ces caractéristiques et avantages de l’invention apparaitront à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple non limitatif, et faite en référence aux dessins annexés, sur lesquels : - la est une vue schématique d’un dispositif de recalage selon l’invention ; et - la est un organigramme d’un procédé recalage selon l’invention, le procédé étant mis en œuvre par le dispositif de recalage de la . On a en effet représenté sur la , un dispositif de recalage 10 selon l’invention apte à déterminer ou à préciser la position d’un dispositif mobile 12 à chaque instant. Le dispositif mobile 12 est apte à se déplacer par exemple selon deux ou trois degrés de liberté par exemple dans un repère terrestre. Selon un exemple de réalisation, le dispositif mobile 12 est un engin mobile, tel que par exemple un véhicule terrestre ou un aéronef, embarquant le dispositif de recalage 10. Dans le cas d’un aéronef, il peut s’agir d’un avion ou d’un hélicoptère pilotable par un opérateur et/ou par un système avionique à partir de celui-ci, ou d’un drone pilotable à partir d’un centre de contrôle distant par un opérateur et/ou par un système avionique adapté. Selon un autre exemple de réalisation, le dispositif mobile 12 est porté par un opérateur ou par tout autre utilisateur et présente par exemple un téléphone portable. Dans ce cas, le dispositif de recalage 10 est par exemple intégré dans le dispositif mobile 12. Le dispositif de recalage 10 est apte à déterminer la position du dispositif mobile 12 en utilisant des mesures de visée prises par rapport aux amers 14-1,…,14-N par au moins certains des capteurs d’une pluralité de capteurs 16. Chaque amer 14-1,…,14-N présente un objet naturel ou artificiel disposé par exemple de manière fixe sur la surface terrestre selon une position connue avec une précision également connue et repérable depuis le dispositif mobile 12. Ainsi, chaque amer présente par exemple un pylône, une éolienne, un phare, un bâtiment, une montagne, un rocher, etc. Chaque amer 14-1,…,14-N est par exemple identifiable par un identifiant qui lui est propre. En outre, selon un mode de réalisation particulier de l’invention, l’identifiant de chaque amer 14-1,…,14-N permet d’obtenir sa position en utilisant par exemple une base de données de positions prévue à cet effet. Une telle base de données peut être publiquement accessible ou alors, présenter un accès restreint. Au moins certains des capteurs 16 sont aptes à déterminer des mesures de visée en relation avec chacun des amers 14-1,…,14-N. Avantageusement selon l’invention, chaque mesure de visée déterminée par les capteurs 16 correspondants correspond à une mesure de gisement ou à une mesure de site, par rapport à un repère de référence. Le repère de référence est par exemple lié au dispositif mobile 12 et dans le cas d’un engin mobile, présente par exemple un axe coïncidant avec le cap de cet engin. Selon un autre exemple, le repère de référence présente un repère terrestre. Comme cela est visible sur la , le dispositif de recalage 10 comprend un module d’entrée 21, un module de calcul 22 et un module de sortie 23. Le module d’entrée 21 est connecté aux capteurs 16 et est apte à acquérir l’ensemble des mesures déterminées par ces capteurs 16. Selon un exemple de réalisation, le module d’entrée 21 est connecté en outre à au moins un système, par exemple un système embarqué, apte à fournir une position connue du dispositif mobile 12. Cette position peut par exemple correspondre à la dernière position connue du dispositif mobile 12 ou la position qui a besoin d’être précisée. Le module de calcul 22 est apte à analyser les données acquises par le module d’entrée 21 afin de fournir un résultat déterminant ou précisant la position du dispositif mobile 12. Enfin, le module de sortie 23 est apte à fournir le résultat déterminé par le module de calcul 22 à tout système intéressé. Par exemple, le module de sortie 23 est apte à fournir à un écran d’affichage la position et/ou l’attitude déterminée(s) du dispositif mobile 12 pour l’afficher. Selon un autre exemple, le module de sortie 23 est apte à fournir cette position par exemple à un dispositif de précision apte à comparer cette position avec celle fournie par exemple par un système de géolocalisation de type GNSS, afin de fournir une position précise. Chacun des modules 21 à 23 se présente par exemple sous la forme d’un logiciel mis en œuvre par un processeur adapté. En variante, au moins l’un de ces modules 21 à 23 se présente au moins partiellement sous la forme d’un circuit logique programmable, tel qu’un circuit de type FPGA (de l’anglais « Field Programmable Gate Array »). Le procédé de recalage sera désormais expliqué en référence à la présentant un organigramme de ses étapes. Initialement, il est considéré que le dispositif mobile 12 se déplace à proximité des amers 14-1, …, 14-N et que dans une position fixe, les capteurs 16 effectuent M mesures de gisement et de site par rapport à au moins M amers 14-1, …, 14-N différents. Par exemple, les capteurs 16 effectuent une mesure de gisement et de site pour chacun des M amers visibles par ces capteurs 16. Avantageusement, selon l’invention, le nombre M peut être supérieur ou égal à 4 et de manière générale, peut correspondre à un nombre quelconque inférieur ou égal à N et supérieur ou égal à 4. Lors d’une étape initiale 110, le module d’entrée 21 acquiert l’ensemble des M mesures de gisement et de site, et les transmet au module de calcul 22. Par la suite, les mesures de gisement sont dénotées par et les mesures de site sont dénotées par , où correspond à l’indice associé à chacun des M amers et varie donc entre 1 et M. Lors de la même étape 110, le module d’entrée 21 acquiert en outre les positions des amers 14-1, …, 14-N ainsi que les précisions associées à ces positions. Ces données sont par exemple issues d’une base de données relative aux amers, telle que définie précédemment. Lors de l’étape 120 suivante, le module de calcul 22 détermine un écart angulaire pour M-1 couples d’amers 14-1, …, 14-N. Ces écarts angulaires sont corrélés mais ne sont pas redondantes. En outre, ces écarts sont dits « orientés » dans la mesure où ils sont déterminés selon les orientations respectives des amers correspondants par rapport au dispositif mobile 12. Pour déterminer une telle orientation pour chaque couple d’amers, le signe de la différence des mesures de gisement de ces amers est utilisé. En particulier, un écart angulaire entre un amer ayant l’indice et un autre amer ayant l’indice peut être déterminé comme suit : lorsque , et lorsque . Lors de l’étape 130 suivante, le module de calcul 22 détermine une approximation initiale de la position du dispositif mobile à partir desdits écarts angulaires et avantageusement à partir des positions connues des amers correspondants à ces écarts angulaires . Pour cela, le module de calcul 22 applique par exemple un algorithme de type de recuit simulé. En particulier, le module de calcul 22 initialise l’algorithme sur le barycentre des amers. Autrement dit, le vecteur initialisant l’algorithme peut être déterminée comme suit : où sont les positions connues des M amers. Le module de calcul 22 réalise également un choix optimisé d’un facteur d’énergie, d’un facteur de refroidissement et d’un espace de recherche aléatoire. Conformément à l’algorithme de type recuit simulé, la fonction d’énergie est choisie comme la valeur du carré des résidus : où est un écart angulaire théorique entre l’amer et l’amer lorsque le dispositif mobile 12 est dans la position . L’algorithme de type recuit simulé est itéré avec les paramètre suivants : - seuil d’énergie minimum (pour terminer un recuit) = 10 -10 ; - seuil d’énergie de sortie (pour terminer les super itérations) = 10 -6 ; - température initiale T0 = 10000° ; - nombre de paliers (itérations) max = kmax = 10000 ; - nombre maximal de tentatives (relance du recuit simulé sur la dernière solution en cas d’échec) Kmax= 10 ; - facteur d’énergie : facteur appliqué à la fonction d’énergie; - facteur de refroidissement : facteur appliqué à la température à chaque nouveau palier. La probabilité de conserver une nouvelle solution x n+1 d’énergie f(x n+1 ) supérieure à f(x n ) est exprimé comme suit : où Δe = f(x n+1 ) – f(x n ). Le module de calcul définit en outre les paramètres de facteur d’énergie et de refroidissement tels que : - au départ (T=T0) : Probabilité d’augmenter en énergie P(T) = Pmax = 0,5 ; - à la fin (T=Tmin au bout de kmax itérations de palier) : P(T) = Pmin = 10 -3 ; avec d’où : ; avec et ; d’où : Le module de calcul 22 définit en outre le voisinage de tirage aléatoire d’une solution voisine comme une zone homothétique de la zone de théâtre de rayon R évoluant entre 1/coef_Rmax = 1/10 du rayon initial Ro de la zone de théâtre (T=T0) vers Rmin = 0,1 m (pour T=Tmin) L’algorithme cœur de recuit simulé consiste à itérer les paliers en abaissant progressivement la température : . Chaque nouvelle solution x n+1 est tirée aléatoirement autour de x n avec un rayon R n = R(T n ) (différencié en horizontal et en vertical). Elle est retenue (x e = x n+1 ) si : - Δe - selon une probabilité P(Δe) = si Δe >0. Le critère d’arrêt de l’algorithme cœur est f(x n ) min et n = kmax. Le critère d’arrêt global est f(x e ) e et K = Kmax. À la fin de cette étape, le module de calcul 22 obtient ainsi une approximation initiale de la position du dispositif mobile 12. Lors de l’étape 140 suivante, le module de calcul 22 détermine une approximation améliorée de la position du dispositif mobile 12 à partir de l’approximation initiale et desdits écarts angulaires . Cette approximation améliorée est dite également par la suite première approximation de la position du dispositif mobile 12. Pour cela, le module de calcul 22 utilise un estimateur du maximum de vraisemblance, de préférence un algorithme de type Gauss-Newton. Cet estimateur a l’expression suivante : où est la nouvelle solution obtenue à partir de la solution précédente en sachant que . L’approximation améliorée de la position du dispositif mobile 12 est obtenue par itération de Gauss-Newton avec un critère d’arrêt sur la différence de deux solutions consécutives ou sur dépassement d’un nombre d’itérations maximum. Dans l’expression de l’estimateur ci-dessus, la matrice est une matrice jacobienne (sans prise en compte l’incertitude des mesures) avec une mesure d’écart angulaire relative et la matrice est une matrice d’écarts entre les écarts mesurés et les écarts théoriques correspondants. Lorsque M=4, la matrice peut prendre la forme suivante : où signifie l’écart angulaire théorique entre l’amer et l’amer selon la solution à l’itération . Lorsque M=4, la matrice peut prendre la forme suivante : Autrement dit, en tenant compte des orientations des écarts : où : et Lors de l’étape 150 suivante, le module de calcul 22 détermine des distances entre le dispositif mobile 12 et chacun M amers. Ces distances peuvent être déterminées comme suit : Ces distances permettent de déterminer les coordonnées sphériques des amers dans un repère du dispositif mobile 12, tel qu’expliqué ci-dessous, et donc les positions mesurées des amers dans ce repère. Lors de l’étape 160 suivante, le module de calcul 22 détermine une première approximation de l’attitude du dispositif mobile 12 à partir desdites positions mesurées des amers dans le repère du dispositif mobile 12 et avantageusement, à partir des positions connues des M amers. Pour cela, le module de calcul 22 détermine par un algorithme de type Kabsch en déterminant une matrice de rotation optimale entre un repère tangent et le repère du dispositif mobile 12. Le repère tangent est par exemple un repère terrestre de type NED (signifiant « North », « Est » et « Down » en anglais ou « Nord », « Est », « Bas » en français). Par la suite, les valeurs relatives à ce repère tangent sont référencées par l’indice NED et les valeurs relatives au repère du dispositif mobile 12 sont référencées par l’indice PTF. En particulier, le module de calcul 22 détermine tout d’abord les barycentres et des positions des amers dans chacun de ces repères. Autrement dit : et Puis, le module de calcul 22 détermine les écarts et de la position de chacun des amers par rapport au barycentre respectivement dans le repère du dispositif mobile 12 et le repère tangent : et La matrice de rotation optimale est déterminée grâce à la décomposition en valeurs singulières de la matrice . En particulier, en dénotant par la décomposition en valeurs singulières de la matrice , la matrice de rotation peut être déterminée comme suit : La matrice de rotation peut également s’écrire sous la forme suivante : où et sont les angles d’Euler correspondant respectivement au tangage, au roulis et au cap du dispositif mobil 12. Ces angles peuvent être déterminés comme suit : ; ; et Ces angles et forment donc une première approximation de l’attitude du dispositif mobile 12. Les étapes 170 et 180 décrites ci-dessous sont optionnellement mises en œuvre par le module de calcul 22. En particulier, lors de l’étape 170, le module de calcul 22 détermine des deuxièmes approximations de la position et de l’attitude du dispositif mobile 12 par estimation conjointe. Cette estimation est non biaisée et plus précise que la première approximation. Elle permet également d’estimer les corrélations entre coordonnées de position et coordonnées d’attitude. Autrement dit, l’étape 170 permet d’affiner les premières estimations de l’attitude et de la position du dispositif mobile 12, obtenues précédemment. Pour mettre en œuvre cette étape 170, le module de calcul 22 utilise par exemple un algorithme de type Levenberg-Marquardt. En particulier, lors de cette étape 170, le module de calcul 22 met en œuvre les itérations suivantes : où est le vecteur d’état à l’itération qui est composé du vecteur correspondant à la position tridimensionnelle du dispositif mobile 12 dans le repère tangent et du vecteur correspondant à l’attitude tridimensionnelle du dispositif mobile 12 ; est la nouvelle solution obtenue à chaque itération à partir de la solution précédente en sachant que avec et correspondant respectivement à la position et à l’attitude du dispositif mobile 12 obtenues lors des étapes précédentes ; est un paramètre d’itération avec , et ; est une matrice jacobienne ; est une matrice d’écarts des mesures réelles de visée et des mesures théoriques correspondantes. La matrice jacobienne dans le cas M=4 peut prendre la forme suivante : La matrice est déterminée comme suit : avec est la mesure vectorielle de visée de l’amer (coordonnées sphériques de cet amer dans le repère du dispositif mobile 12). Dans les expressions des deux matrices et , correspond à la mesure de visée théorique de l’amer dans le repère tangent pour l’état du dispositif mobile 12. Cette expression prend la forme suivante : où est la matrice de rotation définie précédemment ; et est le vecteur unitaire exprimant la direction dans le repère tangent NED. Lors de l’étape 180 suivante, le module de calcul 22 détermine une précision des deuxièmes approximations de l’attitude et de la position du dispositif mobile 12 déterminées lors de l’étape précédente. Cette précision peut par exemple être calculée comme une borne de Cramer-Rao ( ) en utilisant l’expression suivante : où est une matrice jacobienne globale (de dimensions 20x18 lorsque M=4) ; est une matrice de covariances des mesures. La matrice jacobienne globale est composée d’une matrice jacobienne des mesures de positions tridimensionnelles des amers et d’une matrice jacobienne des relèvements. Autrement dit : où La matrice de covariances des mesures s’écrit sous la forme suivante : où est une matrice de covariances des mesures de relèvement (de dimension 8x8 lorsque M=4) et est une matrice de covariances des mesures de position des amers (de dimension 12x12 lorsque M=4). La matrice est calculée en fonction des mesures de visée et expliquées précédemment. On conçoit alors que la présente invention présente un certain nombre d’avantages. Tout d’abord, l’invention permet de déterminer ou de préciser la position et l’attitude tridimensionnelles d’un dispositif mobile en utilisant uniquement des mesures de site et de gisement des amers, sans tout autre dispositif de détermination de l’attitude (tel qu’une centrale inertielle par exemple). En particulier, la solution proposée consiste à décomposer le problème de détermination en deux sous-problèmes : - estimation de la position tridimensionnelle par une méthode de relèvement, sans utilisation d’un dispositif spécifique de détermination de l’attitude (tel qu’une centrale inertielle par exemple) ; - estimation de l’attitude. Enfin, grâce à l’invention, il est possible d’estimer la précision de la solution obtenue. Procédé de recalage sur une pluralité d’amers (14-1, …, 14-N) lors de la géolocalisation d’un dispositif mobile (12), comprenant les étapes suivantes : - acquisition (110) de mesures de gisement et de site d’au moins quatre amers (14-1, …, 14-N) dans un repère du dispositif mobile (12) ; - détermination (120) d’un écart angulaire pour chacun d’au moins trois couples desdits amers (14-1, …, 14-N) ; - détermination (130) d’une approximation initiale de la position du dispositif mobile (12) à partir desdits écarts angulaires ; - détermination (140) d’une approximation améliorée de la position du dispositif mobile (12) à partir de l’approximation initiale et desdits écarts angulaires ; - détermination (150) des distances entre le dispositif mobile (12) et chacun desdits amers (14-1, …, 14-N) ; - détermination (160) d’une première approximation de l’attitude du dispositif mobile (12) à partir desdites distances. Procédé de recalage selon la revendication 1, dans lequel l’écart angulaire de chaque couple d’amers (14-1, …, 14-N) est déterminé selon l’orientation respective de ces amers (14-1, …, 14-N) par rapport au dispositif mobile (12). Procédé de recalage selon la revendication 1 ou 2, dans lequel l’approximation initiale de la position du dispositif mobile (12) est déterminée par un algorithme de type de recuit simulé. Procédé de recalage selon l’une quelconque des revendications précédentes, dans lequel l’approximation améliorée de la position du dispositif mobile (12) est déterminée par un estimateur du maximum de vraisemblance, de préférence en utilisant un algorithme de type Gauss-Newton. Procédé de recalage selon l’une quelconque des revendications précédentes, dans lequel la première approximation de l’attitude du dispositif mobile (12) est déterminée par un algorithme de type Kabsch en déterminant une matrice de rotation optimale entre un repère tangent et un repère du dispositif mobile (12). Procédé de recalage selon l’une quelconque des revendications précédentes, comprenant en outre une étape (170) de détermination d’une deuxième approximation de l’attitude du dispositif mobile à partir de la première approximation de cette attitude et d’une approximation de la position du dispositif mobile (12), avantageusement de l’approximation améliorée de la position du dispositif mobile (12). Procédé de recalage selon la revendication 6, dans lequel la deuxième approximation de l’attitude du dispositif mobile (12) est déterminée par un algorithme de type Levenberg-Marquardt. Procédé de recalage selon l’une quelconque des revendications précédentes, comprenant en outre une étape (180) de détermination d’une précision de la première ou deuxième approximation de l’attitude du dispositif mobile (12) et/ou de l’approximation améliorée de la position de ce celui-ci. Produit programme d'ordinateur comportant des instructions logicielles qui, lorsque mises en œuvre par un équipement informatique, mettent en œuvre le procédé selon l’une quelconque des revendications précédentes. Dispositif de recalage (10) sur une pluralité d’amers (14-1, …, 14-N) lors de la géolocalisation d’un dispositif mobile (12), comprenant des moyens techniques (21, 22, 23) adaptés pour mettre en œuvre le procédé selon l’une quelconque des revendications 1 à 8.