FR 2598808 A1 19871120 FR 8606861 A 19860513 PROCEDE DE REGLAGE DES CORRECTEURS D'HOMOGENEITE DU CHAMP MAGNETIQUE CREE PAR UN AIMANT La présente invention est due à la collaboration du Service National des Champs Intenses (Directeur M. Guy AUBERT), et elle a pour objet un procédé de réglage des correcteurs d'homogénéité du champ magnétique créé par un aimant. Elle trouve son application particulièrement dans le domaine médical o des aimants. sont utilisés pour procéder à des investigations non invasives, en appliquant des méthodes d'imagerie par résonance magnétique nucléaire. Elle peut néanmoins trouver son application dans d'autres domaines, notamment dans celui de la recherche scientifique o des 10 champs intenses produits par des aimants sont utilisés. La résonance magnétique nucléaire est un phénomène d'oscillation du moment magnétique des noyaux des atomes ou molécules cd'un corps à une fréquence qui dépend de l'intensité d'un champ magnétique dans lequel baigne ce corps. Il ressort de cette consta15 tation que si l'intensité du champ magnétique varie, la fréquence du phénomène de résonance varie également. Aussi, pour des raisons technologiques et techniques, dans de telles expérimentations, il est de la plus haute importance que le champ produit par l'aimant soit très homogène dans sa zone d'intérêt. L'homogénéité requise est 20 couramment de l'ordre de quelques parties par million dans le domaine médical, voire de quelques parties par milliard dans le domaine scientifique. Pour arriver à ce but, on essaye de construire des aimants donnant un champ le plus parfaitement homogène. Malheureusement, quel que soit le soin apporté à la constru25 ction des aimants, leur réalisation n'est jamais aussi parfaite que la théorie qui a conduit à leur calcul. Par ailleurs, même si ce défaut peut être éliminé, l'aimant pour être utilisé doit matériellement être placé en un endroit donné. Or aucune des régions de l'espace sur la terre, dans un environnement industriel ou urbain, n'est tota30 J lement exempte cd'éléments magnétiques perturbateurs. Il en résulte qu'une fois installé sur site, le champ produit dans la zone d'intérêt d'un aimant présente des inhomogénéités qu'il faut maintenant tenter de corriger. Le principe de la correction des inhomogénéités de champ est celui de la superposition: on ajoute des bobines, des pièces magnétiques, ou tous autres moyens de produire des champs magnétiques, de manière à corriger les imperfections du champ principal et de manière à produire un champ global homogène. La théorie de 10 la correction a été développée à l'origine par MM. SAUZADE M. et KAN S. K. dans la revue "Advances In Electronics and Electron physics", volume 34, 1973 Academic Press, et une relation de leurs travaux est faite dans la collection des "Techniques de l'ingénieur" chapitre E 4351. Les auteurs sont partis du principe qu'un champ 15 magnétique, en un point quelconque de l'espace, peut être exprimé mathématiquement sous la forme d'une expression dite analytique du type Bz(x,y,z) = f (x,y,z). Dans cette expression Bz représente la composante principale du champ en un point de coordonnée x,y,z, et f est la fonction analytique utilisée. La fonction f que ces auteurs 20 ont explicitée s'écrit sous la forme d'une décomposition polynomiale en x,y,z, cd'ordre, ou de puissance, croissant de zéro à l'infini. Sous certaines restrictions, notamment de proximité des régions intéressées aux moyens de corrections qui en mofidient le champ, les termes d'ordres élevés de cette expression peuvent être négligés. Ces auteurs ont pu tirer de cette décomposition un principe général d'indépendance des corrections. Cette indépendance signifie qu'il est possible de construire des moyens de correction agissant sur des ordres supérieurs, mais pas inférieurs, à un ordre donné. Par exemple pages 4 à 6 de la collection précitée, ces auteurs décrivent 30 une série de douze bobines pour corriger successivement des inhomogénéités de champ jusqu'aux ordres 4 de l'expression analytique en question. Ainsi, par exemple, après avoir corrigé les inhomogénéités de premier ordre du champ principal avec une première bobine conçue dans ce but, on peut corriger des inhomogénéités de - 1 deuxième ordre du champ principal, auquelles se sont ajoutées entre temps des inhomogénéités de deuxième ordre dues à la première bobine, au moyen d'une deuxième bobine de correction. La mission de cette deuxième bobine est de corriger les inhomogénéités de deuxième ordre et de ne pas produire cPd'inhomogénéités de premier ordre. Ainsi de suite, la correction s'effectue successivement pour tous les ordres. Pour le réglage, après chaque correction, avant de passer à la correction suivante, on remesure le champ. Puis on estime la correction à effectuer en comparant les défauts de valeur 10 de champ en des endroits de mesure à des champs qu'on peut obtenir avec la bobine suivante. Cette procédure est laborieuse, elle n'aboutit à de bons résultats qu'après de nombreuses tentatives. Par ailleurs le réglage n'est jamais effectué au mieux. Dans leur complexité, la théorie et le réglage se simplifient par une possibilité de différenciation entre des corrections dites axiales, celles qui correspondent aux variations du champ principal que l'on rencontre quand on se déplace parallèlement à ce champ dans la zone d'intérêt, et des corrections dites radiales qui correspondent aux inhomogénéités que l'on rencontre quand on se déplace 20 perpendiculairement à l'orientation de ce champ principal. Il ressort de cette étude que les corrections axiales et radiales sont indépendantes les unes des autres, au moins si on ne se préoccupe que de la composante du champ principal orientée selon Pl'orientation recherchée de ce champ. Les auteurs précités ont ainsi élaboré, dans un 25 exemple, douze bobines propres à constituer un système de correction complet. Ce système est universellement connu. Il est mis en oeuvre par beaucoup de constructeurs. Plus récemment Mme F. ROMEO et Mr D.I. HOULT ont publié dans la revue MAGNETIC RESONANCE IN MEDICINE 1, 44- 65 de 30 1984 un article intitulé "MAGNET FIELD PROFILING: ANALYSIS AND CORRECTING COIL DESIGN". La démarche de ces auteurs est complémentaire de celle des premiers. Ils montrent en particulier que pour l'expression analytique du champ la décomposition polynomiale à base de coordonnées cartésiennes n'est pas la meilleure. Au contraire, les coordonnées sphériques qui mènent à une expression analytique du champ sous la forme d'une combinaison d'harmoniques sphériques sont plus propices à la compréhension et à la réalisation des structures types de moyens de correction à 5 construire. Les derniers auteurs enseignent en outre une méthodologie de correction. Cette méthodologie comporte les étapes suivantes. Premièrement le champ principal, à corriger, est mesuré. On peut arriver à ce résultat par des moyens classiques. Une fois mesuré on en cherche une expression analytique dans la forme 10 préférée. De cette expression analytique, que l'on compare à une expression analytique idéale (représentant un champ idéal correspondant, par exemple à un Bz constant), on déduit l'expression analytique d'un champ correcteur à produire. Enfin, par des considérations analogues à celles qui ont prévalues pour les précédents 15 auteurs, il est proposé de construire des bobines, ou tout autre moyen de correction, pour arriver à l'homogénéité recherchée du champ principal. Cette méthodologie savante présente cependant un inconvénient majeur: elle n'est pas industrielle. En effet, elle suppose que les systèmes de correction doivent être calculés, 20 fabriqués, et installés spécifiquement pour chaque aimant. Techniquement et financièrement cette méthodologie est inapplicable en dehors du laboratoire. D'une part parce qu'elle correspond à un métier nouveau non encore bien représenté dans l'industrie: celui de calculateur constructeur à la demande de bobines de correction 25 et d'autre part parce que même dans ce cas le coût et le temps passé à de telles réalisations seraient prohibitifs. La présente invention a pour objet de résoudre les inconvénients cités en proposant une méthode qui perfectionne ces deux enseignements pour les rendre plus efficaces et plus industriels. 30 Avec le procédé de réglage de Pinvention, on est capable de régler Phomogénéité dun aimant en quelques heures alors qu'il faut quelques jours avec le premier état de la technique cité. On utilise par ailleurs pour faire le réglage des moyens de correction standard et non pas spécialisés. Ces moyens qui peuvent être des bobines sont seulement indépendants, au sens de l'invention, pour les ordres de correction auxquels ils s'attaquent. Le procédé de l'invention est même tel qu'il est possible de régler le champ d'un aimant alors que cet aimant est muni de moyens de correction dont on ne connait pas la structure et sur lesquels on ne peut agir que par leurs commandes. La présente invention a pour objet un procédé de réglage des correcteurs d'homogénéité du champ magnétique créé par un aimant caractérisé en ce qu'il comporte les phases suivantes: - on mesure le champ réel créé par l'aimant, - on calcule l'expression analytique du champ mesuré, - on calcule les expressions analytiques des champs produits par les moyens de correction en fonction de leur état de réglage, -on calcule l'expression analytique d'un champ idéal à obtenir avec l'aimant, - on procède à une composition algébrique des expressions analytiques des champs de correction pour la rendre égale, et de signe contraire, à la différence des expressions analytiques des champs réel et idéal, - et on règle les moyens de correction en fonction des résultats de 20 cette composition algébrique. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen de la figure qui l'accompagne. Cette description n'est donnée qu'à titre indicatif et nullement limitatif de l'invention. En particulier l'application au domaine médical n'est pas 25 restrictive La figure I unique représente schématiquement en perspective un appareil d'imagerie par résonance magnétique nucléaire du type de ceux utilisés par exemple dans le domaine médical. Un aimant I y produit un champ intense B0 dans un volume d'intérêt 2 ici d'aspect 30 cylindrique circulaire. Des moyens générateurs 3 sous le contrôle d'une unité de commande 4 sont prévus pour alimenter des moyens 5 d'excitation radiofréquence et de différenciation spatiale du moment magnétique de particules placées pour examen dans le volume 2. Les moyens 5, notamment dans leur aspect radiofréquence, sont aussi utilisés pour mesurer une onde de désexcitation représentative du phénomène de résonance magnétique. Ils sont reliés à des moyens 6 de réception et de traitement du signal mesuré. Les moyens 6, également sous le contr8le des moyens de commande 4 peuvent produire des images de coupe du corps examiné. Ces images peuvent être affichées sur des moyens de visualisation 7. Selon l'enseignement préconisé par les deux états de la technique cités, ce système comporte en outre un ensemble 8 de moyens de correction comportant en particulier, dans un exemple, un 10 certain nombre de bobines schématiquement représentées telles que 9,10,11, etc 12. Ces bobines sont alimentées en courant par un générateur 13 par l'intermédiaire de commandes de réglage de l'intensité du courant respectivement 14 à 17. Les commandes 14 à 17 servent à effectuer le réglage des moyens de correction. Comme il a été indiqué ci-dessus les moyens 8, dans l'invention, peuvent être des moyens standards (en général des bobines correctrices au nombre de 12) ou spécialisés. La structure de ces moyens, de ces bobines, peut être connue ou même inconnue de l'opérateur qui va effectuer ler réglage selon l'enseignement de l'invention. 20 Au moyen d'une sonde 18, reliée à un Gaussmrtre 19 et supportée par un dispositif d'exploration 20, on mesure le champ réel créé par l'aimant I dans le volume 2. Le dispositif d'exploration 20 À comporte schématiquement un arbre 21 qui peut coulisser et tourner sur lui-même à l'intérieur d'un palier 22. L'axe de l'arbre 21 est 2. -aligné avec l'axe 23 du volume cylindrique circulaire 2. A un bout de l'arbre 21, est montée une échelle transversale 24 sur laquelle peut être fixée la sonde 18. Suivant la position en rotation de l'arbre 21, réchelle 24 est orientée selon une parmi np directions possibles. La hauteur de la position de la sonde 18 sur l'échelle 24 peut être prise parmi nr positions possibles. Enfin l'arbre 21 peut coulisser dans le 30 support 22 pour que la sonde puisse se trouver successivement à nz abscisses le long de l'axe 23. Un triplet de nombres ir1ip et iz représente ainsi les coordonnées cylindriques d'un point quelconque du volume 2 o on peut placer la sonde 18, et pour lequel on mesure le champ B0 réel. Dans l'invention les mesures effectuées par le Gaussmètre 19 sont introduites dans un ordinateur 27 qui les édite dans une liste 25. En application du procédé de l'invention l'ordinateur 27 peut également éditer une liste 26 d'indications K. des valeurs de réglage des commandes de réglage 14 à 17. Quand les mesures correspondant à la liste 25 sont connues, on peut, au moyen d'un logiciel approprié utilisant en particulier l'enseignement évoqué par le deuxième état de la technique cité, calculer l'expression analytique du champ Bo sous la forme préférée d'une décomposition en harmoniques sphériques. Cette expression analytique peut alors s'écrire de la manière suivante: B Ho + b ([ j - n me H+ k an. m=] + Fr' (nr)npYl( nnl n nn n réduit r/r0 mesure le module de la distance r d'un point quelconque du volume d'intérêt au rayon r0 d'une sphère utile. Les fonctions Pn r)n (cos) et Pm (cos e) sont respectivement les polynomes de Legendre et l es polynomes de Legendre associés d'indices n et _m. Les coefficients Hn Im et sm on n20 représentatifs du champ B: pour un champ B mesuré ils sont rdi /0 mersuenleatodus de cham distancer dun capon qu esrl onqu connus. Dans la pratique, l'ordinateur 27 peut être programmé pour accomplir ce calcul à partir des mesures. On utilise de préférence la correspondance exacte qui existe entre les coordonnées sphériques r, e, et f de. l'expression analytique et le triplet des nombres ir, et iz qui ont servi pour la mesure du champ. Dans lPinvention on a remarqué que s'il était possible de connaître l'expression analytique du champ principal B0, il était aussi possible de connaître l'expression analytique du champ produit par l'un quelconque (9 à 12) des moyens de correction 8. Il est donc 10 possible d'écrire le champ relatif à chacun de ces moyens sous la forme d'une expression analytique similaire à celle vue plus haut et à chaque fois paramétrée par un coefficient Ki: K14 à K17. Ce coefficient, dans le cas de bobines est proportionnel au courant qui les traverse. Par ailleurs, il est aussi possible de calculer 15 l'expression analytique, par exemple aussi selon une décomposition en harmoniques sphériques, d'un champ idéal à obtenir et qui correspondrait par exemple à un Bz constant. La soustraction de l'expression analytique réelle de Pexpression analytique idéale, par exemple correspondant à Bz constant, 20 donne une expression analytique de correction. Cette expression analytique de correction ainsi que les expressions analytiques, paramétrées par Ki, des champs magnétiques de correction peuvent s'écrire de manière concordante. Par exemple tous les coefficients m m Hn, Im, 3 n du champ de correction porteront un indice cor; les 25 mêmes coefficients correspondant aux commandes (14 à 17) comporteront des indices i o i peut valoir de 1 à N, N étant le nombre des commandes. Quand le champ produit par les moyens de correction est égal, et de signe contraire, au champ de correction, on peut écrire que Pl'expression analytique de correction est égale, et 30 de signe contraire, à une composition algébrique des expressions analytiques (à chaque fois multipliées par des coefficients Ki) des champs produits par les bobines de correction. Ceci peut encore s'écrire de la manière suivante: N EAcor. = Ki'.EAi I o EA signifie expression analytique. Dans l'invention on a alors remarqué que cette égalité conduisait à une pluralité d'égalités ou d'équations de la forme: N Hn = Ki À Hn cor. 1 1 N m n cot.r - i ' n i ZIm 10. N im = K. m n cor. i ' n i o n et m peuvent être quelconque (entre 1 et l'infini) puisque les fonctions Pn et P n sont indépendantes. Dans ce système d'équations les inconnues sont les Ki. Pour une valeur de N donnée, (N = 12 dans le premier état de la technique cité), ce système d'équations ne peut admettre de solution si n ou m peuvent être infinis. En effet ceci 20 conduit à une infinité déquations indépendantes avec un nombre limité (N) d'inconnues. Par contre, en se gardant une possibilité de tolérance sur le résultat à atteindre, on peut admettre qu'au delà d'un certain ordre les équations ne seront plus vérifiées. Par exemple on peut admettre qu'au-delà de n = 4 les coefficients de 25 l'expression analytique du champ de correction ainsi que ceux des moyens de correction sont suffisamment faibles pour être négligés. Il en résulte alors que l'on peut se trouver confronté à un système d'équations comportant un nombre donné N d'équations et un même nombre donné N d'inconnues. Un tel système, connu mathéma30 tiquement sous le nom de système de CRAMER, est soluble. La solution qui comporte notamment une opération d'inversion de matrice est typiquement de celles qui peuvent être exécutées par l'ordinateur 27. Ainsi il résulte que, connaissant tous les coefficients de décomposition en harmoniques sphériques liés au champ de correction et liés à chacun (i) des moyens de correction, il est possible de trouver simplement la liste des valeurs Ki de réglage de ces 5 moyens. Dans le cas de bobines les valeurs Ki sont proportionnelles aux courants qui traversent chacune des bobines. On conçoit cependant que cette affirmation ne serait pas vraie si au lieu de disposer d'un nombre N de bobines de correction indépendantes on disposait N fois de la même bobine. Cette impossibilité s'exprime autrement en 10 disant qu'il est nécessaire pour qu'une solution puisse être trouvée selon le procédé de calcul de l'invention que les bobines de correction produisent des champs qui, au regard de la décomposition en harmoniques sphériques (du moins pour les premiers ordres intéressés), puissent être considérés comme algébriquement indé15 - pendants. Cette indépendance a par exemple pour effet qu'il doit être impossible de reproduire le champ d'une quelconque bobine avec une combinaison algébrique des champs produits par les autres bobines de correction, et ceci au moins pour les premiers ordres considérés. Cependant, on peut se rendre compte que cette con20 trainte supplémentaire d'indépendance dans l'invention est très simple à réaliser et est beaucoup moins contraignante que celle correspondant à l'état de la technique cité. Dans celui-ci en effet elle s'exprime différemment, mais on y attend quand même d'une bobine corrigeant des inhomogénéités à un ordre donné qu'elle 25 n'induise pas d'inhomogénéités à un ordre inférieur (ou encore qu'une bobine de correction des inhomogénéités radiales n'induise pas d'ilnhomogénéités axiales et vice versa). La connaissance de la décomposition en harmoniques sphériques du champ de chacune des bobines de correction peut être 30 connue de diverses manières. En connaissant la structure des bobines de correction et en suivant l'enseignement du deuxième état de la technique citée, il est possible de calculer l'expression analytique du champ théorique qu'elles produisent. En fait ces bobines, comme Paimant à corriger, ne produisent pas un champ théorique mais un 1 I., 1 1 -. champ réel légèrement différent du champ théorique. Il en résulte une différence de l'expression analytique qui correspond à ce champ réel. Si on admet que le champ principal connaît des inhomogénéîtés de l'ordre du millième, le champ nominal produit par les moyens de 5 correction doit être de cet ordre de grandeur. Construit avec le même soin que l'aimant, les moyens de correction connaissent alors des défauts du même ordre. En conséquence, les inhomogénéités des champs de correction sont elles-mêmes du millième du millième, c'est-à-dire de l'ordre du millionnième du champ principal. Dans ces 10 conditions, dans une première approche, elles peuvent être négligées et l'expression analytique du champ théorique peut être utilisée. Cependant si tel n'est pas le cas, ou, dans un cas particulièrement intéressant si la structure des moyens de correction n'est pas connue, il est possible de calculer l'expression analytique du 15 champ que chacun d'eux produit. Par exemple pour des bobines on les fait traverser tour à tour par un courant nominal, proportionnel à un coefficient K0, et on mesure le champ réel qu'elles produisent. Puis on calcule de la même manière que pour le champ principal la décomposition analytique à laquelle ce champ correspond. La pré20 sente invention enseigne ainsi que, quelque soient les moyens de corrections existants et pourvu qu'ils soient indépendants au sens de Pinvention, il est possible de remédier aux inhomogénéités dès qu'on en a calculé l'expression analytique. Un autre avantage de la méthode consiste à indiquer avec 25 quelle précision géométrique doit être menée l'exploration du volume 2, ainsi que l'ordre des expressions analytiques de champ dont il faut se servir pour réduire des inhomogénérté jusqu'à une valeur donnée. Le nombre de points de mesure est égal à P o: P = nz(l + (nr- l)np) Dans la pratique, on s'arrange pour que nz soit un nombre impair de manière à ce qu'une des abscisses corresponde au centre du volume d'intérêt 2. Pour nz, nr, et np connus, on peut donc relever le champ en P points et calculer son expression analytique. Cette expression est forcement limitée à un ordre donné. Connaissant cette expression analytique, on calcule le champ théorique auquelle elle correspond. Ce champ théorique est différent du champ réel mesuré puisqu'on a négligé les ordres supérieurs dans l'expression analytique. *' On calcule alors l'écart quadratique moyen entre le champ réel mesuré et le champ théorique calculé. Si cet écart quadratique moyen est supérieur à la précision de mesure du Gaussmètre 19 (ou aussi à la tolérance recherchée sur l'homogénéité), on en déduit que l'expression analytique est incorrecte. Elle peut être incorrecte pour deux raisons. Premièrement parce que le nombre de points de mesure est insuffisant. Deuxièmement, corrélativement, parce que l'ordre de décomposition en harmoniques sphériques est insuffisant également. Il convient dans ce cas d'augmenter nz, nr, np et l'ordre de décomposition jusqu'à ce que lécart quadratique moyen soit de 'Pordre ou inférieur à la précision voulue. La décomposition en harmoniques sphériques nous apprend que m doit être inférieur ou égal à n. Sachant que m y apparaît dans les fonctions propres sous la forme de cos m et sin m, il en résulte que la limitation à une valeur donnée du nombre np limite le nombre m. En effet, l'exploration radiale selon un nombre donné np de directions privilégiées est identique à un échantillonnage, et en appliquant la théorie de Shannon, on sait qu'on ne peut pas valablement définir de fonctions propres sin m et cos m ou m est supérieur à la partie entière de np/2. Des considérations mathé25 matiques et physiques du même ordre conduisent à un choix de n maximum égal à (nr + nz - 2). Ainsi, connaissant nr, np, nz donc P, on en déduit l'ordre supérieur de l'expression analytique. L'augmentation de l'un entraîne celle de l'autre. Dans un exemple on a ainsi découvert pour un volume cylindrique dintérêt, dont le -* 30 diamètre, égal à la longueur, vaut de l'ordre de 50cm, que la connaissance du champ en 451 points (nr = 6, np = 8, nz = 11) permet cfexprimer analytiquement le champ avec des harmoniques sphériques jusqu'à l'ordre n = 15 et m = 4 avec un écart quadratique moyen d'une partie par million. La présente invention quoique décrite dans une solution o le réglage est simple (par réglage de courant dans des bobines de correction) n'est cependant pas limitée à cette application. En particulier il est possible d'utiliser des morceaux d'aimants permanents de formes, de dimensions, et d'aimantation appropriées ainsi que des pièces aimantables (en fer doux) de formes, de dimensions et de susceptibilité magnétique choisies. Ces "cales" magnétiques ou magnétisables sont placées à proximité de l'aimant. Elles exercent leur influence principalement sur le volume d'intérêt 2. Le réglage 1 peut être exercé sous une forme quelconque mais d'une manière 10 préférée il est obtenu en modifiant l'orientation de ces cales. De telles réalisations sont décrites dans une demande de brevet conjointe déposée le même jour par la demanderesse. Dans ce cas, les paramètres Ki doivent être remplacés par des fonctions de para15 métrisation. Ces fonctions de paramétrisation correspondent à des changements d'orientation des cales magnétiques. La notion de composition algébrique des expressions analytiques, revendiquée cidessus, s'étend aussi à ce cas. Elle y signifie en fait composition algébrique des termes des expressions analytiques. En effet ces 20 termes y sont paramétrés en fonction des angles de pointage des cales. REVENDICATIONS 1. Procédé de réglage des correcteurs d'homogénéité du champ magnétique créé par un aimant, caractérisé en ce qu'il comporte les phases suivantes: -on mesure (18-25) le champ réel créé par l'aimant, -on calcule (24) l'expression analytique de ce champ mesuré, - on calcule les expressions analytiques des champs produits par les moyens de correction en fonction de leur état de réglage, - on calcule l'expression analytique d'un champ idéal à obtenir avec raimant, - on procède (24) à une composition algébrique des expressions 10 analytiques des champs de correction pour la rendre égale, et de signe contraire, à la différence des expressions analytiques du champ réel et idéal, - et on règle (19-17) les moyens de correction en fonction des résultats de cette composition algébrique. 15 2. Procédé selon la revendication 1, caractérisé en ce que pour calculer Pexpression analytique des champs produits par les moyens de correction, on mesure le champ que produit chacun de ces moyens dans un état de réglage donné. 3. Procédé selon l'une quelconque des revendications I ou 2, caractérisé en ce que pour calculer les expressions analytiques, on calcule les décompositions en harmoniques sphériques des champs produits. 4. Procédé selon l'une quelconque des revendications 1 ou 2, 25 caractérisé en ce que pour le calcul des expressions analytiques des champs, on calcule les décompositions en harmoniques cylindriques des champs produits. 5. Procédé selon l'un quelconques des revendications I à 4, caractérisé en ce que pour mesurer le champ on déplace une sonde 30 de mesure en un certain nombre de points répartis cylindriquement à rintérieur d'un volume d'intérêt. w t Z, 6. Procédé selon la revendication 5, caractérisé en ce qu'après avoir calculé la décomposition en harmoniques sphériques du champ mesuré, - on calcule le champ magnétique correspondant à cette décom: position, on compare le champ magnétique calculé au champ magnétique mesuré, - on déduit de cette comparaison le nombre minimal de points de mesure et l'ordre minimal de la décomposition compatibles avec la 10 précision de l'appareil de mesure utilisé pour faire la mesure et/ou avec la tolérance à respecter par le champ réel.