FR 2664044 A1 19920103 FR 9008230 A 19900629 Procédé et dispositif de détermination d'une orientation liée à un système mobile, notamment de la- ligne de visée dans un viseur de casque La présente invention se rapporte à la détermination d'une orientation liée à un système mobile, par exemple à un viseur de casque et a plus particulièrement pour objet un procédé, et le dispositif correspondant, de détermination d'une orientation liée à un système mobile permettant d'obtenir une grande précision. Le problème consiste à repérer avec précision la direction d'un axe orienté sans avoir de contact matériel avec lui, en particulier dans les systèmes viseurs de casque utilisés dans les avions de combat ou les hélicoptères. Le fonctionnement d'un tel système est rappelé brièvement ci-après : par un dispositif collimateur annexe, le pilote voit, à travers une glace semi-réfléchissante solidaire du casque et interposée sur un axe de vision, l'image d'un réticule projeté à l'infini en superposition avec la scène extérieure. Lorsqu'il veut désigner une cible, le pilote amène ce réticule à coïncider avec la cible et signale que la coïncidence est réalisée au moyen par exemple d'une commande par bouton poussoir prévue à cet effet. A condition de repérer, au moment de la signalisation, la position exacte du casque, il est possible de déterminer par rapport à l'avion, la direction de visée et de désigner l'objectif à un système d'arme, ou de pointer dans cette direction un système optique ou autre. On connaît notamment par le brevet français nO 79 14441 un dispositif de mesure d'orientation, notamment pour système viseur de casque, dans lequel, pour ramener la ligne de visée déterminée par la position du casque du pilote dans un repère lié à l'aéronef, un capteur de position constitué de trois bobines électromagnétiques orthogonales est placé sur le casque; par affineurs, en un point fixe de la cabine, sont placées trois autres bobines électromagnétiques, formant un émetteur ou "radiateur" selon la terminologie en usage dans le domaine. Le procédé consiste alors à faire passer successivement, dans chaque bobine de l'émetteur (formant le trièdre fixe), un courant électrique.Ces courants font apparaître successivement, trois champs magnétiques y associés qui sont captés par les bobines du capteur (formant le trièdre mobile lié au casque). Pour chaque champ magnétique émis les trois composantes orthogonales correspondantes sont captées dans les bobines du capteur. Les neuf composantes résultantes permettent alors de trouver la transformation (rotation) permettant de passer du trièdre mobile au trièdre fixe. En effet, ces neuf composantes dépendent - de la position du capteur dans l'espace - de l'orientation du capteur. En fait, la transformation résultant d'une telle méthode peut être assez différente de la transformation des coordonnées permettant de passer effectivement du trièdre de référence mobile au trièdre de référence fixe par rapport à l'aéronef, notamment du fait que l'on utilise des champs magnétiques. En effet, ces champs peuvent être notablement perturbés du fait de tous les éléments métalliques de l'environnement, qui créent des réflexions d'ondes magnétiques dans la cabine. De plus, des défauts inhérents à la réalisation des bobines d'émission ou de réception magnétiques créent également des perturbations. La détermination d'une direction en utilisant les lois de l'électromagnétisme en espace libre n'est donc pas possible quand l'espace considéré est la cabine de pilotage d'un aéronef. C'est pourquoi certains dispositifs connus prévoient, avant chaque mesure, la mise en mémoire de mesures préalables et la comparaison de mesures effectuées à un instant donné avec les mesures mises en mémoire afin de déterminer la rotation qui fait passer du trièdre fixe au trièdre mobile. Une telle méthode nécessite la mise en mémoire d'une quantité très importante de mesures qui doivent être reprises intégralement dès que des paramètres tels que, par exemple, la position exacte de l'émet- teur mobile par rapport à la ligne de visée, changent. L'invention a pour objet un procédé, et le dispositif correspondant, de détermination d'une orientation dans un système mobile permettant d'utiliser en fonctionnement le même procédé de mesure que précédemment, mais en évitant les défauts sus-mentionnés. Pour cela, dans une phase préalable, une modélisation prenant en compte les perturbations notamment magnétiques dans la cabine dans laquelle est placé le système de détection, est établie à partir d'orientations particulières connues du capteur mobile ou à partir d'orientations particulières connues de la ligne de visée et pour chaque orientation à partir d'un ensemble de positions de m, origine du repère constitué par le capteur mobile. Cette phase permet d'approcher à partir de mesures réelles les relations complexes entre position et orientation dues aux réflexions des champs magnétiques et aux limites physiques de l'émetteur et du capteur et conduit à la détermination préalable d'une fonction de modélisation globale tenant compte des effets perturbateurs.Cette phase préalable, éventuellement complétée par une initialisation permet alors en fonctionnement de déterminer directement l'orientation associée à une position quelconque du système mobile, à partir des mesures expérimentales réelles. Selon l'invention, un procédé de détermination d'une orientation liée à un système mobile dans un espace déterminé, utilisant un système magnétique de détection de position et d'orientation comportant un émetteur formé de trois bobines orthogonales ayant une position fixe dans l'espace déterminé et un capteur formé de trois bobines orthogonales fixé au système mobile, r étant la distance entre l'origine du repère émetteur et l'origine du repère capteur, un cycle de détection consistant à émettre successivement dans les trois bobines de l'émetteur et à mesurer les composantes de champ générées dans les bobines du récepteur, est caractérisé en ce qu'une phase préalable de modélisation de l'espace déterminé consiste à effectuer des mesures des champs magnétiques dans le capteur pour un ensemble de positions connues du système mobile par rapport à un repère dit repère cartographique, des paramètres de modélisation déduits des mesures de référence en chaque point étant calculés pour en déduire la rotation relative du capteur dans le repère cartographique des mesures de référence, et en ce que, en cours d'utilisation, l'orientation courante du système mobile est déterminée directement à partir des paramètres de modélisation constituant une fonction compacte de modélisation position des effets perturbateurs, fonction de la position. L'invention a également pour objet le dispositif destiné à la mise en oeuvre de ce procédé. L'invention sera mieux comprise et d'autres caractéristiques apparaîtront à l'aide de la description qui suit en référence aux figures annexées. - La figure 1 est un schéma synoptique du dispositif selon l'invention. - La figure 2 illustre les différentes étapes du procédé selon l'invention. La figure 1 représente schématiquement un émetteur 1 qui comporte trois bobines disposées selon un trièdre trirectangle cet émetteur est placé en un point fixe d'un espace donné et, dans l'exemple décrit, cet espace est la cabine de pilotage d'un aéronef ; une unité centrale 10, qui fait partie d'un ensemble de matériel électronique embarqué, délivre un signal électrique d'excitation successivement aux trois bobines de l'émetteur 1, au cours de cycles successifs, via un dispositif amplificateur 20 sur commande à partir d'un panneau de commande 30. Un capteur 2, disposé en phase de vol sur le casque du pilote de l'avion et formé également de trois bobines disposées selon un trièdre trirectangle, est également relié à l'unité centrale 10 via un module de connexion 40 ; ce capteur peut changer de place dans l'espace donné, en phase de vol selon les mouvements de tête du pilote et en phase préalable via un robot, pour occuper des positions prédéterminées. L'unité centrale 10 comporte un organe de mesure qui reçoit du capteur 2, au cours de chaque cycle de mesure, neuf signaux, en 3 fois 3 signaux, qui sont les signaux induits par l'alimentation successive de chacune des bobines de l'émetteur 1, dans chacune des bobines du capteur 2. Les amplitudes de ces neuf signaux constituent les neuf valeurs d'une matrice carrée de dimension 3x3.Le problème est alors de déterminer, à l'aide des lois de l'électromagnétisme, la transformation, i. e. la rotation, qui permet de passer du repère mobile formé par le capteur 2 au repere avion. Par convention l'origine du repère formé par l'émetteur 1 est le point 0, la distance Om de l'émetteur 1 au capteur 2 est appelée r, avec OmOm = u ou uf est le vecteur unité dans la direction Om avec pour projections x, y, z dans le trièdre fixe. Si au point m origine du repère formé par le capteur une première mesure est effectuée, une matrice de neuf composantes P est obtenue ; si une seconde mesure est effectuée, avec une orientation différente du capteur et la même origine m, une seconde matrice M de neuf composantes est obtenue ; la matrice M peut s'écrire M = R.P. où R est une matrice qui correspond à une rotation pure, si bien que tR = R-1 ce qui exprime que la matrice transposée t R d'une matrice rotation est égale à son inverse R ; donc le produit de la transposée tM de M par M s'écrit tMM: tP.tR.R.P = tpR i'R,P = tp p Le produit d'une matrice de mesure transposée avec la matrice elle-même est donc en principe invariant en un point m donné, quelle que soit la rotation du repère ; la matrice produit est de plus symétrique. Théoriquement P et tM.M peuvent être exprimés en fonction de r, x, y, z à l'aide des équations de Maxwell caractérisant les lois de propagation électromagnétique ; en réalité ce n'est pas possible car l'émetteur et le capteur ne sont pas des dipôles parfaits, du fait que l'espace de la cabine n'est pas un espace libre, et que le casque porte des masses perturbatrices mobiles en même temps que lui. Comme indiqué ci-dessus, pour résoudre ce problème, le procédé selon l'invention consiste à établir une modélisation préalable à partir de mesures réelles effectuées pour des positions et orientations prédéterminées - connues du capteur (ou d'une ligne de visée qui lui est liée). Pour cela, dans une phase préalable, l'unité centrale commande des mesures systématiques, effectuées, dans une première étape, en conservant toujours la même orientation à un capteur, 2, et cela selon un maillage constitué par un ensemble de points m de la cabine où peut se trouver le capteur 2 lorsqu'il est placé sur le casque du pilote à son poste de pilotage. Un ensemble de matrices de mesures de dimensions 3 sur 3 associé à l'ensemble des points m est donc obtenu. Ces matrices sont transmises, via un circuit d'interface 50, à la mémoire d'un calculateur 60 dit de modélisation en temps différé pour y être stockées et traitées. Ces mesures peuvent être effectuées avec le récepteur seul, dit alors capteur de cartographie magnétique, déplacé par un robot, c'est-à-dire sans le casque et le pilote.Ces mesures permettent donc d'établir une "cartographie" magnétique du volume dans lequel peut se trouver le capteur 2, pour une orientation particulière de ce capteur. De préférence ce relevé est effectué selon un maillage régulier dans les trois axes, avec un pas de 3 à 4cm. Ensuite d'autres mesures sont effectuées, toujours avec le capteur de cartographie seul, mais avec d'autres orientations mécaniquement connues de ce capteur, en plusieurs points régulièrement espacés. Si P est identifiée à la matrice de mesure obtenue lorsque les axes des bobines du capteur coïncident avec le repère de la cartographie, c'est-à-dire le repère fixe du robot directement lié au repère avion, une fonction globale reliant les composan tes de P aux composantes de tMM est recherchée à partir des matrices de mesures correspondant à la première orientation. Cette fonction traduit donc les effets perturbateurs liés au déplacement du capteur dans la cabine, le capteur étant maintenu dans une orientation fixe. Pour cela, dans un mode de réalisation, cette étape du procédé consiste à rechercher 9 relations polynomiales optimum liant chaque composante de la matrice P (à neuf composantes) aux 6 composantes de la matrice symétrique tMM. Les coefficients des monômes formant les 9 polynômes constituent la fonction F recherchée. La résolution est p facilitée par l'homogénéisation suivante r4P = F (r7 .tMM), où r = Om.En effet, les champs magnétiques décroissant en i un produit matriciel du type tMM décroît en ; pour rechercher une fonction reliant P à tMM, il est intéressant d'homogénéiser en dimension les deux entités à relier, et de leur trouver une dimension commune favorable. La dimension la plus favorable est la distance r. La meilleure approximation de r est donnée par la valeur de la racine 6e de la trace de la matrice tMM, résultant du calcul dans les conditions dipolaires. Après homogénéisation comme indiqué ci-dessus, on identifie chaque composante de la matrice r4P à un polynôme fonction des termes de la matrice r7 tMM, par régression linéaire sur les 6 termes de r7 tM. M en prenant pour régresseurs ek tous les produits de degré au plus égal à 3 calculés en un point m à partir de la matrice de mesure M en ce point. Les produits de degré 1 sont les six termes r7 tMM, les produits de degré 2 sont tous les produits de deux termes de tMM, y compris les termes au carré etc...Jusqu'à tordre 3, 7 pour les 6 termes de r tMM résultant d'une matrice de mesure M au point m, il y a k = 84 régresseurs ek Si n mesures donnant chacune une matrice M ont été établies en n points m, chaque composante de r4 (P), Pj (j = 1 à 9) est identifiée à n polynômes fonction chacun d'un ensemble de régresseurs différent (E) i = (e1, e2, ... ek)i, i = 1 à n calculés à partir des composantes de r7 tMM, et des coefficients (a0, a1, ... ak)j recherchés pour cette composante Pj = a0 + +ae1 + ... akek. Les coefficients sont obtenus en résolvant, au sens des moindres carrés selon la méthode dite de GRAM SCHMIDT, les 9 systèmes indépendants, P = E. A, où E est la matrice ligne définie par les k régresseurs calculés pour chaque point m de la cartographie (E a k colonnes et n lignes). A est la matrice colonne formée des coefficients (a0 ... ak)j de j = 1 à 9 recherchés, et P est la matrice colonne constituée de la valeur mesurée de cette composante Pj répétée n fois. L'exploitation directe des 9 ensembles de coefficients (a0 . . . ak) ... (a0 ... ak) 9 conduit à une première approximation Mc de la matrice cartographie P et donc des trois angles G1 (gisement), S1 (site) et R1 (roulis) caractéristiques de la rotation du capteur par rapport au repère de la cartographie, par la relation M = R.Mc. Puis dans une étape suivante, à partir des mesures effectuées en plusieurs points dans diverses orientations du capteur, une fonction globale F est recherchée telle que (G2, S2, R2) = Fr(G1, S1 R1, MM), (G2, S2, R2) étant des angles vrais. Cette fonction globale traduit les effets perturbateurs induits par la rotation du capteur. Cette étape consiste donc à rechercher 3 relations polynomiales optimum entre d'une part chaque angle vrai (G2, S2, R2) et d'autre part les trois angles approximés correspondants résultant de la mesure (G1, S1, R1) calculés précédemment et les différentes positions, en rotation, du repère capteur, c'est-à-dire les 6 composantes de r7 tM.M. La méthode utilisée est analogue à celle utilisée pour rechercher Fp, c'est-à-dire que l'on procède par régression 7 linéaire sur les 6 termes de r tMM et les 3 termes G1, S1, R en prenant pour regresseurs tous les produits de degré au plus égal à 3, c'est-à-dire une constante et les produits de degré 1, G1, S1, R1 et les 6 termes de tMM, les produits de degré 2 qui sont les précédents élevés au carré et les termes croisés de type G1S1, GlR1 etc... (G1) (tMM)ij ... (R1) ( MM)ij etc. Le degré étant limité à 3, on est conduit à calculer 220 régresseurs. L'exploitation directe des 3 ensembles de coefficients a. formant les paramètres de modélisation issus de cet algorithme de modélisation permet alors d'obtenir, avec une très grande précision, l'orientation du capteur dans le repère de la cartographie, et cela à partir des angles approximés dans une phase précédente, G1, S1, R1. Une fois effectuées les opérations de cartographie décrites ci-dessus dans la cabine du porteur, et après que les mesures de cartographie aient été mises en mémoire dans la mémoire de l'ensemble 60 calculateur de modélisation en temps différé, les paramètres de modélisation issus des algorithmes de modélisation décrits ci-dessus sont mis en mémoire dans l'unité centrale embarquée 10. Puis des mesures sont effectuées après que le casque ait été réglé (réglage de l'optique, réglage de sa fixation sur la tête du pilote etc...). Si M est alors une matrice de mesures correspondant aux composantes de champ magnétique détectées dans le repère formé par le capteur fixé au casque lors de l'émission successive des trois bobines de l'émetteur fixe pour une position du casque et donc du repère mobile, le produit tM. M est calculé et la mesure des angles du capteur dans le repère de cartographie, Mc, peut alors être calculée directement La rotation R du capteur lié au casque par rapport au capteur de cartographie est alors déterminée par R = McM 1. La position du capteur étant alors fixée de manière précise, il est possible d'effectuer rapidement, si l'on connaît la matrice de rotation I de la ligne de visée par rapport au capteur, la transformation permettant de passer du repère capteur au repère cartographique lié au repère avion, notamment pour donner à un système de tir la direction de tir déterminée à partir du viseur de casque. L'équation de la ligne de visée Lv dans le repère lié au porteur est alors obtenue de la manière suivante : Lv = RI 1 où R = McM 1, soit Lv = McM 1I 1 Mais en règle générale une étape d'initialisation de la ligne de visée ou alignement initial, doit être prévue, car les diverses possibilités de réglage offertes à l'utilisateur ne permettent pas de connaître a priori l'orientation de la ligne de visée dans le repère du capteur sur le casque. Pour cela une mesure capteur est acquise par le dispositif quand le pilote vise une direction prédéfinie liée au porteur, par exemple quand le pilote vise un réticule colimaté.Soit Lvo cette direction prédéfinie dite d'initialisation, et soit Mo la mesure capteur lors de cette phase d'initialisation ; la mesure correspondante Mco dans le repère de cartographie est calculée et on écrit la matrice d'initialisation : I 1 = McoMo1. De plus, lorsque le repère de cartographie n'est pas celui du porteur, il est nécessaire d'effectuer un changement de repère, dit "harmonisation", caractérisé par une matrice d'harmonisation H. et on a alors Lv = HRi H Lvo ou encore Lv = H(McM 1) (MoMcolH 1Lvo, Les termes R et I devraient pouvoir être approximés avec précision à l'aide de la fonction F constituée des paramètres de modélisation décrite précédemment. En pratique l'application de cette fonction à cette étape du procédé peut s 'avérer inefficace du fait des éléments perturbants liés au casque, comme par exemple les tubes cathodiques. En conséquence, une étape de compensation des éléments perturbants liés au casque est également mise en oeuvre, selon la même approche que celle décrite ci-dessus (pour tenir compte des éléments perturbants liés à la cabine ou dus aux défauts de l'émetteur ou du capteur). Pour cela, un ensemble de mesures est effectué, avec le casque, en plusieurs points régulièrement espacés et dans diverses orientations, mécaniquement connues, du casque, l'une de ces orientations coïncidant avec l'orientation approximative de la direction de visée dite d'initialisation Lvo. L'application de la première fonction Fp permet de déterminer une première approximation de I et de chaque rotation capteur associée aux diverses orientations, et ainsi permet de définir une première approximation de l'orientation de la ligne de visée par (G1, Si, Roi). A ce niveau de l'algorithme, le résultat, bien qu'imprécis est déjà fortement indépendant de la valeur de I de l'initialisation en utilisation réelle. En effet le produit RI 1, quel que soit le montage du capteur sur le casque est toujours lié à la ligne de visée. La valeur de I trouvée lors de la cartographie correspond à un réglage particulier du casque. Une fonction F r tenant compte de la présence du casque et modélisant directement la ligne de visée théorique (G2, S2, R2) à partir de Gt, S1, R1 et r7 tMM est alors déterminée. La figure 2 illustre de manière synthétique les différentes étapes de la phase de modélisation du procédé de détermination d'orientation selon l'invention dans sa première version indiquée ci-dessus - La première étape est l'étape d'acquisition de matrices de mesure M1... M n pour un ensemble de n positions de m, pour une orientation donnée du capteur - La seconde étape 12 est l'identification de la matrice P à la matrice de mesure (Mc) correspondant à la position du capteur coïncidant avec le repère de la cartographie - La troisième étape 13 est le calcul de r, puis de la matrice r4(P) comme indiqué ci-dessus - La quatrième étape est l'étape de calcul pour toutes les matrices de mesures Mi de (r7 tMiMi) = r7(tMi)(Mi) pour i = 1 àn - La cinquième étape 15 est l'étape de calcul des n ensembles de régresseurs (E.) = (e1 ek)i pour i = 1-à n - La sixième étape 16 est l'étape de résolution des 9 systèmes d'équation P = E.A au sens des moindres carrés, pour déterminer les ensembles de coefficients (aO ... ak) ... (a0 ... 8k)9 formant la fonction de modélisation F p - La septième étape 17 permet de trouver à partir de cette fonction de modélisation, une première approximation de la matrice de cartographie Mc et donc de la rotation capteur rapport au repère de la cartographie [G1, S1, R1i. Cette étape termine le niveau mesure initial. Puis, pour le niveau dit "angulaire", une étape 18 d'acquisition de mesures pour diverses orientations connues du capteur est effectuée. Une étape de calcul 19, permet alors à partir de ces mesures, de rechercher la fonction globale F Une variante du procédé selon l'invention permettant d'éviter la première étape d'acquisition de mesures (capteur dans une orientation fixe déplacé seul) consiste à calculer, pour chaque mesure avec casque, une matrice de mesure virtuelle Mv telle que : Lv = HMvM 1 I 1 H 1 Lvo I peut bien sûr être déterminé à partir d'une cartographie avec capteur en orientation fixe et de la fonction Fp correspondante. Mais il est préférable de le déterminer par une méthode itérative, comme décrit ci-après. A chaque phase de l'itération, la matrice de rotation I de (G, S, R) choisie conduit à une série de mesures virtuelles Mv à partir desquelles est calculée une fonction E'p, par optimisation, pour que MvM ILvo soit le plus proche possible de Lv, par la relation r4 Mv = F' (r7 tMM). p A l'issue de l'itération, on obtient la valeur de I et simultanément les mesures virtuelles Mv associées. Les mesures virtuelles apportent un gain par rapport aux mesures précédentes car elles gomment la distorsion liée à la rotation du capteur et intègrent les déformations dues aux masses perturbantes liées au casque. De plus la recherche itérative de I permet d'éviter la phase de mesures initiales à partir du capteur seul en orientation fixe. Enfin la fonction F'p est également optimisée. Le temps de calcul pour la modélisation est évidemment un peu plus long, tout en restant acceptable, du fait de l'optimisation de I. La dernière étape reste alors le calcul de la fonction F r modélisant directement la ligne de visée théorique à partir de la meilleure fonction F'p, qui est effectuée à partir des mêmes mesures. Le procédé est également adapté à une correction directe des distorsions optiques liées au pare-brise (ou verrière) du porteur, lors de la visée. Pour cela, ces distorsions sont introduites lors de la détermination de Fr, soit par calcul, soit par visées réelles test, automatiquement ou sur commande par le pilote. Un avantage du procédé selon l'invention est que la détermination préalable d'une fonction compacte de modélisation permet de tenir compte des aberrations d'un système magnétique de détection de position et d'orientation lorsqu'un tel système est utilisé dans un environnement métallique perturbateur, type cabine d'hélicoptère. Cette fonction compacte de modélisation est déterminée à partir de mesures expérimentales effectuées pour des positions connues du repère capteur ou directement de la ligne de visée. Ensuite, en cours de vol, l'exploitation du système en temps réel est rapide, conduit à des mesures fiables, et ne nécessite que peu d'espace mémoire dans le système embar qué . L'invention n'est pas limitée aux modes de réalisation précisément décrits ci-dessus. En particulier, au niveau mesure qui aboutit à I'obtention des paramètres de modélisation des mesures de référence, l'identification des composantes de la mesure à un polynôme a été décrite comme étant effectuée par 7 régression linéaire sur les 6 termes de r Pour améliorer cette étape, il est possible d'effectuer une régression linéaire sur un nombre supérieur de termes, par exemple 9, en ajoutant aux 6 termes de r7 tMM, les 3 produits scalaires en ligne de la matrice M, ces termes étant significatifs de la rotation, ou encore 3 termes ,ui définis par r4 (M) .. = r7,ui(tMM)ii, c'est-à-dire des termes reliant les termes diagonaux de M aux termes diagonaux de tMM.Dans ces deux cas la fonction Fp est alors déterminée à partir de 220 régresseurs. Une autre possibilité consiste à utiliser les termes ,u décrits ci-dessus, sauf pour les produits où n'apparaissent aucun terme de la matrice tMM. Pour ces produits, les lli sont remplacés par les produits scalaires en ligne de la matrice M. De plus, il a été indiqué ci-dessus que la recherche des coefficients liant une composante de P aux différents termes de tMM était facilitée par une homogénéisation du type r4P = Fp(r7 tMM) En pratique, d'une manière analogue à celle décrite pour l'optimisation de I dans la variante décrite, il est possible d'effectuer une optimisation itérative (par la méthode de FIBO NACCI par exemple) sur ltexposant de r. L'expérience montre 7 ex que cette optimisation peut conduire à remplacer r par r ou ex a une valeur non entière comprise par exemple entre 6,5 et 7 et qui peut dépendre de l'environnement. Par ailleurs, en ce qui concerne la modélisation des écarts angulaires et le calcul de la ligne de visée qui résulte de la rotation précise, une simplification permettant de gagner en temps le calcul peut être mise en oeuvre. Pour cela la matrice 3x3, ROT = McM 1MoMco 1 est calculée et les 9 termes de cette matrice ROT sont utilisés comme régresseurs pour la modélisation directe d'écarts angulaires sur G, S ou R. F r est alors telle que G (ou S, ou R) = F r (ROT). Une méthode voisine consiste à modéliser les écarts sur les termes de la matrice ROT elle-même, chaque composante de la matrice ROT = McM 1MoMco I étant exprimée à partir d'un polynôme composé des 9 composantes issues de la première étape. Les angles sont alors extraits de la matrice après régression. Le lecteur comprendra aisément que d'autres modifications peuvent être faites, notamment pour effectuer des calculs plus rapides ou pour améliorer la précision des résultats obtenus sans sortir du cadre de l'invention, c'est-à-dire toujours en effectuant une phase préalable de modélisation des relations complexes reliant les mesures réelles à la position du capteur, ou à l'orientation de la ligne de visée qui lui est liée. Dans un exemple de réalisation de la variante utilisant des orientations prédéterminées de la ligne de visée, la modélisation a été obtenue à partir de 85 orientations et pour chacune de ces orientations 27 points de mesure (i. e. 27 positions de m), soit 9 points équirépartis sur 3 plans équidistants formant un maillage cubique dont le pas est de l'ordre de 3 à 4cm.Les mesures ont été faites pour des valeurs de R=0 et G = 1900, 60 , 30 , 0 , -30 300 > 600 > -900) S = [40 , 200, 00, -200, -4003 pour R = 0 et R = (+100, - 10} et G = [900, 450; 0; -450, -9003 et S = [400, 200, 00 > -200, -400j soit 2295 points de mesure m. Le temps de calcul nécessaire pour la modélisation préalable est fonction de la puissance des moyens de calcul. A titre d'exemple pour l'obtention de Fp à partir de 220 régresseurs, la durée de calcul peut être de l'ordre de 3 heures ; à partir de 84 régresseurs, le calcul d'une fonction de modélisation peut être de l'ordre de 5 minutes. REVENDICATIONS 1. Procédé de détermination d'une orientation liée à un système mobile dans un espace déterminé, utilisant un système magnétique de détection de position et d'orientation comportant un émetteur formé de trois bobines orthogonales ayant une position fixe dans l'espace déterminé et un capteur formé de trois bobines orthogonales fixé au système mobile, r étant la distance entre l'origine du repère émetteur (0) et l'origine du repère capteur (m), un cycle de détection consistant à émettre successivement dans les trois bobines de l'émetteur et à mesurer les composantes de champ générées dans les bobines du récepteur, caractérisé en ce qu'une phase préalable de modélisation de l'espace déterminé consiste à effectuer des mesures des champs magnétiques dans le capteur pour un ensemble de positions connues du système mobile par rapport à un repère dit repère cartographique, des paramètres de modélisation déduits des mesures de référence en chaque point étant calculés pour en déduire la rotation relative du capteur dans le repère cartographique des mesures de référence, et en ce que, en cours d'utilisation, l'orientation courante du système mobile est déterminée directement à partir des paramètres de modélisation constituant une fonction compacte de modélisation position (F ) des effets p perturbateurs, fonction de la position. 2. Procédé selon la revendication 1, caractérisé en ce que la phase préalable de modélisation est complétée par une seconde phase pour améliorer la précision de la rotation relative obtenue à l'étape précédente, par modélisation des écarts angulaires entre des rotations vraies (G2, S2, R2) et la rotation capteur approximée dans la phase précédente (G1 > R1), R1) > cette rotation précise permettant de calculer l'orientation courante (Lv). 3. Procédé selon la revendication 1, caractérisé en ce que ltorientation courante est calculée à partir de la rotation capteur et en ce qu'une seconde phase est mise en oeuvre pour améliorer la précision de l'orientation courante (Lv) par modélisation des écarts entre des orientations vraies et les orientations correspondantes calculées. 4. Procédé selon l'une des revendications précédentes caractérisé en ce que le calcul de l'orientation courante (Lv) à partir de la rotation capteur comporte une étape d'alignement initial à partir d'une position connue de l'orientation (Lvg), pour déterminer une matrice d'initialisation (I) caractérisant la rotation de la position connue de l'orientation par rapport au repère cartographique. 5. Procédé selon la revendication 1, caractérisé en ce que les paramètres de modélisation sont déduits des mesures de référence en chaque point par identification par régression linéaire des composantes d'une matrice de mesure (P) prise comme base à des polynômes obtenus à partir de régresseurs calculés à partir des 6 termes des matrices symétriques tMM, eux-mêmes calculés à partir des différentes matrices de mesure M, les produits tMM étant supposés invariants du fait que les matrices de mesures se déduisent les unes des autres par rotation. 6. Procédé selon la revendication 5, caractérisé en ce que les paramètres de modélisation résultent de l'identification des composantes de r4P à Fp (rex tMM) où F p est une fonction polynomiale des composantes de rex tMM, ex étant un exposant compris entre 6, 5 et 7, l'identification étant effectuée par la méthode des moindres carrés. 7. Procédé selon la revendication 6, caractérisé en ce que l'exposant ex est optimisé pour minimiser les écarts. 8. Procédé selon la revendication 6, caractérisé en ce que les paramètres de modélisation résultent de l'identification par régression linéaire des composantes de r4P à des polynômes obtenus à partir de régresseurs calculés à partir des 6 termes des matrices tMM et des produits scalaires en ligne des matrices M. 9. Procédé selon la revendication 5, caractérisé en ce que les paramètres de modélisation résultent de l'identification, par régression linéaire, des composantes de r4P à des polynômes obtenus à partir de régresseurs calculés à partir des 6 termes des matrices symétriques tMM et de termes ,u. définis par r4(M)ii = r711ittE)ii, r4(M)ii étant les termes diagonaux de la matrice r4M et tMM.. les termes diagonaux de la matrice r tMM. les 10. Procédé selon la revendication 1, caractérisé en ce que la fonction compacte de modélisation (F ) est calculée à p partir de matrices de mesures associées à un ensemble d'orientations définies à partir d'une orientation fixe, des matrices de mesures virtuelles (Mv) étant déduites de ces orientations définies (Lv) et des matrices de mesure correspondantes (M). il. Procédé selon la revendication 10, caractérisé en ce qu'un alignement destiné à évaluer la rotation initiale (I) de l'orientation fixe (Lvo) par rapport au repère cartographique est effectué en même temps que le calcul des matrices de mesures virtuelles (Mv) par une méthode itérative permettant d'optimiser une fonction compacte de modélisation F'p telle que r4Mv = F'p(r7 tMM), et MvM-1ILvo étant le plus proche possible de Lv. 12. Procédé selon l'une des revendications 2 et 3, caractérisé en ce que la modélisation des écarts angulaires est effectuée à partir de rotations vraies connues (G2, S2, R2) et des matrices de mesures relevées associées, par identification des composantes rotations vraies connues (G2, S2, R2) à des polynômes constitués à partir de régresseurs calculés sur les termes des matrices r7 (tMM) et sur les angles des rotations approximées (G1, S1, R1), les coefficients résultant formant la fonction de modélisation (Fr) des effets perturbateurs fonction de la rotation. 13. Procédé selon l'une des revendications 2 et 3, caractérisé en qu'une matrice de rotation ROT = McM 1MoMco -1 est calculée à partir de chaque matrice de mesure M, Mc étant la matrice résultant de M après application de la fonction de modélisation position (Fp) et Mo et Mco étant les matrices M et Mc pour une orientation fixe connue d'initialisation (Lvo), et en ce que la modélisation des écarts angulaires est effectuée par identification des rotations vraies connues ((y2 S2, R2) à des polynômes constitués à partir de régresseurs calculés sur les termes des matrices ROT, les coefficients résultant formant la fonction de modélisation (Fr) des effets perturbateurs fonction de la rotation. 14. Procédé selon l'une quelconque des revendications 2 et 3, caractérisé en ce qu'une matrice de rotation ROT = McM 1MoMcol est calculée à partir de chaque matrice de mesure M, Mc étant la matrice résultant de M après application de la fonction de modélisation position (Fp) et Mo et Mco étant les matrices M et Mc pour une orientation fixe connue d'initialisation (Lvo), et en ce que la modélisation des écarts angulaires est effectuée par identification des termes de la matrice ROT à des polynômes constitués à partir de régresseurs calculés sur les termes de la matrice ROT, les coefficients résultant formant la fonction de modélisation (Fr) des effets perturbateurs fonction de la rotation. 15. Dispositif de détermination d'une orientation liée à un système mobile comportant un émetteur (1) formé de trois bobines magnétiques orthogonales ayant une position fixe dans un espace déterminé et un capteur (2) formé de trois bobines orthogonales fixé au système mobile, destiné à la mise en oeuvre du procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte des moyens de déplacement du capteur et de l'orientation qui lui est liée selon une ensemble de positions connues, des moyens d'émission de champ (20) par émetteur (1) et de mesure de champ par le capteur (2) couplés à des moyens de calcul (10), et en ce qu'un calculateur de modélisation (60) couplé aux moyens de calcul (10) via un circuit d'interface (50) contient sous forme logicielle les algorithmes destinés à la modélisation préalable de l'espace correspondant.