Brevet FR2586296A1 — publié le 19870220

PROCEDE DE MODULATION DE L'EFFET DE LA VITESSE DES PARTIES MOBILES D'UN CORPS DANS UNE MESURE DE DENSITE PAR RESONANCE MAGNETIOUE NUCLEAIRE, ET MISE EN OEUVRE DU PROCEDE POUR EN DEDUIRE LA VITESSE DES PARTIES MOBILES CONCERNEES. La présente invention a pour objet un procédé de modulation de l'effet de la vitesse des parties mobiles d'un corps dans une mesure de densité par résonance magnétique nucléaire (RMN) ainsi que la mise en oeuvre du procédé pour en déduire la vitesse des parties mobiles concernées. L'invention trouve son application parti- culièrement dans le domaine médical o les corps examinés sont des corps humains et o les parties mobiles sont des cellules du sang circulant dans des veines ou des artères ou des organes mobiles tels que le muscle cardiaque. Dans cette application, l'invention peut même plus particulièrement être mise en oeuvre avec un procédé d'imagerie pour donner une image représentative de la répartition des vitesses des parties mobiles dans une coupe du corps examine. En effet, l'imagerie par résonance magnétique nucléaire se développe principalement comme moyen de diagnostic médical. Elle permet la visualisation des structures internes tissulaires avec un- contraste et une résolution jamais atteints simultanément par les autres procédés d'imagerie. Pour obtenir une image par résonance magnétique nucléaire d'une coupe d'un corps avec différenciation des caractéristiques tissulaires de ce corps, on utilise la propriété qu'ont certaines particules, comme les protons, d'orienter leur moment magnétique en acquérant de l'énergie lorsqu'ils sont placés dans un champ magnétique principal constant B0. Une zone parti- culière d'un corps contenant des particules présente alors un moment magnétique global que l'on peut faire basculer selon une orientation donnée, perpendiculaire ou parallèle au champ B0, en induisant une résonance par émission d'un champ magnétique radio- fréquence perpendiculaire au champ principal. Toutes les particules possédant alors un moment magnétique, tournant à une vitesse de précession dite de Larmor, tendent à retrouver l'orientation initiale parallèle à B0 en émettant un signal radiofréquence à la fréquence de résonance caractéristique de B0 et de la particule. Ce signal peut être capté par une antenne de réception. La durée de retour à l'équilibre du moment magnétique global d'une région considérée et la décroissance du signal dépendent de deux facteurs importants: l'interaction spin-réseau et l'inter- action spin-spin des particules avec la matière environante. Ces deux facteurs conduisent à la définition de deux temps de relaxation caractéristiques du tissu et dits respectivement T1 et T2. Une région considérée d'un objet émet donc un signal dont l'intensité dépend de T1, de T2, de la densité de particules dans la région, et du temps qui s'est écoulé depuis l'excitation radiofréquence. Si le champ orientateur est parfaitement homogène, des particules mobiles dans une région considérée émettent en réponse un signal identique à celui des particules fixes de cette région. Par contre si le champ orientateur n'est pas homogène, ou plus géné- ralement si pour diverses raisons (en particulier pour faire de l'imagerie) on applique pendant ou après l'excitation magnétique radiofréquence un champ magnétique perturbateur, présentant un gradient d'intensité, on peut montrer que les contributions apportées par les particules mobiles dans le signal global émis sont affectées d'uhe composante de phase dépendant justement de leur vitesse. Ceci peut se comprendre simplement: le signal de résonance émis vibre à une fréquence f0 qui dépend de l'intensité du champ ma- gnétique orientateur B0 et d'un rapport dit gyromagnétique carac- téristique du milieu examiné. Toutes variations de l'intensité du champ B0 provoquent donc une variation correspondante de la fréquence de résonance. En conséquence, une particule fixe soumise, après l'excitation radiofréquence, pendant un premier temps au champ Bo résonne à une fréquence f0, dans un deuxième temps, soumise à un champ plus fort B0 + B0, elle résonne à une fréquence plus élevée f + fO' Dans un troisième temps à nouveau soumise au champ B0 elle vibre à nouveau à la fréquence f0. Au cours du troisième temps le signal émis est alors déphasé par rapport à sa phase pendant le premier temps. Ce déphasage est proportionnel à l'amplitude de la perturbation l Bo et à la durée de cette perturbation. Si toutes les particules du milieu sont fixes ou, si la perturbation qui a atteint dans le temps tout le milieu n'a pas de gradient, il en résulte simplement que le signal global émis est seulement retardé. Il en va différemment par contre pour les particules animées i0 d'une certaine vitesse quand la perturbation présente un gradient. Celles-ci fréquentent au cours des trois périodes (du fait de leur vitesse de déplacement pendant ces périodes) des régions de l'espace o les champs orientateurs et perturbateurs sont différents. Ils sont différents respectivement du fait de l'existence d'inhomogénéités ou du fait de l'existence de gradients. En conséquence la contribution des particules mobiles dans le signal se trouve apportée avec une phase qui dépend non seulement de l'amplitude des perturbations rencontrées (comme pour les particules fixes) mais aussi de la variation d'amplitude de ces perturbations le long du parcours qu'elles ont emprunté. Cette variation qui constitue le gradient est géographiquement imposée. En conséquence, le déphasage du signal des particules mobiles dépend alors de leur vitesse puisque elles fréquentent d'autant plus de régions de l'espace que leur vitesse est grande. Si les vitesses de déplacement, les inhomogénéités, ou les gradients de champ sont trop grands, les phases des différentes contributions peuvent en être à ce point affectées qu'elles finissent par s'opposer. Dans ce cas, ces contributions s'annulent mutuel- lement et le signal global résultant est moins fort. Dans la pratique cet effet est tel qu'il donne souvent l'illusion qu'il n'y a pas de matière dans un corps à l'endroit o circulent des particules mobiles. Pour faire apparaître l'existence des particules mobiles et pour en mesurer les caractéristiques, la densité et éventuellement la vitesse de déplacement, on peut procéder selon une méthode décrite par E.L.HAHN en Février 1960 dans le Journal of GEOPHYSICAL RESEARCH, volume 65, n 2, pages 776 et suivantes. L'auteur y suggère de soumettre le milieu étudié à une séquence d'un gradient particulier: à le coder. Le principe de ce codage consiste à appliquer après le basculement de l'impulsion radiofréquence, un gradient dit bipolaire selon l'axe d'une composante de vitesse que l'on veut reconnaître. Un gradient bipolaire est tel que son intégrale temporelle, du temps correspondant au début de l'impulsion radio- fréquence au temps correspondant au début de la mesure, est nulle. Le moment magnétique de spin d'une particule immobile ne subit dans ce cas qu'un déphasage global nul. En effet le déphasage subit pendant l'application d'une première partie du gradient bipolaire est compensé par l'application de la deuxième partie de ce gradient bipolaire. Par contre une particule mobile, ayant une vitesse posi- tive suivant l'axe du gradient, subit lors de la deuxième partie de l'impulsion un déphasage plus important en valeur absolue que pendant la première partie. Parce qu'elle fréquente pendant cette deuxième partie une région de l'espace o, du fait du gradient le champ magnétique perturbateur est plus fort. En comparant une mesure relevée avec un tel gradient bipolaire à une mesure relevée sans qu'il ait été appliqué, on peut en déduire la vitesse et le nombre des particules mobiles. Ouels que soient les buts poursuivis, mesure simple ou mesure avec image, et quelles que soient les procédures retenues, la sensibilité du phénomène de vitesse au champ magnétique pertur- bateur appliqué est telle, qu'en définitive, la mise en évidence des phénomènes de déplacement n'est possible que tant que les vitesses maximum sont inférieures à une limite. En particulier en imagerie, selon que la composante de vitesse à mettre en évidence est parallèle ou perpendiculaire au plan de la coupe imagée, la sensibi- lité des machines de RMN est actuellement de l'ordre de 1 ra- dian/(cm/s) à 0,2 rd/(cm/s). Ceci signifie qu'une particule qui se déplace de I cm par seconde dans le plan de la coupe contribue au signal global emis avec un déphasage de 1 radian par rapport aux contributions émises par les particules fixes. Dans le corps humain une vitesse nominale de circulation du sang de 50 cm/s est couram- ment atteinte, voire même de plusieurs mètres par seconde dans le coeur. De plus la distribution des vitesses dans un vaisseau est étalée entre zéro, sur les bords du vaisseau, et la vitesse nominale au centre du vaisseau. Donc chaque particule d'un vaisseau contribue au signal avec un déphasage qui peut valoir de 0 à 50 radians! Sachant que des contributions déphasées de Ir radians s'opposent mutuellement, le signal résultant est nul: cela revient à prendre la valeur moyenne d'un signal sinusoïdal sur plusieurs périodes. A titre d'exemple, Paul R. Moran dans un article de la revue Radiology du RSNA de 1985, 154, pages 433 à 441, fait état d'une mesure de vitesse moyenne égale à 0,6 cm par seconde et correspondant à un déphasage d'environ 90 . Au delà de cette limite, la sensibilité des machines est trop forte et les vitesses ne sont plus mesurables. L'invention a pour objet de mettre en évidence l'effet de la vitesse des parties mobiles d'un corps en modifiant la sensibilité des machines d'une manière particulière. La modification de sensibilité, tout en ne modifiant d'aucune manière les signaux émis par les parties fixes, peut avoir pour effet d'annuler la part de déphasage dû à la vitesse. Les particules en mouvement contribuent alors au signal global comme si elles étaient fixes. Le procédé de l'invention permet par ailleurs, plutôt que de tout annuler, de moduler l'effet de la yitesse. En faisant deux mesures, avec des caractéristiques de modulation différentes, et en comparant les deux mesures, on peut supprimer l'incidence des parties fixes pour ne laisser apparaître que celle des parties mobiles. Dans cette comparaison les parties mobiles apparaissent pondérées par les caractéristiques de modu- lation des deux expérimentations. Dans l'invention ces caracté- ristiques sont calculables et l'effet de la vitesse peut être quantifié. Le principe énoncé dans le procédé de l'invention est très. général. Il n'est pas cantoné à une application, plus complète, d'imagerie. Il peut en particulier être mis en oeuvre avec des spectromètres de masse. Par ailleurs, il s'applique quelles que soient les procédures d'excitation radiofréquence retenues: procédures connues souvent sous leurs noms anglo-saxon de saturation recovery, inversion recovery, double saturation recovery, saturation recovery - inversion recovery, spin-echo, etc... De même l'invention, dans le cadre de l'imagerie, est applicable quelle que soit la procédure d'imagerie retenue: par exemple que celle-ci soit une méthode de retroprojection de type P.C. Lauterbur, une méthode par trans- formée de Fourier de type 3DFT ou 2DFT mise au point par A. Kumar et R.R. Ernst ou sa variante connue sous le nom de spin warp (c'est dans cette seconde méthode que l'invention sera décrite), une méthode d'imagerie des volumes sensibles mise au point par W.F. Hinshaw, ou encore une méthode d'acquisition rapide mise au point par P. Mansfield et plus connue sous le nom de "Echo-planar", etc... En effet, dans toutes les situations impliquées par ces procédures, la modulation objet de l'invention est possible puisqu'elle consiste à modifier les champs magnétiques perturbateurs (les gradients de champ) en leur adjoignant des champs magnétiques compensateurs d'allure semblable, et dont les caractéristiques de forme, de durée, et d'amplitude dépendent de ces champs magnétique perturbateurs. L'invention a pour objet un procédé de modulation de l'effet de la vitesse des parties mobiles d'un corps dans une mesure de densité par résonance magnétique nucléaire, mesure pour laquelle - on soumet le corps à un champ magnétique orientateur et constant pour orienter dans une direction unique des moments magnétiques de spin du corps, - on fait subir à ce corps une excitation magnétique radio- fréquence en présence et/ou suivie de l'application d'une séquence d'un champ magnétique perturbateur, - et on relève un signal de résonance magnétique émis en réponse par le corps, ce procédé étant caractérisé en ce qu'on module l'effet de la vitesse des parties mobiles du corps créé par la séquence du champ perturbateur, par l'application, avant le relevé, d'une séquence d'un champ magnétique compensateur, dont l'intégrale calculée sur sa durée est nulle, et dont l'historique et la valeur sont fonction de l'historique et de la valeur du champ perturbateur. L'invention concerne également une mise en oeuvre de ce procédé, caractérisée en ce qu'on le met une première fois en oeuvre pour compenser l'effet de la vitesse, et une autre fois en oeuvre en modifiant une ou plusieurs des caractéristiques du champ magnétique compensateur et en ce qu'on compare les mesures obtenues dans les deux mises en oeuvre pour en déduire la vitesse des parties mobiles concernées. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Comme cela a I0 été évoqué plus haut cette description et ces dessins ne sont donnés qu'à titre indicatif et nullement limitatif de l'invention. En parti- culier l'invention n'est décrite que dans une application d'imagerie. Sur les figures les mêmes repères désignent les mêmes éléments. Elles montrent: - figure 1: une machine pour la mise en oeuvre du procédé selon l'invention; - figure 2 des diagrammes temporels de signaux radiofré- quence d'excitation, de signaux de champ magnétique perturbateur, et de signaux relevés dans une mesure particulière comportant une imagerie de type 2DFT d'une coupe d'un corps sous examen; - figures 3a et 3b: des diagrammes temporels de déphasages résultant, après application de séquences de champs magnétiques perturbateurs, entre les contributions émises par des particules fixes. et des particules en mouvement; - figure 4: une représentation schématique de la réponse d'une partie d'un milieu dont les particules sont en déplacement selon que ce déplacement est parallèle ou perpendiculaire à un plan de coupe imagé. La figure I représente une machine utilisable pour mettre en oeuvre le procédé de l'invention. Elle comporte des moyens symbo- lisés par une bobine 1 pour soumettre un corps 2 à un champ magnétique important et constant B0. Le champ B0 est le champ orientateur. La machine comporte encore des moyens générateurs 3 et bobines 4, reliés aux moyens générateurs 3, pour soumettre le corps ainsi placé à des séquences d'excitation radiofréquence en présence de séquences de champ perturbateur: des gradients de champ orientés selon trois axes de référence de la machine X, Y ou Z (figure 2). Les bobines 4 représentent en même temps les bobines radiofréquence et les bobines de gradient de champ. La machine comporte encore des moyens de réception 5, reliés aux bobines 4, pour recevoir le signal de résonance magnétique. Dans une application d'imagerie, des moyens 6 peuvent permettre de calculer et de mémoriser une première image 1 et une deuxième image 12 d'une coupe 7 du corps 2. Les deux images sont relatives à deux séries d'expérimentation imposées par des commandes symbo- lisées C1 et C2 des moyens générateurs 3. Dans des moyens de comparaison 8 on compare point à point les images Il, I2, et on élabore une troisième image 13. On peut visualiser cette image 13 sur un dispositif de visualisation 50. Dans les images I.1 I2, ou 13, un point d'image, ou pixel, représente classiquement par sa luminosité la densité de particules contenue dans un élément de volume, ou voxel v qui lui correspond dans la tranche 7. Le voxel v comporte de nombreuses particules animées de vitesses V comportant des composantes Vx, Vy et Vz sur chacun des trois axes de référence de la machine. L'invention va permettre de calculer dans le voxel v, et pour tous les autres voxels aussi, la vitesse moyenne et le nombre des particules de ce voxel ainsi en mouvement. L'image 13 est représentative de cette infor- mation de vitesse, par comparaison entre deux images I1 et 12 o l'effet de cette vitesse aura été modulée de deux manières dif- férentes. Dans un premier temps, on va rappeler brièvement la théorie de l'imagerie par une méthode de type 2DFT. Dans cette méthode comme dans toutes les autres il existe des gradients perturbateurs. Dans un deuxième temps on va rappeler pourquoi cette méthode est préférée eu égard à des problèmes d'inhomogénéité du champ orientateur. Dans un troisième temps, on indiquera, dans cet exemple, comment se calculent les champs compensateurs néces- saires pour moduler l'effet de la vitesse. Pour localiser une région d'un milieu il est nécessaire de repérer la nature de son émission en fonction de conditions locales de champ magnétique. Ces conditions locales sont imposées de façon à ce que la fréquence et la phase d'émission soient bien caracté- ristiques de la localisation dans l'espace de cette région du milieu. Pour cela on superpose au champ principal B0 des gradients de champ magnétique pulsés. Ces gradients sont orientés selon des directions X, Y, Z pour définir à tout moment les éléments de volume qui résonnent à des fréquences connues. Pour l'acquisition d'une image entière les conditions locales sont imposées en des séquences programmées (C1, C2). Celles-ci sont mises en mémoire dans un ordinateur pilote. Ces séquences définissent les instants d'application des gradients, les instants d'excitation des particules par les impulsions du champ radiofréquence, ainsi que les instants de lecture ou d'acquisition des données de l'image. La méthode d'imagerie 2DFT permet actuellement d'obtenir la meilleure qualité d'image. Dans cette méthode d'imagerie seul un plan de coupe, le plan 7 est, excité à la fois. La figure 2 montre à cet effet des impulsions radiofréquence 10, 11 et 12 appliquées en présence d'un gradient de champ dit de sélection, par exemple orienté selon l'axe Z, et représenté respectivement par des impul- sions 13, 14 et 15. Quand le gradient de sélection est orienté selon l'axe Z la coupe est transverse, c'est-à-dire selon un plan X, Y. Dans l'imagerie 2 DFT on code en phase les différents signaux acquis. Ceci est obtenu par une impulsion, d'intensité variable d'un gradient dit déphaseur d'axe perpendiculaire à un gradient dit de lecture dont la direction est constante. Par exemple, pour une coupe transverse, le gradient de lecture peut être sur l'axe X et le gradient de déphasage peut être sur l'axe Y. Sur la figure 2, le gradient X est coté 16, le gradient Y est coté 17. Par une double transformé de Fourier spatiale l'image est reconstruite, d'o le nom de la méthode. Une amélioration de cette méthode peut d'ailleurs permettre d'ob- tenir simultanément des images de plusieurs couples parallèles. L'existence du précodage 17, qui prend un certain temps, repousse dans le temps la mesure du signal émis par le corps. Comme ce signal s'atténue très vite, on a imaginé d'en mesurer un écho. Cet écho est produit, au moment de l'application des impul- sions radiofréquence Il puis 12 par une réflexion de la dispersion de phase des contributions apportées par chacune des particules. Pour cette raison ces impulsions sont dessinées sur la figure 2 plus hautes que l'impulsion 10 puisqu'elles font basculer l'orientation des moments magnétiques des spins des particules de 180 alors que l'impulsion 10 ne les avait fait basculer que de 90 . Pour imager une coupe d'un corps il y a lieu de réaliser des séquences d'excitation d'écho de spin en présence de séquences de gradients de champ en un nombre aussi important de fois que la résolution de l'image attendue doit être plus précise. A chaque séquence d'excitation, le gradient déphaseur Y varie par pas successifs partant d'une certaine valeur jusqu'à la même valeur mais de signe opposé. Cette valeur dépend de la forme et de la durée du gradient 16 de lecture. Ce gradient déphaseur permet de faire tourner chaque moment magnétique de spin d'une phase variable, dépendant de son ordonnée suivant l'axe Y et de la valeur de ce gradient. Pour chaque image Il et 12 le gradient Y peut prendre successivement le même nombre de valeurs : d'une manière préférée la définition des deux images sera la même. Pour éviter un précodage parasite imposé par la fin 18 de l'impulsion 13, pour la partie de cette impulsion existant après la fin de l'impulsion radiofréquence 10, et par le début 19 de l'impulsion 16 de lecture il est connu d'appliquer sur ces axes respectivement des impulsions 20 et 21. L'intégrale temporelle de ces impulsions est la même mais de signe opposé, elles neutralisent l'effet parasite de codage. L'intégrale est entendue au sens d'intégrale dans le temps. Les surfaces hachurées 20 et 19 sont donc égales aux surfaces en pointillés 18 et 21 respectivement. Les impulsion 20 et 21 n'inter- fèrent pas dans la sélection de coupe puisqu'au moment o elles sont appliquées l'impulsion 10 n'est plus présente. Les impulsions 14 et 15 sont autoneutralisées puisqu'elles comportent des portions respec- Il tivement 22 et 23, et 24 et 25 antisymétriques par rapport au milieu de l'impulsion radiofréquence, Il et 12 pour lesquelles elles resélec- tionnent la coupe 7 dans le corps 2. La méthode d'imagerie de type 2 DFT est la plus utilisée dans la pratique du fait de sa rapidité d'acquisition (comparée à des méthodes 3D) et à sa robustesse vis-à-vis des imperfections du système physique, en particulier vis-à-vis des inhomogénéités du champ orientateur. Ceci s'apprécie particulièrement en comparaison avec des méthodes par rétroprojection ou avec des méthodes n'uti- lisant pas d'écho de spin. Cette méthode est en contrepartie la plus intrinsèquement sensible au déplacement des particules. L'immunité aux inhomogénéités du champ B0 de la méthode 2 DFT provient du fait que contrairement aux méthodes par rétroprojection le gradient codeur en fréquence est de même direction (le gradient de lec- ture X) et de même amplitude, GL, d'une séquence à l'autre. L'inho- mogénéité du champ B0 ne fait alors que déformer les lignes isofréquences de résonance. Pendant la lecture du signal, l'inhomo- généité du champ devient l'équivalent d'une mauvaise linéarité du gradient de lecture. Les images sont en conséquence distordues mais sans la création de flou et sans la perte de résolution associées que l'on rencontre avec les méthodes par rétroprojection. Dans celles-ci en effet le gradient de lecture change d'orientation à chaque séquence et de ce fait il répartit les inhomogénéités du champ orientateur dans toute l'image pour laquelle il crée ainsi un flou. Les figures 3a et 3b montrent le déphasage à0 résultant dans la contribution de signal émis par deux particules d'un même voxel soumises à des perturbations identiques, quand l'une des particules (trait plein) est immobile alors que l'autre particule (tirets) est animée d'une vitesse V. La figure 3a montre le diagramme temporel d'un signal radiofréquence RF comportant une impulsion 26 pour créer un écho de spin. Selon un axe A un gradient de champ magnétique a été imposé en prenant avant et après l'impulsion 26, par exemple, une valeur G pendant une durée e à chaque fois. En un lieu donné 10 une particule préalablement excitée, immobile, vibre à une fréquence fO au début de l'application de l'impulsion 27. Pendant la durée Ode cette impulsion 27 elle vibre à une fréquence dif- férente. Par exemple cette fréquence est supérieure, de telle manière qu'à la fin de l'impulsion 27 on puisse considérer, puisqu'elle se remet à vibrer à la fréquence fO qu'elle a subi un déphasage t.. Lors de l'application de l'impulsion radiofréquence 26 on peut ad- mettre que tout se passe comme si la phase du signal émis par cette particule s'inversait: - 1' Pendant l'impulsion 28 qui a exac- tement la même allure que l'impulsion 27 le phénomène de dépha- sage continue: la contribution se déphase encore de (1' Mais ce dernier déphasage vient effacer l'effet du premier après qu'il ait été inversé par l'impulsion 26. En conséquence, la particule immobile en 10 retrouve à la fin de l'impulsion 28 la phase qu'elle avait au début de l'impulsion 27. Il est utile de remarquer que le dépha- sage f 1 est proportionnel d'une part à l'amplitude du gradient G et d'autre part à la durée pendant laquelle chaque impulsion du gradient est imposée. Plus généralement qc est proportionnel à l'intégrale de ce gradient pendant la durée d'existence de son impulsion si cette impulsion n'est pas exactement carrée. Il en va tout autrement pour une deuxième particule, toute voisine dans le même voxel de la particule précédente, et qui par contre est animée d'une vitesse V. Sur la figure 3a, on a indiqué que cette particule se trouvait successivement dans les positions 10 à 14 au début et à la fin de la première impulsion 27, au moment de l'impulsion radiofréquence 26, et au début et à la fin de la deuxième impulsion 28 du gradient. Après la première impulsion 27 de gra- dient, le déphasage résultant P2 est proportionnel au gradient G et à la durée O.Mais il est aussi proportionnel à la distance 11 - 10 parcourue dans l'intervalle par la particule. Cette distance est égale à V. 0. Ce qui fait que le temps intervient maintenant au carré: la première portion de déphasage à l'allure d'une parabole (V.G. 82). De la fin de l'impulsion 27 au début de l'impulsion 28, tout se passe comme pour la particule immobile: la contribution possède une phase -)2 au début de l'impulsion 28. Au cours de cette impul- sion 28 l'évolution de la phase est également parabolique comme pour l'impulsion 27. On pourrait s'attendre à ce que les effets de phase s'effacent ici également. Cependant entre la fin de l'impul- sion 27 et le début de l'impulsion 28 la particule animée s'est déplacée de la position 11 à la position 13. Dans la position 13 au moment o l'impulsion 28 apparaît, et du fait de l'existence du gradient de champ, les conditions locales de champ ne sont plus les mêmes que celles qui avaient prévalu quand la particule s'était trouvée entre 10 et 11. En conséquence à la fin de l'impulsion 28 la compensation ne s'opère pas et cette particule animée émet alors un signal déphasé d'une valeur Dû que l'on peut écrire de la manière suivante: = G (). V. T. o r est le rapport gyromagnétique des particules considérées, o G (e)représente le gradient dans les impulsions 27 et 28, et o T est la durée qui sépare l'apparition de ces deux impulsions. Si e n'est pas petit, T mesure la durée entre les milieux de ces impulsions. Plus généralement on peut écrire: A (t) = V G v. d. o' le déphasage est donné maintenant en fonction du temps et o G (t) représente toutes les perturbations appliquées aux parti- cules avant qu'on en mesure à une date t le signal qu'elles émettent. On se souvient simplement que les perturbations considérées sont telles que fG (t') dt= 0 on Car c'est sous cette condition que, d'une manière connue-dans l'état de la technique (figure 2 impulsions 18 à 25), on a compensé les précodages parasites. On ne l'avait bien entendu pas fait pour l'impulsion variable 17 puisque le déphasage qu'elle provoque pour le signal global est, d'une acquisition à l'autre, une des caractéristiques permettant l'imagerie. En l'absence d'une impulsion radiofréquence telle que 26 cré- ant un écho de spin (figure 3b), les impulsion 27 et 28 du gradient bipolaire doivent être remplacées par des impulsions 29, 30 de signes opposés. De cette manière la phase pour les particules fixes évolue jusqu'à 903 en fin de l'impulsion 29 et évolue en sens inverse jusqu'à zéro au cours de l'impulsion inverse 30. Le déphasage tAf entre la contribution des particules mobiles et celles des particules fixes est du même ordre que précédemment. Dans l'expression de A(t), G (e) représente donc toutes les séquences nécessitées par la mesure. Ce qui caractérise l'invention c'est que ces séquences sont maintenues telles quelles mais que leurs effets sur la vitesse sont compensés par l'adjonction de séquences supplémentaires de gradients de champ G' dits compensateurs tels que Ut (G (t).r+ G'(t.r) d= 0 o G' représente des gradients, eux aussi bipolaires, c'est-à-dire tels que ; G(t")dt'= 0 d Autrement dit en agissant ainsi on annule le coefficient multi- plicateur de la vitesse V qui apparaît dans le calcul de Af(t). Du coup, quelle que soit cette vitesse, Af(t) est nul: toutes les contributions sont en phase. Les particules mobiles sont prises en compte comme si elles étaient fixes. En revenant à la figure 2 on voit sur chaque axe X, Y, Z l'adjonction de tels gradients bipolaires compensateurs sous la forme préférée de couples d'impulsions notés 31 et 32, 33 et 34, 35 et 36, 37 et 38, 39 et 40 et 42 et 43. L'incidence de ces impulsions compensatrices est symbolisée dans le dessin par la présence de petites croix. On sait, dans les machines de l'état de la technique imposer de telles impulsions en modifiant les commandes C1 ou C2 des moyens générateurs. L'intégrale dans le temps des impulsions de chaque couple est nulle. D'une manière préférée, puisque l'impulsion de correction de précodage 20 sur l'axe de sélection Z est appliquée entre des dates tl et t2 à la fin de l'impulsion 13, on applique l'impulsion 31 pendant la même période. Dans ce cas l'impulsion 20 "nouvelle" est plus forte. Les impulsions 35 et 37 sont également appliquées pendant cette période. Les impulsions 32, 36 et 38 d'une manière préférée sont appliquées entre la fin t3 de l'impulsion 14 et le début t4 de l'impulsion de lecture 16. Les impulsions 33, 39 et 42 sont appliquées entre la fin t5 de l'impulsion 16 et le début t6 de l'impulsion 15. Les impulsions 34, 40 et 43 sont appliquées entre la fin t7 de l'impulsion 15 et le début t8 d'une impulsion 41 qui sert à la lecture d'un deuxième écho. Sur le plan pratique, la mesure du signal émis est effectuée entre des dates t'4 et t'5 contenues entre les dates t4 et t5 et aussi pendant la durée de l'impulsion 41. Bien entendu les impulsions de compensation appliquées sur un axe, X, Y ou Z ont pour objet d'effacer l'effet des gradients perturbateurs appliqués sur le même axe. L'expression intégrale donne implicitement les gradients G'; il est possible de faire des simplifications. Premièrement tous les termes en G (t).t sont des termes connus (ils appartiennent à la séquence 2DFT): leur intégrale l'est également. Cette intégrale peut donc être sortie de l'expression. D'une manière préférée les formes des impulsions compensatrices 31 à 43 seront des formes connues et, dans la mesure du possible, identiques à des formes déjà pratiquées pour les impulsions normales des séquences pertur- batrices. Compte-tenu des dates retenues pour l'application de ces impulsions, une intégrale unitaire de chacune d'elle peut être calculée et l'expression de G' peut s'écrire d'une autre manière: G (t). t'. d+ A). 'd t:= O Le premier membre est connu il concerne les gradients pertur- bateurs. Le deuxième membre comporte des intégrales de gradients de forme préféréeLeet d'amplitude unité. On sait calculer l'intégrale du deuxième membre puisqu'on connaît la forme choisie f. En conséquence on peut en déduire X. c'est-à-dire l'amplitude des impulsions compensatrices. Dans le cas de l'imagerie conventionnelle, o un gradient au moins est appliqué pendant la période de lecture (t'4, t'5), l'équation implicite contenant G' ne peut être vérifiée pour tout t. Mais on doit au moins la vérifier aux temps tE tE... centre des phases de as mesure. Ces dates sont celles ou l'éecho de spin attendu est le plus fort. Elles sont telles que la durée séparant la date de l'impulsion radiofréquence 10 et tE est le double de la durée entre l'impul- I sion 10 et l'impulsion 11, et telles que la durée entre la date tE et I5 la date tE est égale au double de la durée entre la durée tE et la date de I impulsion radiofréquence 12.... La justification Je ces dates est, d'une part, intuitive et, d'autre part, théorique. L'idée conductrice en ce domaine est que l'on cherche à connaître l'état de la coupe 7 non pas au moment de son excitation (impulsion 10) mais au moment, tE ou tE, o justement on effectue la mesure: au centre du signal émis. 2 Tous les termes de l'équation implicite donnant G' étant con- nus, G' est calculable. Sur la figure 2 on va donner une approche simple de ce calcul. Pour l'axe X on sait que les impulsions 21 et 19 situées de part et d'autre de l'impulsion radiofréquence 11 neutra- lisent mutuellement leurs effets de précodage: l'intégrale de la somme de ces impulsions est nulle. Par contre l'intégrale du produit de ces impulsions par la durée qui les sépare de la date du relevé, tE, n'est pas nulle: l'impulsion 21 étant plus antérieure dans le temps que l'impulsion 19. Tout se passe donc comme si l'impul- sion 21 de "force" donnée est appliquée à un "bras de levier" plus long que l'impulsion 19 de force égale. L'effet sur la phase de l'impulsion 19 est bien entendu inverse de celui de l'impulsion 21 puisqu'elles sont situées de part et d'autre de l'impulsion 11. En conséquence leur effet commun sur la phase s'analyse comme un moment orienté dans le sens de l'impulsion ayant le plus grand bras de levier: l'impulsion 21. Il convient de contrecarrer ce moment résultant par un couple d'impulsions 35, 36 dont l'intégrale totale est nulle (il ne provoque aucun déphasage sur le signal des parties fixes) mais dont le moment global s'oppose au moment global des deux impulsions 21 et 19. Ainsi, l'impulsion 35 étant placée dans le temps à une date plus antérieure que l'impulsion 36 par rapport à la- date tE son effet est prépondérant. Il est bien de sens opposé à l'impulsion 21. De la même façon sur l'axe de sélection Z, des impulsions 18 et 20, et 22 et 23, ce sont à chaque fois les impulsions 18 et 22 qui priment: leur bras de levier est plus long. En conséquence on conçoit bien que l'impulsion 31 doit être dans le même sens que l'impulsion 20. Le cas des impulsions 37 et 38 est un peu différent. En effet elles doivent contrecarrer les effets d'un gradient qui varie d'une acquisition à l'autre. Les amplitudes des impulsions 37 et 38 varieront donc également en correspondance d'une acquisition à l'autre (symbolisé par des flèches). Pour les raisons précédentes quand l'impulsion 17 est positive l'impulsion 37 doit être négative. Pour pouvoir réduire les amplitudes des impulsions compensatrices on cherche bien entendu à donner à leur moment le plus d'efficacité possible: on essaie d'augmenter la durée qui sépare leur application de la date du relevé. On les place le plus en amont possible de la séquence pour la première (31, 35, 37) et le plus en aval pour la deuxième (32, 36, 38). En conclusion elles sont le plus éloignées possible de part et d'autre de l'impulsion radiofréquence 11; et elles sont le plus éloignées possible l'une de l'autre. Lors du premier relevé, en tE, les effets de vitesse des parties mobiles ont été complètement compensés. On compense maintenant au moyen des impulsions 33 et 34 les effets de la deuxième impulsion de sélection 15. La date pivot pour le calcul des moments est maintenant la date tE. Sur l'axe de lecture X on compense également par deux impulsions 42 et 43 les effets de la deuxième moitié 60 de l'impulsion 16 et de la première moitié 61 de l'impulsion 41. La présence des impulsions compensatrices 39 et 40 sur l'axe Y se justifie par le fait que l'impulsion 17 n'est pas bipolaire. Car la compensation de son moment à la date pivot tE par les impulsions 37 et 38 n'implique pas une telle compensation a une autre date pivot telle que tE. Et ceci bien qu'aucune impulsion À 2! supplémentaire de perturbation n'intervienne sur cet axe entre temps. La figure 4 représente un tube d'écoulement 51 sensé se trouver dans la tranche examinée 7 du corps 2. Le haut de la figure représente une situation dans laquelle le tube 51 est coplanaire au plan X, Y de la coupe. Dans la partie basse de la figure la situation est différente: le tube 51 est supposé être perpendiculaire au plan X et Y de la coupe. De part et d'autre sont représentés des voxels notés 52 à 55. Ces voxels ont la forme de parallélépipèdes rectangles ayant une forme de section sensiblement carrée. La section est carrée dans la mesure o, dans l'imagerie, l'échantil- lonnage représentatif est effectué avec la même résolution selon l'axe X et l'axe Y du plan de coupe. La grande longueur du paral- lélépipède rectangle s'étend selon l'axe Z. Cette longueur est liée à l'épaisseur de la coupe 7. Cette épaisseur est un compromis expéri- mental entre la finesse de résolution des coupes et la quantité de signal global mesurable. Plus la coupe est épaisse plus le signal mesurable est élevé, et bien entendu moins l'image est précise. Dans la pratique la grande dimension est de l'ordre de 7 à 15 mm, la section carrée est de l'ordre de 1 mm par I mm. Dans le conduit 51 un fluide non représenté circule avec un diagramme de vitesse 56. Les vitesses 57 des particules situées au centre sont plus élevées que les vitesses 58 des particules situées sur les bords du conduit. Dans les artères la vitesse 57 peut être de l'ordre de 50 cm par seconde. Les voxels 52 et 53 comportent, quand ils sont situés à l'aplomb de la section 59 du conduit 51, néces- sairement des particules animées avec des vitesses très différentes. Avant l'invention l'étalement du spectre de vitesse était tellement important que la contribution des voxels 52 et 53 dans le signal global était nulle: on avait l'impression que le tube était creux. L'image des voxels 54 et 55 pouvait dans certains cas être révé- latrice. Elle l'était dans la mesure o ces voxels se situaient à l'aplomb de la partie de la section du conduit o l'étalement du spectre des vitesses était faible. C'est le cas pour le voxel 55 o les vitesses (au centre du tube) tout en étant élevées sont relativement homogènes entre elles. Leur déphasage correspondant est sensi- blement le même: la moyenne n'est pas nulle. Par contre pour le voxel 54 situé à l'aplomb du bord du conduit 51, là o l'évolution des vitesses est brutale, le signal restitué est nul. Il en résultait que l'image des tubes 51 perpendiculaires au plan de coupe était appré- ciée géométriquement plus petite qu'elle n'était en réalité. Avec l'invention l'image des écoulements est révélée fidèlement. On peut moduler l'effet de la vitesse avec le procédé de l'invention. Il suffit de ne pas annuler totalement le coefficient multiplicateur de V dans l'expansion de Af((t) o les gradients compensateurs sont présents. Dans l'invention on propose de réaliser une deuxième expérimentation dans laquelle la compensation n'est plus complète et dans laquelle en définitive la mesure est sensible à cet effet de vitesse. On compare ultérieurement les résultats obtenus au cours de ces deux expérimentations. On en déduit les caractéristiques de vitesse recherchées. D'une manière préférée on fera varier la modulation selon un axe particulier. Cet axe parti- culier peut être l'axe X, l'axe Y, l'axe Z, ou même un axe différent. Pour faire varier cette modulation on modifie l'amplitude (À) des impulsions de compensation appliquées respectivement sur les axes X, Y, Z ou sur une combinaison de ces axes. Dans la deuxième expérimentation on s'arrange pour que la sensibilité à l'effet de vitesse soit telle que, pour la vitesse nominale attendue à mesurer, l'écart maximum de phase soit de l'ordre ou inférieur à un radian. Dans ces conditions on peut montrer que l'information de vitesse est obtenue simplement en calculant le rapport de la variation relative des mesures, d'une expérience à l'autre, à la mesure de la première (ou de la deuxième) expérience. Autrement dit si S1 est le résultat de la première expérience, et si S2 est le résultat de la deuxième, l'information de vitesse moyenne est représenté par (S1 - S2)/S. Au delà de cette sensibilité d'une part le calcul proposé n'est plus applicable tel quel pour connaître l'information de vitesse, et d'autre part la sensibilité risque d'être tellement forte qu'on retomberait dans les difficultés que l'invention a résolu. C'est un calcul de ce type qu'effectuent les moyens de comparaison 8. Le calcul est entrepris pixel par pixel. D'une expérimentation à l'autre la sensibilité est passée de zero à une certaine valeur calculable par l'intégrale de G (t). et de G' ().'. On peut donc quantifier la vitesse mesurée. REVENDICATIONS 1. Procédé de modulation de l'effet de la vitesse des parties mobiles d'un corps (2) dans une mesure (figures 1-2) de densité par résonance magnétique nucléaire, mesure pour laquelle - on soumet (1) le corps à un champ magnétique orien- tateur (B0) constant pour orienter dans une direction unique (Z) des moments magnétiques de spin du corps, - on fait subir (3, 4) à ce corps une excitation magnétique radiofréquence (10) en présence (13) et/ou suivie (14, 17) de l'ap- plication d'une séquence d'un champ magnétique perturbateur, - et on relève (5) un signal radiofréquence (S) émis en réponse par le corps, caractérisé en ce qu'on module l'effet de la vitesse (V) des parties mobiles du corps créé par la séquence du champ perturbateur par l'application, avant le relevé (tE), d'une séquence (31 - 43) d'un I, champ magnétique compensateur dont l'integrale calculée sur sa durée est nulle, et dont l'historique et la valeur sont fonction de l'historique et de la valeur du champ perturbateur. 2. Procédé selon la revendication 1, caractérisé en ce que la séquence du champ magnétique perturbateur comporte des impul- sions de champ magnétique selon trois axes orthogonaux (X, Y, Z), et en ce que la séquence du champ magnétique compensateur comporte des impulsions de champ magnétique selon ces trois mêmes axes, pour moduler un à un les effets de la vitesse des parties mobiles du corps selon ces trois mêmes axes. 3. Procédé selon la revendication 2, caractérisé en ce que les impulsions de champ magnétique compensateur sont déterminées à priori en forme (ô), en durée, et en position (t1 - t2), et en ce que leur amplitude (')est évaluée pour obtenir la modulation re- cherchée. 4. Procédé selon la revendication 2 ou 3, caractérisé en ce que la séquence des champs magnétiques compensateurs comporte des couples bipolaires d'impulsions. 5. Procédé selon la revendication 2, 3 ou 4, caractérisé en ce que la séquence des champs magnétiques compensateurs comporte des couples d'impulsions, chaque impulsion d'un couple étant de valeur, de forme, de durée, et de signe égal à l'autre impulsion du couple, et ces impulsions se trouvant respectivement avant et après une deuxième impulsion haute fréquence. 6. Procédé selon l'une quelconque des revendications 2 à 5, caractérisé en ce que l'intégrale du produit de la valeur des impulsions de champ perturbateur par la durée qui les sépare du relevé du signal émis est compensée par l'intégrale d'un même produit obtenu avec les impulsions de champ compensateur. 7. Procédé selon la revendication 4 ou 5, caractérisé en ce que les impulsions des couples d'impulsion comportent des impulsions les plus éloignées possible l'une de l'autre. 8. Procédé selon l'une quelconque des revendications I à 7, caractérisé en ce que l'excitation magnétique radiofréquence com- porte une première impulsion (10) radiofréquence pour faire basculer de 90 l'orientation des moments magnétiques des spins, suivie d'au moins une deuxième (11) impulsion magnétique radiofréquence pour faire basculer de 180 l'orientation des moments magnétiques des spins, et en ce qu'on relève (tE) comme signal de résonance magnétique l'écho du signal émis en reponse après qu'il ait été réfléchi par la deuxième impulsion radiofréquence. 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que la mesure comporte la production de l'image d'un volume dont la position et les dimensions sont déterminées par la séquence (13 - 15) du champ magnétique perturbateur. 10. Procédé selon l'une quelconque des revendications I à 8, caractérisé en ce que la mesure comporte la production (6) d'une image d'une coupe (7) du corps dont la position dans le corps est déterminée par la séquence (13 - 15) du champ magnétique pertur- bateur. Il. Procédé selon la revendication 9, caractérisé en ce que la production de l'image est accomplie selon un procédé de type transformée de Fourier inverse à trois dimensions (3DFT). 12. Procédé selon la revendication 8, caractérisé en ce que la production de l'image est accomplie selon un procédé de type transformé de Fourier inverse à deux dimensions (2DFT). 13. Mise en oeuvre du procédé selon l'une quelconque des revendications I à 12, caractérisée en ce qu'on le met en oeuvre une première fois (I1) pour compenser l'effet de la vitesse, et une autre fois (I2) en modifiant une ou plusieurs des caractéristiques du champ magnétique compensateur, et en ce qu'on compare (8, I3) les me- sures obtenues dans les deux mises en oeuvres pour en déduire la vitesse des parties mobiles concernées. 14. Mise en oeuvre selon la revendication 13 et selon la reven- dication 2, caractérisée en ce que la modification porte sur des caractéristiques du champ magnétique compensateur relative à un axe pour en déduire la composante de vitesse des parties mobiles selon cet axe (figure 4). 15. Mise en oeuvre selon la revendication 13 ou la revendi- cation 14, caractérisée en ce que la modification est telle que le déphasage ( Af (t)) correspondant sur le signal émis par les parties mobiles se déplaçant à une vitesse nominale est de l'ordre de 1 radian.