Brevet FR2929395A1 — publié le 20091002

METHODE D'ESTIMATION D'ATTITUDE D'UN SENSEUR STELLAIRE La présente invention concerne les senseurs stellaires destinés à mesurer l'attitude d'un véhicule spatial à partir de mesures de directions d'étoiles. Les mesures effectuées par les senseurs stellaires sont en effet sujettes à des erreurs dont les sources sont multiples, et de nombreux dispositifs et méthodes de correction de ces erreurs de mesure existent à ce jour. Parmi les types d'erreurs existants, on compte : - L'erreur de bruit temporel, qui désigne une erreur décorrélée d'un io pas d'échantillonnage à un autre, également appelée � bruit blanc �, souvent représentée comme un bruit gaussien et due entre autres phénomènes au bruit de la chaîne électronique d'amplification/lecture, ainsi qu'au bruit photonique. - L'erreur � fausse étoile �, qui désigne une perturbation due à un 15 corps étranger dans le champ, soit vu comme une étoile, soit venant perturber la mesure d'une étoile, par exemple un proton, une particule énergétique... - Les erreurs spatiales, qui désignent les erreurs qui, pour une étoile donnée, dépendent de la position de cette étoile dans le champ de vue. On 20 distingue deux types d'erreurs spatiales : - L'erreur dite � pixel �, dont la longueur de corrélation n'excède pas un voire quelques pixels, - L'erreur dite � champ de vue �, dont les longueurs de corrélation sont de l'ordre de grandeur de quelques fractions du champ de vue. 25 Le document US 6 047 226 révèle ainsi une méthode de contrôle d'attitude d'un engin spatial comportant le fait de modéliser des bruits temporels et spatiaux et utiliser un filtre de Kalman se basant sur cette modélisation. Le document US 6 272 432 révèle quand à lui une méthode de correction d'erreur spatiale comportant le fait de comparer l'attitude estimée dans deux 30 référentiels distincts, l'un inertiel l'autre stellaire, appliquer un filtre de Kalman dont les états caractérisent les erreurs spatiales et dont l'une des matrices de recalage caractérise à la fois les erreurs de biais d'un gyroscope et les erreurs d'attitude, recaler les deux estimations d'attitude et minimiser l'erreur sur ces estimations. Ces méthodes sont complexes à mettre en oeuvre et nécessitent de lourds calculs. De plus, un certain nombre de missions ont des exigences plus sévères que celles auxquelles répondent les senseurs stellaires existants, et c'est dans ce cadre que la présente invention se propose de fournir une méthode simple et efficace d'amélioration de la performance des viseurs stellaires visant à rejeter des erreurs ponctuelles sur une étoile telles que les perturbations protons, à calibrer en vol des erreurs spatiales et à réduire le bruit mono-étoile temporel. L'invention se propose de pallier au moins un des inconvénients des io dispositifs de l'art antérieur et a pour but de proposer une méthode et un équipement permettant d'améliorer l'estimation de l'attitude d'un engin spatial muni d'un catalogue d'étoiles embarqué, par des moyens simples et faciles à mettre en oeuvre. Dans le processus qui mène de l'acquisition des images de paysages 15 d'étoiles par le senseur stellaire à la détermination de l'attitude de l'engin spatial sur lequel est installé ledit senseur stellaire, il y a deux étapes qui viennent en amont et en aval de l'utilisation des distances inter-étoiles conformément à la présente invention. En amont du calcul ou de la mesure des distances inter-étoiles, il est 20 nécessaire d'identifier le paysage d'étoiles observé. Les méthodes d'identification sont nombreuses et bien connues de l'homme du métier, et c'est pourquoi nous ne nous attarderons pas à les décrire ici. En aval de l'estimation des erreurs de mesure de position des étoiles, la détermination de l'attitude d'un engin spatial sur lequel est installé un senseur 25 stellaire à partir des directions des étoiles est un procédé connu de l'homme du métier qui comporte le calcul de l'orientation du senseur en déterminant les quaternions d'attitude. L'invention propose donc une méthode d'estimation d'attitude d'un senseur stellaire muni d'un catalogue d'étoiles embarqué, comportant le fait de mesurer en 30 axes relatifs au senseur stellaire des distances angulaires inter-étoiles, calculer lesdites distances à partir d'informations du catalogue d'étoiles embarqué, caractérisée en ce qu'on compare les distances angulaires inter-étoiles mesurées aux distances angulaires inter-étoiles calculées afin d'en déduire une information 15 20 25 30 sur la précision des mesures d'attitude du senseur, ainsi qu'un équipement mettant en oeuvre cette méthode. D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre, et en regard des 5 dessins annexés, donnés à titre d'exemples non limitatifs et sur lesquels : - la figure 1 représente une vue schématique d'un équipement pour engin spatial conforme à la présente invention, - la figure 2 représente un premier mode de réalisation de l'invention permettant de calibrer des erreurs spatiales, - la figure 3 représente une évolution temporelle d'une estimation de biais en radians pour 16 zones du champ, 1600 zones au total, 12 jours de simulation, 12 étoiles dans le champ, une fréquence du senseur stellaire 1 de 10Hz, - la figure 4 représente une évolution temporelle d'une amplitude d'une erreur pixel pour 16 zones du champ, 1600 zones au total, 12 jours de simulation, 12 étoiles dans le champ, une fréquence du senseur stellaire 1 de 10Hz, - la figure 5 représente un second mode de réalisation de l'invention permettant de détecter des biais mono-étoiles, - la figure 6 représente un tableau de performances d'estimation du biais, - la figure 7 représente un tableau de probabilités d'occurrence de rejet d'étoile conforme au second mode de réalisation de l'invention, - la figure 8 représente un autre tableau de probabilités d'occurrence de rejet d'étoile conforme au second mode de réalisation de l'invention, et - la figure 9 représente un dernier tableau de probabilités d'occurrence de rejet d'étoile conforme au second mode de réalisation de l'invention. Comme illustré sur la figure 1, la méthode concerne tout engin spatial sur lequel est installé au moins un équipement 0 comportant au moins un senseur stellaire 1, un ordinateur de bord 5 et un catalogue d'étoiles embarqué 2 qui contient la direction de certaines étoiles dans un repère inertiel. La précision sur les directions de ces étoiles étant très bonne relativement aux précisions de mesure des senseurs usuels, cette information est prise en référence. De cette information absolue de directions, il est possible de tirer des informations relatives de distances angulaires entre les étoiles, appelées io distances inter-étoiles. L'intérêt de cette information relative est qu'elle ne dépend pas de l'attitude 6 du senseur 1 : ainsi, quelle que soit l'attitude 6, la distance angulaire 3 entre deux étoiles mesurées en axes relatifs au senseur stellaire 1 est égale, aux erreurs de mesure près, à la distance angulaire 4 entre ces deux étoiles calculées 15 à partir des informations du catalogue embarqué. Ainsi qu'illustré sur la figure 1, pour chaque couple d'étoiles dont le senseur stellaire 1 mesure la distance inter-étoile, il est possible d'avoir accès à une information sur la précision de la mesure en comparant la distance inter-étoile 3 mesurée par le senseur stellaire 1 à la distance inter-étoile 4 calculée à partir des 20 informations du catalogue, cela via un ordinateur de bord 5. Selon un premier mode de réalisation, la méthode permet une calibration en vol des erreurs spatiales. La première étape 7 peut comporter le fait de modéliser les deux types 25 d'erreurs spatiales. Les sources principales des erreurs champ de vue sont des défauts optiques qui créent des aberrations induisant des biais sur les mesures de positions d'étoiles. Ces erreurs sont modélisées par un vecteur à deux dimensions dont les 30 composantes BX et By sont des fonctions lentement variables dans le champ et continues : �FOV (x, Y) - B y (x, y) La source principale des erreurs pixel est le calcul du barycentre de la tache étoile, erreur de l'ordre du centième de pixel. Ces erreurs sont modélisées par un vecteur à deux dimensions dont les composantes sont des sinus d'amplitude variable dans le champ, de période égale à un pixel, dont la phase est telle que l'erreur est nulle lorsque l'étoile est centrée sur un pixel : Épixel (x, Y) = Ax (x, y)sin 2n Ay (x, y)sin 2n où AX et Ay sont les amplitudes de l'erreur en x et y, px et py les tailles des io pixels en x et y. Ces amplitudes sont supposées lentement variables dans le champ. La seconde étape 8 peut comporter le fait de partitionner le champ de vue en un nombre z2 de zones élémentaires, les paramètres BX, By, AX et Ay prenant des valeurs différentes dans chacune des zones élémentaires. 15 Ce quadrillage permet de discrétiser le calcul qui conduit à l'estimation des erreurs spatiales. La valeur des paramètres ne dépend que de la zone dans laquelle l'étoile se situe. En tout point d'une zone élémentaire, chaque paramètre possède une même unique valeur propre à la zone en question, et ces valeurs sont différentes pour deux zones différentes. 20 Ainsi, pour caractériser complètement les erreurs spatiales il suffit d'identifier avec une bonne précision les z2 valeurs de BX, les z2 valeurs de By, les z2 valeurs de AX et les z2 valeurs de Ay, soit 4z2 paramètres à estimer. L'expérience a montré qu'une valeur de z réalisant un bon compromis entre le nombre de paramètres à estimer, la place mémoire nécessaire et la précision 25 du modèle est z = 40. La troisième étape 9 peut comporter le fait d'exprimer la différence entre les distances inter-étoiles 4 calculées à partir des informations du catalogue d'étoiles embarqué 2 et les distances inter-étoiles 3 mesurées par le senseur stellaire 1. 5 La notion de distance angulaire inter-étoile est exprimée par le produit scalaire entre deux mesures de direction étoile : Yi; La différence entre la mesure et l'expression de ce produit scalaire émanant des informations catalogue peut s'exprimer comme suit : cata log ue mes y i; = ui + L'. avec L' fonction des mesures et des phases pixel : mes u L' = - mes mes mes mes uzi uyj ùuyi uzj - T mes mes mes mes uzi uxj +uxi uzj mes mes mes mes uzi uyi ù uyj uzi mes mes mes mes uzj uxi + uxj uzi1 0 0 0 sincp 0 0 0 0 1 0 0 0 sincp. 0 0 0 0 1 0 0 0 smcpxi 0 0 0 0 1 0 0 0 smcp3 les valeurs Axk, Ayk, Bxk et Byk étant les valeurs des paramètres à estimer qui 10 correspondent à la zone occupée par l'étoile n°k, k allant de 1 au nombre d'étoiles n. La quatrième étape 10 comporte le fait d'utiliser un filtre d'estimation pour recaler un vecteur d'état formé par les amplitudes des erreurs spatiales à partir des distances inter-étoiles calculées 4 et des distances inter-étoiles mesurées 3. 15 Pour cela, un vecteur d'inconnues k formé des amplitudes Ax;, Ay;, Bx;, By;, des erreurs spatiales est constitué et la dérivée partielle des Y; par rapport à X est calculée : mes mes mes mes uzi uyj uyi uzj ) mes mes mes uxj +uxi uz; > mes mes mes mes uzj uyi uyj uzi ) mes mes mes -uZ. uxi mes mes) + uxj uzi mes mes mes mes) U zi uyj uyi uzj ) sin (P xi mes mes (ùuzi uxj mes u Zj mes mes (-uzj uxi mes mes +uxi uzj ) sin (pi mes mes +uxj uzi )sin(pyj mes uyi ùuy1 mes mes u zi ) sin (p xj T aY aX H Le filtre d'estimation proposé peut être un simple filtre de type moindres carrés récursifs, ou Kalman à gain constant. Le vecteur d'état Xk formé des erreurs spatiales est initialisé à 0 et est 5 recalé à chaque pas d'indice. L'équation de recalage au pas k, k allant de 1 à n(n-1)/2, est la suivante : Xk+l = Xk +KkHkYk avec Hk = axk , Ki, = PHk (HkPHk +R) ' Xk étant le vecteur d'état formé des erreurs spatiales, Io Hk la matrice d'observation, P la matrice de covariance de l'état estimé (constante), R la matrice de covariance des distances inter-étoiles, Yk le vecteur de mesure au pas d'indice k composé des Yii calculés à chaque pas d'indice. 15 Le vecteur d'état étant de taille importante, les matrices P et R sont réduites à des scalaires. La matrice R représente la covariance des distances inter-étoiles. Les bruits de mesure étant des bruits mono-étoile gaussiens, indépendants deux à deux et d'écart-type 6, la matrice R se calcule simplement : 20 VAR(Y) = E[(Y - E(Y))2 ] E 1 = E, �=G f reel(3 ` reel (1J `~2 +(u;''''(3) reel (2 `)2 +(u reel (1 ` n Rl; = 1L~u Du ) u J u2 u2 J J 2 reel (2 Llreel (3))2 + reel (3 ulreel (1 2 + ree2l (3 Ll lreel (2))2 + ree2l Ll (1 reel (3))2 + 1 2 2 l reel (2) ureel (3) 2 4 ree2l u (1) reel (1) + u lreel (2 uree2l (2 ) 2 l l En faisant la simplification : u (3) , u ; (3) 1 Ri; = 2.6 2 Ixi uxj) + uyi 11,1)2) La phase des erreurs pixels étant dépendante de la mesure de la position des étoiles par rapport aux pixels, un bruit étoile de 10" (16), qui correspond à environ 1/10 de pixel engendre une erreur de phase de l'ordre de 36°. Il est donc nécessaire de vérifier que la criticité de cette erreur de phase est minime. Des simulations ont ainsi été réalisées en testant différents jeux de réglage, faisant varier plusieurs paramètres dont la durée de la simulation, le io nombre de zones élémentaires et la constante de temps du simulateur. Des hypothèses simplificatrices ont été réalisées, le but étant de vérifier la faisabilité de la méthode. Ainsi, les bruits mono-étoiles sont supposés gaussiens, indépendants deux à deux et d'écart-type 6. La magnitude, la longueur d'onde et la température des étoiles n'ont ainsi pas été prises en compte. 15 Les simulations effectuent, pour un nombre donné de tirages de paysages d'étoiles : - un tirage aléatoire des positions des étoiles dans le champ, - un tirage gaussien de bruits mono-étoiles sur les étoiles considérées, - le calcul et la mesure des distances inter-étoiles, 20 - le recalage de l'état estimé selon les équations précédentes, - des statistiques (moyenne, écart type) des erreurs d'estimation, - des tracés temporels. On peut ainsi choisir, dans l'optique d'obtenir une convergence en 6 jours un réglage tel que le nombre de zones élémentaires soit égal à 40 x 40 = 1600 25 zones, pour un paysage de 12 étoiles dans le champ de vue mesurées à une fréquence de 10 Hz, ces deux derniers paramètres étant représentatifs du cas du senseur stellaire 1 SED 36. Les résultats de ces simulations sont représentés sur les figures 3 et 4. On y voit l'évolution temporelle sur 12 jours de l'estimation des biais et erreurs pixels 5 pour 16 zones du champ parmi 1600. Le terme de covariance P a été réglé à 0.1 �rad pour les termes de biais et d'amplitude de l'erreur pixel. Potentiellement, la performance d'estimation pourrait être encore améliorée en augmentant le temps de convergence, en baissant la valeur de la matrice de io covariance P du filtre par exemple. Le meilleur nombre de zones élémentaires résulte du compromis entre la vitesse de convergence de la méthode et la performance d'estimation. En effet, un nombre important de zones élémentaires garantit une très bonne performance d'estimation mais entraîne par ailleurs une longue durée de 15 convergence ainsi qu'un besoin important en place mémoire. Selon un second mode de réalisation illustré sur la figure 5, la méthode permet de détecter des erreurs de type � fausse étoile �, qui désignent une 20 perturbation due à un corps étranger dans le champ, soit vu comme une étoile, soit venant perturber la mesure d'une étoile, par exemple un proton ou une particule énergétique. Nous adopterons ici pour plus de clarté un formalisme mathématique dit � dans le plan � où les mesures étoiles sont représentées en deux dimensions 25 (x,y) ce qui n'est nullement limitatif par rapport à la méthode exposée. La première étape 14 est le calcul de la différence des carrés des distances inter-étoiles 3 mesurées par le senseur stellaire 1 et des distances inter-étoiles 4 calculées à partir des informations du catalogue d'étoiles embarqué 2 : Mu =dij,m - diu = (Xun -xjm)2+(im-Yjm) -(xi - Xi) - (yi-.yj)2, 30 où dij,m est la distance mesurée par le senseur stellaire, d;i la distance calculée à partir des informations du catalogue d'étoiles embarqué 2, (x;, y;) les coordonnées angulaires cartésiennes de l'étoile i répertoriées dans le catalogue d'étoiles io embarqué 2, (x;m, y;m) les coordonnées angulaires cartésiennes de l'étoile i mesurées dans le repère du senseur stellaire 1. En supposant les biais de mesures c faibles devant les coordonnées étoiles, un développement au premier ordre donne : Adu2 = 2(xi ùXi VY )(�xi -Ex] )+2(yi -yi )(� yi )+o(c ) La seconde étape 15 comporte le fait de calculer à partir de Ad;i2 l'erreur fausse étoile sur une étoile donnée indicée j. Cette erreur se traduit par un biais de mesure anormalement élevé par rapport au bruit habituellement attendu sur l'étoile j. io Soit f)-. = bxj le biais de mesure de l'étoile j à l'instant de mesure t. Les coordonnées mesurées de l'étoile j sont donc Xmj -xi +bxj Ymj ùYj +byi Les n -1 autres étoiles sont supposées non biaisées. Les biais de mesure commis sur les autres étoiles correspondent par conséquent au niveau de bruit 15 mono-étoile. Pour i≠j, nous avons donc : xmi =xi +Exi Ymi = Yi +E yi Soit Ad 2 ,,.2(xi ùXj)(�xi ùbxj+2(yi ù Y;X�yi ù by;)=2bxj(Xmj ùXmi)+2byi(Ymj ymi ) 20 A partir des n - 1 équations précédentes, un système conduisant au calcul du biais de mesure sur l'étoile j est obtenu : A.X. =Y;, avecA; =2 Xmj ù Xml Xmj ù Xm2 Ymj Yml Ymj Ym2 , A de taille (n-1,2), Xmj ù Xmn Ymj Ymn Il xj = b' vecteur d'inconnues de taille (2,1) et Y. = taille (n-1,1). Dans le cas où le nombre d'étoiles n mesurées est supérieur ou égal à 2, la résolution optimale classique par moindres carrés est possible par pseudo inverse, donnant le résultat suivant : =(AT*A)1AT Ce calcul engendre la simple inversion d'une matrice 2*2. Ainsi, en supposant que seule l'étoile n°j est biaisée et que ce biais est grand devant les bruits mono-étoile, il est possible de calculer la valeur du biais io sur l'étoile j grâce aux calculs des distances inter-étoiles. Afin de ne pas pondérer trop faiblement des mesures de distances inter-étoiles d'étoiles trop proches l'une de l'autre, il convient de normaliser les lignes de la matrice A en divisant chaque ligne du vecteur de mesure Y par la norme du vecteur Xmj Ymj Ymi 15 La troisième étape 16 comporte le fait de tester si tous les biais calculés sont inférieurs en norme à un seuil donné au dessus duquel une étoile est considérée comme étant biaisée. Une étoile est dite biaisée lorsqu'elle présente une erreur de mesure supérieure à une certaine limite, en pratique plus grande que la valeur 30" (36) du 20 bruit mono-étoile. Il faut avant de mener à bien la troisième étape déterminer le meilleur seuil de rejet d'étoiles présentant un biais déterminé. La valeur du biais à partir duquel une étoile est rejetée sera prise par exemple égal à 90", valeur correspondant à un pixel du senseur stellaire 1 SED16. Le meilleur seuil de rejet d'étoiles est un 25 seuil permettant d'atteindre les objectifs suivants : - Objectif 1 : Ne pas rejeter d'étoiles lorsqu'aucune étoile du champ n'est réellement biaisée. vecteur de mesure de - Objectif 2 : Détecter une étoile qui présente un biais supérieur à une valeur limite donnée de manière à rejeter sa mesure avec une probabilité suffisante. - Objectif 3 : Ne pas rejeter une autre étoile non biaisée lorsqu'une étoile biaisée est dans le champ (on parle alors de mauvaise détection). - Objectif 4 : Limiter les perturbations induites par une étoile faiblement biaisée, c'est-à-dire une étoile dont le biais est inférieur à 90". Pour atteindre l'objectif 1, la valeur du seuil de rejet d'une étoile doit être supérieure à la performance d'estimation des biais. La performance ao d'estimation des biais étant la précision avec laquelle un biais est déterminé, il est nécessaire de faire des simulations qui fournissent un ordre de grandeur de cette performance. Pour cela, on ne simule qu'une étoile biaisée parmi n. On simule un biais sur cette étoile de 180", ce qui correspond à environ deux pixels du senseur stellaire 1 SED16. On simule des bruits mono-étoile de valeur 30" soit 36 par étoile. On réalise alors 1000 tirages aléatoires de paysages d'étoiles et on étudie la sensibilité des résultats statistiques en fonction du nombre d'étoiles n. La figure 6 représente les performances 6o d'estimation du biais pour différentes valeurs de n. On voit ainsi que pour 10 étoiles, le biais détecté varie entre 180"-19.6" et 180"+19.6" sur l'axe des abscisses et 180"-18.7" et 180+18.7" sur l'axe des ordonnées. Soit une détection du biais en arcsecondes dans l'intervalle [160.4 ; 199.6] suivant l'axe des abscisses et [161.3 ; 198.7] suivant l'axe des ordonnées. Pour atteindre l'objectif 1, il faut donc déterminer la performance 6o d'estimation des biais et choisir une valeur de seuil de rejet supérieure à cette performance. Pour atteindre l'objectif 2, il est nécessaire de choisir une valeur inférieure à la valeur du biais maximal autorisé auquel on soustrait la performance 6o d'estimation des biais. Par exemple, si l'on veut rejeter des étoiles réellement biaisées de 90", il est alors nécessaire de fixer le seuil à une valeur inférieure à 90"-36o. Pour atteindre l'objectif 3, il est nécessaire de calculer l'impact d'un biais étoile sur l'estimation des n-1 autres biais étoiles. En effet, l'impact d'un biais étoile sur les estimations des autres biais est proportionnel au biais lui-même et dépend grandement de la configuration géométrique du paysage d'étoiles. Même s'il est statistiquement très rare de trouver des configurations défavorables rendant très difficile la détection de l'étoile biaisée, des simulations doivent être effectuées pour fixer la meilleure valeur du seuil de rejet. Les simulations effectuent donc pour plusieurs valeurs de seuil de rejet : - un tirage d'un grand nombre de paysages de 10 étoiles, - un biais fixe b = [bx, by] sur une seule étoile j parmi n à chaque tirage de paysage étoiles, - une estimation des n biais étoile. Plusieurs cas sont alors possibles : Cas 1 : Le biais maximum estimé correspond à l'étoile j et ce biais estimé est supérieur au seuil de rejet. Ce cas correspond au cas où la détection d'étoile a été bonne. - Cas 2 : Le biais maximum estimé ne correspond pas à l'étoile j et ce biais est supérieur au seuil de rejet. Ce cas correspond au cas où une mauvaise étoile a été rejetée. - Cas 3 : Tous les biais estimés sont inférieurs au seuil de rejet. Ce cas correspond au cas où aucune étoile n'a été rejetée. Les probabilités d'occurrence des cas 1, 2 et 3 sont ainsi calculées pour plusieurs valeurs de seuil de rejet et plusieurs valeurs de biais et sont représentées sur la figure 7. Pour atteindre l'objectif 4, il est de plus nécessaire de vérifier que des faibles biais ne perturbent pas la détection. Les probabilités d'occurrence des cas 1, 2 et 3 sont calculées pour de faibles valeurs de biais et plusieurs valeurs de seuil de rejet, et sont représentées sur la figure 8. La figure 9 représente en outre les probabilités d'occurrence des cas 1, 2 et 3 lorsqu'il n'y a aucun biais. Les simulations font ainsi apparaître que pour un bruit mono-étoile de 30" (36), en considérant qu'une étoile est biaisée lorsque son biais est supérieur à 30", un paysage de 10 étoiles dans le champ de vue, et un objectif qui est de rejeter efficacement des étoiles présentant un biais supérieur à 90", on constate qu'un choix de seuil de l'ordre de 30" induit trop de rejets intempestifs en présence de biais faibles ou en l'absence de biais. En revanche, des seuils de l'ordre de 50" à 70" permettent d'atteindre tous les objectifs fixés. Le choix du seuil doit donc se faire en fonction des étoiles que l'on souhaite rejeter. Si l'on souhaite rejeter efficacement des étoiles présentant un biais assez faible, de l'ordre de 60" par exemple, il faudra choisir un seuil de rejet inférieur à 50". Il est à noter que plus on baisse ce seuil, plus le risque de rejeter des étoiles non biaisées augmente, ce qui n'est pas un réel problème, puisque l'estimation des n-1 biais qui suit la rejection d'une étoile suite à l'estimation des n biais laisse suffisamment de chance d'éliminer l'étoile biaisée. Une fois le seuil choisi, la troisième étape 16 qui comporte le fait de tester si tous les biais calculés sont inférieurs en norme audit seuil au dessus duquel une étoile est considérée comme étant biaisée peut être menée à bien. La quatrième étape 17 comporte alors le fait de recommencer les étapes suivantes lorsqu'au moins une étoile présente un biais supérieur en norme audit seuil en ayant préalablement éliminé de l'ensemble fini d'étoiles l'étoile ayant le biais le plus grand en norme 18: - calcul des différences Ad2 des carrés des distances inter-étoiles mesurées 3 par le senseur stellaire 1 et des distances inter-étoiles calculées 4 à partir des informations du catalogue d'étoiles embarqué 2, - calcul à partir des Ad2 de l'erreur � fausse étoile � sur une étoile donnée, - tester si tous les biais calculés sont inférieurs en norme audit seuil au dessus duquel une étoile est considérée comme étant biaisée, - recommencer les étapes précédentes lorsqu'au moins une étoile présente un biais supérieur en norme audit seuil en ayant préalablement éliminé de l'ensemble fini d'étoiles l'étoile ayant le biais le plus grand en norme. Plusieurs étoiles biaisées peuvent ainsi être détectées. De plus, même si une étoile non biaisée est rejetée lors de la première boucle de calcul, il est possible d'estimer l'étoile réellement biaisée lors de la seconde boucle, ce qui garantit la détection d'une étoile biaisée. Le nombre de boucles à effectuer est en général limité à deux, voire trois lo dans des cas très rares. REVENDICATIONS1. Méthode (0) d'estimation d'attitude (6) d'un engin spatial muni d'un senseur stellaire (1) et d'un catalogue d'étoiles embarqué (2), comportant le fait de s mesurer en axes relatifs au senseur stellaire (1) des distances angulaires inter-étoiles (3), calculer lesdites distances (4) à partir d'informations du catalogue d'étoiles embarqué (2), caractérisée en ce qu'on compare les distances angulaires inter-étoiles mesurées (3) aux distances angulaires inter-étoiles calculées (4) afin d'en déduire une information sur la précision des mesures d'attitude (6) du 10 senseur (1). 2. Méthode selon la revendication 1, caractérisée en ce que, pour la calibration des erreurs spatiales, elle comporte les étapes (7, 8, 9, 10) consistant à: 15 - modéliser des erreurs spatiales à l'échelle du pixel et du champ de vue du senseur stellaire (1), - partitionner le champ de vue en zones élémentaires, à chaque zone étant associé une valeur d'erreur spatiale à l'échelle du champ de vue et une valeur 20 d'amplitude d'erreur spatiale à l'échelle du pixel, - exprimer la différence des carrés des distances inter-étoiles (3) mesurées par le senseur stellaire (1) et des distances inter-étoiles (4) calculées à partir des informations du catalogue d'étoiles embarqué (2), - utiliser un filtre d'estimation pour recaler un vecteur d'état formé par les 25 amplitudes des erreurs spatiales à partir de la différence des carrés des distances inter-étoiles (4) calculées et des distances inter-étoiles (3) mesurées, les trois premières étapes (7, 8, 9) étant permutables. 30 3. Méthode selon la revendication 2, caractérisée en ce que les erreurs spatiales à l'échelle du champ de vue sont de la forme : �FOV (x, Y) = où Bx et By sont des fonctions continues. 4. Méthode selon au moins l'une des revendications 2 ou 3, caractérisée en w ce que les erreurs spatiales à l'échelle du pixel sont de la forme : Ax (x, y) sin(2n x ) px Ay (x, y) sin(2Tr Y ) py où Ax et Ay sont les amplitudes de l'erreur en x et y et px , py sont les tailles en x et y des pixels, x et y représentant la position des étoiles dans le pixel. 15 5. Méthode selon au moins l'une des revendications 2 à 4, caractérisée en ce que le filtre d'estimation (10) utilisé est un filtre de type moindres carrés récursifs ou Kalman à gain constant. 6. Méthode selon au moins l'une des revendications 2 à 5, caractérisée en ce 20 que l'équation de recalage au pas k, k allant de 1 au nombre de distances inter-étoiles, est : Xk+l =Xk +KkHkYk avec Hk = aYk ,Kk =PHkT(HkPHkT +R) 1 où Xk est le vecteur d'état formé par les erreurs spatiales, Yk la matrice composée des différences des carrés des distances inter-25 étoiles mesurées (3) et calculées (4), Hk la matrice d'observation, Bx (x,Y) By(x, y) �pixel (x,Y) =P la matrice de covariance de l'état estimé, R la matrice de covariance des distances inter-étoiles, Kk la matrice de Kalman. 7. Méthode selon la revendication 1, caractérisée en ce que, pour la détection des biais ponctuels mono-étoile, elle comporte les étapes (14, 15, 16, 17, 18, 20) consistant à : to - calculer pour un ensemble fini d'étoiles la différence des carrés des distances inter-étoiles mesurées (3) par le senseur stellaire (1) et des distances inter-étoiles (4) calculées à partir des informations du catalogue d'étoiles embarqué (2), - calculer à partir de la différence des carrés desdites distances calculées (4) 15 et mesurées (3) le biais ponctuel mono-étoile sur chaque étoile dudit ensemble fini d'étoiles, tester si tous les biais calculés sont inférieurs en norme à un seuil donné au dessus duquel une étoile est considérée comme étant biaisée, recommencer les étapes ci-dessus lorsqu'au moins une étoile présente un 20 biais supérieur en norme au seuil en ayant préalablement éliminé de l'ensemble fini d'étoiles l'étoile ayant le biais le plus grand en norme. - calculer l'attitude (6) du senseur (2) lorsque tous les biais calculés sont inférieurs en norme au seuil. 25 8. Méthode selon la revendication 7, caractérisée en ce que le seuil est tel que la probabilité de rejeter une étoile non biaisée est faible. 9. Méthode selon au moins l'une des revendications 7 et 8, caractérisée en ce que le seuil est tel que l'impact d'un biais étoile sur l'estimation des autres biais 30 étoiles est faible.5 10. Méthode selon l'une au moins des revendications 7 à 9, caractérisée en ce que le seuil est tel que des faibles biais ne perturbent pas la détection d'une étoile biaisée. 11. Méthode selon l'une au moins des revendications 7 à 10, caractérisée en ce que la valeur du seuil pour rejeter des étoiles présentant un biais supérieur à 90" avec une valeur de bruit mono-étoile de 30" est comprise entre 50" et 70". 12. Equipement (0) pour engin spatial, caractérisé en ce qu'il comporte des 10 moyens aptes à mettre en oeuvre une méthode d'estimation d'attitude (6) conforme à au moins une des revendications précédentes. 13. Engin spatial caractérisé en ce qu'il comporte au moins un équipement (0) conforme à la revendication 12. 15