Brevet FR3114164A1 — publié le 20220318

Procédé de contrôle d’intégrité du recalage sur une pluralité d’amers, produit programme d'ordinateur et dispositif de contrôle d’intégrité associés Ce procédé comprenant les étapes de détermination (120) d’un vecteur d’état estimé à partir d’un vecteur de mesures comprenant des mesures de position des amers, des mesures de direction des amers et des mesures de distance; de détermination (130) d’une variable de test ; de détermination (140) d’une valeur de décentralisation de la loi de distribution de ladite variable, détermination (150) d’un biais minimum détectable à partir de la valeur de décentralisation, de la variable de test et d’un vecteur caractéristique associé à ce type de défaillance ; et de détermination (160) d’un rayon de protection à partir du ou de chaque biais minimum détectable. Figure pour l'abrégé : Figure 3 Procédé de contrôle d’intégrité du recalage sur une pluralité d’amers, produit programme d'ordinateur et dispositif de contrôle d’intégrité associés La présente invention concerne un procédé de contrôle d’intégrité du recalage sur une pluralité d’amers. La présente invention concerne également un produit programme d'ordinateur et un dispositif de contrôle d’intégrité associés à ce procédé. 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 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 (azimuts 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. Les méthodes de géolocalisation par amers restent ainsi utilisables lorsque par exemple la géolocalisation par un système GNSS (de l’anglais « Global Navigation Satellite System ») est difficile ou non-souhaitable, ou lorsque par exemple il est nécessaire de déterminer la position du dispositif mobile en utilisant plusieurs sources redondantes. En outre, en cas d’utilisation par exemple des mesures inertielles pour géolocaliser le dispositif mobile, il est possible de recaler ces mesures en utilisant la géolocalisation par amers. Comme toutes les méthodes de géolocalisation se basant sur des mesures extérieures, la géolocalisation par amers présente un risque de perte d’intégrité. Cette perte d’intégrité est notamment due aux erreurs et bruits accumulées lors de la prise des mesures, et conduit généralement à une géolocalisation erronée du dispositif mobile. Toutefois, contrairement par exemple aux méthodes de géolocalisation par satellites GNSS, les méthodes de géolocalisation par amers connues ne contrôlent pas l’intégrité de la solution fournie. Ainsi, une telle solution peut comporter des erreurs qui ne seront pas détectées lors la mise en œuvre de ces méthodes. La présente invention a pour but de proposer une méthode permettant de contrôler l’intégrité des mesures utilisées lors de la géolocalisation d’un dispositif mobile par amers et en cas d’une perte d’intégrité, en avertir l’opérateur exploitant le dispositif mobile ou tout autre système. À cet effet, l’invention a pour objet un procédé de contrôle d’intégrité du recalage sur une pluralité d’amers lors de la géolocalisation d’un dispositif mobile suite à une ou plusieurs défaillances de types prédéterminés. Le procédé comprenant les étapes suivantes : - acquisition d’une première valeur de probabilité correspondant à la probabilité de détection d’une défaillance et d’une deuxième valeur de probabilité correspondant à la probabilité d’une fausse alerte ; - détermination d’un vecteur d’état estimé à partir d’un vecteur de mesures comprenant des mesures de position des amers, des mesures de direction des amers et des mesures de distance jusqu’aux amers à partir du dispositif mobile ; - détermination d’une variable de test correspondant à une distance entre le vecteur de mesures et des mesures théoriques associées au vecteur d’état estimé, la variable de test suivant une loi de distribution prédéterminée ; - détermination d’une valeur de décentralisation de la loi de distribution pour les première et deuxième valeurs de probabilité; - pour le ou chaque type de défaillance, détermination d’un biais minimum détectable à partir de la valeur de décentralisation, de la variable de test et d’un vecteur caractéristique associé à ce type de défaillance ; - détermination d’un rayon de protection à partir du ou de chaque biais minimum détectable, l’intégrité des mesures dans ledit rayon de protection étant garantie avec la première valeur de probabilité en absence d’une fausse alerte avec la deuxième valeur de probabilité. Suivant d’autres aspects avantageux de l’invention, le procédé de contrôle d’intégrité comprend des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - le vecteur d’état estimé est déterminé à partir d’un modèle linéarisé ; - le vecteur d’état estimé est déterminé en utilisant un estimateur du maximum de vraisemblance ; - chaque type de défaillance est choisi dans le groupe comprenant : - biais scalaire identique sur toutes les mesures de direction ; - biais scalaire identique sur toutes les mesures de distance ; - biais vectoriel sur les mesures de position d’au moins un amer dans un repère prédéterminé ; - pour le type de défaillance correspondant au biais vectoriel sur les mesures de position d’au moins un amer, un biais minimum détectable scalaire est déterminé pour chaque axe dudit repère ; - ladite distance entre le vecteur de mesures et les mesures théoriques associées au vecteur d’état estimé est le carré de la distance de Mahalanobis ; - ladite loi de distribution prédéterminée est la loi de chi 2 avec degrés de liberté ; - le nombre correspond au nombre total de mesures diminué de où est le nombre d’amers ; - la détermination de la valeur de décentralisation comprend : - détermination d’un seuil d’alerte en intégrité à partir de la deuxième valeur de probabilité et du nombre et en utilisant la fonction de répartition de ladite loi de distribution ; - détermination de la valeur de décentralisation à partir de la première valeur de probabilité, du nombre , du seuil d’alerte en intégrité et en utilisant la fonction de répartition de ladite loi de distribution ; - une étape de déclenchement d’une alerte lorsque le rayon de protection dépasse un seuil de protection ; - lors de l’étape de déclenchement d’une alerte, une alerte est déclenchée en outre lorsque la valeur de décentralisation dépasse le seuil d’alerte en intégrité. 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 contrôle d’intégrité du recalage sur une pluralité d’amers lors de géolocalisation d’un dispositif mobile suite à une ou plusieurs défaillances de types prédéterminés, 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 système de géolocalisation par amers d’un dispositif mobile, le système de géolocalisation comprenant un dispositif de contrôle d’intégrité selon l’invention ; - la est une vue schématique du dispositif de contrôle d’intégrité de la ; et - la est un organigramme d’un procédé de contrôle d’intégrité selon l’invention, le procédé étant mis en œuvre par le dispositif de contrôle d’intégrité de la figure 1. On a en effet représenté sur la , un système de géolocalisation 10 par amers d’un dispositif mobile 12. 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 système de géolocalisation 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 système de géolocalisation 10 est par exemple intégré dans le dispositif mobile 12. Le système de géolocalisation 10 est apte à déterminer la position du dispositif mobile 12 en utilisant notamment amers 14-1,…,14-N. 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. En référence à la , le système de géolocalisation 10 comprend une pluralité de capteurs 16, un dispositif de traitement 18 et un dispositif de contrôle d’intégrité 20. Les capteurs 16 sont aptes à déterminer des données de mesure en relation avec chacun des amers 14-1,…,14-N. Avantageusement selon l’invention, chaque mesure déterminée par les capteurs 16 correspond à l’un de type de mesure choisi dans le groupe : - mesure de position ; - mesure de direction ; - mesure de distance. En particulier, chaque mesure de position correspond à une mesure de position bidimensionnelle ou tridimensionnelles de l’un des amers 14-1,…,14-N. Une telle mesure est par exemple déterminée par l’un des capteurs 16 en déterminant l’identifiant de l’amer correspondant et en identifiant la position correspondant à cet identifiant dans la base de données de positions définie ci-dessus. Dans ce cas, le capteur correspondant se présente par exemple sous la forme d’un capteur d’images associé à une unité de traitement apte à accéder à la base de données de positions. Selon un autre exemple de réalisation, l’indentification d’au moins certains des amers s’effectue manuellement par l’opérateur. Dans ce cas, l’opérateur peut par exemple introduire manuellement l’identifiant de l’amer correspondant ou alors sa position dans le dispositif de contrôle d’intégrité 20. Selon un autre exemple de réalisation, la position d’au moins l’un des amers 14-1,…,14-N est déterminée en utilisant des informations transmises par l’amer lui-même en utilisant par exemple des signaux radioélectriques. Chaque mesure de direction correspond à une mesure d’azimut ou à une mesure de site ou à une mesure de gisement relative à l’un des amers 14-1,…,14-N. Une telle mesure est par exemple déterminée par l’un des capteurs 16 en utilisant des mesures optiques et dans le cas d’une mesure d’azimut, en comparant ces mesures avec le cap du dispositif mobile 12 ou avec une direction fixe dans l’espace. Chaque mesure de distance correspond à une mesure de distance entre le dispositif mobile 12 et l’un des amers 14-1,…,14-N. Une telle mesure est réalisée par l’un des capteurs 16 comprenant par exemple un radar ou un laser ou tout autre dispositif de télémétrie connu en soi. Le dispositif de traitement 18 est connecté aux capteurs 16 et apte à analyser l’ensemble des mesures fournies par ces capteurs 16 pour déterminer la position du dispositif mobile 10 en mettant en œuvre une méthode de géolocalisation. Une telle méthode de géolocalisation est connue en soi et ne sera pas décrite en détail par la suite. Selon un autre exemple de réalisation, le dispositif de traitement 18 est apte à déterminer ou à préciser la position du dispositif mobile 10 en utilisant en outre des données issues d’autres systèmes, comme par exemple un système de géolocalisation par satellites (système GNSS), selon des méthodes également connues en soi. Le dispositif de traitement 18 est connecté en outre à un écran permettant par exemple d’afficher la position déterminée. Selon un autre exemple de réalisation, le dispositif de traitement 18 est connecté à tout autre système apte à utiliser la position déterminée pour par exemple piloter le dispositif mobile 10 lorsqu’il présente un engin mobile. Le dispositif de contrôle d’intégrité 20 permet de contrôler l’intégrité de la position fournie par le dispositif de traitement 18, en mettant en œuvre le procédé de contrôle d’intégrité selon l’invention. En particulier, le dispositif de contrôle d’intégrité 20 est apte détecter une perte d’intégrité de la position fournie par le dispositif de traitement 18, suite à une ou plusieurs défaillances de types prédéterminés. Chaque type de défaillance est choisi dans le groupe comprenant : - biais scalaire identique sur toutes les mesures de direction (résultant par exemple de l’apparition d’un biais sur le cap) ; - biais scalaire identique sur toutes les mesures de distance (résultant par exemple d’un biais occasionné par un défaut de fonctionnement du capteur de télémétrie) ; - biais vectoriel sur les mesures de position d’au moins un amer (résultant par exemple d’une erreur de déréférencent de l’amer). Pour identifier chaque type de défaillance par la suite, l’indice variant de à est associé à chaque type de défaillance. Lorsque , il s’agit d’un biais scalaire identique sur toutes les mesures de direction. Lorsque , il s’agit d’un biais scalaire identique sur toutes les mesures de distance. Lorsque , il s’agit d’un biais vectoriel sur les mesures de position d’au moins un amer. En référence à la , le dispositif de contrôle d’intégrité 20 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. Le module d’entrée 16 est apte en outre à acquérir des paramètres de contrôle présentant des valeurs de paramétrage du procédé de contrôle d’intégrité. Selon un exemple de réalisation, ces paramètres de contrôle comprennent : - une première valeur de probabilité correspondant à la probabilité de détection d’une défaillance de l’un des types prédéterminés, la valeur signifiant donc le risque de perte d’intégrité ; - une deuxième valeur de probabilité correspondant à la probabilité d’une fausse alerte signifiant une perte d’intégrité ; - un seuil de protection signifiant un rayon de protection maximal dans lequel l’intégrité est garantie avec la première valeur de probabilité en absence d’une fausse alerte avec la deuxième valeur de probabilité . Les paramètres de contrôle sont par exemple transmis au module d’entrée 21 via une interface homme-machine ou par tout autre système connecté à ce module. 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 de contrôle d’intégrité de la position déterminée par le dispositif de traitement 18, comme cela sera expliqué plus en détail par la suite. Enfin, le module de sortie 23 est apte à fournir le résultat de contrôle d’intégrité à tout système intéressé. Par exemple, lorsque la position du dispositif mobile 10 est affichée sur un écran, le module de sortie 23 est apte afficher sur le même écran un indicateur d’intégrité de cette position. En variante ou en complément, le module de sortie 23 est apte à lever une alerte en cas d’une perte d’intégrité. 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 contrôle d’intégrité 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 10 se déplace dans l’espace par exemple selon deux degrés de liberté à proximité des amers 14-1,…,14-N dont les positions bidimensionnelles sont connues. Il est également considéré que les capteurs 16 déterminent l’ensemble des mesures relatives à ces amers, telles que définies précédemment, et que le dispositif de traitement 18 détermine la position du dispositif mobile 10 en utilisant ces mesures. Lors de l’étape initiale 110, le module d’entrée 21 acquiert l’ensemble des mesures déterminées par les capteurs 16 ainsi que les paramètres de contrôle, tels que définis précédemment. En particulier, il est considéré par la suite que le module d’entrée 21 acquiert mesures. Dans ce cas, le nombre est déterminé selon l’expression suivante : où est le nombre d’amers et donc le nombre de positions bidimensionnelles de ces amers 14-1,…,14-N ; est le nombre de mesures de direction (par exemple des mesures d’azimut) des amers 14-1,…,14-N distincts ; et est le nombre de mesures de distance jusqu’aux amers 14-1,…,14-N distincts ; est le nombre de mesures conjointes, c’est-à-dire le nombre à la fois de mesures de direction et de mesures de distance sur un même amer 14-1,…,14-N. À la fin de l’étape 110, le module d’entrée 21 transmet l’ensemble des données acquises au module de calcul 22. Lors de l’étape suivante 120, le module de calcul 22 estime un vecteur d’état pour obtenir un vecteur d’état estimé à partir d’un vecteur de mesures . Cette estimation est faite en utilisant un modèle linéarisé selon lequel : avec où est un biais produit par une défaillance (de l’un des types prédéterminés) ; sont les mesures de direction ; sont les mesures de distance ; sont les mesures de position ; est une matrice jacobienne ; est un bruit aléatoire supposé être un vecteur aléatoire gaussien ; est la position bidimensionnelle du dispositif mobile 12 ; sont des postions bidimensionnelles des amers 14-1,…,14-N. Le vecteur d’état estimé est déterminé par application d’un estimateur du maximum de vraisemblance au modèle linéarisée et se présente sous la forme suivante : où est une matrice pseudo-inverse pondérée de ( ) ; est une matrice de covariance des mesures ; est l’écart-type des mesures de position ; est l’écart-type des mesures de distance ; et est l’écart-type des mesures de direction. Lors de l’étape suivante 130, le module de calcul 22 détermine une variable de test correspondant à une distance entre le vecteur de mesures et des mesures théoriques associées au vecteur d’état estimé . La variable de test dépend du bruit ainsi que du biais . Dans l’exemple décrit, ladite distance est le carré de la distance de Mahalanobis déterminée selon l’expression suivante : où est le vecteur des résidus linéarisés égal à Ainsi, la variable de test prend la forme suivante : où La variable de test suit une loi de distribution prédéterminée qui, dans l’exemple décrit, est la loi à degrés de liberté et décentrée selon une valeur de décentralisation . Le nombre est égal à . L’espérance mathématique de la variable de test est avec . Lors de l’étape 140 suivante, le module de calcul 22 détermine la valeur de décentralisation de la loi de distribution pour les première et deuxième valeurs de probabilité. Pour cela, le module de calcul 22 détermine d’abord un seuil d’alerte en intégrité à partir de la deuxième valeur de probabilité et du nombre et en utilisant la fonction de répartition de ladite loi de distribution. En particulier, et Puis, le module de calcul 22 détermine la valeur de décentralisation à partir de la première valeur de probabilité , du nombre , du seuil d’alerte en intégrité et en utilisant la fonction de répartition de ladite loi de distribution. En particulier, Lors de l’étape 150 suivante, pour chaque type de défaillance et chaque mesure , le module de calcul 22 détermine un biais minimum détectable à partir de la valeur de décentralisation , de la variable de test et d’un vecteur caractéristique associé à ce type de défaillance. Lorsque ou , il s’agit d’un biais minimum détectable scalaire qui est déterminé en résolvant l’équation en utilisant une méthode connue en soi, comme par exemple la méthode de résolution par dichotomie. Dans ce cas, lors d’une défaillance de type , le biais minimum détectable scalaire sur une mesure peut être exprimé sous la forme suivante : où le vecteur binaire caractéristique associé au type de défaillance pour une mesure . Lors d’une défaillance avec un biais inferieur à , l’erreur maximale de position a pour expression : Lorsque , il s’agit d’un biais minimum détectable vectoriel . Ce biais sur la mesure de position peut être considéré comme une somme pondérée d’un biais sur l’abscisse et d’un biais sur l’ordonnée . Ce biais peut être alors identifié en utilisant l’expression suivante : où et sont des vecteurs caractéristique pour respectivement l’abscisse et l’ordonnée . Les biais et peuvent être déterminés en utilisant, tout comme dans le cas précédent, l’expression : En particulier, le biais peut être identifié en fonction du biais par l’expression suivante : où est l’élément de la matrice associé en ligne à la valeur et en colonne à la valeur ; est l’élément de la matrice associé en ligne à la valeur et en colonne à la valeur ; est l’élément de la matrice associé en ligne à la valeur et en colonne à la valeur . Le biais vectoriel minimum détectable de la mesure de position qui impacte le plus fortement l’estimation peut être déterminé en recherchant le biais qui maximise le biais d’estimation avec : avec étant ensuite déterminé selon la formule précédente. Lors de l’étape suivante 160, le module de calcul 22 détermine un rayon de protection dans lequel l’intégrité est garantie avec la première valeur de probabilité en absence d’une fausse alerte avec la deuxième valeur de probabilité . Ce rayon est déterminé selon l’expression suivante : où est la projection du biais minimal détectable (pour la mesure lors de la défaillance de type ) sur l’abscisse de l’estimateur de position ; est la projection du biais minimal détectable sur l’ordonnée. À la fin de cette étape 160, le module de calcul 22 transmet au module de sortie 23 la valeur qui l’affiche par exemple sur l’écran correspondant et la transmet à tout autre système intéressé. Lors de l’étape suivante 170 optionnelle, le module de calcul 22 déclenche une alerte lorsque le rayon de protection dépasse le seuil de protection ou lorsque la valeur de décentralisation dépasse le seuil d’alerte en intégrité . Dans ce cas, la solution fournie par le dispositif de traitement 18 n’est alors plus considérée comme intègre. Procédé de contrôle d’intégrité du recalage sur une pluralité d’amers (14-1, …, 14-N) lors de la géolocalisation d’un dispositif mobile (12) suite à une ou plusieurs défaillances de types prédéterminés, comprenant les étapes suivantes : - acquisition (110) d’une première valeur de probabilité ( ) correspondant à la probabilité de détection d’une défaillance et d’une deuxième valeur de probabilité ( ) correspondant à la probabilité d’une fausse alerte ; - détermination (120) d’un vecteur d’état estimé ( ) à partir d’un vecteur de mesures ( ) comprenant des mesures de position des amers (14-1, …, 14-N), des mesures de direction des amers (14-1, …, 14-N) et des mesures de distance jusqu’aux amers (14-1, …, 14-N) à partir du dispositif mobile (12) ; - détermination (130) d’une variable de test ( ) correspondant à une distance entre le vecteur de mesures ( ) et des mesures théoriques associées au vecteur d’état estimé ( ), la variable de test ( ) suivant une loi de distribution prédéterminée ; - détermination (140) d’une valeur de décentralisation ( ) de la loi de distribution pour les première et deuxième valeurs de probabilité ( ) ; - pour le ou chaque type de défaillance, détermination (150) d’un biais minimum détectable ( ) à partir de la valeur de décentralisation ( ), de la variable de test ( ) et d’un vecteur caractéristique ( ) associé à ce type de défaillance ; - détermination (160) d’un rayon de protection ( ) à partir du ou de chaque biais minimum détectable ( ), l’intégrité des mesures dans ledit rayon de protection ( ) étant garantie avec la première valeur de probabilité ( ) en absence d’une fausse alerte avec la deuxième valeur de probabilité ( ). Procédé selon la revendication 1, dans lequel le vecteur d’état estimé ( ) est déterminé à partir d’un modèle linéarisé. Procédé selon la revendication 2, dans lequel le vecteur d’état estimé ( ) est déterminé en utilisant un estimateur du maximum de vraisemblance. Procédé selon l’une quelconque des revendications précédentes, dans lequel chaque type de défaillance est choisi dans le groupe comprenant : - biais scalaire identique sur toutes les mesures de direction ; - biais scalaire identique sur toutes les mesures de distance ; - biais vectoriel sur les mesures de position d’au moins un amer (14-1, …, 14-N) dans un repère prédéterminé. Procédé selon la revendication 4, dans lequel pour le type de défaillance correspondant au biais vectoriel sur les mesures de position d’au moins un amer (14-1, …, 14-N), un biais minimum détectable scalaire ( ) est déterminé pour chaque axe dudit repère. Procédé selon l’une quelconque des revendications précédentes, dans lequel ladite distance entre le vecteur de mesures ( ) et les mesures théoriques associées au vecteur d’état estimé ( ) est le carré de la distance de Mahalanobis. Procédé selon l’une quelconque des revendications précédentes, dans lequel ladite loi de distribution prédéterminée est la loi de chi 2 avec degrés de liberté. Procédé selon la revendication 7, dans lequel le nombre correspond au nombre total de mesures ( ) diminué de où est le nombre d’amers (14-1, …, 14-N). Procédé selon la revendication 7 ou 8, dans lequel la détermination de la valeur de décentralisation ( ) comprend : - détermination d’un seuil d’alerte en intégrité ( ) à partir de la deuxième valeur de probabilité ( ) et du nombre et en utilisant la fonction de répartition ( ) de ladite loi de distribution ; - détermination de la valeur de décentralisation ( ) à partir de la première valeur de probabilité ( ), du nombre , du seuil d’alerte en intégrité ( ) et en utilisant la fonction de répartition ( ) de ladite loi de distribution. Procédé selon l’une quelconque des revendications précédentes, comprenant en outre une étape (170) de déclenchement d’une alerte lorsque le rayon de protection ( ) dépasse un seuil de protection ( ). Procédé selon les revendications 9 et 10, dans lequel, lors de l’étape (170) de déclenchement d’une alerte, une alerte est déclenchée en outre lorsque la valeur de décentralisation ( ) dépasse le seuil d’alerte en intégrité ( ). 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 contrôle d’intégrité (20) du recalage sur une pluralité d’amers (14-1, …, 14-N) lors de la géolocalisation d’un dispositif mobile (12) suite à une ou plusieurs défaillances de types prédéterminés, comprenant des moyens techniques (21, 22, 23) adaptés pour mettre en œuvre le procédé selon l’une quelconque des revendications 1 à 11.