Brevet FR2616936A1 — publié le 19881223

PROCEDE DE PRISE EN COMPTE, DANS UNE IMAGE, DES MOUVEMENTS D'UN OBJETr. La présente invention a pour objet un procédé de prise en compte, dans une image, des mouvements d'un objet. Elle trouve plus partieulièrement son application dans le domaine médical o on cherche à représenter les mouvements du coeur, des poumons ou l'écoulement du sang, ou au contraire à éliminer dans l'image les perturbations résultant d'artéfacts dits de bougé. La description de l'invention va être faite dans le cadre d'une application d'imagerie par résonnance magnétique nucléaire (RMN) sans pour autant être spécifiquement restreinte à cette application. De même elle sera étudiée dans le cadre de mouvements cycliques o elle peut être facilement mise en oeuvre. Mais son adaptation à des mouvements qui ne se produisent qu'une fois peut cependant être aussi faite. La représentation des objets en mouvement passe d'une manière classique par la saisie de clichés, clichés étant entendu au sens large, successifs régulièrement répartis dans le temps au cours du mouvement à étudier. La visualisation de ces clichés, à un rythme ralenti, permet de bien étudier le mouvement de l'objet. Les procédés de télévision modernes permettent, avec des caméras, de mettre en oeuvre cette technique dans la mesure o les caméras sont suffisamment rapides pour acquérir des images très proches les unes des autres dans le temps et dans la mesure o on visualise par la suite ces images en ralentissant le rythme de projection. Ces techniques générales malheureusement ne sont pas applicables facilement aux procédés d'imagerie par résonance magnétique nucléaire, ni de même d'ailleurs aux procédés de représentation des coupes tomodensitométriques obtenues avec des tomodensitomètres. 26 1 6936 En effet, l'acquisltion d'images de type tomographique avec de tels appareils n'est obtenue qu'après un traitement de reconstruction des images des coupes. Par exemple en RMN, les particules du corps d'un patient à étudier, situées dans la coupe, sont soumises à une séquence d'excitation au cours de l'expeérimentation et restituent à l'issue un signal qui comporte une information représentative de leur nature. Malheureusement toutes les particules de la coupe réémettent un signal en même temps. Avec les procédés de reconstruction 0 d'image connus on tourne la difficulté de cette réponse simultanée en réitérant un certain nombre de fois les expérimentations et en changeant, d'une expérimentation à une autre, les conditions d'acquisition du signal. Typiquement en imagerie par résonance magnétique nucléaire on effectue 256 expérimentations successives pendant une durée totale qui, au plus court, avoisine la demie minute. Dans le domaine des images tomodensitométriques la durée totale d'acquisition de l'ensemble des vues est de l'ordre de la seconde. Or dans un cas comme dans l'autre les cycles de fluctuation du corps humain, par exemple à l'endroit du coeur, sont du même ordre. Il en résulte que l'ensemble des séquences acquises, ou que l'ensemble des vues détectées, s'effectue pendant une durée voisine de celle de la pulsation cardiaque: à l'endroit du coeur les images sont brouillées. Pour retnédier à cet inconvénient de bougé on a, dans un premier temps, effectué une syrnchronisation. Cependant, en particulier en résonance magnétique nucléaire, les procédés d'excitation les plus rapides (dits SSFP Steady State Free Precession) qui correspondent une précession libre en régime d'équilibre des moments magnétiques des particules doivent être appliqués à un rythme propre. Ce rythme propre dépend du temps de relaxation spin-spin dit T2 des moments magnétiques des particules étudiées. Ce rythme propre s'oppose bien entendu à la synchronisation. Il n'y a en effet aucune raison pour que ce rythme propre soit un multiple du rythme cardiaque. Si on synchronise (et donc si on abandonne l'idée du mode SSFP), et qu'en plus on veuille faire une image dite en T2 ( pour montrer les temps de relaxation spin-spin), il faut alors attendre entre chaque expérimentation une durée égale à trois ou quatre fois la durée T1 (temps de relaxation spin-réseau) caractéristique des particules étudiées. Autrement dit chaque expérimentation dure de l'ordre de 2 secondes. L'acquisition des 256 séquences conduit alors à une durée, pour une seule image synchronisée pendant le cycle, d'environ 10 minutes compte tenu des pertes de temps lors des synchronisations. Or les médecins, pour pouvoir bien comprendre le fonctionnement du coeur, demandent de disposer d'environ 32 images pour un cycle cardiaque complet. Ceci conduit à une durée totale d'examen de 5 heures. Cette durée est tout à fait insupportable pour les patients. En outre, même si ces patients la supportaient, le cycle cardiaque n'est pas constant. De sorte que la synchronisation, effectuée par rapport à une date précise du cycle, des autres instants du cycle ne correspond pas toujours à un même état du coeur. L'avance ou le retard de phase par rapport à un état prédéterminé du coeur dépend de l'accélération ou du ralentissement du rythme cardiaque par rapport à un rythme cardiaque nominal estimé. Il en résulte que les images présentent des artefacts. Dans un mode d'excitation en variante, dit T2 rapide avec faible basculement, et en synchronisant, on peut aller beaucoup plus vite puisque typiquement un cycle en T2 rapide dure environ 50 ms. De sorte qu'on peut acquérir 32 phases du cycle cardiaque sur une durée sensiblement égale à deux fois ce cycle. En imagerie avec résolution normale (256 séqences par phase), ceci conduit à une durée totale de l'expérimentation de minutes. Mais l'imagerie T2 rapide ne peut pas être menée systématiquement dans une même tranche: le signal de RMN s'y amortit trop, on perd en contraste. On utilise alors une technique multicoupe: pendant des temps morts on effectue des mesures dans d'autres coupes. Pratiquement l'expérimentation 26 1 6936 dure alors quatre fois plus longtemps: 40 minutes. C'est aussi bien trop long. De plus le même problème de synchronisation que cl-dessus se présente. Lorsque la fréquence cardiaque n'est pas constante (pendant 40 minutes l), les images présentent aussi des artefacts. Une autre technique a été imaginée par Michael L. WOOD et R. Marc HENKELMAN du ONTARIO CANCER INSTITUTE AND DEPARTMENT OF MEDICAL BIOPHYSICS, UNIVERSITY OF TORONTO, TORONTO, CANADA. Elle a été publiée dans Medical Physics, Volume 13, numéro 6 de Novembre/Décembre. 1986. Elle consiste pour l'essentiel à effectuer une moyenne sur plusieurs images successives, en vue d'éliminer des artefacts de bougé. Cette technique d'imagerie, qui s'intéresse A l'élimination des effets des mouvements des parties mobiles, n'est bien entendu pas susceptible de restituer l'image mobile du mouvement que l'on cherche à rendre. Elle n'est de plus pas totalement efficace et certains parasites peuvent rester dans l'image. La solution préconisée par l'invention est une solution générale qui permet, selon les besoins, de représenter les mouvements d'un objet mobile avec un ralenti aussi finement décomposé qu'on le désire, sans entraîner une multiplication correspondante des durées d'acquisitlon. En outre, dans les cas o les parties en mouvement sont considérées comme des perturbations, le procédé de l'invention permet de donner une image des parties fixes de l'objet exemptes des composantes relatives au mouvement. Pour l'essentiel, dans l'invention, on considère qu'un mouvement d'un objet s'analyse, dans une image, comme une variation temporelle de la luminosité des éléments d'image (on dit aussi pixels) de cette image. Et on avance que cette variation temporelle de luminosité d'un élément d'image peut être remplacée par un spectre de composantes fréquentielles (d'un signal élémentaire de base lorsque le mouvement étudié est cyclique). La décomposition spectrale est de préférence une décomposition dite de Fourler. Le signal temporel de variation de luminosité de chaque élément d'image est alors remplacé pour chaque élément d'image par la transformée de Fourier de cette luminosité. Dans un cas simple, quand le phénomène étudié est cyclique, on peut considérer que les luminosités des éléments d'image comportent une composante continue, une composante fondamentale de variation cyclique, ainsi que des composantes harmoniques de variation cyclique. Le cycle étant supposé le même pour tout le corps étudié, la décomposition du mouvement du corps en série de Fourier permet, élément d'image par élément d'image, de disposer d'une série d'images représentatives ehacune respectivement, de l'image fixe, de l'image de variation au rythme fondamental, et des images de variation au rythme des harmoniques du fondamental. On peut ensuite réobtenir l'image globale. du mouvement en utilisant autant de générateurs d'images plus un que d'harmoniques d'image à prendre en compte. On combine, élément d'image par élément d'image, les participations de chacune des images harmoniques, et on visualise l'image ainsi combinée. Pour obtenir un ralenti aussi ralenti qu'on le désire il suffit de faire évoluer les générateurs d'image aussi lentement qu'on le désire. En définitive, on effectue ainsi une analyse harmonique de l'image. On remarque que l'image des parties mobiles seules peut aussi être obtenue en éliminant la composante continue d'image et que, par contre, l'image exempte d'artéfact de bougé peut être obtenue en utilisant que cette image de composante continue d'image seule. L'invention a donc pour objet un procédé de prise en compte, dans une image, des mouvements d'un objet caractérisé en ce que: - on détermine les composantes de la transformée analytique du mouvement, - on élabore des images correspondant à chacune des composantes de la transformée analytique de ce mouvement, 26 1 6936 - et on compose l'image de l'objet mobile en combinant lea images des composantes. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Celles-ci ne sont données qu'à titre indicatif et nullement limitatif de l'invention. Les figures montrent: - Figure 1: un dispositif d'imagerie par résonance magnétique nucléaire utilisable pour la mise en oeuvre du procédé selon l'invention; - Figures 2a à 2d: des diagrammes temporels de signaux de codage des expérimentions en vue de permettre la reconstruction des images; - Figure 3: la modification d'un de ces signaux de codage pour mettre en oeuvre un exemple du procédé de l'invention; - Figure 4: un dispositif utile à la composition de l'image combinée. La figure 1 représente un dispositif d'imagerie par résonance magnétique nucléaire utilisé, dans un exemple, pour la mise en oeuvre du procédé de prise en compte selon l'invention. Cet appareil comporte essentiellement des moyens, symbolisés par une bobine 1, pour produire un champ magnétique uniforme intense B0 dans une région o est placée le corps 2 d'un patient, au-dessus d'un plateau 3. Soumis- à cette aimantation, les moments magnétiques des particules du corps 2 s'orientent selon la direction du champ uniforme B0. Dans cette situation on soumet en outre le corps 2 à une excitation électromagnétique radiofréquence transmise par une antenne, par exemple à barres rayonnantes 4 à 7, alimentée au travers d'un circuit résonnant 8 par un générateur d'excitation 9. A l'issue de l'excitation, une antenne de réception 10 permet de prélever le signal de déséxcitation qui marque le retour à l'équilibre des moments magnétiques des particules du corps 2. Ce signal est reçu et est traité dans un récepteur 1l en vue de la visualisation, sur un dispositif de visualisation 12, d'une image d'une coupe du corps 26 1 6936 2. De manière à discriminer dans le signal reçu les contributions de chacun des éléments de volume de la coupe, une série d'expérimentations est entreprise au cours de chacune desquelles le codage magnétique de l'espace o est situé le corps 2 est fixé a une valeur différente par des champs magnétiques supplémentaires appliqués par des bobines dites de gradient représentées schématiquement en 13 et alimentées par I0 un générateur d'impulsions de gradient 14. L'ensemble de ces moyens fonctionne sous la conduite d'un séquenceur 15. Dans la pratique on retiendra que l'image peut être reconstruite selon un procédé préféré de reconstruction dit 2DFT dont la figure 2 donne un aperçu des phases principales. Une excitation radiofréquence 16 est appliquée en présence d'une impulsion 17 de gradient de sélection de coupe. Une impulsion radiofréquence 18 dite d'écho de spin est émise ultérieurement en présence d'une impulsion 19 de gradient de sélection de coupe pour sélectionner la même coupe et pour éliminer des défauts d'homogénéité du champ B0. Pendant la réception du signal 20 de désexcitation, une impulsion 21 d'un gradient dit de lecture est appliquée (axe X). Entre l'impulsion 16 et l'impulsion 18 un codage de phase 22 du signal de RMN est appliqué. D'une expérimentation à l'autre (au bout d'une durée Tr), ce codage de phase évolue: il atteint par exemple le niveau 23 pendant l'expérimentation suivante. Le calcul de reconstruction d'image est connu. Il comporte deux transformées de Fourier: une sur k1 dependant du gradient de lecture et une sur k2 dependant du gradient de codage de phase. Au cours d'une expérimentation de rang I le gradient de codage de phase vaut k2.1i. Le corps 2 n'est bien entendu pas immobile pendant toutes ces expérimentations, en particulier à cause des pulsations cardiaques du patient. En conséquence les images acquises sont altérées par ces mouvements. Pour connaître 26 16936 l'influence de ces mouvements, dans l'invention, on cherche à connaître la transformée analytique du mouvement. Dans un premier temps, pour simplifier l'explication, cette transformée analytique peut être considérée comme étant la transformée de Fourier. L'analyse dépend en fait de l'espace des fonctions utilisées: Içi les fonctions de Fourier. Dans ce but on dispose à proximité des parties en mouvement, en particulier Ici près du coeur, un détecteur de mouvements cardiaques 24. Ce détecteur 24, du même type que ceux utilisés en synchronisation cardiaque classique, est relié à un récepteur 25 qui produit un signal 26 représentatif des mouvements de ce coeur. Le récepteur 25 est en relation avec un analyseur 27 susceptible de donner en sortie le cycle de base du mouvement. Du fait que le mouvement cardiaque est cyclique, l'analyse de ce mouvement se résume à la mesure de la durée de ce cycle. Les harmoniques d'image selon l'invention sont ensuite calculées selon une méthode décrite ci-dessous. Considérons un moment magnétique M (y,t) d'une particule dont l'amplitude est fonction du lieu y de la particule et du temps t de la manière suivante M (y,t) = a0 (y) ejA sin wt I Dans cette expression, a0 indique l'amplitude de la valeur de ce moment magnétique en fonction du temps, j est tel que j = -1, et w représente la pulsation d'un signal sinusoïdal (ici il s'agit donc d'un signal sinusoïdal pur). L'expression de M peut encore s'écrire: M(y,t) = aO(y) + aO(y) J -ejwt - aO(Y)jA e-jwt ou encore en première approximation M(y,t) = aO(y) + a1(y) ewt + a 1(y) ejwt Au cours d'une excitation de RMN, les moments 26 1 6936 magnétiques des particules situées à une ordonnée y ne contribuent pas seuls au signal de RMN reçu par l'antenne 10. Ils y contribuent en commun avec ceux des particles voisines situées à une même abscisse x mais à une autre ordonnée y, ainsi qu'avec ceux situés à d'autres abscisses x et même à d'autres altitudes z (voir les axes de référence sur la figure 1). Si on suppose que les codages provoqués au cours des séquences par le générateur 14 d'impulsions de gradient ont eu pour effet de cantonner la naissance d'un signal de RMN utilisable à une tranche prédéterminée du corps 2 (on fait une image à deux dimensions avec une sélection de coupe), et si un gradient de lecture 21 a été appliqué au moment de la détection du signal par l'antenne 10, on peut, après filtrage fréquentiel ou après calcul d'une transformé de Fourier, extraire du signal reçu un signal qui peut s'écrire de la manière suivante s(x,i.tr) = e.k2i'y M(y,i.Tr) dy Cette expression signifie que le signal détecté est 2 prélevé au cours d'une séquence d'excitation-mesure de rang i. Comnme les séquences ont chacune une durée Tr, le signal est donc détecté à la date i.Tr. Le signal de RMN à cet instant est alors l'intégrale (dy) du moment magnétique de toutes les particules se trouvant le long de Y. Pour chaque particule se trouvant dans un élément de volume d'ordonnée y, une fois qu'a été codé le signal de résonance par un gradient de codage de phase d'amplitude k2.i (puisqu'il est appliqué au cours de la séquenee de rang i), la phase relative vaut k2. i. y. Dans cette expression on peut remplacer M par sa valeur fonction de a0, al, et a_1. Si de plus on considère, dans un cas particulier qu'ultérieurement on généralisera, qu'un signal de RMN a été émis en Y0 seulement, le signal s devient (en négligeant maintenant l'influence de x qui de ce point de vue ne modifie pas le sens des expressions) 26 1 6936 s(1. Tr) = a0(yO) ejk2'l'Yo + a1(Y()eJk2 i(Y0-wTr/k2)+ a On peut alors considérer que la première partie (en a0) du deuxième membre de l'égalité concerne la contribution des particules fixes à l'endroit de l'élément de volume y0, tandis que ce qui suit (en a1 et a1) concerne les particules mobiles dans l'élément de volume yO. Supposons maintenant qu'avec la même valeur de k2.1 on effectue un nombre NE d'expérimentations successives Identiques. En définitive au lieu d'acquérir une série de 2N séquences (i allant de -N à N-1; dans un exemple type o 2N = 256,1 va de -128 à +127) on peut acquérir NE (par exemple NE =8 pour fixer les idées) séquences identiques pour chaque valeur de 1. On acquiert donc NE x 2N séquences au total (voir Fig. 3). Les NE séquences sont identiques à cela près que les instants d'acquisition ne sont pas les mêmes, et sont bien entendu décalés dans le temps. Supposons, en outre, qu'au lieu d'un signal sinusoïdal pur on ait affaire en fait à un signal périodique de pulsation fondamental w mais possédant en plus certains harmoniques. Si on appelle p(t) le signal représentant la phase de la pulsation fondamentale de variation du signal périodique, on peut écrire: M(Yo,l) a0(yO) - a, (y.) e jp(l) + a_ (yo) e -jp(l) + a2 (yo) e j2p(l)+ a_2() e -j2p() I + a3 (y0) e j3p(l)+ a 3(y0) e -i3p(l) + _3 Dans cette expression 1 désigne le rang d'une des expérimentations prises parmi les NE retenues. Il indique en fait la date 1. Tr A laquelle est mesuré le signal de RMN. Dans 26 1 6936 1! la pratique 1 vaut de O à NE-i. Cette dernière expression contient en substance toutes l'essence de l'invention. M (y0,1) représente la luminosité à afficher d'un point d'ordonnée y0 à une date 1. La valeur a0(yO) représente la composante continue de luminosité à cet endroit de l'image: elle est indépendante du temps puisqu'elle est une composante continue. Les termes en a1 et a 1 représentent, en fonction du fondamental p(1) et du temps 1, l'influence des coefficients pondérateurs de ce fondamental à l'endroit de l'ordonnée y0. Les termes en a2 et a_2 représentent les mêmes influences pour l'harmonique 2 (2p(1)), et ainsi de suite. Autrement dit, pour un lieu y0, connaître la valeur du moment magnétique M (ou ce qui est équivalent la luminosité d'un point correspondant sur l'image) revient à déterminer les coefficients a0 et an et an, jusqu'à n aussi grand qu'on le désire, de manière à bien représenter le mouvement. On sait que la décomposition en série de Fourier permet, pour représenter -la modulation du signal de départ, de se satisfaire d'un nombre restreint d'harmoniques: ceux dont les coefficients an ont encore une -valeur significative. Dans l'explication qui va suivre on se limitera à n = 3 bien que cela ne constitue pas une obligation. Il y a lieu cependant de se rappeler l'application du théorème de Shannon selon lequel la connaissance d'un signal à une fréquence nw n'est atteinte que si pendant une durée 1/nw le nombre d'échantillons pris en compte est supérieure à 2. On peut montrer que dans le cas présent, o on a choisi NE = 8, la connaissance du troisième harmonique est donc totalement envisageable (8 est plus grand que 2 x 3+1). Si l'on effectue, ainsi qu'on l'a prévu, les NE = 8 mesures, on est alors confronté à la situation selon laquelle on dispose de 8 valeurs du signal s(l.Tr) alors qu'on ne cherche en fait que 7 inconnues qui sont a0, a1, a 1,... a_3. La solution de ce système d'équations, peut être apportée par une méthode d'approximation, sur le plan pratique par une méthode d'approximation par des moindres carrés. 12 2616936 Dans un cas particulièrement simple, ou les NE mesures sont effectuées pendant une durée exactement égale à la période fondamentale du phénomène cyclique examiné, et dans lequel la pulsation du fondamental est constante dans le temps, on peut arriver après calcul à la solution simple suivante 1/7 a0 s.(l.Tr) a = 1/8> 1=0 1 =7 -,. 21T.1.1 a1 1_/8 A e-J8.s(l.Tr) l 0 17C.1.3 a3 =l/84 e _8À. s(l.Tr) 1=0 Les deux dernières expressions nous permettent de montrer qu'en fait ah est le coefficient pondérateur de la composante de Fourler de rang h de la décomposition en série de Fourier du phénomène cyclique en question. Mals dans ce cas o la fréquence n'est pas constante, ou bien dans le cas o les NE mesures ne sont pas reparties exactement sur un cycle du mouvement, les fonctions ejp(l) ne sont plus orthogonales. Il faut généraliser le traitement du calcul des ah sous forme, par exemple, d'un calcul des moindres carrés. Ce calcul est rendu nécessaire du fait que, dans ces cas, la décomposition de la fonction temporelle d'aimantation n'est plus effectuée sur une base de fonctions analytiques orthogonales. Autrement dit la base des fonctions analytiques choisies n'a pas besoin d'être orthogonale sous reserve qu'une approximation des ah soit entreprise. La dernière formulation II de M(yo, 1) conduit, en assimilant M à s dans la mesure o seul un élément de volume d'ordonnée yO du corps 2 émet une réponse & l'issue des expérimentations, à la formulation matricielle suivante: I e jp(O) e2jp(O) ...........................e3j()aS) I e jP(0) e3jp(O) e23P !.ej....... s(i) ejp(7)........................................... e3jp(7) a 3 s(7) e a3 s(7) Dans cette equation le signe = a un sens qui doit donc être Interpreté dans le sens d'un calcul par des moindres carrés. Cette équation peut encore s'écrire IPI I"a"-IS Cette formulation matricielle rend compte des valeurs des différents signaux s (LTr) en fonction des harmoniques (np(1)) du mouvement considéré. En définitive ce qui est inconnu est le vecteur "a", tandis que ce qui est connu est le vecteur S ainsi que la matrice P. La matrice P est connue puisque on connaît le cycle fondamental 1.p(l), les harmoniques h.p (1), et puisqu'on connaît également les instants 1.Tr auxquels sont pris les échantillons. La théorie du calcul matriciel permet d'écrire pour mettre en oeuvre le calcul par les moindres carrés: Dans cette expression P* reprsente la matrice conugue de la Dans cette expression P* représente la matrice conjuguée de la 14 2 6 1 6 9 3 6 matrice P, telle que chaque élément qli de P* est tel que qJ = Pji* ai PJi est un élément de P. L'exposant -1 correspond à une inversion de matrice dont la définition est telle que PP I o I est la matrice Identlté. Avec cette dernière méthode on peut donc calculer le vecteur "a".Ce calcul est le calcul des moindres carrés. Il importe maintenant de généraliser l'explication au cas o il n'y a pas seulement un point y0 qui émet un signal utile. En effet dans une expérimentation d'imagerie par résonance magnétique nucléalre, dans une coupe au moment de la lecture, toutes les partleules situées sur une même colonne (donc d'ordonnée y différente et de même abscisse x) contribuent au signal. On peut écrire, en supposant que les a n sont connus, qu'en un endroit d'ordonnée yet à une date 1; le i5 moment magnétique correspondant a une amplitude de la forme h= +n M (y,l) = = ah (y).e jhp(l) II1 h -n déjà vue précédemment mals généralisée cli à toutes les coordonnées et à toutes les harmoniques. Comme tous les signaux situés à une même abcisse x et d'ordonnée y participent en même temps au signal reçu, on peut écrire S(=) e k.2i.y ah (y) e jhp(l) ly - n h qui est une expression Intégrale de laquelle on peut sortir l'expression suivante: h (+n /y jhp(l) s(l) = f eJk2iy.ah (y) dy e h -n 26 1 6936 Dans cette dernière expression, on peut admettre que l'intégrale puisse s'écrire sous la forme Ah. Cette formulation en fait repose sur un certain nombre de suppositions. Premièrement il est supposé que, pendant les NEx2N mesures effectuées les coefficients ah sont constants. Ceci signifie en définitive que le phénomène est stationnaire (mals pas que son cycle est de durée constante) . Si son rythme n'est pas constant, par contre sa forme du mouvement doit pouvoir se déduire, d'un cycle à l'autre, par compression ou dilatation temporelle uniquement. Dans cette expression i représente le codage i. k2 tandis que h dépend de l'harmonique considéré. Le vecteur précédent "a", à sept composantes pour trois harmoniques, devient maintenant une matrice généralisée appelée A, à (2n+1).2N dimensions. Dans l'exemple, elle a 7 x 256 dimensions. Le chiffre 256 représente la résolution désirée de l'image ainsi que le nombre désiré de différentes valeurs de k2.i. Cette matrice peut s'écrire 0 I A256 A AO..................... 0 i A A A256 A256 A0 A256 A 3: i 256 A................................................ A0 1 A 256 _3A13--.................... A 3 26 1 6936 La matrice A correspond polir toutes les séquences 1, aux contributions de toutes les particules d'une même colonne à l'image fixe (A0), aux images du fondamental (A1, A 1), ou aux images des harmoniques (An, A n). Pour chaque séquence, pour chaque valeur de i.k2, compte tenu du fait qu'on a acquis 1 signaux, il est possible de calculer les coefficients Ah par moindres carrés. Autrement dit, A chaque série de NE mesures pour une valeur i. k2 déterminé, il est possible de connaître les coefficients A' en question. Mais ce qui nous intéresse est en fait la connaissance des coefficients ah (y). Il importe donc d'effectuer des transformées de Fourier, en utilisant le fait qu'on dispose de multiples expérimentations iL.k2 (et dans ce cas on dit qu'on fait une transformée de Fourier "sur k2") pour retrouver les coefficients ah (y). Ces transformées de Fourier sur k2 sont identiques à celle pratiquée dans un procédé de reconstruction d'image du type 2DFT classique. La seule particularité de l'invention réside dans le fait que ces dernières transformées de Fourier sont effectuées, ici, 2n+1 fois étant donné qu'on cherche à connaître, pour chaque élément d'image de l'image, non pas une seule valeur de luminosité, mais les 2n+1 valeurs a0, a1, a 1,... a n des harmoniques d'image. La figure 4 montre comment, ultérieurement, est recombinée l'image deséléments mobiles. On dispose en pratique de 2n+1 mémoires d'image numérotées 101 à 104. Ces mémoires d'image sont destinées a être lues simultanément sous la commande d'un dispositif de remise à zéro 105. Leurs contenus sont additionnés, élément d'image, par élément d'image, dans un additionneur 106 en vue d'alimenter une autre mémoire d'image 107 qui, elle, peut être lue par le dispositif de visualisation 12. Dans l'invention on multiplie les valeurs, élément d'image par élément d'image, des mémoires d'images par une fonction cos h wt + J sin h.wt avant l'addition. Des multiplicateurs 108 à 110 sont ainsi montrés à cet effet. La valeur t n'est pas nécessairement prise au rythme o elle évolue en temps réel. 26 1 6936 Dans le cas o on veut un ralenti trés finement représenté on fait évoluer très lentement la valeur t tout en lisant la mémoire d'image 107 à un rythme compatible avec une lecture vidéo classique (25 images par seconde). En résumé dans l'exemple décrit on n'a acquis que huit fois plus d'informations (NE = 8) et on peut avoir un ralenti aussi finement représenté qu'on le désire (en faisant faiblement varier la valeur de _t). Mais l'invention apporte en plus un autre avantage. Si on considère que tout ce qui bouge n'est pas intéressant, on peut l'éliminer en ne prenant en compte que les coefficients de la matrice a0 contenus dans la mémoire d'image 101. De cette manière on élimine simplement les artéfacts de bougé. Il reste maintenant à traiter plusieurs problèmes. Le premier concerne la variation du cycle de base w du mouvement étudié au cours de toutes Ies mesures. En effet si la durée de celle-ci dure de l'ordre de cinq minutes, le rythme cardiaque ne peut pas être considéré comme constant pendant une telle durée. Cependant, avec le détecteur 16, le récepteur 17, et l'analyseur 18, il est possible de connaître non seulement le cycle cardiaque mais même les variations de ce cycle. On peut admettre dans un premier temps que les variations de ce cycle seront faibles de sorte que le nombre d'échantillons à prendre en compte, eu égard au nombre d'harmoniques qu'on cherche à représenter, respectera toujours les contraintes opposées par Shannon. SI tel n'est pas le cas il suffit tout simplement d'augmenter d'une quantité nécessaire le nombre des échantillons. On se retrouve alors dans la situation suivante: les instants 1, o les mesures sont effectuées, sont connus dans l'absolu puisqu'ils sont imposés par le séquenceur. Mais leur position relative par rapport au cycle fondamental évolue. Cela a seulement pour effet de compliquer légèrement le calcul des e]hp(l) puisque, dans ces conditions, la matrice P varie d'une séquence à l'autre. La phase p (1) peut s'écrire p (l) =l Tr (t) dt Dans l'analyseur 27 en fait on ne mesure pas la forme du mouvement si ce mouvement est cyclique: on se contente de mesurer le temps d'apparition d'un événement caractéristique du cycle mesuré. Par définition la phase de la pulsation fondamentale lors de l'apparition de cet évenement est 251 (en fait k.2 rythme, la diastole s'accélère relativement plus que la systole. Ainsl, dans ce cas particulier, la détermination des composants de la transformée analytique (de Fourier) du mouvement se résume à la mesure de l'arrivée d'un évenement caractéristique et à la mesure de la phase relative par interpolation (linéaire ou non). On peut aussi étudier le cas o le mouvement n'est pas seulement créé par un seul mouvement périodique mais par deux mouvements périodiques imbriqués. Dans la pratique on pense ainsi à la simultanéité de la pulsation cardiaque et du rythme respiratoire. Il n'y a aucune raison pour que l'un soit le multiple de l'autre. En pratique tout ce qui a été dit jusqu'à présent peut être généralisé puisque la combinaison de deux mouvements peut s'analyser comme une modulation par multiplication d'un mouvement par l'autre. Le spectre de la décomposition en série de Fourier qu'il y a lieu de prendre en compte est alors équivalent au spectre de la pulsation la moins rapide déplacée des valeurs du fondamental et des harmoniques du spectre de la pulsation la plus rapide. Sur la figure 3 on a représenté. des séries d'expérimentations avec des mêmes Valeurs de k2.1 qui se suivent consécutivement. Il est cependant possible d'interlacer des valeurs de k2 1 égales avec d'autres valeurs de k2. Cette hypothèse se justifie d'autant plus que le procédé de l'invention est compatible avec un mode d'excitation SSFP dont le rythme propre peut être très faible (bientôt quelques millisecondes). Dans un cycle cardiaque d'une seconde il devient alors possible d'effectuer plusieurs trains de sequences à des valeurs i.k2 différentes mals interlacées entre elles: par exemple un train comportant des suites, une séquence i.k2 une séquence (i + 1). k2, suivi d'un train comportant des suites une séquence (i + 2). k2 - une séquence (i + 3).k2, etc.... La décomposition en série de Fourier du mouvement est par ailleurs uniquement justifiée dans la mesure o, dans les multiplicateurs 108 à 110, on utilise, ce qui paraît le plus simple, des fonctions sinusoïdales. On peut cependant envisager de décomposer le mouvement selon d'autres fonctions analytiques dont l'intérêt pourrait se justifier pour telle ou telle raison. Pour cette raison, dans la revendication principale il est évoqué une transformée analytique, cette transformée analytique étant une transformée de Fourier, lorsque que l'analyse du mouvement est une analyse de Fourier. En résumé, dans le cadre d'une application d'imagerie par RMN le procédé de l'invention comporte les phases suivantes. Pour chaque expérimentation on excite le corps, on code magnétiquement l'espace, et on entreprend la lecture du signal réémis par Ies particules. La présence d'un gradient de lecture pendant cette lecture permet, en effectuant une transformée de Fourier sur l'axe de lecture (en général l'axe x) de disposer, pour NE mesures avec un gradient de codage de phase i.k2, et pour toutes les valeurs possibles de ce gradient codeur de phase, d'un ensemble de valeurs de signaux détectés fonction de leur date d'acquisition (1) et du codage qu'ils ont subi (i.k2). Pendant la mesure des signaux réémis, on effectue également une mesure du rythme du mouvement des parties en 2 6 1 6 9 3 6 mouvement (par exemple du rythme cardiaque). On détermine alors par rapport à ee rythme mesuré, des évolutions p (1) utiles en tenant compte, en particulier, des variations possibles de ce rythme au cours de l'ensemble des mesures. On dispose donc de la matrice P utilisable pour chaque série de NE séquences. On calcule alors IP*.PI 1l.P* pour une valeur de i.k2, c'est-à-dire pour NE valeurs de signal mesurées. On calcule ainsi les A" relatifs à la séquence i.k2, dans un h 2' exemple pour h appartenant à -3,+3. Et on réitère ce caleul autant de fois qu'il y a de valeurs i.k2 différentes. Ensuite, en disposant de l'ensemble des valeurs A on calcule par 2n + 1 transformées de Fourier sur k2 l'ensemble des valeurs ah (y) affecté à chacune des ordonnées y des points d'une colonne d'abcisse x. On effectue le même travail pour toutes les abcisses x du plan, et on dispose ensuite des mémoires d'images représentées sur la figure 4. La recombinaison de l'image s'effectue ensuite en fonction, de ce que l'on cherche à montrer. Comme l'élaboration des images est faite par caleul en nombres complexes (j) on peut décider, d'une manière connue, de montrer l'image réelle, l'image imaginaire, ou même l'image en module. On peut même si les besoins de la visualisation le requièrent moduler par un coefficient de contraste certaines des images, fixes ou mobiles, de façon à souligner ou à effacer certaines parties de l'image recombinée. 26 1 69D6 REVENDICATIONS 1. Procédé de prise en compte, dans une image, des mouvements dans un objet caractérisé en ce que: - on détermine les composantes (24-27) de la transformée analytique de ces mouvements, - on élabore (Fig. 3) des Images correspondant à chacune de ces composantes de la transformée analytique de ces mouvements, - et on compose l'image de l'objet en combinant les images des composantes (Fig. 4). 2. Procédé selon la revendication 1 caractérisé en ce que l'analyse est une analyse de Fourier, les transformées analytiques étant des transformées de Fourier. 3. Procédé selon la revendication 2, caractérisé en ce que pour élaborer les composantes de la transformée de Fourier on décompose le mouvement en série de Fourier. 4. Procédé selon la revendication 3 caractérisé en ce que pour élaborer les images des composantes on pratique une inversion de matrice généralisée pour mettre en oeuvre un algorithme de moindres carrés. 5. Procédé selon la revendication 3 ou 4, caractérisé en ce que pour décomposer on mesure la période d'un mouvement cyclique de base de l'objet. 6. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que pour composer on combine certaines ou même une seule, des images des composantes. 7. Procédé selon l'une quelconque des revendications 1 à 6 caractérisé en ce que pour combiner les images des composantes on les affecte de coefficients de contraste. 8. Utilisation du procédé selon l'une quelconque des revendications 1 à 7 dans un procédé d'imagerie par résonance magnétique nuléare. magnéitque nucléaire.