FR 3002328 A1 20140822 FR 1351396 A 20130219 Domaine de l'invention La présente invention a pour objet un procédé et un dispositif de correction de champ magnétique, pouvant être utilisés notamment dans la technique permettant de créer des spectres et des images par résonance magnétique nucléaire (RMN), cette technique étant également connue sous le nom d'imagerie par résonance magnétique (IRM). L'invention concerne également un système d'imagerie par résonance magnétique mettant en oeuvre un tel procédé de correction de 10 champ magnétique. Art antérieur L'IRM et la RMN reposent sur l'utilisation de champs magnétiques, 15 dont un champ magnétique dit « pnncipal» qui doit être aussi uniforme que possible dans la région examinée ou zone d'intérêt ZI. L'usage a consacré le qualificatif « homogène » pour désigner ce caractère uniforme. Ce champ magnétique de grande homogénéité est généré par des aimants dont les plus répandus aujourd'hui sont constitués de 20 bobinages supraconducteurs qui transportent les courants électriques générant le champ sans aucune dissipation d'énergie, à condition qu'ils soient maintenus à très basse température. Un tel dispositif d'aimant a généralement l'apparence extérieure d'un tunnel cylindrique dans lequel on introduit l'objet ou le patient à imager. 25 L'analyse d'échantillons anisotropes, par exemple des solides, par RMN nécessite de faire tourner l'échantillon autour d'un axe orienté à l'angle dit "magique" (arctan(V2) - 54,7°). L'échantillon est généralement de forme cylindrique, et sa longueur est souvent largement plus grande que son diamètre (entre un facteur 2 et un facteur 10). Ce rapport de 30 forme est aussi retrouvé dans la RMN d'échantillons isotropes, par exemple des liquides, où, typiquement, l'échantillon est contenu dans un tube de 5 mm de diamètre et la hauteur d'échantillon dans le tube se situe plutôt de l'ordre du centimètre ou plus. La spectroscopie RMN, qu'elle soit de l'état anisotrope (solide) ou de l'état isotrope (liquide), nécessite un champ magnétique ambiant extrêmement uniforme dans l'espace. L'analyse des échantillons se fait en effet par l'analyse du spectre RMN, qui est constitué par la réponse en fréquence de l'échantillon lorsqu'il est excité par des impulsions radiofréquence (RF). Cette réponse dépend directement de la valeur locale du champ magnétique. La fréquence de réponse locale (c'est-à-dire celle d'un noyau atomique donné) est appelée fréquence de Larmor (f est donnée par 'Larmor, e, fLarmor = \IBO, (1) où y est le rapport gyromagnétique fréquentiel du noyau et Bo est le module du champ magnétique statique local. Cette valeur locale peut être affectée par la composition chimique de l'échantillon, ce qui permet de révéler des informations cruciales sur la nature, la composition et les propriétés de l'échantillon. L'ordre de grandeur des interactions qui affectent le spectre RMN est la partie par million (ppm). Ceci implique que le champ magnétique ambiant doit lui-même être bien plus homogène que 1 ppm sur toute l'étendue de l'échantillon. On conçoit pour cela des aimants spécifiquement configurés qui produisent un champ extrêmement homogène, mais rarement suffisamment pour permettre la spectroscopie sans ajustements supplémentaires. Ceci est d'autant plus vrai que l'échantillon lui-même est souvent à l'origine de distorsions du champ magnétique à cause de sa susceptibilité magnétique intrinsèque. On est donc obligé d'avoir recours à des bobines spécifiques de correction de champ magnétique, dites "bobines de shirn", qui permettent 3002 32 8 3 de compenser les dernières imperfections de l'aimant afin d'obtenir l'uniformité nécessaire. Chaque aimant de RMN "standard", que ce soit pour le liquide ou le solide, est donc équipé de ce que l'on appelle un "fourreau de shims" qui n'est autre qu'un ensemble de bobines de 5 correction de champ magnétique (bobines de shim) dont les courants peuvent être pilotés indépendamment. L'ajustement du courant permet de contrôler l'effet de chaque bobine sur la distribution spatiale du champ. Le fourreau de shims est généralement un objet cylindrique d'épaisseur modérée afin d'occuper le minimum de place dans le trou de l'aimant. Il 10 est bien sûr coaxial avec le trou de l'aimant. On a représenté de façon schématique sur la figure 7 un exemple de spectromètre de RMN mettant en oeuvre des bobines de correction de champ magnétique. Un tel spectromètre comprend une unité d'expérimentation 1, une 15 unité d'activation 2 comprenant un ensemble de composants électroniques et une unité de commande 3 comprenant un ordinateur ou processeur de traitement. L'unité d'expérimentation 1 comprend, à l'intérieur du trou cylindrique d'un aimant (non représenté) de création du champ 20 magnétique principal, d'axe z, un échantillon 17, autour duquel sont disposées des bobines radiofréquence 15, elles-mêmes entourées de bobines de gradient 14 et de bobines 16 de correction de champ magnétique (bobines de shim). L'unité d'activation 2 comprend une unité 21 d'alimentation des 25 bobines de shim 16, une unité 22 d'alimentation des bobines de gradient 14, une unité 23 d'émission de signaux RF vers les bobines RF 15 et de réception des signaux RF émis par les bobines RF 15. L'unité de commande 3 comprend un module 31 de détermination des valeurs des signaux à fournir par l'unité 21 d'alimentation des bobines de shim 16, un module 32 de détermination des valeurs des signaux à fournir par l'unité 22 d'alimentation des bobines de gradient 14, un module 33 d'envoi d'impulsions RF vers l'unité 23 reliée aux bobines RF 15 et un module 34 de réception de signaux de RMN radiofréquence fournis par l'unité 23 reliée aux bobines RF 15. L'obtention d'un champ magnétique hautement uniforme est une tâche complexe. Afin de simplifier le problème, on a recours à une approximation majeure. Celle-ci consiste à considérer que, dans un champ homogène, les variations du module du champ magnétique sont dominées par celle de la composante principale du champ. Dans le cas des aimants supraconducteurs, le champ est généralement vertical dans le référentiel du laboratoire. On a donc pris l'habitude d'assimiler la quantité Bo à Bz, la composante verticale du champ magnétique. On peut alors montrer que les variations de Bo sont dominées par celles de Bz, tant que ces variations sont faibles comparées à la valeur du champ à l'origine. Cette approximation permet donc de ne considérer qu'une composante cartésienne du champ magnétique. Il est possible de montrer que, dans une région vide de sources de champ magnétique, cette composante peut être développée en série d'harmoniques sphériques. On peut ainsi écrire de manière générale, dans le système de coordonnées sphériques('"'), et dans la sphère intérieure aux sources : , B, = Z0 +Ir"[Z'Pn(cosa) + 1(X cosmço +17: sin inço)Pnm (cos,/ m=i (2) où ("'Ç°) sont les coordonnées sphériques du point considéré x = r sin Ocos 9 dans un repère d'axe Oz tel que z = r cost9 et avec y = r sin a sin ço . Les Pn sont les polynômes de Legendre de degré n et ni les Pn sont les polynômes de Legendre associés de degré n et d'ordre in. Ce développement est unique et valable à l'intérieur de la plus grande boule de centre 0 magnétiquement vide. Si la configuration des sources de champ est axisymétrique d'axe Oz, on obtient la forme simple suivante : B z = Z0 + Z 'Pn (cos ,9) n =1 (2.1) C'est la forme appropriée pour un ensemble de spires coaxiales d'axe . Elle se réduit à une forme encore plus simple pour B(z) sur l'axe, soit : B z = Z0 + Ez' Zn n =1 (2.2) Les termes Zn, Xnm, et Y0n1 sont des termes définis par la géométrie des sources de champ. L'équation fondamentale (2) susmentionnée fournit l'outil nécessaire à la résolution du problème de l'homogénéité. En effet, à l'intérieur de la sphère intérieure (i.e. la plus grande sphère ne contenant aucune source de champ), les variations du champ dues à un terme Cn (c'est-à-dire Zn, Xnm, ou Ynni ) de degré n sont en Cnrn. Or on montre (par analyse dimensionnelle ou calcul direct) que Co varie comme 1/r00. Si rs est la plus grande distance d'un point de l'échantillon à l'origine, la contribution du terme Cn au champ variera comme Cn(rs/ro). Pour un échantillon de taille donnée (rs Définition et objet de l'invention La présente invention vise à remédier aux inconvénients précités et à permettre la réalisation de façon simplifiée d'un dispositif de correction d'homogénéité d'un champ magnétique pour un système de spectroscopie ou d'imagerie par résonance magnétique. L'invention vise encore à définir un procédé de réalisation d'un tel dispositif qui soit simplifié et permette néanmoins une optimisation de 10 l'homogénéisation du champ magnétique créé dans le volume d'intérêt. Ces buts sont atteints, conformément à l'invention, grâce à un procédé de correction de champ magnétique d'un système d'imagerie par résonance magnétique, ledit système comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone 15 d'intérêt ZI, un dispositif de support d'un échantillon, avec une orientation d'une dimension principale de l'échantillon selon un angle 00 différent de zéro par rapport à ladite direction Oz, des bobines de gradient, et des bobines radiofréquence et la correction étant effectuée à l'aide de bobines de correction disposées autour dudit dispositif de support de l'échantillon, 20 caractérisé en ce qu'il comprend les étapes suivantes : - définir un système de coordonnées incliné Ox'y'z' attaché à l'échantillon avec un axe principal Oz' correspondant à ladite dimension principale de l'échantillon orientée selon un angle 00 différent de zéro par rapport à ladite direction Oz du champ magnétique homogène, 25 - déterminer par un développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon la forme des bobines de correction à partir d'une fonction flux F, chaque bobine de correction correspondant à un terme du développement en harmoniques sphériques, - tracer pour chaque bobine de correction des iso-contours de la fonction flux F espacés de manière régulière entre des bornes de la fonction flux F sur un cylindre, et - former à partir des iso-contours du cylindre de chaque bobine de correction des éléments de bobinage pour définir une bobine cylindrique de correction positionnée autour du dispositif de support de l'échantillon, l'axe de ladite bobine cylindrique étant confondu avec ladite direction Oz. Selon un mode de réalisation avantageux, ledit angle 00 différent de zéro correspond à un angle dit magique égal à arctan (V2) - 54,7°. Selon un mode de réalisation particulier qui est à la fois simple et performant, on définit des première, deuxième et troisième bobines de correction correspondant respectivement à des termes axiaux Z'1, r2, r3 du développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon. Avantageusement, on alimente chacune des bobines de correction à l'aide d'un courant dont la valeur est ajustable. L'invention concerne également un dispositif de correction de champ magnétique d'un système d'imagerie par résonance magnétique, ledit système comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone d'intérêt ZI, un dispositif de support d'un échantillon, avec une orientation d'une dimension principale de l'échantillon selon un angle 00 différent de zéro par rapport à ladite direction Oz, des bobines de gradient et des bobines radiofréquence, ledit dispositif de correction de champ magnétique comprenant un ensemble de bobines de correction positionnées autour du dispositif de support de l'échantillon, caractérisé en ce que chaque bobine de correction présente un axe confondu avec la direction Oz et comprend des éléments de bobinage réalisés à partir des iso-contours d'une fonction flux F espacés de manière régulière entre des bornes de la fonction flux F sur un cylindre, la forme des iso-contours étant déterminée à partir d'un développement en harmoniques sphériques dans un système de coordonnées incliné Ox'y'z' attaché à l'échantillon avec un axe principal Oz' correspondant à ladite dimension principale de l'échantillon orientée selon un angle 00 différent de zéro par rapport à ladite direction principale Oz du champ magnétique homogène, chaque bobine de correction correspondant à un terme du développement en harmoniques sphériques. Avantageusement, ledit angle eo différent de zéro correspond à un angle dit magique égal à arctan (V2) - 54,70. Selon un mode de réalisation préférentiel, l'ensemble de bobines de correction comprend des première, deuxième et troisième bobines de correction correspondant respectivement à des termes axiaux Z'1, r2, r3 du développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon. L'invention concerne encore un système d'imagerie par résonance magnétique, comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone d'intérêt ZI, un dispositif de support d'un échantillon, avec une orientation d'une dimension principale de l'échantillon selon un angle 00 différent de zéro par rapport à ladite direction Oz, des bobines de gradient et des bobines radiofréquence, caractérisé en ce qu'il comprend un dispositif de correction de champ magnétique tel que défini ci-dessus. Ladite direction Oz peut notamment être verticale ou horizontale, selon les applications envisagées. On notera que favoriser les termes le long de l'axe à l'angle magique, comme selon des modes préférentiels de réalisation de l'invention, est d'autant plus pertinent que la rotation de l'échantillon dans une expérience de rotation à l'angle magique (en anglais « Magic Angle Spinning » ou « MAS ») atteint des vitesses de l'ordre de 10 kHz ou plus, ce qui a pour effet d'annuler l'influence, même importante, des termes non axiaux. Par conséquent, se placer dans le référentiel incliné pour traiter les termes du DHS de la composante de champ magnétique Bz est extrêmement important et efficace. Selon des modes particuliers de réalisation, les éléments de bobinage des bobines de correction peuvent comprendre par exemple des pistes ou fils conducteurs sur des supports isolants. Brève description des dessins D'autres caractéristiques et avantages de l'invention ressortiront de la description suivante de modes particuliers de réalisation de l'invention donnés à titre d'exemples en référence aux dessins annexés, sur lesquels : La Figure 1 représente une vue schématique de la position d'un échantillon pour un système de spectroscopie RMN, auquel est applicable un dispositif de correction de champ magnétique selon l'invention, La Figure 2 montre une vue schématique en coupe axiale, d'un exemple de dispositif de support d'échantillon permettant une inclinaison de cet échantillon et le support d'un dispositif de correction de champ magnétique selon l'invention, La Figure 3 montre, sous une forme développée dans un plan, un exemple des iso-contours de la fonction flux pour une bobine cylindrique correctrice générant un profil de champ dominé par le terme axial Z'i. d'un DHS et dont la puissance dissipée pour une valeur donnée de Z'1 est minimum, La Figure 4 montre, sous une forme développée dans un plan, un exemple des iso-contours de la fonction flux pour une bobine cylindrique correctrice générant un profil de champ dominé par le terme axial Z'2 d'un DHS et dont la puissance dissipée pour une valeur donnée de r2 est minimum, - La Figure 5 montre, sous une forme développée dans un plan, un exemple des iso-contours de la fonction flux pour une bobine cylindrique correctrice générant un profil de champ dominé par le terme axial r3 d'un DHS et dont la puissance dissipée pour une valeur donnée de r3 est minimum, - La Figure 6 est un schéma-bloc simplifié d'un exemple de système de spectroscopie ou d'imagerie par résonance magnétique mettant en oeuvre un dispositif de correction de champ magnétique selon l'invention, La Figure 7 représente une vue schématique d'un exemple de système de spectroscopie et d'imagerie par résonance magnétique connu, - La figure 8 montre un exemple de masque pour la fabrication d'un exemple de bobine de correction de gradient Gx, - La figure 9 montre un exemple de masque pour la fabrication d'un exemple de bobine de correction de gradient Gy, - La figure 10 montre un exemple de masque pour la fabrication d'un exemple de bobine de correction de gradient Gz, et La Figure 11 représente une vue schématique de la position d'un échantillon pour un système connu de spectroscopie RMN. Description détaillée de modes de réalisation préférentiels Selon la présente invention, on propose des bobines de correction de champ magnétique ou bobines de « shim », destinées à prendre en compte de façon préférentielle dans un développement en harmoniques sphériques (DHS) respectivement les termes axiaux r1, r2 et Z'3, le 3002 32 8 12 dispositif de correction comprenant lesdites bobines étant conçu spécifiquement pour l'application à un système d'imagerie par résonance magnétique comprenant un échantillon incliné par rapport à la direction du champ magnétique principal en étant orienté selon un angle 00 qui est 5 de préférence égal à l'angle magique (54,7°). Ces bobines s'inscrivent toujours sur un cylindre coaxial avec le trou de l'aimant créant le champ principal, et peuvent avoir les mêmes dimensions que les fourreaux contenant des bobines de correction (dits « fourreaux de shirn ») déjà en service, permettant leur utilisation directe 10 dans des installations existantes. Elles corrigent directement les termes du DHS attachés au référentiel incliné. On réduit ainsi de manière radicale le travail de l'opérateur qui doit procéder à des corrections des irrégularités du champ magnétique principal créé par son aimant. Il est fait référence ici à des termes axiaux ri., Z'2 et r3 jusqu'au 15 degré 3, car ils sont considérés comme les plus importants, mais l'invention n'est bien sûr aucunement limitée à ce nombre et le procédé selon l'invention permet de calculer les caractéristiques de bobines compensant d'autres termes, de degré plus élevé ou même des termes non axiaux. 20 Un schéma simplifié du contexte d'application de l'invention est présenté sur la Figure 1. On observe l'échantillon 117 placé le long de l'axe z' qui se trouve à l'angle magique par rapport à l'axe z du champ magnétique Bo. L'échantillon 117 peut tourner autour de cet axe z' avec une fréquence cor afin de moyenner les interactions anisotropes et les 25 inhomogénéités résiduelles du champ. Les bobines de correction 116 (Figure 2) sont dessinées pour effectuer des corrections sur les termes du DHS rattaché au système de coordonnées incliné Ox'y'zi. On décrira ci-dessous un exemple de calcul des caractéristiques de bobines de correction 116 selon l'invention. On commence par restreindre la zone d'existence des courants à une surface cylindrique de rayon a et de longueur 2b. Il s'agit ici d'une situation statique et l'on a donc ^-^^---^,, où j est la densité de courant. On a de plus, de nouveau du fait du régime statique, aux limites =0. où n est le vecteur normal à la surface. La densité de courant est donc surfacique et sera baptisée k dans ce qui suit. De par la géométrie imposée ici, on adopte le système de coordonnées cylindriques (p, p, z), avec pour origine le centre du cylindre. Le cylindre est ainsi délimité par deux plans de cotes z = b et z = -b et l'axe du cylindre passe par l'origine en constituant son axe de symétrie. L'axe du cylindre est celui de l'aimant de RMN, c'est à dire Oz dans le référentiel du laboratoire. La distribution de courant prend donc la forme générale On peut ensuite appliquer la loi de Biot et Savart pour trouver la forme générale de la distribution de champ magnétique induite par les courants. On s'intéresse ici à la composante 8, du champ, qui est la composante dominante dans un aimant de RMN standard. On peut écrire, en utilisant les coordonnées cartésiennes (x0,y0,z0), où Oz est le long de l'axe de l'aimant, pour repérer le point de calcul du champ, et les coordonnées cylindriques pour l'intégration sur la surface portant les courants : P f O J°0 ),s0 - 'ro, Yo, = - 47T I ark)(1:: Ia- (6) Il importe ensuite de définir de manière appropriée le repère (Ox'yZ) dans lequel le DHS est exprimé. Dans le cas d'un échantillon cylindrique 117 tournant, on a tout intérêt à fixer Oz' comme l'axe de révolution de l'échantillon. Il est bien connu que la rotation rapide autour de l'axe Oz' permet d'annuler l'effet des termes non axiaux dans le signal RMN. On se simplifie ainsi la tâche en ne devant ajuster que des termes axiaux. On pourrait parler de repère "propre" de l'échantillon. Les termes axiaux ne sont autres que les n-dérivées le long de l'axe Oicalculées à l'origine. On a ainsi , 1 Z" n! dz'n O (7) L'indice alpha de B indique la composante de B le long d'un axe d'orientation arbitraire alpha. Alpha peut être par exemple x, y ou z. Dans cette application spécifique où le champ principal pointe le long de l'axe z et tronque les composantes transverses, alpha peut donc être z. Il est aisé d'exprimer uz, dans le repère (x,y,z) en fonction de son angle d'inclinaison 8 et en supposant cette inclinaison dans le plan x0z. On a ainsi = cos 0 uzo + sin 8 uxo. (8) On peut donc écrire la dérivée le long de Oi comme df + sin() = cos() a Zo (9) Dans le cas de l'angle magique, on trouve ainsi : b ad0c17 47r fbrn.- -- Zin (10) = a k 0(z,0) (a2 + z2)2 (3az + -\12(2a 2- Z2)COS0) 5 (a2 + z2)2 ko(z,0) (12) 3 = 1 a(2z2 3a2) 2-\,(4za2 z3) cos 0 2a(3a2 2z2) cos2 2 (a2 +z2) (13) / 5az(9a 2 r- - 22'2) (63a 2Z2 -36a- 9 9 COS (a2 + z2)2 (a2 + z2)2 30az(5a 2 - 2z2)COS2 q5 + 10a 2 (4Vie 2 - 3-\/2z2) 9 9 COS, k0(Z,0) (14) (a2 z2)2 (a2 z2)2 Comme cela a été vu plus haut, on a : 3 Ainsi que (16) 16 Ces équations impliquent que k est le rotationnel d'un vecteur F tel que F est appelée fonction flux. On peut alors montrer que l'on obtient la distribution de courant minimisant la puissance dissipée pour une valeur d'un terme axial Z'n donné lorsque la fonction F est solution d'une équation de Poisson. Le second membre de cette équation est donné par les contraintes fixées sur le profil de champ, c'est à dire les valeurs relatives désirées des termes axiaux Z'n (y compris l'annulation de certains termes si nécessaire). La fonction flux F minimise ainsi une fonction objectif tout en respectant des contraintes. La fonction objectif P' peut être considérée 20 comme proportionnelle à la puissance P dissipée par effet Joule. On a ainsi les relations suivantes : a P, cre = ff (k2 + k2) clçodz Z (19) (20) 3002 3 2 8 17 où e est l'épaisseur de la couche mince conductrice de conductivité électrique a . Il faudra transformer cette expression pour qu'elle ne fasse intervenir que la fonction flux et il faudra faire de même pour les 5 expressions des contraintes. Ces dernières sont les expressions de coefficients du DHS ou de coefficients Z", qui doivent prendre des valeurs fixées ou ne pas dépasser en valeur absolue une limite donnée. Comme ils sont relatifs à la composante B, du champ produit, ils ne dépendent pas de k, et dépendent linéairement de k, avec la forme générique 10 suivante : = ffs k,,,fi(ço,z)dçodz (21) Une fois la fonction flux F trouvée, il suffit de tracer des isocontours de cette fonction Fespacés de manière régulière entre les bornes 15 de F (en incluant zéro si il y a changement de signe) sur le cylindre pour obtenir, soit la position de boucles décrites par un fil conducteur ou une piste conductrice, soit la position de découpes dans une plaque de conducteur (par exemple du type tôle de cuivre). Les figures 3 à 5 montrent des exemples de bobines cylindriques 20 116A, 116B, 116C générant chacune un profil de champ magnétique dominé par un terme axial rn. Pour chaque figure, Les abscisses correspondent à la direction parallèle à l'axe z du cylindre et les ordonnées à la position angulaire sur le cylindre. Il suffit donc d'enrouler ce dessin autour d'un cylindre de rayon 25 approprié pour obtenir la bobine. Le schéma est à l'échelle et seule la proportion entre le rayon a du cylindre et sa longueur 2b doit être conservée constante pour conserver les propriétés calculées (à part l'intensité du terme généré par unité de puissance qui décroît lorsque le rayon a augmente). La figure 3 montre un exemple de bobine cylindrique 116A générant un profil de champ dominé par Z'1 et dont la puissance dissipée 5 pour une valeur donnée de Z'l est minimum. La figure 4 montre un exemple de bobine cylindrique 116B générant un profil de champ dominé par r2 et dont la puissance dissipée pour une valeur donnée de r2 est minimum. La figure 5 montre un exemple de bobine cylindrique 116C 10 générant un profil de champ dominé par ZI3 et dont la puissance dissipée pour une valeur donnée de r3 est minimum. Pour la réalisation des bobines cylindriques de correction, on peut utiliser par exemple du fil de cuivre isolé de section constante, circulaire ou rectangulaire, collé sur les contours de la fonction de flux. Comme on 15 le voit sur les exemples des figures 3 à 5, il existe un ensemble de différents iso-contours imbriqués les uns dans les autres, à la manière de courbes de niveau. Il convient de passer d'un iso-contour à un iso-contour voisin en ouvrant la boucle et en utilisant un segment de fil droit dont l'emplacement est choisi pour ne pas contribuer au champ principal du 20 dispositif de correction. Ainsi, les boucles fermées des fils conducteurs superposés aux iso-contours sont connectées en série et de préférence les segments de liaison sont disposés parallèlement à l'axe Oz pour qu'ils ne créent pas de champ additionnel dans cette direction. L'invention permet également de réaliser des bobines correctrices 25 de gradient Gx, Gy, et G, à partir des iso-contours d'une fonction flux, comme illustré par exemple sur les figures 8 à 10, relatives respectivement à une bobine correctrice de gradient Gx, à une bobine correctrice de gradient Gy et à une bobine correctrice de gradient G. Les figures 8 à 10 représentent respectivement des masques 114A, 114B et 114C pour la réalisation de pistes conductrices sur une face d'un circuit imprimé de manière à constituer des éléments de bobinage pour les gradients Gx, Gy, et G. On voit en particulier sur ces masques les passages pour le courant entre les pistes correspondant à des iso-contours voisins, afin de définir des connexions en série. Il est possible d'utiliser un masque sur chacune des faces d'un circuit imprimé afin de doubler l'efficacité de la bobine. Le courant passe d'une face à l'autre par des vias placés au centre des contours du centre. Les bobines correctrices de gradient permettent d'augmenter la linéarité de la variation du champ dans une direction fixée Oz, dans une région d'intérêt, de même que des bobines de correction telles que les bobines de correction 116A, 116B et 116C visent à rendre aussi invariable que possible dans la région d'intérêt le champ magnétique dans la direction Oz. On a représenté sur la figure 2 un exemple de dispositif de spectroscopie RMN comprenant un boîtier 180 inséré dans le tunnel d'un aimant (non représenté) qui crée dans une zone d'intérêt ZI un champ magnétique homogène Bo ayant une composante axiale orientée selon un axe z 161 du laboratoire. Un dispositif de mesure 140 comprend un carter 143 relié par des 20 éléments de support 142 au boîtier 180. Le carter 143 contient un échantillon 117 de forme allongée orienté selon un axe z' 141 formant un angle 0 par rapport à l'axe z 161 du champ magnétique principal Bo. L'échantillon 117 peut être entraîné en rotation autour de son axe z' (mouvement rotatif 151) par un dispositif 170 d'entraînement en rotation. 25 On a représenté de façon schématique sur la figure 2 des bobines RF 115 (qui entourent le carter 143 et sont coaxiales avec l'échantillon 117 orienté selon l'axe z' 141), des bobines de gradient 114 (qui ont pour axe l'axe z 161 du laboratoire) et des bobines cylindriques 116 de correction de champ magnétique (de rayon a et de longueur 2b) dont les caractéristiques sont déterminées de la façon mentionnée plus haut en prenant en considération un repère Ox'y'z' lié à l'échantillon 117 et qui ont pour axe l'axe z 161 du laboratoire. La figure 6 représente de façon schématique l'ensemble d'un 5 système de spectroscopie et d'imagerie par résonance magnétique auquel est applicable l'invention. Une unité d'expérimentation 101 comprend, à l'intérieur du trou cylindrique d'un aimant 118 de création d'un champ magnétique principal comprenant une composante Bz orientée selon un axe z, de l'extérieur 10 vers l'intérieur, des bobines 116 de correction de champ magnétique coaxiales avec l'axe z), des bobines de gradient 114 (également coaxiales avec l'axe z) et des bobines RF 115 placées au plus près de l'échantillon 117. L'échantillon 117 de forme allongée selon un axe z' est lui-même 15 incliné d'un angle prédéterminé, par exemple l'angle magique, par rapport à l'axe z, de même que les bobines RF 115 qui entourent l'échantillon 117. Une unité d'activation 102 alimente les différentes bobines de l'unité d'expérimentation 101 et reçoit également en retour les signaux RF modulés des bobines RF 115. 20 Une unité de commande 103 qui peut être constituée par un ordinateur comprend un module 136 de communication entre une unité centrale de traitement 139 et l'unité d'activation 102, des unités de mémoires vive 137, des unités de mémoire morte 138 et une interface utilisateur 135. Les valeurs des différents signaux à fournir par l'unité 25 d'activation 102 sont déterminées par l'unité de commande 103. En résumé, l'espace disponible pour une source de champ magnétique aux caractéristiques données est souvent très limité dans certaines directions et conduit à mettre en oeuvre des distributions de courant sur une surface géométrique imposée. En pratique, ces distributions superficielles sont réalisées de façon approchée soit en disposant des conducteurs filiformes sur la surface, soit par des découpes appropriées de feuilles conductrices minces ou par des techniques de circuits imprimés. Par exemple, les machines d'imagerie par résonance magnétique (IRM) doivent être munies de sources de gradients de la composante principale du champ magnétique dans trois directions Gx, Gy, et Gz, aussi homogènes que possible. Dans la plupart des machines, les sources de gradients doivent être placées à l'intérieur du cylindre circulaire libre de l'aimant principal et y occuper le minimum de place, ce qui les confine à un espace cylindrique annulaire de faible épaisseur. On peut les réaliser au moyen de bobinages en fil de cuivre de formes appropriées (en hélices pour Gz et en selles de cheval pour Gx ou Gy) mais également au moyen de tubes (ou de portions de tube) minces en cuivre dans lesquels on pratique des découpes pour créer des canaux de circulation du courant. Dans ces mêmes machines ainsi que dans tous les aimants destinés à produire un champ magnétique très homogène dans une région donnée, il est nécessaire de prévoir des dispositifs correcteurs des imperfections du champ résultant d'une réalisation imparfaite des sources (bobinages 20 parcourus par des courants ou matériaux aimantés) ou de perturbations tant externes à l'aimant qu'internes (échantillon et ses supports). Si ces dispositifs appelés « shims » sont réalisés par des courants, la recherche d'un encombrement minimum conduit également à faire appel à des distributions superficielles générant un champ magnétique dont la 25 composante suivant la direction Oz du champ principal présente un profil donné dans la région d'intérêt. L'invention, qui permet la détermination de densités superficielles de courant portées par un cylindre circulaire générant un profil donné de la composante du champ magnétique suivant l'axe Oz du cylindre, peut REVENDICATIONS1. Procédé de correction de champ magnétique d'un système d'imagerie par résonance magnétique, ledit système comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone d'intérêt ZI, un dispositif de support d'un échantillon (117), avec une orientation d'une dimension principale de l'échantillon (117) selon un angle 90 différent de zéro par rapport à ladite direction Oz, des bobines de gradient (114), et des bobines radiofréquence (115) et la correction étant effectuée à l'aide de bobines (116) de correction disposées autour dudit dispositif de support de l'échantillon (117), caractérisé en ce qu'il comprend les étapes suivantes : - définir un système de coordonnées incliné Ox'y'z' attaché à l'échantillon 15 (117) avec un axe principal Oz' correspondant à ladite dimension principale de l'échantillon (117) orientée selon un angle 00 différent de zéro par rapport à ladite direction Oz du champ magnétique homogène, - déterminer par un développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon (117) la 20 forme des bobines (116) de correction à partir d'une fonction flux F, chaque bobine (116) de correction correspondant à un terme du développement en harmoniques sphériques, - tracer pour chaque bobine (116) de correction des iso-contours de la fonction flux F espacés de manière régulière entre des bornes de la 25 fonction flux F sur un cylindre, et - former à partir des iso-contours du cylindre de chaque bobine (116) de correction des éléments de bobinage pour définir une bobine (116) cylindrique de correction positionnée autour du dispositif de support de 3002 3 2 8 24 l'échantillon (117), l'axe de ladite bobine cylindrique (116) étant confondu avec ladite direction Oz. 2. Procédé selon la revendication 1, caractérisé en ce que la fonction flux F est telle que l'on obtient une distribution de courant minimisant la 5 puissance dissipée pour une valeur d'un terme axial Z'n donné lorsque la fonction Fest solution d'une équation de Poisson. 3. Procédé selon la revendication 1 ou la revendication 2, caractérisé en ce que ledit angle 80 différent de zéro correspond à un angle dit magique égal à arctan (V2) - 54,70. 10 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que l'on définit des première, deuxième et troisième bobines (116) de correction correspondant respectivement à des termes axiaux Z'1, r2, du développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon (117). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que l'on alimente chacune des bobines (116) de correction à l'aide d'un courant dont la valeur est ajustable. 6. Dispositif de correction de champ magnétique d'un système d'imagerie par résonance magnétique, ledit système comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone d'intérêt ZI, un dispositif de support d'un échantillon (117), avec une orientation d'une dimension principale de l'échantillon (117) selon un angle eo différent de zéro par rapport à ladite direction principale Oz, des bobines de gradient (114) et des bobines radiofréquence (115), ledit dispositif de correction de champ magnétique comprenant un ensemble de bobines (116) de correction positionnées autour du dispositif de support de l'échantillon (117), caractérisé en ce que chaque bobine (116) de correction présente un axe confondu avec la direction Oz et comprend des éléments de bobinage réalisés à partir des iso-contours d'une fonctionflux F espacés de manière régulière entre des bornes de la fonction flux F sur un cylindre, la forme des iso-contours étant déterminée à partir d'un développement en harmoniques sphériques dans un système de coordonnées incliné Ox'y'z' attaché à l'échantillon (117) avec un axe principal Oz' correspondant à ladite dimension principale de l'échantillon (117) orientée selon un angle 80 différent de zéro par rapport à ladite direction principale Oz du champ magnétique homogène, chaque bobine (116) de correction correspondant à un terme du développement en harmoniques sphériques. 7. Dispositif selon la revendication 6, caractérisé en ce que ledit angle 80 différent de zéro correspond à un angle dit magique égal à arctan (V2) - 54,7°. 8. Dispositif selon la revendication 6 ou la revendication 7, caractérisé en ce que l'ensemble de bobines (116) de correction comprend des première, deuxième et troisième bobines (116) de correction correspondant respectivement à des termes axiaux Z'1, r2, r3 du développement en harmoniques sphériques dans ledit système de coordonnées incliné Ox'y'z' attaché à l'échantillon (117). 9. Dispositif selon l'une quelconque des revendications 6 à 8, caractérisé 20 en ce que les éléments de bobinage des bobines (116) de correction comprennent des pistes ou fils conducteurs sur un support isolant. 10. Système d'imagerie par résonance magnétique, comprenant un dispositif de création d'un champ magnétique homogène selon une direction Oz, dans une zone d'intérêt ZI, un dispositif de support d'un 25 échantillon (117), avec une orientation d'une dimension principale de l'échantillon (117) selon un angle 80 différent de zéro par rapport à ladite direction Oz, des bobines de gradient (114) et des bobines radiofréquence (115), caractérisé en ce qu'il comprend un dispositif de correction de champ magnétique selon l'une quelconque des revendications 6 à 9. 11. Système selon la revendication 10, caractérisé en ce que ladite direction Oz est verticale. 12. Système selon la revendication 10, caractérisé en ce que ladite direction Oz est horizontale.