FR 3037415 A1 20161216 FR 1555296 A 20150610 1 PROCÉDÉ DE TRAITEMENT DES MESURES D'AU MOINS UN CAPTEUR ELECTRONIQUE LOGE DANS UN APPAREIL PORTABLE [1] L'invention concerne un procédé de traitement des mesures d'au moins un capteur électronique logé dans un appareil portable par un utilisateur afin d'en déduire les valeurs d'une donnée à des instants successifs. L'invention a également pour objet un support d'enregistrement d'informations, une unité électronique et un appareil pour la mise en oeuvre de ce procédé. [2] De tels procédés connus de traitement comporte : - une phase active lors de laquelle, à chaque instant tk, un premier calculateur exécute un premier algorithme d'estimation d'une nouvelle valeur PIk de la donnée, cette exécution comprenant : - la réception d'une nouvelle mesure d'une grandeur physique, réalisée par le capteur électronique, cette nouvelle mesure étant plus récente que la mesure reçue lors de la précédente exécution du premier algorithme, puis - l'estimation de la nouvelle valeur PIk de la donnée à partir de la nouvelle mesure reçue, - à un instant t1, l'arrêt de la phase active et le passage dans une phase d'indisponibilité qui dure plusieurs instants tk et pendant laquelle l'exécution par le 20 premier calculateur du premier algorithme est inhibée, puis - après un instant t2, le retour à la phase active. [3] Typiquement, le premier calculateur est connu sous le terme anglais de « host processor » ou « application processor » ou « context hub ». [4] Le passage dans la phase d'indisponibilité peut correspondre à une mise en 25 veille du premier calculateur ou à une mise en attente de l'exécution du premier algorithme car le premier calculateur doit exécuter d'autres tâches plus prioritaires. Après la phase d'indisponibilité, le procédé revient à la phase active et le premier algorithme est de nouveau exécuté. Toutefois, si rien n'est fait, les nouvelles mesures de la grandeur physique réalisées pendant la phase d'indisponibilité sont perdues et il 30 est alors impossible d'estimer une valeur P''' de la donnée à un instant tw compris dans l'intervalle de temps ]t1 ; t2]. [5] Pour résoudre ce problème, les déposants ont envisagés d'enregistrer toutes les nouvelles mesures de la grandeur physique réalisées pendant la phase d'indisponibilité. Ensuite, à l'issue de la phase d'indisponibilité, le premier calculateur 35 exécute d'abord le premier algorithme sur les données enregistrées avant de traiter les nouvelles mesures reçues depuis le retour à la phase active. Toutefois, la quantité de mesures enregistrées pendant la phase d'indisponibilité peut être importante. Dans ce cas, le traitement par le premier calculateur des mesures enregistrées peut 3037415 2 retarder de façon significative le traitement des nouvelles mesures reçues après l'instant t2. Ceci n'est pas souhaitable. [6] Les déposants ont aussi envisagé de faire exécuter le premier algorithme par un second calculateur pendant la phase d'indisponibilité. Plus précisément, il a été 5 envisagé d'utiliser en tant que second calculateur, le calculateur connu sous le terme anglais de « sensor processor » ou « sensor hub ». En effet, un tel second processeur est généralement déjà présent dans tout appareil incorporant un capteur électronique. Il n'est donc pas nécessaire d'ajouter un nouveau calculateur, en plus de ceux déjà existants, dans l'appareil. Typiquement, ce second calculateur est utilisé 10 pour gérer le transfert des mesures entre le capteur électronique et le premier calculateur. Ce second calculateur est aussi souvent utilisé pour mettre en forme les nouvelles mesures ou filtrer les nouvelles mesures avant de les transmettre au premier calculateur. Toutefois, la puissance de calcul de ce second calculateur est souvent beaucoup moins importante que celle du premier calculateur. Ainsi, ce 15 second calculateur est souvent incapable d'exécuter assez rapidement le premier algorithme. [7] Dans ce contexte, l'invention vise à proposer un tel procédé de traitement des mesures d'un capteur électronique : - qui permet à la fois d'estimer la valeur P''' de la donnée à au moins un instant tw 20 compris dans l'intervalle ]ti ; t2] tout en retardant moins le redémarrage de la phase active après la phase d'indisponibilité, et - qui ne nécessite pas obligatoirement l'usage d'un second calculateur aussi puissant que le premier calculateur. [008] Elle a donc pour objet un procédé de traitement des mesures d'au moins un 25 capteur électronique, logé dans un appareil portable par un utilisateur, afin d'en déduire les valeurs d'une donnée à des instants successifs dans le temps, dans lequel : - pendant la phase d'indisponibilité, un second calculateur : - acquière les nouvelles mesures de la grandeur physique réalisées par le 30 capteur électronique pendant cette phase d'indisponibilité, - relève des caractéristiques du signal mesuré par le capteur électronique à partir de ces nouvelles mesures reçues pendant cette phase d'indisponibilité, et - enregistre ces caractéristiques, le nombre de caractéristiques relevées 35 pendant la phase d'indisponibilité étant strictement inférieur au nombre de nouvelles mesures acquises pendant cette phase d'indisponibilité, - à un instant t2, l'arrêt de la phase d'indisponibilité et le passage dans une phase de réveil actif lors de laquelle le premier calculateur exécute un second algorithme d'estimation d'une valeur P''' de la donnée à un instant tw compris dans l'intervalle ]t1 ; t] à 3037415 3 partir des caractéristiques relevées et enregistrées pendant la phase d'indisponibilité par le second calculateur, puis - l'arrêt de la phase de réveil actif et le retour à la phase active. [009] Dans le procédé revendiqué, le second calculateur est utilisé pour relever et 5 enregistrer des caractéristiques du signal mesuré qui sont moins nombreuses que les nouvelles mesures réalisées pendant la phase d'indisponibilité. Pour faire un tel traitement, le second calculateur n'a pas besoin d'être aussi puissant que le premier calculateur. Ensuite, avant de retourner à la phase active, le premier calculateur exécute, lors d'une phase de réveil actif, un second algorithme qui permet d'estimer la 10 valeur P''' de la donnée à un instant tw compris dans l'intervalle ]t1 ; t2] à partir des caractéristiques relevées par le second calculateur et sans utiliser directement les nouvelles mesures de la grandeur physique réalisée pendant la phase d'indisponibilité. Etant donné que les caractéristiques relevées sont moins nombreuses que les nouvelles mesures réalisées pendant la phase d'indisponibilité, 15 l'exécution de ce second algorithme est plus rapide que s'il était exécuté sur l'ensemble de ces nouvelles mesures. Le retard pour redémarrer la phase active est donc limité. [0010] Les modes de réalisation de ce procédé peuvent comporter une ou plusieurs des caractéristiques suivantes : - lors de la phase active, l'estimation de la nouvelle valeur Pik de la donnée est réalisée à partir de la nouvelle mesure reçue et également à partir de la précédente valeur Pik_i de la donnée estimée lors de la précédente exécution du premier algorithme à l'instant tk_i, et lors de la première exécution du premier algorithme après l'arrêt de la phase de réveil actif, la précédente valeur PIk_i de la donnée nécessaire pour estimer la nouvelle valeur Pik de la donnée est construite à partir de la valeur P''' fournie par l'exécution du second algorithme d'estimation pendant la phase de réveil actif ; - - le capteur électronique comporte une centrale inertielle embarquée dans l'appareil portable et la grandeur physique mesurée par cette centrale inertielle est représentative de la direction de déplacement de l'appareil et de l'amplitude de ce déplacement entre un instant précédent tk_i et l'instant tk, - la donnée est la localisation de l'appareil portable à l'intérieur d'un espace et chaque valeur de la donnée comporte au moins une position Pik possible de l'appareil à l'intérieur de cet espace, - le premier et le second algorithmes d'estimation sont des algorithmes d'estimation d'au moins une position possible Pik de l'appareil à l'intérieur de cet espace, - à chaque instant tk, l'exécution du premier algorithme, par le premier calculateur, génère au moins une position possible Pik à partir d'au moins une 3037415 4 position possible Pik-i estimée à l'aide de ce premier algorithme à l'instant précédent tk-1, - avant l'instant t1, au moins une dernière position possible Pit1 est construite à partir de la position possible Pik générée à l'instant tk précédent immédiatement 5 l'instant t1, puis enregistrée, - pendant la phase d'indisponibilité, les caractéristiques du signal mesuré par la centrale inertielle qui sont relevées comportent une suite temporellement ordonnée de caractéristiques d'un trajet parcouru par l'appareil pendant la phase d'indisponibilité, cette suite de caractéristiques constituant un relevé du 10 trajet parcouru par l'appareil pendant la phase d'indisponibilité, - pendant la phase de réveil actif, l'exécution du second algorithme génère au moins une position initiale Pit2, pour cela la position initiale Pit2 est estimée à partir de la dernière position construite Pit1 et du trajet parcouru relevé pendant la phase d'indisponibilité, et 15 - lors de la première exécution du premier algorithme après l'instant t2, la position précédente Pik-i est construite à partir de la position initiale Pit2; - - à chaque instant tk, l'exécution du premier algorithme, par le premier calculateur, génère plusieurs positions possibles Pik à partir de plusieurs positions possibles Pik-i estimées à l'aide de ce premier algorithme à l'instant 20 précédent tk-i, - avant l'instant t1, plusieurs positions possibles Pit1 différentes sont construites à partir des positions possibles Pik générées à l'instant tk précédent immédiatement l'instant t1, puis enregistrées, - pendant la phase de réveil actif, l'exécution du second algorithme génère 25 plusieurs positions initiales Pit2, pour cela chaque position initiale Pit2 est estimée à partir d'une dernière position construite Pit1 respective et du trajet parcouru relevé pendant la phase d'indisponibilité, et - lors de la première exécution du premier algorithme après l'instant t2, chaque position précédente Pik-i est construite à partir d'une position initiale possible 30 Pit2 respective ; - la construction desdites plusieurs dernières positions initiales Pit1 comporte : - la sélection de M positions les plus probables de l'appareil à partir de poids associés aux positions estimées Pik générées à l'instant tk précédent immédiatement l'instant t1, où M est un entier supérieur à deux et strictement 35 inférieur au nombre de positions possibles Pik générées par l'exécution du premier algorithme à l'instant tk précédent immédiatement l'instant t1, et chaque poids associé à une position possible Pik respective représente la probabilité que l'appareil se trouve dans cette position Pik, puis 3037415 5 - l'enregistrement en tant que dernières positions possibles P't1 des seules M positions les plus probables identifiées ; - l'exécution du second algorithme comporte : - la correction du trajet relevé pendant la phase d'indisponibilité à l'aide d'un 5 jeu de contraintes prédéfinies sur les déplacements de l'appareil dans l'espace, cette correction comportant des modifications des distances entre des points du trajet relevé et des modifications des amplitudes des changements de direction au niveau de ces points du trajet relevé pour obtenir un trajet corrigé qui n'enfreint aucune des contraintes prédéfinies, chaque 10 modification étant réalisée en tenant compte aussi bien des points situés en amont qu'en aval du point affecté par cette modification de la distance ou de l'amplitude du changement de direction, puis - la construction de l'estimation de la position initiale P't2 comporte : - le positionnement de l'origine du trajet corrigé en fonction de la dernière 15 position P't1 estimée par l'exécution du premier algorithme, puis - l'obtention de la position P't2 à partir de la position du point d'arrivée du trajet corrigé lorsque l'origine de ce trajet corrigé a été ainsi positionnée ; - - à chaque instant tk, le premier algorithme calcule également la valeur d'un facteur correctif d'un biais de mesure affectant l'un des capteurs embarqués à 20 l'intérieur de l'appareil, et - le procédé comporte la correction du trajet parcouru relevé pendant la phase d'indisponibilité pour obtenir un trajet corrigé, cette correction étant réalisée en appliquant aux mesures reçues de ce capteur embarqué le facteur correctif calculé à l'instant tk précédent immédiatement l'instant t1, puis 25 - la position P't2 est estimée en utilisant ce trajet corrigé ; - - le capteur électronique comporte une centrale inertielle embarquée dans l'appareil portable et la grandeur physique mesurée par cette centrale inertielle est représentative de la direction de déplacement de l'appareil et de l'amplitude de ce déplacement entre un instant précédent tk_i et l'instant tk, 30 - la donnée est la localisation de l'appareil portable à l'intérieur d'un espace et chaque valeur de la donnée comporte au moins une position PIk possible de l'appareil à l'intérieur de cet espace, - le premier et le second algorithmes d'estimation sont des algorithmes d'estimation d'au moins une position possible PIk de l'appareil à l'intérieur de 35 cet espace, - pendant la phase d'indisponibilité, les caractéristiques du signal mesuré par la centrale inertielle qui sont relevées comportent des points caractéristiques 3037415 6 d'un trajet parcouru par l'appareil pendant cette phase d'indisponibilité, le relevé des caractéristiques du signal mesuré comportant à cet effet: - l'identification, à partir des nouvelles mesures de la centrale inertielle reçues pendant la phase d'indisponibilité, de points caractéristiques du 5 trajet moins nombreux que les nouvelles mesures de la centrale inertielle reçues pendant la phase d'indisponibilité, et - pour chaque paire de points caractéristiques identifiés successifs le long du trajet, l'approximation du trajet de l'appareil entre ces deux points caractéristiques par un polynôme dont les coefficients sont déterminés ; 10 - l'identification de points caractéristiques du trajet comporte l'identification, en tant que point caractéristique du trajet, de chaque point où l'amplitude d'un changement de direction de déplacement de l'appareil dépasse un seuil prédéterminé ; - l'identification de points caractéristiques du trajet comporte également 15 l'identification d'un point caractéristique supplémentaire à chaque fois qu'au moins l'une des conditions suivantes est satisfaite : - la distance parcourue depuis le dernier point caractéristique identifié est supérieure à un seuil prédéterminé, - l'amplitude d'une variation de pression mesurée par un capteur de pression 20 embarqué à l'intérieur de l'appareil est supérieure à un seuil prédéterminé, - l'amplitude d'un déplacement de l'appareil par rapport au corps d'une personne qui transporte manuellement l'appareil est supérieure à un seuil prédéterminé, ou - le temps écoulé depuis le dernier point caractéristique identifié est supérieur 25 à un seuil prédéterminé ; - pendant toute la durée de la phase d'indisponibilité, le premier calculateur est mis en veille pour économiser de l'énergie ; - lorsqu'il est exécuté par le premier calculateur, le premier algorithme estime la position Pik de l'appareil en utilisant un jeu de contraintes prédéfinies sur les 30 déplacements de l'appareil dans l'espace ; - le premier algorithme est un filtre particulaire. [0011] Ces modes de réalisation présentent en outre les avantages suivants : - L'utilisation de la valeur P''' fournie par l'exécution du second algorithme d'estimation pour initialiser la valeur Pik-i du premier algorithme permet d'initialiser 35 plus précisément cette valeur Pik_i car il est alors tenu compte de ce qui s'est passé pendant la phase d'indisponibilité. - Après l'instant t2, initialiser la position Pik-i à partir de la dernière position construite P'ti avant le passage dans la phase d'indisponibilité et du trajet parcouru relevé 3037415 7 pendant la phase d'indisponibilité, permet d'initialiser plus précisément cette position P'k-1- - Générer plusieurs positions initiales possibles en sortie de la phase d'indisponibilité à partir de plusieurs positions respectives enregistrées au moment où le calculateur 5 bascule dans la phase d'indisponibilité permet de reprendre, à l'issue de la phase d'indisponibilité, les mêmes hypothèses sur la position de l'appareil que celles formulées juste avant le passage dans la phase d'indisponibilité. On assure ainsi la continuité de ces hypothèses malgré l'existence d'une phase d'indisponibilité. Autrement dit, on limite la perte d'informations liée à la phase d'indisponibilité et donc 10 on améliore la localisation de l'appareil au moment où l'on revient dans la phase active. - Enregistrer seulement les M positions les plus probables juste avant le basculement dans la phase d'indisponibilité permet de diminuer la quantité d'informations stockées et le nombre d'hypothèses à traiter sans pour autant décroître la précision de la 15 localisation de l'appareil au moment de la sortie de la phase d'indisponibilité. - L'estimation de la position initiale à la sortie de la phase d'indisponibilité en utilisant un second algorithme qui corrige la position de chaque point du trajet en prenant en compte tous les autres points du trajet aussi bien situés en amont qu'en aval de ce point, permet d'obtenir un trajet corrigé plus précis. 20 - Corriger le trajet relevé pendant la phase d'indisponibilité à l'aide du facteur correctif d'un biais de mesure estimé juste avant le basculement dans la phase d'indisponibilité permet d'obtenir un trajet corrigé plus précis. On obtient alors une estimation plus précise de la position de l'appareil en sortie de la phase d'indisponibilité. 25 - L'enregistrement de points caractéristiques du trajet à la place de toutes les positions de l'appareil, permet de limiter la quantité d'information enregistrée pendant cette phase d'indisponibilité et donc d'économiser de l'espace mémoire. De plus, le fait de limiter la quantité d'informations enregistrée permet d'accélérer l'exécution du second algorithme qui estime la position actuelle de l'appareil à la sortie de la phase 30 d'indisponibilité. - Identifier des points caractéristiques supplémentaires en plus de ceux où un changement de direction important se produit, permet d'améliorer la précision du trajet relevé. [0012] L'invention a également pour objet un support d'enregistrement d'informations 35 comportant des instructions pour exécuter le procédé ci-dessus lorsque ces instructions sont exécutées par un calculateur électronique. [0013] L'invention a également pour objet une unité électronique de traitement des mesures d'au moins un capteur électronique afin d'en déduire les valeurs d'une donnée à des instants successifs dans le temps, cette unité comportant au moins un 40 premier et un second calculateurs électroniques, dans lequel : 3037415 8 - le premier calculateur est programmé pour : - exécuter, pendant une phase active, à chaque instant tk, un premier algorithme d'estimation d'une nouvelle valeur PIk de la donnée, cette exécution comprenant : 5 - la réception d'une nouvelle mesure d'une grandeur physique, réalisée par le capteur électronique, cette nouvelle mesure étant plus récente que la mesure reçue lors de la précédente exécution du premier algorithme, puis - l'estimation de la nouvelle valeur PIk de la donnée à partir de la 10 nouvelle mesure reçue, - à un instant t1, arrêter de la phase active et passer dans une phase d'indisponibilité qui dure plusieurs instants tk et pendant laquelle l'exécution par le premier calculateur du premier algorithme est inhibée, - le second calculateur est programmé pour, pendant la phase d'indisponibilité : 15 - acquérir les nouvelles mesures de la grandeur physique réalisées par le capteur électronique pendant cette phase d'indisponibilité, - relever des caractéristiques du signal mesuré par le capteur électronique à partir de ces nouvelles mesures reçues pendant cette phase d'indisponibilité, et 20 - enregistrer ces caractéristiques, le nombre de caractéristiques relevées pendant la phase d'indisponibilité étant strictement inférieur au nombre de nouvelles mesures acquise pendant cette phase d'indisponibilité, et - le premier calculateur est programmé pour : - à un instant t2, arrêter la phase d'indisponibilité et passer dans une phase de 25 réveil actif lors de laquelle il exécute un second algorithme d'estimation d'une valeur P''' de la donnée à un instant tw compris dans l'intervalle ]t1 ; t] à partir des caractéristiques relevées et enregistrées pendant la phase d'indisponibilité par le second calculateur, puis - arrêter de la phase de réveil actif et retourner à la phase active. 30 [0014] Enfin, l'invention a également pour objet un appareil directement transportable à la main par un être humain, cet appareil comportant au moins un capteur électronique apte à mesurer une grandeur physique, dans lequel l'appareil comporte aussi l'unité électronique revendiquée, cette unité électronique étant raccordée à ce capteur électronique. 35 [0015] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif, et faite en se référant aux dessins sur lesquels : 3037415 9 - la figure 1 est une illustration schématique en coupe verticale d'un bâtiment à l'intérieur duquel est mis en oeuvre un procédé de localisation d'un appareil ; - la figure 2A est une illustration schématique d'un appareil de localisation ; - la figure 2B est une illustration schématique d'un mode de réalisation d'une unité de 5 calcul de l'appareil de la figure 2B ; - la figure 3 est une illustration schématique d'un jeu de contraintes utilisé pour localiser l'appareil de la figure 2A dans le bâtiment de la figure 1 ; - la figure 4 est un organigramme d'un procédé de localisation mis en oeuvre par l'appareil de la figure 2A; 10 - la figure 5 est une illustration schématique d'un exemple de graphe des trajets possibles utilisé par l'appareil de la figure 2A; - les figures 6 et 7 sont des exemples de trajets, respectivement, relevé et corrigé par l'appareil de la figure 2A; - la figure 8 est une illustration schématique d'une étape de correction d'un trajet 15 relevé à l'aide du graphe de la figure 5, - la figure 9 est un chronogramme de différentes phases de fonctionnement de différents calculateurs utilisés dans l'appareil de la figure 2A, et - la figure 10 est un organigramme d'un procédé de traitement des mesures d'au moins un capteur logé dans un appareil tel que l'appareil de la figure 2A. 20 [0016] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. Dans la suite de cette description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détail. [0017] L'exemple de mode de réalisation de l'invention est décrit pas suite dans le cas particulier d'un procédé de localisation d'un appareil à l'intérieur d'un bâtiment. 25 Toutefois, la solution décrite dans ce cas particulier s'applique à tout contexte où le problème exposé dans l'introduction de cette demande de brevet se rencontre. [0018] La figure 1 représente un ensemble comportant un bâtiment 2, à l'intérieur duquel un piéton 4 peut se déplacer librement en marchant. Le bâtiment 2 se divise en plusieurs étages. Ici, seul un rez-de-chaussée 6 et un premier étage 8 sont 30 représentés. Les étages sont reliés entre-eux par des zones de changement d'étage telles qu'un escalier ou un ascenseur. Chaque étage comporte des pièces et des couloirs délimités par des murs infranchissables par le piéton 4. Le piéton 4 peut rentrer à l'intérieur d'une pièce uniquement en franchissant une porte. Ici, l'intérieur d'une pièce peut également comporter des obstacles infranchissables par le piéton 4 35 comme, par exemple, des piliers ou autres éléments de construction du bâtiment 2. [0019] Pour aider le piéton 4 à se localiser à l'intérieur du bâtiment 2, celui-ci transporte, par exemple directement à la main, un appareil 10 de localisation. L'appareil 10 est capable de se localiser sur une carte du bâtiment 2 seulement à l'aide d'un capteur de position relative. En particulier, l'appareil 10 peut repérer sa position à l'intérieur du bâtiment 2 sans utiliser un système de navigation faisant 3037415 10 appel à des balises externes de repérage, par exemple des balises implantées dans l'environnement du bâtiment 2. Ces balises externes peuvent être des satellites ou des émetteurs d'ondes radio fixés au bâtiment 2. Par conséquent, l'appareil 10 peut repérer sa position sans utiliser un système GPS (« Global Positionning System »). 5 [0020] La figure 2A représente plus en détail l'appareil 10. L'appareil 10 comporte une unité électronique 11 de calcul qui capable, en outre, de localiser l'appareil 10. Cette unité 11 comporte une mémoire 12 ainsi qu'un calculateur électronique principal 14 et un calculateur électronique 15 dédié au traitement des mesures. Les calculateurs 14 et 15 sont des calculateurs programmables capables chacun 10 d'exécuter des instructions enregistrées dans une mémoire. Le calculateur 14 est connu sous le terme anglais de « host processor » ou « application processor ». Il est configuré pour exécuter les différentes opérations nécessaires au fonctionnement général et normal de l'appareil 10. A cet effet, le calculateur 14 est raccordé à une mémoire 12. La mémoire 12 est connue sous le terme anglais de « host memory ». 15 La mémoire 12 comporte les instructions nécessaires pour exécuter le procédé de la figure 4. De plus, la mémoire 12 comprend au moins : - une carte 16 du bâtiment 2, - des instructions pour exécuter un premier algorithme 18 d'estimation de la position de l'appareil 10 à l'intérieur du bâtiment 2, et 20 - des instructions pour exécuter un second algorithme 20 d'estimation de la position de l'appareil 10 à l'intérieur du bâtiment 2. [0021] L'unité 11 et la carte 16 sont décrits plus en détail en référence, respectivement, aux figures 2B et 3. [0022] L'algorithme 18 est un algorithme « en ligne » ou un « algorithme 25 d'apprentissage incrémental » plus connu sous le terme anglais de « online algorithm ». L'algorithme 18 est également désigné sous le terme « d'algorithme d'estimation temps réel » ou de « real-time algorithm » en anglais. Il s'agit d'un algorithme capable d'estimer, à chaque instant tk, la position actuelle PAk de l'appareil 10 à l'intérieur du bâtiment 2 en utilisant seulement les mesures acquises par une 30 centrale inertielle 22. Les instants tk sont les instants auxquels l'algorithme 18 est exécuté pour estimer la position actuelle PAk de l'appareil 10. Typiquement, deux instants tk_i et tk immédiatement successifs sont séparés l'un de l'autre par une durée supérieure à 0,01s ou 0,1s. Ici, l'algorithme 18 est un filtre particulaire, plus connu sous le terme anglais de « Particle filter ». Les filtres particulaires sont bien connus. 35 Par exemple, le lecteur peut se référer aux documents suivants : - les demandes de brevet W02012158441 et US8548738B1, - l'article O. Woodman et Al, « Pedestrian localisation for indoor environments », ACM, 2008. - la la thèse de J. Straub, « Pedestrian indoor localisation and tracking using a 40 particule filter combined with a leaming accessibility map », thèse, août 2010, 3037415 11 Université technique de Munich. Cette thèse est téléchargeable à l'adresse suivante : http://people.csail.m it. edu/jstraub/down load/Straub1OPedestrian Local ization. pdf. Par la suite, cette thèse est référencée sous le terme « Straub2010 ». [0023] Ainsi, le fonctionnement d'un filtre particulaire est considéré comme connu. 5 Par exemple, dans ce mode de réalisation, l'algorithme 18 est celui décrit dans la demande française déposée sous le numéro FR1455575 le 18 juin 2014 par le Commissariat à l'énergie atomique et aux énergies alternatives. [0024] On rappelle ici que dans cette demande FR1455575, les coordonnées de chaque particules Si sont mise à jour à l'aide d'une loi de déplacement. La loi de 10 déplacement permet de calculer, à partir de la direction ek et l'amplitude lk du déplacement de l'appareil 10, acquis à l'instant tk, le déplacement de la particule Si depuis sa position précédente Pik-i jusqu'à sa nouvelle position Pik. Ce déplacement est directement corrélé à celui de l'appareil 10. Typiquement, ce déplacement entre les positions Pik-i et Pik est identique ou très proche de celui de l'appareil 10 entre les 15 instants tk_i et tk. Par la suite, l'exposant « i » est l'identifiant de la particule et l'indice « k » est le numéro d'ordre de l'instant tk où est exécuté l'algorithme 18. [0025] Lorsque le piéton 4 se déplace en marchant sur le plancher de l'étage 8, une loi de déplacement est donnée par les relations suivantes : = vk * At * cosek ; 20 yik = vk * At *sinek, où : - (xik, yik) et (xik_i, Yu) sont les coordonnées, dans le plan du plancher, des positions Pi', et Pik_i de la particule S'; - At est l'intervalle de temps tk-tk_i, et 25 - vk est la vitesse du déplacement entre les instants tk et tk_i. Elle est donnée par la rapport I k/At. [0026] Dans le cas d'un procédé de localisation de l'appareil 10 mettant en oeuvre un filtre particulaire, il est intéressant d'explorer le plus grand nombre possible de trajectoires avec les particules. Ainsi, classiquement, le déplacement de chaque 30 particule est perturbé de façon aléatoire. Par exemple, pour cela, la loi de déplacement utilisable est la suivante : Xik = Xik-i+ vk *At * cosek+pix; yik= vk * At *sinek+ ; où pix et piy sont des variables aléatoires. 35 [0027] A chaque instant tk et pour chaque particule Si, les valeurs de ces variables pix et piy sont tirées au sort en fonction d'une loi prédéfinie de probabilité centrée, c'est-à-dire caractérisée par une espérance mathématique nulle. Ainsi, la moyenne des valeurs de chaque variable aléatoire pix et piy aux différents instants successifs tk tend vers zéro lorsque k augmente. Par exemple, cette loi prédéfinie de probabilité est la 40 même pour les variables aléatoires pix et piy et pour toutes les particules Si. Par la 3037415 12 suite, elle est notée Lpxy. Cette loi Lpxy est caractérisée par un écart type an, prédéterminé. Ici, l'écart type an, est constant et indépendant des mesures de la centrale inertielle 22 pour une mise à jour à chaque pas. Par exemple, l'écart type o-xy est supérieur à 5 cm ou 10 cm et, de préférence, inférieur à 35 cm. Par exemple, la 5 loi Lpxy est une distribution uniforme ou une distribution gaussienne. [0028] Dans la réalité, il peut aussi exister un biais de mesure, appelé biais de direction, sur la mesure de la direction ek. Un tel biais de direction peut provenir d'un défaut sur les capteurs de la centrale inertielle 22. Ce biais de direction peut aussi être causé par le fait que le piéton 4 ne pointe pas l'appareil 10 dans sa direction de 10 déplacement. De façon similaire, il peut aussi exister un biais de mesure, appelé ici biais de pas, sur la mesure de l'amplitude lk du déplacement de l'appareil 10. Ce biais de pas peut provenir d'un défaut des capteurs de la centrale inertielle 22 ou d'une erreur de modélisation. Typiquement, ces biais sont constants au moins pendant un intervalle de temps assez long pour pouvoir les estimer et les corriger comme décrit 15 par la suite. [0029] Ici, pour compenser et corriger ces biais de direction et de pas, la loi de déplacement utilisée intègre des facteurs correctifs, respectivement a' et ci associés à chaque particule Si. Par exemple, la loi de déplacement est donnée par les relations suivantes : 20 xik = + vk * At * (1-Fcik) * cos(ek+ aik) + ; y'k= + vk *At * (1-Fcik) * sin(ek+ aik) + Eik = £'k-1 PiE ; ak = ak_l + pa ; où : 25 - Eik et Eik-i sont les valeurs, respectivement aux instants tk et tk_i, du facteur correctif c' utilisé pour corriger le biais de pas, et - ail, et aik_i sont les valeurs, respectivement aux instants tk et tk_i, du facteur correctif a' utilisé pour corriger le biais de direction, - pi, et pi, sont des variables aléatoires. 30 [0030] Les variables aléatoires pi, et pi, sont utilisées pour les mêmes raisons et de la même manière que les variables pix et piy précédemment introduites. Ainsi, une nouvelle valeur des variables pi, et pi, est tirée au sort à chaque nouvel instant tk et pour chaque particule Si en fonction, respectivement, d'une loi prédéfinie de probabilité Lp, et d'une loi prédéfinie de probabilité Lpa. Typiquement, ces lois Lp, et 35 Lpa sont les mêmes pour toutes les particules Si. Ici, les espérances mathématiques des lois Lp, et Lpa sont égales à zéro. Par conséquent, comme pour les variables aléatoires pix et piy, la moyenne des valeurs de chaque variable aléatoire pi, et pi, aux différents instants successifs tk tend vers zéro lorsque k augmente. [0031] De plus, la fonction des variables pi, et pi, est seulement de venir légèrement 40 perturber les précédentes valeurs Eik_i et aik_1 des facteurs correctifs ci et a' pour que 3037415 13 les valeurs des facteurs correctifs c' et a' restent stables dans le temps. A cet effet, les écarts types a, et aa, respectivement, des lois Lp, et Lpa, ne permettent pas une variation rapide des valeurs des facteurs correctifs c' et a'. Ici, a cet effet, l'écart type a, est choisi suffisamment petit pour que le ratio Zo-Ek /T soit inférieur à 10 %/s et, de 5 préférence, inférieur à 5 %/s ou 1%/s, où : - o-Ek est l'écart type de la loi Lp, lors de la k-ième itération de l'étape 126 (figure 4), - Zo-Ek est la somme des écarts types o-Ek entre la q-ième itération et la p-ième itération de l'étape 126, où q est un entier strictement inférieur à p, - T est la durée en seconde de l'intervalle de temps qui s'est écoulé entre la q-ième et 10 la p-ième itération de l'étape 126. [0032] Ici, l'écart type a, est constant. Ainsi, le précédent ratio peut aussi s'écrire : (p-q)o-E/T. Dans ce cas, quels que soient p et q, le ratio est constant. La différence p-q est généralement assez grande pour couvrir une période de temps supérieure à 1s ou 4s et, généralement, inférieure à 10 min ou 5 min ou 1 min. Par exemple, cette 15 différence entre p et q est constante quel que soit p. [0033] De façon similaire, l'écart type a, est choisi suffisamment petit pour que le ratio Zo-ak/T soit inférieur à 10°/s et, de préférence, inférieur à 5°/s ou 1°/s, où : - o-ak est l'écart type de la loi Lpa lors de la k-ième itération de l'étape 126, - Zo-ak est la somme des écarts types o-ak entre la q-ième itération et la p-ième itération 20 de l'étape 126, où q est un entier strictement inférieur à p, - T est la durée en seconde de l'intervalle de temps qui s'est écoulé entre la q-ième et la p-ième itération de l'étape 126. [0034] Ici, l'écart type a, est aussi constant. Ainsi, le précédent ratio peut aussi s'écrire : (p-q)o-a/T. 25 [0035] Comme pour les variable p'x et p'y, les variables p', et p', permettent d'explorer un grand nombre de valeurs possibles pour les facteurs correctifs c' et a'. [0036] Cette loi de déplacement fonctionne particulièrement bien dans les situations où le piéton 4 marche sur un sol horizontal. Par contre, la gestion des changements d'étages est réalisée, par exemple, comme décrit dans Straub2010. 30 [0037] L'algorithme 20 est un algorithme hors-ligne plus connu sous le terme anglais de « offline algorithm ». Contrairement à un algorithme temps réel, l'algorithme hors-ligne utilise toutes les mesures préalablement acquises entre deux instants t1 et t2 pour estimer la position PA, à plusieurs intants tw compris entre t1 et t2. Par exemple, l'algorithme 20 est un algorithme de « graph-matching » tel que celui décrit dans la 35 thèse suivante : Ivan Spassov : « Algorithm for Map-Aided Autonomous Indoor Pedestrian Positionning and Navigation », THÈSE N° 3961 (2007), présentée le 23 novembre 2007 à la faculté de l'environnement naturel, architectural et construit, Laboratoire de topométrie, ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE. Par la suite, cette thèse est référencée sous le terme « Spassov2007 ». 3037415 14 [0038] Le calculateur 14 est généralement un calculateur plus puissant que le calculateur 15, c'est-à-dire capable d'exécuter un plus grand nombre d'instructions par seconde. Ce calculateur 14 est également apte à basculer entre une phase active et une phase d'indisponibilité. Pendant la phase active, il a accès à un très grand 5 nombre d'informations dont notamment la carte 16 et il exécute des applications de haut niveau tel que les algorithmes 18 et 20. Au contraire, dans la phase d'indisponibilité, le calculateur 14 est incapable d'exécuter les algorithmes 18 et 20. La suite de la description de ce mode de réalisation est faite dans le cas particulier où la phase d'indisponibilité correspond à une phase où le calculateur 14 est mis en 10 veille pour économiser de l'énergie. Dans ce cas particulier, typiquement, le calculateur 14 bascule de la phase active vers la phase d'indisponibilité automatiquement en absence d'interaction avec le piéton 4 pendant un délai prédéterminé. Le piéton 4 interagit avec l'appareil 10 par l'intermédiaire d'une interface homme-machine. Le calculateur 14 peut aussi être mis en veille en réponse 15 à une commande transmise par le piéton 4. [0039] Le calculateur 15 est programmé pour gérer, en outre, le fonctionnement et l'acquisition des mesures de la centrale inertielle 22. Le calculateur 15 peut également effectuer des traitements limités sur les mesures acquises. Par contre, il est incapable d'exécuter les algorithmes 18 et 20 car ceux-ci nécessitent une 20 puissance de calcul dont le calculateur 15 ne dispose pas. Par contre, le calculateur 15 transmet les mesures, éventuellement traitées localement, au calculateur 14 si celui-ci est dans sa phase active. A cause de sa puissance de calcul plus limitée, le calculateur 15 consomme beaucoup moins d'énergie que le calculateur 14 pour fonctionner. Ainsi, ici, le calculateur 15 est en permanence alimenté tant que l'appareil 25 10 est allumé. Par exemple, il fonctionne en permanence tant qu'il est nécessaire de localiser l'appareil 10 dans le bâtiment 2. Le calculateur 15 a aussi accès à moins d'informations que le calculateur 14. Par exemple, le calculateur 15 ne peut pas accéder à la carte 16 enregistrée dans la mémoire 12. [0040] L'appareil 10 comporte la centrale inertielle 22. La centrale inertielle 22 logée 30 à l'intérieur de l'appareil 10 mesure la direction de déplacement et l'amplitude du déplacement de cette appareil dans une direction depuis une position précédente de ce même appareil. Une centrale inertielle est un capteur de position relative car elle permet seulement d'estimer une nouvelle position de l'appareil par rapport à sa précédente position. Ainsi, pour estimer correctement la position de l'appareil à un 35 instant tk, il faut absolument connaître la précédente position occupée à l'instant tk_i. A l'inverse, un capteur de position, tel qu'un capteur GPS (Global Positionning System), est un capteur de position absolue car il est capable d'estimer la position de l'appareil sans avoir aucune connaissance sur la position précédente de ce même appareil. Ici, la centrale inertielle 22 comporte un accéléromètre 24, un gyromètre 25 et un 40 magnétomètre 26. Les capteurs peuvent être des capteurs triaxes. De plus, dans ce 3037415 15 mode de réalisation, la centrale inertielle 22 comporte aussi un baromètre 27 pour mesurer l'altitude de l'appareil 10. [0041] La localisation de l'appareil 10 à partir des mesures de la centrale inertielle pose des problèmes particuliers qu'on ne rencontre pas avec l'utilisation de capteurs 5 de position absolue. En particulier, il est généralement complexe d'utiliser directement et seulement les mesures d'un capteur de position relative pour localiser l'appareil car, dans ce cas, les erreurs sur la localisation de l'appareil s'accumulent au cours du temps. La localisation ainsi obtenue devient alors rapidement inexploitable. Il faut donc être capable de régulièrement corriger cette localisation sans pour autant avoir 10 recours à un capteur de position absolue. Différentes solutions ont été proposées. Par exemple, pour corriger et améliorer l'estimation de la position de l'appareil, il est connu d'exploiter le fait qu'il existe des contraintes prédéfinies sur les déplacements de l'appareil à l'intérieur du bâtiment 2. Par exemple, une contrainte typique est qu'un déplacement ne peut pas traverser un mur. Ici, c'est la carte 16 qui contient ces 15 contraintes prédéfinies sur les déplacements de l'appareil 10 à l'intérieur du bâtiment 2. [0042] La centrale inertielle 22 transmet au calculateur 15, par l'intermédiaire d'un bus 19 de transmission d'informations, des mesures représentatives de la direction dans laquelle l'appareil 10 se déplace et de l'amplitude du déplacement dans cette 20 direction depuis la dernière position relevée pour cet appareil 10. [0043] L'appareil 10 est équipé d'une interface homme-machine comportant notamment un écran 28 permettant d'afficher une représentation graphique 30 de la carte 16 et, sur cette représentation graphique, un point PAk représentant la position actuelle de l'appareil 4 à l'intérieur du bâtiment 2. Ce point PAk est donc situé dans la 25 représentation graphique 30 à l'emplacement de la carte 16 correspondant à la position actuelle de l'appareil 10 et donc du piéton 4. [0044] Par exemple, dans ce mode de réalisation, l'appareil 10 est un téléphone ou une tablette électronique programmée pour exécuter le procédé de la figure 4. [0045] La figure 2B représente plus en détail un mode de réalisation possible de 30 l'unité 11 de calcul. L'unité 11 comprend une unité de traitement de déplacement désignée ici sous l'acronyme MPU (Motion Processing Unit) 40. Le MPU 40 comprend des capteurs internes. Ici, ces capteurs internes correspondent à ceux de la centrale inertielle 22. Un capteur interne est typiquement un MEMS (MicroElecroMechanical System) intégré sur la puce qui forme le MPU 40 ou intégré 35 dans le même enrobage isolant que celui qui enrobe les autres éléments du MPU 40. [0046] Le calculateur 14, la mémoire 12 et le MPU 40 sont raccordés les uns aux autres par le bus 19. Par exemple, le bus 19 est un bus PCIe (Peripherical Component Interconnect express), un bus USB (Univervsal Serial Bus), un bus série UART (Universal Asynchronous Receiver/Transmitter (UART), une interface AMBA 40 (Advanced Microcontroller Bus Architecture), un bus I2C (Inter-Integrated Circuit), un 3037415 16 bus SDIO (Serial Digital Input Output) ou autre bus similaire. De plus des signaux de signalisation additionnelle peuvent être utilisés en plus du bus 19 tel qu'une ligne d'interruption ou similaire. Des détails supplémentaires sur la réalisation du calculateur 14 et de la mémoire 12 peuvent être trouvés dans la demande déposée le 5 6 juin 2007 sous le numéro US11/774488, et dans la demande US12/106,921, déposée le 11 avril 2008. Des détails supplémentaires sur la réalisation de l'unité 40 sont disponibles auprès de la société InvenSense Inc en Californie. [0047] Dans ce mode de réalisation, le MPU 40 comprend, en plus des capteurs 24 à 27, le calculateur 15 et une mémoire interne 44. La mémoire 44 stocke des 10 instructions exécutables par le calculateur 15 pour la mise en oeuvre du procédé de la figure 4. La mémoire 44 est aussi utilisable pour enregistrer les données mesurées par les capteurs 24 à 27. Par exemple, la mémoire 44 comporte une mémoire tampon FIFO (first in, first out) pour enregistrer un jeu des mesures réalisées par les capteurs 24 à 27. Le MPU 40 comprend aussi un bus 46 qui raccorde entre eux le calculateur 15 15, les capteurs 24 à 27 et la mémoire 44. Le bus 46 est par exemple similaire ou identique au bus 19. La mémoire 44 stocke, par exemple, un seul échantillon des mesures au même instant de mesure. Cette échantillon est ensuite récupérer par le calculateur 14 avant d'être remplacé par une nouvelle échantillon contenant des mesures plus récentes. Pour cela, par exemple, le MPU 40 signale au calculateur 14 20 qu'un nouvel échantillon de mesures est disponible en envoyant au calculateur 14 un signal d'interruption. Dans ce mode de réalisation, la fréquence à laquelle le signal d'interruption est transmis au calculateur 14 peut dépendre de la période d'échantillonnage des mesures. Ainsi, un signal d'interruption est transmis au calculateur 14 après l'enregistrement de chaque nouvel échantillon dans la mémoire 25 44 ou après l'enregistrement d'un nombre prédéterminé supérieur à deux de nouveaux échantillons dans la mémoire 44. Dans ce dernier cas, la mémoire 44 peut stocker simultanément plusieurs échantillons successifs de mesures. Le mode de fonctionnement qui vient d'être décrit dans le cas particulier des capteurs internes 24 à 27 peut aussi s'appliquer aux mesures de capteurs externes 48 raccordés par un 30 bus 50 au MPU 40. Les capteurs externes 48 peuvent être situés et intégrés à l'intérieur de l'unité 11 ou situé sur un autre appareil distant. L'appareil distant peut être un bracelet, une montre ou un autre appareil distant portable par l'utilisateur. Cet appareil distant est raccordé à l'appareil 10 par une liaison filaire ou une liaison sans-fil. Le nombre de capteurs externes 48 peut être quelconque. Il peut s'agir d'un 35 accéléromètre, d'un gyroscope, d'un magnétomètre, d'un baromètre, d'un hydromètre, d'un thermomètre, d'un microphone, d'un capteur de proximité ou d'un capteur d'intensité lumineuse ou autres. [0048] La figure 3 représente graphiquement un exemple de contenu de la carte 16 pour l'étage 8 du bâtiment 2. Ce qui va maintenant être décrit pour l'étage 8 du 40 bâtiment 2 s'applique à chaque étage de ce bâtiment et au rez-de-chaussée 6. 3037415 17 [0049] Le plan de l'étage 8 est typiquement horizontal. La carte 16 est similaire à celle décrite dans la demande française déposée sous le numéro FR1455575. Ainsi, seuls les détails nécessaires à la compréhension de l'invention sont donnés. [0050] Le repère XYZ est un repère orthogonal dans lequel les directions X et Y sont 5 horizontales et la direction Z est verticale. [0051] Sur la figure 3, la périphérie de chaque pièce est délimitée par des murs représentés en traits fins. Ces murs sont des obstacles infranchissables par le piéton 4. De plus, chaque pièce comporte au moins une ouverture pour accéder à l'intérieur de cette pièce. Sur la figure 3, les ouvertures sont situées entre les extrémités des 10 murs marquées par des points. Une ouverture est typiquement une porte. Une pièce peut aussi contenir d'autres obstacles infranchissables par le piéton 4. Par exemple, un obstacle infranchissable est une cloison intérieure ou un pilier ou tout autre élément du bâtiment 2 qui ne peut pas être traversé par le piéton 4. Chaque objet infranchissable constitue une contrainte sur le déplacement de l'appareil 10 à 15 l'intérieur du bâtiment 2. [0052] Dans cette carte 16, la position et les dimensions de chaque objet infranchissable sont codées par un segment horizontal contenu dans le plan du plancher. Ainsi, chaque identifiant d'obstacle est associé à une paire de points Ejd et EJf. Les points Ejd et EJf marquent respectivement, le début et la fin du segment [E,d ; 20 EJf], où j est l'identifiant de l'obstacle infranchissable. Les coordonnées des points Ejd et E,f, dans le plan du plancher, sont connues et contenues dans la carte 16. Sur la figure 3, quelques points Ea à ER sont représentés. [0053] Le fonctionnement de l'appareil 10 va maintenant être décrit en référence au procédé de la figure 4 et à l'aide des figures 5 à 8. 25 [0054] A un instant to, lors d'une étape 100, le piéton 4 déclenche manuellement l'exécution par l'unité 11 du procédé de localisation de l'appareil 10 à l'intérieur du bâtiment 2. Le calculateur 14 est alors dans sa phase active ou bascule dans sa phase active. [0055] Ensuite, lors d'une étape 102, la centrale inertielle 22 mesure, à une 30 fréquence plus élevée que les instants tk, les grandeurs physiques telles que l'accélération, la direction du champ magnétique, la vitesse de rotation de l'appareil 10 et la pression atmosphérique. [0056] Lors d'une étape 104, le calculateur 15 acquiert, à cette fréquence plus élevée que celle des instants tk, les mesures de la centrale inertielle 22. Puis, dans ce 35 mode de réalisation, il exécute un algorithme de détection d'un pas du piéton 4 à partir des mesures de l'accéléromètre. En effet, à chaque fois que le piéton 4 pose un pied par terre, cela correspond à une évolution temporelle caractéristique de l'accélération mesurée. Le calculateur 15 utilise cette évolution temporelle caractéristique pour détecter l'instant tp auquel le piéton 4 pose un pied par terre. Par 40 exemple, le lecteur peut se référer à la demande de brevet déposée sous le numéro 3037415 18 FR14050950, le 5 juillet 2014 par le déposant et intitulée « Procédé et dispositif de comptage de pas ». [0057] Lors de l'étape 108, le calculateur 15 établit les mesures de la direction ek et de l'amplitude lk du déplacement de l'appareil 10 à l'instant tp puis les enregistre dans 5 la mémoire 44. Par la suite, on considère que la différence entre l'instant tp où les nouvelles mesures sont acquises et enregistrées dans la mémoire 44 et l'instant tk où l'algorithme 18 est exécuté par le calculateur 14 est négligeable. Dès lors, les instants tp et tk sont considérés comme étant égaux et seule la notation « tk » est aussi utilisée pour désigner l'instant où un pas est détecté. 10 [0058] Lors de l'étape 108, le calculateur peut établir également d'autres grandeurs physiques supplémentaires qui sont ensuite exploitées par les algorithmes 18 et 20 pour accroître la précision de la localisation de l'appareil 10. Par exemple, le calculateur 15 établit la valeur d'un indicateur Manip(tk). L'indicateur Manip(tk) est égal à « 1 » si l'appareil 10 a été déplacé par rapport au corps du piéton 4 à l'instant tk, et 15 égal à « 0 » dans le cas contraire. Pour cela, le lecteur peut se référer à la demande de brevet déposée sous le numéro FR1363723 le 31 décembre 2013 par le déposant et intitulée « Procédé et dispositif de détection de manipulation d'un appareil portable ». [0059] Toujours lors de l'étape 108, le calculateur 15 établit aussi une grandeur 20 physique P(tk) dont la valeur est égale à la différence entre les pressions mesurées aux instants tk-1 et tk. [0060] Ensuite, lors d'une étape 110, le calculateur 15 vérifie si le calculateur 14 est dans sa phase active. Dans l'affirmative, il rend disponible pour le calculateur 14 les mesures réalisées puis il retourne ensuite directement à l'étape 102 sans exécuter 25 l'étape 112 suivante. Par exemple, les mesures préparées sont enregistrées dans la mémoire 44 et un signal d'interruption est transmis au calculateur 14 pour lui indiquer qu'il peut récupérer les nouvelles mesures enregistrées dans la mémoire 44. Dans la négative, c'est-à-dire si le calculateur 14 est dans sa phase d'indisponibilité, alors, avant de retourner à l'étape 102, le calculateur 15 procède à l'étape 112. 30 [0061] Lors de l'étape 112, le calculateur 15 relève des caractéristiques du trajet parcouru par l'appareil 10 pendant que le calculateur 14 est dans sa phase d'indisponibilité. Par exemple, lors de cette étape 112, à chaque instant tk, le calculateur 15 construit la position PAk de l'appareil 10 uniquement à partir de la direction ek et l'amplitude lk mesurées et de la dernière position PAk_i construite. Dans 35 ce mode de réalisation, c'est cette suite de positions qui constitue le trajet relevé par le calculateur 15 pendant la phase d'indisponibilité du calculateur 14. Le trajet relevé est enregistré dans la mémoire 44. [0062] Lors de l'étape 112, le calculateur 15 cherche uniquement à relever le trajet parcouru par l'appareil 10 et non pas à le corriger. En effet, seul le calculateur 14 a 40 accès à la carte 16 pour corriger la position relevée de l'appareil 10. Ainsi, lors de 3037415 19 cette étape 112, le calculateur 15 ne corrige pas les positions construites pour éviter l'accumulation des erreurs de mesure de la centrale inertielle 22. En particulier, le calculateur 15 n'utilise pas la carte 16 ou tout autre source d'informations stockant des contraintes sur les déplacements du piéton 4 dans le bâtiment 2 pour éviter 5 l'accumulation de telles erreurs. Le trajet relevé lorsque le calculateur 14 est dans sa phase d'indisponibilité peut donc comporter des erreurs importantes. En particulier, les biais de mesure des capteurs de la centrale inertielle 22 s'accumulent au fur et à mesure que le temps s'écoule. [0063] Par contre, ici, lors de l'étape 112, le calculateur 15 procède à une opération 10 114 de relevé de caractéristiques du trajet parcouru pendant la phase d'indisponibilité. Cette opération 114 vise à limiter le nombre d'informations à traiter par le calculateur 14, lors de son retour à la phase active, pour obtenir une estimation de la position de l'appareil 10 à partir des mesures réalisées pendant la phase d'indisponibilité. Pour cela, le calculateur 15 identifie des points caractéristiques dans 15 le trajet parcouru par l'appareil 10 et enregistre ces points caractéristiques dans la mémoire 44. De plus, pour chaque paires de points caractéristiques successifs le long du trajet parcouru par l'appareil 10, le calculateur 15 détermine les coefficients d'un polynôme qui approxime le trajet de l'appareil 10 entre ces deux points caractéristiques. Par exemple, ici, ce polynôme est un polynôme de degré un. Ainsi, 20 le trajet de l'appareil 10 entre deux points caractéristiques immédiatement successifs est approximé par un segment rectiligne qui relie directement ces deux points caractéristiques. Les positions intermédiaires de l'appareil 10 entre deux points caractéristiques relevées par le calculateur 15 sont, par exemple, effacées. Pour obtenir moins de points caractéristiques qu'il y a de positions PAk construites, les 25 points caractéristiques sont, par exemple, les points du trajet où l'amplitude du changement de la direction de déplacement de l'appareil 10 est supérieure à un seuil Se prédéterminé. Par exemple, un changement de direction est identifié quand lek - Capl > Se, où Cap est la direction de déplacement de l'appareil 10 obtenue à partir des précédentes mesures. Le seuil Se est par exemple supérieur à 500 ou 80°. Ici, le 30 seuil Se est pris égal à 85°. [0064] Un exemple de méthode possible pour déterminer ces points caractéristiques où la direction de déplacement change beaucoup, est décrit en détail dans la thèse Spassov2007 au chapitre 3.2.1. La description de cette méthode n'est donc pas reprise ici. 35 [0065] De plus, dans ce mode de réalisation, à ces points caractéristiques correspondant à des changements importants dans la direction de déplacement, le calculateur 15 ajoute d'autres points caractéristiques supplémentaires. Ces points caractéristiques supplémentaires ne correspondent pas nécessairement à un changement de direction important, mais à un changement important dans une autre 40 grandeur physique mesurée par la centrale inertielle 22. Par exemple, le calculateur 3037415 20 15 ajoute un point caractéristique supplémentaire à chaque fois qu'au moins l'une des conditions suivantes est satisfaite : 1) La distance AD parcourue depuis le dernier point caractéristique identifié est supérieure à un seuil SD prédéterminé, 5 2) La différence de pression P(tk) est supérieure à un seuil prédéterminé Sp, 3) L'indicateur Manip(tk) vient de changer de valeur, et 4) Si la différence entre tk - T est supérieure à un seuil prédéterminé ST, où T est l'instant où le dernier point caractéristique a été identifié. [0066] Par exemple, le seuil SD est supérieur ou égal à 50m ou 100m, et inférieur ou 10 égal à 500m ou 250m. Ce seuil SD peut être ajusté en fonction des caractéristiques du bâtiment 2. Par exemple, si le bâtiment 2 comporte de larges espaces ouverts, le seuil SD sera plus grand que si le bâtiment 2 comporte seulement des petites pièces. Ici, le seuil SD est pris égal à 100m. Le seuil Sp est supérieur à 0,05 hPa et, généralement, supérieur à 0,1 hPa. Ici, le seuil Sp est pris égal à 0,1 hPa. Le seuil ST 15 est supérieur à 5s ou 10s. Ici, le seuil ST est égal à 100s. [0067] La position de chaque point caractéristique et de chaque point caractéristique supplémentaire est enregistrée, par le calculateur 15, dans la mémoire 44. [0068] En parallèle, immédiatement après l'instant to, le calculateur 14 s'apprête à exécuter pour la première fois l'algorithme 18. Il n'existe donc aucune position 20 précédemment estimée pour l'appareil 10 à l'aide de cet algorithme 18. Ainsi, le calculateur 14 commence par exécuter une étape 120 d'initialisation de l'algorithme 18. Pour cela, le calculateur 14 génère alors un jeu initial de No particules S'. Le nombre No de particules S' dépend notamment des connaissances initiales que l'on a sur la position de l'appareil 10 dans le bâtiment 2 et de la surface de ce bâtiment. 25 Typiquement, No est supérieur à 10 ou 100. No est également généralement inférieur à 5000 ou 1000. Chaque particule S' est associée initialement : - à une position initiale P'0 à l'intérieur du bâtiment 2, cette position P'0 comportant les coordonnées x'0, y'o et z'o de la particule S' dans le repère XYZ, - à une valeur initiale w'o du poids w' représentant la probabilité que l'appareil 10 se 30 trouve à l'emplacement de cette particule S' à l'instant to, - à une valeur initiale a'0 pour le facteur correctif a', et - à une valeur initiale Elo pour le facteur correctif c'. [0069] Lors de l'étape 120, la position initiale P'0 et les valeurs initiales w'o, a'0 et Elo sont initialisées. De nombreuses méthodes pour initialiser les positions P'0 et les 35 valeurs w'o de chaque particule sont connues. Par exemple, si la position initiale de l'appareil 10 est connue à plus ou moins 1 m près, les valeurs P'0 de toutes les particules S' sont tirées au hasard à l'intérieur d'un cercle centré sur la position connue et de rayon égale à 1 m. Si la position initiale de l'appareil 10 est totalement inconnue, alors les positions P'0 sont, par exemple, uniformément réparties dans 3037415 21 l'ensemble du bâtiment 2. A titre d'illustration, chaque valeur w'o est prise égale à 1/N0, où No est le nombre initial de particules générées. [0070] Chaque valeur a'o est tirée au hasard de manière à ce que la distribution des valeurs initiales a'o suive une loi Lpao de probabilité prédéterminée telle qu'une 5 distribution uniforme ou gaussienne ou autre. La loi Lpao n'est généralement pas la même que la loi Lpa utilisée pour obtenir les valeurs de la variable aléatoire p',. C'est la connaissance a priori que l'on a sur la distribution du biais de direction qui permet de choisir la loi de probabilité des valeurs initiales a'o qui ressemble le plus à celle observée dans la réalité. Par exemple, c'est aussi cette connaissance à priori sur la 10 répartition des biais de direction qui permet de fixer la valeur de l'écart type cyao de la loi Lpao. Par exemple, l'écart type o-ao est choisi égale à 360° si l'on a aucune information sur le biais de direction. Dans un autre exemple, l'écart type o-ao est choisi inférieur à 45° si l'on a quelques informations sur le biais de direction. Contrairement au cas des variables aléatoires précédentes, cette loi de probabilité Lpao n'a pas 15 nécessairement une moyenne nulle. [0071] Chaque valeur initiale E'o est choisie comme décrit précédemment pour les valeurs initiales a'o sauf que l'on utilise une loi prédéfinie LpEo de probabilité pour le biais de pas à la place de la loi Lpao. De plus, la loi de probabilité LpEo n'est pas nécessairement identique à la loi Lpao. En effet, généralement, les biais de direction et 20 de pas ne sont pas corrélés. Typiquement, l'écart type o-E0 de la loi LpEo est égale 30 % à plus ou moins 5 % près si l'on a aucune information sur le biais de pas. Dans un autre exemple, l'écart type o-E0 est choisi inférieur à 20 % si l'on a quelques informations sur le biais de pas. [0072] Ensuite, le calculateur 14 exécute une phase 122 de localisation de l'appareil 25 10 en exécutant l'algorithme 18. [0073] Plus précisément, lors d'une étape 124, le calculateur 14 reçoit la direction ek et l'amplitude lk précédemment établies par le calculateur 15 pour l'instant tk. A partir de ces mesures, le calculateur 14 calcule la vitesse vk du déplacement actuel de l'appareil 10 depuis sa dernière position. Ensuite, le calculateur 14 estime les 30 nouvelles positions Pik en exécutant les étapes suivantes. [0074] Lors d'une étape 126, chaque particule est déplacée de manière corrélée avec le déplacement mesuré de l'appareil 10 depuis sa précédente position Pik-i jusqu'à une nouvelle position Pik. Pour cela, les coordonnées de chaque particule S' sont mises à jour en fonction : 35 - de la direction ek et l'amplitude lk reçues, - des coordonnées de la précédente position Pik-i de cette particule, et - de la loi de déplacement utilisée. [0075] A chaque nouvelle exécution de l'étape 126, les nouvelles valeurs pour les variables p'x, p'y, p', et p', sont tirées au sort à l'aide des lois, respectivement, Lpxy, Lpa 40 et LpE. 3037415 22 [0076] Lors d'une étape 128, le calculateur 14 met à jour les poids w' de chaque particule S'. Plus précisément, le calculateur 14 diminue le poids w' de la particule S' si son dernier déplacement de la position PIk-i à la position P'k a enfreint des contraintes prédéfinies associées à la carte 16. 5 [0077] Typiquement, dans cette exemple de mode de réalisation, pour chaque particule S', le calculateur 14 vérifie la contrainte suivante : - contrainte 1) : le segment [P'k_i ; P'k] ne doit pas traverser un obstacle infranchissable de la carte 16 [0078] De façon générale, ici on définit une contrainte sur le déplacement de la 10 particule S' comme étant une condition qui, si elle est vérifiée par la particule S', est utilisée pour augmenter le poids w' de cette particule S' par rapport au poids des particules qui ne vérifient pas cette condition. A l'inverse si cette condition n'est pas vérifiée par la particule S', alors elle est utilisée pour diminuer le poids w' de cette particule S' par rapport aux poids des particules qui vérifient cette condition. 15 [0079] Pour évaluer la contrainte 1), pour chaque particule S', le calculateur 14 recherche s'il existe une intersection entre le segment [P'k_i ; P'k] et chaque obstacle infranchissable de la carte 16. Ici, puisque chaque obstacle est codé par un segment, cette recherche d'intersection revient à rechercher une intersection entre deux segments. 20 [0080] S'il existe une intersection, alors, une valeur très faible ou nulle est affectée au poids w'. Une valeur très faible est une valeur inférieure à 0,2 ou 0,1 par exemple. Dans le cas contraire, la valeur du poids w' reste inchangée. [0081] Lors d'une étape 130, le calculateur 14 procède à la normalisation des poids w' de toutes les particules S' pour que la somme de tous ces poids soit égale à un. Par 25 exemple, le calculateur 14 calcule la somme W de tous les poids w' puis divise chaque poids w' par la somme W. [0082] Lors d'une étape 132, le calculateur 14 ré-échantillonne les particules S'. Cette étape de ré-échantillonnage consiste à supprimer les particules dont les poids w' sont devenus trop faibles pour les remplacer par des particules associées à des 30 paramètres dont les poids sont plus forts. [0083] De nombreuses techniques de ré-échantillonnage sont connues. Par exemple, ici, on applique la méthode SIR (Sequential Importance Resampling) décrite dans le livre suivant : B. Ristic, S. Arulampalam, N. Gordon, "Beyond the Kalman Filter, particle filter for tracking applications", Artech House, 2004. 35 [0084] En résumé, cette méthode de ré-échantillonnage consiste d'abord à classer les particules en deux groupes : les particules à régénérer, c'est-à-dire toutes les particules dont le poids est inférieur à un seuil prédéterminé et les particules survivantes, c'est-à-dire les autres particules. [0085] Ensuite, pour chaque particule à régénérer, le calculateur 14 supprime 40 l'ancienne particule puis génère une nouvelle particule pour la remplacer. Pour 3037415 23 générer une nouvelle particule, le calculateur 14 tire aléatoirement une particule dans le groupe des particules survivantes. Ce tirage au sort est réalisé de préférence de telle sorte que la probabilité d'être tiré soit proportionnelle au poids de la particule survivante. 5 [0086] Ensuite, la position PIk et les valeurs a'k et E'k de la particule survivante tirée sont affectées à la nouvelle particule générée. Lors de cette étape, de préférence, le calculateur ajoute une perturbation aléatoire sur la position et sur les valeurs des facteurs correctifs en utilisant, par exemple, les variables aléatoires px, py, pa et pE. Toutefois, ces perturbations sont calibrées de manière à ce que les valeurs a'k et E'k 10 affectées à la particule générée restent très proches des valeurs actuelles des facteurs correctifs associés à la particule survivante tirée. Typiquement, la moyenne des valeurs a'k affectées aux particules générées est plus proche de la moyenne des valeurs a'k des particules survivantes que de la moyenne des valeurs a'k des particules supprimées. Il en est de même pour les valeurs E'k affectées aux particules 15 générées. [0087] Un nouveau poids est également affecté à chaque particule générée et, éventuellement, à la particule survivante dont elle est issue. [0088] Après le ré-échantillonnage, lors d'une étape 134, pour chaque particule S' survivante, le calculateur 14 enregistre dans la mémoire 12 la position Pik, les valeurs 20 a'k et E'k et son poids w'. [0089] Ensuite, lors d'une étape 136, le calculateur 14 estime la position PAk de l'appareil 10 à partir des positions PIk et des poids w' de toutes les particules S'. De nombreuses méthodes sont possibles pour faire cela. Par exemple, la position PAk est prise égale à celle de la particule S' ayant le poids w' le plus fort. Dans un autre 25 mode de réalisation, la position PA, est prise égale à la moyenne des positions PIk des particules S' en pondérant chaque position PIk par le poids w'. Une autre méthode est également décrite dans la demande W02012158441. [0090] Enfin, lors d'une étape 138, un point est affiché à la position PAk sur la représentation graphique 30 de la carte 16 pour indiquer au piéton 4 sa position sur 30 cette carte et donc sa position à l'intérieur du bâtiment 2. [0091] Les étapes 124 à 138 sont réitérées en boucle tant que le calculateur 14 reste dans sa phase active. Au fur et à mesure de ces itérations, la détérioration de la précision de l'estimation de la position de l'appareil 10 est limitée au maximum car seules les positions PIk les plus probables sont conservées lors de l'étape 132 de ré- 35 échantillonnage. [0092] De plus, puisque les valeurs des facteurs correctifs a' et c' convergent vers les véritables valeurs des biais de direction et de pas au fur et à mesure des itérations des étapes 124 à 138, les facteurs correctifs corrigent de plus en plus précisément ces biais de direction et de pas. Dès lors, la présence de ces facteurs correctifs dans 3037415 24 la loi de déplacement permet d'augmenter sensiblement la précision de l'estimation de la position PA, même en présence de tels biais de direction et de pas. [0093] A un instant t1, le calculateur 14 bascule, par exemple automatiquement, vers sa phase d'indisponibilité. Par exemple, ce basculement vers la phase 5 d'indisponibilité est provoqué par l'absence d'interaction entre le piéton 4 et l'appareil 10 pendant une durée supérieure à un seuil prédéterminé. L'exécution de l'algorithme 18 est donc interrompue. Lorsque le calculateur 14 bascule vers sa phase d'indisponibilité, par exemple, il en informe le calculateur 15. Ici, ce basculement du calculateur 14 vers sa phase d'indisponibilité permet d'économiser de l'énergie et 10 donc de limiter la décharge de la batterie de l'appareil 10. Par la suite, les dernières positions Pik, les derniers valeurs \mil, des poids wi et les dernières valeurs ail, et Eik enregistrées pour chaque particule Si juste avant l'instant t1, sont notées, respectivement, Pin, Mi et aiti et Eiti. [0094] Pendant la phase d'indisponibilité, le calculateur 15 continue d'exécuter les 15 étapes 102 à 108 et exécute en plus l'étape 112. Typiquement, la phase d'indisponibilité dure plusieurs intervalles At. Par exemple, elle dure au moins une 1 s et souvent plus de 1 min ou 5 min. [0095] A un instant t2, le calculateur 14 redevient disponible pour exécuter l'un quelconque des algorithmes 18 et 20. La phase d'indisponibilité s'interrompt donc à 20 cet instant. Par exemple, l'arrêt de la phase d'indisponibilité est provoqué par une interaction du piéton 4 avec l'interface homme-machine de l'appareil 10. A cet instant, l'algorithme 18 doit de nouveau être initialisé, puisque entre les instants t1 et t2, l'appareil 10 a pu être déplacé à l'intérieur du bâtiment 2. Pour cela, le calculateur 14 exécute d'abord une phase 148 de réveil actif. Dans le cas particulier décrit ici, cette 25 phase 148 vise à initialiser pour chaque particule Si sa position Pit2, son poids wit2 et les valeurs initiales ait2 et Eit2 des facteurs correctifs a' et ci en utilisant les données enregistrées par le calculateur 15 pendant la phase d'indisponibilité. Par la suite, pour simplifier la description, on suppose ici que les paramètres wit2, ait2 et Eit2 sont pris égaux, respectivement, aux paramètres Mi, aiti et Eiti. 30 [0096] Pour déterminer la valeur de la position Pit2 à l'instant t2, le calculateur 14 exécute, entre les instants t2 et t2 + At2, l'algorithme 20. [0097] Dans ce mode de réalisation, pour corriger le trajet relevé entre les instants t1 et t2, l'algorithme 20 utilise un graphe 150 (Figure 5) préenregistré dans la mémoire 12. Le graphe 150 code tous les trajets possibles à l'intérieur du bâtiment 2. Pour 35 simplifier la figure 5, seule une partie du graphe 150 codant les trajets possibles à l'intérieur d'une seule pièce 152 a été représentée. Cette pièce 152 est délimitée par des murs 154 à 157 et comporte également, à l'intérieur, deux obstacles infranchissables 160 et 162. [0098] Dans le graphe 150, chaque trajet possible comporte une succession de 40 noeuds successifs Nm reliés les uns aux autres par des segments rectilignes A. Par la 3037415 25 suite, le premier noeud et le dernier noeud d'un trajet possible sont appelés, respectivement « noeud de départ » et « noeud d'arrivée » du trajet. De façon similaire, pour chaque segment AJ, qui s'étend entre deux noeuds, le noeud le plus proche du noeud de départ le long du trajet est appelé « noeud amont ». Tandis que le 5 noeud le plus éloigné de ce noeud de départ est appelé « noeud aval ». Chaque noeud est situé à une position PNm qui peut être occupée par l'appareil 10. La position PNm de chaque noeud Nm est exprimée dans le repère XYZ. Chaque position PNm est connue. Par exemple, ici, ces positions PNm sont choisies pour que la distance entre deux noeuds immédiatement consécutifs dans un trajet possible soit toujours égale à une 10 constante D. Ici, cette constante D est égale à 1m. Les segments AJ ne traversent jamais un obstacle infranchissable. Ainsi, l'existence d'un segment entre deux noeuds indique que l'appareil 10 peut passer librement d'un de ces noeuds à l'autre. Chaque segment AJ est associé à un poids wAJ. Par exemple, le poids wA, indique la longueur du segment A. Dans cet exemple, les poids wA, sont tous égaux. 15 [0099] Sur la figure 5, seuls huit noeuds Ni à N8 ont été représentés, reliés les uns aux autres par huit segments Al à A8. [00100] L'algorithme 20 comporte les étapes suivantes pour d'obtenir une position initiale Plt2 associée au poids w't2 et aux valeurs a't2 et E't2. [00101] Éventuellement, lors d'une étape 198, le calculateur 14 sélectionne un 20 nombre M limité de particules parmi les particules S' et applique les étapes suivantes uniquement aux particules S' sélectionnées. Par exemple, le calculateur 14 sélectionne uniquement les particules S' dont le poids w'ti est supérieur à un seuil prédéterminé S. [00102] Ensuite, lors d'une étape 200, le calculateur 14 corrige le trajet, relevé 25 pendant la phase d'indisponibilité puis extrait de la mémoire 44, à l'aide des facteurs correctifs a' et c' associés à la particule S'. Cette étape est maintenant décrite dans le cas particulier du trajet représenté sur la figure 6. La figure 6 représente un exemple simplifié d'un trajet relevé 202 comportant seulement trois points caractéristiques 204 à 206. Les points 204 et 206 sont, respectivement, les points de départ et d'arrivée du 30 trajet 202. Les points 204 et 205 sont connectés l'un à l'autre par un segment 208, et les points 205 et 206 sont connectés par un segment 210. La direction de parcours du trajet 202 est représentée par des flèches dans la figure 6. [00103] Pour commencer, le calculateur 14 corrige la longueur du segment 208 en le multipliant par (1 + On tient compte ainsi du biais de pas estimé pour cette 35 particule S' lors de la phase 122. Ensuite, le calculateur 14 fait tourner d'un angle a'n le segment 208 autour d'un axe vertical passant par le point 204. Ainsi, le calculateur 14 corrige le biais de direction estimé pour cette même particule S' lors de la phase 122. [00104] Le segment corrigé ainsi obtenu porte la référence 212 sur la figure 7. Le 40 segment corrigé 212 pointe, depuis le point 204 vers un nouveau point caractéristique 3037415 26 214. Le point 214 est identique au point 205, sauf que sa position a été corrigée à l'aide des facteurs correctifs cin et ain. [00105] Le même traitement est ensuite appliqué au segment suivant 210. On obtient ainsi le segment corrigé 216 de la figure 7. Ce segment 216 pointe depuis le point 5 caractéristique corrigé 214 vers un nouveau point caractéristique 218. Le point 218 est identique au point 206, sauf que sa position a été corrigée à l'aide des valeurs a in et cin. Les points 206, 214 et 218 et les segments 212 et 216 forment un trajet 220 corrigé. Lors de l'étape suivante, c'est ce trajet 220 qui est utilisé et non pas le trajet non corrigé 202. 10 [00106] On notera que selon la particule Si traitée, le trajet corrigé obtenu n'est pas nécessairement le même, puisque les valeurs ain et cin peuvent être différentes d'une particule à l'autre. [00107] Ensuite, lors d'une étape 228, le calculateur 14 corrige une nouvelle fois le trajet relevé mais cette fois à partir des obstacles infranchissables associés à la carte 15 16. Pour cela, il procède aux opérations suivantes. [00108] La position Pin ne coïncide généralement pas exactement avec la position PNm d'un noeud du graphe 150. Ainsi, lors d'une opération 230, le calculateur 14 associe la particule Si à un ou plusieurs noeuds du graphe 150 les plus proches. Cette étape est illustrée à l'aide de la figure 5. Sur la figure 5, en plus du graphe 150, 20 on a représenté la position Pin de la particule Si. La position Pin ne coïncide avec aucune des positions PNm des noeuds du graphe 150. Le calculateur 14 associe alors la particule Si aux seuls noeuds les plus proches de la position Pin. Ici, si plusieurs noeuds se trouvent à moins de X mètres de la position Pin, chacun de ces noeuds est associé à la particule Si. Par exemple, la distance X est constante et inférieure à 2m 25 ou à 1m. Ici, la distance X est prise égale à 1m. Dans le cas où aucun noeud ne se trouve à moins de X mètres de la position Pin, la particule Si est associée au noeud le plus proche. Dans l'exemple représenté sur la figure 5, la particule Si est associée aux noeuds Ni et N2 car elle est à moins d'un mètre de ces deux noeuds. [00109] Lors d'une opération 232, le calculateur 14 affecte un poids à chaque noeud 30 associé à la particule Si. Ici, la somme des poids affectés aux noeuds associés à la même particule Si est égale au poids \Afin. De plus, de préférence, le poids affecté à chaque noeud est fonction de la distance qui le sépare de la position Pin. Ici, les poids WiN1 et WiN2, respectivement des noeuds Ni et N2, sont pris égaux, respectivement à (d1*wi1)/(d1+d2) et (d2*wi1)/(d1+d2), où dl et d2 sont les distances qui séparent la 35 position Pin, respectivement, des positions des noeuds Ni et N2. [00110] Lors d'une opération 234, pour chaque noeud associé à la particule Si, le calculateur 14 identifie le ou les noeuds avals les plus probables. Cette étape est illustrée à l'aide de la figure 8. Cette étape va maintenant être décrite dans le cas particulier du noeud Ni et du segment 212, sachant que ce qui est décrit dans ce cas 40 particulier s'applique à l'identique aux autres noeuds associés à la particule Si. 3037415 27 [00111] Le segment 212 est déplacé pour que son point amont soit confondu avec la position PNi. On obtient ainsi un segment 236 représenté sur la figure 8. Le point aval du segment 236 porte la référence 238 et occupe une position P238. Le calculateur 14 associe le point 238 au(x) noeud(s) le plus proche. Par exemple, il procède de façon 5 identique à ce qui a été décrit lors de l'opération 230. Ici, cela conduit à associer le point 238 aux noeuds N3, N4 et N7. On considère ici que le noeud N8 est à plus de 1 m de la position P238. En associant plusieurs noeuds avals au même noeud amont, on définit plusieurs chemins hypothétiques susceptibles d'avoir été suivis par l'appareil 10. Par exemple, dans le cas illustré ici : 10 - le noeud N3 correspond à un premier chemin hypothétique égal à la réunion des segments Al et A2, - le noeud N4 correspond à un deuxième chemin hypothétique égal à la réunion des segments Ai, A2 et A3, et - le noeud N7 correspond à un troisième chemin hypothétique égal à la réunion des 15 segments Ai, A5 et A7. [00112] Ici, pour simplifier, le calculateur 14 sélectionne uniquement le chemin hypothétique le plus direct entre les noeuds amont et aval. [00113] Lors d'une opération 240, le calculateur 14 affecte un poids à chaque noeud aval. Par exemple, le calculateur 14 procède comme décrit lors de l'opération 232, 20 sauf que la somme des poids des noeuds avals est égale au poids du noeud amont. Ainsi, la somme des poids affectés aux noeuds N3, N4 et N7 est égale au poids du noeud Ni. On procède aussi de la même façon pour les noeuds avals obtenus en sélectionnant le noeud N2 comme étant le noeud amont du segment 212. [00114] Ensuite, les opérations 234 et 240 sont réitérées en prenant comme noeud 25 amont chacun des noeuds avals identifiés lors de la précédente exécution de l'opération 234 et en utilisant le segment suivant 216 du trajet corrigé à la place du segment précédent 212. [00115] Les opérations 234 et 240 sont réitérées jusqu'à ce que la totalité des segments du trajet corrigé ait été utilisée successivement et dans l'ordre où ils ont été 30 relevés pendant la phase d'indisponibilité. [00116] Une fois la réitération des opérations 234 à 240 terminée, on dispose généralement de plusieurs chemins hypothétiques susceptibles d'avoir été suivis par l'appareil 10 en partant de la position P't1. Lors d'une opération 250, le calculateur 14 sélectionne alors le chemin hypothétique le plus probable. Dans le cas décrit ici, le 35 chemin le plus probable est celui dont le noeud d'arrivée est associé au poids le plus élevé. Ce chemin sélectionné est alors le nouveau trajet corrigé obtenu à l'issue de l'étape 228. [00117] Enfin, lors d'une étape 252, le calculateur 14 construit la position initiale P't2 de la particule S' à partir de sa position P't1 et du trajet corrigé obtenu à l'issue l'étape 40 228. Par exemple, le calculateur déplace la particule S' en suivant le trajet corrigé 3037415 28 depuis la position P't1 jusqu'à son point d'arrivée. La position initiale P't2 est alors prise égale aux coordonnées, dans le repère XYZ, de ce point d'arrivée. [00118] On remarque qu'en procédant ainsi, on conserve la diversité des positions P'ti. En effet, il est très peu probable que deux particules ayant des positions P't1 5 différentes permettent d'obtenir des positions initiales P't2 identiques. [00119] Les étapes 200 à 252 sont appliquées à chacune des particules sélectionnées lors de l'étape 198. [00120] A l'instant t2 + At2, l'exécution de l'algorithme 20 est terminée et la phase 148 se termine. Le calculateur 14 dispose maintenant de positions initiales P't2 pour 10 chaque particule sélectionnée. Si le nombre de particules sélectionnées est insuffisant, alors le calculateur 14 peut exécuter l'étape de ré-échantillonnage sur les particules sélectionnées pour accroître leur nombre. Cette étape de ré-échantillonnage est par exemple identique à l'étape 132. [00121] Ensuite, le procédé retourne à la phase active 122. Pour cela, le calculateur 15 14 retourne à l'étape 124 et utilise en tant que précédentes positions P'kl les positions Plt2 obtenues à l'issue de la phase 148. [00122] La phase 148 est exécutée à chaque fois que le calculateur 14 bascule de sa phase d'indisponibilité vers sa phase active. [00123] La figure 9 représente l'ordonnancement dans le temps du fonctionnement 20 des calculateurs 14 et 15 lors de l'exécution du procédé de la figure 4. Plus précisément, sur cette figure, les flèches du haut et du bas représentent les phases de fonctionnement, respectivement, des calculateurs 14 et 15. [00124] Entre les instants to et t1 puis après l'instant t2, le calculateur 14 est dans sa phase active. A l'inverse, entre les instants t1 et t2, le calculateur 14 est dans sa phase 25 d'indisponibilité. La phase d'indisponibilité est représenté par une zone hachurée 260. Sur cette figure, à l'intérieur de chacune des flèches, les références numériques des étapes exécutées par, respectivement, les calculateurs 14 et 15 sont indiquées dans l'intervalle temporelle où elles sont exécutées. Les notations « 102-108 » et « 102- 112 » signifie qu'il s'agit, respectivement, des étapes 102 à 108 et des étapes 102 à 30 112. [00125] Différentes variantes du mode de réalisation principal vont maintenant être décrites. [00126] Variantes sur la structure de l'appareil 10 : [00127] La mémoire 12 peut être un ensemble de plusieurs mémoires physiques 35 distinctes. [00128] Dans une variante, l'appareil 10 comporte des capteurs supplémentaires capables de délivrer des informations supplémentaires sur la position de l'appareil 10. Par exemple, l'appareil 10 peut comporter un capteur GPS qui délivre une position absolue de l'appareil 10. Il peut aussi comporter un capteur capable de mesurer la 40 position absolue de cet appareil 10 par rapport à des bornes fixes fixées sans aucun 3037415 29 degré de liberté à l'intérieur du bâtiment 2. Ces bornes sont par exemple des émetteurs-récepteurs Wi-fi. Dans ce cas, le capteur mesure la puissance des ondes émises par ces bornes et calcule sa position par rapport à ces bornes, par exemple, par triangulation. Ensuite, cette position par rapport aux bornes est convertie en une 5 position de l'appareil 10 dans la carte 16 à partir de la connaissance de la position de ces bornes à l'intérieur du bâtiment 2. Dans ce cas, l'algorithme 18 et/ou l'algorithme 20 peuvent, en plus des mesures de la centrale inertielle 22, utiliser des mesures de ces capteurs supplémentaires pour estimer la position PAk. Par exemple, si une mesure approximative de la position de l'appareil 10 est disponible en utilisant 10 d'autres informations, alors le poids w' est augmenté si la position Pik est proche de cette position approximative et, au contraire, diminuée si la position Pik est éloignée de cette position approximative. [00129] Dans un autre mode de réalisation, la centrale inertielle comporte des capteurs plus sophistiqués qui fournissent une valeur moyenne de la variable 15 mesurée à un instant k ainsi qu'un écart-type sur cette mesure (voir par exemple Straub2010). Dans ce cas, les valeurs mesurées de l'angle ek et de l'amplitude lk sont obtenues en tirant au hasard les valeurs en utilisant une loi de probabilité gaussienne dont la moyenne et l'écart-type sont égaux à ceux transmis par les capteurs. Cela permet notamment de prendre en compte le bruit de mesure. 20 [00130] La centrale inertielle peut aussi comporter des capteurs supplémentaires tels qu'un gyroscope, un baromètre, un hydromètre, un thermomètre, un microphone, un capteur de proximité ou un capteur d'intensité lumineuse ou autres. [00131] L'appareil 10 peut prendre la forme de n'importe quel appareil qui peut être déplacé dans l'espace par un utilisateur et dont la position relative peut être mesurée. 25 De préférence, l'appareil 10 est un appareil transportable, directement à la main, par l'utilisateur. Il peut s'agir, en plus des exemples précédemment donnés, d'un téléphone filaire, d'un assistant digital personnel plus connu sous l'acronyme « PDA » (Personal Digital Assistant), d'une console de jeux, d'un dispositif de navigation, d'un dispositif de suivi de l'activité ou de la forme de l'utilisateur tels que ceux incorporés 30 dans un bracelet intelligent, d'une montre intelligente, et de tout autre appareil électronique portable par l'utilisateur. Par exemple, il peut aussi s'agir d'un appareil d'accès à internet mobile connu sous l'acronyme MID (Mobile Internet Device), d'un dispositif de navigation personnel connu sous l'acronyme PND (Personal navigation device), d'un appareil photo numérique, d'une caméra vidéo, de jumelles, d'un objectif 35 de longue focale (« telephoto lens » en anglais), d'un appareil multimédia portable pour jouer de la musique ou des vidéos ou pour écouter la radio, d'une télécommande ou d'une combinaison des appareils précédemment cités. [00132] L'appareil 10 peut être un appareil autonome ou, au contraire, fonctionner uniquement en combinaison avec un autre appareil portable ou non. Par exemple, 40 l'autre appareil est un ordinateur portable ou de bureau, un serveur informatique ou 3037415 30 autre. L'autre appareil communique alors avec l'appareil 10 par l'intermédiaire d'un réseau filaire ou sans-fil. Si une communication filaire est utilisée, celle-ci peut être de n'importe quel type tel qu'une liaison série, une liaison parallèle, une transmission par parquets, etc. De façon similaire, si une liaison sans-fil est utilisée, celle-ci peut être 5 de n'importe quel type comme une transmission sans fil via des ondes électromagnétiques ou des ondes infrarouges ou autres. Il peut aussi s'agir d'une liaison mettant en oeuvre une combinaison de technologies filaire et sans-fil. [00133] En variante, le calculateur 14 n'est pas plus puissant que le calculateur 15. Dans une autre variante, les calculateurs 14 et 15 sont intégrés à l'intérieur d'un 10 même composant électronique, voire même, sur un même substrat en silicium. Le calculateur 14 peut se présenter sous la forme d'un regroupement de plusieurs microprocesseurs distincts. Dans ce dernier cas, l'un de ces microprocesseurs distincts est utilisé pour exécuter l'algorithme 18 et un autre pour exécuter l'algorithme 20. 15 [00134] Le calculateur 15 n'est pas nécessairement intégré à l'intérieur du MPU 40. En variante, le calculateur 15 est une calculateur externe au MPU 40 et raccordé par un bus externe au MPU 40. Dans une autre variante, le calculateur 15 est composé de plusieurs microprocesseurs. Dans ce dernier cas, l'un de ces microprocesseur peut être intégré au MPU 40 pour traiter les mesures de la centrale 22 et les 20 enregistrer dans la mémoire 44 tandis qu'un autre microprocesseur, externe au MPU 40, réalise le reste des opérations comme l'opération 114. [00135] Dans un autre mode de réalisation, l'algorithme 20 est exécuté par un serveur distant et non pas par l'unité 11. Dans ce cas, l'unité 11 est raccordée par l'intermédiaire d'un réseau sans fil à ce serveur distant, et transmet à ce serveur 25 distant le trajet relevé pendant la phase d'indisponibilité. En réponse, le serveur distant exécute l'algorithme 20 puis renvoie les positions P't2 à l'appareil 10. Cette architecture est intéressante si le serveur distant exécute plus rapidement l'algorithme 20 que le calculateur 14. [00136] Dans une autre variante, le calculateur 15 bascule également entre une 30 phase active et une phase d'indisponibilité. Dans la phase active, le calculateur 15 acquiert les mesures de la centrale inertielle 22 et les traite comme décrit précédemment. Dans la phase d'indisponibilité, il n'acquiert aucune mesure de la centrale inertielle 22, et ne réalise donc pas les opérations de traitement précédemment décrites. Par contre, même dans ce cas, le calculateur 15 reste dans 35 sa phase active, tant que l'appareil 10 doit être localisé et, en particulier, pendant les phases d'indisponibilité du calculateur 14. Par exemple, le calculateur 15 bascule dans sa phase d'indisponibilité, uniquement lorsque aucune application de localisation de l'appareil 10 n'est en cours d'exécution par le calculateur 14. Le calculateur 15 peut aussi basculé automatiquement vers sa phase d'indisponibilité 40 dès que les capteurs de la centrale inertielle 22 sont éteints. 3037415 31 [00137] Variantes des étapes 102 à 112: [00138] En variante, les instants tk ne correspondent pas aux instants où l'utilisateur pose un pied parterre. Par exemple, les instants tk sont des instants séparés les uns des autres par une période prédéterminée. Dans ce cas, la détection des pas lors de 5 l'étape 104 peut être omise. [00139] Dans une autre variante, les pas sont détectés à partir des mesures d'un ou plusieurs autres capteurs comme, par exemple, les mesures du gyromètre. [00140] Lors de l'étape 108 d'autres grandeurs physiques supplémentaires sont établies puis exploitées ultérieurement pour augmenter la précision de la localisation. 10 Par exemple le calculateur 15 identifie un changement de vitesse lors d'un déplacement en ligne droite ou un piétinement de l'utilisateur. Dans une variante simplifiée, l'étape 108 est omise. [00141] Les étapes 108 et 110 peuvent être exécutées, au moins pour partie, en parallèle. 15 [00142] Lors de l'opération 114, l'ajout de points caractéristiques supplémentaires peut être omis. Cette opération 114 peut aussi être simplifiée en utilisant moins de conditions que celles précédemment décrites pour ajouter des points caractéristiques supplémentaires. Par exemple, seulement une seule ou deux de ces conditions choisies parmi l'ensemble des conditions précédemment décrites sont utilisées. Il est 20 aussi possible d'ajouter d'autres points caractéristiques supplémentaires en plus de ceux déjà cités en utilisant pour cela les mesures de tout type de capteur embarqué dans l'appareil 10. Par exemple, un point caractéristique supplémentaire est ajouté dès que la variation de la direction du champ magnétique mesurée par le magnétomètre 26 franchit un seuil prédéterminé. Dans un autre exemple, un point 25 caractéristique supplémentaire est enregistré à chaque fois qu'un changement d'étage est détecté. Ce changement d'étage est détecté, par exemple, à partir de l'accélération verticale et du baromètre de la centrale inertielle 22 et non pas à partir d'un simple changement de la direction de déplacement. [00143] D'autres méthodes sont possibles pour détecter un point caractéristique où 30 un changement important de direction se produit. Par exemple, un changement important de direction peut être détecté par le fait que la vitesse angulaire, mesurée par un gyromètre de l'appareil 10, dépasse un seuil prédéterminé. Un changement de direction peut aussi être détecté à partir d'une variation dans la direction du champ magnétique mesuré par un magnétomètre de l'appareil 10. 35 [00144] Le trajet entre deux points caractéristiques peut aussi être approximé, non pas par un segment rectiligne, mais par un polynôme de degré strictement supérieur à un, tel que par une Spline. Par exemple, le trajet entre deux points caractéristiques successifs est approximé par un polynôme de degré supérieur ou égale à deux, trois ou plus. On notera que le cas d'un polynôme de degré un correspond au mode de 40 réalisation précédemment décrit. Le degré du polynôme peut être fixé une fois pour 3037415 32 toute ou, au contraire, choisi dynamiquement pour minimiser une erreur d'approximation. Par exemple, le degré du polynôme utilisé pour approximer le trajet relevé entre deux points caractéristiques successifs est choisi pour que l'erreur d'approximation soit inférieure à un seuil prédéterminé. L'erreur d'approximation entre 5 le polynôme et le trajet relevé est par exemple calculée comme une erreur quadratique. Dans ce cas, le degré du polynôme utilisé entre deux points caractéristiques donnés n'est pas nécessairement le même que celui utilisé entre deux autres points caractéristiques. Le polynôme qui approxime le trajet relevé entre deux points caractéristiques successifs peut être enregistré dans la mémoire 44 en 10 enregistrant ses coefficients ou plusieurs des points qu'il interpole. [00145] En variante, l'opération 114 est remplacée par une autre méthode connue d'extraction de caractéristiques dans un signal relevé entre les instants t1 et t2. Dans une autre variante, le trajet relevé n'est pas compressé et l'opération 114 est omise. [00146] Le trajet compressé généré par l'opération 114 est enregistré dans la 15 mémoire 44 en plus du trajet non-compressé relevé lors de l'étape 112 ou, au contraire, remplace dans la mémoire 44 le trajet non-compressé relevé. Dans ce dernier cas, le trajet non-compressé relevé est effacé de la mémoire 44 au moment où le trajet compressé est enregistré dans cette mémoire. [00147] Si le trajet compressé et/ou le trajet relevé occupent trop de place dans la 20 mémoire 44, il est possible de prévoir que le calculateur 15 réveille le calculateur 14 pour qu'il bascule dans sa phase active et traite les données enregistrées dans la mémoire 44. Ensuite, le calculateur 14 retourne automatiquement dans sa phase d'indisponibilité et les données traitées par le calculateur 14 sont effacées de la mémoire 44. Ainsi, de l'espace mémoire est libéré pour l'enregistrement de la suite du 25 trajet relevé et/ou du trajet compressé. [00148] Dans une autre variante, l'étape 112 est mise en oeuvre en permanence, et pas seulement lorsque le calculateur 14 est dans sa phase d'indisponibilité. Dans un mode de réalisation particulier de cette variante, le relevé du trajet est systématiquement mis en oeuvre mais l'opération 114 n'est exécutée que si le 30 calculateur 14 est dans sa phase d'indisponibilité. [00149] Variantes de l'étape 120 et de la phase active 122: [00150] L'étape 120 d'initialisation de l'algorithme 18 peut être réalisée différemment. Par exemple, la position initiale de l'appareil 10 dans le bâtiment 2 est initialisée en tenant compte de la dernière position absolue mesurée à l'aide d'un GPS ou à partir 35 de bornes Wifi. [00151] L'algorithme 18 peut être un autre algorithme qu'un filtre particulaire. Par exemple, l'algorithme 18 peut être un filtre de Kalman utilisé pour estimer la position soit de chaque particule, soit directement de l'appareil 10. L'algorithme 18 peut aussi être remplacé par un algorithme IMM (Interactive Multiple Models). 3037415 33 [00152] Dans le cas où l'algorithme 18 est un filtre particulaire, d'autres modèles de déplacement sont possibles. Par exemple, la distance parcourue entre deux positions P'k_i et P'k peut être obtenue en intégrant l'accélération mesurée pour obtenir la vitesse de déplacement entre ces deux positions, puis en intégrant cette vitesse entre 5 les instants tk_i et tk. Il existe également d'autres façons d'implémenter le filtre particulaire. Par exemple, on pourra se référer à ce sujet à la thèse Straub2010 ou aux demandes WO 2012158441 et US 8548738. [00153] D'autres contraintes sur les déplacements de l'appareil 10 sont utilisables que celles précédemment décrites. Par exemple, ces autres contraintes sont identiques à 10 celles décrites dans la demande FR1455575. Ces autres contraintes sont aussi bien utilisables lors de la phase 122 que lors de l'étape 228 pour corriger le trajet relevé. Par exemple, en variante, la carte 16 ne code pas la position des murs à l'intérieur du bâtiment 2. Au lieu de cela ou en plus, la carte 16 contient : - une cartographie de la puissance des ondes Wifi mesurables dans différentes zones 15 du bâtiment 2, - une cartographie de la fréquentation mesurée dans différentes zones du bâtiment 2, - une cartographie des niveaux sonores ou des ambiances sonores présents dans différentes zones du bâtiment 2, - une cartographie de l'amplitude des champs magnétiques mesurables dans 20 différentes zones du bâtiment 2. [00154] De tels cartes sont utilisables comme décrit précédemment dans le cas particulier de la carte 16. Par exemple, l'appareil 10 mesure la puissance des ondes Wifi ou la fréquentation ou le niveaux sonores ou l'ambiance sonore ou l'amplitude du champ magnétique avant et après avoir été déplacé d'un pas. Si d'après ces mesures 25 ce déplacement a fait passer l'appareil 10 d'une zone de la carte à une autre, alors le poids associé à toutes les particules S' pour lesquels les positions P'k_i et P'k sont restées dans la même zone est diminué. A l'inverse le poids des particules S', pour lesquels les positions Pik_i et P'k se trouvent dans des zones différentes de la carte susceptibles de correspondre au déplacement mesuré par l'appareil 10, est soit 30 augmenté soit laissé inchangé. [00155] Lors de l'étape 124, le calculateur 14 peut recevoir des données supplémentaires relatives au déplacement de l'appareil 10 pour calculer le déplacement de chaque particules S'. Par exemple, le calculateur 14 peut recevoir plusieurs angles utiles pour déterminer la direction ek tels que : 35 - un angle d'orientation de l'appareil 10, - un angle de correction entre l'orientation de l'appareil 10 et la direction de la marche, plus connu sous le terme anglais de « Misalignement angle », ou - un angle, mesuré à l'aide du magnétomètre, entre un repère dans lequel est calculé le trajet relevé et un repère de référence. 3037415 34 [00156] Lors de l'étape 124, le déplacement peut aussi être donné sous la forme d'une vitesse vx , vy ou sous la forme d'incréments de position Ax , y. [00157] Des informations additionnelles peuvent être transmises au calculateur 14 lors de cette étape 124 tels qu'un indicateur de changement de position de l'appareil 5 10, un indicateur de piétinement, la pression atmosphérique, etc. [00158] Lors de l'étape 126, les lois Lpxy, Lpa et Lp, peuvent être modifiées en fonction des informations additionnelles transmises par le calculateur 15 lors de l'étape 124. Par exemple, si un changement de position de l'appareil 10 est détecté, la dispersion de la loi Lpa est augmentée. 10 [00159] D'autres méthodes de ré-échantillonnage peuvent être mises en oeuvre lors de l'étape 132. Par exemple, toutes les particules S' sont ré-échantillonnées lors de cette étape et non pas seulement une partie. Lors du ré-échantillonnage des particules, de nombreux autres algorithmes pour déterminer la position initiale des particules régénérées sont connus et peuvent être utilisés à la place de celui décrit 15 précédemment. Par exemple, l'algorithme KLD (Kullbak-Leibler-Divergence) peut être utilisé. L'étape de ré-échantillonnage peut aussi être omise. [00160] Variantes sur le passage entre les phases 122 et 148: [00161] En variante, c'est le calculateur 14 qui indique au calculateur 15 qu'il va passer dans sa phase d'indisponibilité. Dans ce cas, le calculateur 15 n'a pas besoin 20 d'interroger le calculateur 14 pour déterminer s'il est ou non dans sa phase d'indisponibilité. [00162] D'autres causes que la mise en veille peuvent provoquer l'interruption de l'exécution de l'algorithme 18 par le calculateur 14 et le passage dans la phase d'indisponibilité. Quelle que soit la cause qui provoque l'interruption momentanée de 25 l'exécution de l'algorithme 18 par le calculateur 14, cette cause peut être traitée comme décrit ci-dessus. Ainsi, lorsque cette cause disparaît et que le calculateur 14 est de nouveau capable d'exécuter l'algorithme 18, il peut être correctement initialisé en exécutant la phase 148 de réveil actif. Par exemple, dans un mode de réalisation particulier, le calculateur 14 bascule vers la phase d'indisponibilité lorsqu'il est trop 30 occupé pour exécuter l'algorithme 18. Dès lors, lorsque le calculateur 14 est de nouveau disponible pour exécuter l'algorithme 18, la phase 148 est exécutée et l'algorithme 18 est correctement initialisé. De façon similaire, lorsque l'algorithme 18 est exécuté sur un serveur distant, la cause qui provoque l'arrêt de l'exécution de l'algorithme 18 peut être une interruption de la liaison entre l'appareil 10 et ce serveur 35 distant. [00163] Variantes de la phase 148 de réveil actif : [00164] En variantes, l'algorithme 20 est identique à l'algorithme 18, sauf que l'algorithme 20 est exécuté hors ligne et qu'il traite le trajet compressé et non pas directement toutes les mesures réalisées pendant la phase d'indisponibilité. Dans ce 40 cas, il est réglé pour tenir compte du fait que les distances entre les points 3037415 35 caractéristiques sont généralement plus grandes que la distance d'un pas. Par exemple, le nombre initial No de particules est beaucoup plus grand lors de l'exécution de l'algorithme 20, que lorsque l'algorithme 18 est exécuté. [00165] En variante, le calculateur 14 n'enregistre pas, avant l'instant t1, la position P'k, 5 la valeur w'k du poids w' et les valeurs a'k et E'k pour toutes les particules S' mais seulement pour un nombre limité N1 de particules pour limiter l'espace mémoire occupé par cet enregistrement. Par exemple, le calculateur 14 remplace chaque amas de particules par une seule particule S'n dont la position est égale au barycentre des positions PIk de chaque particule de cet amas. Le barycentre est calculé à partir 10 des positions PIk de chaque particule de cet amas pondéré par le poids w'k de cette particule. Il est aussi possible de pondérer chaque position P Ik par le même poids w. Un amas de particules est un groupe de particules, toutes situées à une distance inférieure à un seuil Sa du barycentre de cet amas. Il est donc possible de déterminer si une particule appartient ou non à un amas et ainsi d'identifier toutes les particules 15 qui appartiennent au même amas. Par exemple, le seuil Sa est inférieur ou égal à 2m ou 1m ou 50cm. La valeur w'ti poids w' de la particule S'n est calculé à partir des poids w' de toutes les particules S' appartenant à cet amas. Par exemple, le poids w'ti est égal à la moyenne des valeurs w'k des poids des particules appartenant à cet amas. [00166] Dans un mode de réalisation très simplifié, le calculateur 14 enregistre 20 seulement la position PAk estimée immédiatement avant l'instant t1 en tant que position P't1. Dans ce cas, il n'existe donc qu'une seule position P'n. [00167] En variante, les poids w' sont réinitialisés à la sortie de la phase d'indisponibilité. Par exemple, on affecte la même valeur w't2 à tous les poids w'. Typiquement, cette valeur est prise égale à 1/N2, où N2 est le nombre de particules S' 25 à l'instant t2 + At2. La valeur w't2 peut aussi être calculée en fonction de la valeur w'ti et du poids du noeud d'arrivée du chemin le plus probable utilisé pour construire la position PIt2. [00168] À la sortie de la phase d'indisponibilité, il est aussi possible de réinitialiser les valeurs des facteurs correctifs a' et c' comme décrit lors de l'étape 100. 30 [00169] L'étape 198 peut être omise. Dans ce cas, toutes les particules S' enregistrées à l'instant t1 sont utilisées. Dans un autre mode de réalisation, l'étape 198 est réalisée juste avant l'étape 134 d'enregistrement. Puis, lors de l'étape 134, seule les données des particules sélectionnées lors de l'étape 198 sont enregistrées. Dans ce cas, l'étape 198 n'est généralement pas réalisée une nouvelle fois au début 35 de la phase 148. [00170] Pour simplifier le procédé, l'étape 200 peut être omise. Dans ce cas, le trajet relevé pendant la phase d'indisponibilité est corrigé uniquement par l'exécution de l'étape 228. [00171] Dans un autre mode de réalisation, lors de l'étape 200, seul l'un des facteurs 40 correctifs a' et c' est utilisé pour corriger le trajet relevé et non pas les deux. Dans une 3037415 36 autre variante, le trajet est corrigé en tenant compte de toutes les valeurs a'n et c'ti de toutes les particules S'. Par exemple, le trajet est corrigé à partir de la moyenne de toutes les valeurs a'n et de la moyenne de toutes les valeurs c'ti. Cette moyenne peut être pondérée avec les poids w'ti de chaque particule S'. Dans ce cas, le trajet corrigé 5 est le même pour toutes les particules S'. [00172] Lors de l'étape 250, plutôt que de sélectionner seulement le chemin le plus probable, le calculateur 14 sélectionne les Ne chemins les plus probables, où Ne, est un entier strictement supérieur à 1 et généralement inférieur à 5 ou 3. Ensuite, l'étape 252 est appliquée pour chacun de ces Ne chemins les plus probables. On initialise 10 ainsi les positions P't2 de Nd, particules à partir d'une seule particule S, dont la position P't1 a été enregistrée. [00173] La position P't2 d'une particule à l'instant t2 n'est pas nécessairement confondue avec la position du point d'arrivée du trajet relevé. Par exemple, elle peut aussi être construite en ajoutant les variables aléatoires p'x et p'y aux coordonnées de 15 la position du point d'arrivée de ce trajet. Dans un autre mode de réalisation, on peut aussi ajouter les variables aléatoires p'x et p'y aux coordonnées du point de départ du trajet utilisé. [00174] Dans une autre variante de la phase 148, pour obtenir les positions P't2, le calculateur 14 fait parcourir à chaque particule S' dont la position a été enregistrée 20 juste avant l'instant t1, le trajet relevé pendant la phase d'indisponibilité sans exécuter l'étape 228. Ainsi, lors de ce parcours, le calculateur 14 ne tient pas compte des obstacles infranchissables. L'algorithme 20 n'utilise donc pas le graphe 150 des trajets possibles. Les positions de ces particules S' après qu'elles aient parcouru le trajet relevé pendant la phase d'indisponibilité sont alors considérées comme étant 25 les positions P't2 nécessaires à l'initialisation de l'algorithme 18 à l'instant t2 + At2. [00175] Dans une autre variante, les noeuds du graphe 150 ne sont pas uniformément répartis dans le bâtiment 2. Par exemple, le graphe 150 comporte un seul noeud par pièce ou par intersection entre deux pièces. [00176] En variante, pour limiter le nombre d'opérations réalisées lors de l'exécution 30 de l'algorithme 20, celui-ci est modifié pour : - arrêter d'explorer un chemin hypothétique dès que le poids associé à un noeud aval de ce chemin hypothétique devient inférieur à un seuil prédéterminé et cela avant que tous les segments du trajet corrigé aient été utilisés, et/ou - limiter le nombre de noeuds avals que le calculateur 14 peut identifier, par exemple, 35 à un seul noeud aval, et/ou - remplacer dans le trajet corrigé les segments dont la longueur est inférieure à un seuil prédéterminé par un seul point caractéristique, par exemple situé au milieu de ce segment court, et/ou - appliquer l'algorithme 20 uniquement aux particules dont le poids est strictement 40 supérieur à un seuil prédéterminé, et/ou 3037415 37 - appliquer l'algorithme 20 uniquement aux particules qui sont espacées les unes des autres d'une distance supérieure à un seuil prédéterminé, et/ou - fusionner plusieurs segments du trajet corrigé en un seul segment. [00177] La phase 148 peut aussi être exécutée en parallèle de la phase 122, c'est-à5 dire en même temps que le calculateur 14 est dans sa phase active et exécute l'algorithme 18. Dans ce cas, les positions estimées par l'exécution de l'algorithme 20 sont combinées aux positions estimées pour le même instant tk par l'algorithme 18 de manière à améliorer la précision de l'estimation de la position de l'appareil 10. Par exemple, la phase 148 est exécutée à intervalle régulier en parallèle de la phase 122. 10 [00178] Autres variantes : [00179] Le lancement de l'exécution du procédé de la figure 4 peut être déclenché automatiquement, par exemple, en réponse à la perte d'un signal GPS ou en réponse à la mesure d'un signal Wifi. [00180] La position estimée PA, n'est pas nécessairement affichée mais, en variante, 15 seulement transmise à un autre appareil pour traitement de cette information. [00181] Tout ce qui a été précédemment décrit s'applique également au cas où le piéton 4 est remplacé par une personne équipée d'un véhicule roulant comme un vélo, un chariot ou une poussette ou encore un fauteuil roulant. L'appareil 10 est alors soit transporté par cette personne soit directement déposé sur le véhicule roulant. 20 Dans ce cas, l'accéléromètre 24 peut être remplacé par un compte-tour qui mesure le nombre de tours de roue du véhicule roulant au fur et à mesure que cette personne se déplace dans le bâtiment 2. [00182] Le procédé précédemment décrit s'applique à tout type d'espace tridimensionnel ou bidimensionnel où l'appareil 10 peut être déplacé. Par exemple, il 25 peut aussi s'agir d'un espace situé en plein air et en dehors de tout bâtiment. Par exemple, ce procédé est utile pour localiser une personne dans un lieu où une localisation par GPS ou à partir de relais téléphonique est impossible. [00183] L'utilisation de points caractéristiques supplémentaires, telle que décrite précédemment, peut être mise en oeuvre indépendamment des différentes étapes 30 décrites pour initialiser l'algorithme 18 à la sortie de la phase d'indisponibilité. Par exemple, l'algorithme 20 peut être utilisé dans un appareil mobile dont les calculateurs ne basculent jamais dans une phase d'indisponibilité. [00184] Autres applications : [00185] On notera que les enseignements précédemment décrits dans le cadre d'un 35 appareil et d'un procédé de localisation de cet appareil déplacé à l'intérieur d'un espace tridimensionnel peuvent également être mis en oeuvre dans le cadre d'autres applications. Notamment, ces enseignements s'appliquent à tout appareil muni d'un ou de plusieurs capteurs délivrant des mesures, l'appareil étant associé à au moins un calculateur de mesure et un calculateur d'applications. Le calculateur de mesure, 40 généralement embarqué dans l'appareil et directement connecté aux capteurs, est 3037415 38 programmé pour fonctionner dans un mode « temps réel » et, en alternance, dans un mode de stockage. Dans le mode temps réel, le calculateur de mesure transfert, au fur et à mesure qu'elles sont réalisées, les mesures des capteurs à destination du calculateur d'application. Dans le mode de stockage, le calculateur de mesure relève 5 et enregistre, dans une mémoire qui lui est associée, des caractéristiques extraites des signaux mesurés par le ou les capteurs pour limiter le nombre d'informations à traiter par le calculateur d'applications. [00186] Le calculateur d'application bascule, automatiquement ou en réponse à une commande transmise par un utilisateur, entre une phase active et une phase 10 d'indisponibilité. Dans sa phase active, le calculateur d'application est apte à traiter les nouvelles mesures transmises par le calculateur de mesure au fur et à mesure qu'elles sont reçues. A l'inverse, dans sa phase d'indisponibilité, le calculateur d'application n'est pas capable de traiter les nouvelles mesures au fur et à mesure qu'elles lui sont transmises. De plus, le calculateur d'applications est programmé 15 pour, lorsqu'il est dans sa phase active : - exécuter un premier algorithme qui le conduit à recevoir les mesures transférées en temps réel par le calculateur de mesure et à les traiter au fur et à mesure de leur réception pour estimer en temps réel une nouvelle valeur d'une donnée, et, en alternance, 20 - exécuter un second algorithme qui le conduit à traiter les caractéristiques extraites, par le calculateur de mesure, pendant la phase d'indisponibilité, pour en déduire une valeur de cette même donnée mais à partir des mesures réalisées pendant cette phase d'indisponibilité. [00187] Dans ce cas plus général, le procédé de traitement des mesures s'articule 25 comme décrit en référence à la figure 10. Initialement, lors d'une étape 280, le calculateur de mesure acquière les nouvelles mesures réalisées par le ou les capteurs et les enregistrent dans une mémoire telle que la mémoire 44. [00188] Ensuite, lors d'une étape 282, le calculateur de mesure détermine si le calculateur d'applications est dans sa phase d'indisponibilité ou non. S'il est 30 déterminé que le calculateur d'applications est dans sa phase d'indisponibilité, alors le calculateur de mesure procède à une étape 284. Dans le cas contraire, le calculateur procède à une étape 286. [00189] Lors de l'étape 284, le calculateur de mesure bascule ou reste dans son mode de stockage des caractéristiques extraites des mesures du ou des capteurs. En 35 effet, il ne peut pas dans ce cas transmettre les nouvelles mesures au calculateur d'application au fur et à mesure qu'elles sont réalisées. [00190] Lors de l'étape 286, le calculateur de mesure détermine si le calculateur d'application vient juste de sortir de sa phase d'indisponibilité et revient à sa phase active. Dans l'affirmative, il procède à une étape 288. Dans le cas contraire, il procède 40 à une étape 290. 3037415 39 [00191] Lors de l'étape 288, le calculateur de mesure transmet les caractéristiques enregistrées lors de l'étape 284 au calculateur d'application. Ensuite, toujours lors de cette étape 288, le calculateur d'application exécute une phase de réveil actif lors de laquelle il traite les caractéristiques reçues de façon à estimer une valeur d'une 5 donnée à un instant tw compris dans l'intervalle de temps ]t1 ; t2]. L'intervalle de temps ]t1 ; t2] est l'intervalle de temps compris entre les instants t1 et t2 et qui comprend l'instant t2 mais qui ne comprend pas l'instant t1. Comme précédemment, les instants t1 et t2 sont, respectivement, l'instant où le calculateur d'application bascule dans sa phase d'indisponibilité et l'instant où le calculateur d'application 10 revient dans sa phase active. On notera que contrairement au mode de réalisation décrit précédemment, l'instant tw n'est pas nécessairement égale à l'instant t2. Par exemple, il peut être utile d'avoir la position de l'appareil 10 à un instant intermédiaire entre les instants t1 et t2. [00192] Ensuite, une fois la phase de réveil actif terminée, le procédé retourne à 15 l'étape 280. [00193] Lors de l'étape 290, le calculateur est dans sa phase active et la phase de réveil actif est terminée. Dans ce dernier cas, le calculateur de mesure est dans son mode temps réel et le calculateur d'application exécute un algorithme qui traite les nouvelles mesures au fur et à mesure qu'elles sont reçues. 20 [00194] A titre d'exemple, les capteurs sont des capteurs permettant le suivi de l'activité humaine. Par exemple, ces capteurs mesurent les déplacements et l'orientation du corps de la personne et/ou les déplacements des membres de cette personne les uns par rapport aux autres. Dans ce cas la donnée dont les valeurs successives doivent être estimées à partir des mesures des capteurs est une donnée 25 représentative d'une activité humaine telle que s'asseoir, marcher, courir, dormir ou autres. La mise en oeuvre du procédé de la figure 10 dans ce contexte permet au calculateur d'application d'estimer l'activité de la personne en prenant en compte des caractéristiques relevées pendant sa phase d'indiponibilité.De plus, dans ce type d'application, le passage de la phase d'indisponibilité à la phase active du calculateur 30 d'application peut être déclenchée par la détection, par le calculateur de mesure, d'une caractéristique associée à un mouvement de la personne reflétant son réveil, par exemple un mouvement d'un bras d'une amplitude dépassant un seuil prédéterminé. [00195] Dans un autre exemple d'application, le capteur électronique est un capteur 35 audio qui mesure un signal audio audible dans l'environnement immédiat de l'appareil. Dans ce cas, les premier et second algorithmes sont, par exemple, des algorithmes de reconnaissance vocale et les valeurs de la donnée qui doivent être estimées correspondent à des mots successivement prononcés par un être humain. Le premier algorithme effectue la reconnaissance vocale au fur et à mesure que les 40 nouvelles mesures du capteur sont reçues. Le second algorithme est quant à lui 3037415 40 conçu pour identifier des mots à partir de certaines caractéristiques, comme par exemple, un spectre de fréquences ou l'intensité du signal audio à certaines fréquences prédéterminées. Dans ce cas, lors de la phase d'indisponibilité, le calculateur de mesure est programmé pour extraire du signal audio mesuré les 5 caractéristiques nécessaires à la reconnaissance vocale par le second algorithme. Par exemple, ces caractéristiques sont les intensités du signal audio aux fréquences prédéterminées. A l'issue de la phase d'indisponibilité, pendant la phase de réveil actif, le calculateur d'application exécute le second algorithme pour identifier des mots prononcés juste avant ou au moment où ce calculateur d'application est revenu 10 dans sa phase active. Ensuite, suite au retour dans la phase active, le calculateur d'application exécute à nouveau le premier algorithme pour reconnaître les mots suivants à partir des mesures qui lui sont transmises par le calculateur de mesure. [00196] Dans une autre variante, le premier algorithme utilise les mots précédemment reconnus pour reconnaître les mots suivants. Dans ce cas, les mots prononcés 15 pendant la phase d'indisponibilité et identifiés pendant la phase de réveil actif suite à l'exécution du second algorithme sont utilisés, lors de la première exécution suivante du premier algorithme, pour aider à la reconnaissance des mots dans le signal audio mesuré après le retour à la phase active. Par exemple, le premier algorithme utilise, en plus des nouvelles mesures du signal audio, la connaissance des mots 20 précédemment identifiés par le second algorithme pour affiner et augmenter la fiabilité de la reconnaissance des mots suivants. En effet, il est connu qu'il est plus facile d'identifier un mot lorsqu'on connaît le contexte dans lequel il a été prononcé plutôt qu'en absence de toute information sur ce qui s'est dit juste avant. REVENDICATIONS1. Procédé de traitement des mesures d'au moins un capteur électronique, logé dans un appareil portable par un utilisateur, afin d'en déduire les valeurs d'une 5 donnée à des instants successifs dans le temps, dans lequel : - pendant une phase active, à chaque instant tk, un premier calculateur exécute (122 ; 290) un premier algorithme d'estimation d'une nouvelle valeur Pik de la donnée, cette exécution comprenant : - la réception (124) d'une nouvelle mesure d'une grandeur physique, réalisée 10 par le capteur électronique, cette nouvelle mesure étant plus récente que la mesure reçue lors de la précédente exécution du premier algorithme, puis - l'estimation (126 - 132) de la nouvelle valeur Pik de la donnée à partir de la nouvelle mesure reçue, - à un instant t1, l'arrêt de la phase active et le passage dans une phase 15 d'indisponibilité qui dure plusieurs instants tk et pendant laquelle l'exécution par le premier calculateur du premier algorithme est inhibée, caractérisé en ce que : - pendant la phase d'indisponibilité, un second calculateur : - acquière (102; 280) les nouvelles mesures de la grandeur physique réalisées 20 par le capteur électronique pendant cette phase d'indisponibilité, - relève (114; 284) des caractéristiques du signal mesuré par le capteur électronique à partir de ces nouvelles mesures reçues pendant cette phase d'indisponibilité, et - enregistre (114; 284) ces caractéristiques, le nombre de caractéristiques 25 relevées pendant la phase d'indisponibilité étant strictement inférieur au nombre de nouvelles mesures acquises pendant cette phase d'indisponibilité, - à un instant t2, l'arrêt de la phase d'indisponibilité et le passage dans une phase de réveil actif lors de laquelle le premier calculateur exécute (148) un second algorithme d'estimation d'une valeur P', de la donnée à un instant tw compris dans l'intervalle ]t1 ; 30 t] à partir des caractéristiques relevées et enregistrées pendant la phase d'indisponibilité par le second calculateur, puis - l'arrêt de la phase de réveil actif et le retour à la phase active. 2. Procédé selon la revendication 1, dans lequel : 35 - lors de la phase active, l'estimation (126 - 132) de la nouvelle valeur P ik de la donnée est réalisée à partir de la nouvelle mesure reçue et également à partir de la précédente valeur PIk-i de la donnée estimée lors de la précédente exécution du premier algorithme à l'instant tk_i, et 3037415 42 - lors de la première exécution du premier algorithme après l'arrêt de la phase de réveil actif, la précédente valeur Pik-i de la donnée nécessaire pour estimer la nouvelle valeur Pik de la donnée est construite à partir de la valeur P', fournie par l'exécution du second algorithme d'estimation pendant la phase de réveil actif. 3. Procédé selon la revendication 2, dans lequel : - le capteur électronique comporte une centrale inertielle embarquée dans l'appareil portable et la grandeur physique mesurée par cette centrale inertielle est représentative de la direction de déplacement de l'appareil et de l'amplitude de ce 10 déplacement entre un instant précédent tk_i et l'instant tk, - la donnée est la localisation de l'appareil portable à l'intérieur d'un espace et chaque valeur de la donnée comporte au moins une position Pik possible de l'appareil à l'intérieur de cet espace, - le premier et le second algorithmes d'estimation sont des algorithmes d'estimation 15 d'au moins une position possible Pik de l'appareil à l'intérieur de cet espace, - à chaque instant tk, l'exécution (122) du premier algorithme, par le premier calculateur, génère au moins une position possible Pik à partir d'au moins une position possible Pik-i estimée à l'aide de ce premier algorithme à l'instant précédent tk_i, - avant l'instant t1, au moins une dernière position possible P't1 est construite à partir 20 de la position possible Pik générée à l'instant tk précédent immédiatement l'instant t1, puis enregistrée, - pendant la phase d'indisponibilité, les caractéristiques du signal mesuré par la centrale inertielle qui sont relevées comportent une suite temporellement ordonnée de caractéristiques d'un trajet parcouru par l'appareil pendant la phase 25 d'indisponibilité, cette suite de caractéristiques constituant un relevé du trajet parcouru par l'appareil pendant la phase d'indisponibilité, - pendant la phase de réveil actif, l'exécution (148) du second algorithme génère au moins une position initiale P't2, pour cela la position initiale P't2 est estimée à partir de la dernière position construite P't1 et du trajet parcouru relevé pendant la phase 30 d'indisponibilité, et - lors de la première exécution du premier algorithme après l'instant t2, la position précédente PIk-i est construite à partir de la position initiale P't2 . 4. Procédé selon la revendication 3, dans lequel : 35 - à chaque instant tk, l'exécution (122) du premier algorithme, par le premier calculateur, génère plusieurs positions possibles Pik à partir de plusieurs positions possibles Pik-i estimées à l'aide de ce premier algorithme à l'instant précédent tk_i, 3037415 43 - avant l'instant t1, plusieurs positions possibles Pin différentes sont construites à partir des positions possibles Pik générées à l'instant tk précédent immédiatement l'instant t1, puis enregistrées, - pendant la phase de réveil actif, l'exécution (148) du second algorithme génère 5 plusieurs positions initiales Pit2, pour cela chaque position initiale Pit2 est estimée à partir d'une dernière position construite Pin respective et du trajet parcouru relevé pendant la phase d'indisponibilité, et - lors de la première exécution du premier algorithme après l'instant t2, chaque position précédente Pik-i est construite à partir d'une position initiale possible Pit2 10 respective. 5. Procédé selon la revendication 4, dans lequel la construction desdites plusieurs dernières positions initiales P't1 comporte : - la sélection (198) de M positions les plus probables de l'appareil à partir de poids 15 associés aux positions estimées Pik générées à l'instant tk précédent immédiatement l'instant t1, où M est un entier supérieur à deux et strictement inférieur au nombre de positions possibles Pik générées par l'exécution (122) du premier algorithme à l'instant tk précédent immédiatement l'instant t1, et chaque poids associé à une position possible Pik respective représente la probabilité que l'appareil se trouve dans cette 20 position Pik, puis - l'enregistrement en tant que dernières positions possibles P't1 des seules M positions les plus probables identifiées. 6. Procédé selon la revendication 3 ou 4, dans lequel l'exécution (148) du second 25 algorithme comporte : - la correction (228) du trajet relevé pendant la phase d'indisponibilité à l'aide d'un jeu de contraintes prédéfinies sur les déplacements de l'appareil dans l'espace, cette correction comportant des modifications des distances entre des points du trajet relevé et des modifications des amplitudes des changements de direction au niveau 30 de ces points du trajet relevé pour obtenir un trajet corrigé qui n'enfreint aucune des contraintes prédéfinies, chaque modification étant réalisée en tenant compte aussi bien des points situés en amont qu'en aval du point affecté par cette modification de la distance ou de l'amplitude du changement de direction, puis - la construction (252) de l'estimation de la position initiale P't2 comporte : 35 - le positionnement de l'origine du trajet corrigé en fonction de la dernière position P't1 estimée par l'exécution du premier algorithme, puis - l'obtention de la position P't2 à partir de la position du point d'arrivée du trajet corrigé lorsque l'origine de ce trajet corrigé a été ainsi positionnée. 3037415 44 7. Procédé selon la revendication 3 ou 4, dans lequel : - à chaque instant tk, le premier algorithme calcule (126) également la valeur d'un facteur correctif d'un biais de mesure affectant l'un des capteurs embarqués à 5 l'intérieur de l'appareil, et - le procédé comporte la correction (200) du trajet parcouru relevé pendant la phase d'indisponibilité pour obtenir un trajet corrigé, cette correction étant réalisée en appliquant aux mesures reçues de ce capteur embarqué le facteur correctif calculé à l'instant tk précédent immédiatement l'instant t1, puis - la position P't2 est estimée en utilisant ce trajet corrigé. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel : - le capteur électronique comporte une centrale inertielle embarquée dans l'appareil portable et la grandeur physique mesurée par cette centrale inertielle est 15 représentative de la direction de déplacement de l'appareil et de l'amplitude de ce déplacement entre un instant précédent tk_i et l'instant tk, - la donnée est la localisation de l'appareil portable à l'intérieur d'un espace et chaque valeur de la donnée comporte au moins une position PIk possible de l'appareil à l'intérieur de cet espace, 20 - le premier et le second algorithmes d'estimation sont des algorithmes d'estimation d'au moins une position possible PIk de l'appareil à l'intérieur de cet espace, - pendant la phase d'indisponibilité, les caractéristiques du signal mesuré par la centrale inertielle qui sont relevées comportent des points caractéristiques d'un trajet parcouru par l'appareil pendant cette phase d'indisponibilité, le relevé (114) des 25 caractéristiques du signal mesuré comportant à cet effet: - l'identification, à partir des nouvelles mesures de la centrale inertielle reçues pendant la phase d'indisponibilité, de points caractéristiques du trajet moins nombreux que les nouvelles mesures de la centrale inertielle reçues pendant la phase d'indisponibilité, et 30 - pour chaque paire de points caractéristiques identifiés successifs le long du trajet, l'approximation du trajet de l'appareil entre ces deux points caractéristiques par un polynôme dont les coefficients sont déterminés. 9. Procédé selon la revendication 8, dans lequel l'identification de points 35 caractéristiques du trajet comporte l'identification, en tant que point caractéristique du trajet, de chaque point où l'amplitude d'un changement de direction de déplacement de l'appareil dépasse un seuil prédéterminé. 3037415 45 10. Procédé selon la revendication 9, dans lequel l'identification de points caractéristiques du trajet comporte également l'identification d'un point caractéristique supplémentaire à chaque fois qu'au moins l'une des conditions suivantes est satisfaite : 5 - la distance parcourue depuis le dernier point caractéristique identifié est supérieure à un seuil prédéterminé, - l'amplitude d'une variation de pression mesurée par un capteur de pression embarqué à l'intérieur de l'appareil est supérieure à un seuil prédéterminé, - l'amplitude d'un déplacement de l'appareil par rapport au corps d'une personne qui 10 transporte manuellement l'appareil est supérieure à un seuil prédéterminé, ou - le temps écoulé depuis le dernier point caractéristique identifié est supérieur à un seuil prédéterminé. 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel 15 pendant toute la durée de la phase d'indisponibilité, le premier calculateur est mis en veille pour économiser de l'énergie. 12. Procédé selon l'une quelconque des revendications 3 à 7, dans lequel, lorsqu'il est exécuté par le premier calculateur, le premier algorithme estime (126 - 132) la 20 position PIk de l'appareil en utilisant un jeu de contraintes prédéfinies sur les déplacements de l'appareil dans l'espace. 13. Procédé selon l'une quelconque des revendications précédentes, dans lequel le premier algorithme est un filtre particulaire. 25 14. Support (12) d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions pour exécuter un procédé conforme à l'une quelconque des revendications précédentes, lorsque ces instructions sont exécutées par un calculateur électronique. 30 15. Unité électronique de traitement des mesures d'au moins un capteur électronique afin d'en déduire les valeurs d'une donnée à des instants successifs dans le temps, cette unité comportant au moins un premier et un second calculateurs électroniques (14, 15), dans lequel : 35 - le premier calculateur (14) est programmé pour : - exécuter, pendant une phase active, à chaque instant tk, un premier algorithme d'estimation d'une nouvelle valeur PIk de la donnée, cette exécution comprenant : 3037415 46 - la réception d'une nouvelle mesure d'une grandeur physique, réalisée par le capteur électronique, cette nouvelle mesure étant plus récente que la mesure reçue lors de la précédente exécution du premier algorithme, puis - l'estimation de la nouvelle valeur Pik de la donnée à partir de la nouvelle mesure 5 reçue, - à un instant t1, arrêter de la phase active et passer dans une phase d'indisponibilité qui dure plusieurs instants tk et pendant laquelle l'exécution par le premier calculateur du premier algorithme est inhibée, caractérisé en ce que : 10 - le second calculateur est programmé pour, pendant la phase d'indisponibilité : - acquérir les nouvelles mesures de la grandeur physique réalisées par le capteur électronique pendant cette phase d'indisponibilité, - relever des caractéristiques du signal mesuré par le capteur électronique à partir de ces nouvelles mesures reçues pendant cette phase d'indisponibilité, et - enregistrer ces caractéristiques, le nombre de caractéristiques relevées pendant la phase d'indisponibilité étant strictement inférieur au nombre de nouvelles mesures acquise pendant cette phase d'indisponibilité, et - le premier calculateur est programmé pour : - à un instant t2, arrêter la phase d'indisponibilité et passer dans une phase de réveil actif lors de laquelle il exécute un second algorithme d'estimation d'une valeur P''' de la donnée à un instant tw compris dans l'intervalle ]t1 ; t] à partir des caractéristiques relevées et enregistrées pendant la phase d'indisponibilité par le second calculateur, puis - arrêter de la phase de réveil actif et retourner à la phase active. 16. Appareil directement transportable à la main par un être humain, cet appareil comportant au moins un capteur électronique (22) apte à mesurer une grandeur physique, caractérisé en ce que l'appareil comporte aussi une unité électronique (11) conforme à la revendication 15, cette unité électronique étant raccordée à ce capteur électronique.