FR 2574551 A1 19860613 FR 8418997 A 19841212 PROCEDE DE GENERATION ET DE TRAITEMENT DE SIGNAUX POUR L'OBTENTION PAR RESONANCE MAGNETIQUE NUCLEAIRE D'UNE IMAGE EXEMPTE DE DISTORSIONS A PARTIR D'UN CHAMP DE POLARISATION INHOMOGENE La présente invention se rapporte aux techniques d'imagerie par la méthode dite de résonance magnétique nucléaire (RMN) actuellement en voie de développement rapide, notamment pour l'examen du corps humain, dont elle permet d'obtenir des images d'une qualité et d'une précision inconnues jusqu'ici avec les méthodes de la radiologie classique aux rayons X. Elle s'applique à la phase de construction de l'image obtenue par codage spatial et calcul à l'aide d'un ordinateur. La méthode d'imagerie RMN utilise la propriété de résonance magnétique nucléaire de certains noyaux présents dans le corps humain, essentiellement les pro- tons répartis dans tout l'organisme, et quelques noyaux d'intérêt biologique mais en abondance beaucoup moins grande tels que le phosphore 31p, le potassium 39K le sodium 23Na. Elle permet, en dressant la carte de leur concentration en chaque point du volume examiné, de fabriquer des images des tissus vivants, notamment à partir de l'hydrogène contenu dans l'eau et les lipides qui sont deux composants essentiels de toute matière vivante. Une installation de RMN médicale se compose principalement d'un aimant destiné à produire dans toute la zone du corps à examiner un champ magnétique de pola- risation Bo, statique et uniforme, auquel on superpose, à l'aide d'une bobine auxiliaire, des impulsions de champ de radiofréquence dans un plan perpendiculaire à la direction du champ o précédent. Les images sont obtenues en faisant résonner le plus souvent des noyaux d'hydrogène ou protons, contenus dans les tissus biologiques. Cette résonance est possible parce que chaque proton se comporte comme un microscopique aimant, et que le champ de radiofré- quence est équivalent à deux champs tournant en sens inverse dont celui qui tourne dans le sens de précession des spins est susceptible de se coupler avec eux. Lorsque l'on applique le champ statique Bo, les spins s'orientent tous parallèlement à l'axe du champ. Ils n'ont donc plus que deux orientations possibles: dans le sens du champ ou en sens contraire. Lorsqu'on applique en plus le champ de radiofréquence à une fréquence quelconque, il ne se passe rien. Mais si cette fréquence est choisie égale ou très voisine de fo (fréquence de Larmor), telle que 2wfo=yBO, formule dans laquelle y est une constante physique caractéristique des noyaux que l'on veut mettre en résonance et dénommée rapport gyromagnétique, le-couplage des spins et du champ tournant atteint une valeur telle que ceux-ci entrent en résonance. Le signal de résonance détecté au cours de leur retour à l'équilibre en précession libre est proportionnel à l'aimantation t des noyaux placés dans le champ magnétique Bo de polarisation. L'observation du phénomène de résonance- requiert la présence d'un champ magnétique élevé (de 0,1 à l tesla) et-de très grande uniformité dans tout le volume à imager. On rappellera maintenant un certain nombre de phénomènes connus et l'on fixera un certain nombre de conventions qui aiderontà la compréhension de la suite du présent texte. Tout d'abord on repérera le volume à imager par rapport à un trièdre trirectangle O,x,y,z dans lequel l'axe Oz sera, par convention, parallèle au champ continu de polarisation Bo. Dans cette hypothèse, le vecteur aimantation macroscopique M du milieu à imager sera, au repos, parallèle à cet axe Oz et le champ de radiofréquence B1, tournant à la vitesse angulaire wo sera situé dans le plan x,o,y qui est également celui dans lequel on recueille le signal de résonance. Pour décrire ces divers phénomènes, on a souvent recours à un référentiel x', O, y', z tournant autour de Oz à la pulsation wo (Fig. 1 et 2). Dans ce référentiel, B1 est aligné sur Ox'. Une impulsion de radiofréquence à la fréquence de Larmor fo dirigée selon l'axe Ox' a la propriété de faire tourner le vecteur aimantation macroscopique d'un angle e autour de Ox dans le plan y'Oz.tel que e = yBlt, formule dans laquelle y est le rapport gyro- magnétique, B1 le module du champ de radiofréquence et t la durée d'application de ce champ. Dans les techniques de RMN, on utilise surtout les impulsions dans lesquelles e = 900 et e = 1800 dites respectivement impulsion à 900 et. impulsion à . Ces deux types d'impulsion ont respectivement les propriétés suivantes. Les impulsions de 900 selon Ox' font basculer le vecteur M selon l'axe Oy' comme on le voit sur la Fig. 1. A partir de cette position initiale, le vecteur M se met à précesser dans le plan xoy, également à la pul- sation wo = 2wfo comme B1 et c'est ce mouvement qui induit précisément le signal de résonance à la fréquence de Larmor que l'on détecte à l'aide d'une bobine non représentée placée dans ce même plan xOy. Au moment de l'application selon Ox' d'une impulsion à 900, tous les moments magnétiques des diffé- rents noyaux sont mis en phase selon Oy' et le vecteur M a son module maxiinmal-M selon Oy'. Au fur et à mesure de l'écoulement du temps, une dispersion des moments ma- gnétiques selon par exemple deux cposantes M et M2 tournant en sens inverse dans le plan x'Oy' par rapport à Oz se produit sous la double influence de l'inhomogénéité du champ B = B0 + E (x,y,z) et des interactions entre noyaux. C'est le phénomène de la relaxation transverse qui conduit le moment résultant M à décroître avec une constante de temps T2 de façon telle que IMI = Me-t/T2 et le signal de résonance à s'éteindre. Les impulsions de 180' ont la double propriété de réaliser une inversion de population des spins en plaçant le vecteur aimantation macroscopique Mo dans sa position antiparallèle et d'effectuer une transfor- mation symétrique par rapport à la direction de cette impulsion, à la manière de l'image dans un miroir plan, des vecteurs aimantation existant au moment de l'impulsion. Ceci équivaut par conséquent, dans le cas o l'impulsion de radiofréquence à 1800 est dirigée selon Oy', à changer en -.0 les phases 0 des vecteurs M1 et M2 qui s'échangent alors mutuellement,-et à changer 0 en (i - 0) dans le cas o cette impulsion à 180 serait dirigée comme B1 selon Ox', M1 et 2 prenant alors les places de 3.etM4. Les différentes considérations précédentes ayant été exposées, on expliquera maintenant les techniques connues pour obtenir par le calcul une image de RMN d'un organisme vivant. Si le champ continu de polarisation B est parfaitement homogène, tous les vecteurs aimantation issus des différents points du volume à imager sont en phase et il est impossible, dans le signal global reçu: s(t) (x,y,z) e dx dydz xyz d'attribuer à chaque point une fréquence et une intensité particulières. Pour parvenir à identifier la contri- bution à ce signal global de chaque volume élémentaire, on réalise un codage en fréquence de l'espace à imager, en superposant au champ Bo des champs additionnels variant linéairement suivant l'une des coordonnées, de sorte que la composante Bz a des gradients spatiaux Bz, y z B z constants selon x, y et z, satisfaisant les trois égalités suivantes: VBz 0 B z(x) = x x x Bz(y) = ZBz y. -y = Gy.y Bz(z) = DBz z = Gz.z dans lesquelles Bo+Bz est la composante du champ continu de polarisation suivant Oz et Gx, Gy, Gz les trois gradients constants de l'induction résultante selon Ox, Oy et Oz. L'usage a prévalu en imagerie par RMN d'appeler, de façon impropre, ces champs additionnels "gradients linéaires" ("linear gradients") et nous nous conformerons à cet usage dans la suite du présent texte. On comprend facilement le but recherché à l'aide de ces gradients linéaires: la fréquence ou pulsation de Larmor w = yB étant, en tout point du volume à imager, proportionnelle à l'induction en ce point, on réalise ainsi un codage spatial de ce volume en phase et en fréquence et on peut montrer que la transformée de Fourier du. signal temporel s(t) ou spectre de fréquences s(w) de ce même signal est l'image tridimensionnelle elle-même de l'échantillon en densité d'aimantation puisqu'en presence d'un gradient linéaire la fréquence représente la coordonnée selon laquelle est appliqué ce gradient. A titre illustratif, la Fig. 3 montre la répartition des vecteurs Bz dans le cas d'un gradient linéaire Gx selon la direction Ox. C'est la méthode de construction de l'image dite méthode de Fourier directe; elle s'applique aussi bien d'ailleurs à l'imagerie tridimensionnelle (en anglais "3D, FT") que bidimensionnelle (en anglais "2D, FT") et on en rappel- lera maintenant de façon succincte le principe, car c'est à cette méthode d'imagerie que s'applique spécia- lement bien l'invention. Méthode de Fourier directe & trois dimensions (3D, FT) On a représenté sur la Fig. 4 l'évolution dans le temps représenté en abscisses des différents champs appliqués au volume à imager et du signal de résonance magnétique des protons reçu. Le graphe 1 montre l'impulsion de radiofréquence à 904 qui déclenche une séquence d'analyse du volume à imager en déclenchant le processus de relaxation trans- verse de l'aimantationM et son mouvement de précession libre à la fréquence de Larmor autour de la direction Oz parallèle au champ continu Bo de polarisation. Pendant un temps t suivant cette impulsion, on applique comme représenté sur les graphes 2 et 3 les gradients linéaires Gz et Gy qui réalisent un codage en phase dans ces deux directions en chaque point du volume, les déphasages ainsi acquis par l'aimantation M(x, y,z) pouvant s'écrire iyGz.z.T iYGy.y.t en notation complexe e et. Pratiquement d'ailleurs, les différents gradients Gz et Gy sont appliqués par valeurs croissantes de -G à +G et, pour chaque valeur d'un gradient Gz on applique successivement tous les gradients Gy de -Gy à +Gy, ce qui donne si N et N' sont les nombres respectifs de gradients Gz et G appliqués, un nombre de mesures égal à N x N'. A la fin de chaque codage de phase en z et en y, on applique pendant un temps TX le gradient Gx ou gradient de lecture, comme représenté. sur le graphe 4 et l'on recueille alors le signal de résonance repré- senté sur le graphe 5. Ce signal, en supposant le champ de polarisation homogène, est de la forme: s(t) = fffm(x,y,z)e t/T2eiYGYyT eiTGz*z*T xyz. eiG .x.t.dx.dy.dz formule dans laquelle les deux premières exponentielles représentent les déphasages acquis par les vecteurs aimantations M pendant le codage de phase et la troisième le déphasage acquis pendant l'acquisition du signal. C'est sur ce signal s(t) que l'on effectue ensuite, dans un ordre quelconque les trois transformées de Fourier par rapport à Gz à G et au temps t pour y obtenir, de façon connue,à l'aide d'un calculateur, la répartition spatiale de l'aimantation des différents points du volume à imager sous forme d'une image complexe tridimensionnelle, comportant un volume "réel" et un volume "imaginaire". Méthode de Fourier directe à deux dimensions (2D, FT): On a représenté sur la Fig. 5 l'évolution dans le temps, représenté en abscisses, des mêmes grandeurs que sur la Fig. 4, à l'aide des mêmes graphes 1, 2, 3, 4, 5. Dans cette méthode toutefois, on utilise le premier gradient Gz pour définir, dans le volume à imager, une coupe ou tranche, perpendiculaire à Oz, d'épaisseur Az = 2w/yGzTz centrée sur la cote zo = (2wfo - yB)Gz et dans laquelle on travaille ensuite selon les deux coordonnées x et y. Cette coupe sélective est obtenue en associant pendant le temps TZ (Fig. 5) une impulsion de radiofréquence à 90 , RF, à la fréquence moyenne fo et le gradient linéaire Gz = -z comme le montrent les graphes 1 et 2. Ensuite interviennent, comme dans la méthode tridimensionnelle, une période de codage de phase avec le gradient Gy de durée Ty et une période de codage en fréquences avec le gradient Gx de durée Tx, pendant laquelleon effectue également l'acquisition du signal s(t) =J M(xy)e-t/T2 eiYGy.y'Ty iyGx.x.t dx.dy Xy o les deux exponentielles représentent les déphasages acquis pendant le codage de phase et le codage de fréquence par les différents vecteurs aimantations M de la tranche examinée. Les deux transformations de Fourier par rapport à Gy et par rapport au temps t interviennent ensuite et permettent d'obtenir un signal analytique complexe comportant une image "réelle" et une image "imaginaire" de la répartition spatiale de l'aimantation dans la tranche sélectionnée, et, par là-même, une image tomo- graphique de cette même tranche. La théorie complète des principes de la résonance magnétique nucléaire et de son application à l'imagerie médicale sont notamment décrites dans les ouvrages suivants: Shaw, D. Fourier Transform NMR Spectroscopy, Elsevier, Amsterdam, 1976.; Farrar, T.C. and Becker, E.C. Pulse and Fourier Transform NMR, Academic Press, 1971.; Abragam, A. The Principles of Nuclear Magnetism, Oxford University Press (Clarendon), London, 1961; Mansfield, P. and Morris, P.G. NMR Imaging in Biomedicine, Academic Press, 1982; Kaufman, L., Crooks, L.E. and Margulis, A.R. Nuclear Magnetic Resonance Imaging in Medicine Igaku-Shoin, NewYork, Tokyo, 1981, Chapters 1 and 2. La-théorie connue de l'imagerie RMN qui vient d'être rappelée brièvement Suppose toutefois une homogénéité-parfaite du champ de polarisation Bo appliqué, ce qui, dans la pratique, n'est jamais réalisé complé- tement et conduit à des'défauts graves et coûteux pour les installations de RMN réalisées à ce jour. En effet, si l'on représente le champ de pola- risation par la fonction B(x,y,z) = Bo + E(x,y,z) dans laquelle E(x,y,z) représente l'inhomogénéité du champ au point (x,y,z) par rapport à Bo souhaité, on doit compléter l'écriture précédente du signal s(t) reçu par deux exponentielles en facteur avec les précédentes et qui s'écrivent: eiyE(x'y'z)T. e iyE(x,y,z)t pour tenir compte des phases acquises par les différents vecteurs aimantation, dues à l'inhomogénéité du champ, durant les périodes de codage de phase et de codage de fréquence respectivement. La première de ces deux exponentielles est une constante qui se conserve en l'état après le codage de phase, mais la seconde est, malheureusement à la fois une fonction de x,y,z et du temps t. Elle conduit,quand on utilise les transformées de Fourier,à des distorsions stigmatiques aggravées d'un brouillage de l'information des images "réelle" et "imaginaire" qui mélange les termes d'absorptionet de dispersion et rend ces défauts non corrigibles si l'on ne connait pas, à l'avance, la fonction E(x,y,z) et la carte du champ. Ces défauts se traduisent, lors de la reconstruction de l'image par ordinateur, par une affectation spatiale erronée de l'intensité d'aimantation de chaque volume élémentaire constituant l'échantillon. Ces distorsions stigmatiques conduisent à des déforma- tions importantes de l'image et sont surtout gênantes lorsque l'on fait de la microimagerie et que l'on veut représenter des détails très fins. Par ailleurs, l'inhomogéneité du champ en RMN affecte également la sensibilité des mesures (rapport signal/bruit). En effet, si l'inhomogênéité de B se traduit sur le-volume de l'échantillon par une dispersion AB, la limite de résolution spatiale est approximativement AB/G en présence d'un gradient G. La limite de résolution étant fixée, la dimension de l'échantillon étant donnée, l'intervalle spectral sur lequel s'étalent les signaux de résonance est d'autant plus élevé que AB est importante. L'énergie du signal par intervalle spectral en est d'autant plus faible et le travail en champ inhomogène est mauvais pour la sensibilité. On peut évidemment chercher à réaliser des aimants perfectionnés dont le champ soit homogène à un haut degré de précision près. Toutefois, l'obtention de AB/B de l'ordre de 10-6 ou mieux sur des volumes avoisinant 10 litres, nécessite un excellent alignement des éléments et/ou bobines qui créent le champ magnétique. Le prix d'aimants de haute résolution est, pour ces raisons techniques, très élevé (ex:un cryoaimant donnant un champ de 2,3 Teslas avec une inhomogénéité de 10-6 dans une sphère de 30 cm de diamètre conte approximativement 1,5 MF). Or l'imagerie du proton (eau b 80 % du corps) ne pose pas de problème quant à la sensibilité, car les signaux protoniques sont très intenses et nombreux. Ainsi, une méthode d'imagerie -4 en champ inhomogène (AB/B 6 10), qui s'affranchirait des distorsions dues à la non uniformité du champ ma- gnétique de polarisation, peut s'envisager dans le cas du proton (imageurs médicaux essentiellement), sans grande perte au niveau de la sensibilité. I1 en résulterait une substantielle économie sur le coQt -des aimants, donc le coût général de l'appareillage. On signalera enfin pour mémoire, les méthodes qui consistent à dresser au préalable la carte E(x,y,z) du champ et à entrer les valeurs mesurées dans le calculateur qui construit l'image. Il est bien évident que ces méthodes sont peu satisfaisantes, car longues et délicates d'emploi, et par là-même, très coQteuses. 1! La présente invention a précisément pour objet un procédé de génération et de traitement de signaux pour l'obtention par RMN d'une image exempte.de distorsions à partir d'un champ de polarisation inhomogène qui permet, à l'aide d'un simple doublement du nombre de mesures effectuées, de s'affranchir totalement des inhomogénéités du champ de polarisation, y compris de celles introduites par la présence du corps à imageret qui fournit, de surcroît, comme sous-produit, la carto- graphie du champ de l'aimant de polarisation. Ce procédé, qui utilise les méthodes connues de codage de l'image en deux ou trois dimensions (en terminologie anglosaxonne "2DFT" et "3DFT") comportant l'application d'une impulsion de radiofréquence à 90 dans le plan xOy, puis un codage de phase pendant un temps T suivi d'un codage de fréquence par application selon Ox d'un gradient de lecture pendant l'acquisition du signal de précession libre, signal auquel on fait subir ensuite des transformations de Fourier par rapport au temps et par rapport au(x.gradient(s)..appliqué(s) et que l'on recueille sous forme complexe, se caractérise pr les étapes suivantes: a) au cours d'une première série de séquences d'acquisition des signaux comportant un codage de phase et un codage de fréquence, on recueille de la façon. connue précédente les n signaux complexes si (T, t); b) au cours d'une deuxième série de séquences comportant l'application successive des mêmes gradients, on soumet le système en précession libre à une seconde impulsion de radiofréquence de 180 au début de l'éta- blissement de chaque gradient de lecture, de façon à changer en - 0 {ou X - 0 selon le cas) les phases 0 acquises par les différents moments magnétiques durant le codage de phase et l'on recueille les n signaux imaginaires s2 (T, t), les étapes a) et b) précédentes pouvant être successives ou imbriquées dans le temps; c) on forme les signaux s (T,t) conjugués * 2 de s2 (T,t) qui sont tels que s2 (T,t) = s1 (T, -t), et on les associe un par un avec leur homologue s1(T,t) de façon à former n signaux complexes s(T,t) dans lesquels t varie de -T à +T, T étant le temps d'acqui- sition de chaque signal compté à partir du début de l'application du gradient de lecture; d) on effectue sur ces signaux s(T,t) pour lesquels T appliqué(s), ce qui fournit une image complexe multi- dimensionnelle dont les parties réelle R et imaginaire I ne comportent que des termes d'absorption pure; e) on effectue sur chaque point x,y,z, les corrections de distorsion selon Ox et d'affectation d'intensité d'aimantation en calculant successivement la valeur correcte de la phase 0 = Arc cos - la fonction E(x,y,z) représentant la cartographie du champ par la relation 0(x,y,z) = yE(x,y,z)T, puis par la relation Gx.x + E(x,y,z)-Gx. x' = 0, la valeur x' de x à laquelle doit être affectée l'intensité d'aiman- tation trouvée. Le fait de doubler le nombre des mesures effectuées par rapport aux méthodes connues jusqu'à ce jour en faisant précéder une fois sur deux l'éta- blissement du gradient de lecture Gx d'une impulsion de 180 permet donc de définir un signal de résonance qui évolue dans le temps de -T à +T et d'obtenir, après les deux transformées de Fourier, une image, -certes distordue dans la direction du gradient de lecture, mais exempte de brouillage, c'est-à-dire dont les parties réelle R et imaginaire I selon la direction de ce même gradient de lecture ne comportent que des termes d'absorption pure: ceci permet de corriger ensuite la distorsion, par calcul en chaque point de la phase exacte, de l'inhomogénéité E(x,y,z) du champ, puis enfin de la valeur exacte x' de l'ordonnée x à laquelle doit être attribuée l'intensité d'aimanta- tion trouvée au point (x,y,z). Dans les explications qui vont suivre, on raisonnera, pour plus de simplicité dans les exposés mathématiques,en imagerie bidimensionnelle (x, y) avec sélection de coupe selon Oz. Il doit être bien entendu toutefois que ceci n'est nullement limitatif du domaine d'application de l'invention et que la généralisation au cas d'une imagerie tridimensionnelle est tout à fait à la portée de l'homme de métier. Par ailleurs, le choix particulier qui est fait dans tout le présent texte des rôles particuliers des trois coordonnées x,y,z vis-à-vis du champ de polarisation et des directions des gradients appliqués résulte de la pratique internationale habi- tuelle en ce domaine, mais ne saurait nullement cons- tituer une limitation de portée du procédé objet de la présente invention. On donnera maintenant l'explication mathé- matique de l'efficacité du procédé selon l'invention, en se référant aux Fig. 6 à 16 suivantes sur lesquelles: - les Fig. 6 et 7 illustrent les deux séquences nécessaires, dans le cadre de l'invention, à la construc- tion d'une image RMN par le procédé 2DFT avec sélection de tranche perpendiculaire à Oz; - les Fig. 6a et 7a illustrent les deux séquences nécessaires dans le cadre de l'invention, à la construction d'une image RMN par le procédé 3DFT; - les Fig. 8, 9 et 10 illustrent les distorsions de forme et d'intensité d'une image & deux dimensions dues à l'inhomogénéité du champ de polarisation Bo que l'invention permet de corriger; - les Fig. 11 à 16 représentent les termes d'absorption et de dispersion du spectre complexe de l'image obtenue par transformée de Fourier dans différentes hypothèses envisagées. Pour tenir compte de l'inhomogénéité du champ de polarisation, on écrira la valeur de ce champ sous la forme B(x,y) = Bo + E(x,y), E(x,y) représentant le défaut d'homogénéité ou carte du champ au point (x,y). Dans cette hypothèse, au cours de la première série de séquences d'acquisition de signaux sans impulsion de 180 (Fig. 6), on recueille les signaux: sl(try) = fjM(xY) eiYGY'YTY.eiYE(X'Y)TY xy eiyGx.x.t eiyE(x, y)tdxdy formule dans laquelle on néglige par commodité l'amortis- sement e t/T2 et la pulsation eiw t du signal, ce qui est légitime pour ce dernier élément car on effectue une détection synchrone. Dans cette formule, les deux premières expo- nentielles représentent les phases acquises par chaque vecteur aimantation M (x,y) pendant le codage de phase et dues respectivement à l'application du gradient Gy et au défaut d'homogénéité du champ E(x,y). Pour une valeur de y donnée, ces phases sont des constantes. Les deux dernières exponentielles représentent les phases, acquises par ces mêmes vecteurs M (x,y) pendant le temps t de lecture, et dues respectivement à l'application du gradient Gx et au défaut d'homogé- néité du champ E(x,y). Ces phases dépendent du temps t et, pour la dernière, iyE(x,y)t, aàla fois de x et de y. C'est précisément ce terme qui, à cause de sa double dépendance de x et de y serait très gênant et "obscurcirait" l'image de façon irrécupérable si on en restait là. Au cours de la deuxième série de séquences d'acquisition de signaux, représentée Fig. 7, on fait précéder l'application de chaque gradient de lecture Gx d'une deuxième impulsion de radiofréquence de 180 à la fréquence wo, impulsion que l'on supposera, dans cette partie de la démonstration, en quadrature avec la première, c'est-à-dire selon l'axe Oy' de la Fig. 2. Il en résulte donc, conformément aux explications fournies dans l'introduction du présent texte, un changement de signe des phases 0 acquises par chaque vecteur aimantation M(x,y) pendant le codage de phase et les signaux s2(Ty,t) recueillis ont la forme générale: sx (te y) i8GY Y TY i*(X.Y)Iy iyx.x.t iyE(xy)t xS 2(trY)y e e e dxdy Si l'on forme maintenant les signaux imagi- naires conjugués s* (t,T), ils s'écrivent sous la 2 y forme générale: s. (ty) =M(xy)eiyGy-y.Ty.eiyE(xy)TY -iyGxx.t e-iyE(x,y)tdxdy 2JJY)=fMxy)e.e. e x etdxdy xy et l'on constate facilement l'égalité: s2 (t, Ty) =Sl(-t,Ty) C'est là la caractéristique essentielle du procédé selon l'invention qui permet, en associant deux par deux les 2n signaux globaux reçus de 0 à +T de définir n signaux s (t, Ty) dans lesquels t varie de -T à + T. Il suffit, pour ce faire, d'associer chaque signal s1(t, Ty) de la première séquence avec *y son homologue s2(t, Ty) de la seconde séquence obtenu après application d'une impulsion à 180 et conjugaison du signal. C'est cette possibilité d'intégrer les signaux de -à à + au cours des transformées de Fourier ou, plus exactement de -T à T, car l'infini n'existe pas pour le physicien, qui permet, bien que les signaux soient réellement lus de 0 à T, de supprimer le brouillage de l'image dont les composantes, grâce au procédé objet de l'invention, ne comportent plus que des termes d'absorption pure. Ce résultat découle d'un théorème de mathématiques connu, d'après lequel: r+o J e-ix dx = A(A)+iB(X), X paramètre quelconque, alors que J eix dx = A(X), A et B fonction de X, c'est-à-dire que l'intégration de - à ± permet de supprimer la partie imaginaire du résultat de l'inté- gration. Cette propriété sera démontrée plus loin dans un cas particulier. Les signaux sl(t, Ty) précédents, peuvent s'écrire également sl(Gy, t) puisqu'en fait, c'est, lors de chaque mesure Gy qui varie, pendant le même temps Ty de -G à +G. Dans ces conditions, la double transformée de Fourier par rapport à Gy et à t s'écrit: A(yo6xx) (= iyGy.yo.T -iyGx.x't A(yoGX.X')(Sy,t)e-.e dGy.dt La première transformée de Fourier par rapport à Gy est égale à: 1 (T[ Mx, o)iyE(xyO)t eiy{Gx.x+E(x,yo)-Gx.x'} tdx lITJM (xyo)e e.dX formule dans laquelle x' est l'abscisse corrigée de la distorsion à laquelle il faut attribuer l'intensité trouvée pour le point x,y. La deuxième transformée de Fourier par rapport à t est égale à: A(yo,Gx. x') = M(xoyo)eiYE(XOyO)T (1) formule qui donne finalement les deux tableaux de points (x,y) affectés de leur intensité M(x,y) constituant cequ'il est convenu d'appeler les parties "réelle" et "imaginaire" de l'image complexe et dans laquelle xo et i sont donnés par les équations: Gx.xo + E(xo,yo) - Gx.x' 0 (2) et - -Jà7x- Gx.x + E(xyot - Gx.x'] x = xo + 3) Les couples de valeurs M, x et y donnés par l'équation (1) sont entachés d'une erreur systématique due à l'inhomogénéité E(x,y) qui se traduit par une distorsion de l'image selon la direction x du gradient de lecture. En d'autres termes, il est encore nécessaire pour corriger cette distorsion de l'image de chercher, pour une valeur donnée de y, la valeur x' exacte de x à laquelle on doit attribuer l'intensité d'aimantation M. Ceci est réalisé de la façon suivante. La comparaison des parties réelle R et imaginaire I de l'équation (1) permet de calculer la phase 0 du signal au point x,y par exemple par la formule: 0 = Arc cos R et la valeur de E(xo,yo). Cette valeur, reportée dans (2), permet de déduire xo car Gx et Gx.x' sont connus. Enfin, la connaissance de E(xo,yo) pour toute valeur de x permet de calculer $. On parvient ainsi finalement à la détermination de M(x,y) recherché et on obtient E(x,y) ou cartographie du champ à titre de sous-produit de la méthode. On voit donc que c'est la possibilité d'effec- tuer les transformées de Fourier sur un signal qui varie de -T à +T qui permet, selon le procédé objet de l'invention, de construire une image corrigée de ses distorsions par calcul de la phase exacte en chaque point et de l'abscisse exacte x' pour chaque intensité M(x,y) mesurée. Par ailleurs, le fait que le procédé conduise à l'obtention, à titre de sousproduit, de la cartographie exacte du champ E(x,y) -est aussi l'un des avantages majeurs de l'invention. L'extension de ce qui précède à une imagerie à trois dimensions ne pose aucun problème de principe. En effet, les séquences nécessaires à une telle imagerie RMNNdu type "3DFT" appliquant le procédé objet de l'invention sont représentées schématiquement sur les Fig. 6a (impulsion à 90 seule) et 7a (impulsions à 900 et à 180 ) qui se lisent d'elles-mêmes par compa- raison avec les Fig. 6 et 7. On voit simplement que les graphes 2 et 3 correspondent à des codages de phase selon Oy et Oz, les graphes 4 au codage de fréquence suivant Ox, les impulsions de radiofréquence (graphes 1 et 3') et la lecture (graphes 5) se déroulant comme dans l'exemple des Fig. 6 et 7. Conformément à l'invention, les signaux obtenus à l'aide des deux séquences des Fig. 6a et 7a sont associés comme expliqué à propos de l'exemple des Fig. 6 et 7, de façon à définir un signal dans les intervalles [-T, +T], [-Gz +Gz] et [-G, +Gy] L'analyse mathématique du procédé mis en oeuvre est alors la suivante: soit A(x,y,z) = yE(x,y,z) l'écart de pulsation au point (x,y,z) par rapport à la valeur nominale de la pulsation; A(x,y,z) mesure donc l'inhomogénéité en (x,y, z). Le signal recueilli est de la forme: PPi s(Gyz,t = Gdx dy dz M(x,y,z) ei xyz)T xyz eiyGyyT eiyGzzr eiyGxxt eiA(x,y,z)t Les trois transformations de Fourier classiques sont effectuées par rapport à G G et t dans un ordre indifférent. Pour un signal S, ces transformations de Fourier sont: - par rapport à G i t. S e-iG dG y 2w - par rapport à Gz: 1f+GZ -iyGzz , 1 J+Gz S e dGz 2 J-G- z - par rapport à t: +T -ikt 1 f -S e at i J-T d'o le résultat i -M du champ. Comme dans le cas de l'imagerie bidimen- sionnelle, on corrige ces erreurs de distorsion de l'image en calculant successivement pour chaque point: - la phase correcte 0 = Arc cos R, R et I étant respectivement les parties réee +e et imaginaire de l'image complexe selon la direction x; - la fonction E(x,y,z) par la relation: 0(x,y,z) = yE(x,y,z)T; - la valeur x' de x à laquelle doit -tre affectée l'intensité M(x,y,z) trouvée pour le point considéré, par la relation: Gx.x + E(x,y,z) - Gx.x' = 0. Il est utile, au stade actuel de l'exposé, de faire les remarques suivantes. D'abord, le calcul précédent suppose que l'équation générale Gx.x + E(x,y, z) - Gx.x' = 0 n'ait qu'une seule solution. Ceci est heureusement réalisé sans difficulté si Gx est suffisamment grand pour que Gx.x + E(x,y,z)soit une fonction monotone, ce qui est en pratique toujours réalisé. Ensuite, les calculs précédents sont effectués dans la réalité sur des échantillonnages discrets des fonctions étudiées, puisque tous les signaux sont digitalisés pour traitement par l'ordinateur et le signal analytique complexe recueilli, qui résulte comme toujours en RMN d'une détection en quadrature (partie "réelle" ou "en cosinus" et partie "imaginaire" ou "en sinus") se compose en fait d'une double liste d'ensembles de points discrets. Ceci n'est pas un obstacle à l'exécution pratique des calculs formels exposés. Enfin, selon une variante du procédé objet de l'invention, la seconde impulsion de radiofréquence de 1800 est en phase avec la première appliquée au début du codage de phase et l'on obtient des signaux S2 (Tt) = -Sl(T,-t), que l'on change de signe avant de les associer à leur homologue de la première séquence. En effet, si la seconde impulsion de radio- fréquence de 1800 est appliquée, selon Ox' en phase avec la première de 900, les phases acquises pendant le codage de phase se changent instantanément en (r - 0) comme expliqué sur la Fig. 2. Il en résulte que les exponentielles de l'expression s,(t,T) concernant la phase 0 acquise pendant le codage de phase se trans- formnt de eo. eni (ir-O) = - forment de ei9 enei( -) = -e -. Il en résulte que, dans ce cas, s* (T,t) = -Sl(T,-t). Il est donc nécessaire, dans ce cas, de changer de signe les signaux s2 (T,t) avant de les associer a leur homologue de la première séquence pour calculer les intégrales de Fourier de -T à +T. On donnera maintenant, -à titre de précisions complémentaires sur le procédé objet de l'invention, un certain nombre d'explications et-de justifications sur les raisons d'être et le mode de fonctionnement de ce procédé qui permettront, dans le cas parti- culier d'une image bidimensionnelle, de mieux comprendre les données théoriques qui ont été exposées jusqu'ici.- Pour l'image complexe calculée par l'ordina- teur, l'inhomogénéité du champ Bo de polarisation conduit, aussi bien pour la partie réelle que pour la partie imaginaire à deux types de distorsions: a) distorsion de forme (ex: une croix droite (Fig. 8) voit sa branche verticale devenir inclinée (Fig. 9), car les distorsions de forme se manifestent selon les directions Oy perpendiculaires à celles Ox du gradient de lecture Gx qui est horizontal); b) distorsion d'intensité (ex: une région dont le niveau de gris est normalement uniforme, est alors striée de bandes alternativement plus intenses ou plus sombres, selon des directions perpendiculaires aux équipotentielles de champ magnétique, qui sont ici horizontales (Fig. 10) (cas d'un gradient résiduel selon Oy pour E(x,y)). L'image imaginaire présente des distorsions de type a) et b) analogues. Toutefois, la grande différence entre le procédé selon l'invention et l'art antérieur réside dans le fait que ces distorsions s'effectuent sans brouillage de l'information; ce sont des distorsions sans astigmatisme alors qu'en procédure usuelle elles sont astigmatiques, donc de véritables abréviations qu'il n'est pas possible de corriger (sauf si on connait au préalable la cartographie du champ et que l'on en tient compte dans le processus de construction de l'image). L'exposé qui suit permet de comprendre pourquoi. Pour le comprendre, il est nécessaire de se rapporter aux propriétés particulières de la transformation de Fourier, transformations opérées sur des signaux dont l'amplitude décroît avec la valeur du paramètre d'acquisition. En imagerie par 2DFT, il y a deux paramètres d'acquisition: -2574551 - le temps (s o si ItI > Tm, Tm de l'ordre de T, en raison de la relaxation transversale, - le gradient Gy; si Gy est élevé, le déphasage des spins nucléaires est très rapide et le signal recueilli petit (s o si IGYl > Gm, Gm de l'ordre de G valeur maximale absolue du gradient de codage utilisé, -G.. Gy Gy occupe toujours un domaine comportant des valeurs positives et des valeurs négatives, en général opposées entre elles et qui correspondent toutes à ce qu'on appelle l'échantillonnage de la valeur du gradient. Pour simplifier l'analyse, on considérera la dépendance d'un signal par rapport au temps (par rapport à Gy le problème est absolument le même, donc un seul raisonnement suffit). Pour simplifier encore, considérons un signal "monochromatique", ce qui n'enlève aucune généralité au raisonnement puisque le signal reçu peut être considéré comme la super- position (l'addition) de signaux monochromatiques. Le raisonnement portera donc sur un signal: s(t) = a e-t/T2 e jeot 1 T2 caractérise l'amortissement transverse du signal (en RMN on l'appelle temps de relaxation spin-spin apparent, et il gouverne la valeur de la largeur de raie). La transformée de Fourier de s(t) (ou spectre complexe) est définie comme suit: + S(in = s(t) e jwt dt Si le signal est défini de 0 à (s(t) 0 si t Le spectre est complexe. La partie réelle Ro = aT2 1+(Wo-W) 2T22 représente le terme d'absorption (notation usuelle en spectroscopie) représenté sur la Fig. 11; La partie imaginaire représente le terme de dispersion: Io = aT22(mo-m) 1+ (wo-w)ZTZ représenté sur la Fig. 12. Si le signal est défini de -- +, on peut écrire:+ S(1T+ (wo-)t S(W) = a e (-1/T-+j(wo-w)t dt, mais ceci est faux, car il faut tenir compte du fait que si t ±- s(t) ±> 0. Ceci conduit aux égalités suivantes S(m)= a C ) = a f f - iti/T2 eJ (WO-W tdt e(-1/T2+j ( wo-w)tdt {-e o S( = a{j-oeItI/T2ej(wo-w)t dt + Pe(-1/T2+j(Wo-W)tdtl S () = a{ f e-t/T2e-j ((O-W) t(dt)+ t e( 1/T2+j (o-w) idt} -1/T2-j(Wo-W) -1/T2+j(Wo-W) S(w) = 2a + T2 1+ (wo-w) 'T2 Le spectre est réel et ne comporte qu'un seul terme qui est un terme d'absorption égale à 2Ro(w). Dans le cas d'un champ inhomogène, il faut considérer un signal qui présente un argument supplémentaire 0 tel que: s'(t) = s(t) ejg, ej0 étant une constante pour chaque spin. Dans ce cas, cette phase 0 supplémentaire est celle qui est introduite par l'inhomogénéité du champ. En effet, si tous les noyaux ne voient pas le méme champ B, ils précessionnent à des vitesses angulaires différentes de celles que leur imposent: - le codage de phase Gy; - le codage en fréquence Gx lors de la lecture. Dans le cas idéal 0 est donc nul (pas d'inhomogénéités), dans le cas non idéal 0 est fonction de la position spatiale des spins. La transformation de Fourier donne S'(w) = S(W)eJ0-au lieu de S(w) étudié précédemment. a) si le signal n'est défini que de 0 à ± eJi = cos 0 + j sin 0 et S'(w) = (Ro(w) + j Io(w)) (cos 0 + j sin 0) = Ro(w)cos 0 - Io(w) sin 0 + j [Ro(w) sin 0 + Io(w) cos 0] On constate que la partie réelle et la partie imaginaire sont des combinaisons linéaires des termes d'absorption et de dispersion représentées sur les fig. 13 et 14 sur lesquelles: R(w) = Ro(w) cos Â - Io(w) sin f et I(w) = Ro(w) sin t + Io(w) cos Le spectre est donc obscurci par la présence de ces mélanges et il en résulte une perte notable de la résolution, c'est-à-dire un brouillage irrécupérable de ces images. B) si le signal est défini de -ô à S' (w) = 2Ro cos 0 + j2Ro sin 0 c'està-dire que les parties imaginaire I = 2Ro.sinT et réelle R = 2Ro cos Y sont des termes d'absorption purs modulés par cos 0 et sin 0 (voir Fig. 15 et 16). La valeur de 0 n'a donc pas besoin d'être connue pour que l'on effectue une correction de phase, car connaissant R et I, on en déduit: 0 = arc cos ( R ce qui, dans le procédé Me l'invention, revient à déterminer l'inhomogénéité de B en chaque point. Ces résultats mathématiques s'appliquent immédiatement au procédé d'obtention d'images RMN correctes en champ inhomogène comme expliqué maintenant. La variation de phase 0 (x,y) qui sera corrigée apparaît au cours de la période de codage de phase par le gradient Gy et s'écrit: 0(x,y) = yE(x,y) T o E(x,y) est l"'erreur" sur B au point x,y et T la durée d'application du gradient Gy. De plus, pendant la période de lecture (O-t), il apparaît un glissement de fréquence yE(x,y)/2e qui dépend malheureusement à la fois de x et de y pour le point (x,y) par rapport à la valeur observable en champ parfaitement uniforme. La phase (généralisée) acquise pendant Gx due à l'inhomogénéité E(x,y) est yE(x,y)t, et à chaque instant, introduit donc i) si l'acquisition a lieu de 0 à Y, un mélange des termes d'absorption et des termes de dispersion qui est irrécupérable si on ne connait pas E(x,y) à l'avance-et cecidans les parties réelle et imaginaire des transormées de Fourier par rapport au temps. ii) si l'acquisition a lieu de -ô a +, une partie réelle et une partie imaginaire débarrassées des termes de dispersion, et qui permettent de déter- miner l'écart par rapport à la valeur nominale de B en tout point x,y. Toutefois, le glissement en fréquences persiste, ce qui explique que les images sont déformées selon la direction ox du gradient de lecture. Aussi, faudra-t-il opérer en deux temps pour effectuer la correction. On construit tout d'abord l'image, tout en conservant les deux informations: réelle et imaginaire. On dispose ainsi de l'image réelle et de l'image imaginaire. Seule la direction x est affectée par des distorsions de forme, sur chaque image. 1) Repérons un point (x,y) sur cette image complexe. En ce point, la comparaison des valeurs fournies par l'image réelle (R) et l'image imaginaire (I) donne: E(x,y) par la relation: yE(x,y)T = arc cos R R2 + I2 * et l'intensité en ce point est R2 + 12 = Ro 2) La position y est correcte, mais non la position x qui est x', tel que: Gx + E(x,y) - Gx' = 0 o x' est la valeur de l'abscisse à laquelle il faudra affecter l'intensité Ro. Cette intensité-est donc affectée définitivement au couple (x',y) et la forme correcte de l'image, avec son signal image correct, est restituée. Quelques remarques ou explications complé- mentaires seront données pour terminer l'exposé de l'invention et de ses principales applications. a) Le procédé de l'invention demande un temps d'expo- sition double du temps que nécessiterait l'imagerie ordinaire. Cela peut être considéré comme un incon- vénient, qui semble relativement mineur en comparaison du progrès apporté dans la qualité de l'image. b) La valeur de l'erreur E(x,y,z), ou inhomogéneité, est déterminée par la méthode de restitution et non par une cartographie préalable. C'est un "sous-produit" du procédé et il peut avoir son utilité si l'on cherche à déterminer l'inhomogénéité introduite par la présence de l'échantillon. Par exemple, en imagerie RMN microscopique (objets inférieurs à 1 mm), on peut ainsi observer (tout en s'en débarrassant) les déformations des lignes du champ magnétique provoquées par la distribution de domaines dont les susceptibilités diamagnétiques sont différentes. On en déduit une méthode d'imagerie de la répartition des susceptibilités diamagnétiques; en outre, le procédé de l'invention évite de recourir à des valeurs de gradients très importantes dans les techniques de microscopie, valeurs justement imposées actuellement pour atténuer l'hété- rogénéité des objets observés. c) Un petit inconvénient persiste: le fait que le déphasage soit défini à 2w près, si bien que la connaissance de E s'effectue modulo Eo,o yEot = 2w. Toutefois, par continuité, il est pratiquement toujours possible de déterminer la valeur correcte de E. d) Il est de plus, possible, d'opérer en présence de gradients non constants, à condition que la variation spatiale du champ imposée par le gradient soit monotone (c'est-à-dire si elle ne crée pas de ddoublements de la valeur du champ lors de leur utilisation). En effet, appliquer un gradient non constant revient à superposer une inhomogénéité supplémentaire E'(x,y). Cependant, la détermination correcte de E(x,y) à partir de la valeur obtenue Eo(x,y) = E(x,y) + E'(x,y) impose la connaissance de E'(x,y). e) Le procédé s'applique également à toute opération d'imagerie effectuée dans un champ considéré comme "uniforme par les constructeurs. AB/Bo étant donné par le constructeur, la limite de résolution spatiale dans l'image peut être estimée égale à Al = AB/G ou G est la valeur absolue du gradient de lecture utilisée (Gx). Ainsi, il semble possible de reculer la limite de résolution spatiale accessible couramment dans le plan de l'image à une valeur nettement inférieure sans utiliser des valeurs de gradient nettement plus élevées que celles qui sont employées. Il en résulte une nette amélioration dans la manière d'appliquer ces gradients: alimentations en courant plus -faibles, consommation électrique réduite. Enfin, le procédé selon l'invention s'applique à des nombreux types d'imageries RMN, notamment: - les imageries 2DFT et 3DFT; - l'imagerie microscopique; - les séquences "saturation-récupération", inversion récupération, spin écho, les séquences spécifiques à l'étude des mouvements liquidiens et les séquences réalisant le contraste eau/graisses en imagerie des protons sur échantillons biologiques. REVENDICATIONS 1. Procédé de génération et de traitement de signaux pour l'obtention d'une image de résonance magnétique nucléaire (RMN) exempte de distorsions à partir d'un champ de polarisation inhomogène. B=Bo + E(x,y, z), utilisant les méthodes connues de codage de l'image en deux ou trois dimensions (en terminologie anglosaxonne "2DFT" et "'3DFT") comportant l'application d'une impulsion de radiofréquence à 90 dans le plan xOy, puis un codage de phase pendant un temps T suivi d'un codage de fréquence par application selon Ox d'un gradient de lecture pendant l'acquisition du signal de précession libre, signal auquel on fait subir ensuite des transformations de Fourier par rapport au temps et par rapport au(x) gradient(s) appliqué(s) et que l'on recueille sous forme complexe, caractérisé par les étapes suivantes: a) au cours d'une première série de séquences d'acquisition des signaux comportant un codage de phase et un codage de fréquence, on recueille de la façon connue précédente les n signaux complexes S1(T,t); b) au cours d'une deuxième série de séquences comportant l'application successive des mêmes gradients, on soumet le système en précession libre à une seconde impulsion de radiofréquence de 180 au début de l'éta- blissement de chaque gradient de lecture, de façon à changer en -0 (ou r 0 selon le cas) les phases 0 acquises par les différents moments magnétiques durant le codage de phase et l'on recueille les n signaux imaginaires s2 (T, t), les étapes a) et b) -précédentes pouvant être successives ou imbriquées dans le temps; c) on forme les signaux s2 (T,t) conjugués de s2 (t,t) qui sont tels que s2 (t,t) = s1(T, -t), et on les associe un.par un avec leur homologue sl (T,t) de façon à former n signaux complexes s(t,t) dans lesquels t varie de -T à +T, T étant le temps d'acquisition de chaque signal compté à partir du début de l'application du gradient de lecture; d) on effectue sur ces signaux s{(T, t), pour. lesquels--T appliqué(s), ce qui fournit une-image complexe multi- dimensionnelle dont les parties réelle R et imaginaire I selon la direction x ne comportent que des termes d'absorption pure; e) on effectue sur chaque point x,y,z, les corrections de distorsion selon Ox et d'affectation d'intensité d'aimantation en calculant successivement la valeur correcte de la phase: R = Arc cos- R 42 + 2 R+ I la fonction E(x,y,z) représentant la cartographie du champ par la relation 0(x,y,z) = yE(x,y,z)T puis par la relation Gx.x + E(x,y,z) - Gx. x' = 0, la valeur x' de x à laquelle doit être affectée l'intensité d'aiman- tation trouvée. 2. Procédé d'obtention et de traitement de signaux selon la revendication 1, caractérisé en ce que la seconde impulsion de radiofréquence de est en quadrature avec la première appliquée au début du codage de phase et l'on obtient des signaux s2 (T, t) = S l(T; - t). 3. Procédé d'obtention et de traitement de signaux selon la revendication 1, caractérisé en ce que la seconde impulsion de radiofréquence de 1800 est en phase avec la première appliquée au début du codage de phase et l'on obtient des signaux s2 (T,t) = - S1(T, -t), que l'on change de signe avant de les associer à leur homologue de la première séquence.