FR 2853080 A1 20041001 FR 0303868 A 20030328 Procédé et installation d'imagerie multidimensionnelle par résonance magnétique en champ inhomogène. La présente invention est relative aux procédés 5 d'imagerie multidimensionnelle par résonance magnétique en champ inhomogène. Plus particulièrement, l'invention concerne un procédé d'imagerie multidimensionnelle par résonance magnétique d'un échantillon placé dans un gradient de champ 10 magnétique principalement monodirectionnel le long d'une direction de gradient, dans lequel on soumet l'échantillon à une séquence d'impulsions radiofréquences, et on détecte l'énergie émise par ledit échantillon. En imagerie par résonance magnétique nucléaire, on utilise des gradients de champ magnétique pulsés pour coder la fréquence de résonance des noyaux possédant un spin en fonction de leur position spatiale. Plutôt que de générer un gradient de champ magnétique avec un équipement 20 spécifique, il est connu et appliqué de se servir directement du champ de fuite d'un aimant comme gradient de champ magnétique, par exemple dans le cadre de la technique " STRAFI " (STRAy Field Imaging) présentée dans " Stray Field Imaging of Solids (STRAFI) ", Samilenko, Zick, Bruker 25 Report 1 (1990) 40-41. Néanmoins, dans ce cas, un fort gradient de champ magnétique n'est obtenu que dans une direction. Pour obtenir une image tridimensionnelle de l'échantillon à imager, on doit alors faire pivoter l'échantillon dans chacune des trois directions de 30 l'espace, par exemple mécaniquement à l'aide d'une sonde telle que celle décrite dans le document GB 2 268 275, prendre une image monodirectionnelle de l'échantillon le long de chacune de ces directions, et retrouver par calcul la densité de spin volumique de l'échantillon, ce qui 35 constitue un procédé trop long et fastidieux pour être mis en oeuvre de manière efficace. La présente invention a notamment pour but de fournir un procédé simple et rapide permettant de fournir des images multidimensionnelles (2D ou 3D) d'un échantillon 5 placé dans un fort gradient de champ magnétique monodimensionnel. A cet effet, selon l'invention, un procédé du genre en question est essentiellement caractérisé en ce qu'on fait tourner ledit échantillon autour d'un axe formant un 10 angle non nul avec ladite direction de gradient, de manière à ce que ladite direction de gradient soit successivement orientée, relativement à l'échantillon, le long d'au moins une première et une deuxième directions distinctes d'échantillon, et on synchronise ladite séquence d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence de résonance en fonction de la position le long des au moins première et deuxième directions d'échantillon. Grâce à ces dispositions, on peut coder au cours 20 d'une même période de rotation la fréquence de résonance des noyaux possédant un spin en fonction de leur position spatiale dans un gradient de champ magnétique ayant une direction principale fixe, tel le champ de fuite d'un aimant, et obtenir une image de Fourier multidimensionnelle 25 de l'échantillon. Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes: - on synchronise ladite séquence d'impulsions avec 30 la rotation de l'échantillon en émettant une impulsion radiofréquence au moins quand la première et la deuxième direction d'échantillon sont orientées le long de ladite direction de gradient; - on soumet l'échantillon à une séquence 35 d'impulsions comportant au moins un temps d'acquisition indirecte au cours duquel la première direction d'échantillon est globalement orientée le long de la direction de gradient, et un temps d'acquisition directe, les temps d'acquisition étant séparés par un temps de 5 rotation pendant lequel la rotation de l'échantillon amène ladite deuxième direction d'échantillon le long de la direction de gradient et pendant lequel l'acquisition est stoppée; pendant au moins un temps de rotation, on soumet 10 l'échantillon à au moins une impulsion radiofréquence synchronisée avec la rotation de l'échantillon soit pour stocker soit pour refocaliser l'aimantation de l'échantillon; - on détecte la vitesse de rotation de 15 l'échantillon, et on utilise ladite vitesse détectée pour synchroniser la séquence d'impulsions avec la rotation de l'échantillon; - un signal électrique relatif à la position temporelle détectée de l'échantillon est utilisé pour 20 commander la séquence d'impulsions; - au cours de rotations successives de l'échantillon, on fait varier la durée des temps d'acquisition indirects de manière à coder l'énergie détectée en fonction des temps d'acquisition; - on calcule, à partir de l'énergie détectée codée, une image de Fourier de l'échantillon; - l'échantillon a été placé en une région du champ magnétique dans laquelle le gradient de champ magnétique est globalement constant; - l'acquisition d'image présente une erreur systématique connue liée au positionnement de l'échantillon dans un gradient non constant du champ magnétique et, au cours du calcul de l'image, on redresse l'image à l'aide d'une carte de champ de l'espace; le procédé comprend en outre une étape préliminaire au cours de laquelle on évalue ladite carte de champ en prenant une image de calibration d'un objet de géométrie connue, et en calculant la carte de champ nécessaire à appliquer à l'image de calibration pour représenter ladite géométrie connue; - l'échantillon tourne sensiblement par rapport à la direction de gradient pendant au moins un temps d'acquisition indirecte ou directe, ce qui induit une erreur systématique de l'image et, au cours du calcul de 10 l'image, on redresse numériquement l'image; - on soumet l'échantillon à une séquence d'impulsion comportant: * une première impulsion avant un temps d'acquisition indirecte au cours duquel la première 15 direction d'échantillon est globalement orientée le long de la direction de gradient, * un temps de rotation pendant lequel la rotation de l'échantillon amène la deuxième direction d'échantillon le long de la direction de gradient, et 20 pendant lequel l'acquisition est stoppée, et * un temps d'acquisition directe au cours duquel la deuxième direction d'échantillon est globalement orientée le long de la direction de gradient; - on fait tourner ledit échantillon autour d'un 25 axe formant un angle non nul avec ladite direction de gradient, de manière à ce que ladite direction de gradient soit successivement orientée, relativement à l'échantillon, le long de la première, d'une troisième et de la deuxième directions non coplanaires d'échantillon, et on synchronise ladite séquence d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence de résonance en fonction de la position le long des première, deuxième et troisième directions d'échantillon; - on sélectionne au moins une plage énergétique de 35 l'énergie détectée pour obtenir au moins une coupe de l'échantillon; - on soumet l'échantillon à une séquence d'impulsion comportant: * une impulsion sélective au cours de laquelle 5 la première direction d'échantillon est globalement orientée le long de la direction de gradient, * un premier temps de rotation pendant lequel la rotation de l'échantillon amène la troisième direction d'échantillon le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, * une première impulsion avant un temps d'acquisition indirecte au cours duquel ladite troisième direction d'échantillon est globalement orientée le long de la direction de gradient, 15. un deuxième temps de rotation pendant lequel la rotation de l'échantillon amène la deuxième direction d'échantillon le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, et * un temps d'acquisition directe au cours 20 duquel la deuxième direction d'échantillon est globalement orientée le long de la direction de gradient; - on soumet l'échantillon à une séquence d'impulsion comportant: * une première impulsion avant un temps 25 d'acquisition indirecte au cours duquel la première direction d'échantillon est globalement orientée le long de la direction de gradient, * un temps de rotation pendant lequel la rotation de l'échantillon amène la troisième direction 30 d'échantillon le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, * une deuxième impulsion avant un temps d'acquisition indirecte au cours duquel la troisième direction d'échantillon est globalement orientée le long de 35 la direction de gradient, un temps de rotation pendant lequel la rotation de l'échantillon amène la deuxième direction d'échantillon le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, 5. un temps d'acquisition directe au cours duquel la deuxième direction d'échantillon est globalement orientée le long de la direction de gradient; - les directions d'échantillon sont orthogonales deux à deux; - on soumet l'échantillon à une rotation autour d'un axe formant un angle d'environ 54,70 avec la direction de gradient. Selon encore un autre aspect, l'invention concerne une installation d'imagerie multidimensionnelle par 15 résonance magnétique d'un échantillon placé dans un gradient de champ magnétique principalement monodirectionnel le long d'une direction de gradient, comprenant: une sonde adaptée pour soumettre l'échantillon à 20 une séquence d'impulsions radiofréquences, et pour détecter l'énergie émise par ledit échantillon, un dispositif moteur adapté pour faire tourner ledit échantillon autour d'un axe formant un angle non nul avec ladite direction de gradient, de manière à ce que 25 ladite direction de gradient soit successivement orientée, relativement à l'échantillon, le long d'au moins une première et une deuxième directions distinctes d'échantillon, et une unité centrale adaptée pour synchroniser ladite 30 séquence d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence de résonance en fonction de la position le long des au moins première et deuxième directions d'échantillon. Dans des modes de réalisation préférés de 35 l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes - l'installation comprend un aimant adapté pour générer le champ magnétique présentant un gradient principalement monodirectionnel; - l'échantillon est adapté pour être placé dans une région de champ magnétique homogène, et l'installation comprend en outre au moins une bobine de gradient adaptée pour générer le gradient de champ magnétique le long d'au moins une direction de gradient au niveau de 10 l'échantillon; - la sonde est une sonde MAS (" magic angle spinning ") standard. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description 15 suivante, donnée à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins - les figures la à lc sont des vues schématiques en perspective représentant un mode de réalisation d'une 20 installation, la figure 2 est une vue schématique plane des lignes de champ magnétique à proximité d'un aimant, - les figures 3a et 3b sont des diagrammes représentant très schématiquement l'invention, et - les figures 4a, 4b, 5 et 6 sont des séquences radiofréquences utilisées dans le cadre de l'invention. Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. On décrit ici l'invention en référence à des 30 techniques d'imagerie par résonance magnétique nucléaire, mais l'invention est susceptible d'application dans tout type d'imagerie par résonance magnétique, comme par exemple la résonance paramagnétique électronique, ou autre. Les figures la à lc représentent un premier mode de 35 réalisation de l'invention. On utilise par exemple un aimant 1 supraconducteur présentant principalement la forme d'un cylindre évidé centralement en un espace interne 2, à l'air libre, également principalement cylindrique. L'aimant produit un champ magnétique B de l'ordre par exemple de 5 quelques Tesla à quelques dizaines de Tesla au niveau de l'espace interne 2. Cet aimant est typiquement un électroaimant formé par une bobine 18 de fils supraconducteurs, embobinés autour d'un axe z par exemple vertical, et plongés dans un liquide réfrigérant, tel que par exemple de 10 l'hélium liquide. La bobine, parcourue par un courant électrique i génère le champ magnétique dans l'espace interne 2. Un tel champ magnétique est représenté sur la figure 2. Celle-ci représente les lignes de champ 15 magnétique générées par l'aimant 1 au niveau de l'espace intérieur 2. Le champ magnétique B est relativement uniforme dans une portion centrale 16 de l'espace intérieur, et cette portion centrale 16 est la région dans laquelle sont généralement placés les échantillons à imager 20 en imagerie par résonance magnétique traditionnelle. Dans la technique d'imagerie en champ inhomogène " STRAFI ", on positionne classiquement l'échantillon à imager dans une région excentrée 17, décalée le long de l'axe longitudinal de l'espace interne 12, par rapport à la portion centrale 25 16, par exemple de quelques centimètres. Cette région excentrée présente un fort gradient G de champ magnétique dans la direction de cet axe. On peut estimer par exemple un gradient de 17 Hz/pm pour un champ magnétique central de 11,75 T et un échantillon placé à 7 cm du centre de 30 l'aimant. Il est connu qu'il existe, pour le type d'aimant de la figure 1, une position particulière dans laquelle le gradient G de champ magnétique est quasiment constant en fonction de la côte suivant l'axe z, et on peut éventuellement, mais pas nécessairement, placer l'échantillon dans cette position particulière. A cet effet, on dispose en outre d'une sonde 3 (figure lb), destinée à être introduite dans l'espace interne 2 au sein de l'aimant 1, et présentant à cet effet une forme globalement complémentaire de l'espace interne 2, 5 telle par exemple une sonde " MAS " (" magic angle spinning ") classique. La sonde 3 est ainsi de forme cylindrique allongée le long d'un axe longitudinal qui est confondu avec l'axe z de l'aimant, une fois la sonde insérée dans l'espace interne 2. La sonde 3 présente une 10 portion de réception 4 destinée à recevoir l'échantillon à imager 5, un dispositif moteur 12, un dispositif de détection de la rotation 13, et une ou plusieurs bobines radiofréquences 14. La portion de réception 4 présente par exemple un ou plusieurs éléments de fixation dans lesquels 15 va être monté un porte-échantillon 6 contenant l'échantillon à imager 5. Cet échantillon peut par exemple être un fragment de roche, de tissu biologique, d'aliment ou autre. La figure lc représente en gros plan ce porte20 échantillon. Il s'agit par exemple d'un tube cylindrique creux, par exemple en plastique ou autre, de quelques millimètres de longueur et de diamètre. Le porteéchantillon recevant l'échantillon de manière fixe (par exemple à l'aide d'éléments de maintien placés à 25 l'intérieur du tube cylindrique) s'étend longitudinalement le long d'un axe de porte-échantillon. Il comporte également, par exemple sur une extrémité, une portion de rotor 7 constituée d'une pluralité d'ailettes 8, par exemple en matière plastique, et une portion signal 9, par 30 exemple située en une autre extrémité. La portion signal 9 comporte par exemple sur tout ou partie de sa périphérie une alternance connue de régions claires 10 et de régions foncées 11 lisible par exemple par moyens optiques, ou tout autre moyen permettant de facilement coder la vitesse de 35 rotation du porte- échantillon. Le porte-échantillon 6 est placé dans le logement de réception 4 de sorte que son axe longitudinal forme un angle Gr avec la direction principale du gradient de champ magnétique, et dans notre exemple avec l'axe longitudinal 5 de la sonde 3. On peut par exemple choisir Gm égal à l'" angle magique " classiquement utilisé en résonance magnétique, à savoir un angle d'environ 54,7 , ou tout angle 0 racine de l'équation suivante: 3 cos20 - 1 = 0. Un des intérêts de cette orientation est que des rotations 10 successives de 1200 de l'échantillon autour de cet axe permettent une permutation périodique dans le repère du laboratoire des axes de coordonnées (Xe, Ye, Ze) d'un repère orthonormé lié à l'échantillon. Le porte-échantillon 6 est également disposé de 15 manière à ce que le dispositif moteur 12 puisse provoquer la rotation du porte-échantillon autour de son axe longitudinal. Le dispositif moteur peut classiquement être constitué par une valve de sortie d'air 12 capable de générer un flux d'air, par exemple sur la portion de rotor 20 du porte-échantillon, de manière principalement normale aux ailettes 8. Un tel dispositif est classiquement utilisé en résonance magnétique nucléaire et ne sera pas décrit plus en détail. De plus, n'importe quel dispositif susceptible de générer la rotation de l'échantillon, à une vitesse 25 sensiblement constante de l'ordre de quelques hertz à quelques centaines de hertz, autour de l'axe longitudinal du porte-échantillon peut être utilisé dans le cadre de l'invention. Par exemple, et si les vitesses de rotation utilisées le permettent, on pourra utiliser un moteur pas à 30 pas, ou toute autre technique adaptée. En regard de la portion signal 9 du porteéchantillon, on peut disposer, fixé au corps de la sonde 3, d'un dispositif de détection de rotation 13, tel par exemple un émetteur/détecteur optique capable de détecter 35 une onde optique réfléchie par la portion signal 9 ou tout autre moyen permettant de lire la vitesse de rotation de l'échantillon. Ce dispositif de détection de rotation 13 peut être relié à une unité centrale 15 adaptée pour mesurer une vitesse de rotation du porte-échantillon. 5 L'unité centrale peut par exemple comporter une unité de stabilisation raccordée au dispositif moteur 12 de manière à stabiliser la vitesse de rotation de l'échantillon à une valeur sensiblement constante. La sonde 3 comporte également une ou plusieurs 10 bobines radiofréquences 14 par exemple capable d'exciter les résonances magnétiques nucléaires des différents spins (par exemple de proton) dans les champs magnétiques considérés. Ces bobines, contenant des éléments inductifs et capacitifs, sont disposées pour générer, selon une 15 séquence d'impulsions, un champ magnétique B1 de manière orthogonale au champ magnétique principal Bo. Ces impulsions font tourner l'aimantation globale de l'échantillon. La séquence d'impulsions peut classiquement comporter des impulsions 90 , qui sont dimensionnées de 20 manière à faire tourner le vecteur aimantation nette de l'échantillon de 90 , des impulsions 1800, qui sont dimensionnées de manière à faire tourner le vecteur aimantation nette de l'échantillon de 1800, ou tout autre type d'impulsions nécessaires. Les bobines radiofréquences 25 14 sont éventuellement reliées à l'unité centrale, de manière à recevoir l'information de la vitesse de rotation du porte-échantillon délivrée par le dispositif de détection de rotation 13, et à lui synchroniser la séquence d'impulsions. Ces bobines sont de plus aptes à détecter 30 l'énergie dégagée par la relaxation de l'aimantation de l'échantillon. A des fins d'imagerie, l'échantillon doit être placé dans un gradient de champ magnétique qui permet de coder la fréquence de résonance des noyaux présentant des 35 spins (par exemple des noyaux 1H) en fonction de leur position spatiale. Dans le cadre de la technique " STRAFI ", on utilise le gradient de fuite de l'aimant 1, comme représenté sur la figure 2, plutôt que de générer par un équipement spécifique un tel gradient. Les figures 3a et 3b représentent schématiquement le principe de la réalisation de l'invention. La figure 3a est un diagramme temporel sur une période tr de révolution de l'échantillon, partagé en trois temps de rotation ta et trois temps d'acquisition t1, t2, t3. Pendant le temps 10 d'acquisition t1, l'échantillon est principalement orienté de manière à ce qu'un repère (Xe, Ye, Ze) orthonormé lié à l'échantillon a par exemple son axe Ze orienté dans la direction du gradient G, et son axe xe orienté avec B, (gauche de la figure 3b) . La fréquence de résonance est 15 codée en fonction de la position selon l'axe Ze des noyaux comportant des spins par une impulsion radiofréquence, par exemple de 90 , émise par la bobine radiofréquence juste avant le temps d'acquisition t1. Une impulsion avec un champ radiofréquence d'amplitude 125 kHz durant environ 2 20 ps peut facilement être insérée entre les temps de rotation ta et d'acquisition t1. Pendant le temps d'acquisition ti, qui est relativement court par rapport à la période de rotation de l'échantillon, l'échantillon tourne peu autour de l'axe 25 longitudinal du porte-échantillon. De ce fait, la position spatiale de l'échantillon ne varie que légèrement pendant le temps t1, de sorte que l'échantillon peut être considéré comme immobile. Cette variation de position est d'autant plus faible que le temps d'acquisition est court. On stoppe alors l'acquisition pendant un temps de rotation t,, pendant lequel l'échantillon effectue une rotation autour de l'axe longitudinal du porte-échantillon, de sorte que pendant le temps d'acquisition suivant t2, les positions des axes du repère orthonormé lié à l'objet aient 35 été interchangées. On effectue par exemple une rotation d'environ 120 si l'échantillon est disposé le long d'un axe formant un angle égal à l'" angle magique " avec la direction du gradient G. Par exemple, à l'instant t2, c'est Xe qui est globalement aligné avec B1 et ye avec G, comme 5 représenté au centre de la figure 3b. La fréquence est maintenant en outre codée en fonction de la position selon l'axe Ye des noyaux comportant des spins par une nouvelle impulsion radiofréquence émise avant le temps d'acquisition t2. On coupe alors à nouveau l'acquisition pendant un temps ta nécessaire à amener l'échantillon dans une position dans laquelle c'est au tour de l'axe x, d'être amené parallèle au gradient G (droite de la figure 3b) . La fréquence est maintenant en outre codée en fonction de la 15 position selon l'axe xe des noyaux comportant des spins par une nouvelle impulsion radiofréquence émise avant le temps d'acquisition t3. Puis on stoppe l'acquisition pendant un temps suffisant pour ramener l'échantillon dans la position de 20 départ (gauche de la figure 3b). Pendant les temps de rotation de l'échantillon, l'évolution de l'aimantation nette est bloquée, par exemple, mais pas exclusivement, par des techniques classiques de résonance magnétique nucléaire du solide sous 25 angle magique, telle que l'utilisation de filtres-z, par impulsions d'inversion synchronisées avec la rotation de l'échantillon, ou autres. On peut effectuer un grand nombre de fois cette suite périodique d'opérations pour obtenir une image de 30 Fourier de l'échantillon. A cet effet, pour chaque nouvelle révolution de l'échantillon, on peut être amené à faire varier les temps t1 et t2 pour ne réaliser qu'une acquisition directe pendant t3. Le signal ainsi obtenu est modulé par les variations de t1 et t2, ce qui permet 35 d'obtenir l'image de Fourier de l'échantillon par analyse spectrale du signal obtenu. Il est ainsi essentiel de synchroniser avec précision la rotation de l'échantillon avec l'échéance et la durée des impulsions radiofréquences. La vitesse de 5 rotation de l'échantillon peut être mesurée par le dispositif de détection de rotation 13. Cette donnée est transmise à l'unité centrale, qui peut en conséquence adapter le fonctionnement du dispositif moteur 12, par exemple en régulant le débit de gaz atteignant les ailettes 10 8, pour garantir une vitesse de rotation uniforme. La vitesse de rotation et la phase correspondant à la position de l'échantillon au cours de la rotation peuvent également ou alternativement être transmises à la ou les bobines radiofréquences de sorte que l'émission des impulsions 15 radiofréquences soit contrôlée par la position et/ou la vitesse de rotation de l'échantillon, ou vice-versa. La portion signal 9 du porte-échantillon comporte par exemple une unique zone claire et une unique zone sombre, et le signal correspondant mesuré par le dispositif 20 de détection de rotation est un créneau transmis à l'unité centrale pour déclencher automatiquement la séquence d'impulsions en fonction de la position de l'échantillon. Pendant les différents temps d'acquisition, l'échantillon étant mobile, les axes Xe, Yer ze du repère 25 qui lui est lié tournent légèrement autour de l'axe de rotation de l'échantillon. Pour des temps d'acquisition donnés, on pourra être tenté de réduire la vitesse de rotation, de manière à ce que le déplacement de l'échantillon pendant le temps d'acquisition soit minimal, 30 et donc la résolution spatiale de l'image améliorée. Une vitesse de rotation très lente rallonge néanmoins la durée de la séquence d'impulsions nécessaire pour l'acquisition totale de l'image. Pour certains types d'échantillons, le vecteur aimantation peut posséder en outre un temps de 35 relaxation trop court pour que le signal puisse être mesuré efficacement à de faibles vitesses de rotation. Un compromis peut donc être nécessaire entre la résolution souhaitée et l'intensité du signal mesuré. En particulier dans le cas de temps d'acquisition 5 relativement longs par rapport à la période de rotation, on pourrait obtenir une distorsion de l'image. Cette distorsion constitue toutefois une erreur systématique de l'installation, de sorte qu'elle peut être calculée numériquement en fonction des paramètres d'acquisition et 10 de rotation de l'échantillon. L'image peut alors être corrigée en conséquence. Une autre erreur systématique pouvant conduire à une distorsion de l'image peut être due au fait que le gradient de champ magnétique n'est pas constant à l'endroit 15 o le porte-échantillon est placé. On peut s'en affranchir, par exemple par une étape préliminaire de calibration, qui peut par exemple être effectuée avant chaque image, périodiquement, ou une fois pour toute, par exemple au moment de la mise en service de l'installation. Au cours de 20 cette étape, on peut par exemple détecter l'image d'un fantôme de géométrie connue, tel que par exemple le porteéchantillon lui-même rempli d'eau ou autre. A l'aide de l'image obtenue et de la connaissance de la géométrie de l'échantillon, on peut générer numériquement un champ de 25 déformation à appliquer aux différents pixels de l'image détectée pour obtenir une image fidèle de la géométrie 3D de l'échantillon. Ce champ de déformation est par exemple mémorisé dans l'unité centrale de l'installation qui l'applique à chaque image prise ultérieurement dans les 30 mêmes conditions. La figure 4a décrit une séquence d'impulsions radiofréquence synchronisée avec la rotation de l'échantillon selon l'" angle magique " pour une image tridimensionnelle de l'échantillon selon l'invention. Cette 35 séquence reprend principalement la figure 3a, dans laquelle les temps de détection indirects t, et t2 sont encadrés par deux impulsions radiofréquences 900. La détection directe de l'écho tout entier peut par exemple être effectuée au cours des deux périodes entières suivantes de rotation du 5 porte-échantillon. A cet effet, et de manière classique en résonance magnétique nucléaire, une impulsion 180 est générée après un temps techo sensiblement égal à une période de révolution, et la détection de tout l'écho est effectuée pendant un temps d'acquisition directe t3 consécutif. Le temps d'acquisition techo peut également être choisi très court, par exemple sensiblement plus court que la période de rotation, de manière à ce que l'angle de rotation puisse être négligeable pendant techo. Le signal est alors mesuré en t3 à partir de son sommet apparaissant 15 en un temps techo après l'impulsion 1800. La séquence ne dure alors plus qu'une période de rotation. Pour une rotation lente de l'échantillon, on peut également mettre en oeuvre une séquence telle que représenté sur la figure 4b, au cours de laquelle la détection directe 20 de la moitié du signal est effectuée pendant un temps d'acquisition t3. Il peut exister d'autres techniques de détection pouvant être appliquées dans le cadre de l'invention. Plutôt qu'une image 3D, on peut également 25 facilement acquérir des projections planes de l'échantillon, toujours en le faisant tourner autour d'un axe formant un angle égal à l'" angle magique " avec la direction principale du gradient. On peut par exemple mettre en oeuvre la séquence de la figure 5, représentée 30 uniquement à titre illustratif dans le cadre d'une détection par écho. Après par exemple un premier temps d'acquisition tl, on peut attendre un temps équivalent à la rotation d'environ 1200 de l'échantillon autour de son axe (n=l sur la séquence de la figure 5), ou un temps 35 équivalent à la rotation d'environ 240 autour de son axe (n=2 sur la séquence de la figure 5) pour' effectuer l'acquisition directe. Ainsi on peut obtenir une projection sur le plan (Yeze) pour n=1 et une projection sur le plan (XeZe) pour n=2, avec Ze, Ye, et Xe, les directions de 5 l'objet alignées avec la direction du gradient G aux instants tl, t2, et t3 respectivement. En imagerie bi- ou tridimensionnelle, en introduisant un délai A avant chaque séquence d'acquisition, on peut même obtenir une série de 10 projections de l'échantillon décalées angulairement d'un angle correspondant à l'angle de rotation parcouru par l'échantillon durant le temps A. Cette caractéristique peut être particulièrement utile en imagerie bidimensionnelle pour visualiser l'échantillon sous plusieurs angles. Pour de l'imagerie purement bidimensionnelle, il n'est pas nécessaire de faire tourner l'échantillon autour d'un axe orienté selon l'" angle magique ". Il suffit de faire tourner l'échantillon autour d'un axe de l'espace de sorte que deux directions orthogonales d'un repère lié à 20 l'échantillon soient successivement alignées avec l'axe portant le gradient G de champ magnétique, et de synchroniser de façon adéquate la séquence d'impulsions pour permettre une acquisition aux deux instants auxquels le repère lié à l'échantillon voit une de ses directions 25 alignées avec la direction du gradient. Par exemple, une rotation autour d'un axe perpendiculaire à la direction du gradient est envisageable. On obtient alors une projection bidimensionnelle dans le plan des deux directions de l'échantillon susnommées. Il n'est même pas nécessaire que les deux ou trois directions de l'échantillon forment un repère orthonormé. On peut tout à fait, si on le souhaite, obtenir une image bidimensionnelle (respectivement tridimensionnelle) en faisant tourner l'échantillon de manière à aligner avec la 35 direction du gradient G deux (respectivement trois) directions de l'échantillon formant un repère du plan souhaité (respectivement de l'espace), quitte à reconstruire par calcul l'image dans un plan souhaité (respectivement l'espace) orthonormé. Plutôt que d'obtenir des projections de l'échantillon, on peut également facilement mettre en oeuvre l'invention pour obtenir des coupes à travers l'échantillon, en sélectionnant une tranche de résonance, par exemple en utilisant des impulsions carrées de durée 10 plus longue, par des impulsions sélectives comme représenté sur la figure 6, ou autres. On émet une impulsion sélective, représentée par une Gaussienne sur la figure 6, quand par exemple une direction Ze d'échantillon est orientée avec la direction 15 de gradient de champ magnétique. On peut alors obtenir, dans le plan (xe, ye) des directions d'échantillon alignées avec la direction de gradient aux temps de détection indirecte t2 et directe t3, une coupe correspondant à la tranche de résonance sélectionnée. On peut également augmenter la résolution de l'image en éloignant encore l'échantillon de la portion centrale de l'espace interne de l'aimant vers une région o le gradient est encore plus fort. Il peut alors être nécessaire de réduire la durée des impulsions 25 radiofréquences si on souhaite exciter uniformément toute la gamme des fréquences de résonance. On a ici présenté différentes séquences d'impulsions utilisées pour mettre en oeuvre l'invention. D'autres séquences d'impulsions, non décrites, et pouvant 30 par exemple servir à améliorer la qualité du signal, ou autres, comme il en existe dans l'état de la technique, pourraient également être appliquées dans le cadre de l'invention. L'invention a été présentée en référence à un 35 gradient de champ magnétique généré par le champ de fuite d'un aimant, dans le cadre du procédé " STRAFI ". Il n'est pas forcément nécessaire de travailler dans un tel cadre, l'invention étant applicable à n'importe quel gradient de champ magnétique dirigé principalement selon une direction 5 de gradient donné. On peut par exemple tout à fait utiliser une unique bobine de gradient conventionelle générant un gradient de champ magnétique dans une seule direction pour un échantillon situé dans un champ magnétique globalement homogène. On peut ainsi générer une image multidimensionnelle en n'utilisant qu'une seule bobine de gradient, ce qui est un avantage vis-à-vis des techniques antérieures qui en nécessitent une par dimension d'image. Egalement, l'invention peut tout à fait être mise en oeuvre en bénéficiant du fort gradient de champ 15 magnétique obtenu à la surface d'un aimant ouvert, ce qui permet d'appliquer l'invention à des techniques d'imagerie bi- ou tridimensionnelle ex-situ. La sonde 3 peut par exemple être alors posée directement à la surface de l'aimant ouvert, en une position o le gradient de champ 20 magnétique est adéquat, l'échantillon placé directement à proximité de l'aimant. REVENDICATIONS 1. Procédé d'imagerie multidimensionnelle par résonance magnétique d'un échantillon (5) placé dans un 5 gradient de champ magnétique principalement monodirectionnel le long d'une direction de gradient (z), dans lequel on soumet l'échantillon (5) à une séquence d'impulsions radiofréquences, et on détecte l'énergie émise par ledit échantillon, caractérisé en ce qu'on fait tourner ledit échantillon (5) autour d'un axe formant un angle non nul avec ladite direction de gradient (z), de manière à ce que ladite direction de gradient (z) soit successivement 15 orientée, relativement à l'échantillon, le long d'au moins une première et une deuxième directions distinctes (xe, ye, Ze) d'échantillon, et on synchronise ladite séquence d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence 20 de résonance en fonction de la position le long des au moins première et deuxième directions d'échantillon. 2. Procédé selon la revendication 1, dans lequel on synchronise ladite séquence d'impulsions avec la rotation de l'échantillon en émettant une impulsion radiofréquence 25 au moins quand la première et la deuxième direction d'échantillon sont orientées le long de ladite direction de gradient. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel on soumet 30 l'échantillon à une séquence d'impulsions comportant au moins un temps d'acquisition indirecte (tl, t2) au cours duquel la première direction d'échantillon (Ze, Ye) est globalement orientée le long de la direction de gradient (z), et un temps d'acquisition directe (t3), les temps 35 d'acquisition étant séparés par un temps de rotation (ta) pendant lequel la rotation de l'échantillon amène ladite deuxième direction d'échantillon le long de la direction de gradient et pendant lequel l'acquisition est stoppée. 4. Procédé selon la revendication 3 dans lequel, 5 pendant au moins un temps de rotation, on soumet l'échantillon (5) à au moins une impulsion radiofréquence synchronisée avec la rotation de l'échantillon soit pour stocker soit pour refocaliser l'aimantation de l'échantillon. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel on détecte la vitesse de rotation de l'échantillon, et dans lequel on utilise ladite vitesse détectée pour synchroniser la séquence d'impulsions avec la rotation de l'échantillon. 6. Procédé selon la revendication 5, dans lequel un signal électrique relatif à la position temporelle détectée de l'échantillon (5) est utilisé pour commander la séquence d' impulsions. 7. Procédé selon l'une quelconque des revendications précédentes dans lequel, au cours de rotations successives de l'échantillon, on fait varier la durée des temps d'acquisition indirects (tl, t2) de manière à coder l'énergie détectée en fonction des temps d'acquisition. 8. Procédé selon la revendication 7, dans lequel on calcule, à partir de l'énergie détectée codée, une image de Fourier de l'échantillon. 9. Procédé selon la revendication 8, dans lequel l'échantillon a été placé en une région (17) du champ 30 magnétique dans laquelle le gradient de champ magnétique est globalement constant. lO.Procédé selon la revendication 8, dans lequel l'acquisition d'image présente une erreur systématique connue liée au positionnement de l'échantillon (5) dans un 35 gradient non constant du champ magnétique, et dans lequel, au cours du calcul de l'image, on redresse l'image à l'aide d'une carte de champ de l'espace. ll.Procédé selon la revendication 10, comprenant en outre une étape préliminaire au cours de laquelle on évalue 5 ladite carte de champ en prenant une image de calibration d'un objet (6) de géométrie connue, et en calculant la carte de champ nécessaire à appliquer à l'image de calibration pour représenter ladite géométrie connue. 12.Procédé selon l'une quelconque des revendications 8 à 11, dans lequel l'échantillon tourne sensiblement par rapport à la direction de gradient pendant au moins un temps d'acquisition indirecte (tl, t2) ou directe (t3) , ce qui induit une erreur systématique de l'image, et dans lequel, au cours du calcul de l'image, on 15 redresse numériquement l'image. 13.Procédé d'imagerie bidimensionnelle selon l'une quelconque des revendications précédentes, dans lequel on soumet l'échantillon à une séquence d'impulsion comportant: - une première impulsion avant un temps d'acquisition indirecte (tl, t2) au cours duquel la première direction d'échantillon (ze, Ye) est globalement orientée le long de la direction de gradient (z), - un temps de rotation (ta) pendant lequel la 25 rotation de l'échantillon amène la deuxième direction d'échantillon (xe> le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, et un temps d'acquisition directe (t3) au cours duquel la deuxième direction d'échantillon est globalement 30 orientée le long de la direction de gradient. 14.Procédé selon l'une quelconque des revendications 1 à 12, dans lequel on fait tourner ledit échantillon (5) autour d'un axe formant un angle non nul avec ladite direction de gradient (z), de manière à ce que 35 ladite direction de gradient (z) soit orientée, relativement à l'échantillon, le long de la première (ze)j d'une troisième (ye) et de la deuxième (xe) direction non coplanaires d'échantillon, et dans lequel on synchronise ladite séquence 5 d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence de résonance en fonction de la position le long des première, deuxième et troisième directions d'échantillon. 15.Procédé d'imagerie bidimensionnelle selon la 10 revendication 14, dans lequel on sélectionne au moins une plage énergétique de l'énergie détectée pour obtenir au moins une coupe de l'échantillon (5). 16.Procédé d'imagerie bidimensionnelle selon la revendication 15, dans lequel on soumet l'échantillon à une 15 séquence d'impulsion comportant: une impulsion sélective au cours de laquelle la première direction d'échantillon (Ze) est globalement orientée le long de la direction de gradient (z), - un premier temps de rotation (ta) pendant lequel 20 la rotation de l'échantillon amène la troisième direction d'échantillon (xe) le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, - une première impulsion avant un temps d'acquisition indirecte (t2) au cours duquel ladite 25 troisième direction d'échantillon (xe) est globalement orientée le long de la direction de gradient (z), - un deuxième temps de rotation (ta) pendant lequel la rotation de l'échantillon amène la deuxième direction d'échantillon (Ye) le long de la direction de 30 gradient (z), et pendant lequel l'acquisition est stoppée, et - un temps d'acquisition directe (t3) au cours duquel la deuxième direction d'échantillon est globalement orientée le long de la direction de gradient. 17.Procédé d'imagerie tridimensionnelle selon la revendication 14, dans lequel on soumet l'échantillon à une séquence d'impulsion comportant: une première impulsion avant un temps d'acquisition indirecte (ti) au cours duquel la première 5 direction d'échantillon (Ze) est globalement orientée le long de la direction de gradient (z), - un temps de rotation (ta) pendant lequel la rotation de l'échantillon amène la troisième direction d'échantillon (xe) le long de la direction de gradient, et 10 pendant lequel l'acquisition est stoppée, - une deuxième impulsion avant un temps d'acquisition indirecte (t2) au cours duquel la troisième direction d'échantillon (xe) est globalement orientée le long de la direction de gradient (z), - un temps de rotation (ta) pendant lequel la rotation de l'échantillon amène la deuxième direction d'échantillon (ye) le long de la direction de gradient, et pendant lequel l'acquisition est stoppée, et - un temps d'acquisition directe (t3) au cours duquel la deuxième direction d'échantillon est globalement orientée le long de la direction de gradient. 18.Procédé selon l'une quelconque des revendications précédentes, dans lequel les directions d'échantillon (xe, ye, Ze) sont orthogonales deux à deux. 19.Procédé d'imagerie selon l'une quelconque des revendications précédentes, dans lequel on soumet l'échantillon à une rotation autour d'un axe formant un angle d'environ 54,7 avec la direction de gradient (z). 20.Installation d'imagerie multidimensionnelle par 30 résonance magnétique d'un échantillon (5) placé dans un gradient de champ magnétique principalement monodirectionnel le long d'une direction de gradient, comprenant: une sonde (3) adaptée pour soumettre l'échantillon 35 (5) à une séquence d'impulsions radiofréquences, et pour détecter l'énergie émise par ledit échantillon, un dispositif moteur (12) adapté pour faire tourner ledit échantillon (5) autour d'un axe formant un angle non nul avec ladite direction de gradient, de manière à ce que 5 ladite direction de gradient soit successivement orientée, relativement à l'échantillon, le long d'au moins une première et une deuxième directions distinctes d'échantillon, et une unité centrale adaptée pour synchroniser ladite 10 séquence d'impulsions avec la rotation de l'échantillon de manière à coder la fréquence de résonance en fonction de la position le long des au moins première et deuxième directions d'échantillon. 21.Installation selon la revendication 20 15 comprenant un aimant (1) adapté pour générer le champ magnétique présentant un gradient principalement monodirectionnel. 22.Installation selon la revendication 20, dans laquelle l'échantillon est adapté pour être placé dans une 20 région de champ magnétique homogène, et comprenant en outre au moins une bobine de gradient adaptée pour générer le gradient de champ magnétique le long d'au moins une direction de gradient au niveau de l'échantillon (5). 23.Installation selon l'une quelconque des 25 revendications 20 à 22, dans laquelle la sonde (3) est une sonde MAS standard.