LU81529A1 - Procede et appareil de determination des densites de noyaux par echo de spin pour elaboration de tomogrammes - Google Patents

Description

D. 50.617 -»GRAND-DUCHÉ DE LUXEMBOURG Brevet N°__________^ .....^.......2 || du ...19 j.uillet...l9.7.9 _ (ffîWXty Monsieur le Ministre , . rïvietSj de l’Économie' Nationale et des Classes Moyennes

Titre delivre : ........................................ ÎSK

i Service de la Propriété Industrielle

' LUXEMBOURG

Demande de Breyet d’invention I. Requête .....THE RKCTNTS OF TTTR TTNTVKRfitTY Q.R....CM.XEQSMIA f......22o.O ...Ulllyerslty.....m -„Aweme.,......à ..BERKELEYr ...Etat de ..Calif ornie.,.....Etats-r.U.ixis.....d.î.Âiaër.iq.ue/ ...

.....représenté.....par....Monsieur...Jacqu.es...de..M.uyser..,......agissan.t....en.....qiialir·. (2) ......té...da..mandataire...........................................................................................................................................................................................................

dépose ce .....dixrr.neu£...juillet .. 19 oo....solxante.rr.dlx-jieuf--------------------------(3) à______1.5........... .... heures, au Ministère de l’Économie Nationale et des Classes Moyennes, à Luxembourg : 1. la présente requête pour l’obtention d’un brevet d’invention concernant : ......fî-Procédé.....et.appareil.....de..déterminafciondes....densités.....de....noyaux... (4) _________.par.....é.dtia—de....apla....pour.....ëlaboratlon....de.....tomogrammes"..._______________________________________ déclare, en assumant la responsabilité de cette déclaration, que l’(es) inventeur (s) est (sont) : .....voir .. .au ...verso..................................................................................................................................................................................................................(5) 2. la délégation de pouvoir, datée de........................................................................ le ......................................................................

3. la description en langue.........française.................................. de l’invention en deux exemplaires ; 4.............9................ planches de dessin, en deux exemplaires ; 5. la quittance des taxes versées au Bureau de l’Enregistrement à Luxembourg, \ le 19 juillet....1.9..7..9.....................................................................................................................................-......................................................................

revendique pour la susdite demande de brevet la priorité d’une (des) demande(s) de (6)............................brevet........................................déposée® φ/(Ί)......auX;...Eta.ts-Unis....d..,Amérique................

le 2o.. .juillet .19.7 δ...........(No* 92 6 /-571)...................................................................................................................................(8) au nom de .s.....inventeurs...........................................................................................................................................................................................(9) élit domicile pour lui (elle) et, si désigné, pour son mandataire, à Luxembourg _______________________________________ ......35T.....bld».....Royal............................................................................................................................................................................................................(10) sollicite la délivrance d’un brevet d’invention pour l’objet décrit et représenté dans les annexes susmentionnées, — avec ajournement de cette délivrance à ................6.___________________________mois.

Π. Procès-verbal de Dépôt •La susdite demande de brevet d’invention a été déposée au Ministère de l’Économie Nationale et des Classes Moyennes, Service de la Propriété Industrielle à Luxembourg, en date du :

Gp //tG 1.9 juillet 1979 T\ .........···. Pr. le Ministre fi £ Us> à..........15........... heures / '·. de l’Économie Nationale et des Classes Moyennes, A 68007 '.........“ 1. - Lawrence Eugene CROOKS, 5439 Sacramento Avenue, à RICHMOND,

Conté de Contra Costa, Etat de Californie, Etats-Unis d'Amérique 2. - John Charles HOENNINGER III, 2o Roanoke Road, à BERKELEY, Conté de Alameda, Etat de Californie, Etats-Unis d'Amérique 3. - Mitsuaki ARAKAWA, 112o East Hillsdale Blvd., FOREST CITY, Conté de San Mateo, Etat de Californie, Etats-Unis d'Amérique +.

* y— D 50 bl7

REVENDICATION Pif LA PRIORITE

de la demande de brevet / du modèle d'utilité En ETATS-UNIS D'AMERIQUE Du 20 JUILLET 1978 Mémoire Descriptif déposé à l'appui d'une demande de

BREVET D’INVENTION

au

Luxembourg au nom de: THE REGENTS OF THE UNIVERSITY OF CALIFORNIA.

A.

- pour : " Procédé et appareil de détermination des densités de noyaux par écho de spin pour élaboration de tomogrammes." « * ·» ' φ c -

La présente invention concerne un procédé et un appareil de formation d'une image d'une coupe représentant des densités nucléaires relatives dans un objet. Plus précisément, elle concerne la mise en oeuvre des techniques de résonance magnétique 5 nucléaire pour le relevé in vivo des lignes de densités de noyaux à la résonance, dans un être humain ou un autre animal.

Les techniques actuellement utilisées pour l'examen non traumatisant d'un corps mettent en oeuvre des rayons X (par exemple dans la tomographie assistée par ordinateur ou tacographie 10 et des ultrasons. Un autre procédé non traumatisant met en oeuvre la résonance magnétique nucléaire RMN pour la formation d'une image d'une coupe donnant les densités des noyaux à l'intérieur d'un corps. Le numéro atomique moyen Z des noyaux présents dans les tumeurs a tendance à différer nettement de celui des tissus normaux. 15 Les densités de noyaux d'hydrogëne détectés par résonance magnétique nucléaire sont actuellement considérées comme un indicateur puissant des changements du numéro atomique moyen Z dans les tissus puisque l'hydrogène est le noyau de nombre de masse impair le plus abondant présent dans les tissus vivants. La résonance 20 magnétique nucléaire convient parfaitement au relevé des densités relatives des noyaux d'hydrogène dans .les tissus vivants.

Les principes de la résonance magnétique nucléaire sont bien connus dans leur généralité. Tous les noyaux ayant un nombre impair de protons et de neutrons (c'est-à-dire un nombre de masse 25 impair) se comportent en fait comme de petits aimants. Lorsqu'ils sont disposés dans un champ magnétique externe permanent, l'axe magnétique des noyaux (et en conséquence des atomes) présente une précession correspondant à un certain angle par rapport à l'axe du champ imposé et il apparaît un moment magnétique ou "spin" à 30 une fréquence dite de "Larmor". La fréquence de Larmor f est reliée au champ magnétique Bq au niveau du noyau suivant l'équation f = 2 TT y Bq dans laquelle y est une constante, le rapport gyroma-gnëtique caractéristique d'un type particulier de noyau.

On sait que la direction du moment angulaire résultant 35 ou "spin" d'un groupe de noyaux (et ainsi leur axe magnétique résultant) peut être réorientée par rapport au champ magnétique externe par des signaux électromagnétiques ayant une fréquence égale à la fréquence de Larmor. Le signal électromagnétique provoque 2 * % ' k la nutation (réorientation) du spin résultant des noyaux à la résonance (fréquence de Larmor) d'une quantité correspondant à. son amplitude et sa durée. Le sens de la nutation est fonction de la phase du signal électromagnétique par rapport à la préces-, 5 sion nucléaire autour du champ magnétique imposé.

c Si l'on suppose que les spirs nucléaires sont initiale ment alignés et sont ensuite réorientés transversalement à la direction initiale, les orientations des axes magnétiques des noyaux qui changent provoquent l'induction d'un signal caractéris-10 tique à haute fréquence dans une bobine reliée à un récepteur de signaux à haute fréquence. La repolarisation (nouvel alignement par rapport au champ imposé) se manifeste ensuite avec un temps de relaxation (longitudinale) spin-réseau T1·

Lorsqu'il existe un gradient de champ magnétique dans 15 un objet, ou lorsqu'il existe des défauts d'homogénéité dans le champ magnétique imposé, des noyaux ayant le même rapport gyro-magnétique possèdent des fréquences différentes de Larmor suivant leurs positions dans l'objet. Un volume contenant des noy;aux dans l'objet peut ainsi être considéré comme ayant une plage de fré-20 quences de prêcession, centrée autour d'une fréquence donnée de Larmor.

Il est commode de considérer le procédé nucléaire décrit dans la suite à partir d'un système rotatif de référence qui tourne à la fréquence de Larmor, si bien qu'un moment magnétique 25 nucléaire ayant une précession exactement à la fréquence de T,qrmpr paraît fixe. Dans un tel système de référence, l'aimantation macroscopique M est alignée nominalement sur la direction du champ magnétique imposé Bq (jusqu'à la réorientation) .

Lorsqu'un nouvel alignement nucléaire est effectué, 30 les phases relatives des s pins individuels (la phase étant l'angle formé par la projection du spin sur un plan dans le système rotatif de référence, et un axe, disposé dans ce plan et passant par l'axe de rotation) commencent à diverger lorsque certains noyaux présentent une précession plus rapide et certains une précession 35 plus lente que la fréquence centrale de Larmor. Ainsi, il apparaît un "déphasage" progressif des spins nucléaires individuels et une perte correspondante de la cohérence de phase. Initialement, après réorientation, une tension relativement intense est induite dans * % - JS.

3 les bobines réceptrices et diminue progressivement en amplitude étant donné les échanges d'énergie entre spins (constante de , temps de relaxation spin-spin T2) et le déphasage des spins, ces deux phénomènes étant caractérisés cumulativement par le temps , 5 de relaxation T2 . Ce signal représente la décroissance libre de toute induction ou simplement"décroissance libre*!

On sait qu'un "écho de spin" ou une représentation ultérieure de la décroissance libre peut être créé par remise des spins en cohérence de phase.

10 Par exemple, lorsque, un temps après la réorienta tion des spins nucléaires (par exemple à 90° par rapport à une direction initiale), par une première impulsion électromagnétique de fréquence, d'amplitude et de durée convenables (appelée dans xa suite impulsion à 90°) , un autre signal électromagnétique 15 de fréquence, d'amplitude et de durée convenables est appliqué afin que les spins nucléaires subissent une nutation de 180° (ce signal étant appelé impulsion à 180° dans la suite du présent mémoire), chaque spin individuel tourne en fait de 180° (dans le système rotatif de référence) . Cela signifie que la phase est 20 alors opposée à la phase accumulée avant l'impulsion à 180°.

L'accumulation d'autres déphasagesdes spins nucléaires individuels est la même que précédemment et en conséquence, au temps 2 c (après la perturbation initiale) tous les spins individuels re trouvent à nouveau la cohérence de phase (la phase négative annule i 25 la phase nouvellement accumulée). De cette manière un "écho de spin" de la décroissance libre est créé. L'amplitude de crête de l'écho de spin dépend de la constante de temps de relaxation transversale T2* L'écho de spin est en fait une image spéculaire et l'écho de la décroissance libre centrée sur un temps postérieur du temps 2 30 à la perturbation initiale.

Il faut noter que l'écho de spin présente toujours son maximum à une période suivant l'application de l'impulsion à 180° qui est égale à l'intervalle de temps compris entre l'application de la perturbation initiale (l'impulsion à 90° créant la décroissance 35 libre) et l'application de l'impulsion à 180°. Ce phénomène est appelé dans la suite la "règle des temps égaux".

On peut se reporter à l'ouvrage de Farrar et Becker "Pulse and Fourier Transform NMR Introduction to Theory and Methods", ► ¥ 4

Academie Press, New York, 1971 pour la description plus détaillée des principes fondamentaux de la résonance magnétique nucléaire.

Bien qu'on ait utilisé depuis longtemps les techniques de résonance magnétique nucléaire pour la mesure des champs magné-5 tiques et en analyse chimique, ce n'est que récemment que la technique a été appliquée à la formation d'image dans le domaine médical. De façon générale, les techniques de formation d'image par résonance magnétique nucléaire reposent sur l'hypothèse qu'on dispose volontairement un échantillon dans un 10 champ magnétique variable avec la position (un champ ayant une intensité qui varie en fonction de la position), les fréquences de Larmor des noyaux occupant des positions différentes présentant des différences correspondantes. Ainsi, une discrimination de fréquences est assurée entre les spins des atomes occupant des 15 positions différentes, et la densité de spins dans un volume unitaire ou élémentaire compris dans le volume excité de noyaux est représentée par une composante à une fréquence particulière de la décroissance libre.

Les techniques de formation d'image par résonance magnë-20 tique nucléaire entrent dans l'une des cinq catégories suivantes : la formation d'image à partir de projection, la technique FONAR, la formation d'image par points sensibles, la formation d'image suivant les transformations de Fourier, et la formation d'image * par irradiation sélective.

25 La formation d'image par la technique de projection nécessite la formation d'une multiplicité de projections suivant différentes orientations, par exemple par création d'un gradient linéaire de champ dans l'objet et l'enregistrement d'une projection unidimensionnelle de la densité nucléaire dans la direction 30 définie par le gradient. Une image est alors reconstituée à partir des projections, par des techniques mathématiques analogues à celles qui sont utilisées en tomographie mettant en oeuvre des rayons X, ce procédé étant décrit par exemple dans l'article de Lauterbur, Nature, 242:190, mars. 1973.

35 La technique FONAR met en oeuvre des champs magnétiques de configuration convenable, appliqués transversalement à l'objet afin que seule une petite fenêtre à la résonance, à l'intérieur de l'échantillon, donne un signal de résonance magnétique nucléaire.

"* * 5

La région sensible balaye alors l'objet/ par exemple par déplacement physique. La technique FONAR est décrite par exemple dans l'article de "Focusing Nuclear Magnetic Résonance (FONAR)" Visua-lization of a Tumor in a Live Animal, Science, Vol. 194, p. 1430-5 1432, décembre 1976 et dans le brevet des Etats-Unis d'Amérique n° 3 789 832.

La technique de formation d'image par points sensibles, aussi appelée relevé de spin, est un procédé selon lequel le signal de résonance magnétique nucléaire provenant de volumes élémentaires 10 particuliers est enregistré successivement. Un gradient de champ magnétique, alternant à une basse fréquence prédéterminée (de l'ordre de 50 Hz) est créé suivant un premier axe de l'objet. Le signal de résonance magnétique nucléaire de tous les éléments dans l'objet est alors modulé à la fréquence du changement de gradient, 15 à l'exception des protons se trouvant dans le plan nul du gradient (plan de zéro). Des gradients alternatifs analogues peuvent être appliqués à des fréquences asynchrones suivant les axes transversaux afin qu'un point nul ou noeud soit formé en fait dans l'objet à l'intersection des plans nuis du gradient ou des noeuds. 20 Un filtrage passe-bas convenable donne alors une indication correspondant au signal de résonance magnétique nucléaire provenant du point d'intersection des trois plans nuis. Un balayage de l'objet, suivant une trame, est obtenu par variation des gradients i relatifs. Une telle technique de formation d'image par points 25 sensibles est décrite dans l'article de Hinshaw, Journal of Applied Physics, vol. 47, n° 8, août 1976.

Un procédé par points sensibles multiples, mettant en oeuvre deux gradients alternatifs orthogonaux pour la formation d'une droite nulle et une chaîne d'impulsions cohérentes à haute 30 fréquence, de phase régulièrement répartie et alternant, est suggéré dans l'article de Andrew et Col., "NMR" Images ly Multiple Sensitive Point Method, Application to Larger Biological Systems", Phys. Med. Biol. 1977; vol. 22, n° 5, 971-974, 1977. Cet article indique que la transformation dicrète de Fourier du signal reçu 35 entre les impulsions à. haute fréquence est utilisée pour la formation d'indices relatifs à la densité des protons suivant la droite d'intersection des deux plans de zéro de gradients alternatifs.

Les techniques de formation d'image par transformation * 6 de Fourier mettent en oeuvre en général une impulsion initiale à haute fréquence pour la réorientation des spins des protons dans l’objet, de 90°. Pendant le signal résultant de décroissance libre, l'objet est soumis à des gradients successifs appliqués 5 consécutivement et très rapidement le long de trois axes cartésiens principaux du système. Le signal de décroissance libre est échantillonné en présence du dernier gradient appliqué, et une transformation tridimensionnelle de Fourier est réalisée afin qu'une image tridimensionnelle soit formée. On connaît aussi deux 10 procédés mettant en oeuvre une transformation bidimensionnelle de Fourier. L'article de Kumar et col., "NMR Fourier Zeugmatogra-phy" Journal of Magnetic Résonance 18:69-83 (1975) décrit ces techniques mettant en oeuvre les transformations de Fourier dans le cas de la résonance magnétique nucléaire.

15 Les techniques de formation d'image par irradiation sélective impliquent l'utilisation d'une séquence d'impulsions électromagnétiques ayant des spectres prédéterminés de fréquences Un premier gradient magnétique est appliqué suivant un axe donné et l'objet est irradié par une séquence d'impulsions électromagné-20 tiques ayant un spectre combiné de fréquencesd'intensité égale à toutes les fréquences de Larmor dans l'objet, à l'exception d'une bande étroite. A la suite de 1'irradiation, tous les noyaux présents dans l'objet, à l'exception de ceux de la bande étroite, sont saturés. Les atomes à saturation ne sont donc plus sensibles 25 à d'autres signaux électromagnétiques pendant un temps de l'ordre de la constante de temps de relaxation spin-réseau T^. Le premier gradient magnétique est remplacé par un gradient en direction orthogonale et l'objet est à nouveau irradié par une séquence d'impulsions électromagnétiques ayant cette fois une largeur de 30 bande correspondant à une bande élémentaire particulière dans le plan non saturé. La seconde séquence d'impulsions provoque la nutation des spins des atomes dans la bande prédéterminée d'un angle de 90° avec création d'une décroissance litre. Cette décroissance libre est alors enregistrée en présence d'un gradient magné-35 tique dans la troisième direction orthogonale (suivant la direction de la bande) et une transformation de Fourier est utilisée afin que la répartion de densités nucléaires suivant cette ligne droite soit déterminée. Le brevet des Etats-Unis d'Amérique n° 4 021 726 » * 7 décrit plus en détail la formation d'image par irradiation sélective.

On peut aussi se référer, pour la description des * techniques de formation d'image par résonance magnétique nu-5 cléaire citées précédemment ou différentes aux articles suivants deP.C. Lauterbur et al., "Magnetic Résonance Zeugmato-graphy" 18th Amper. Conf. 1974; de P. Mansfield, P.K. Grannel & A.A. Maudsley, "Diffraction and Microscopy in Solids hy NMR", 18th Amper. Conf. 1974, p. 431-432 ; de P.C. Lauterbur, ''Magnetic 10 Résonance Zeugmatography", de P.C. Lauterbur, "Flow Measurements by NMR Zeugmatography", 24 octobre 1973, de P.C. Lauterbur, "Stable Isotope Distributions ly NMR", Proc. First International Conf. on Stable Isotopes Conf. 730525, 9-18 mai 1973, p. 255-260, de Lauterbur, "Image Formation hy Induced Local Interactions : 15 Exemples Employing Nuclear Magnetic Résonance", Nature, Vol. 242, 16 mars 1973, p. 190-191, de P.C. Lauterbur et col., "ESR Zeugmatography - Distributions of Unpaired Electrons Within Objects",

Gordon Conf. 12-16 août 1974, de P.C. Lauterbur et col., "In Vivo Studies of Cancer ty NMR Zeugmatography", Gordon Conf. 12-16 août 20 1974, de P.C. Lauterbur, "Reconstruction in Zeugmatography - The

Spatial Resolution of Magnetic Résonance Signals", Intl. Workshop on 3-D Image Reconstruction Techniques, 16-19 juillet 1974, de A.N. Garrovay, "Velocity Profile Measurements hy NMR", l8th Amper. Conf. 1974, p. 435-436, de W.S. Hinshaw, "The Application of Time 25 Dépendent Field Gradients to NMR Spin Mapping", 18 th Amper.

Conf. 1974, p. 433-434, de J.M.S. Hutchinson, J.R. Mallard & C.C. Goll, "In Vivo Imaging of Body Structures Using Proton Résonance", I8th Amper. Conf. 1974, p. 283-284, de P. Mansfield & A.A. Maudsley, "Line Scan Proton Spin Imaging in Biological Structures hy NMR", 50 Phys, in Medicine and Biology 21 n° 5 (1976), p. 847-852, de P.K. Grannel, "NMR Body Images", Physics Bulletin, mars 1976, p. 95-96, de P.C. Lauterbur, D.M. Krammer, W.V. House, C. Chen, "Zeugmatography High Resolution NMR Spectroscopy, Images of Chemical Inhomogeneity Within Macroscopie objects", American Chemical Society 35 Journal, 97:23, 12 novembre 1975, de P. Mansfield & P.K. Grannel, "Diffraction and Microscopy in Solids and Liquids by NMR", Physical Review B, vol. 12, n° 9,· 1er novembre 1975, p. 3618-3634, de P. Mansfield, A.A. Maudsley & T. Baines, "Fast Scan Proton Density Imaging by NMR", J. of physics E, vol. 9, 1976, ——-- « 8 p. 271-278, de P.C. Lauterhir, "Bibliography on Magnetic Résonance Zeugmatography", 3 juin 1975, de A.N. Garroway, P.K. Grannel & P. Mansfield, "Image Formation in NMR ly a Sélective Irradiative Process", J. Phys. C : vol. 7, 1974, p. 457-462, de A. Kumar, D.

5 Welt & R. Ernst, "NMR Fourier Zeugmatography", J. Mag. Res. 18, 69-83 (1975), de P. Mansfield & A.A. Maudsley, "Medical Imaging ly NMR", British Journal of Radiology 50, 188-194 (1977), de D.I. Hoult, "Zeugmatography : A criticism of the Concept of a Sélective Pulse in the presence of a Field Gradient", J. Mag.

10 Res. 26, 165-167 (1977), de P. Mansfield & A.A. Maudsley, "Planar Spin Imaging ly NMR", J. of Physics C. vol. 9, 1976, p. L409-412, de P. Mansfield, "Proton Spin Imaging ly Nuclear Magnetic Résonance", Contemporary Physics, vol. B', n° 6, 1976, p. 553-576, de R. Dama-dian et col., "Field Focusmg Nuclear Magnetic Résonance (FONAR) : 15 Visualization of a Tumor in a Live Animal", Science, vol. 194, 24 décembre 1976, p. 1430-1431, E.R. Andrew, "Zeugmatography", IVth Amper. Summer School, septembre 1976, de W.S. Hinshaw, "Image Formation by nuclear Magnetic Résonance : The Sensitive-Point Method", J. of Applied Physics, vol. 47, n° 8, août 1976, 20 de R. Damadian, M. Goldsmith & L. Minkoff, "NMR in Cancer : XVI FONAR Image of the Live Human Body", Physiol. Chem, and Phys. 9, (1977) p. 97-108, de G.N. Holland & P.A. Bottomley, "A Colour Display Technique for NMR Imaging", J. of Physics E : 10 (1977), p. 714-716, de T, Baines & P. Mansfield, "An Improved Picture Dis-25 play for NMR Imaging", Journal of Physics E : Scientific Instruments 9 (1976), p. 809-811, de E.R. Andrew et col. "NMR Images ly the Multiple Sensitive Point method : Application to Larger Bio-logical Systems", Physics in Medicine and Biology 22, n° 5, 971-974 (1977), de L. Minkoff, R. Damadian, T.E. Thomas, N. Hu, M.

30 Goldsmith, J. Koutcher & M. Stanford, "NMR in Cancer : XVII.

Devar for a 53-Inch Superconducting NMR Magnet", Physiol. Chem, and Phys. 9 (1977), p. 101-109, de Ros Herman, "NMR Makes Waves in Medical Equipment Companies", New Scientist, 12 janvier 1978, de L.E. Crooks, T.P. Grover, L. Kaufman & J.R. Singer, "Tomographie 35 Imaging with Nuclear Magnetic Résonance", Investigative Radiology, 13, 63, janvier-février 1978, de W.S. Hinshaw, P.A. Bottomley & G.N. Holland, "Radiographie Thin-Section Image of the Human Wrist fcy Nuclear Magnetic Résonance", Nature, vol. 270, n° 22, ê 9 29 décembre 1977, p. 722-723, et de T.C. Farrar & E.D. Becker, "Pulse and Fourier Trans form NMR - Introduction to Theory and Methods", Académie Press, 1971, New York, p. 1-33.

Les brevets des Etats-Unis d'Amérique n° 3 975 675, 5 4 021 726, 4 015 196, 4 034 191, 3 789 832, 3 932 805, 3 651 396 et 3 999 118, concernent aussi de telles techniques.

Il faut noter que chacune des techniques précitées présente des inconvénients dans différents domaines. Par exemple, les techniques de formation d'une image à partir de projections 10 nécessitent un traitement mathématique très important des données. La technique FONAR nécessite apparemment soit une installation extrêmement complexe pour le talayage du champ magnétique, soit un dispositif quelconque permettant un mouvement relatif entre le champ et le sujet.

15 Les techniques mettant en oeuvre une transformation de

Fourier de type tridimensionnel nécessitent que tous les plans soient balayés simultanément plusieurs fois afin que les données soient suffisante pour que celles qui correspondent aux divers plans puissent être séparées mathématiquement. Dans les techniques 20 mettant en oeuvre une transformation de Fourier de type bidimensionnel, la fréquence de répétition est limitée par le temps de relaxation spin-réseau T1 des noyaux car chaque irradiation affecte la totalité du système des spins. En outre, une capacité importante de mémoire est nécessaire dans un calculateur.

25 Les techniques de formation d'image par irradiation sélective selon lesquelles la totalité de l'objet est saturée à l'exception d'un plan unique présente des inconvénients car les installations correspondantes ne peuvent pas être adaptées à un talayage séquentiel rapide de plusieurs plans. Ainsi., avant qu'un 30 second plan puisse être adressé, un temps suffisant doit s'écouler pour que l'objet cesse d'être saturé.

L'invention concerne une technique d'irradiation sélective de l'objet par des impulsions électromagnétiques permettant la création d'échos de spin (et non la détection des signaux de 35 décroissance libre) pour un balayage rapide et facile de plusieurs plans.

Le phénomène des échos de spin est bien connu comme indiqué précédemment. Cependant, jusqu'à présent, on a utilisé 10 essentiellement l'écho de spin pour la mesure de la constante de temps de relaxation transversale T2 d'un échantillon. Un exemple de système mettant en oeuvre les échos de spin pour la mesure du temps de relaxation T2 lors du relevé des formations terrestres 5 traversées par un trou de forage, est décrit dans le brevet des Etats-Unis d'Amérique n° 3 128 425. De manière analogue, le brevet des Etats-Unis d'Amérique n° 3 213 355 décrit un système de mesure des dimensions d'un récipient, mettant en oeuvre les échos de spin pour la détermination de la constante de temps de 10 relaxation transversale T2· L'article de Mansfield and Maudsley, "Planar Spin Imaging ty NMR", J. Phys. C : Solid State Physics, vol. 9, 1976, cité précédemment, paraît indiquer que, après la décroissance libre, différents arrangement de refocali.ation de signaux 15 (impulsions sélectives de 180 et 90° et diverses combinaisons avec inversion du gradient de champ) peuvent être utilisés pour le rappel du signal en vue de la formation d'une moyenne. Le mécanisme particulier des arrangements de refocalisation n'est cependant pas décrit dans l'article. De manière analogue, le bre-20 vet des Etats-Unis d'Amérique n° 3 781 650 paraît décrire un spectromètre de résonance magnétique nucléaire dans lequel les signaux de décroissance libre et les échos de spin sont combinés en vue de la réduction des interférences.

La communication précitée de Hoult (Journal of Magnetic 25 résonance, 26 : 165-167 (1977)) indique que les techniques d'irradiation sélective violent le"principe d'incertitude" à moins que des défauts de linéarité, présents dans le système de résonance magnétique nucléaire, soient exploités. Cet auteur indique qu'une impulsion sélective rectangulaire comprend un spectre large 30 de fréquences et que, pendant le temps d'application de l'impulsion, les spins perturbés présentent un déphasage. Cependant, cet article indique aussi que cette situation n'est pas irréversible car, si le gradient de champ est inversé après l'impulsion, un écho de celle-ci est formé et, au milieu de l'écho pour les impul-35 sions à petit déphasage, tous les spins sont en phase.

Cette communication de Hoult implique que la configuration de la région choisie est essentiellement identique à celle du spectre d'irradiation sélective . Cependant, on note selon * 11 l'invention que, comme la réponse des spins nucléaires au champ magnétique à haute fréquence n'est pas linéaire, la configuration du volume choisi ne correspond pas exactement à celle du spectre du champ magnétique excitateur à haute fréquence. Par exemple, 5 un spectre ayant un bloc parfaitement rectangulaire de fréquences provoque l'excitation d'un volume qui recouvre une plage de fréquences légèrement plus large que ce bloc de fréquences, et les * bords du volume excité sont inclinés et non verticaux. La confi guration du volume excité peut être calculée à l'aide des équa-10 tions de Block, décrites dans l'article précité de Farrar et Becker, p. 7 & 8, avec un champ magnétique à haute fréquence dépendant du temps de façon appropriée, ayant le spectre de fréquences considéré. L'utilisation des équations de Block détermine aussi le déphasage de spin qui se manifeste pendant l'impulsion 15 à haute fréquence. On peut ainsi montrer que la suggestion de Hoult de l'inversion du gradient peut être utilisée même pour les angles de perturbation qui ne sont pas petits. Pour un angle de perturbation de 90°, l'inversion de gradient nécessaire à l'obtention d'un signal maximal à partir des spins du volume 20 choisi, est un gradient inversé de même intensité que le gradient original mais dont la durée est à peu près égale à la moitié de la durée de l'irradiation. La durée du gradient inverse dépend de la configuration de l'impulsion à haute fréquence. L'effet de ce gradient inverse est de regrouper la plupart des spins qui ont 25 présenté un déphasage pendant l'irradiation sélective. Comme le déphasage pendant cette irradiation n'est pas linéaire, la remise en phase n'est pas parfaite mais elle est importante. Après la fin du gradient inversé, on dispose d'un signal qu'on considère comme étant un signal de décroissance libre bien que Hoult l'appelle 30 "écho". L'autre angle important de perturbation qui est fréquemment utilisé est 180°. Une irradiation sélective pour cette valeur ne nécessite pas de correction de la phase. La raison en est que le déphasage de spin pendant les 90° initiaux de la perturbation est annulé par la remise en phase dans la seconde partie de 90° '35 de la perturbation.

L'inversion du gradient après une irradiation sélective constitue l'un des types de correction de phase nécessaire à la mise en oeuvre des techniques de relevé linéaire décrites dans * 12 le présent mémoire. L'application d'un gradient inverse pendant une période égale à peu près à la moitié de la durée de l'irradiation sélective est appelée, dans la suite du présent mémoire, une correction de phase de type I. En pratique, la surface com-5 prise sous la forme d'onde représentant le gradient de correction en fonction du temps est le facteur primordial. Lorsque le gradient de correction a une intensité double, il doit être appliqué seulement pendant la moitié du temps. Cette caractéristique s'applique à tous les types de corrections de phase et la description 10 n'indique l'exemple des intensités égales que par raison de simplicité. Un second type de correction de phase, appelé correction de type II, est la correction immédiate pour une impulsion de gradient qui vient juste de se terminer ; l'étalement de la phase qui apparaît dans les parties de l'objet qui ne sont pas soumise.. 15 à une irradiation sélective pendant l'application de cette impulsion de gradient, est corrigé par l'application immédiate d'un gradient égal et de sens opposé, pendant le même temps que l'impulsion originale de gradient. Une extension de la correction de phase de type II est l'introduction d'un temps de retard avant 20 application du gradient de correction. Des événements tels que des échos de spin peuvent être observés pendant ce temps de retard. Une telle correction de phase dans laquelle le gradient de polarité opposée est appliqué après un temps de retard, est e appelée correction de phase de type III. Un autre type de correc-25 tion de phase, appelé type IV, est analogue à la correction de type III, mais une impulsion à haute fréquence à 180° est appliquée au volume intéressant pendant la période comprise entre l'application du gradient original et l'application du gradient de correction si bien que la polarité du gradient de correction 30 est la même que celle du premier gradient. Le gradient de correction a la même polarité que le premier gradient étant donné que l'impulsion intermédiaire à haute fréquence de 180° rend les phases négatives.

L'invention concerne un procédé de formation d'image 35 par irradiation sélective dans lequel l'analyse directe de l'écho de spin (différente de l'analyse de la décroissance libre) permet un balayage successif rapide de plusieurs plans (volumes de forme étalée) à l'intérieur de l'objet. Des plans parallèles 13 en nombre prédéterminé sont excités successivement d’une manière sélective afin que les spins des atomes (noyaux) placés à l'intérieur se réorientent d'environ 90°. Des plans transversaux en nombre prédéterminé/ a l'intérieur de l'objet, sont alors 5 excités sélectivement afin que les spins des atomes contenus se réorientent de 180'. Les noyaux placés aux intersections des plans réorientés à 90° et des plansréorientés à 180° produisent ensuite des signaux d'écho de spin à des moments correspondant à la règle des temps égaux. Les signaux d'écho de spin sont 10 ainsi créés pendant une séquence de périodes correspondantes. La mesure des échos de spin en présence d'un champ magnétique variant avec la position le long de la ligne d'intersection permet la détermination des densités de spin des volumes unitaires à chaque intersection individuelle, par transformation de Fourier des 15 signaux correspondants d'écho de spin.

Lorsque les jeux de lectures répétées d'une ligne particulière sont effectués trop souvent, les effets du temps de relaxation réduisent l'intensité des signaux provenant des volumes unitaires. Ainsi, l'intensité du signal indique les effets 20 combinés de la densité de spin et de la constante de relaxation T^. Les volumes élémentaires correspondant à des constantes importantes T1 ont des densités "apparentes" de spin inférieures à la densité normale. Les éléments de volumes ayant des constantes T1 courtes, présentent très peu de variations par rapport à leur 25 densité véritable de spin. Cet effet permet la construction d'une image de constante par formation de plusieurs images de densité "apparente" de spin, par mise en oeuvre de temps de plus en plus courts entre les jeux de lectures répétées. La réduction de la 30 densité apparente en fonction du temps entre les lectures répétées peut être utilisée pour le calcul de la constante de temps de chaque élément de volume. Les constantes de tous les éléments de volume peuvent alors être utilisées pour la formation d'une image de constante de l'objet.

35 D'autres caractéristiques et avantages de l'invention ressortiront mieux de la description qui va suivre, faite en' référence aux dessins annexés sur lesquels : les figures IA à 1D sont des schémas en perspective d'un objet soumis à des séquences de détermination des signaux 14 d'écho de spin par balayage, avec une seule raie, selon l'invention ; la figure 2 est un diagramme des temps illustrant a la séquence de signaux d'écho de spin sélectifs des figures IA 5 à 1D ; la figure 3 est un graphique schématique représentant la séquence de signaux d'écho de spin indiquée en référence aux figures IA à 1D et 2 ; les figures 4A à 4C sont des perspectives schématiques d'un objet soumis à une séquence de balayages donnant des signaux d'écho de spin à plusieurs raies ; la figure 5 est un diagramme des temps représentant une séquence de balayages donnant des signaux d'écho de spin à deux raies ; 15 les figures 6 et 7 sont des graphiques schématiques représentant d'autres séquences de signaux d'écho de spin à plusieurs raies selon l'invention ; les figures 8, 9 et 10 sont des graphiques schématiques représentant des séquences de balayages donnant des signaux 20 d'écho de spin selon l'invention, et mettant en oeuvre une copie ; les figures lia et 11b forment ensemble un diagramme synoptique, constituant la figure 11, d'un appareil destiné à créer les séquences de balayages de signaux d'écho de spin selon l'invention ; 25 la figure 12 est une perspective étalée d'une partie de l'appareil des figures lia et 11b ; et la figure 13 est un- diagramme synoptique d'un démodulateur et d'un circuit générateur de déphasage et synthétiseur de fréquence, du type représenté sur les figures lia et 11b formant 30 la figure 11.

On considère maintenant, en référence aux figures IA à 1D et 2, une séquence de balayage avec formation de signaux d'écho de spin, permettant la détermination de la densité normale ou modifiée de constante des spins nucléaires (appelée dans 35 le présent mémoire "densité de spin") dans un volume élémentaire 11 contenu dans un objet 10. Un appareil qui convient à la mise en oeuvre de cette séquence de balayage de signaux d'écho de spin est décrit dans la suite en référence aux figures lia et 11b, 12 15 4 et 13.

La première étape dans la séquence de balayage, comme dans les autres techniques de formation d'image en résonance magnétique nucléaire, est l'alignement initial des s pins nucléaires 5 dans l'objet 10. A cet effet, un champ magnétique intense Bq est créé par exemple suivant la direction Z d'un système nominal de coordonnées cartésiennes centrées dans l'alignement de l'objet 10 (figure IA). Comme indiqué précédemment, les spins nucléaires des atomes ont des fréquences de Larmor proportionnelles au champ 10 magnétique agissant sur eux et onttendance à s'aligner sur le champ magnétique Bq.

Ensuite, un volume particulier interne à l'objet 10 est excité afin que les spins des atomes compris dans le volume présentent une nutation vers une orientation transversale, de préférence à 90°. Il faut se rappeler que les noyaux ayant des spins à une fréquence donnée de Larmor ne sont sensibles qu'aux signaux électromagnétiques extérieurs pratiquement à cette même fréquence de Larmor. Ainsi, si l'objet 10 est soumis à un champ magnétique qui varie avec la position, par exemple suivant un 20 gradient magnétique dans la direction X, tel que g = SB /Sx, les noyaux placés dans les divers plans Y-Z (volumes étalés), ayant des positions différentes le long de l'axe χ, ont des fréquences différentes de Larmor. Un plan unique Y-Z peut ainsi être adressé par irradiation de l'objet 10 (en présence du champ 25 variant avec la position X), à l'aide d'un signal électromagnétique dont le spectre de fréquences correspond à la fréquence particulière de Larmor de ce plan. Evidemment, en pratique, le volume adressé contenant les atomes a une dimension X finie et en conséquence comprend les spins ayant les fréquences de Larmor 50 d'une largeur de bande prédéterminée.

Ainsi, après alignement initial des spins par le champ Bq, un champ magnétique variant dans la direction X, avantageusement avec un gradient g (figure 2) est créé dans l'objet 10 et donne une différence de fréquences de Larmor suivant l'axe X. En 35 présence du gradient g , un volume étalé particulier 12 (figure IB) à l'intérieur de l'objet 10 est adressé par irradiation de cet objet avec une impulsion électromagnétique X^ à 90°, ayant un spectre de fréquences comprenant pratiquement les composantes correspondant aux fréquences de Larmor des atomes placés dans 16 4 ce volume étalé ou plan 12.

La relation entre l'impulsion X. à 90° et le gradient • 1 g est indiquée sur la figure 2. Il faut noter que cette figure

A

5 (ainsi que la figure 5) représente la présence ou l'absence du gradient et non la configuration ou forme du gradient lui-même.

Pour simplifier la description, on adopte les conventions suivantes. Le processus d'irradiation de l'objet par un signal électromagnétique de spectre prédéterminé de fréquences, 10 en présence d'un champ variant avec la position, pour l'excitation d'un volume prédéterminé de noyaux, est appelé processus d'irradiation sélective dans la suite du présent mémoire. De manière analogue, les volumes étalés dans lesquels les spins nucléaires sont réorientés de 90 ou 180° sont appelés plans à 90° 15 et plans à 180° respectivement.

Comme indiqué précédemment, les spins nucléaires, dans le volume 12, après l'impulsion à 90° et la correction de phase de type I -g , sont initialement alignés et à 90° respec-

X

tivement par rapport à l'orientation originale et induisent une 20 tension relativement intense dans une bobine placée autour de l'axe X. Cependant, la tension induite diminue au cours du temps lorsque les phases des spins, pour les différentes positions dans le volume 12, s'étalent. Le signal de décroissance libre DL créé par les noyaux présents dans le volume 12 à la suite de 25 l'impulsion à 90° X., est représenté sur la figure 2.

« X

Ensuite, une correction de phase de type IV est appliquée au gradient g . Ce gradient est appliqué afin que Z1 les phases s'étalent suivant la direction Z, de manière qu'elles puissent être refocalisées par le gradient g pendant la première Z2 30 moitié de la lecture de l'écho de spin. Le gradient g a la même Z1 surface que la première moitié du gradient g . La figure 2 repré- Z2 sente le gradient g qui ne recouvre pas un autre gradient, pour Z1 une simplification de la représentation du signal de décroissance litre DL. Il est possible que les gradients g_ et -g se re- Z1 35 couvrent, sans effet nuisible sur l'écho de spin, et avec une réduction du temps nécessaire pour la séquence.

L'étape suivante de la séquence de ' balayage d'écho 17 de spin est une rotation de 180° des spins nucléaires dans un volume transversal au volume 12 et recoupant ce dernier de ma-. nière que le volume élémentaire 11 soit commun auxdeux volumes.

Un champ magnétique variant avec la position dans la direction 5 Y, avantageusement un gradient g = S B /5y, est donc créé dans Y z l'objet 10, et un volume étalé ou plan prédéterminé X-Y portant la référence 14 sur la figure IC est excité par une impulsion à 180° ayant un spectre de fréquences qui correspond à la bande ; des fréquences de Larmor des noyaux présents dans le volume 14.

10 Ainsi, les phases des spins des noyaux à l'intérieur du volume 14 sont inversées et le volume 14 devient un plan à 180°.

L'effet de l'inversion de phase sur les noyaux communs aux deux plans 12 et 14 (appelés volume 16 d'intersection sur la figure 1D) est un regroupement des phases des spins. Ainsi, un 15 écho de spin (figure 2) est créé par le volume 16 d'intersection.

Si l'on se réfère brièvement à la figure 3, on note que l'impulsion Y^ à 180° constitue en fait une réflexion du signal 22 de décroissance libre. Lorsque les phases commencent à se regrouper, une réflexion du signal de décroissance est formée 20 et un signal maximal est créé au temps de 2Vît, v étant égal au temps compris entre la création de l'impulsion X-j^ à 90° et la création de l'impulsion X^ à 180°, suivant la règle des temps égaux. Les phases des spins s'étalent ensuite et le signal d'écho de spin décroît comme décroît le signal de décroissance libre, i 25 Comme décrit dans la suite, la règle des temps égaux peut être utilisée dans un système à balayage à plusieurs raies pour la formation de signaux d'écho de spin sans interférence correspondant aux différentes raies. En outre, lorsque le signal d'écho de spin a une amplitude inférieure au signal de décroissance linéaire 30 étant donné la constante de temps de relaxation transversale, l'effet de réflexion donne une période d'échantillonnage double pour la collecte des données.

Le signal d'écho de spin est enregistré en présence d'un champ magnétique variant avec la position, l'intensité du 35 champ magnétique variant en fonction de la position le long de la courbe d'intersection. Comme indiqué en référence aux figures IA à 1D et 2, les densités individuelles de spins des volumes unitaires dans le volume intersection 16 sont déterminées par * 18 > application à l'objet 10 d'un champ magnétique variant avec la position suivant la direction Z, avantageusement avec un gradient

g = <$ B/$ z (figure 2) pendant la période dans laquelle le signal * Z Z

24 d'écho de spin est formé. Le signal d'écho de spin est échan-5 tillonné et la transformée de Fourier du signal échantillonné est déterminée afin que l'intensité des diverses composantes de fréquences du signal d'écho de spin soit mesurée. La densité de spin du volume unitaire particulier 11 est ainsi représentée par l'intensité de la composante de fréquences correspondant 10 à la fréquence particulière de Larmor de ce volume élémentaire.

Lorsque le gradient gz est tel que sa contribution au champ Bz est nulle au centre du volume 16 d'intersection, les volumes élémentaires à des distances Z égales de part et d autre du centre ont des fréquences de Larmor qui correspondent 15 à des différences égales au-dessus et au-dessous de la fréquence centrale de Larmor. Les composantes démodulées du signal provenant de ces volumes élémentaires sont à des fréquences égales et opposées et peuvent donc être séparées par les techniques classiques de détection de quadrature, comme décrit dans la suite. 20 La séquence de balayage de signaux d'écho de spin décrite précédemment est particulièrement avantageuse car elle convient facilement à un balayage rapide à plusieurs raies. L'adressage sélectif d'une séquence de plans à l'aide d'impulsions à 90° puis l'adressage successif de plans transversaux choi-25 sis à l'aide d'impulsions à 180° permettent l'utilisation d'une séquence synchronisée de signaux d'écho de spin provenant de plusieurs raies ou lignes d'intersection. La règle précitée des temps égaux permet une discrimination dans le temps entre les échos respectifs de spins. On considère à titre illustratif 30 une séquence à deux lignes ou raies en référence aux figures 4A à 4C et 5.

Comme dans la séquence de balayage à une seule ligne, un premier plan 12 (figure 4A) est adressé par une impulsion X^ à 90° en coopération avec un gradient X g (figure 5). Ensuite,

. X

35 les opérations successives d'irradiation sélective sont mises en oeuvre afin qu'elles provoquent des invers ions de phase dans les plans transversaux 14A et 14B (figure 4B). Comme indiqué sur la figure 5, un champ magnétique variant avec la position Y, 19 a avantageusement de gradient g^, est créé dans l'objet 10 et une impulsion Y.^ à 90° ayant un spectre convenable de fréquences , est appliquée afin qu'elle excite sélectivement le volume 14A à 180°. Un second plan X-Z 14B est excité par une opération 5 d'irradiation sélective mettant en oeuvre une seconde impulsion Ï2 à 180° et un gradient Y Sy2> suivant une relation temporelle prédéterminée par rapport aux impulsions et Y^. Le plan X-Z 14B est choisi a une position Y différente de celle du plan 14A par décalage convenable du spectre de fréquences de l'impulsion 10 Yj à 189° par une fréquence convenahLe AF, par modification du champ magnétique principal Bq d'une quantité prédéterminée ÀB (le déplacement AB provoque un déplacement de la fréquence relative de Larmor dans l'objet), ou par une combinaison des deux.

Le procédé mettant en oeuvre AB est schématiquement indiqué iur 15 la figure 5.

Les noyaux communs au plan à 90° et au plan à 180°, c'est-à-dire des volumes intersection X·^ et ΧχΥ2 (figure 4C) donnent des signaux respectifs d'ëchcs de spin (aussi désignés par les références X-j^ et X-j^ sur la figure 5) à des moments 20 correspondants à la règle des temps égaux. Les signaux d'échos de spin sont enregistrés en présence de champs magnétiques variant avec la position, de gradients g ^ et gz2 (figure 5) respectivement afin qu'une différence de fréquences de Larmor soit imposée le long de la ligne d'intersection. Les transformées 25 de Fourier des signaux d'échos de spin enregistrés sont alors * effectuées afin que des valeurs des densités de spins, modifiées ou non par la constante Tl7 des volumes unitaires des volumes intersection soient déterminées.

Il faut noter que les échos respectifs de spin Xj^ 30 et X^^2 apparaissent à des périodes qui suivent l'impulsion Y^ et l'impulsion Y2 respectivement, égales aux périodes comprises entre les impulsions Y^ et X.^ et entre les impulsions Y2 et X^. Ainsi, la synchronisation convenable des impulsions Y^ et Y2 permet l'apparition des échos de spin X^Y.^ et X-j^ avec un inter -’ 35 valle voulu de temps si bien que la discrimination entre les raies peut être obtenue.

Les phases des spins ; des noyaux autres que ceux qui se trouvent dans les plans choisis de manière instantanée s'étalent 20

du fait de l'exposition aux gradients magnétiques correspondants. Un étalement parasite de phases dû aux gradients et aux gradins de champ magnétique AB est corrigé par application de différents types de corrections de phase. Par exemple,comme indiqué sur la 5 figure 5, le gradient g est la correction de phase de type I

’ Λ nécessaire après une irradiation sélective à 90°. Le gradient g^ provoque un étalement des phases dans le plan 14B à 180° avant l'impulsion Y g et Ie gradient g^ provoque un étalement des phases dans le plan 14A à 180° après l'impulsion . La cor -10 rection de cet étalement parasite des phases doit être réalisée avant la création d'un signal d'écho de spin a partir des noyaux qui se trouvent dans le plan affecté. Ainsi, pour éviter des erreurs dues à un étalement erroné de phases d'un signal d'écho de spin X-jY.^, ou applique un gradient négatif -g 2 15 juste après la fin du gradient g 2* Cependant, le gradient négatif ~gy2 affecte le plan 14B après l'impulsion Y2 et doit lui-même être corrigé. La correction des effets des:gradients gyl et g^ sur le plan 14B et ainsi l'écho X^2 est assurée par application d'un gradient g g de correction de type IV et d'un 20 gradient g ^ de correction de type III respectivement. La première moitié du gradient de lecture gzl doit refocaliser les phases de spins pour la formation d'un écho. Cette opération nécessite que les spins soient dêfocalisés au début de la séquence. Le gradient gzg remplit cette fonction. La surface due 25 au gradient gzg est égale à la moitié de celle du gradient gzl si bien que la refocalisation est complète au centre de l'écho.

La première moitié du gradient g ^ est en fait une correction de phase de type IV pour le gradient gzg· Ce dernier peut être appliqué en même temps que le gradient -g .

30 Les effets d'étalement des phases de la seconde moitié du gradient g . sur le plan 14B sont corrigés par application d'un gradient négatif -g . (correction de phase de type II). L'écho de spin est refocalisé pendant cette correction et des échantillons prélevés pendant cette période peuvent être utilisés pour la for-35 mation d'une moyenne du signal.

La figure 6 représente schématiquement une séquence de balayage d'écho de spin à trois lignes ou raies. Une impulsion Χχ à 90° est appliquée au temps tg (en présence d'un gradient X) # 21 et une séquence d'impulsions Y (en présence de gradients Y) irradie l'objet à des temps 3ΐ/2, 2t et 5^/2 respectivement.

Ainsi, un écho de spin X-^Y.^ est créé à un temps 3τ par les noyaux disposés en commun dans les plans 90° associés à l'impulsion Υχ.

5 De manière analogue, les échos de spin X^Y2 et X1Y3 sont créés par les noyaux communs aux plans associés aux impulsions X^ et Y2 et aux impulsions X^ et , aux temps 4f et δτ' respectivement.

La correction de phase non représentée peut être réalisée d'une manière . analogue à celle qu'on a décrit en référence à la figure 5. Les 10 échos de spin sont enregistrés en présence d'un gradient Z et les transformées de Fourier sont exécutées pour la détermination des valeurs des densités de spins relatives normales ou modifiées par la constante 1^.

Comme indiqué schématiquement sur la figure 7, plusieurs 15 plans Y -Z peuvent être excités séquentiellement en coopération avec un ou plusieurs plans Y-Z pour la formation d'un balayage à plusieurs lignes ou raies. Par exemple, si l'on suppose des plans Y-Z 12A, 12B et 12C (figure 4A) excités à des orientations à 90° par un processus d'irradiation sélective comprenant l'irra-20 diation de l'objet par les impulsions successives à 90° X^ X2 et Xg en présence de gradients X et cpmme noté précédemment, les plans suivant l'axe X (et suivant l'axe Y) sont déterminés de façon unique par sélection des spectres convenables de fréquences pour les impulsions respectives à 90°, par décalage du champ ma-. 25 gnétique Bq ou par une combinaison de ces deux techniques. Après excitation à 90° des plans Y-Z 12A, 12B et 12C, les plans X-Z choisis 14A et 14B sont excités à 180° par une opération similaire d'irradiation sélective à l'aide d'impulsions à 180° Y^ et Y2·

Les intersections des plans respectifs Y-Z et X-Z provoquent la 30 formation de plusieurs lignes d'intersection comme indiqué sur la figure 4C, Χ.^, Χ2Υχ' Χ3Υ17 X1Y2' X2Y2 et X3Y2* Chac!ue ügne d'intersection crée un signal d'écho de spin repéré sur la figure 7 par la ligne associé d'intersection, suivant la règle précitée des temps égaux.

35 La sélection convenable des intervalles de tempsentre les impulsions permet la création des échos de spin suivant une séquence ne présentant pas d'interférence. Par exemple, si l'on suppose que les impulsions à 90° X^, X2 et X^ sont créées à des 22 temps 0,x et 2τ respectivement et les impulsions à 180° Y^ et Y2 sont créées à des temps lie/4 et 3f respectivement, l'écho de spin provenant de la ligne d'intersection X^Y^ est créé au temps 3-r/2, celui de la ligne X3Y2 est créé au temps 4<r, celui de 5 la ligne X2Yi est cr^e au temps 9ττ/2, celui de la ligne X2Y2 est créé au temps 5-r , celui de la ligne X^Y^ est créé au temps ll-r/2, et celui de la ligne Χ·^Υ2 est créé au temps 6 t. La correction de ' l'étalement de phase dû aux effets des gradients respectifs est à nouveau réalisée comme décrit précédemment en référence à la 10 figure 5. Les échos de spin sont alors enregistrés en présence d'un gradient Z et les transformées de Fourier sont effectuées pour la détermination des valeurs des densités de spins normales ou modifiées par la constante des volumes élémentaires individuels dans les lignes d'intersection. Ainsi, l'utilisation d'une 15 séquence d'échos de spin selon l'invention permet le balayage de nombreuses lignes dans le volume 10 suivant une séquence rapide, à des fréquences qui ne sont pas limitées par le temps de relaxation spin-réseau de l'objet, comme dans les techniques de balayage mettant en oeuvre une saturation.

20 Dans certains cas, il est souhaitable que plusieurs lectures de l'écho de spin soient réalisées pour une ligne donnée d'intersection, afin qu'une moyenne des signaux puisse être réalisée. Comme indiqué précédemment, l'article de Mansfield et * ‘ Maudsley indique que, lorsque la décroissance libre a disparu, 25 différents arrangements de refocalisation de signaux peuvent être utilisés pour le rappel du signal en vue de la formation de la moyenne. Un processus analogue peut être utilisé pour le rappel du signal d'écho de spin. Comme indiqué sur la figure 8, un écho de spin X1Y1 est créé par irradiation sélective initiale d'un 30 premier plan (volume plan) avec une impulsion X^ à 90° puis par irradiation sélective d'un plan transversal avec une impulsion Y^ à 180°. Si l'on suppose que l'intervalle de temps compris entre les impulsions X-j^ et est égal à l'écho de spin Xj^ est créé avec une valeur de crête au temps après l'irradiation • 35 de l'impulsion X1 à 180°, suivant la règle des temps égaux. Lors -qu'une impulsion supplémentaire à 180° Y^, ayant un spectre de fréquences comprenant celui de l'impulsion X-^, est créée au- temps après la valeur de crête de l'écho de spin X.jY , les phases 23 des spins nucléaires créant l'écho de spin X^Y^ sont inversées si tien que les phases convergent et provoquent la création d'une copie ou réplique X^Y^ ' de l'écho de spin X^Y^.

La règle des temps égaux s'applique au phénomène de 5 copie ou réplique. Ainsi, la réplique X^Y^' apparaît un temps ^2 après la création d'une impulsion à 180° y '. D'autres répliques peuvent être formées par application d'autres impulsions à 180°. Cependant, les amplitudes globales des signaux diminuent avec le temps de relaxation transversale T2· Cette décroissance repré -5 10 sente une limite pratique au nombre de répliques qui peuvent être produites.

Il faut noter que, si l'impulsion Y1' à 180° comprend des composantes spectrales qui ne sont pas contenues dans l'impulsion Ylf une inversion de phase est réalisée dans lesparties 15 du plan original à 90° qui ne sont pas affectées par l'impulsion Y1. Ainsi, ces parties du plan original à 90° forment un écho parasite de spin P à un temps 2-r^ + <xr après la création de l'impulsion Y^'. L'écho parasite P peut être évité lorsque le spectre de fréquences de l'impulsion Y1' à 180° est identique à celui 20 de l'impulsion Y^.

Des répliques de plusieurs échos de spin peuvent être créées par impulsion d'une seule impulsion à 180° ayant un spectre de fréquences recouvrant les spectres individuels des * - impulsions originales à 180°, c'est-à-dire recouvrant les plans 25 respectifs à 180°. Les figures 9 et 10 illustrent schématiquement * une telle procédure. Sur la figure 9, un premier plan est irradié sélectivement par l'impulsion X^ à 90° puis des plans transversaux sont irradiés sélectivement successivement par des impulsions à 180° Yj_, Y2 et Y3 respectivement, appliquées à des temps 3t/2, 30 2*t et 5ττ/2. Ainsi, les échos de spin provenant des lignes d'in tersection respectives X.^, χχΥ2 et xiY3 sont créés aux temps 3t, 4t et 5τ' respectivement.

L'application d'une impulsion Y' à 180° à spectre large, par exemple au temps 11^/2, provoque la création de ré-. 35 pliques des échos de spin dans l'ordre inverse. Selon la règle des temps égaux, la réplique Χ,Υ-,' est créée au temps ôt* {<Z/2 après application de l'impulsion Y' à intervalle de temps égal au temps écoulé entre l'apparition de l'écho de spin X-j^ à 5^ « 24 « et l'impulsion Y'). De manière analogue, la réplique Y1Y2 apparaît au temps 7<t· et la réplique X^^1 au temps 8-f.

L'application d'une seconde impulsion à 180° Y" à spectre large par exemple au temps 17t/2 provoque la création 5 d'autres répliques des échos de spin, cette fois dans le même ordre que les échos originaux. Les répliques X.^", X1Y2" et X^" respectivement apparaissent aux temps 9^, 10·^ et ll-f .

Il faut cependant noter qu'une impulsion Y' à 180° à spectre large provoque des inversions de phases dans des parties 10 du plan à 90° qui n'ont pas subi une inversion de phasespar l'impulsion Y^ si tien qu'un écho parasite P de spin est créé en conséquence. Selon la règle des temps égaux, cet écho parasite est créé ä un temps lltr. Ainsi, l'écho parasite P rend pratiquement inutilisable la réplique Y^X-^ ' .

15 L'écho parasite peut être réduit par réduction des parties du plan à 90° qui ne subissent pas l'inversion de phase par l'impulsion Y^. Les effets des échos parasites peuvent aussi être réduits ou étalés entre diverses répliques par modification des synchronisations relatives des impulsions Y' et Y" à 180° à 20 spectre large.

La figure 10 représente une variante de séquence de balayage d'écho de spin mettant en oeuvre une telle copie ou réplique. Dans ce cas, plusieurs plans à 90° sont créés par irradiation sélective avec des impulsions à 90° X^, X2 et X^, 25 suivies par une impulsion à 180° Y. ayant un spectre prédéterminé » * de fréquences. Les lignes résultantes d'intersection créent des échos de spin correspondants Χ^Υ^, X2Yi et xlYi suivant la règle des temps égaux. Ensuite, la copie est réalisée par irradiation sélective à l'aide d'impulsions respectives à 180° Y^' et Y^' 30 de même spectre de fréquences que l'impulsion Y^. Il faut noter que, comme le spectre de fréquences des impulsions à 180° est le même, aucun écho parasite n'est créé.

On se réfère maintenant aux figures lia et 11b et 12 représentant un appareil convenant à la mise en oeuvre des sé-35 quences de balayage d'échos de spin selon l'invention. Les figures lia et 11b représentent sous forme synoptique l'ensemble de l'appareil, alors que l'appareil utilisé pour la création des divers champs magnétiques (à l'exception du circuit de commande) est « 25 représenté en perspective éclatée sur la figure 12.

L'objet 10 est placé dans une bobine 22 à haute fréquence disposée autour de l'axe nominal X d'un système de coordonnées cartésiennes. La bobine 22 à haute fréquence est 5 utilisée pour l'irradiation de l'objet, pour le prélèvement des signaux d'écho de spin créé par l'objet 10 et en outre, elle assure avantageusement le support de l'objet 10.

La bobine 22 est elle-même disposée entre les extrémités 24 et 26 des pièces polaires d'un aimant principal 28.

10 Celui-ci est utilisé pour la création du champ magnétique principal Bq le long de la direction Z du système de coordonnées nominales (transversale à l'axe de la bobine 22) . L'aimant 28 comprend avantageusement un noyau classique de verre autour duquel sont enroulés des conducteurs de cuivre qui coopèrent 15 avec une alimentation non représentée, et qui sont avantageusement refroidis par circulation d'eau afin que la température reste constante.

Des champs variant avec la position Z sont appliqués sélectivement par des bobines de gradient Z 34 et 36 (g ). Ces 20 bobines 34 et 36 sont placées sur les faces respectives des extrémités des pièces polaires 24 et 26 et ont avantageusement la configuration géométrique de bobine de Maxwell. Ainsi, les bobines sont circulaires et concentriques aux extrémités des pièces polaires, le rayon des bobines étant déterminé en fonc-25 tion de la distance qui les sépare. Les bobines 34 et 36 sont montées électriquement en série si bien que les champs magnétiques créés sont en opposition et s'annulent à l'origine du système de coordonnées.

Les champs magnétiques présentant un gradient de posi-30 tion X sont créés sélectivement par les bobines 38 et 40 créant le gradient X (g ). Les bobines 38 et 40 ont avantageusement une

A

configuration rectangulaire et sont placées sur les bouts 24 et 26 des pièces polaires respectivement afin qu'elles simulent des conducteurs infinis placés dans la direction X. Les bobines 35 38 et 48 sont reliées électriquement afin qu'elles créent des champs magnétiques opposés qui s'annulent aussi à l'origine du système de coordonnées.

Les champs magnétiques variant avec la position Y

« 26 sont créés sélectivement par des bobines 42 et 44 au gradient Y (g^.) . Les bobines 42 et 44 ont avantageusement la même configuration et la même surface que les bobines 38 et 40 de gradient X mais sont disposées sur les extrémités des pièces polaires 24 5 et 26 afin qu'elles simulent les conducteurs infinis disposés dans la direction Y. Les bobines 42 et 44 sont aussi reliées électriquement afin qu'elles créent des champs opposés qui 5 s'annulent à l'origine du système de coordonnées.

Des changements du champ principal Bq, indépendant 10 de la position, peuvent aussi être réalisés le cas échéant par des bobines Δβ 32A et 32B disposées autour de la partie externe des bouts des pièces polaires 24 et 26 et elles sont reliées électriquement en série afin qu'elles créent un champ magnétique additif (de façon mutuelle). Le brevet des Etats-Unis d'Amérique 15 n° 4 015 196 décrit plus en détail des bobines convenant à la formation de gradients magnétiques.

La bobine 22 à haute fréquence est reliée électriquement par un circuit convenable 50 d'adaptation d'impédance à une première borne d'un circuit directionnel ou circulateur 52 20 (circuit à raccord magique en T). Le circuit 50 d'adaptation d'impédance comprend avantageusement deux condensateurs variables destinés à assurer l'accord de la bobine 22 à haute fréquence.

Les condensateurs variables sont de préférence en matière non 4 magnétique, par exemple de cuivre ou de laiton, et sont placés 25 avantageusement aussi près que possible de la bobine 22 à haute fréquence.

Le circuit 52 à raccord magique en T assure le couplage sélectif de la bobine 22 à haute fréquence (d'une manière s'excluant mutuellement) à un émetteur 54 et un préamplificateur 56.

30 On sait que le circuit 52 assure la connexion de l'émetteur 54 à la bobine lorsque l'émetteur 54 est en fonctionnement, et relie la bobine 22 au préamplificateur 56 lorsque l'émetteur 54 n'est pas en cours d'émission.

Le préamplificateur 56 est avantageusement de type 35 classique possédant un gain élevé et une grande largeur de bande. Il est souhaitable que le préamplificateur soit commuté à l'état de repos pendant les périodes d'émission ou intermédiaires afin que l'isolement par rapport à l'émetteur 54 soit excellent. Ainsi, 27 on peut utiliser un préamplificateur déclenché. Les signaux de sortie provenant du préamplificateur 56 sont transmis par un détecteur 62 à diodes et un filtre passe-bas 64 le cas échéant en vue du contrôle pendant l'accord du système, et ils sont 5 transmis par un circuit tampon 58 à un démodulateur 60.

Le démodulateur 60 est de préférence un démodulateur en quadrature si bien que non seulement l'information de décalage de fréquences (par rapport à la fréquence centrale comme décrit dans la suite) et la formation d'amplitudes sont transmises mais 10 le signe de décalage de fréquences peut aussi être déterminé et permet la détermination des composantes de fréquences par rapport aux positions dans l'objet 10, de part et d'autre de l'origine.

En outre, l'utilisation d'un démodulateur en quadrature évite les erreurs de phase dues aux retards des circuits, ^.e démodulais teur 60 peut recevoir des signaux représentatifs du signal émis (ainsi que du signal reçu) et il crée des signaux de sortie en phase I et en quadrature Q (le signal de sortie Q est déphasé de 90° par rapport au signal de sortie I). Ces signaux I et Q comprennent des composantes représentatives de la somme et de la 20 différence de la fréquence du signal reçu et de la fréquence du signal émis. Le signal d'écho de spin (moins la fréquence porteuse) est restitué par détermination de la racine carrée de la somme des carrés des signaux I et Q. On décrit plus en détail 4 * un démodulateur convenable 60 en référence à la figure 13.

25 Les signaux de sortie I et Q du démodulateur sont transmis par des filtres passe-bas correspondants 68 à un convertisseur analogique-numérique 70 à deux canaux. Ce convertisseur est lui-même relié par un circuit 72 de couplage à une mémoire à accès direct, à un calculateur convenable 96.

30 L'émetteur 54 est avantageusement un amplificateur de classe A ayant une largeur de bande qui suffit pour couvrir la largeur de bande voulue des fréquences de Larmor, et il est piloté par des signaux provenant d'un modulateur 74. Celui-ci est avantageusement un mélangeur équilibré (avec filtrage) et • 35 il reçoit un signal de réglage de forme d'impulsion d'une commande convenable 76 à micro-ordinateur (comme décrit dans la suite) et un signal à une fréquence porteuse voulue de Larmor (fréquence centrale).

28

Le signal de fréquence porteuse de Larraor est avantageusement crée à partir du signal de sortie d'un oscillateur piézoélectrique 78, par un circuit 80 générateur de déphasage et synthétiseur de fréquence, verrouillés en phase par rapport , 5 à l'oscillateur 78 et un circuit convenable 84 formant filtre passe-bande et circuit de réglage de niveau. Ce filtre 84 crée une onde sinusoïdale ayant une enveloppe prédéterminée constante à partir du signal de sortie du synthétiseur 80 de fréquence. On décrit plus en détail dans la suite en référence à la figure 13 10 un générateur de déphasage et synthétiseur de fréquence 80. La création du signal à la fréquence centrale de Larmor est réalisée en fonction des signaux de commande provenant de la commande 76 à micro-ordinateur.

La commande 76 règle essentiellement la séquence 15 d'événements dans toute l'installation de résonance magnétique nucléaire : le couplage d'un ordinateur 96 avec l'installation, l'affichage des données par exemple avec un terminal 100 à tube à rayons cathodiques, la création des gradients de champ ainsi que la synchronisation, l'amplitude, la fréquence et la phase 20 des signaux électromagnétiques émis. La commande 76 est avantageusement à base d'un micro-ordinateur tel que LSI-11. Un tel micro-ordinateur peut être modifié afin qu'il ait un fonctionnement plus rapide (par exemple par rapport à la commutation i * du signal de porteuse, à la sélection de direction du gradient 25 et à la sélection de la phase) par addition de circuits câblés de couplage d'application particulière, comme décrit dans l'article de "An NMR Sequencer for Imaging", de J. Hoenninger et L. Crooks, actuellement sous presse.

Lors de la création des gradients de champ magnétique, 30 la commande 76 crée des signaux de commande représentatifs des valeurs voulues pour les gradients et pour les sens particuliers des gradients. Un signal de commande analogue représentatif des valeurs particulières de ΔΒ et/ou Δί est créé le cas échéant. Les signaux de réglage de valeurs de gradients et/ou le signal 35 AB sont appliqués à des circuits classiques 90 formant convertisseurs tension-courant et amplificateurs qui eux-mêmes appliquent les signaux de gradients aux bobines convenables g 34, 36,

Z

g„ 38, 40, et g 42 et 44, les signaux AB étant transmis aux x y 29 » bobines 32A et 32B.

La commande 76 transmet des signaux de commande repré-, sentatifs de la forme voulue pour l'impulsion à haute fréquence, au modulateur 74 afin que l'amplitude et durée des signaux élec-5 tromagnétiques soient ainsi réglées, et ainsi l'angle de réorientation du spin (nutation) soit 90 ou 180°, réalisé par le signal. La configuration de l'impulsion électromagnétique et le facteur d'échelle de l'amplitude, ayant les valeurs voulues, sont conservés en mémoire et utilisés sélectivement pour la création des 10 signaux de commande.

Au cours de l'opération d'irradiation sélective, il est souhaitable que les impulsions électromagnétiques aient un spectre à bande étroide de fréquences. Ainsi, on utilise une configuration d'impulsion (sin t)/t (donnant un spectre de fré-15 quences presque rectangulaire). On considère aussi l'utilisation d'une impulsion de forme gaussienne (donnant un spectre gaussien de fréquences).

La commande 76 permet aussi le réglage et l'échantillonnage par le convertisseur analogique-numérique 70 et la com-20 mande de la transmission des données au calculateur 96 par le circuit 72 de couplage. Sous la commande des signaux provenant de la commande 76, le convertisseur 70 prélève un nombre prédéterminé d'échantillons de signaux démodulés et transmet les données au calculateur principal 96 (la mémoire 98) par l'inter-, 25 médiaire du circuit 72. Lorsqu'un groupe d'échantillons est reçu par le calculateur principal 96, il est conservé à des emplacements convenables de la mémoire 98. La programmation de la commande 76 à micro-ordinateur et du calculateur principal 96 est coordonnée de manière que les données reçues à partir 30 du circuit 72 puissent être convenablement interprétées par rapport à la séquence. Ensuite, une transformation de Fourier des données est réalisée, et la transformée de Fourier des balayages respectifs de ligne est affichée sur le dispositif 101 d'affichage à tubes à rayons cathodiques. Le cas échéant, des 35 servo-mécanismes de positionnement de l'échantillon (objet 10) par rapport aux diverses bobines peuvent aussi être commandés.

On se réfère maintenant à la figure 13 pour la description détaillée d'un circuit démodulateur en quadrature con- 30 venable 60 et d'un circuit convenable 80 générateur de déphasage et synthétiseur de fréquence. Ce circuit 80 et le démodulateur 60/ afin qu'ils réduisent au minimum les fuites, fonctionnent avantageusement essentiellement à des fréquences intermédiaires 5 autres que la fréquence réellement émise par l'émetteur 54 et reçue par le préamplificateur 56. A cet égard, on peut se référer au brevet des Etats-Unis d'Amérique n° 3 651 396. Ainsi, l'oscillateur piézoélectrique 78 transmet un signal de sortie de forme rectangulaire, ayant une fréquence prédéterminée (10 MHz) 10 au circuit 80. Dans l'hypothèse où la fréquence voulue de Larmor est de l'ordre de 15 MHz, le signal de l'oscillateur à 10 MHz est alors transmis convenablement à un circuit 102 diviseur de fréquence (par.un facteur 5) afin qu'un signal à 2 MHz soit formé. Ce signal est appliqué convenablement à un circuit 104 15 synthétiseur de fréquence qui est programmable et à une boucle 106 à verrouillage de phase. Le synthétiseur programmable 104 peut recevoir des signaux de commande indiquant la fréquence voulue à la sortie, provenant de la commande 76. Le signal de sortie du synthétiseur 104 est transmis par une porte 108, 20 d'après les signaux de la commande 76 et du filtre passe-bas 110, vers un mélangeur équilibré classique 112.

La boucle 106 à verrouillage de phase comprend avantageusement un oscillateur commandé en tension fonctionnant autour d'une fréquence centrale prédéterminée telle que 52 MHz, 25 et verrouilléeen phase sur le signal à 2 MHz. Le signal de sortie à 52 MHz de la boucle 106 verrouillée en phase est convenablement transmis à un générateur classique 116 de quatre signaux déphasés, créant des signaux à la fréquence intermédiaire de 13 MHz, avec des phases de 0, 90, 180 et 270°. Les signaux de 30 sortie du générateur 116 sont transmis à un circuit sélecteur 118 à porte commandé par les signaux de la commande 76. Le signal choisi à la fréquence intermédiaire est transmis par un filtre passe-bas 120 au mélangeur 112 qui crée une composante d'un signal différent ä la fréquence voulue de la porteuse de Larmor.

' 35 Si l'on suppose que la fréquence voulue est de 15 MHz, le synthétiseur programmable 104 est réglé par la commande 76 afin qu'il transmette un signal à 28 MHz si bien que la composante obtenue par différence, formée par le modulateur 112 (28 - 13 MHz) « 31 « correspond à la fréquence voulue de 15 MHz. Cette composante différence est extraite par le circuit 84 formant filtre passe-bande et le circuit de réglage de niveau afin qu'une onde sinusoïdale d'enveloppe constante soit formée à la frë-5 quence voulue de la porteuse de Larmor.

Le signal de sortie du synthétiseur 104 est aussi transmis à un amplificateur tampon 122 incorporé au démodulateur 60. Le signal traité est transmis par une porte 124/ à un mélangeur équilibré 126 en fonction des signaux de la com-10 mande 76. Ce mélangeur 126 peut recevoir aussi le signal reçu et amplifié par le préamplificateur 56 et le circuit tampon 58. Dans 11 exemple pour lequel la fréquence porteuse émise est de 15 MHz et le signal de sortie du synthétiseur est à 28 MHz, la composante différence du signal de sortie du mélangeur 126 15 est a la fréquence intermédiaire de 13 MHz. La composante différence est extraite du signal de sortie du mélangeur par un filtre passe-bas 128 et parvient à un amplificateur convenable 130 à fréquence intermédiaire, accordé à 13 MHz. Les signaux de sortie de l'amplificateur 130 parviennent à des mélangeurs équi-20 librés 132 et 134. Ceux-ci peuvent recevoir les signaux intermédiaires ayant les phases 0 et 90° provenant du générateur 116 de déphasage incorporé au circuit 80. Ainsi, les mélangeurs 132 et 134 transmettent des signaux à la fréquence acoustique en phase et en quadrature, représentatifs des composantes respectives 25 d'aimantation du signal d'écho de spin.

Comme indiqué précédemment, les signaux de sortie I et Q sont filtrés et échantillonnés puis conservés dans la mémoire 98 du calculateur 96. Celui-ci calcule en fait la somme vectorielle des composantes I et Q et effectue une transformation de Fourier 30 sur la somme vectorielle afin qu'il crée des valeurs correspondant aux amplitudes des composantes respectives des fréquences des signaux des échos de spin.

Il faut noter que, comme le démodulateur 60 et le circuit 80 générateur de déphasage et synthétiseur de fréquence fonction-35 nent essentiellement à des fréquences intermédiaires autres que la fréquence porteuse de Larmor, les fuites de l'émetteur dans le démodulateur 60 sont notablement réduites.

Par exemple, on considère le fonctionnement de l'appa- « 32 reil des figures lia, 11b, 12 & 13 pendant un balayage d'une seule ligne. Les séquences particulières voulues d'irradiation sélective sont introduites (ou rappelées depuis la mémoire) et la séquence est déclenchée après un signal de préparation prove-, 5 nant du calculateur principal 87 et destiné à la commande 76.

L'aimant principal 28 est commandé afin que les spins nucléaires soient alignés initialement.

La commande 76 calcule ou rappelle dans la mémoire ou dans une table de consultation, la configuration d'impulsions, 10 les amplitudes et durées nécessaires aux nutationsde spin de 90 et 180°, les intervalles de temps voulus entre les impulsions et les valeurs voulues pour les gradients.

Les signaux convenables de commande sont créés afin qu'un volume particulier soit irradié sélectivement par une im-15 pulsion à 90° (plan Y-Z 12 de la figure 10). Un signal convenable est appliqué aux bobines g 38 et 40 afin qu'un gradient X soit formé dans l'objet 10. Un signal sous forme d'une impulsion à haute fréquence d'amplitude et durée convenant à l'angle voulu de rotation (c'est-à-dire 90°) est créé et appliqué au rnodula-20 teur 74.

Simultanément, des signaux convenables de commande sont transmis au circuit 80 générateur de déphasage et synthétiseur de fréquence afin qu'un signal de porteuse soit formé à la fréquence de Larmor associée au plan 12. Le synthétiseur 25 programmable 104 de fréquence (figure 13) est chargé par un code convenable de fréquence provenant de la commande 76 à microordinateur, et le mode de création de signaux commence.

Il est parfois souhaitable, pour la formation de la moyenne des signaux, que les spins subissent une nutation sélec-50 tive de l'angle choisi (90, 180°) dans des directions différentes (par exemple que les rotations des spins soient réalisées vers les axes Y ou -Y ou les axes X ou -X). Ainsi, la commande 76 crée aussi un signal convenable de commande transmis à la porte 118 de sélection du circuit 30 générateur et synthétiseur (fi-35 gure 13) afin que le signal à fréquence intermédiaire de phase convenable soit transmis au mélangeur 112, le signal de la porteuse ayant une phase telle que l'impulsion électromagnétique émise assure la rotation dans le sens voulu. La porteuse 108 * 33 (figure 13) est alors fermée afin que le signal à la fréquence de la porteuse soit appliqué au modulateur 74.

Le modulateur 74 transmet donc une impulsion conformée à haute fréquence, à la fréquence voulue de Larmor (avec la 5 phase voulue) à l'émetteur 54. Celui-ci crée en conséquence un signal par l'intermédiaire du circuit 52, du commutateur 53 à haute fréquence et du circuit d'adaptation 50 de 50 ohms, vers la bobine 22 à haute fréquence, et l'objet 10 est ainsi irradié. Ainsi, un plan Y-Z 12 dans l'objet 10 est excité à 90° par une 10 opération d'irradiation sélective.

La commande 76 crée alors des signaux convenables de réglage de gradient de polarité opposée pendant une période prédéterminée, afin qu'un gradient négatif soit créé et permette la correction de la phase, le signal de gradient étant alors 15 supprimé des bobines g 38, 40 si bien que le gradient X est supprime.

Un signal convenable est appliqué aux bobines g 34 et 36 afin qu'un gradient Z soit formé dans l'objet 10. Après un temps prédéterminé, ces signaux sont supprimés, si bien que le 20 gradient Z disparaît. La synchronisation des événements peut être modifiée afin que le gradient Z soit appliqué pendant un temps qui recouvre celui pendant lequel le gradient X utilisé pour la correction de phase est appliqué.

Des signaux convenables de commande sont alors créés 25 afin qu'ils assurent l'irradiation sélective du plan X-Z 14 (figure IC) par une impulsion à 180°. Un signal sous forme d'une impulsion à haute fréquence, ayant l'amplitude et la durée correspondant à un angle de rotation de 180°, est appliqué au modulateur 74 et les signaux convenables de commande sont trans-30 mis au synthétiseur 104 et au sélecteur 116 de déphasage afin que la fréquence et la phase du signal de porteuse soient déterminées. La porte 108 est préparée et le signal modulé (impulsion à 180°) est transmis à la bobine afin que le plan 14 soit irradié sélectivement. Simultanément, un signal convenable est 35 transmis aux bobines g^ 42 et 44 afin qu'un gradient Y soit formé dans l'objet 10. Après l'écoulement du temps approprié, l'émetteur est arrêté et le gradient est supprimé, par l'intermédiaire de la porte 103, suivant le signal donnant la valeur du gradient « 34 transmis par la commande 76.

Après une période correspondant à la règle des temps égaux, la ligne d'intersection 16 (figure 1D) des plans transversaux (12 et 14) provoque la création d'un écho de spin. Juste 5 avant la création prévue de l'écho de spin, des signaux convenables de commande sont utilisés pour la transmission du signal correspondant à la valeur du gradient aux bobines g 34 et 36 si bien qu'un gradient Z est formé dans l'objet 10. La porte 124 est préparée afin qu'elle provoque la mise en fonctionnement du 10 démodulateur.

Les signaux d'écho de spin sont induits dans la bobine 22 et sont transmis par le circuit 52 au préamplificateur 56 puis à un circuit tampon 58 vers le démodulateur 60 en quadrature.

Ce dernier, en coopération avec le filtre passe-bas 68, crée un 15 signal à fréquence acoustique qui est échantillonné sélectivement à une fréquence prédéterminée par le convertisseur analogique-numérique 70 en fonction de signaux provenant de la commande 76.

Les échantillons numérisés sont alors chargés dans la mémoire 98 par l'intermédiaire du circuit 72 de couplage et du calcula-20 teur 96.

La séquence se répète un nombre prédéterminé de fois et la moyenne convenable des signaux est réalisée dans le calculateur 96. Ce dernier assure aussi le calcul des transformées de Fourier et les résultats sont présentés sur le dispositif 101 25 d'affichage à tube à rayons cathodiques.

Lorsqu'un balayage à plusieurs lignes doit être réalisé, les irradiations sélectives successives de plans parallèles sont réalisées soit par modification de la fréquence porteuse, par chargement convenable du synthétiseur programmable 104 avec la 30 fréquence de Larmor +Af, soit par changement approprié du champ magnétique de base par application de signaux convenables aux bobines B 32A et 32B, soit par ces deux opérations.

Il faut noter que l'appareil décrit est un exemple de différents types d'appareils qui peuvent être utilisés pour la • 35 mise en oeuvre des séquences de balayages de signaux d'écho de spin selon l'invention. Il faut noter que d'autres appareils conviennent aussi. Par exemple, l'appareil décrit dans le brevet des Etats-Unis d'Amérique n° 4 021 726 peut être adapté à la » 35 mise en oeuvre d'une séquence de balayage de signaux d'écho de spin. En outre, il faut noter que, bien que les divers conducteurs représentés entre les différents éléments des figures lia, 11b, 12 et 13 soient représentés par des lignes uniques, il ne 5 s'agit pas d'une caractéristique limitative, chaque trait pouvant présenter plusieurs connexions, de manière bien connue dans la technique.

Il est bien entendu que l'invention n'a été décrite et représentée qu'à titre d'exemple préférentiel et qu'on pourra 10 apporter toute équivalence technique dans ses éléments constitutifs sans pour autant sortir de son cadre.

* «·

Claims (11)

1. 36
2. 1. Procédé de détermination des densités relatives des noyaux dans un objet, les noyaux ayant des spins nucléaires aux fréquences de Larmor qui dépendent de l’intensité d'un champ 5 magnétique appliquée, ledit procédé étant caractérisé en ce qu'il comprend la création sélective de signaux d'écho de spin provenant de noyaux disposés dans un volume prédéterminé de l'objet, et la détection de l'intensité des signaux d'écho de spin afin qu'une mesure de la densité normale ou de la densité 10 modifiée par la constante de temps T des noyaux dans le volume prédéterminé soit obtenue.
3. 2. Procédé selon la revendication 1, caractérisé en ce que la création sélective des signaux d'écho de spin par les noyaux placés dans un volume prédéterminé comprend l'excitation 15 sélective d'un premier volume étalé ou plan contenant le volume prédéterminé afin que les spins nucléaires des noyaux contenus dans ce premier volume plan subissent une nutation d'environ 90°, et l'excitation sélective d'un second volume plan, transversal au premier, un intervalle prédéterminé de temps après, 20 le volume prédéterminé se trouvant à l'intersection du premier et du second volume plan, afin que les spins nucléaires des noyaux contenus dans le second volume plan subissent une nutation de 180°, si bien que les noyaux qui se trouvent à l'intersection créent ensuite des signaux d'écho de spin à des moments qui 25 dépendent de l'intervalle prédéterminé de temps.
4. 3. Procédé selon l'une des revendications 1 et 2, caractérisé en ce que la détection de l'intensité des signaux d'écho de spin comprend la réception de ces signaux, la création d'un champ magnétique dans le volume prédéterminé, avec une intensité 30 qui varie en fonction de la position, pendant la réception des signaux d'écho de spin, et la création de valeurs représentant les intensités relatives des composantes à des fréquences différentes, présentes dans les signaux d'écho de spin.
5. 4. Appareil destiné à donner des valeurs relatives aux 35 densités nucléaires relatives normales ou modifiées par la constante de temps , des noyaux qui se trouvent dans un objet, ces noyaux ayant des spins nucléaires ayant des fréquences de Larmor correspondant à l'intensité du champ magnétique appliqué, ledit 37 appareil étant caractérisé en ce qu'il comprend un dispositif d'alignement initial des spins nucléaires des noyaux, un dispositif d'excitation sélective successive des noyaux disposés dans un nombre prédéterminé de premiers volumes sensiblement plans 5 et parallèles qui se trouvent dans l'objet afin que les spins nucléaires des noyaux qui se trouvent dans les premiers volumes sensiblement plans soient réorientés d'un premier angle prédéterminé, un dispositif d'excitation sélective et successive des noyaux qui se trouvent dans ce nombre prédéterminé de seconds 10 volumes sensiblement plans transversaux aux premiers volumes sensiblement plans et recoupant ces premiers volumes sensiblement plans à l'intérieur de l'objet afin que les spins nucléaires des noyaux soient réorientés d'un second angle prédéterminé dans les seconds volumes sensiblement plans, si bien que les spins 15 nucléaires des noyaux qui se trouvent aux intersections respectives des premiers et seconds volumes créent des signaux respectifs d'écho de spin pendant des périodes successives de temps, et un dispositif destiné à recevoir les signaux d'écho de spin et à déterminer à partir de ceux-ci la densité normale ou modi-20 fiée par la constante des noyaux qui se trouvent aux intersections respectives.
6. 5. Appareil selon la revendication 4, caractérisé en ce que le dispositif de réception des signaux d'écho de spin respectifs et de détermination de la densité normale ou modifiée 25 par la constante de temps des noyaux aux intersections respectives, comprend un dispositif destiné à créer, pendant les périodes successives respectives de temps, des champs magnétiques d'intensité variant avec la position dans les intersections respectives, et un dispositif destiné à former la transformée 30 de Fourier des signaux d'écho de spin.
7. 6. Procédé de détermination des densités nucléaires relatives normales ou modifiées par la constante de temps Τι7 à l'intérieur d'un objet comprenant des noyaux ayant des spins nucléaires à des fréquences de Larmor correspondant à l'inten- 35 sité du champ magnétique appliqué, le procédé étant du type qui comprend la création d'une séquence de champs magnétiques variant avec la position afin que l'intensité des champs magnétiques varie successivement dans l'objet en fonction de la posi- i 38 tion dans des directions transversales, l'irradiation de l'objet avec une séquence de signaux ayant des spectres prédéterminés de fréquences, en présence de champs associés variant avec la position, afin que les spins nucléaires des noyaux ayant des 5 fréquences de Larmor comprises dans lesdits spectres de fréquences soient réorientés sélectivement, et la détection des signaux créés par les noyaux, ledit procédé étant caractérisé en ce qu'il comprend la création d'une séquence supplémentaire de champs magnétiques variant avec la position, chaque champ 10 successif étant créé après un premier champ magnétique correspondant variant avec la position, et la variation de l'intensité du champ magnétique dans l'objet dans une direction convenable par rapport à la direction du premier champ correspondant variant avec la position afin que l'étalement parasite des phases dans 15 l'objet dû au premier champ correspondant soit corrigé.
8. 7. Procédé de détermination des densités nucléaires normales ou modifiées par la constante , dans un objet, les noyaux ayant des spins nucléaires à des fréquences de Larmor dépendant de l'intensité du champ magnétique appliqué., ledit 20 procédé étant caractérisé en ce qu'il comprend : la disposition de l'objet dans un premier champ magnétique intense afin que les spins nucléaires des noyaux soient alignés sur une première direction prédéterminée, * ‘ la création, pendant une première période prédéterminée, ‘ 25 d'un premier champ magnétique variant avec la position afin que l'intensité du champ magnétique varie dans l'objet en fonction de la position suivant une seconde direction prédéterminée si bien que les fréquences de Larmor des noyaux varient en fonction de la position relative suivant la seconde position prédéterminée, 30 l'irradiation de l'objet pendant la première période prédéterminée avec un signal ayant un spectre de fréquences comprenant sensiblement des composantes correspondant à une bande prédéterminée de fréquences nucléaires de Larmor, afin que les noyaux ayant des fréquences de Larmor comprises dans cette 35 bande prédéterminée soient excités, les spins nucléaires des noyaux se trouvant dans un premier plan prédéterminé dans l'objet étant réorientés d'un premier angle prédéterminé par rapport à leur orientation, k 39 » la création, pendant des secondes périodes prédéterminées suivant la première période, d'un second champ magnétique variant avec la position afin que le champ magnétique varie dans l'objet suivant la position le long d'une troisième direc-5 tion prédéterminée, si bien que les fréquences de Larmor des noyaux se trouvant dans l'objet varient en fonction de la position relative des noyaux suivant la troisième direction prédéterminée, cette troisième direction étant transversale par rapport à la seconde, 10 l'irradiation de l'objet pendant les secondes périodes prédéterminées, à l'aide de signaux ayant des spectres de fréquences comprenant sensiblement des composantes correspondant aux bandes prédéterminées de fréquences de Larmor afin que les spins nucléaires des noyaux qui se trouvent dans les plans sup-15 plémentaires respectifs se réorientent d'un second angle prédéterminé par rapport ä la première orientation, ces plans supplémentaires ne se recoupant pas mutuellement et recoupant chacun le premier plan dans l'objet suivant des lignes uniques d'inter sections, des signaux d'écho de spin étant ensuite créés pen-20 dant des troisièmes périodes respectives, la création, pendant les troisièmes périodes respectives, de champs magnétiques supplémentaires variant avec la position afin que l'intensité du champ magnétique varie dans l'objet en fonction de la position suivant les lignes d'inter-25 sections, si bien que les fréquences de Larmor des noyaux se trouvant le long des lignes d'intersections varient en fonction de la position relative le long des lignes d'intersections, et la détection des signaux d'écho de spin et la détermination de l'intensité relative des composantes de fréquences 30 de ces signaux, chaque composante de fréquences des signaux d'écho de spin étant représentative de la densité de spin normale ou modifiée par la constante de temps T.^ des noyaux occupant une position associée le long des lignes d'intersections.
9. 8. Procédé de détermination des densités relatives nor- 35 males ou modifiées par la constante de temps de noyaux placés dans un objet, les noyaux ayant des spins nucléaires à des fréquences de Larmor qui dépendent de l'intensité du champ magnétique appliqué , ledit procédé étant caractérisé en ce qu'il 40 comprend : la disposition de l'objet dans un premier champ magnétique intense afin que les spins nucléaires des noyaux soient alignés suivant une première direction prédéterminée, 5 la création, pour une première période prédéterminée, d'un premier champ magnétique variant avec la position afin que l'intensité du champ magnétique varie dans l'objet suivant la position le long d'une seconde direction prédéterminée, les fréquences de Larmor des noyaux variant en fonction de la posi-10 tion relative le long de la seconde direction prédéterminée, l'irradiation de l'objet pendant la première période prédéterminée avec un signal ayant un spectre de fréquences comprenant sensiblement des composantes correspondant à une bande prédéterminée de fréquences nucléaires de Larmor, afin que des 15 noyaux ayant des fréquences de Larmor dans la bande prédéterminée soient excitées sélectivement, les spins nucléaires des noyaux disposés dans le premier plan prédéterminé dans l'objet étant réorientés d'un premier angle prédéterminé par rapport à leur orientation, 20 la création, pendant des secondes périodes prédéterminées successives, postérieures à la première période, d'un second champ magnétique variant avec la position afin que le champ magnétique varie dans l'objet en fonction de la position le long * ’ d'une troisième direction prédéterminée, les fréquences de Larmor *' 25 des noyaux de l'objet variant en fonction de la position relative des noyaux le long de la troisième direction prédéterminée, cette troisième direction prédéterminée étant transversale par rapport à la seconde direction prédéterminée, l'irradiation de l'objet pendant les secondes périodes 30 prédéterminées successives avec des signaux ayant un spectre de fréquences comprenant sensiblement des composantes correspondant à une bande prédéterminée de fréquences de Larmor afin que les spins nucléaires des noyaux ayant des fréquences de Larmor comprises dans la bande prédéterminée de fréquences de Larmor - 35 se réorientent d'un second angle prédéterminé par rapport à leur première orientation, la modification cumulée de l'intensité du premier champ magnétique d'une quantité prédéterminée dans chaque * 41 seconde période prédéterminée afin que les noyaux particuliers ayant des fréquences de Larmor comprises dans la bande prédéterminée de fréquences de Larmor soient modifiés pendant chaque seconde période prédéterminée, si bien que les spins nucléaires, 5 dans des plans supplémentaires à l'intérieur de l'objet se réorientent du second angle prédéterminé, les plans supplémentaires recoupant le premier plan suivant des lignes uniques * d'intersection, les noyaux placés autour des lignes d'intersections créant des signaux d'écho de spin ultérieurement pendant 10 des troisièmes périodes successives, la création, pendant les troisièmes périodes, de champs magnétiques supplémentaires variant avec la position afin que l'intensité du champ magnétique varie dans l'objet suivant la position le long des lignes d'intersections, si bien que les 15 fréquences de Larmor des noyaux se trouvant le long des lignes d'intersectionsvarient avec la position relative le long de la ligne d'intersections, et la détection des signaux d'écho de spin et la détermination de l'intensité relative des composantes de fréquences de 20 ces signaux, chacune des composantes de fréquences des signaux d'écho de spin étant représentative de la densité normale ou modifiée par la constante de temps des noyaux occupant une position associée le long de la ligne d'intersections créant « l'écho de spin. , 25 9. Procédé selon la revendication 7, caractérisé en ce qu'il comprend : la création, pendant des quatrièmes périodes suivant les troisièmes périodes, d'un quatrième champ magnétique variant avec la position afin que l'intensité du champ magnétique variant 30 dans l'objet suivant la position le long de la troisième direction prédéterminée, l'irradiation de l'objet pendant une partie au moins de la quatrième période prédéterminée avec un troisième signal ayant un spectre de fréquences qui comprend toutes les fréquences 35 des irradiations des secondes périodes afin que des signaux de réplique des signaux d'écho de spin soient créés ultérieurement, pendant des cinquièmes périodes prédéterminées respectives, la création, pendant les cinquièmes périodes prédéter- £ 42 minées, d'un quatrième champ magnétique variant avec la position afin que l'intensité du champ magnétique varie dans l'objet suivant la position le long de la ligne d'intersection, a la détection des répliques, et . 5 le traitement des signaux d'écho de spin et des ré pliques en combinaison afin que le signal traité obtenu soit représentatif de la densité des noyaux.
10. 10. Appareil destiné à donner des valeurs correspondant aux densités nucléaires relatives dans des volumes élémentaires 10 d'un objet, les noyaux ayant des spins nucléaires à -des fréquences de Larmor dépendant de 1'intensité du champ magnétique appliqué, ledit appareil étant caractérisé en ce qu'il comprend : un dipositif commandé par des signaux de commande qu'il reçoit et destiné à créé un premier champ magnétique intense 15 dans l'objet, afin que les spins des noyaux s'alignent dans une première direction prédéterminée, un premier dispositif générateur, pendant une première période prédéterminée, d'un premier champ magnétique variant avec la position, faisant varier l'intensité du champ magnétique 20 dans l'objet en fonction de la position dans la seconde direction prédéterminée si bien que les fréquences de Larmor des spins nucléaires varient en fonction des positions relatives des noyaux suivant la seconde direction prédéterminée, un second dispositif d'irradiation de l'objet pendant 25 la première période prédéterminée avec un signal ayant un spectre de fréquences comprenant sensiblement des composantes correspondant à une bande prédéterminée de fréquences de Larmor afin que les noyaux ayant des fréquences de Larmor comprises dans cette bande prédéterminée soient excités sélectivement, les spins des 30 noyaux se trouvant dans un premier plan prédéterminé dans l'objet se réorientant d'un premier angle prédéterminé par rapport à leur orientation, un troisième dispositif générateur, pendant des secondes périodes prédéterminées séquentielles respectives suivant la 35 première période, d'un second champ magnétique variant avec la position, faisant varier le champ magnétique dans l'objet d'après la position suivant une troisième direction prédéterminée si bien que les fréquences de Larmor des noyaux présents dans l'objet fc 43 4 varient en fonction de la position relative des noyaux dans la troisième direction prédéterminée/ cette troisième direction étant transversale par rapport à la seconde, un quatrième dispositif destiné à irradier l'objet 5 pendant les secondes périodes prédéterminées successives avec des signaux ayant un spectre de fréquences qui comprend essentiellement les composantes de fréquences correspondant à une bande ' prédéterminée de fréquences de Larmor afin que les spins des noyaux ayant des fréquences de Larmor comprises dans cette bande 10 prédéterminée soient réorientés d'un second angle prédéterminé par rapport à leur première orientation, un cinquième dispositif générateur de signaux de commande destiné au premier dispositif et destiné à faire varier sélectivement l'intensité du premier champ magnétique, par des 15 quantités prédéterminées, pendant chaque seconde période prédéterminée afin que les noyaux particuliers ayant des fréquences de Larmor dans la bande prédéterminée soient changés pendant chaque seconde période prédéterminée, des plans supplémentaires de l'objet étant ainsi réorientés du second angle prédéterminé, 20 ces plans supplémentaires recoupant le premier plan suivant des lignes uniques d'intersection, les noyaux placés le long des lignes d'intersection créant ensuite des signaux d'écho de spin pendant des troisièmes périodes successives respectives, un sixième dispositif générateur, pendant les troi- « 25 sièmes périodes respectives, de champs magnétiques supplémen-taires variant avec la position, faisant varier l'intensité du champ magnétique dans l'objet suivant la position le long des lignes d'intersection si bien que les fréquences de Larmor des noyaux se trouvant le long de ces lignes d'intersection 30 varient en fonction de la position relative le long de la ligne d'intersection, un septième dispositif destiné à détecter les signaux d'écho de spin et à déterminer les intensités relatives des composantes de fréquences de ces signaux, chacune des composantes 35 de fréquences des signaux d'écho de spin étant représentative de la densité des atomes occupant une position associée le long de la ligne d'intersection créant l'écho de spin, et un huitième dispositif coopérant avec le premier, le 44 troisième et le cinquième dispositif et destiné à faire varier sélectivement l'intensité du champ magnétique dans l'objet suivant la position le long de la seconde et de la troisième direction et le long des lignes d'intersection, pendant 5 diverses périodes prédéterminées afin que l'étalement parasite de la phase des signaux de spin soit corrigé.
11. 11. Appareil de détermination des densités nucléaires relatives dans des volumes élémentaires d'un objet, les noyaux ayant des spins nucléaires à des fréquences de Larmor dépendant 10 de l'intensité du champ magnétique appliqué, les noyaux étant sensibles uniquement à des signaux à leur fréquence de Larmor, l'objet étant considéré par rapport à un système nominal de coordonnées cartésiennes ayant des directions orthogonales X, Y et Z, ledit appareil étant caractérisé en ce qu'il comprend : 15 un premier dispositif commandé par des signaux de commande qui lui sont appliqués et destiné à créer de façon réglée un champ magnétique indépendant de la position, dans l'objet et le long de la direction Z, un second dispositif commandé par des signaux de 20 commande qui lui sont appliqués et destiné à créer de façon réglée des champs magnétiques respectifs variant avec la position, dans l'objet, si bien que l'intensité du champ magnétique dans la direction Z, dans l'objet, varie en fonction de la position X, „ Y et Z respectivement, avec variation des fréquences de Larmor » 25 des spins nucléaires dans l'objet en fonction des positions X, Y et Z respectivement, un troisième dispositif commandé par des signaux de commande qu'il reçoit et destiné à créer un signal de porteuse à une fréquence prédéterminée de porteuse, 30 un quatrième dispositif commandé par des signaux de commande qu'il reçoit et destiné à créer de façon réglée un signal de modulation ayant une forme d'onde, une amplitude et une durée prédéterminées et choisies, un cinquième dispositif commandé par le signal de modu-35 lation et le signal de porteuse et destiné à moduler ce dernier par le signal de modulation et à irradier l'objet avec le signal modulé, un sixième dispositif destiné à créer un premier jeu * 45 de signaux de commande destinés au second, au troisième et au quatrième dispositif afin que, pendant un nombre prédéterminé de premières périodes successives, un champ variant avec la position X soit créé et irradie l'objet par des premiers signaux 5 ayant une première amplitude et une première durée prédéterminées, et ayant des spectres prédéterminés de fréquences correspondant à des plages particulières de fréquences de Larmor si bien que des spins nucléaires, dans les plans respectifs Y-z de l'objet, se réorientent de 90° environ, ce sixième dispositif 10 étant aussi destiné à créer un second jeu de signaux de commande destiné au second, au troisième et au quatrième dispositif et destiné à provoquer, pendant un nombre prédéterminé de secondes périodes successives suivant les premières périodes, la création d'un champ magnétique variant avec la position Y, et 1'irradiais tion de l'objet par des seconds signaux ayant une seconde amplitude et une seconde durée prédéterminées et ayant des spectres de fréquences correspondant à des plages particulières de fréquences de Larmor si bien que les spins nucléaires compris dans des plans X-Z respectifs se réorientent de 180° environ, les plans 20 respectifs X-Z recoupant les plans respectifs Y-Z le long de lignes uniques respectives d'intersection dans l'objet, l'intersection provoquant ultérieurement la création de signaux d'écho de spin, ce sixième dispositif étant aussi destiné à créer un ♦ * troisième jeu de signaux de commande destiné au second dispositif * 25 et à créer des gradients Z respectifs à des moments dépendant des intervalles de temps entre les irradiations de l'objet par les premiers et les seconds signaux si bien que des signaux d'écho de spin sont créés en présence des gradients Z, un septième dispositif destiné à recevoir les signaux 30 d'écho de spin et à créer des signaux représentatifs correspondants , un démodulateur en quadrature, commandé par les signaux représentatifs des signaux d'écho de spin et des signaux respectifs en phase et en quadrature représentatifs du signal de por-35 teuse, et destiné à créer des signaux respectifs de sortie représentatifs des différences de phases et de fréquences entre les signaux d'écho de spin et le signal en phase, et les signaux d'écho de spin et le signal en quadrature, 46 un dispositif de calcul des sommes vectorielles des signaux de sortie du démodulateur en quadrature, un dispositif de calcul des transformées de Fourier des sommes vectorielles, destiné à former des valeurs corres-5 pondant aux composantes de fréquences, et un huitième dispositif destiné à transmettre des signaux de commande au second dispositif afin que les champs magnétiques variant avec la position soient appliqués à l'objet pendant diverses périodes pëdêterminées, de manière que l'étale- 10 ment parasite des phases des signaux de spin soit corrigé. -A $ *>