FR 2976353 A1 20121214 FR 1154915 A 20110607 PROCÉDÉ D'ESTIMATION SIMPLIFIE DE L'ORIENTATION D'UN OBJET ET CENTRALE D'ATTITUDE METTANT EN OEUVRE UN TEL PROCÉDÉ La présente invention se rapporte à un procédé d'estimation de l'orientation d'un objet dans l'espace, animé ou non d'une accélération propre et soumis ou non à une perturbation magnétique (variations inconnues autour d'un champ principalement uniforme dans l'environnement proche et constant dans le temps), et à un dispositif mettant en oeuvre un tel procédé. L'estimation de l'orientation permet de résoudre toute une classe de problèmes. L'obtention de l'orientation nécessite généralement la mise en oeuvre de plusieurs modalités de capteurs, faisant partie d'un ensemble désigné par 1 o dispositif de capture de mouvement, également désigné par centrale d'attitude. Des capteurs MEMS (« Micro-Electro-Mechanical Systems » ou microsystèmes électromécaniques) peuvent être utilisés pour constituer cette centrale, ceux-ci présentent l'avantage d'être peu encombrants et de coût 15 réduit. Le capteur offrant une telle combinaison de capteurs MEMS est souvent désigné sous le nom de centrale inertielle ou IMU en anglais (Inertial Measurement Unit). Une centrale d'attitude est une centrale inertielle augmentée des moyens de traitement nécessaires pour obtenir l'orientation. L'emploi de tels capteurs MEMS permet d'envisager l'utilisation de centrales 20 d'attitude dans des domaines d'applications variés, notamment le domaine biomédical, pour la surveillance des personnes âgées à domicile, la rééducation fonctionnelle, dans le domaine sportif, pour l'analyse des mouvements des sportifs, dans les domaines automobile, robotique, de la réalité virtuelle, de l'animation en trois dimensions et plus généralement dans 25 tout domaine dans lequel on cherche à déterminer ou à observer un mouvement. Cependant, par rapport à des capteurs non MEMS (notamment des capteurs de fabrication classique, non issus des micro-technologies et utilisés par exemple, dans le domaine de la navigation), ces capteurs MEMS sont très 30 peu onéreux mais ont pour inconvénients d'être relativement peu performants, bruités et biaisés. Actuellement il existe plusieurs méthodes pour obtenir, à partir de mesures fournies par des accéléromètres, des magnétomètres et des gyromètres, une estimation de l'orientation de l'objet. Il existe notamment des méthodes mettant en oeuvre un observateur, qui permettent de fusionner une information double : celle provenant des mesures fournies par les capteurs et celle provenant d'un modèle d'évolution, et ceci tout en maintenant un temps de calcul compatible avec une implémentation en temps réel sur des processeurs relativement puissant (PC, Smartphone). Les méthodes connues mettant en oeuvre un observateur reposent principalement sur l'utilisation d'un filtre de Kalman. L'avantage de cette technique est de permettre la fusion des données tout en tenant compte de la qualité de l'information apportée par les mesures fournies par les capteurs et de la qualité du modèle d'évolution des états cinétiques. Parmi les différents types de filtre de Kalman, on utilise particulièrement le filtre de Kalman étendu (« Extended Kalman Filter » ou EKF en terminologie anglo-saxonne) ; celui-ci est rapide et simple à mettre en oeuvre, une de ses applications à la capture de mouvement est notamment décrite dans le document « Quaternion-based extended Kalman filter for determination orientation by inertial and magnetic sensing », SABATINI A.M., IEEE Transactions on Biomedical Engineering, 2006, 53(7). Outre le choix du filtre, la qualité des mesures injectées dans le filtre a une grande importance, et notamment la confiance qu'on accorde à leur valeur. En effet, les mesures comportent une partie informative directement liée à l'orientation de l'objet en mouvement et une partie perturbatrice dont la nature dépend du capteur considéré. Au premier ordre, il s'agit des accélérations propres pour les mesures fournies par les accéléromètres, des perturbations magnétiques pour les mesures délivrées par le magnétomètre et le biais pour les gyromètres. Il faut également prendre en compte le bruit de mesure. Cependant celui-ci est classiquement traité dans le filtre. Il existe plusieurs méthodes pour traiter les perturbations. On peut notamment ne pas ternir compte des mesures perturbées, mais dans le cas où les mesures de plusieurs capteurs présentent simultanément une perturbation, l'observateur n'a plus assez d'informations pour proposer une estimation correcte de l'orientation. On peut également inclure une représentation des perturbations dans le vecteur et le modèle d'état du filtre de Kalman pour les estimer conjointement avec les paramètres de l'attitude du solide. Toutefois, l'estimation conjointe des perturbations et de l'orientation s'avère délicate du fait d'un manque d'observabilité. Elle nécessite en outre le réglage d'un grand nombre de paramètres, ce qui complexifie sa mise en oeuvre. Ces inconvénients ont été surmontés par la demande de brevet internationale publiée sous le n cWO2010/007160, dont la demanderesse est co-titulaire, qui offre un procédé d'estimation de l'orientation fournissant une estimation précise de l'orientation, en présence ou non d'accélérations propres et de perturbations magnétiques, et ceci de manière simplifiée par rapport aux procédés existants. Le but énoncé précédemment est atteint par un procédé d'estimation de l'orientation sur la base de mesures selon les trois axes de l'espace de l'accélération, du champ magnétique et de la vitesse de rotation, comportant : - une étape de prétraitement de ces mesures pour détecter l'existence d'une perturbation et estimer des mesures non perturbées, et - une étape d'estimation de l'orientation sur la base des valeurs de mesures issues de l'étape de prétraitement. Ce procédé présente cependant l'inconvénient de nécessiter l'utilisation 20 d'une puissance de calcul importante qui n'est pas disponible dans un microcontroleur. C'est le but de la présente invention de procurer un procédé d'estimation des paramètres d'attitude d'un solide en mouvement, instrumenté par des capteurs de type accéléromètre, gyromètre et magnétomètre, qui remplace la 25 plupart des étapes de calcul du filtre de Kalman par une prise en compte directe de mesures desdits capteurs, corrigées ou non, selon que des perturbations sont détectées ou pas, dans le vecteur d'état. A cet effet, la présente invention divulgue un procédé d'estimation de 30 l'orientation d'un objet dans l'espace à un instant k utilisant au moins une mesure d'au moins un champ physique sensiblement uniforme dans le temps et dans l'espace et des mesures de la vitesse de rotation dudit objet selon trois axes de l'espace, ledit procédé comprenant une première étape de détection de la présence de perturbations dans les mesures dudit au moins 35 un champ physique et une deuxième étape de calcul d'un quaternion d'orientation dudit objet dans l'espace déduit d'au moins une représentation des mesures dudit au moins un champ physique en sortie de la première étape, ledit procédé étant caractérisé en ce que, lorsque en sortie de ladite première étape est détectée une perturbation de ladite au moins une mesure, ladite mesure perturbée en sortie de la première étape est combinée en entrée de la deuxième étape avec un quaternion d'orientation dudit objet calculé en sortie d'une troisième étape d'intégration d'une sortie corrigée de la mesure de vitesse de rotation de l'objet à partir d'un quaternion d'orientation dudit objet en sortie de la deuxième étape à l'instant k-1. Avantageusement, le au moins un champ physique est la gravité terrestre. Avantageusement, le au moins un champ physique est le champ magnétique terrestre. Avantageusement, le procédé de l'invention utilise des mesures d'au moins deux champs physiques, la gravité terrestre et le champ magnétique terrestre. 20 Avantageusement, le procédé de l'invention utilise des mesures d'un seul champ physique et des valeurs d'un champ synthétique construites à partir dudit seul champ physique, ledit champ synthétique étant défini tel qu'il ait au moins une composante orthogonale audit champ physique et que son produit vectoriel avec ledit champ physique soit non nul. 25 Avantageusement, ledit champ physique est la gravité GO et ledit champ synthétique est un vecteur constitué par une combinaison linéaire de ladite gravité GO et d' un vecteur non nul choisi dans un plan orthogonal à GO. 30 Avantageusement, ledit champ physique est le champ magnétique terrestre HO et ledit champ synthétique est un vecteur constitué par une combinaison linéaire dudit champ magnétique terrestre HO et d'un vecteur non nul choisi dans un plan orthogonal à HO.15 Avantageusement, le procédé de l'invention comprend en outre dans la première étape, une sous-étape de correction de la vitesse de rotation consistant à retrancher aux mesures de ladite vitesse de rotation un biais moyen déterminé lors d'une étape de position au repos et à lui substituer une vitesse de rotation corrigée. Avantageusement, la première étape de détection d'une perturbation dans les mesures dudit au moins un champ consiste à déterminer si la différence de la norme des mesures dudit au moins un champ et de leurs valeurs de 1 o référence est supérieure à une valeur de seuil choisie. Avantageusement, la deuxième étape de calcul d'un quaternion d'orientation dudit objet s'effectue par composition d'au moins deux mesures de champs physiques en sortie de la première étape et utilise un algorithme de type 15 TRIAD. Avantageusement, la deuxième étape de calcul d'un quaternion d'orientation dudit objet s'effectue par composition de mesures de plus de deux mesures de champs physiques en sortie de la première étape et utilise un algorithme 20 de type QUEST. Avantageusement, la troisième étape d'intégration d'une sortie corrigée de la mesure de vitesse de l'objet à partir d'un quaternion d'orientation dudit objet en sortie de la deuxième étape à l'instant k-1 est effectuée par 25 résolution de l'équation différentielle : w 0 L(t)- 1 Gq_ -2OL~G q 30 Avantageusement, le quaternion estimé de l'objet en sortie de la deuxième étape est normalisé. La présente invention divulgue également une centrale d'attitude d'un objet dans l'espace comprenant des moyens de mesure d'au moins un champ physique et de la vitesse de rotation dudit objet selon trois axes de l'espace, ladite centrale comprenant un premier module de prétraitement à l'instant k pour détecter la présence du au moins un champ physique et un deuxième module de calcul d'un quaternion d'orientation d'au moins une représentation de l'orientation dudit objet calculée à partir d'au moins une représentation des mesures dudit au moins un champ physique en sortie du premier module, ladite centrale étant caractérisée en ce que, lorsque en sortie dudit premier module est détectée une perturbation de ladite au moins une mesure ladite mesure perturbée en sortie du premier module est remplacée en entrée du deuxième module par une mesure combinée avec un quaternion d'orientation dudit objet calculé en sortie d'un troisième module d'intégration 1 o d'une sortie corrigée de la mesure de vitesse de l'objet à partir d'un quaternion d'orientation dudit objet en sortie du deuxième module à l'instant k-1. L'invention présente est utilisable pour toutes les combinaisons de capteurs 15 comportant : - au moins un capteur apte à restituer une information représentative de la vitesse angulaire (par exemple : un gyroscope, un couple d'accéléromètres séparés par un bras de levier, ...) - et au moins un capteur apte à restituer une information représentative 20 d'un champ physique Po principalement uniforme et constant dans le temps et l'espace Le dispositif de l'invention sera préférablement composé d'un gyroscope, d'un accéléromètre (apte à mesurer le champ terrestre gravitationnel Go) et d'un magnétomètre (apte à mesurer le champ terrestre magnétique Ho). 25 L'invention présente l'avantage de permettre une estimation de la vraie accélération.et des vraies perturbations magnétiques : MesuresAccelero=MatriceOrientation*(AccVraie-Go(0,0,1)) AccVraie=Transpose(MatriceOrientation)*MesuresAccelero+Go(0,0,1). L'invention présente en outre l'avantage de permettre le choix des capteurs 30 dont les mesures sont à privilégier en fonction de différents cas d'emploi. Notamment, dans le cas où une perturbation est détectée sur l'accéléromètre ou le magnétomètre, on privilégie les sorties de l'autre capteur et on réalise un traitement adapté qui est compatible avec une faible capacité de calcul et avec une faible capacité mémoire. 35 La présente invention sera mieux comprise à l'aide de la description qui va suivre et des dessins en annexe, sur lesquels : - la figure 1 représente une architecture fonctionnelle de mise en oeuvre du procédé selon l'art antérieur pour l'échantillon k ; - la figure 2 représente une architecture fonctionnelle de mise en oeuvre du procédé de l'invention pour l'échantillon k, dans un de ses modes de réalisation ; La figure 1 représente une architecture fonctionnelle de mise en oeuvre du 1 o procédé selon l'art antérieur pour l'échantillon k. Cette figure illustre la mise en oeuvre du procédé divulgué par la demande internationale publiée sous le n cWO2010/007160, dont les notations et certaines étapes sont utilisées dans le cadre de la mise en oeuvre de la présente invention. 15 On souhaite obtenir l'orientation d'un objet se déplaçant dans l'espace, par exemple l'orientation d'une personne (si le capteur est disposé sur le buste ou à la ceinture) ou d'un segment d'une personne (si le capteur est disposé sur la partie du corps à laquelle on s'intéresse). Pour cela, on utilise une centrale d'attitude comportant des capteurs aptes à fournir des mesures de 20 l'accélération totale, de champ magnétique et de la vitesse de rotation selon les trois axes de l'espace. Les capteurs sont avantageusement des capteurs MEMS offrant un prix de revient réduit et un encombrement limité. Ce procédé de l'art antérieur comporte une étape 110 d'initialisation de la centrale d'attitude, une étape 120 de prétraitement des mesures fournies par 25 les capteurs et une troisième étape de traitement 130 par un observateur, de type Extended Kalman Filter ou EKF. Selon l'invention, on peut diversifier les capteurs ; par exemple, pour la mesure de l'accélération, le système peut être constitué d'un accéléromètre 30 tri axe ou de trois accéléromètres mono axe fournissant une mesure sur chacun des axes. Notre invention décrit principalement le procédé utilisant les trois modalités de mesures précitées mais le raisonnement est le même pour un autre triplet de capteurs. Les propriétés minimales du dispositif doivent, de manière privilégiée, être les suivantes : - Un capteur de mesure yG à au moins un axe sensible, apte à restituer une mesure liée à la vitesse angulaire co par un modèle identifiable ; - Un capteur de mesure yA à au moins un axe sensible apte à restituer une mesure d'un champ de référence Go fixe dans le repère de référence et identifiable à un moment donné, auquel peut s'ajouter une variation a de même nature que Go mais non connue a priori ; Préférablement, ces capteurs sont respectivement : un gyromètre et un 1 o accéléromètre ou un magnétomètre. Il faut noter que dans le procédé s'appuyant sur le précédent dispositif, une étape de construction d'une mesure d'un champ physique synthétique orthogonal à Go sera alors réalisée. Un dispositif préférable un peu plus coûteux mais plus performant possède 15 les propriétés minimales suivantes : - Un capteur de mesure yG à au moins un axe sensible, apte à restituer une mesure liée à la vitesse angulaire co par un modèle identifiable ; - Un capteur de mesure yA à au moins un axe sensible apte à restituer une mesure d'un champ de référence Go fixe dans le repère de référence et 20 identifiable à un moment donné, auquel peut s'ajouter une variation a de même nature que Go mais non connue a priori ; - Un capteur de mesure yM à au moins un axe de mesure apte à restituer une mesure d'un champ de référence Ho fixe dans le repère de référence et identifiable à un moment donné, auquel peut s'ajouter une 25 variation m de même nature que Go mais non connue a priori. Préférablement, ces capteurs sont respectivement : un gyromètre, un accéléromètre et un magnétomètre. Il faut noter que ces trois types de capteurs fournissent des informations complémentaires en vue d'une meilleure estimation de l'orientation, 30 notamment si l'on dispose de capteurs tri-axiaux : - Le gyromètre, si on l'intègre, apporte des informations angulaires à trois degrés de liberté, mais en relatif, c'est-à-dire avec une propension à la dérive si ledit gyromètre est utilisé seul ; - L'accéléromètre apporte deux degrés de liberté angulaires absolus (roll et pitch), qui sont accessibles en statique si ledit accéléromètre est utilisé seul ; - Le magnétomètre apporte deux degrés de liberté angulaires (dont le 5 lacet ou yaw), qui sont accessibles dans un contexte non perturbé magnétiquement si ledit magnétomètre est utilisé seul. On peut généraliser la présente invention qui reste valable pour tout capteur axial apte à mesurer une direction uniforme de l'espace pour remplacer les 1 o capteurs de mesure du champ physique apportant une information additionnelle au gyroscope. On pense notamment aux mesures de DOA (Direction of Arrivai) pour tout type d'ondes électromagnétiques, en particulier les champs d'ondes émis par les balises que les opérateurs de téléphonie mobile utilisent pour mailler 15 l'espace de couverture de réseau, ou des cellules photoélectriques selon des modalités divulguées par la demande de brevet européenne publiée sous le n° EP1984696). Il peut s'agir aussi d'un dispositif utilisant des ondes ultrasonores planes provenant d'une source lointaine et dont la direction est mesurée et joue le rôle du champ physique. 20 Un des avantages du couple accéléromètre-magnétomètre est la capacité de chacun de ces deux capteurs à fournir ces deux DDL complémentaires. On peut cependant très bien utiliser d'autres capteurs vectoriels qui se complémentent de la même manière. L'accéléromètre est cependant nécessaire si on veut estimer l'accélération propre, la vitesse linéaire ou la 25 position... Pour les mesures de champ magnétique, les capteurs privilégiés peuvent être un magnétomètre tri axe ou trois magnétomètres mono axe. Pour la mesure de la vitesse de rotation, les capteurs privilégiés peuvent être trois gyromètres mono axe, avantageusement deux gyromètres bi axe ou un 30 gyromètre tri axe. Les triaxes peuvent être alignés ou non, dans ce dernier cas l'orientation relative entre les axes doit être connue. Dans la suite de la description, à des fins de simplicité, nous désignerons le ou les accéléromètres, par un accéléromètre, le ou les magnétomètres par 10 15 20 25 un magnétomètre et le ou les gyromètres par un gyromètre. Ces capteurs sont attachés à l'objet dont on souhaite connaître l'orientation. Nous disposons des mesures Y, que nous modélisons en notations 5 vectorielles par : YA=GR(GGO+'a)+ VA YM=G'R(GHo+'h)+vM (1) YG=Lw+b+vG Où les préfixes L et G précisent si le vecteur est exprimé dans le repère global {G} fixe ou dans le repère local {L} mobile. avec : YA : mesure tri-axe de l'accélération totale fournie par l'accéléromètre, YM : mesure tri-axe du champ magnétique fournie par le magnétomètre, YG : mesure tri-axe de la vitesse de rotation fournie par le gyromètre, notée également co : est la matrice de passage du repère global au repère local (matrice de rotation), Go : champ de gravité terrestre (vecteur 3x1), exprimé dans le repère fixe en général, soit (0,0,1), Ho : champ magnétique terrestre (vecteur 3x1), également exprimé dans le repère fixe en général, soit (0.5 ; 0 ; Vj 2 ) w : vitesse angulaire, a : accélérations propres ou variations perturbatrices d'accélération. h : perturbations magnétiques, b : biais du gyromètre, vA : bruit de mesure de l'accéléromètre, vM : bruit de mesure du magnétomètre, vG : bruit de mesure du gyromètre. L'orientation est estimée par rapport à un repère de référence, entièrement défini par la donnée des vecteurs Go et Ho. Par exemple, 5 le repère géocentrique est défini par les vecteurs Go (0;0; 1) etHo(0.5;0; 2 ) A des fins de simplicité, nous ne distinguerons pas les mesures selon les trois directions de l'espace. Comme cela apparaît clairement au vu de la définition 1 o mathématique des mesures (voir par exemple : Harold D. Black. A passive system for determining the attitude of a satellite. American Institute of Aeronautics and Astronautics Journal, 2(7):1350-1351, July 1964), chacune de ces mesures comporte respectivement une première partie « RG0 », « R.Ho » et w, qui contient l'information permettant d'obtenir une estimation 15 de l'orientation, une deuxième partie a, d et b qui représente les perturbations éventuelles pouvant apparaître, de façon aléatoire dans les mesures, et enfin une troisième partie VA ,VM, vG représentant le bruit de mesure au niveau de chaque capteur. On commence par définir ci-dessous les conventions de notation 20 utilisées dans la suite des développements . On utilise de manière préférentielle les quaternions pour la représentation des grandeurs. Il existe plusieurs notations de quaternions. La notation académique est celle de Hamilton. Nous utiliserons plutôt celle qui privilégie les multiplications de quaternions dans un ordre naturel. Elle a été proposée 25 par le JPL (Jet Propulsion Laboratory) en 1999 (voir par exemple : W. G. Breckenridge, "Quaternions - Proposed Standard Conventions," JPL, Tech. Rep. INTEROFFICE MEMORANDUM. 10M 343-79-1199, 1999 ; Nikolas Trawny et Stergios I. Roumeliotis, "Indirect Kalman Filter for 3D Attitude Estimation", A Tutorial for Quaternion Algebra, Department of 30 Computer Science & Engineering, University of Minnesota. Technical Report, 2005-002, Rev. 57, March 2005). Selon cette notation, le quaternion sera noté q et défini comme suit : = q4 +q1i+q2J+q3k partie partie réelle imaginaire En notation matricielle, on peut aussi écrire : _ T q1 q2 q3 q4 Si q et q4 remplissent les conditions ci-après, ce qui sera le cas pour tout notre exposé, q est dit quaternion unitaire. Il représente alors une rotation d'angle 0 autour du vecteur unitaire k : kx sin(O / 2) q = ky sin(8 / 2) = k sine / 2) et q4 = cos(8 / 2) (4) kz sin(O/2) Un vecteur p dans l'espace (3 coordonnées) peut être représenté sous forme de quaternion : q= q q4 (2) (3) =P 0 (5) Nous utiliserons les notations suivantes pour préciser les repères de référence dans lesquels les vecteurs s'expriment : LP=LC(q )GP (6) Où : - - q=cLq - GC(q) est une matrice de rotation (3x3) qui exprime le repère global {G} dans le repère local {L}. Avec cette notation, on peut passer du système de coordonnées du repère 15 local (repère capteur mobile) au global (repère fixe) et inversement par l'expression suivante : LP-G'q ' x L2 -~ (7) Nous manipulons des vitesses angulaires du système de coordonnées local {L} communément celui du capteur en mouvement par rapport au repère global {G}, la dérivée temporelle du quaternion s'exprime à partir de l'expression générale (voir Trawny et Roumeliotis, référence citée ci-dessus) : q (t) = lim 1 (L(t+o) q_L(t+At q (t)) (8) oto At Or: L(t+ot)--L(t+ot)-O L(t) Gq- L(t)q Gq (9) Où la rotation de {L(t) } à {L(t+dt)} s'exprime de la façon suivante: L(t+ot)- - L(t) - k sin(9/ 2) cos(9/2) k-6/2 1 z Se 1 At,- (10) -Le vecteur S9 est de même direction que l'axe de la rotation faisant passer 1 o de {L(t) } à {L(t+dt)} et a une magnitude égale à l'angle de rotation. Or, la définition de la vitesse angulaire nous donne : w=limse ot,o At Donc, d'après les expressions (8),(9),(10) et (11) : w 0 OL~G q (12) Nous aurons besoin d'intégrer un quaternion pour estimer les angles à partir des mesures co du gyro. Intégrer un quaternion est équivalent à résoudre 15 l'équation différentielle du 1 er ordre ci-dessus. Dans Trawny et Roumeliotis, (référence citée ci-dessus), on trouve la solution à l'ordre 0: /cos/ W l j bu -At + 1 sin/ _ At - SL(uo) q(tk ) 2 w 2 (13) 0 - û2 coy - (ox n((o)_ (14) Où (13) a une expression plus stable lorsque w-0 L- I col G q (tk+1 ) 4x4 + Ot n((o)1 G q (t ) (15) Notons que (13) peut se réécrire aussi : w . w -.sin - 1u) 2 L- Gq tk+l _ iW cos Ot 2 (16) où l'on voit que la solution à tk+1 n'est autre que la rotation du repère précédent pris à tk selon un axe défini par co avec un angle de rotation égal à Iwlot, ce qui correspond à l'hypothèse à l'ordre 0 qui veut que co est constant entre tk et tk+1. Un homme du métier de la capture de mouvement pourra consulter la littérature, notamment Trawny et Roumeliotis, (référence citée ci-dessus) pour établir les solutions aux ordres supérieurs (accélération angulaire constante, jerk angulaire constant, etc...). Le cadre de référence des notations permettant de décrire la présente invention est maintenant posé. Nous pouvons à présent décrire les étapes du procédé de traitement à l'aide des notations précédemment définies.15 La figure 2 représente une architecture fonctionnelle de mise en oeuvre du procédé de l'invention pour l'échantillon k, dans un de ses modes de réalisation. Avec les mêmes nommages de variables que dans la référence Black citée plus haut, on peut écrire les modèles de mesure pour les différents capteurs du dispositif considéré (deux capteurs de champs Go et Ho uniformes dans le repère absolu, par exemple un accéléromètre et un magnétomètre, et un capteur représentatif de la vitesse angulaire, par exemple un gyromètre) sous la forme : /Yc v w+b /c YA = g0(Go+a)0q-1 + vA (17) GYM/ ~q 0(Ho+h)0. -1/ wM/ A to, on peut considérer que le dispositif est dans sa position neutre ou de référence et qu'il est immobile, de sorte que seuls les champs Go et Ho sont mesurés, sans accélération a ni perturbation h. Le repère capteur se confond avec le repère fixe de référence et, si on néglige les bruits capteurs : /Yc (to )' / o YA(to) _ 'Go (18) c ~YMt°J H°J 15 Cette équation montre qu'on peut enregistrer par la mesure les champs de référence Go et Ho à to ou, si on le désire, à chaque fois que a et h sont nuls. A t>to, considérons que le dispositif a quitté sa position neutre et se met en mouvement. Le principe général des traitements selon l'invention est le suivant : 20 - Les mesures sont pré-traitées soit pour éliminer le biais (cas du gyromètre), soit pour détecter les perturbations dont elles sont éventuellement affectées ; ces pré-traitements sont similaires à ceux effectués selon l'art antérieur ; ils sont indiqués sur la figure 2 par les références 210, 220 et 230 ; 25 - Lorsque les mesures en sortie de l'accéléromètre et du magnétomètre (ou d'autres capteurs remplissant les mêmes fonctions, comme indiqué plus haut) ne sont pas perturbées,10 l'orientation est fournie par un algorithme classique représenté sur la figure 2 par la référence 240 ; - Lorsqu'une des mesures en sortie d'un de ces capteurs est perturbée, elle est remplacée en entrée de l'algorithme de calcul de l'orientation par une mesure reconstruite à partir du quaternion calculé par un algorithme référencé 250 sur la figure 2, auquel sont fournis la vitesse angulaire en sortie du gyromètre et le quaternion en sortie de l'algorithme de calcul 240 à l'instant k-1. 1 o Les mesures sont donc d'abord pré-traitées : - Le traitement 210 permet de filtrer le gyromètre pour éliminer son biais b. Ce traitement consiste à évaluer sa variance sur une fenêtre glissante d'une durée choisie T (par exemple 1 seconde). Si la variance reste en deça d'un certain seuil, on peut considérer que le gyromètre ne subit pas 15 de vitesse angulaire. On établit la moyenne sur la fenêtre, ce qui donne une estimée du biais b sur chaque axe de mesure. 210 s'exprime donc ainsi: sivar(YG(tk-T+1),...,YG(tk))aG YG(tk)=YG(tk)-(YG)T r l r 11 r l_ r l (19) si var(YG (tk-T+1 ),..., YG (tk )) i aG YG (tk ) = YG (tk ) on applique ce procédé préferentiellement à l'initialisation en veillant à ce que le mobile observe une période stable supérieure à T ; 20 néanmoins, le biais peut être calculé en permanence sachant qu'il ne sera rafraîchi que si les conditions de seuil (19) sont respectées; dans (19) , le test est effectué pour chaque axe du gyroscope, ce qui permet de corriger éventuellement un axe lorsqu'aucun mouvement rotatoire n'est observé par ledite axe ; néanmoins, on peut remplacer YG par sa 25 norme pour ce test de la variance ; En ce cas, il faut que le mobile observe une position stable absolue (complète) pour que les corrections s'appliquent (par conséquent sur tous les axes) ; on confondra parfois la notation tk avec la notation k pour désigner l'instant. 30 Le traitement 220 permet de tester les mesures de l'accéléromètre (ou d'un capteur équivalent de champ uniforme) ; on applique le test proprement dit de l'existence ou non d'une accélération propre. Pour ce faire, on compare la norme de la mesure YA,k à la norme du champ gravitationnel (pour rappel, on travaille en multiple de Go), on fait donc une comparaison par rapport à 1 : Si YA k -1 0 YA,k = YA si test230 > 0 YM,k YM (20) - Si les tests précédents 220 et/ou 230 sont négatifs, il y a des perturbations, on construit alors, via le filtre CNP (Construction Non 25 Perturbée) une mesure artificielle basée sur le ou les champs de référence et sur le quaternion estimé au même instant par le filtre 250, dont le fonctionnement est décrit plus loin. La construction de cette mesure non perturbée s'exprime comme suit (valable pour YA et YM) : Yk =' k OO Go 0'7k (21) L'équation ci-dessus s'applique pour Ak autant que pour Mk , avec le même 30 formalisme (en utilisant Ho au lieu de Go). Les mesures Ak et Mk , sont ensuite directement mises en entrée d'un algorithme de calcul traditionnel d'orientation de type TRIAD (TRI-axial Attitude Determination) par exemple (cf : htt~:ten. ikidia.orlvikiTriad Method ou publication 1964 de Harold Black déjà citée). Cet algorithme permet d'estimer un quaternion de rotation à partir de la mesure des deux directions de référence fournies par les données Ak et Mk . Mais tout autre algorithme permettant de calculer l'orientation d'un objet dans un repère de référence à partir de sa position par rapport aux vecteurs champ gravitationnel et champ 1 o magnétique terrestres (ou de tout autre champ uniforme) peut remplacer l'algorithme TRIAD. En particulier, dans le cas de plus de deux capteurs, un algorithme de type QUEST (Quaternion ESTimator), décrit notamment dans la publication http:l,www.dept.aoe.vt.edulhcdhalllcourseslaoe4140lattde.pdf et, en comparaison avec TRIAD dans 15 htt~:/;t.malcolmdshuster.comlcomm 1981 a J TRAD-QUEST..df Le quaternion qk résultant est la réponse de notre système pour fournir l'orientation du dispositif. Cette estimation qk est ensuite ré-introduite dans le système d'intégration 250. Ce traitement permet de considérer l'attitude estimée actuelle qk et de 20 mettre à jour à l'instant k+1 grâce à la mesure de vitesse angulaire du gyromètre et la rotation élémentaire qui en découle entre k et k+1. D'après (16), le nouveau quaternion est donné par : L c q tk+i = (22) Les exemples décrits ci-dessus sont donnés à titre d'illustration de modes de 25 réalisation de l'invention. Ils ne limitent en aucune manière le champ de l'invention qui est défini par les revendications qui suivent. REVENDICATIONS1. Procédé d'estimation de l'orientation d'un objet dans l'espace à un instant k utilisant au moins une mesure (yA, yM) d'au moins un champ physique sensiblement uniforme dans le temps et dans l'espace et des mesures de la vitesse de rotation (yG) dudit objet selon trois axes de l'espace, ledit procédé comprenant une première étape de détection (220, 230) de la présence de perturbations dans les mesures dudit au moins un champ physique et une deuxième étape (240) de calcul d'un quaternion (qk) d'orientation dudit objet dans l'espace déduit d'au 1 o moins une représentation des mesures dudit au moins un champ physique en sortie de la première étape, ledit procédé étant caractérisé en ce que, lorsque en sortie de ladite première étape est détectée une perturbation de ladite au moins une mesure (yA, yM), ladite mesure perturbée en sortie de la première étape est combinée en entrée de la 15 deuxième étape avec un quaternion d'orientation dudit objet calculé en sortie d'une troisième étape (250) d'intégration d'une sortie corrigée de la mesure de vitesse de rotation de l'objet à partir d'un quaternion d'orientation dudit objet en sortie de la deuxième étape à l'instant k-1 20 2. Procédé selon la revendication 1, caractérisé en ce que le au moins un champ physique est la gravité terrestre. 3. Procédé selon la revendication 1, caractérisé en ce que le au moins un champ physique est le champ magnétique terrestre. 4. Procédé selon la revendication 1, caractérisé en ce qu'il utilise des mesures d'au moins deux champs physiques, la gravité terrestre et le champ magnétique terrestre. 30 5. Procédé selon la revendication 1, caractérisé en ce qu'il utilise des mesures d'un seul champ physique et des valeurs d'un champ synthétique construites à partir dudit seul champ physique, ledit champ synthétique étant défini tel qu'il ait au moins une composante 25orthogonale audit champ physique et que son produit vectoriel avec ledit champ physique soit non nul. 6. Procédé selon la revendication 5, caractérisé en ce que ledit champ physique est la gravité GO et ledit champ synthétique est un vecteur constitué par une combinaison linéaire de ladite gravité GO et d' un vecteur non nul choisi dans un plan orthogonal à GO. 7. Procédé selon la revendication 5, caractérisé en ce que ledit champ 1 o physique est le champ magnétique terrestre HO et ledit champ synthétique est un vecteur constitué par une combinaison linéaire dudit champ magnétique terrestre HO et d'un vecteur non nul choisi dans un plan orthogonal à HO. 15 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce qu'il comprend en outre dans la première étape, une sous-étape de correction (210) de la vitesse de rotation (yG,k) consistant à retrancher aux mesures de ladite vitesse de rotation un biais moyen (bnoyen) déterminé lors d'une étape de position au repos et à lui substituer une 20 vitesse de rotation corrigée (ÿG,k ). 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce que la première étape de détection (220, 230) d'une perturbation dans les mesures dudit au moins un champ consiste à déterminer si la différence 25 de la norme des mesures (yA, yM) dudit au moins un champ et de leurs valeurs de référence (Go, Ho) est supérieure à une valeur de seuil choisie (aA, aM). 10. Procédé selon l'une des revendications 1 à 9, caractérisé en ce que la 30 deuxième étape de calcul d'un quaternion (qk) d'orientation dudit objet s'effectue par composition d'au moins deux mesures de champs physiques en sortie de la première étape et utilise un algorithme de type TRIAD. 11. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que la deuxième étape de calcul d'un quaternion (qk) d'orientation dudit objet s'effectue par composition de mesures de plus de deux mesures de champs physiques en sortie de la première étape et utilise un algorithme de type QUEST. 12. Procédé selon l'une des revendications 1 à 11, caractérisé en ce que la troisième étape d'intégration d'une sortie corrigée de la mesure de vitesse de l'objet à partir d'un quaternion d'orientation dudit objet en 1 o sortie de la deuxième étape à l'instant k-1 est effectuée par résolution de l'équation différentielle : L~c = 2 0 L~c 13. Procédé selon l'une des revendications 1 à 11, caractérisé en ce que le 15 quaternion estimé de l'objet en sortie de la deuxième étape est normalisé. 14. Centrale d'attitude d'un objet dans l'espace comprenant des moyens de mesure d'au moins un champ physique (yA, yM) et de la vitesse de 20 rotation (yG) dudit objet selon trois axes de l'espace, ladite centrale comprenant un premier module de prétraitement à l'instant k pour détecter la présence du au moins un champ physique (yA, yM) et un deuxième module de calcul d'un quaternion (qk) d'orientation d'au moins une représentation de l'orientation dudit objet calculée à partir 25 d'au moins une représentation des mesures dudit au moins un champ physique en sortie du premier module, ladite centrale étant caractérisée en ce que, lorsque en sortie dudit premier module est détectée une perturbation de ladite au moins une mesure (yA, yM), ladite mesure perturbée en sortie du premier module est remplacée en entrée du 30 deuxième module par une mesure combinée avec un quaternion d'orientation dudit objet calculé en sortie d'un troisième module d'intégration d'une sortie corrigée de la mesure de vitesse de l'objet à partir d'un quaternion d'orientation dudit objet en sortie du deuxième module à l'instant k-1.