EP0208601B1 - Procédé d'imagerie par résonance magnétique nucléaire - Google Patents

Procédé d'imagerie par résonance magnétique nucléaire Download PDF

Info

Publication number
EP0208601B1
EP0208601B1 EP86401423A EP86401423A EP0208601B1 EP 0208601 B1 EP0208601 B1 EP 0208601B1 EP 86401423 A EP86401423 A EP 86401423A EP 86401423 A EP86401423 A EP 86401423A EP 0208601 B1 EP0208601 B1 EP 0208601B1
Authority
EP
European Patent Office
Prior art keywords
sequences
image
axis
excitation
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
EP86401423A
Other languages
German (de)
English (en)
Other versions
EP0208601A1 (fr
Inventor
Eric C/O Thomson Csf Scpi Breton
Denis C/O Thomson Csf Scpi Le Bihan
Patrick C/O Thomson Csf Scpi Le Roux
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Electric CGR SA
Original Assignee
General Electric CGR SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from FR8509824A external-priority patent/FR2584188B1/fr
Priority claimed from FR8512352A external-priority patent/FR2586296B1/fr
Application filed by General Electric CGR SA filed Critical General Electric CGR SA
Publication of EP0208601A1 publication Critical patent/EP0208601A1/fr
Application granted granted Critical
Publication of EP0208601B1 publication Critical patent/EP0208601B1/fr
Expired legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion imaging

Definitions

  • the present invention relates to a nuclear magnetic resonance imaging method.
  • the use of this method is particularly envisaged in the medical field to represent sections of organs of the human body.
  • Nuclear magnetic resonance imaging is mainly developed as a means of medical diagnosis. It allows visualization of internal tissue structures with a contrast and resolution never achieved simultaneously by other imaging methods.
  • To obtain a nuclear magnetic resonance image of an organ section with differentiation of the tissue characteristics of this organ we use the property that certain atomic nuclei, such as protons, have of orienting their magnetic moment by acquiring energy when placed in a constant main magnetic field B o .
  • a particular zone of an object containing nuclei then presents a global magnetic moment which can be tilted according to a given orientation, perpendicular or parallel to the field Bo, by inducing a resonance by emission of a radiofrequency field perpendicular to the main field .
  • This signal can be received by a receiving antenna.
  • the duration of return to equilibrium of the global magnetic moment of a region considered and the decrease of the signal depend on two important factors: the spin-network interaction and the spinspin interaction of the nucleus with the surrounding matter. These two factors lead to the definition of two relaxation times called T i and T z respectively .
  • a region considered of an object therefore emits a signal whose intensity depends on T i , ⁇ 2 , the proton density of the region, and the time that has elapsed since the radio frequency excitation.
  • phase dispersion As the displacements of the molecules are in all directions, the resulting phase dispersion has the effect that certain contributions come to oppose each other. The signal received is then weaker. This loss of sensitivity somehow measures the diffusion characteristic of a medium: a strongly diffusing medium sees the intensity of its magnetic resonance signal decrease very rapidly with the echo time used.
  • the object of the present invention is therefore to propose images where the parameter highlighted is precisely the characteristic of molecular diffusion in the tissues studied in order to improve their differentiation.
  • the present invention makes it possible to establish a diffusion image while avoiding these drawbacks.
  • the images are of constant thickness, which makes acquisitions on humans possible, even with a multi-cut process.
  • the sensitivity to diffusion is good: it is due to the use of fairly long echo times and to gradients effective in their intensity and position.
  • the exact determination of the diffusion coefficients is obtained without standard substance.
  • the absolute effect of the diffusion could be calculated from the acquisition parameters.
  • FIG. 1 represents a device for implementing the imaging method according to the invention.
  • This device comprises means symbolized by a coil 1 for subjecting a body 2 to a large and constant magnetic field B o .
  • This device also comprises generator means 3 and coils 4 for subjecting the body thus placed to spin echo sequences to one or more echoes in the presence of sequences of field gradients (FIG. 2 a).
  • the coils 4 represent the radio frequency coils and the field gradient coils.
  • It also comprises reception means 5 connected to the coils 4 for receiving the magnetic resonance signal, and means 6 for calculating and memorizing a first image 1 1 and a second image 1 2 relating to two series of experiments imposed by commands symbolized Ci and C 2 of the generating means 3.
  • comparison means 7 we compare re point to point the images 1 1 , 1 2 , by calculating the logarithm of the ratio of the values representative of the magnetic resonance signals. The means 7 then develop a third image 1 3 in which two regions 8 and 9 from the middle to the place of the imaged cut have responses in different diffusion value whereas they would have presented identical responses in conventional magnetic resonance image . These images can be viewed on a viewing device 50.
  • the process for calculating the images l 1 or 1 2 is conventional.
  • the imaging method used by the means 6 is a so-called 2DFT method.
  • This imaging method currently provides the best image quality.
  • only one section plane is excited at a time, by means of radio frequency excitations (90 ° or 180 ° ) of particular shape, in the presence of a so-called selection gradient.
  • the selection gradient can be oriented along the Z axis to select a transverse section, that is to say along an X, Y plane.
  • the principle of 2DFT imaging is to phase code the different signals acquired.
  • phase shifter of axis perpendicular to a reading gradient whose direction is constant.
  • the reading gradient could be of gradient X and the phase shift gradient could be gradient Y.
  • a double spatial Fourier transform the image is constructed, hence the name of the method.
  • a description of this mode of imagery can be found in the book: "Magnetic resonance imaging” by M. LE BIHAN published by Editions MASSON, Paris, MARCH 1985. An improvement of this method can moreover make it possible to obtain simultaneously images of several parallel sections.
  • FIG. 2 a typically represents sequences of field gradient along the three axes Z, X, Y as well as the instants of application of radio frequency excitations tending to tilt the spins of the nuclei by 90 ° for small excitations and by 180 ° for great excitement.
  • To image a section of a body it is necessary to carry out sequences of excitations with spin echo in the presence of sequences of field gradient, in as large a number of times as the resolution of the expected image must be more precise.
  • the phase shift gradient, Y varies in successive steps from a certain value to the same value but of opposite sign. This value depends on the shape and the duration of the reading gradient.
  • This phase shift gradient makes it possible to rotate each spin by a variable phase, dependent on its ordinate along the Y axis and on the value of this gradient.
  • the gradient Y can successively take the same number of values: in a preferred manner the definition of the two images will be the same.
  • the number of radiofrequency excitation sequences with spin echo and / or the intensity and the shape of the gradient sequences is linked to the way in which the images have been acquired and to the information content which these images contain.
  • N excitations 12 said to be 180 ° succeed a single excitation 11 at 90 ° .
  • the number N is greater than or equal to one.
  • the second sequences of spin echo excitations there is only one 180 ° pulse which follows the 90 ° excitation.
  • the duration of the echo time T E of the second sequence of spin echo excitations is equal to N times the duration of the echo time T E ' of the first sequence of spin echo excitations.
  • the successive echo times are equal to each other. It is easier, in fact, to divide a given duration T E into an integer of equal elementary durations T E ' . But if the successive echo times T E ' of the first spin echo sequences are not equal to each other, the invention can still work: it is important that the duration between the 90 ° 11 pulse and the measurement date 10 in both cases is the same.
  • the first sequence with spin echoes comprises four equal echo durations T E ′ in the middle of which each time a 180 ° pulse takes place.
  • T E for the same duration T E , a single pulse at 180 ° 26 is interposed.
  • the duration T E is worth 112 milliseconds
  • the duration T E is worth 28 milliseconds.
  • the repetition time TR which is the time which separates for one or the other image, each sequence with echo of spin during which the gradient Y takes a different value, is worth one second. This value of 112 milliseconds is chosen to be large enough to obtain good sensitivity to the effect of the diffusion and a reasonable signal-to-noise ratio. It is known by commands C 1 and C 2 to cause such sequences of excitation-measures. These commands depend on the machines used.
  • p is the density of protons at the point considered
  • y is the gyromagnetic ratio of the molecules of the medium in this place
  • G ij and dij correspond to a pair of compensated homologous gradients as presented in FIGS. 2b or 2c .
  • the sums relate to the number m of homologous pairs compensated for gradient pulses and to the number N of spin echoes.
  • N the compensated homologous gradient pairs corresponding to a normal imaging sequence.
  • NT E the second sequence of excitations
  • NT E ' is replaced by T E
  • N is replaced by 1.
  • Gij, dij, and my are then different.
  • the signal is said to be S 1 (T E ).
  • the contributions of a point of the medium are brought in the magnetic resonance signal in common with the contributions of the other points of the section examined in the medium.
  • the 2DFT imaging method ultimately makes it possible to assign to each point of the images 1 1 or 1 2 a value which is representative of these contributions.
  • FIG. 3a represents by the envelope curve 13 the theoretical signal resulting from the existence of a relaxation time T 2 when all the other parasitic effects, including the diffusion, can be neglected.
  • the 3a comprise an excitation tending to cause the nuclei of the body to tilt 90 ° (excitation 15) and four successive excitations tending to cause these nuclei to tilt 180 ° (16 to 19).
  • the dotted signal 20 is the nuclear magnetic resonance signal after the first 90 ° excitation.
  • the signal S N of the preceding formula corresponds to the amplitude of the peak 24.
  • FIG. 3a also shows the envelope 25 of these peaks when the body is subjected to diffusing sequences.
  • Curve 25 is also reported in Figure 3b.
  • a diffusing sequence a single excitation at 180 ° 26 follows the pulse at 90 ° 15.
  • the signal measured at the end of the echo time T E is rated S 1 .
  • the value of this signal is given by a formula identical to the previous one: the only difference comes from the fact that in this case N is worth 1.
  • N is worth 1.
  • the curve 14 approaches the curve 13. This makes it possible to represent an image where the effect 28 of the diffusion alone can be even better highlighted. It is useful to note that the third image obtained with the method of the invention differs from the teaching given in the aforementioned document. With this teaching, the signal is proportional to a difference 29 between the effects of diffusion in the body studied and in a standard body whose response in diffusion would be written in the dashed curve 30.
  • the ratio of the signals S N and S 1 is carried out and the logarithm of this ratio is taken.
  • FIG. 2a An additional gradient is shown in Figure 2a below the graph of radio frequency excitations. It can be applied to any of the X, Y or Z axes. It can for example be applied to the phase coding or selection axis. Applied on the axis of the selection, it does not interfere in the thickness of the selected slice in the body insofar as its application, additional, is made outside the times of application of the radio frequency excitations during which the gradient is present. Indeed, if the selection gradient receives additional value during the 90 ° selection and 180 ° reversal excitations, the inhomogeneity it causes will be greater. For the same given passband of the excitation, the thickness of the slice excited in the body 2 will be thinner. The signal that this thinner slice will restore will be weaker.
  • the additional gradient will be applied in addition to the reading gradient.
  • it will be arranged so that the gradients added are added in periods as distant as possible from those in which the tilting excitations, the excitations 12, are applied. -.
  • FIG. 3b gives an idea of the durations during which the field gradient supplements 41 to 48 can be applied in a preferred manner. These are applied outside the dates on which the peaks of the nuclear magnetic resonance signal are noted.
  • the gradients 41 and 48 are dyssimetric. But as they are homologous compensated with respect to the tilting pulse 26 their integrals with respect to time are equal.
  • the large pulses 41 and 48 can be used alone. They are as far as possible from the excitation 26.
  • the choice of the axis on which the additional pulses of the gradients are applied may not be indifferent. If we consider that molecular diffusion is an isotropic phenomenon this is indeed the case. On the other hand, if we consider that molecular diffusion can, in certain cases, be anisotropic, it may be interesting to choose the axis on which we will favor the diffusion effect.
  • certain tissues of the human body have a privileged orientation by their arrangement in the body. This privileged orientation results from a privileged form of the cells which constitute them. These cells which have no reason to be globally symmetrical on three axes then have different molecular diffusion coefficients according to each of these three axes.
  • the measurement of the diffusion given in the third image is a quantitative measurement starting from the calculation of the coefficients b km (k being worth 1 or N according to the sequences) and therefore of factor K.
  • the imperfections inherent in the imaging systems can involve an overall attenuation of the signals S N or Si at each point of the image, for example slightly different cutting thickness, poor rephasing of the cutting selection of one of the two images, ...
  • This overall attenuation can be calibrated to l using a standard substance 33 placed along the body studied so that it appears on the edges of the image field.
  • Se N and Se i are the magnetic resonance signals in the standard body, corresponding to the signals S N and Si in the body studied, respectively at the end of the weakly diffusing and diffusing sequences, we can write that the mollecular diffusion coefficient d '' a point of the cut is worth: where De is the known molecular diffusion coefficient of the standard body. Under these conditions the calibration is obtained simply.
  • the invention is implemented by modulating the effect of the speed of the moving parts (the blood) of the tissue according to a modulation method, the description of which is appended hereafter.
  • the modulation method described below makes it possible to modulate the effect (on the resonance signal) of these microcirculations.
  • this modulation process it is known, by adding so-called compensating bipolar gradients, to ensure that in the image the molecules of the blood are taken into account as if they were fixed. Under these conditions the diffusion phenomenon appears alone.
  • This modulation process was previously described in a French patent application filed on AUGUST 13, 1985 under M * 85 12352. The content of this prior application is now an integral part of the present invention.
  • the micro-circulations (they are slow displacements) disturb the diffusing sequences because of the sufficiently strong sensitization gradients that these sequences contain. It is therefore preferably the diffusing sequences which are compensated.
  • the inventors of the present invention believe that this non-neutralization is due to a temporally inconsistent distribution of molecular diffusion displacements.
  • micro-circulations are coherent displacements since their speeds are almost constant. In other words, the reduction of the resonance signal measured during a diffusing sequence (compared to that measured during a non-diffusing sequence) is now due only to the diffusion and no longer also to the micro-circulations whose effect has been neutralized.
  • the object of the improvement is a method of modulating the effect of the speed of the moving parts of a body in a density measurement by nuclear magnetic resonance (NMR) as well as the implementation of the method to deduce the speed of the parts. affected mobiles.
  • NMR nuclear magnetic resonance
  • This improvement finds its application particularly in the medical field where the bodies examined are human bodies and where the moving parts are blood cells circulating in veins or arteries or mobile organs such as the heart muscle. In this application, it can even more particularly be implemented with an imaging method to give an image representative of the distribution of the speeds of the moving parts in a section of the body examined.
  • the resonance signal emitted vibrates at a frequency fo which depends on the intensity of the orienting magnetic field Bo and on the gyromagnetic ratio y characteristic of the medium examined.
  • the magnetic moment of spin of an immobile particle undergoes in this case only a zero global phase shift. Indeed, the phase shift suffered during the application of a first part of the bipolar gradient is compensated by the application of the second part of this bipolar gradient.
  • a mobile particle having a positive speed along the axis of the gradient, undergoes during the second part of the pulse a phase shift greater in absolute value than during the first part. Because it frequents during this second part a region of space where, because of the gradient the disturbing magnetic field is stronger.
  • the sensitivity of the speed phenomenon to the disturbing magnetic field applied is such that, ultimately, the detection of displacement phenomena n It is possible that as long as the maximum speeds are below a limit.
  • the sensitivity of NMR machines is currently around 1 radian / (cm / s) at 0, 2 rd / (cm / s). This means that a particle which moves 1 cm per second in the plane of the section contributes to the global signal emitted with a phase shift of 1 radian compared to the contributions emitted by the fixed particles.
  • the object of improvement is to highlight the effect of the speed of the moving parts of a body by modifying the sensitivity of machines in a particular way.
  • the modification of sensitivity while not modifying in any way the signals emitted by the fixed parts, can have the effect of canceling the part of phase shift due to speed.
  • the moving particles then contribute to the overall signal as if they were fixed.
  • the modulation process also allows, rather than canceling everything, to modulate the effect of speed.
  • the principle stated in the modulation process is very general. It is not confined to a more complete application of imagery. It can in particular be implemented with mass spectrometers. Furthermore, it applies regardless of the radio frequency excitation procedures selected: procedures often known by their Anglo-Saxon names of saturation recovery, inversion recovery, double saturation recovery, saturation recovery - inversion recovery, spin-echo, etc. .. Likewise, the invention, within the framework of imaging, is applicable whatever the imaging procedure chosen: for example that it is a rear projection method of the PC Lauterbur type, a method by Fourier transform of the 3DFT or 2DFT type developed by A.
  • the improvement also relates to an implementation of this method, characterized in that it is used for the first time to compensate for the effect of speed, and another time for implementing it by modifying one or more of the characteristics of the magnetic field. compensator and in that the measurements obtained in the two implementations are compared in order to deduce therefrom the speed of the moving parts concerned.
  • FIG. 4 represents a machine usable for implementing the improvement process. It includes means symbolized by a coil 1 for subjecting a body 2 to a large and constant magnetic field B o .
  • Field B o is the orienting field.
  • the machine also includes generator means 3 and coils 4, connected to the generator means 3, for subjecting the body thus placed to radio frequency excitation sequences in the presence of disturbing field sequences: field gradients oriented along three reference axes of machine X, Y or Z ( Figure 5).
  • the coils 4 represent the radio frequency coils and the field gradient coils at the same time.
  • the machine also includes receiving means 5, connected to the coils 4, for receiving the magnetic resonance signal.
  • means 6 can be used to calculate and store a first image 1 1 and a second image 1 2 of a section 70 of the body 2.
  • the two images relate to two series of experiments imposed by symbolized commands C 1 and C 2 of the generating means 30.
  • comparison means 80 the images 1 1 , 1 2 are compared point to point, and a third image l 3 is produced . This image is can be viewed on a viewing device 50.
  • an image point, or pixel p conventionally represents by its luminosity the density of particles contained in an element of volume, or voxel y which corresponds to it in the section 70.
  • the voxel v comprises numerous animated particles of velocities V comprising components V x , Vy and V z on each of the three reference axes of the machine.
  • the invention will make it possible to calculate in the voxel v, and for all the other voxels also, the average speed and the number of the particles of this voxel thus in motion.
  • Image 1 3 is representative of this speed information, by comparison between two images 1 1 and 1 2 where the effect of this speed will have been modulated in two different ways.
  • FIG. 5 shows for this purpose radiofrequency pulses 100, 110 and 120 applied in the presence of a so-called selection field gradient, for example oriented along the Z axis, and represented respectively by pulses 130, 140 and 150.
  • the selection gradient is oriented along the Z axis the section is transverse, that is to say along an X, Y plane.
  • the different acquired signals are phase coded. This is obtained by a pulse, of variable intensity of a gradient called phase shifter with an axis perpendicular to a so-called reading gradient whose direction is constant.
  • the reading gradient can be on the X axis and the phase shift gradient can be on the Y axis.
  • the X gradient is rated 160
  • the Y gradient is rated 170
  • An improvement of this method can moreover make it possible to obtain images of several parallel couples simultaneously.
  • precoding 170 which takes some time, postpones the measurement of the signal emitted by the body. As this signal fades very quickly, we imagined to measure an echo. This echo is produced, at the time of application of the radiofrequency pulses 110 and then 120 by a reflection of the phase dispersion of the contributions made by each of the particles. For this reason these pulses are drawn in FIG. 5 higher than the pulse 100 since they cause the orientation of the magnetic moments of the spins of the particles to be tilted by 180 ° whereas the pulse 100 had only tilted them by 90 ° .
  • phase shift gradient Y varies in successive steps starting from a certain value up to the same value but of opposite sign. This value depends on the shape and the duration of the reading gradient 160.
  • This phase shift gradient makes it possible to rotate each magnetic spin moment of a variable phase, depending on its ordinate along the Y axis and on the value of this gradient.
  • the gradient Y can successively take the same number of values: in a preferred manner the definition of the two images will be the same.
  • pulses 200 and 210 respectively to these axes.
  • the time integral of these pulses is the same but of opposite sign, they neutralize the parasitic coding effect.
  • the integral is understood in the sense of integral in time.
  • the hatched areas 200 and 190 are therefore equal to the dotted areas 180 and 210 respectively.
  • the pulses 200 and 210 do not interfere in the cut selection since at the time when they are applied the pulse 100 is no longer present.
  • the pulses 140 and 150 are self-neutralized since they comprise portions 220 and 230 respectively, and 240 and 250 which are antisymmetrical with respect to the middle of the radiofrequency pulse, 110 and 120 for which they reselect the cut 70 in the body 2.
  • the type 2 DFT imaging method is the most used in practice because of its speed of acquisition (compared to 3D methods) and its robustness against imperfections of the physical system, in particular vis- vis-à-vis the inhomogeneities of the orienting field. This is particularly appreciated in comparison with rear projection methods or with methods not using spin echo. In return, this method is the most intrinsically sensitive to the movement of particles.
  • the immunity to inhomogeneities of the B o field of the 2 DFT method stems from the fact that, unlike rear projection methods, the frequency encoder gradient is in the same direction (the reading gradient X) and of the same amplitude, G L , of a sequence to another. The inhomogeneity of the Bo field then only distorts the isofrequency lines of resonance.
  • the inhomogeneity of the field becomes the equivalent of a bad linearity of the reading gradient.
  • the images are consequently distorted but without the creation of blurring and without the associated loss of resolution which one meets with the methods by rear projection.
  • the reading gradient changes orientation with each sequence and therefore it distributes the inhomogeneities of the orienting field throughout the image for which it thus creates a blur.
  • Figures 6a and 6b show the phase shift ⁇ resulting in the signal contribution emitted by two particles of the same voxel subjected to identical disturbances, when one of the particles (solid line) is stationary while the other particle (dashes) ) is driven by a speed V.
  • Figure 6a shows the time diagram of an RF radio frequency signal comprising a pulse 260 to create a spin echo. According to an axis A a magnetic field gradient was imposed by taking before and after the pulse 260, for example, a value G for a duration 8 each time. In a given place, a previously excited, immobile particle vibrates at a frequency fo at the start of the application of pulse 270. During the duration e of this pulse 270 it vibrates at a different frequency.
  • this frequency is higher, so that at the end of the pulse 270 it can be considered, since it resumes vibrating at the frequency f o that it has undergone a phase shift ⁇ 1.
  • the radio frequency pulse 260 it can be assumed that everything happens as if the phase of the signal emitted by this particle were reversed: - ⁇ 1 .
  • the phase shift phenomenon continues: the contribution is still dephased by ⁇ 1 . But this last phase shift erases the effect of the first after it has been reversed by pulse 260. Consequently, the motionless particle finds it at the end of pulse 280 the phase it had at the start of the pulse 270.
  • phase shift ⁇ 1 is proportional on the one hand to the amplitude of the gradient G and on the other hand to the duration during which each pulse of the gradient is imposed. More generally ç1 is proportional to the integral of this gradient during the duration of existence of its pulse if this pulse is not exactly square.
  • the pulses 270 and 280 of the bipolar gradient must be replaced by pulses 290, 300 of opposite signs.
  • the phase for the fixed particles evolves to ⁇ 3 at the end of the pulse 290 and evolves in opposite direction to zero during the inverse pulse 300.
  • the phase shift ⁇ between the contribution of the moving particles and those fixed particles is of the same order as before.
  • G (t) therefore represents all the sequences required by the measurement. What characterizes the improvement is that these sequences are maintained as they are but that their effects on the speed are compensated for by the addition of additional sequences of field gradients G 'called compensators such as where G 'represents gradients, also bipolar, i.e. such that
  • the micro-circulation phenomenon is particularly present in angiomas and tumors. Hence the interest in removing its effects to correctly measure the diffusion.
  • the intensity of the micro-circulations can be quantified by comparison of a diffusion image not compensated for in speed with a diffusion image compensated in speed. Indeed the first image is an image of diffusion and micro-circulation while the second is only an image of pure diffusion. The subtraction of the two images, if their characteristics are equal, leads to an image of micro-circulation.
  • each axis X, Y, Z the addition of such compensating bipolar gradients in the preferred form of pairs of pulses denoted 310 and 320, 330 and 340, 350 and 360, 370 and 380, 390 and 400 and 420 and 430.
  • the incidence of these compensating impulses is symbolized in the drawing by the presence of small crosses. It is known, in the machines of the state of the art, to impose such pulses by modifying the commands C 1 or C 2 of the generating means. The time integral of the pulses of each couple is zero.
  • the pulse 310 is applied during the same period. In this case the "new" pulse 200 is stronger.
  • the pulses 350 and 370 are also applied during this period.
  • the pulses 320, 360 and 380 are preferably applied between the end t 3 of the pulse 140 and the start t 4 of the read pulse 160.
  • the pulses 330, 390 and 420 are applied between the end ts of pulse 160 and the start t 6 of pulse 150.
  • the pulses 340, 400 and 430 are applied between the end t 7 of pulse 150 and the start t 8 of a pulse 410 which is used for reading a second echo.
  • the measurement of the signal emitted is carried out between dates t ' 4 and t's contained between the dates t 4 and t s and also during the duration of the pulse 410.
  • the compensation pulses applied on an axis , X, Y or Z are intended to erase the effect of the disturbing gradients applied on the same axis.
  • the integral expression implicitly gives the gradients G '; it is possible to make simplifications.
  • G ( ⁇ ) . ⁇ are known terms (they belong to the 2DFT sequence): their integral is also. This integral can therefore be taken out of expression.
  • the forms of the compensating pulses 310 to 430 will be known forms and, as far as possible, identical to forms already practiced for the normal pulses of the disturbing sequences. Given the dates used for the application of these pulses, a unit integral of each of them can be calculated and the expression of G 'can be written in another way:
  • the first member is known, it concerns disturbing gradients.
  • the second member has integrals of gradients of preferred shape r and of unit amplitude. We know how to calculate the integral of the second member since we know the chosen form r. Consequently, we can deduce ⁇ , that is to say the amplitude of the compensating pulses.
  • the effect on the phase of the pulse 190 is of course opposite to that of the pulse 210 since they are located on either side of the pulse 110. Consequently their common effect on the phase is analyzed as a moment oriented in the direction of the pulse having the largest lever arm: pulse 210.
  • This moment should be countered by a couple of pulses 350, 360 whose total integral is zero (it does not causes no phase shift on the signal of the fixed parts) but whose global moment is opposed to the global moment of the two pulses 210 and 190.
  • the pulse 350 being placed in time on a date earlier than the pulse 360 by compared to the date t E1 its effect is preponderant. It is indeed of opposite direction to the impulse 210.
  • pulses 180 and 200, and 220 and 230 it is each time the pulses 180 and 220 which take precedence: their lever arm is longer. Consequently, it is well understood that the pulse 310 must be in the same direction as the pulse 200.
  • the case of the pulses 370 and 380 is a little different. In fact, they must counteract the effects of a gradient which varies from one acquisition to another. The amplitudes of the pulses 370 and 380 will therefore also vary in correspondence from one acquisition to another (symbolized by arrows). For the above reasons when the pulse 170 is positive the pulse 370 must be negative.
  • FIG. 7 represents a flow tube 510 supposed to be in the examined section 70 of the body 2.
  • the top of the figure represents a situation in which the tube 510 is coplanar with the plane X, Y of the section. In the lower part of the figure the situation is different: the tube 510 is assumed to be perpendicular to the plane X and Y of the section.
  • On both sides are represented voxels denoted 520 to 550. These voxels have the shape of rectangular parallelepipeds having a shape of substantially square section.
  • the section is square insofar as, in the imagery, the representative sampling is carried out with the same resolution along the X axis and the Y axis of the section plane.
  • the great length of the rectangular parallelepiped extends along the axis Z. This length is linked to the thickness of the section 70. This thickness is an experimental compromise between the resolution resolution of the sections and the quantity of overall measurable signal. The thicker the cut, the higher the measurable signal, and of course the less precise the image. In practice, the large dimension is of the order of 7 to 15 mm, the square section is of the order of 1 mm by 1 mm.
  • a fluid not shown circulates with a speed diagram 560.
  • the speeds 570 of the particles located in the center are higher than the speeds 580 of the particles located on the edges of the conduit.
  • the speed 570 can be of the order of 50 cm per second.
  • Voxels 520 and 530 include, when they are located directly above section 590 of conduit 510, necessarily animated particles with very different speeds.
  • the spreading of the speed spectrum was so important that the contribution of the voxels 520 and 530 in the overall signal was zero: one had the impression that the tube was hollow.
  • the image of voxels 540 and 550 could in some cases be revealing.
  • the inventors have in fact realized, from experience, that the compensation of the speed of the blood in the capillaries linked to the micro-circulation, in order to show only the phenomenon of diffusion alone, was not always effective because the velocities of micro-circulation in the tissue were not constant in value and in direction during the measurement time (echo time).
  • the constancy in direction is related to the geometry of the capillary network in the studied tissue: for example the direction of the microcirculation in a muscle can be constant (rectilinear vessels) over a distance greater than the distance traveled by the blood during the measurement time.
  • the micro-circulation (we also say sometimes the perfusion) is oriented in all directions.
  • the tangle of capillary networks is even such that, despite the speed of NMR measurements (an excitation-measurement sequence can typically last 100 ms), the orientations of the molecules' displacement speeds change direction during the measurement: which can be interpreted as a macro-diffusion on the scale of a voxel.
  • the representation in the image of the diffusion parameter D of the tissue is distorted: the coefficient represented contains a component linked to the value of the true diffusion (D) for the static part of the tissue and a component (D ') linked to the presence of microcirculation in a given voxel.
  • S (p, Ti, T 2 ) and S '(p, Ti, T 2 ) represent the NMR signal respectively for the tissue (the brain tissue) and for the fluid of the active capillaries.
  • the function F equal to 1 for sequences with little scattering, has a value less than 1 in the case of scattering sequences. It depends on the gradients, the speed of the blood in the microcirculation, and the geometry of the capillaries.
  • F is a cardinal sinus of the gradient and the speed of microcirculation.
  • Micro-circulation is assimilated to macro-diffusion.
  • This expression is an equation with two unknowns, D and f, of which we can find a solution by repeating the experiments so as to measure two values of Dm for two different values of b k m. From the system of two equations with two unknowns obtained we extract D and f. The two values of b k m rete in addition, they comply with the indications given above relating to the need to make fF 1 negligible.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Vascular Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

  • La présente invention a pour objet un procédé d'imagerie par résonance magnétique nucléaire. L'utilisation de ce procédé est particulièrement envisagée dans le domaine médical pour représenter des coupes d'organes du corps humain.
  • 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'organe avec différenciation des caractéristiques tissulaires de cet organe, on utilise la propriété qu'ont certains noyaux atomiques, 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 Bo. Une zone particulière d'un objet contenant des noyaux 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 Bo, en induisant une résonance par émission d'un champ radiofréquence perpendiculaire au champ principal.
  • Tous les noyaux possédant alors un moment magnétique, tournant à une vitesse de précession dite de Larmor, tendent à retrouver l'orientation initiale parallèle à Bo en émettant un signal radio fréquence à la fréquence de résonance caractéristique de Bo et du noyau. 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'interaction spinspin du noyau avec la matière environnante. Ces deux facteurs conduisent à la définition de deux temps de relaxation dits respectivement Ti et Tz. Une région considérée d'un objet émet donc un signal dont l'intensité dépend de Ti, de Ï2, de la densité protonique de la région, et du temps qui s'est écoulé depuis l'excitation radiofréquence.
  • Pour localiser une région de l'organe 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 de l'organe. Pour cela on superpose au champ principal Bo des gradients de champs magnétiques pulsés. Ces gradients sont orientés selon des directions X, Y, Z pour définir à tous moments 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. 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 noyaux par les impulsions du champ radiofréquence, ainsi que les instants de lecture ou d'acquisition des données de l'image.
  • Un autre facteur intervient pour modifier l'intensité du signal capté lorsque les noyaux retournent à leur orientation d'équilibre. Cet autre facteur dépend de la diffusion moléculaire du milieu. La diffusion moléculaire concerne les déplacements que subissent les molécules d'un milieu en fonction du temps. L'in- homogénéité du champ magnétique dans lequel se trouve ces molécules a alors pour conséquence que la fréquence de résonance magnétique de ces molécules évolue. En effet, cette fréquence est liée par le rapport gyromagnétique de ces molécules à l'intensité de ce champ. Aussi, au cours d'une expérimentation de résonance magnétique, en particulier une expérimentation avec une séquence d'échos de spin, le signal capté est inférieur à l'intensité attendue. En effet, la fréquentation par les molécules de régions de l'espace où les champs magnétiques sont différents a eu pour conséquence de modifier les phases relatives des contributions apportées par chacune de ces molécules au signal de résonance magnétique global capté. Comme les déplacements des molécules sont dans tous les sens, la dispersion de phase qui en résulte a pour effet que certaines contributions viennent mutuellement s'opposer. Le signal capté est alors moins fort. Cette perte de sensibilité mesure en quelque sorte la caractéristique de diffusion d'un milieu : un milieu fortement diffusant voit l'intensité de son signal de résonance magnétique décroître très rapidement avec le temps d'écho utilisé.
  • Il existe dans les organes humains des tissus pathologiques par exemple les angiomes et les tumeurs, qui présentent des signaux de résonance magnétique nucléaire classiques souvent identiques. Autrement dit des images de ces organes mettant en exergue les temps de relaxation Ti ou T2 ne permettent pas de discriminer ces conformations. Il en résulte que l'examen de l'image classique ne permet pas de prendre une décision thérapeutique. La présente invention a alors pour but de proposer des images où le paramètre mis en exergue est justement la caractéristique de diffusion moléculaire dans les tissus étudiés afin d'améliorer leur différenciation.
  • On connaît par le livre Biomedical Magnetic Résonance, édité par Radiology Research and Education Foundation (San Francisco) en 1984, un article dû à Messieurs WESBEY George et al, et intitulé "Translational Molecular Self-Diffusion in Magnetic Resonance Imaging : Effects and Applications". Cet article suggère de mesurer la constante de diffusion des régions d'un milieu en comparant l'effet relatif de la diffusion sur le milieu étudié et sur un corps étalon lors de séquences d'excitations magnétiques différentes. Dans la description qu'il donne de cette technique, il apparaît un inconvénient. En effet, ces séquences sont obtenues en incrémentant l'intensité d'un gradient de sélection de coupe, ce qui modifie l'épaisseur de la coupe étudiée. Cette méthode ne s'applique alors qu'à des objets plus fins que l'épaisseur de coupe la plus fine obtenue par les séquences utilisées. Elle n'est donc pas applicable chez l'homme. De plus, la sensibilité de cette méthode à la diffusion étant relativement faible (temps d'écho faibles, gradients inéfficacement placés dans la séquence), les auteurs utilisent plusieurs acquisitions pour atteindre une précision raisonnable sur la mesure. D'autre part, le recours à un même étalon est indispensable pour servir de référence aux mesures.
  • La présente invention permet d'établir une image de diffusion en évitant ces inconvénients. En particulier, les images sont à coupe d'épaisseur constante, ce qui rend les acquisitions sur l'homme réalisables, même avec un procédé multicoupe. La sensibilité à la diffusion est bonne : elle est due à l'utilisation de temps d'écho assez long et à des gradients efficaces par leur intensité et leur position. De plus, la détermination exacte des coefficients de diffusion est obtenue sans substance étalon. Dans l'invention, l'effet absolu de la diffusion a pu être calculé à partir des paramètres d'acquisition.
  • L'invention a pour objet un procédé d'imagerie par résonance magnétique nucléaire caractérisé en ce qu'il comporte les étapes suivantes pour donner une image de la diffusion moléculaire d'un corps étudié :
    • - on place le corps dans un champ magnétique constant Bo ;
    • - on soumet le corps ainsi placé à une première pluralité de premières séquences d'excitations à écho de spin en présence de premières séquences de gradients de champ ; ces premières séquences à écho de spin comportant un nombre entier N supérieur ou égal à 1 d'excitations où les moments magnétiques des noyaux du corps sont basculés de 180° après une excitation où ces moments ont été basculés de 90° pour faire ainsi des séquences à N échos peu diffusantes ;
    • - on relève les signaux de résonance magnétique en fin de ces premières séquences et on calcule une première image correspondant à l'écho N de ces signaux en attribuant à chaque point de l'image une valeur correspondant au signal de résonance magnétique du point qui lui correspond dans le corps ;
    • - on soumet le corps ainsi placé à une deuxième pluralité de deuxièmes séquences d'excitations à échos de spin en présence de deuxièmes séquences de gradients de champ ; ces deuxièmes séquences à échos de spin comportant une excitation à 180° après une excitation à 90° pour faire des séquences à un écho diffusantes, la durée d'écho des deuxièmes séquences d'excitations étant égale à la durée totale d'écho des premières séquences d'excitations-;
    • - on relève les signaux de résonance magnétique en fin de ces deuxièmes séquences d'excitations et on calcule une deuxième image correspondant aux échos de ces signaux en attribuant à chaque point de l'image une valeur correspondant au signal de résonance magnétique du point qui lui correspond dans le corps ;
    • - on compare point à point les valeurs attribuées pour la première image aux valeurs attribuées pour la deuxième image, pour élaborer une troisième image représentative de la diffusion moléculaire en chaque point du corps
  • L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Sur ces figures les mêmes repères désignent les mêmes éléments. Elles représentent :.
    • - figure 1 un dispositif apte à la mise en oeuvre du procédé d'imagerie de l'invention ;
    • - figure 2a des diagrammes temporels de séquences d'excitations et de gradients de champ appliqués au corps pour élaborer les premières et deuxièmes images ;
    • - figures 2b et 2c des diagrammes temporels de gradients dits homologues compensés, utilisés pour augmenter la sensibilité à l'effet de la diffusion ;
    • - figures 3a et 3b les allures des signaux relevés respectivement à l'issue des premières et deuxièmes séquences d'excitations ;
    • - figure 4 une machine pour la mise en oeuvre dans l'invention d'un procédé de modulation selon un perfectionnement ;
    • - figure 5 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 6a et 6b 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 7 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 1 représente un dispositif pour la mise en oeuvre du procédé d'imagerie selon l'invention. Ce dispositif comporte des moyens symbolisés par une bobine 1 pour soumettre un corps 2 à un champ magnétique important et constant Bo. Ce dispositif comporte encore des moyens générateur 3 et bobines 4 pour soumettre le corps ainsi placé à des séquences d'écho de spin à un ou plusieurs échos en présence de séquences de gradients de champ (figure 2 a). Les bobines 4 représentent les bobines radiofréquence et les bobines de gradient de champ. Il comporte encore des moyens de réception 5 reliés aux bobines 4 pour recevoir le signal de résonance magnétique, et des moyens 6 pour calculer et mémoriser une première image 11 et une deuxième image 12 relatives à deux séries d'expérimentations imposées par des commandes symbolisées Ci et C2 des moyens générateurs 3. Dans des moyens de comparaison 7 on compare point à point les images 11, 12, en calculant le logarithme du rapport des valeurs représentatives des signaux de résonance magnétique. Les moyens 7 élaborent alors une troisième image 13 dans laquelle deux régions 8 et 9 du milieu à l'endroit de la coupe imagée présentent des réponses en valeur de diffusion différentes alors qu'elles auraient présenté des réponses identiques en image de résonance magnétique classique. Ces images peuvent être visualisées sur un dispositif de visualisation 50.
  • Le processus de calcul des images l1 ou 12 est classique. Dans un exemple la méthode d'imagerie employée par les moyens 6 est une méthode dite 2DFT. Cette méthode d'imagerie permet actuellement d'obtenir la meilleure qualité d'image. Dans cette méthode d'imagerie seul un plan de coupe est excité à la fois, par le moyen d'excitations radiofréquence (90° ou 180°) de forme particulière, en présence d'un gradient dit de sélection. Dans la figure 1, le gradient de sélection peut être orienté selon l'axe Z pour sélectionner une coupe transverse, c'est-à-dire selon un plan X, Y. Le principe de l'imagerie 2DFT est de coder 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 de lecture dont la direction est constante. Par exemple, pour une coupe transverse, le gradient de lecture pourra être de gradient X et le gradient de déphasage pourra être le gradient Y. Puis par une double transformée de Fourier spatiale, l'image est construite, d'où le nom de la méthode. On pourra trouver une description de ce mode d'imagerie dans le livre : "Imagerie par résonance magnétique" de M. LE BIHAN publié par les Editions MASSON, Paris, MARS 1985. Une amélioration de cette méthode peut d'ailleurs permettre d'obtenir simultanément des images de plusieurs coupes parallèles.
  • La figure 2 a représente typiquement des séquences de gradient de champ selon les trois axes Z, X, Y ainsi que les instants d'application des excitations radiofréquence tendant à faire basculer les spins des noyaux de 90° pour les petites excitations et de 180° pour les grandes excitations. Pour imager une coupe d'un corps il y a lieu de réaliser des séquences d'excitations à écho de spin en présence de séquences de gradient 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'excitations 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 10 de lecture. Ce gradient déphaseur permet de faire tourner chaque spin d'une phase variable, dépendante de son ordonnée suivant l'axe Y et de la valeur de ce gradient. Pour chaque image h 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.
  • Ce qui diffère dans l'invention entre la première image 11 et la deuxième image 12 est le nombre des séquences d'excitations radiofréquence à écho de spin et/ou l'intensité et la forme des séquences de gradient. Le calcul effectué par les moyens 7 de comparaison est lié à la manière dont les images ont été acquises et au contenu d'information que ces images renferment. Dans les premières séquences d'excitations à échos de spin, N excitations 12 dites à 180° succèdent à une seule excitation 11 à 90°. Le nombre N est plus grand ou égal à un. Dans les deuxièmes séquences d'excitations à écho de spin il n'y a qu'une seule impulsion à 180° qui suit l'excitation à 90°. Par ailleurs la durée du temps d'écho TE de la deuxième séquence d'excitations à échos de spin est égale à N fois la durée du temps d'écho TE' de la première séquence d'excitations à échos de spin. Pour des raisons de simplicité dans les premières séquences à échos de spin, les temps d'écho successifs sont égaux entre eux. Il est plus facile en effet de diviser une durée donnée TE en un nombre entier de durées élementaires TE' égales. Mais si les temps d'écho successifs TE' des premières séquences à écho de spin ne sont pas égaux entre eux, l'invention peut quand même fonctionner : il importe que la durée entre l'impulsion 90° 11 et la date de mesure 10 dans l'un et l'autre cas soit la même.
  • Dans l'exemple qui apparaît sur la figure 3a la première séquence à échos de spin comporte quatre durées d'écho égales TE' au milieu desquelles chaque fois une impulsion de 180° prend place. Sur la figure 3b, pendant une même durée TE, une seule impulsion à 180° 26 vient s'interposer. Dans un exemple la durée TE vaut 112 millisecondes, et la durée TE, vaut 28 millisecondes. Dans cet exemple le temps de répétition TR, qui est le temps qui sépare pour l'une ou l'autre image, chaque séquence à écho de spin au cours de laquelle le gradient Y prend une valeur différente, vaut une seconde. Cette valeur de 112 millisecondes est choisie suffisamment grande pour obtenir une bonne sensibilité sur l'effet de la diffusion et un rapport signal sur bruit raisonnable. On sait par des commandes C1 et C2 provoquer de telles séquences d'excitations-mesures. Ces commandes dépendent des machines utilisées.
  • Dans l'invention, on a trouvé que la contribution d'une partie du milieu au signal de résonance magnétique relevé au bout du temps d'écho N.TE' dans une première séquence d'excitations a la forme suivante :
    Figure imgb0001
    où
    Figure imgb0002
  • Dans ces expressions p est la densité de protons à l'endroit considéré, y est le rapport gyromagnétique des molécules du milieu dans cet endroit, et Gij et dij correspondent à une paire de gradients homologues compensés tels que présentés sur les figures 2b ou 2c. Une paire de gradients homologues compensés correspond à une paire de gradients d'intensité Gij et G'ij et de durées dij et d'ij tels que Gij.dij = G'ij.d'ij. Si les gradients sont placés de part et d'autre d'une impulsion d'excitation radiofréquence à 180°; les intensités Gij et G'ij sont de même signe. Si les gradients sont placés d'un même côté, elles sont de signe contraire. Les sommes portent sur le nombre m de paires homologues compensées d'impulsions de gradient et sur le nombre N d'échos de spin. Par exemple, sur la figure 3b, il y a quatre paires homologues compensées d'impulsions de gradient auxquelles pour calculer m il faut rajouter les paires de gradient homologues compensés correspondant à une séquence d'imagerie normale. Pour la deuxième séquence d'excitations N.TE' est remplacé par TE, et N est remplacé par 1. Les paramètres Gij, dij, et m y sont alors différents. Dans ce cas, le signal est dit S1(TE). Les contributions d'un point du milieu sont apportées dans le signal de résonance magnétique en commun avec les contributions des autres points de la tranche examinée dans le milieu. La méthode d'imagerie 2DFT permet en définitive d'affecter à chaque point des images 11 ou 12 une valeur qui est représentative de ces contributions.
  • La figure 3a représente par la courbe enveloppe 13 le signal théorique résultant de l'existence d'un temps de relaxation T2 quand tous les autres effets parasites, y compris la diffusion, peuvent être négligés. La courbe 14 montre ce qu'il advient de cette courbe théorique quand on la mesure avec des séquences dites peu diffusantes : c'est-à-dire quand l'effet de la diffusion se fait peu sentir. Cet effet se fait peu sentir si le nombre N est suffisamment grand pour un temps TE (= NTE') donné : par exemple ici il vaut 4. En effet dans la formule de la contribution, plus N est grand plus les coefficients dij et Iij sont petits et plus l'effet de la diffusion est faible. Les excitations radiofréquence notées 15 à 19 de la figure 3a comportent une excitation tendant à faire basculer à 90° les noyaux du corps (excitation 15) et quatre excitations successives tendant à faire basculer de 180° ces noyaux (16 à 19). Le signal en pointillé 20 est le signal de résonance magnétique nucléaire après la première excitation de 90°. Les sommets 21 à 24 des pics de ce signal, mesuré après des périodes de temps d'échos TE' successives, suivent la courbe 14. Le signal SN de la formule précédente correspond à l'amplitude du pic 24.
  • La figure 3a présente encore l'enveloppe 25 de ces pics lorsque le corps est soumis à des séquences diffusantes. La courbe 25 est par ailleurs reportée sur la figure 3b. Dans une séquence diffusante une seule excitation à 180° 26 suit l'impulsion à 90° 15. Le signal mesuré au bout du temps d'écho TE est coté S1. La valeur de ce signal est donnée par une formule identique à la précédente : la seule différence vient du fait que dans ce cas N vaut 1. On comprend à l'examen de ces deux figures que l'on va pouvoir en comparant une image l1, relative à des pics qui suivent l'enveloppe 14, à une image 12, relative à des pics qui suivent une courbe 25, représenter une image 13 relative à l'écart 27 entre ces deux courbes. Dans la mesure où la séquence peu diffusante peut être de moins en moins diffusante la courbe 14 se rapproche de la courbe 13. Cela permet de représenter une image où l'effet 28 de la diffusion seule peut être encore mieux mis en évidence. Il est utile de noter que la troisième image obtenue avec le procédé de l'invention diffère de l'enseignement donné dans le document précité. Avec cet enseignement, le signal est proportionnel à un écart 29 entre les effets de la diffusion dans le corps étudié et dans un corps étalon dont la réponse en diffusion serait inscrite dans la courbe en tirets 30.
  • Dans l'invention, on effectue le rapport des signaux SN et S1 et on prend le logarithme de ce rapport. On peut ainsi élaborer un signal représentatif (point à point dans l'image) et correspondant à la formule suivante :
    Figure imgb0003
    avec
    Figure imgb0004
  • Cette formule montre que le traitement ainsi effectué sur SN et S1 rend compte, à un facteur K près, de la mesure de la constante de diffusion à l'endroit considéré dans l'image. Mais ce rapport ne donne le résultat ci-dessus indiqué que dans la mesure où, d'une part, les instants temps de prélèvement des signaux, TE et N.TE' , sont égaux entre eux et où, d'autre part, le temps de répétition TR est le même pour les séquences à N échos et pour celles à un écho. La fonction f prend alors la même valeur pour les deux séquences et elle disparaît dans le calcul de l'image 13 où n'apparaît plus que l'effet de la diffusion.
  • Pour que l'effet de la diffusion, courbe 25, se fasse bien sentir il est intéressant d'avoir un temps d'écho TE dans les deuxièmes séquences qui soit suffisamment long et/ou d'utiliser des gradients plus forts, plus longs, que dans une séquence classique. Notamment des gradients 41 à 48 supplémentaires, écartés au maximum des impulsions d'excitation radiofréquence à 180°, peuvent être utilisés en dehors des périodes pendant lesquelles les excitations radiofréquence sont appliquées ou reçues. Par contre pour que l'effet de la diffusion, dans la séquence non diffusante, soit minimisé, on peut prendre un nombre N suffisamment important : N = 4 semble suffisant.
  • Un gradient supplémentaire est montré sur la figure 2a en dessous du graphe des excitations radiofréquence. Il peut être appliqué sur l'un quelconque des axes X, Y ou Z. On peut par exemple l'appliquer sur l'axe de codage de phase ou de sélection. Appliqué sur l'axe de la sélection, il n'interfère pas dans l'épaisseur de la tranche sélectionnée dans le corps dans la mesure où son application, supplémentaire, est faite en dehors des temps d'application des excitations radiofréquence pendant lesquelles le gradient de sélection est lui présent. En effet si le gradient de sélection reçoit pendant les excitations de sélection 90° et de retournement 180° un supplément de valeur, l'inhomogénéité qu'il provoque sera plus forte. Pour une même bande passante donnée de l'excitation, l'épaisseur de la tranche excitée dans le corps 2 sera plus fine. Le signal que restituera cette tranche plus fine sera plus faible. Si on veut qu'il soit aussi fort, il faut alors augmenter la bande passante des excitations et ce procédé nuit alors à la simplicité de la méthode. De manière préférée, le gradient supplémentaire sera appliqué en supplément du gradient de lecture. D'une manière préférée, toutes choses étant égales par ailleurs, on s'arrangera pour que les gradients ajoutés le soient en des périodes les plus éloignées possible de celles où les excitations de basculement, les excitations 12, sont appliquées. -.
  • La figure 3b donne une idée des durées pendant lesquelles on peut appliquer d'une manière préférée les suppléments 41 à 48 de gradient de champ. Ceux-ci sont appliqués en dehors des dates auxquelles on relève les pics du signal de résonance magnétique nucléaire. Dans la représentation les gradients 41 et 48 sont dyssimétriques. Mais comme ils sont homologues compensés par rapport à l'impulsion de basculement 26 leurs intégrales par rapport au temps sont égales. Les grandes impulsions 41 et 48 peuvent être utilisées seules. Elles sont le plus éloignées possible de l'excitation 26.
  • Le choix de l'axe sur lequel les impulsions supplémentaires des gradients sont appliquées peut ne pas être indifférent. Si on considère que la diffusion moléculaire est un phénomène isotrope c'est effectivement le cas. Par contre si on considère que la diffusion moléculaire peut, dans certains cas, être anisotrope, il peut être intéressant de choisir l'axe sur lequel on va favoriser l'effet de diffusion. En particulier certains tissus du corps humain ont par leur disposition dans le corps une orientation privilégiée. Cette orientation privilégiée résulte d'une forme privilégiée des cellules qui les constituent. Ces cellules qui n'ont aucune raison d'être globalement symétriques sur trois axes présentent alors des coefficients de diffusion moléculaires différents selon chacun de ces trois axes. En comparant les images de diffusion moléculaire obtenues selon le procédé de l'invention avec des suppléments de gradients de champ appliqués selon un axe puis selon un autre axe il peut être ultérieurement possible de déterminer à quelle nature de tissu on a à faire.
  • La mesure de la diffusion donnée dans la troisième image est une mesure quantitative à partir du calcul des coefficients bkm (k valant 1 ou N selon les séquences) et donc du facteur K. Toutefois, les imperfections inhérentes aux systèmes d'imagerie peuvent entraîner une atténuation globale des signaux SN ou Si en chaque point de l'image, par exemple épaisseur de coupe légèrement différente, mauvais rephasage de la sélection de coupe d'une des deux images, ... Cette atténuation globale peut être calibrée à l'aide d'une substance étalon 33 disposée le long du corps étudié de façon à ce qu'elle apparaisse sur les bords du champ de l'image. Si SeN et Sei sont les signaux de résonance magnétique dans le corps étalon, correspondant aux signaux SN et Si dans le corps étudié, respectivement à la fin des séquences peu diffusantes et diffusantes, on peut écrire que le coefficient de diffusion molléculaire d'un point de la coupe vaut :
    Figure imgb0005
    où De est le coefficient de diffusion moléculaire connu du corps étalon. Dans ces conditions la calibration est obtenue simplement.
  • En fait dans un tissu vivant deux phénomènes de déplacement se produisent. Un premier phénomène est dû à la diffusion moléculaire comme évoquée jusqu'ici. Un deuxième phénomène est dû aux micro-circulations dans les tissus. Ces micro-circulations relèvent essentiellement de la vascularisation sanguine. Ces micro-circulations viennent perturber l'image de diffusion moléculaire. Aussi, dans une variante, l'invention est mise en oeuvre en modulant l'effet de la vitesse des parties mobiles (le sang) du tissu selon un procédé de modulation dont la description est annexée ci-après. Le procédé de modulation décrit ci-après permet de moduler l'effet (sur le signal de résonance) de ces micro-circulations. Dans ce procédé de modulation on sait, par adjonction de gradients dits bipolaires compensateurs, faire en sorte que dans l'image les molécules du sang soient prises en compte comme si elles étaient fixes. Dans ces conditions le phénomène de diffusion apparaît seul. Ce procédé de modulation a été préalablement décrit dans une demande de brevet français déposée le 13 AOUT 1985 sous le M* 85 12352. Le contenu de cette demande préalable fait maintenant partie intégrante de la présente invention.
  • Pour le perfectionnement ainsi apporté, on peut expliquer que les micro-circulations (elles sont des déplacements lents) perturbent les séquences diffusantes à cause des gradients de sensibilisation suffisamment forts que ces séquences comportent. C'est donc d'une manière préférée les séquences diffusantes qui sont compensées. La question se pose de savoir comment cette compensation, recherchée par ce qu'elle neutralise l'effet indésirable des micro-circulations, ne neutralise pas en même temps l'effet de la diffusion. Les inventeurs de la présente invention pensent que cette non-neutralisation est dûe à une distribution temporellement incohérente des déplacements de diffusion moléculaires. Par contre les micro-circulations sont des déplacements cohérents puisque leurs vitesses sont quasiment constantes. Autrement dit, la réduction du signal de résonance mesurée au cours d'une séquence diffusante (par rapport à celui mesuré au cours d'une séquence non diffusante) est dûe maintenant uniquement à la diffusion et non plus aussi aux micro-circulations dont l'effet a été neutralisé.
  • Le perfectionnement 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. Ce perfectionnement trouve son application particuliè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, il peut même plus particulièrement être mis 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 examiné.
  • Au cours d'une expérience de résonance, si le champ orientateur Bo 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'une composante de phase dépendant justement de leur vitesse. Ceci peut se comprendre simplement : le signal de résonance émis vibre à une fréquence fo qui dépend de l'intensité du champ magnétique orientateur Bo et du rapport gyromagnétique y caractéristique du milieu examiné. Toutes variations de l'intensité du champ Bo 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 fo, dans un deuxième temps, soumise à un champ plus fort Bo + à Bo, elle résonne à une fréquence plus élevée fo + Afo. Dans un troisième temps à nouveau soumise au champ Bo elle vibre à nouveau à la fréquence fo. 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 à 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 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 mutuellement 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 ra-- - diofré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 positive 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.
  • Quels 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 perturbateur 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 sensibilité des machines de RMN est actuellement de l'ordre de 1 radian/(cm/s) à 0,2 rd/(cm/s). Ceci signifie qu'une particule qui se déplace de 1 cm par seconde dans le plan de la coupe contribue au signal global émis 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 couramment 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 radians s'opposent mutuellement, le signal résultant est nul : cela revient à prendre la valeur moyenne d'un signal sinusoidal 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.
  • Le perfectionnement 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 modulation permet par ailleurs, plutôt que de tout annuler, de moduler l'effet de la vitesse. 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 modulation 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 modulation est très général. Il n'est pas cantonné à 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 transformé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 le perfectionnement sera décrit ), 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 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.
  • Le perfectionnement 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.
  • Le perfectionnement 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.
  • La figure 4 représente une machine utilisable pour mettre en oeuvre le procédé du perfectionnement. Elle comporte des moyens symbolisés par une bobine 1 pour soumettre un corps 2 à un champ magnétique important et constant Bo. Le champ Bo 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 5). 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 11 et une deuxième image 12 d'une coupe 70 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 30. Dans des moyens de comparaison 80 on compare point à point les images 11 , 12, et on élabore une troisième image l3. On peut visualiser cette image is sur un dispositif de visualisation 50.
  • Dans les images 11 , 12, ou l3, un point d'image, ou pixel p, représente classiquement par sa luminosité la densité de particules contenue dans un élément de volume, ou voxel y qui lui correspond dans la tranche 70. 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 information de vitesse, par comparaison entre deux images 11 et 12 où l'effet de cette vitesse aura été modulée de deux manières diffé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'inhomo- généité du champ orientateur. Dans un troisième temps, on indiquera, dans cet exemple, comment se calculent les champs compensateurs nécessaires 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 Bo 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, Cz). 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 5 montre à cet effet des impulsions radiofréquence 100 , 110 et 120 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 impulsions 130 , 140 et 150. 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 5, le gradient X est coté 160, le gradient Y est coté 170. 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'obtenir simultanément des images de plusieurs couples parallèles.
  • L'existence du précodage 170, 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 impulsions radiofréquence 110 puis 120 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 5 plus hautes que l'impulsion 100 puisqu'elles font basculer l'orientation des moments magnétiques des spins des particules de 180° alors que l'impulsion 100 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 160 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 11 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 180 de l'impulsion 130, pour la partie de cette impulsion existant après la fin de l'impulsion radiofréquence 100, et par le début 190 de l'impulsion 160 de lecture il est connu d'appliquer sur ces axes respectivement des impulsions 200 et 210. 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 200 et 190 sont donc égales aux surfaces en pointillés 180 et 210 respectivement. Les impulsion 200 et 210 n'interfèrent pas dans la sélection de coupe puisqu'au moment où elles sont appliquées l'impulsion 100 n'est plus présente. Les impulsions 140 et 150 sont autoneutralisées puisqu'elles comportent des portions respectivement 220 et 230, et 240 et 250 antisymétriques par rapport au milieu de l'impulsion radiofréquence, 110 et 120 pour lesquelles elles resélectionnent la coupe 70 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 Bo 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 lecture X) et de même amplitude, GL, d'une séquence à l'autre. L'inhomogénéité du champ Bo ne fait alors que déformer les lignes isofréquences de résonance. Pendant la lecture du signal, l'inhomogé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 6a et 6b montrent le déphasage Δφ 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 6a montre le diagramme temporel d'un signal radiofréquence RF comportant une impulsion 260 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 260, par exemple, une valeur G pendant une durée 8 à chaque fois. En un lieu donné la une particule préalablement excitée, immobile, vibre à une fréquence fo au début de l'application de l'impulsion 270. Pendant la durée e de cette impulsion 270 elle vibre à une fréquence différente. Par exemple cette fréquence est supérieure, de telle manière qu'à la fin de l'impulsion 270 on puisse considérer, puisqu'elle se remet à vibrer à la fréquence fo qu'elle a subi un déphasage φ1. Lors de l'application de l'impulsion radiofréquence 260 on peut admettre que tout se passe comme si la phase du signal émis par cette particule s'inversait : - φ1. Pendant l'impulsion 280 qui a exactement la même allure que l'impulsion 270 le phénomène de déphasage 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 260. En conséquence, la particule immobile en la retrouve à la fin de l'impulsion 280 la phase qu'elle avait au début de l'impulsion 270. Il est utile de remarquer que le déphasage φ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 ç1 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 6a, on a indiqué que cette particule se trouvait successivement dans les positions la à 14 au début et à la fin de la première impulsion 270, au moment de l'impulsion radiofréquence 260, et au début et à la fin de la deuxième impulsion 280 du gradient. Après la première impulsion 270 de gradient, le déphasage résultant φ2 est proportionnel au gradient G et à la durée e Mais il est aussi proportionnel à la distance h - Io parcourue dans l'intervalle par la particule. Cette distance est égale à V. e. Ce qui fait que le temps intervient maintenant au carré : la première portion de déphasage à l'allure d'une parabole (V.G.e2). De la fin de l'impulsion 270 au début de l'impulsion 280, tout se passe comme pour la particule immobile : la contribution possède une phase - φ2 au début de l'impulsion 280. Au cours de cette impulsion 280 l'évolution de la phase est également parabolique comme pour l'impulsion 270. On pourrait s'attendre à ce que les effets de phase s'effacent ici également. Cependant entre la fin de l'impulsion 270 et le début de l'impulsion 280 la particule animée s'est déplacée de la position h à la position l3. Dans la position 13 au moment où l'impulsion 280 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 Io et h. En conséquence à la fin de l'impulsion 280 la compensation ne s'opère pas et cette particule animée émet alors un signal déphasé d'une valeur Δφ que l'on peut écrire de la manière suivante :
    Figure imgb0006
    où y est le rapport gyromagnétique des particules considérées, où G (e) représente le gradient dans les impulsions 270 et 280, et où T est la durée qui sépare l'apparition de ces deux impulsions. Si 0 n'est pas petit, T mesure la durée entre les milieux de ces impulsions. Plus généralement on peut écrire :
    Figure imgb0007
    où le déphasage est donné maintenant en fonction du temps et où G(τ) représente toutes les perturbations appliquées aux particules avant qu'on en mesure à une date le signal qu'elles émettent. On se souvient simplement que les perturbations considérées sont telles que :
    Figure imgb0008
  • Car c'est sous cette condition que, d'une manière connue dans l'état de la technique (figure 5 impulsions 180 à 250), on a compensé les précodages parasites. On ne l'avait bien entendu pas fait pour l'impulsion variable 170 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 260 créant un écho de spin (figure 6b), les impulsion 270 et 280 du gradient bipolaire doivent être remplacées par des impulsions 290 , 300 de signes opposés. De cette manière la phase pour les particules fixes évolue jusqu'à φ3 en fin de l'impulsion 290 et évolue en sens inverse jusqu'à zéro au cours de l'impulsion inverse 300. Le déphasage Δφ entre la contribution des particules mobiles et celles des particules fixes est du même ordre que précédemment.
  • Dans l'expression de Δφ (t), G (t) représente donc toutes les séquences nécessitées par la mesure. Ce qui caractérise le perfectionnement 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
    Figure imgb0009
    où G' représente des gradients, eux aussi bipolaires, c'est-à-dire tels que
    Figure imgb0010
  • Autrement dit en agissant ainsi on annule le coefficient multiplicateur de la vitesse V qui apparaît dans le calcul de Δφ (t). Du coup, quelle que soit cette vitesse, Δφ (t) est nul : toutes les contributions sont en phase. Les particules mobiles sont prises en compte comme si elles étaient fixes. Il faut cependant remarquer que seules les modifications du signal dues à des vitesses V, constantes pendant la mesure, sont ainsi annulées. Ceci est en particulier le cas pour les micro-circulations de cellules de sang. Par contre, avec la diffusion, la vitesse de protons au cours d'une même séquence (NTE', TE) n'est pas constante. Donc le terme V ne peut être sorti de l'expression de Δφ. La compensation de G par G' ne peut pas se produire. Dans ces conditions, l'atténuation du signal de résonance, due à la diffusion, persiste. Elle est mesurée et sa mesure permet la différenciation du milieu étudié. Il est à noter que le phénomène de micro-circulation est particulièrement présent dans les angiomes et les tumeurs. D'où l'intérêt de supprimer ses effets pour mesurer correctement la diffusion. Par ailleurs, on peut quantifier l'intensité des micro-circulations par comparaison d'une image de diffusion non compensée en vitesse à une image de diffusion compensée en vitesse. En effet la première image est une image de diffusion et de micro-circulation alors que la deuxième n'est qu'une image de diffusion pure. La soustraction des deux images, si leurs caractéristiques sont égales, conduit à une image de micro-circulation.
  • En revenant à la figure 5 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 310 et 320 , 330 et 340 , 350 et 360 , 370 et 380 , 390 et 400 et 420 et 430. 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 200 sur l'axe de sélection Z est appliquée entre des dates t1 et t2 à la fin de l'impulsion 130, on applique l'impulsion 310 pendant la même période. Dans ce cas l'impulsion 200 "nouvelle" est plus forte. Les impulsions 350 et 370 sont également appliquées pendant cette période. Les impulsions 320 , 360 et 380 d'une manière préférée sont appliquées entre la fin t3 de l'impulsion 140 et le début t4 de l'impulsion de lecture 160. Les impulsions 330 , 390 et 420 sont appliquées entre la fin ts de l'impulsion 160 et le début t6 de l'impulsion 150. Les impulsions 340 , 400 et 430 sont appliquées entre la fin t7 de l'impulsion 150 et le début t8 d'une impulsion 410 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's contenues entre les dates t4 et ts et aussi pendant la durée de l'impulsion 410. 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 (τ).τ 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 310 à 430 seront des formes connues et, dans la mesure du possible, identiques à des formes déjà pratiquées pour les impulsions normales des séquences perturbatrices. 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 :
    Figure imgb0011
  • Le premier membre est connu il concerne les gradients perturbateurs. Le deuxième membre comporte des intégrales de gradients de forme préférée r et d'amplitude unité. On sait calculer l'intégrale du deuxième membre puisqu'on connaît la forme choisie r. En conséquence on peut en déduire λ, 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's), l'équation implicite contenant G' ne peut être vérifiée pour tout t. Mais on doit au moins la vérifier aux temps tEi , tE2... centre des phases de mesure. Ces dates sont celles où l'écho de spin attendu est le plus fort. Elles sont telles que la durée séparant la date de l'impulsion radiofréquence 100 et tEi est le double de la durée entre l'impulsion 100 et l'impulsion 110, et telles que la durée entre la date tE1 et la date tE2 est égale au double de la durée entre la durée tE1 et la date de l'impulsion radiofréquence 120.... La justification de 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 70 non pas au moment de son excitation (impulsion 100) mais au moment, tE1 ou tE2, où justement on effectue la mesure : au centre du signal émis.
  • Tous les termes de l'équation implicite donnant G' étant connus, G' est calculable. Sur la figure 5 on va donner une approche simple de ce calcul. Pour l'axe X on sait que les impulsions 210 et 190 situées de part et d'autre de l'impulsion radiofréquence 110 neutralisent 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é, tEi, n'est pas nulle : l'impulsion 210 étant plus antérieure dans le temps que l'impulsion 190. Tout se passe donc comme si l'impulsion 210 de "force" donnée est appliquée à un "bras de levier" plus long que l'impulsion 190 de force égale. L'effet sur la phase de l'impulsion 190 est bien entendu inverse de celui de l'impulsion 210 puisqu'elles sont situées de part et d'autre de l'impulsion 110. 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 210. Il convient de contrecarrer ce moment résultant par un couple d'impulsions 350 , 360 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 210 et 190. Ainsi, l'impulsion 350 étant placée dans le temps à une date plus antérieure que l'impulsion 360 par rapport à la date tE1 son effet est prépondérant. II est bien de sens opposé à l'impulsion 210.
  • De la même façon sur l'axe de sélection Z, des impulsions 180 et 200 , et 220 et 230, ce sont à chaque fois les impulsions 180 et 220 qui priment : leur bras de levier est plus long. En conséquence on conçoit bien que l'impulsion 310 doit être dans le même sens que l'impulsion 200. Le cas des impulsions 370 et 380 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 370 et 380 varieront donc également en correspondance d'une acquisition à l'autre (symbolisé par des flèches). Pour les raisons précédentes quand l'impulsion 170 est positive l'impulsion 370 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 (310 , 350 , 370) et le plus en aval pour la deuxième (320 , 360 , 380). En conclusion elles sont le plus éloignées possible de part et d'autre de l'impulsion radiofréquence 110 ; et elles sont le plus éloignées possible l'une de l'autre.
  • Lors du premier relevé, en tEi, les effets de vitesse des parties mobiles ont été complètement compensés. On compense maintenant au moyen des impulsions 330 et 340 les effets de la deuxième impulsion de sélection 150. La date pivot pour le calcul des moments est maintenant la date tE2- Sur l'axe de lecture X on compense également par deux impulsions 420 et 430 les effets de la deuxième moitié 600 de l'impulsion 160 et de la première moitié 610 de l'impulsion 410. La présence des impulsions compensatrices 390 et 400 sur l'axe Y se justifie par le fait que l'impulsion 170 n'est pas bipolaire. Car la compensation de son moment à la date pivot tE1 par les impulsions 370 et 380 n'implique pas une telle compensation à une autre date pivot telle que tE2. Et ceci bien qu'aucune impulsion supplémentaire de perturbation n'intervienne sur cet axe entre temps.
  • La figure 7 représente un tube d'écoulement 510 sensé se trouver dans la tranche examinée 70 du corps 2. Le haut de la figure représente une situation dans laquelle le tube 510 est coplanaire au plan X , Y de la coupe. Dans la partie basse de la figure la situation est différente : le tube 510 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 520 à 550. 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'échantillonnage 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 parallélépipède rectangle s'étend selon l'axe Z. Cette longueur est liée à l'épaisseur de la coupe 70. Cette épaisseur est un compromis expérimental 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 1 mm.
  • Dans le conduit 510 un fluide non représenté circule avec un diagramme de vitesse 560. Les vitesses 570 des particules situées au centre sont plus élevées que les vitesses 580 des particules situées sur les bords du conduit. Dans les artères la vitesse 570 peut être de l'ordre de 50 cm par seconde. Les voxels 520 et 530 comportent, quand ils sont situés à l'aplomb de la section 590 du conduit 510, nécessairement 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 520 et 530 dans le signal global était nulle : on avait l'impression que le tube était creux. L'image des voxels 540 et 550 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 550 où les vitesses (au centre du tube) tout en étant élevées sont relativement homogènes entre elles. Leur déphasage correspondant est sensiblement le même : la moyenne n'est pas nulle. Par contre pour le voxel 540 situé à l'aplomb du bord du conduit 510, là ou l'évolution des vitesses est brutale, le signal restitué est nul. Il en résultait que l'image des tubes 510 perpendiculaires au plan de coupe était appréciée géométriquement plus petite qu'elle n'était en réalité. Avec le perfectionnement l'image des écoulements est révélée fidèlement.
  • Un autre perfectionnement a été apporté à l'invention. Les inventeurs se sont rendus compte en effet, à l'expérience, que la compensation de vitesse du sang dans les capillaires liée à la micro-circulation, pour ne faire apparaître que le phénomène de diffusion seul, n'était pas toujours efficace parce que les vitesses de micro-circulation dans le tissu n'étaient pas constantes en valeur et en direction pendant le temps de mesure (temps d'écho). La constance en direction est liée à la géométrie du réseau capillaire dans le tissu étudié : par exemple la direction de la micro-circulation dans un muscle peut-être constante ( vaisseaux rectilignes) sur une distance supérieure à la distance parcourue par le sang pendant le temps de mesure. Par contre pour d'autres tissus, comme dans le cerveau humain, la micro-circulation (on dit aussi quelques fois la perfusion) est orientée en tous sens. L'enchevêtrement des réseaux de capillaires y est même tel que, malgré la rapidité des mesures RMN (une séquence excitation-mesure peut durer typiquement 100 ms) les orientations des vitesses de déplacement des molécules changent de direction pendant la mesure : ce qui peut-être interprété comme une macro-diffusion à l'échelle d'un voxel. Autrement dit la représentation dans l'image du paramètre de diffusion D du tissu est faussée : le coefficient représenté contient une composante liée à la valeur de la diffusion vraie (D) pour la part statique du tissu et une composante (D') liée à la présence de la micro-circulation dans un voxel donné.
  • Si on appelle f la proportion de fluide circulant dans les capillaires et (I-f) la proportion du reste du tissu dans un voxel donné, on peut dire que la contribution de signal mesurée relative à ce voxel peut s'écrire :
    • - pour une séquence peu diffusante (peu sensible à la diffusion et au phénomène de micro-circulation) SN=(1-f)S(ρ,T1,T2).A + f.S'(p,T1,T2).A.FN(avecNechos,Ne 1)
    • - et pour une séquence diffusante:
    • Si = (1-f)S(p,Tt,T2).e-bkmD+ f.S'(p.Ti.T2).Fi.e-bkm.D
  • Dans ces formules S (p,Ti, T2) et S'(p,Ti, T2) représentent le signal de RMN respectivement pour le tissu (le tissu du cerveau) et pour le fluide des capillaires actifs. La valeur A = ebNmD représente l'effet de la diffusion pendant les premières séquences (N , m). Comme les premières séquences sont peu diffusantes A vaut 1 et est par la suite négligé dans les calculs. La fonction F, égale à 1 pour les séquences peu diffusantes, a une valeur inférieure à 1 en cas de séquences diffusantes. Elle dépend des gradients, de la vitesse du sang dans la micro-circulation, et de la géométrie des capillaires. Quand les capillaires sont rectilignes (sur une distance fonction du temps de mesure) mais avec des orientations différentes au sein d'un même voxel, F est un sinus cardinal du gradient et de la vitesse de micro-circulation. Quand la distance moyenne pour laquelle un capillaire peut être considéré comme rectiligne est inférieure à la distance parcourue pendant le temps de mesure on peut écrire :
    Figure imgb0012
  • La micro-circulation est assimilée a une macro-diffusion.
  • Dans l'autre perfectionnement ainsi apporté on s'est rendu compte que Fi, lié à la micro-circulation, pouvait être bien plus petit que e - bkm.D lié à la diffusion D vraie à connaître. On s'arrange alors pour que bkm, qui dépend de la forme et de la force des gradients dans la séquence diffusante, soit suffisant pour rendre l'atténuation e-bkm.D sensible (au point qu'on peut mesurer D) et soit à ce point fort que f.Fi puisse être considéré comme négligeable. Dans un exemple préféré bkm correspond à un k = 1 : les séquences diffusantes comportent de préférence un seul écho de spin (m correspondant au nombre de paires de gradients homologues compensés, aussi dits bipolaires, de la séquence). Dans ce cas on peut dire que le coefficient de diffusion mesuré Dm est, en application de la formule précédente donnée au début de cette description, tel que :
    Figure imgb0013
  • En supposant que S est voisin de S', on peut simplifier par S (ou S'), et l'expression de Dm devient :
    Figure imgb0014
  • Cette expression est une équation à deux inconnues, D et f, dont on peut trouver une solution en réitérant les expériences de manière à mesurer deux valeurs de Dm pour deux valeurs de bkm différentes. Du système de deux équations à deux inconnues obtenu on extrait D et f. Les deux valeurs de bkm retenues sont conformes par ailleurs aux indications données ci-dessus relatives à la nécessité de rendre f.F1 négligeable.
  • Autrement dit le nouveau procédé, ainsi mis en oeuvre pour procurer les fruits de cet autre perfectionnement, doit comporter :
    • - un (premier) ensemble de séquences d'un premier type (des premières séquences peu diffusantes),
    • - un (deuxième) ensemble de séquences d'un deuxième type (des deuxièmes séquences diffusantes),
    • - et un (troisième) ensemble de séquence du deuxième type (des deuxièmes séquences, là aussi diffusantes, mais dont les caractéristiques de diffusion, bkm, sont changées). La première valeur de Dm est tirée de la comparaison du signal de RMN résultant du premier ensemble de séquences au signal de RMN résultant du deuxième ensemble de séquences. La deuxième valeur de Dm est tirée de la comparaison du signal de RMN résultant du premier ensemble de séquences au signal de RMN résultant du troisième ensemble de séquences. Il est inutile en effet de réiterer le premier ensemble de séquences. La valeur vraie D de la diffusion en chaque point est extraite du système des deux équations a deux inconnues qui vaut pour ce point. On sait résoudre un tel système d'équations. On peut aussi par le procédé de l'invention connaître f, ou I-f ce qui revient au même. Il est à noter que la connaissance de ce paramètre peut aussi être intéressante pour établir un diagnostic. On obtient finalement deux images : une de diffusion vraie D, l'autre de perfusion f (en fait de la proportion de capillaires actifs dans le tissu, dans chaque voxel). Le cas échéant, les séquences diffusantes peuvent aussi être compensées en vitesse si une composante rectiligne de la vitesse dans les capillaires est remarquée. Dans ce cas, du fait de la compensation, F1 devient égal à 1.
  • On remarque en outre que pour le deuxième et le troisième ensemble de séquences (diffusantes) f. Fi est négligeable. En effet pour ces séquences la valeur (forme et force) des gradients est adaptée à cette contrainte. On peut alors comparer point à point le signal de RMN relatif à ces séquences. Soit Si et Si le signal relatif, en chaque point, respectivement aux deuxième et troisième séquences, on peut écrire:
    Figure imgb0015
    Figure imgb0016
    puisque la contribution de la micro circulation dans chaque pixel est négligeable. Il en resulte que :
    Figure imgb0017
  • D'où on extrait D vrai. Dans ces conditions, on élimine les parasites de micro circulation, tout en ne faisant que deux séquences, mais deux séquences diffusantes (Si, S'i), avec des caractéristiques de diffusion différentes. En conséquence quand on dit, en particulier dans les revendications ; que le procédé de l'invention comporte des séquences peu diffusantes et des séquences diffusantes il faut entendre en fait, toutes choses étant égales par ailleurs, des séquences moins diffusantes et des séquences plus diffusantes.
  • En outre trois cas peuvent se présenter :
    • - premier cas - les capillaires du tissu sont rectilignes et parallèles-on obtient alors l'image de diffusion vraie avec une séquence non diffusante et une séquence diffusante compensée.
    • - deuxième cas - les capillaires ne sont pas rectilignes - on obtient alors l'image de diffusion vraie D, ou l'image de perfusion f, avec les trois séquences (dont deux diffusantes) et il ne sert à rien de compenser.
    • - troisième cas - les capillaires sont rectilignes mais dans tous les sens - on obtient l'image de diffusion vraie comme dans le premier cas ou comme dans le deuxième cas. Le deuxième cas nécessite une acquisition supplémentaire de l'image (en faisant varier bkm) ce qui est plus long mais qui donne l'avantage de proposer en outre l'image de perfusion f.

Claims (24)

1. Procédé d'imagerie par résonance magnétique nucléaire (1-9, 50) comportant les étapes suivantes pour donner une image (13) de la diffusion moléculaire d'un corps étudié :
- on place le corps (2) dans un champ magnétique constant (Bo);
- on soumet (3-4, Ci) le corps ainsi placé à une première pluralité de premières séquences d'excitations à écho de spin en présence de premières séquences de gradient de champ ; ces premières séquences à écho de spin comportant un nombre entier N supérieur ou égal à 1 d'excitations (16-19) où les moments magnétiques des noyaux sont basculés de 180° après une excitation (15) où ces moments ont été basculés de 90° pour faire ainsi des séquences à N échos peu diffusantes ;
- on relève (4,5) les signaux de résonance magnétique en fin de cespremières séquences d'excitations et on calcule (6) une première image (11) correspondant aux échos N de ces signaux en attribuant à chaque point de l'image une valeur (SN) correspondant au signal de résonance magnétique du point qui lui correspond dans le corps;
- on soumet (3-4, C2) le corps ainsi placé à une deuxième pluralité de deuxièmes séquences d'excitations à écho de spin en présence de deuxièmes séquences de gradient de champ ; ces deuxièmes séquences à écho de spin comportant au moins une excitation à 180° (26) après l'excitation à 90° (15) pour faire des séquences à au moins un écho diffusantes ; la durée totale d'écho (TE) des deuxièmes séquences d'excitations étant égale à la durée totale d'écho (N.TE') des premières séquences d'excitations ;
- on relève (4,5) les signaux de résonance magnétique en fin de ces deuxièmes séquences d'excitations et on calcule (6) une deuxième image (12) correspondant aux échos de ces signaux en attribuant à chaque point d'image une valeur correspondant au signal de résonance magnétique du point qui lui correspond - - dans le corps.;
- on compare (7) point à point les valeurs attribuées pour la première image aux valeurs attribuées pour la deuxième image, pour élaborer une troisième image (13) représentative de la diffusion moléculaire en chaque point du corps.
2. Procédé selon la revendication 1, caractérisé en ce que les deuxièmes séquences comportent des impulsions plus longues et/ou plus fortes (41-48) de gradients de champ hors présence des excitations d'écho de spin et orientés selon un axe pour constituer des suppléments de gradient de champ.
3. Procédé selon la revendication 2, dans lequel les séquences de gradient de champ comportent des impulsions de gradient de champ orienté selon trois axes (X,Y,Z) : l'axe (Z) du champ constant (Bo), ou axe de sélection, et deux axes (X,Y) orthogonaux à cet axe appelés respectivement axe de codage de phase (Y) et axe de lecture (X), dans lequel la méthode de calcul des images est du type 2DFT, caractérisé en ce que l'axe des suppléments est l'axe de lecture.
4. Procédé selon la revendication 2, dans lequel les séquences de gradient de champ comportent des impulsions de gradient de champ orienté selon trois axes (X,Y,Z) : l'axe (Z) du champ constant (Bo), ou axe de sélection, et deux axes (X,Y) orthogonaux à cet axe appelés respectivement axe de codage de phase (Y) et axe de lecture (X), dans lequel la méthode de calcul des images est du type 2DFT, caractérisé en ce que l'axe des suppléments est l'axe de sélection.
5. Procédé selon la revendication 2, dans lequel les séquences de gradient de champ comportent des impulsions de gradient de champ orienté selon trois axes (X,Y,Z) : l'axe (Z) du champ constant (Bo), ou axe de sélection, et deux axes (X,Y) orthogonaux à cet axe appelés respectivement axe de codage de phase (Y) et axe de lecture (X), dans lequel la méthode de calcul des images est du type 2DFT, caractérisé en ce que l'axe des suppléments est l'axe de codage de phase.
6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il est mis en oeuvre à plusieurs reprises pour réaliser des troisièmes images multicoupes du corps étudié.
7. Procédé selon l'une quelconque des revendications 2 à 6, caractérisé en ce que les instants d'application des suppléments de gradient sont le plus éloignés possible dans le temps des excitations radio- fréquence à 180°.
8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'on place à côté du corps un corps étalon (33) pour calibrer les calculs de la troisième image.
9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les durées de temps d'écho (TE') des premières séquences d'écho de spin sont toutes égales entre elles.
10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que N vaut 4.
11. Procédé selon la revendication 2, caractérisé en ce qu'on réitère les étapes relatives à la deuxième image avec des suppléments de gradient orientés selon un autre axe pour élaborer une autre troisième image pour déterminer la nature des régions imagées.
12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que pour la comparaison point par point, on calcule le logarithme du rapport desdites valeurs.
13. Procédé selon l'une quelconque des revendications 1 à 12, ou 21 ou 22, caractérisé en ce qu'on module l'effet de la vitesse (V) des parties mobiles du corps créé par des séquences de gradient de champ dits perturbateurs par l'application, avant le relevé (tEi), d'une séquence (310 - 430) 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 des champs perturbateurs.
14. Procédé selon la revendication 13, caractérisé en ce que la séquence du champ magnétique perturbateur comporte des impulsions 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.
15. Procédé selon la revendication 14, 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, et en ce que leur amplitude est évaluée pour obtenir la modulation recherchée.
16. Procédé selon la revendication 13 ou 14 ou 15, caractérisé en ce que la séquence des champs magnétiques compensateurs comporte des couples bipolaires d'impulsions.
17. Procédé selon l'une quelconque des revendications 13 à 16, 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 excitation radio-fréquence.
18. Procédé selon l'une quelconque des revendications 13 à 17, 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.
19. Procédé selon la revendication 16 ou 17, caractérisé en ce que les impulsions des couples d'impulsion comportent des impulsions les plus éloignées possible l'une de l'autre.
20. Procédé selon l'une quelconque des revendications 13 à 19, caractérisé en ce qu'on module l'effet de la vitesse des parties mobiles du corps créé seulement par les deuxièmes séquences de gradient de champ.
21. Procédé selon la revendication 1 caractérisé en ce que :
- on soumet le corps ainsi placé à une troisième pluralité de troisième séquences d'excitation à écho de spin en présence de troisièmes séquences de gradient de champ ; les troisièmes séquences à écho de spin comportant au moins une excitation à 180 ° après l'excitation à 90 ° pour faire des séquences à au moins un écho diffusantes ; la durée totale d'écho des troisièmes séquences d'excitation étant égale à la durée totale d'écho des premières séquences d'excitation ; les troisièmes séquences de gradient de champ étant différentes des deuxièmes séquences de gradient de champ ;
- on relève les signaux de résonance magnétiques en fin de ces troisièmes séquences d'excitation et on calcule une quatrième image correspondant aux échos de ces signaux en attribuant à chaque point d'image une valeur correspondant au signal de résonance magnétique du point qui lui correspond dans le corps ;
- on compare point à point les valeurs attribuées pour la première image aux valeurs attribuées pour la quatrième image pour élaborer une cinquième image représentative de la diffusion moléculaire en chaque point du corps ;
- on compare point à point les valeurs attribuées pour la troisième image aux valeurs attribuées pour la cinquième image pour élaborer une sixième image représentative de la diffusion moléculaire vraie dans le corps et exempte de parasites de micro-circulation.
22. Procédé selon la revendication 1 caractérisé en ce que :
- on soumet le corps ainsi placé à une troisième pluralité de troisièmes séquences d'excitation à écho de spin en présence de troisièmes séquences de gradient de champ : les troisièmes séquences à écho de spin comportant au moins une excitation à 180° après l'excitation à 90° pour faire des séquences à au moins un écho diffusantes ; la durée totale d'écho des troisièmes séquences d'excitation étant égale à la durée totale d'écho des premières séquences d'excitation; les troisièmes séquences de gradient de champ étant différentes des deuxièmes séquences de gradient de champ ;
- on relève les signaux de résonance magnétiques en fin de ces troisièmes séquences d'excitation et on calcule une quatrième image correspondant aux échos de ces signaux en attribuant à chaque point d'image une valeur correspondant au signal de résonance magnétique du point qui lui correspond dans le corps ;
- on compare point à point les valeurs attribuées pour la première image aux valeurs attribuées pour la quatrième image pour élaborer une cinquième image représentative de la diffusion moléculaire en chaque point du corps ;
- on compare point à point les valeurs attribuées pour la troisième image aux valeurs attribuées pour la cinquième image pour élaborer une septième image repésentative d'un phénomène de perfusion dans le corps.
23. Procédé selon les revendications 21 et 22 caractérisé en ce qu'on élabore simultanément la sixième et la séptième image.
24. Procédé selon l'une quelconque des revendications 1 à 23 caractérisé en ce que les séquences diffusantes sont des séquences à un seul écho de spin.
EP86401423A 1985-06-27 1986-06-27 Procédé d'imagerie par résonance magnétique nucléaire Expired EP0208601B1 (fr)

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
FR8509824 1985-06-27
FR8509824A FR2584188B1 (fr) 1985-06-27 1985-06-27 Procede d'imagerie par resonance magnetique nucleaire
FR8512352 1985-08-13
FR8512352A FR2586296B1 (fr) 1985-08-13 1985-08-13 Procede de modulation de l'effet de la vitesse des parties mobiles d'un corps dans une mesure de densite par resonance magnetique nucleaire, et mise en oeuvre du procede pour en deduire la vitesse des parties mobiles concernees
US82352286A 1986-01-29 1986-01-29
US823522 1992-01-21

Publications (2)

Publication Number Publication Date
EP0208601A1 EP0208601A1 (fr) 1987-01-14
EP0208601B1 true EP0208601B1 (fr) 1989-10-25

Family

ID=27251293

Family Applications (1)

Application Number Title Priority Date Filing Date
EP86401423A Expired EP0208601B1 (fr) 1985-06-27 1986-06-27 Procédé d'imagerie par résonance magnétique nucléaire

Country Status (3)

Country Link
US (1) USRE33391E (fr)
EP (1) EP0208601B1 (fr)
DE (1) DE3666632D1 (fr)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2621693A1 (fr) * 1987-10-13 1989-04-14 Thomson Cgr Procede d'imagerie de mouvements intravoxels par rmn dans un corps
US5093620A (en) * 1990-08-09 1992-03-03 General Electric Encoding for nmr phase contrast flow measurement
WO1995004940A1 (fr) * 1993-08-06 1995-02-16 The Government Of The United States Of America, Represented By The Secretary Of The Department Of Health And Human Services Procede et systeme servant a mesurer le tenseur de diffusion et a generer des images correspondantes
US6614226B2 (en) 2000-03-31 2003-09-02 The General Hospital Corporation Diffusion imaging of tissues
NL1025819C2 (nl) * 2004-03-26 2005-09-27 Lely Entpr Ag Werkwijze en inrichting voor het melken van een melkdier.
US8540821B2 (en) * 2004-03-26 2013-09-24 Maasland N.V. Teat cup cleaning device and method
SE0702063L (sv) * 2007-05-31 2009-01-13 Colloidal Resource Ab Metod, system, datoravläsbart medium och användning för magnetisk resonanstomografi
US9097778B2 (en) * 2009-03-30 2015-08-04 Hitachi, Ltd. Magnetic resonance device with diffusion gradient phase variation positionally corrected
US8729893B2 (en) * 2010-10-19 2014-05-20 Baker Hughes Incorporated Nuclear magnetic resonance 1H and 13C multiphase flow measurements, estimating phase selected flow rates from velocity distributions, volume fractions, and mean velocity
WO2012059581A1 (fr) 2010-11-04 2012-05-10 Micronic Mydata AB Procédé et dispositif utilisant un déflecteur acousto-optique (aod) perfectionné et un balai dense de spots mobiles

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4307344A (en) * 1979-01-25 1981-12-22 Emi Limited Imaging systems
US4616180A (en) * 1983-11-14 1986-10-07 Technicare Corporation Nuclear magnetic resonance imaging with reduced sensitivity to motional effects
FR2586296B1 (fr) * 1985-08-13 1988-06-17 Thomson Cgr Procede de modulation de l'effet de la vitesse des parties mobiles d'un corps dans une mesure de densite par resonance magnetique nucleaire, et mise en oeuvre du procede pour en deduire la vitesse des parties mobiles concernees
US4683431A (en) * 1985-08-16 1987-07-28 Picker International, Inc. Magnetic resonance imaging of high velocity flows

Also Published As

Publication number Publication date
EP0208601A1 (fr) 1987-01-14
USRE33391E (en) 1990-10-16
DE3666632D1 (en) 1989-11-30

Similar Documents

Publication Publication Date Title
LU81529A1 (fr) Procede et appareil de determination des densites de noyaux par echo de spin pour elaboration de tomogrammes
FR2546642A1 (fr) Procede de generation d'images du temps de relaxation transversale en resonance magnetique nucleaire, par l'utilisation de sequences d'echos de spin multiples
EP0188145B1 (fr) Procédé de génération et de traitement de signaux pour l'obtention par résonance magnétique nucléaire d'une image exempte de distorsions à partir d'un champ de polarisation inhomogène
EP0210038A2 (fr) Imagerie du déplacement chimique avec des corrections de l'inhomogénéité de champ
EP3295201B1 (fr) Procede et dispositif d'imagerie par resonance magnetique
EP0208601B1 (fr) Procédé d'imagerie par résonance magnétique nucléaire
FR2520117A1 (fr) Procede de compensation d'une oscillation acoustique dans des mesures de resonance magnetique nucleaire
EP0312427B1 (fr) Procédé d'imagerie de mouvements intravoxels par RMN dans un corps
Ropele et al. Super‐resolution MRI using microscopic spatial modulation of magnetization
EP0291387A1 (fr) Procédé de mesure des flux dans une expérimentation de résonance magnétique nucléaire
EP0316312B1 (fr) Procede d'excitation radiofrequence dans une experimentation de rmn
FR2630825A1 (fr) Procede et dispositif pour la correction d'une degradation de phase provoquee par des courants parasites
JP4679968B2 (ja) 磁気共鳴イメージング装置
EP0297935B1 (fr) Procédé de prise en compte, dans une image, des mouvements d'un objet
FR2586296A1 (fr) Procede de modulation de l'effet de la vitesse des parties mobiles d'un corps dans une mesure de densite par resonance magnetique nucleaire, et mise en oeuvre du procede pour en deduire la vitesse des parties mobiles concernees
JP3173612B2 (ja) 磁気共鳴分光法および映像法における、またはそれに関連する改良
EP0216657B1 (fr) 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
EP0290554A1 (fr) Procede de representation des parties en mouvement dans un corps par experimentation de resonance magnetique nucleaire
FR2838195A1 (fr) Procede d'imagerie rapide par resonnance magnetique nucleaire
EP0307433B1 (fr) Procede d'imagerie rapide par resonance magnetique nucleaire
JPH09164125A (ja) １回のスキャンで水画像と脂肪画像とを分離するのに必要な情報を収集するための方法
Sýkora K-space formulation of MRI
FR2584188A1 (fr) Procede d'imagerie par resonance magnetique nucleaire
FR2623907A1 (fr) Procede d'acquisition rapide d'images tridimensionnelles par experimentation rmn
FR2649213A1 (fr) Procede d'imagerie tridimensionnelle par resonance magnetique nucleaire

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): DE FR GB IT NL

17P Request for examination filed

Effective date: 19870206

17Q First examination report despatched

Effective date: 19890308

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: GENERAL ELECTRIC CGR SA

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): DE FR GB IT NL

ITF It: translation for a ep patent filed
GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)
REF Corresponds to:

Ref document number: 3666632

Country of ref document: DE

Date of ref document: 19891130

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed
PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 19920526

Year of fee payment: 7

Ref country code: DE

Payment date: 19920526

Year of fee payment: 7

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 19920527

Year of fee payment: 7

ITTA It: last paid annual fee
PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 19920630

Year of fee payment: 7

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Effective date: 19930627

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Effective date: 19940101

NLV4 Nl: lapsed or anulled due to non-payment of the annual fee
GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 19930627

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FR

Effective date: 19940228

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Effective date: 19940301

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED.

Effective date: 20050627