FR 2630825 A1 19891103 FR 8905604 A 19890427 263082S PROCEDE ET DISPOSITIF POUR LA CORRECTION D'UNE DEGRADATION DE PHASE PROVOQUEE PAR DES COURANTS PARASITES La présente invention a trait à l'imagerie à résonance magnétique (MRI) et plus particulièrement à la réduction d'images fant8mes provoquées par des cou- rants parasites dans de telles images. Les gradients utilisés dans des systèmes d'ima- gerie à résonance magnétique (MRI) engendrent des cou- rants parasites. Les courants parasites modifient la phase de chaque point d'image d'une manière inconnue. Par conséquent, lorsqu'on mesure des différences de phase, telles que dans certaines séquences d'écoulement, il n'est pas possible de savoir de combien les différen- ces de phase sont dues à l'écoulement et de combien elles sont dues à des courants parasites. Certaines séquences de balayage utilisent des gradients d'une première polarité suivis de gradients de polarité opposée en tant que parties de la séquence pour quantifier l'écoulement par exemple. Les gradients produisent chacun des courants parasites qui ne sont pas nécessairement adaptés et par conséquent ne sont pas annulés par des procédés de soustraction. Par conséquent, dans ces séquences de balayage, les différences de phase dues à des courants parasites provoquent des mesures erronnées d'écoulement. Par conséquent, l'homme de l'art a tenté et tente encore de développer des systèmes et/ou des procé- dés pour corriger une dégradation de l'angle de phase provoquée par des courants parasites. 263082S Par conséquent, un des buts de la présente in- vention est de proposer un moyen et un appareil en vue de corriger une dégradation de l'angle de phase provoquée par des courants parasites. Selon un large aspect de la présente invention, un procédé en vue de minimiser des variations d'angle de phase provoquées par des courants parasites dans des données d'imagerie à résonance magnétique (MRI) est pro- posé, ledit procédé comportant les opérations consistant à: acquérir des données à faible résolution en- gendrées par des spins qui ont été perturbés dans un plan transversal, ladite acquisition de données se pro- duisant dans le domaine temporel en utilisant une séquen- ce de balayage multidimensionnelle, effectuer la transformation de Fourier des don- nées à faible résolution du domaine temporel acquises pour obtenir des valeurs d'amplitude et de phase pour des points d'une images, lesdites valeurs comprenant des variations de valeur de phase dûes à des courants para- sites et à un mouvement des spins, effectuer une opération sur les données du do- maine temporel acquises pour ne laisser passer sensible- ment que des données avec des variations de phase dûes à des courants parasites, effectuer une transformation de Fourier sur les données soumises à ladite opération pour obtenir des va- leurs de phase et d'amplitude pour chaque point d'image dans lesquelles les valeurs de phase sont celles provo- quées par des courants parasites, et retrancher l'angle de phase de chaque point d'image obtenu par transformation de Fourier des don- nées ayant subi ladite opération de l'angle de phase de chaque point d'image obtenu en effectuant une transfor- mation de Fourier sur les données du domaine temporel n'ayant pas subi ladite opération pour ainsi fournir une image limitée à des variations de phase provoquées par un mouvement des spins sans les changements de phase provoqués par des courants parasites. - Une autre caractéristique de la présente inven- tion est de proposer un procédé de correction de distor- tions de phase provoquées par des courants parasites et utilisant les angles de phase mesurés corrigés pour déterminer un écoulement, ledit procédé comportant des opérations consistant à: acquérir une image à faible résolution non compensée d'un fluide en écoulement, acquérir une image à faible résolution compen- sée d'écoulement dudit fluide en écoulement, obtenir une image non compensée à résolution normale dudit fluide'en écoulement, obtenir une image compensée d'écoulement à ré- solution normale dudit fluide en écoulement, lesdites images comprenant des données d'angle de phase, retrancher les données d'angle de phase des images à faible résolution des données d'angle de phase correspondantes des images à résolution normale pour produire des images compensées et non compensées à angle de phase corrigé, retrancher lesdites images compensées à angle de phase corrigé desdites images non compensées à angle de phase corrigé pour obtenir une image de l'écoulement avec une distortion d'angle de phase rendueminimale, et utiliser l'angle de phase de ladite image de l'écoulement avec une distortion d'angle de phase mini- male pour obtenir la vitesse d'écoulement pour chaque point d'image de ladite image de l'écoulement. Les caractéristiques et objets qui précèdent, ainsi que d'autres, de la présente invention seront mieux compris à la lecture de la description qui va suivre d'un aspect large d'invention faite en regard des dessins an- nexés. La Figure 1 est un schéma synoptique du système inventif. La Figure 2a est un schéma représentant une "fenêtre" ou une opération de "multiplication" du système de la Figure 1 utilisé pour séparer les images d'angle de phase provoquées par des courants parasites. La Figure 2b représente le traitement des si- gnaux reçus pour une image de résolution normale en vue d'obtenir l'image d'angle de phase totale. La Figure 2c représente la soustraction de l'image d'angle de phase de la Figure 2a de l'image d'an- gle de phase de la Figure 2b pour obtenir une image d'an- gle de phase qui est proportionnelle à la vitesse sans la distortion d'angle de phase provoquée par les courants parasites. La Figure 3 représente une séquence de balayage en vue d'acquérir des données non compensées. La Figure 4 représente une séquence de balayage en vue d'acquérir des données compensées. Le système de formation d'images à résonance magnétique 11 représenté sur la Figure 1 comprend un ai- mant 12 de grande dimension en vue de recevoir un malade pour des études de résonance magnétique. Le champ magnéti- que statique important est engendré par un générateur de champ magnétique 13, désigné par HO. Le champ magné- tique statique important, ainsi qu'il est bien connu, aligne des noyaux possédant un nombre impair de protons et/ou de neutrons; c'est-à-dire des particules nucléai- res possédant à la fois un moment angulaire de spin et un moment de dipôle magnétique. Ces particules atomiques, ayant ces deux caractéristiques sont désignées ci-après par "spins". Un champ magnétique à radio fréquence (RF) est appliqué sous la forme d'impulsions transversalement au champ magnétique statique important. Les impulsions RF appliquent de l'énergie aux spins, les perturbant ou les "basculant" dans un plan transversal, par exemple. Après l'excitation par l'impulsion de champ magnétique RF, les spins se dispersent (se déphasent) progressive- ment et retournent également en alignement avec le champ statique. En le faisant, les spins restituent de l'énergie sous la forme de signaux d'amortissement d'induction li- bre (FID) faibles mais détectables. Ces signaux FID sort détectés et utilisés par un système commandé par un dispositif de commande-calcu- lateur pour produire des images d'une manière bien con- nue. La fréquence d'excitation des impulsions RF et la fréquence FID sont définies par ce qui est connu par re- lation de Larmor qui établit que la fréquence angulaire (S0) de la précession des spins est le produit du champ magnétique B agissant sur les spins et du rapport gyro- magnétique y, une constante physique nucléaire fondamen- tale pour chaque type de noyau. Mathématiquement: f = yB/2 T. Par conséquent, en superposant un champ de gradient linéaire au champ uniforme statique, des spins dans un plan choisi peuvent être excités par choix approprié de la fréquence d'excitation RF. Les techniques de résonance magnétique sont allées de mesures de spec- troscopie à la formation d'une image et plus récemment à la formation d'une image du type qui peut à la fois représenter et quantifier un débit sanguin. La mesure d'un écoulement en utilisant des techniques de résonance magnétique a été divulguée par le docteur Hahn dans le Journal of Geophysical Research, Vol. 65, N 2, pages 776 (1960). Ainsi qu'il est bien connu, lorsqu'on utilise une séquence d'échos de spins, les spins sont tout d'abord basculés de 90 par une im- pulsion RF à 90 dans le plan transversal et ensuite après un certain déphasage sont ramenés (rephasés) par une impulsion RF à 1800. Un écho apparait à un instant (T) après l'application de l'impulsion à 180 . Le temps 1 est égal au temps entre l'impulsion à 90 et l'impul- sion à 1800. Ainsi, un écho est normalement produit à un instant 2 r après l'application de l'impulsion de 90 , si les spins n'ont aucun mouvement, c'est-à-dire si ds/dt et/ou des ordres plus élevés sont égaux à zéro. Si les spins sont en mouvement, c'est-à-dire, s'ils possè- dent leur propre vitesse, par exemple, un écho apparai- tra encore à l'instant 2t sauf que l'angle azimutal auquel les spins reconvergent dans le cadre tournant sera 263082e décalé d'une quantité l = (G V)t 2; O G est le vecteur gradient magnétique linéaire, et V la vitesse de l'écoulement. Ainsi, la différence ú entre l'angle azimutal dans une image sensible à un mouvement et une image insensible à un mouvement par exemple, permet idéalement un calcul pour obtenir la vitesse de l'objet dont l'image a été formée tel qu'un fluide corporel pour chaque point de l'image. La quantité indéterminéede la dégradation de phase due à des courants parasites cependant réduit sé- rieusement la précision du calcul. L'effet défavorable du courant parasite sur la phase est accru dans des expériences de mesure de mouve- ment utilisant des gradients de codage de mouvement du fait que pour obtenir des images sensibles à un mouvement il a été constaté qu'il était nécessaire d'utiliser des gradients importants. Voir par exemple, le document 4.516.075 qui enseigne l'utilisation de gradients de di- mensions relativement grandes. Les champs de gradient magnétique enseignés dans ce brevet, provoquent des cou- rants parasites intenses et par conséquent une distortion de phase indéterminée importante. Ainsi, une mesure de la vitesse déterminée en fonction de la phase inclut norma- lement des erreurs dûes aux courants parasites. De nombreuses tentatives ont été effectuées pour surmonter ce problème. Une telle tentative est dé- crite dans un article publié dans MR Research Note Se- lected Abstracts présenté au 6ème SMRM, (Phillips Medi- cal Systems Publication) à la page 16 (1987). L'article décrit un procédé utilisant une multiplication complexe conjointement avec des techniques itératives fournissant une correction de phase spatialement non linéaire pour corriger une dégradation de phase provoquée par des courants parasites. Le procédé est indûment complexe. Une autre technique en vue de surmonter les erreurs de phase provoquées par les gradients utilisés pour obtenir une image rendue sensible à un déplacement utilise des bobines de gradients blindées pour minimiser les effets de courants parasites, voir par exemple, l'article intitulé "Multi-echo Magnetic Resonance An- giography" par C.L. Dumoulin et al, publié dans Magnetic Resonance in Medicine, Vol. 5, pages 47-457 (1987). La technique unique décrite ici n'exige pas des bobines de gradient blindées ni une manipulation complexe telle que les techniques itératives et obtenant les produits cartésiens de polyn8mes de Tchebychev. A la place, un système d'acquisitions de données d'image à faible résolution relativement simple conjointement avec mise en ouvre de fenêtres sont utilisés ainsi qu'il sera expliqué dans la description de la Figure 1. L'emplacement des signaux FID (ou échos) est déterminé en utilisant des champs magnétiques à gradient. Pour obtenir les champs magnétiques à gradient, un géné- rateur de gradient de champ magnétique 14 produisant un gradient de champ dans la direction X, un générateur de champ magnétique 16 délivrant un gradient de champ dans la direction Y et un générateur de gradient de champ ma- gnétique 17 qui crée un gradient de champ magnétique dans la direction Z sont représentés. Dans le mode d'émission, une impulsion RF est tout d'abord émise à travers une bobine, non représentée, dans l'aimant principal 12. L'impulsion RF est prélevée du générateur RF 18 qui envoie l'impulsion à radio fré- quence à un modulateur 19 en vue d'une mise en forme. Le signal mis en forme modulé est émis par un émetteur 21 par l'intermédiaire d'un duplexeur 22. Dans le mode réception, le signal reçu (FID ou écho) est prélevé par une bobine non représentée dans l'aimant principal 12 et traverse le duplexeur 22 pour être appliqué à un récepteur 23. La bobine réceptri- ce peut être la même bobine utilisée pour émettre les impulsions à radio fréquence. Le signal reçu est envoyé à un démodulateur 24. A la fois le modulateur et le dé- modulateur sont alimentés avec la fréquence de modulation par le générateur de fréquence de modulation 20. Le si- gnal démodulé conformément à l'invention de la présente demande est envoyé à la fois à un convertisseur analogi- que numérique 26 et à un opérateur de système 27 pour laisser passer seulement une partie des données du do- maine temporel. Le résultat de l'opérateur 27 est une image à faible résolution. L'opérateur 27 dans un mode de réalisation est une "fenêtre". Sur la Figure 2a, la fenêtre est repré- sentée en tant que fonction échelon. Aux fins de la pré- sente invention, toute fenêtre qui laisse passer seule- ment la partie centrale du signal du domaine temporel reçu peut être utilisée. Une telle fenêtre laisse passer seulement les signaux de fréquence spatiale. Les signaux ayant subi une dégradation de phase provoquée par des courants parasites sont des signaux de fréquence spatiale relativement basse et par conséquent, ils traversent la fenêtre 27. La variation d'angle de phase provoquée par un déplacement de spins sont des signaux à fréquence spatiale relativement élevée qui ne traversent pas la - fenêtre 27. A la fois les signaux qui ont traversé la fe- nêtre 27 et le signal qui n'a pas été soumis à une opéra- tion par la fenêtre 27 sont traités et subissent une transformation de Fourier bidimensionnelle comme indiqué par les blocs de traitement de transformation de Fourier 29 et 31 respectivement. Les sorties des blocs effectuant la transformation de Fourier sont des données d'amplitude et de phase qui sont placées dans des mémoires 32, 33 zone par zone, chaque zone correspondant à un point de l'image. La mémoire 33, qui contient des données à fai- ble résolution provenant de la fenêtre, contient sensi- blement uniquement l'information de dégradation de phase provoquée par les courants parasites. La mémoire 32 contient des données provenant de signaux qui n'ont pas été soumis à l'opération de fenêtre et contient ainsi la valeur de la phase telle qu'affectée par à la fois le mouvement de spins et la dégradation provoquée par les courants parasites. Par conséquent, pour minimiser la dégradation de phase provoquée par les courants parasites, la mémoire 33 est retranchée de la mémoire 32 dans un amplificateur de soustraction 34. La sortie de l'amplifi- cateur de soustraction 34 est l'image à phase corrigée qui est placée dans la mémoire 36 en vue d'une visuali- sation dans l'ensemble de visualisation 37. La visualisa- tion dans l'ensemble de visualisation 37 est par consé- quent pratiquement exempte d'images fantômes dies à des variations d'angle de phase provoquées par des courants parasites. Le même processus est prévu à la fois pour des signaux reçus compensés et non compensés. Les signaux compensés tels qu'utilisés ici signifient que la phase est indépendante de la vitesse de spins de sorte que tous les spins - mobiles et stationnaires - possèdent la même phase. Des données non compensées signifient que la phase est fonction de la vitesse des spins. La Figure 2a représente graphiquement la "fe- nêtre" ou l'opération de "multiplication'" de l'opérateur 27. Le signal reçu représenté en 41 est tout d'abord sou- mis à une opération par une fenêtre pour obtenir une image à faible résolution laissant passer seulement des signaux à fréquence spatiale faible. Le signal reçu est multiplié par ou passé à travers une fenêtre telle qu'une fonction échelon 42 qui se traduira par un signal 43 comprenant seulement des fréquences spatiales faibles. Comme indi- qué précédemment d'autres formes de fenêtre c'est-à-dire gaussiennes, trapézoidales, etc, pourraient également être utilisées dans le cadre de l'invention. Le critère est le passage de fréquences spatiales faibles seule- ment. Le signal 43 représente les variations de phase provoquées par les courants parasites. Ainsi, il con- tient l'erreur de phase qui provoque les images fantô- mes dans des expériences d'écoulement et que l'homme de l'art a essayé d'éliminer depuis un temps relativement long. Le signal 43 est dans le domaine temporel de sorte qu'il subit la transformation de Fourier indiquée en 44 pour fournir des données de phase pour des zones de la matrice 46 qui correspondent à des points d'image. Le signal reçu 41 qui n'a pas traversé une fe- nêtre comme représenté sur la Figure 1 subit directement une opération de transformation de Fburier comme repré- senté en 47 pour fournir des données de phase pour une mémoire 48. Les valeurs de dcnnées de phase contiennent des variations de phase dûes à un déplacement de spins ainsi que la distortion de phase provoquée par les cou- rants parasites. Les données de la mémoire 46 qui con- tiennent seulement les variations de phase provoquées par les courants parasites sont retranchées des données de la mémoire 48 pour fournir une image exempte des variations de phase provoquées par le courant parasite. Ceci est indiqué dans la mémoire 49 qui est une mémoire des données de phase sans dégradation provoquée par les courants parasites. Le système en vue de minimiser des variations d'angle de phase provoquées par des courants parasites trouve présentement une application utile dans la mesure de débit sanguin ou d'autres fluides corporels, par exem- ple. Ainsi, si des données d'image compensées et non compensées sont acquises et soumises à une opération par le système de la Figure 1 comme représenté sur la Figure 2, alors les données d'écoulement finales ne souffriront pas d'imprécisions dues à des variations d'angle de phase provoquées par des courants parasites. La séquence de la Figure 3 est une séquence en vue d'acquérir des données non compensées et la sé- quence de la Figure 4 en vue d'acquérir des données d'écoulement compensées. D'autres séquences pourraient être utilisées à la place de la séquence particulière représentée. Par exemple, des séquences expliquées dans la demande de brevet déposée en Israël le 29 Janvier 1988, N de série 85259, au nom du présent inventeur pourraient également être utilisées à la place des séquences re- présentées sur les Figures 3 et 4. Les Figures 3 et 4 sont basées sur la séquence classique d'acho de spin. Sur la Figure 3, une impulsion à 90 51 est appliquée durant l'application d'un gradient de sélection de tranche 52. Le gradient de sélection de tranche possède la partie négative habituelle 53. L'im- pulsion 51 bascule les spins dans le plan transversal. Les spins qui sont basculés sont sélectionnés par l'im- pulsion de gradient de sélection de tranche. Ensuite une des impulsions de codage de phase des impulsions de codage de phase 54 est appliquée. Ensuite, une impulsion RF à 180 56 est appli- quée pour remettre en phase les spins déphasés dans le plan transversal. A un instant après l'application de l'impulsion RF à 180 qui est éga] au temps entre l'im- pulsion RF à 90 et l'impulsion RF à 120 , un écho appa- raît. Le signal d'écho est représenté en 57. Il se pro- duit durant l'application d'une impulsion de gradient de lecture 58. Une séquence compensée est représentée sur la Figure 4. La différence principale entre la séquence des Figures 3 et 4 est les impulsions à gradient bipo- laire 71, 72. Les impulsions à gradient bipolaire 71, 72 sont utilisées pour coder des spins en mouvement en créant une différence de phase dans des spins en déplace- ment seulement. Du fait qu'elles sont bipolaires, elles n'ont pratiquement aucun effet sur des spins stationnaires. La phase des spins en déplacement est utilisée pour obtenir la vitesse de l'équation suivante (voir la demande de brevet déposée en Israël le 29 Janvier 1988, N de série 85259, au nom du présent demandeur): t = Y [ Z0Gz(t)dt - V Gz(t)dt] t0 to o: Z0 est la position de spin à t = % (t), et V est la vitesse de spin (supposée constante). Cependant, comme indiqué ci-dessus, cet angle de phase à n'est normalement pas déterminé avec préci- sion par suite d'une dégradation de phase provoquée par les courants parasites. En utilisant le système de la Figure 1 comme représenté sur la Figure 2, la dégradation de phase due à des courants parasites peut être suppri- mée de sorte que la phase représente une fonction linéai- re vraie de la vitesse. Par conséquent, en fonctionnement, le système représenté sur la Figure 1 peut fournir une image de vitesse d'écoulement vraie à partir de la soustraction des deux images. Les deux images sont: 1) l'image de l'écoulement comprenant des variations de phase dans les points d'image provoquées par le courant parasite; et 2) une image dans laquelle sensiblement toutes les va- riations d'angle de phase provoquées par le courant para- site ont été éliminées. Ensuite, en retranchant la phase de la seconde image de celle de la première image, l'image résultante est une image d'angle de phase due à la vites- se des spins sans la dégradation provoquée par les cou- rants parasites. Bien que l'invention ait été décrite en réfé- rence à des modes de réalisation particuliers, on com- prendra que ces modes de réalisation sont seulement exem- plaires et n'ont pas pour objet de limiter la portée de l'invention. R E V E N D I C A T I ON S 1. Procédé de correction de distorsions de phase provoquées par des courants parasites dans les images à résonance magnétique, ledit procédé caractérisé par les opé- rations consistant à: acquérir des données à résolution normale d'i- mages à résonance magnétique (MRI). dans le domaine tem- porel, acquérir des données à faible résolution d'i- mages à résonance magnétique (MRI) dans le domaine tem- porel, effectuer la transformation de Fourier des don- nées MRI dans le domaine temporel à résolution normale, effectuer la transformation de Fourier des don- nées MRI à faible résolution, obtenir les images de phase dans le domaine de fréquence après les opérations de transformation de Fourier, et retrancher l'image de phase des données à fai- ble résolution de l'image de phase des données à résolu- tion normale pour fournir une-image de phase à distorsion de phase provoquée par des courants parasites corrigée. 2. Procédé de correction de distorsions de phase provoquées par des courants parasites selon la revendi- cation 1, caractérisé en ce que les opérations d'acquisition de données MRI dans le domaine temporel à résolution norma- le comportent l'utilisation d'une séquence d'acquisition de données normales pour obtenir des données MRI dans le domaine temporel originales, et dans lequel l'opération d'acquisition des données MRI à faible résolution comporte l'exécution d'une opération sur lesdites don- nées MRI à résolution normale pour obtenir des données de fréquence spatiale faible uniquement. 3. Procédé selon la reverdication 2, caractérisé en ce çue ladite exécution d'une opération sur lesdites données MRI à résolution normale comporte un filtrage en fréquen- ce spatiale desdites données MRI du domaine temporel ori- ginales. 4. Procédé selon la revendication 2, caractérisé ence que ladite exécution de l'opération sur lesdites données MRI à résolution normale comporte le passage à travers une fenêtre desdites données MRI originales avec une fonction en vue d'obtenir des données de fréquence spatiale faible uniquement. 5. Procédé de correction d'une distorsion de phase provoquée par des courants parasites et d'utilisation des angles de phase mesurés corrigés pour déterminer un écou- lement, ledit procre- caractérisé arles opérations consistant à: acquérir une image à faible résolution non com- pensée d'un liquide en mouvement, acquérir une image à faible résolution compen- sée d'écoulement dudit fluide en mouvement, obtenir une image non compensée à résolution normale dudit fluide en mouvement, obtenir une image compensée d'écoulement dudit fluide en mouvement, lesdites images comprenant des données d'angle de phase, retrancher les données d'angle de phase de l'i- mage à faible résolution des données d'angle de phase cor- respondantes de l'image à résolution normale pour fournir des images compensées et non compensées à angle de phase corrigé, retrancher lesdites images compensées et non compensées à angle de phase corrigé pour obtenir une ima- ge d'écoulement avec une distorsion d'angle de phase ren- due minimale, et utiliser l'angle de phase de ladite image d'é- coulement avec une distorsion d'angle de phase minimale pour obtenir la vitesse d'écoulement pour chaque point de ladite image de l'écoulement. 6. Procédé selon la reveiOication 5, ceractérisé en ce que la vitesse est déterminée en obtenant une solution pour V dans l'équation suivante: AL =Ty[ ZoGz(t)dt - VGz(t)dt] 00o o: Z0 est la position du spirn à l'instant t=+ (to) et V est la vitesse de spin. 7. Système de correction de distorsions de phase provoquées par des courants parasites dans les images à résonance magnétique, ledit système caractfrise par des moyens (26) en vue d'acquérir des données d'image à résonance magnétique (MRI) dans le domaine temporel à résolution normale, des moyens (27) en vue d'acquérir des données d'image à résonance magnétique (MRI) dans le domaine temporel à faible résolution, des moyens (29) pour effectuer une transformation de Fourier des données MRI dans le domaine temporel à réso- lution normale, des moyens (31) pour effectuer une transformation de Fourier des données MRI à faible résolution, des moyens (37) pour obtenir les images de phase dans le domaine de fréquence après les opérations de trans- formation de Fourier, et des moyens (34) en vue de retrancher I'image de phase des données à faible résolution de l'image de phase des don- nées à résolution normale pour fournir une image de phase à distorsion de phase provoquée par des courants parasi- tes corrigée. 8. Système de correction de distorsions de phase provoquées par des courants parasites selon la revendi- cation 7,caractérisé en ce cue les royens (26,29) envued'acqué- rir des données MRI dans le domaine temporel à résolution normale comportent des moyens pour utiliser une séquence d'acquisition de données normales pour obtenir des don- nées MRI dana le domaine temporel originales, et dans le- quel les moyens (27,28,31) en vue d'acquérir les données MRI à faible résolution comportent des moyens pour effec- tuer une opération sur lesdits données MRI originales en vue d'obtenir des données de fréquence spatiale faible uniquement. 9. Système selon la revendication S,carectgrisé en ce que lesdits moyens en vue d'effectuer une opération sur les- dites données MRI originales comportent des moyens de filtrage (27) en vue d'un filtrage de fréquence spatia- le desdites données MRI dans le domaine temporel-origi- nales. 10. Système selon la revendication 8,caractérisé en ce que lesdits moyens en vue d'effectuer une opération sur les- dites données MRI originales comportent des moyens de fenêtre pour filtrer (27) lesdites données MRI origina- les avec une fonction pour obtenir des données de fré- quence spatiale faible uniquement. 11. Système de correction d'une distorsion de pha- se provoquée par des courants parasites et d'utilisation des angles de phase mesurés corrigés pour déterminer un écoulement, ledit système zcractêrisé par: des moyens (26) en vue d'acquérir une image à faible ré- solution non compensée d'un fluide en mouvement, des moyens (27) en vue d'acquérir une image à faible ré- solution compensée d'écoulement dudit fluide en mouve- ment, des moyens pour obtenir une image non compensée à résolu- tion normale dudit fluide en écoulement, des moyens pour obtenir une image compensée d'écoulement à résolution normale dudit fluide en mouvement, lesdites images comprenant des données d'angle de phase, des moyens (34) en vue de retrancher des données d'angle de phase de l'image à faible résolution des données d'an- gle de phase correspondantes de l'image à résolution normale pour fournir des images compensées et non compen- sées à angle de phase corrigé, des moyens pour retrancher lesdites images compensées et non compensées à angle de phase corrigé pour obtenir une image de l'écoulement avec une distorsion d'angle de phase minimale, et des moyens utilisant l'angle de phase de ladite image d'écoulement à distorsion d'angle de phase minimale pour obtenir la vitesse d'écoulement pour chaque point de la- dite image d'écoulement. 12. Système selon la revendication.ll,caiactérisé en ce que la vitesse est déterminée par des moyens en vue d'obtenir u-ne solution pour V dans l'équation suivante: tt =[T ZoGz(t)dt VGz(t)dt] t t o o o Z0 est la position de spin à l'instant t =4(to) et V est la vitesse de spin.