METHOD OF AND DEVICE FOR THE PHASE CORRECTION OF MR INVERSION RECOVERY IMAGES

Abstract: During known MRI measurements, for example spin echo and inversion recovery measurements, phase errors occur in pixels, inter alia due to eddy currents caused by the application or interruption of magnetic field gradients during a measurement. There is proposed an inversion recovery method combined with a method where the magnetization is not inverted in advance in order to determine for each pixel the phase error in an image obtained by means of an inversion recovery method, so that inversion recovery images can be obtained which do not contain phase errors.

Description: The invention relates to a method of determining a nuclear magnetization distribution in a region of a body which is situated in a steady, uniform magnetic field, spind resonance signals being generated in the region of the body in two sequences by high-frequency excitation pulses and being measured in the presence of magnetic field gradients superposed on the uniform magnetic field, the excitation pulse in one sequence being preceded by a high-frequency inversion pulse which inverts the nuclear magnetization, using a transformation a first image of complex image values being formed from said one sequence whilst a second image of complex image values is formed from the other sequence, each pixel in the first and the second image being represented by an image amplitude and an image phase. The invention also relates to a device for deter­ mining a nuclear magnetization distribution in a region of a body, including: means for generating a steady, uniform magnetic field, means for generating high-frequency excitation pulses, means for generating magnetic field gradients which are superposed on the uniform magnetic field, means for measuring spin resonance signals, processing means for processing spin resonance signals measured, control means for controlling the above means, said control means being suitable for the measurement of different sequences, said processing means being suitable for combining measure­ ment results from two different sequences and for deter­ mining images from the sequences. Such a method and device are known from EP 0.145.276. Devices for determining a nuclear magnetization distribution in a region of a body and the principles on which the operation of such devices is based are known, for example from the article "Proton NMR tomography" in Philips Technical Review, Volume 41, 1983/84, No. 3, pages 73-88. Reference is made to this article for the description of their construction and principles. In a method described in EP 0.145.276, a body to be examined is exposed to a steady, uniform magnetic field Bo during a so-called saturation recovery measure­ ment wherein a selective 90° excitation pulse is followed by a non-selective 180° pulse in order to obtain an echo signal, said steady, uniform magnetic field Bo having a direction which coincides with, for example the z-axis of a cartesian coordinate system (x,y,z). The steady magnetic field Bo causes a slight polarization of the spins present in the body and enables spins to perform a precessional motion about the direction of the magnetic field Bo. After the establishment of the magnetic field Bo, a magnetic field gradient which acts as a selection gradient is applied and at the same time a 90° rf pulse is generated which rotates the magnetization of the nuclei present in a selected slice through an angle of 90°. After termination of the 90° pulse, the spins will perform a precessional motion about the field direction of the magnetic field Bo, thus generating a resonance signal (FID signal). After the 90° pulse, there are simul­ taneously applied field gradients G z , G x and G y whose field direction coincides with that of the magnetic field Bo and whose gradient directions extend in the z, the x, and the y direction, respectively. The field gradient G z then opposes the G z applied during the 90° pulse and serves for rephasing the excited spins in a slice perpendicular to the z direction. The field gradients G x and G y serve for encoding in the x direction and the y direction, respectively. After termination of the three above field gradients and after application of a 180° echo pulse, a field gradient G x is applied and an echo resonance signal of the original FID signal is sampled, after which an image is reconstructed after Fourier transformation. In the cited European Specifi­ cation EO 0.145.276 it is noted that the phase of the pixels generated by means of this so-called saturation recovery measurement can be disturbed by inter alia eddy currents which arise in conductive parts of an MRI apparatus due to the application or termination of mag­ netic field gradients during a measurement. In order to eliminate or at least mitigate phase disturbances in pixels generated by means of a saturation recovery measure­ ment, the cited Patent Specification EP 0.145.276 pro­ poses an extension of the method in the form of an in­ version recovery measurement, including a non-selective 180° pulse in order to obtain an echo signal. The pulse sequence thereof has the same appearance as that of a saturation recovery measurement, be it that for an in­ version recovery measurement the pulse sequence is in­ creased to include a high-frequency electromagnetic 180° pulse which is generated some time before the high-­ frequency electromagnetic 90° pulse. It will be apparent that the time sequence in which the saturation recovery measurement and an inversion recovery measurement are performed is irrelevant. By determining the quotient per pixel of the images obtained from an inversion recovery measurement and a saturation recovery measurement, res­ pectively, ideally an image without phase errors will be obtained, provided that the two type of measurement have both been performed under the same experimental conditions so that the phase errors cancel one another. It is a drawback of the known, extended method that said quotient in the cited Patent Specification EP 0.145.276 contains an imaginary term which represents the difference between the phase errors in the two sequences and which is not always negligibly small in practice, so that corrected pixels will still contain phase errors which give rise to annoying differences in intensity in a corrected image. It is a further draw­ back that the inversion recovery measurement does not produce a pure inversion recovery image, but rather, in conjunction with the saturation recovery measurement, a proportional image. It is the object of the invention to provide a method which enables the determination of pixels of an inversion recovery image which do not contain phase errors, so that an image is obtained which does not contain annoying intensity differences. To this end, a method in accordance with the invention is characterized in that for each corresponding pixel the value of a sign function is determined from the phase difference between the pixels of the first image and the second image. If only the complex image values of said first image were used to form an amplitude image by taking the modulus of each complex image value, the sign information of the first image would be lost. When each image value in said amplitude image is multiplied by the sign function whose value depends on said phase difference and when a suitable sign function is chosen, the sign function of the first image will not be lost. One version of a method in accordance with the invention is characterized in that the sign function has the value one if - π /2 < the phase different < π /2, and otherwise has the value minus one. When the phase er­ rors are of the same order of magnitude, this sign func­ tion which discriminates equally as regards a positive or a negative sign, will be a suitable choice. A different variation of the sign function could be considered when the phase errors deviate substantially. It is to be noted that phase errors may be due not only to systematic causes, such as for example eddy currents, but also to non-systematic causes, for example spins which are subject to flow. A further version of a method in accordance with the invention where the variation of the magnetic field gradients is identical in the two sequences for each pair or resonance signals, is characterized in that in the other sequence a dummy gradient is applied at the location which corresponds to that at which the high-frequency inversion pulse occurs in said one sequence. This dummy gradient, not having a function per se, merely serves to make the gradient identical in the sequence. This dummy gradient and all other gradients in the other sequence should have the same amplitude and duration and should be switched at instants which correspond to those in the one sequence. It is thus achieved that the phase errors of systematic origin will be at least sub­ stantially equal. A further version of a method in accordance with the invention is characterized in that the other sequence is generated with a different resolution and that by interpolation and transformation an image is obtained therefrom having the same resolution as an image from said one sequence. By generating the other sequence with a different, preferably lower resolution, it is achieved that the total experimentation time is substantially reduced, which is favourable both as regards the rf dose for, for example a patient to be examined as well as regards rf load for a high-frequency pulse transmitter. A comparatively short experimentation time also reduces the risk of image artefacts due to motions of, for example, a patient to be examined. A further version of a method in accordance with the invention, where a third image is formed from the first and the second complex image, is characterized in that each pixel of the third image is derived from the product of the image amplitude of the corresponding pixel of the first image and the associated value of the sign function. An image thus obtained represents complex image values of the first image wherefrom phase errors have been removed as well as possible, even in the case of substantial phase errors. A further version of a method in accordance with the invention, where for each pixel there if deter­ mined a ratio which depends on the longitudinal relaxation time T₁, is characterized in that the ratio is the quotient of the image amplitude of the corresponding pixels of the first and the second image, multiplied by the value of the sign function for these pixels. Inter alia from EP 0.140.276 it is known that the local longitudinal relaxation time T₁ can be cal­ culated from the quotient per pixel of the images ob­ tained from an inversion recovery measurement and a saturation recovery measurement, respectively. As has already been explained, this quotient will contain an imaginary term which cannot be ignored in practice. Due to phase errors, this will lead to incorrect deter­ mination of the local relaxation time T₁. The said ratio offers a better basis for the determination of T₁ per pixel, and hence for the formation of a T₁ image. A device in accordance with the invention is characterized in that for corresponding pixels the pro­ cessing means determine the value of a sign function from the phase difference between the pixels of the first and the second image. The invention will be described in detail hereinafter with reference to pulse sequences of an in­ version recovery measurement with spin-echo reading and a spin-echo measurement as shown in the figures 1 and 2 respectively. The image intensity per pixel of an image obtained by means of NMR techniques generally depends on the proton density ρ locally present in an excited layer, on the longitudinal relaxation time T₁ and on the transverse relaxation time T₂; these are important para­ meters which influence the nuclear magnetization dis­ tribution in the layer. Inversion recovery measurements are suitable for obtaining images whose intensity depends mainly on the density ρ locally present in the excited layer and on the relaxation time T₁. Spin-echo measure­ ments are usually performed for obtaining an image inten­ sity which depends mainly on the density ρ and the relaxation time T₂. In practice, inversion recovery measure­ ments and spin-echo measurements are performed during one experiment in order to determine the nuclear magnetization distribution by means of images from an exited layer which are mainly ρ - T₁ or ρ = T₂ dependent. One problem of NMR measurements in general is the occurrence of eddy currents. These eddy currents themselves cause magnetic field gradients which intro­ duce phase errors in an image to be obtained from a measure­ ment. This is because a given number of nuclei which have the same precessional phase and the same precessional frequency will make a given signal contribution to the ultimate signal. Because the magnetic field gradients caused by eddy currents will also influence the precessional phase and the precessional frequency of that given number of nuclei, the number of nuclei generating a given signal contribution will change. As a result, the intensity of that signal contribution will also change. Figure 1 shows a normal inversion recovery pulse sequence with spin-echo reading. A clockwise ro­ tating orthonormal coordinate system will be considered in which a static, uniform magnetic field is directed along the z-axis, a high-frequency electromagnetic field is directed perpendicularly to the direction of the uni­ form magnetic field, and the field directions of gradient fields coincide with the field direction of the uniform field, whilst the gradient directions extend perpendicularly with respect to one another. During an interval 1, bounded by the instants t = t a and t = t b , a magnetic field gra­ dient G is applied and at the same time a high-frequency electromagnetic 180° pulse is generated. As a result, the spins will be excited in a plane Z = Z o , which means that the magnetization generated by the spins in the plane Z = Z o encloses an angle of 180° with respect to the positive z-axis, i.e. is directed along the negative z-axis. Because the 180° pulse has a given bandwidth, spins will be excited not only in the plane Z = Z o , but also in the immediate vicinity thereof. By variation of the frequency content of the 180° pulse, spins in another plane, parallel to the plane Z = Z o , can be excited. After the instant t = t b , the spins will reorient in the direction of the static, uniform magnetic field, because as from that instant they will no longer be subject to a high-frequency electromagnetic pulse. This relaxation process continues until the magnetization is directed in the original positive z-direction again. the magnetization will then successively decrease along the negative z-­ axis, become zero and increase along the positive z-axis until the original value is reached again, as has already been stated. During an interval 3, bounded by the instants t = t c and t = d d , a high-frequency electromagnetic 90° pulse is generated, the magnetic field gradient G being present, in order to cause a precessional motion of the magnetization in the excited layer at an angle of 90° about the z-axis. During an interval 4, bounded by the instants t = t d and t = t e , three magnetic field gradients G, G and G are applied. In order to save time, these three field gradients are applied during the same time interval. However, because no high-frequency electromag­ netic pulses are generated therein, the effects of the field gradients can be separately considered. The nega­ tive gradient G rephases the spins in the z-direction so that the spins which have been selectively excited in and in the immediate vicinity of the plane Z =Z o are no longer subject to phase discrimination in the z-direc­ tion. The magnetic field gradients G and G serve for deliberate dephasing of the spins in the x-direction and the y-direction, respectively. Encoding is thus re­ alized in both directions. Between the instants t = t f and t = t h there is generated a resonance signal which depends on the magnetization locally present in the z-­ direction immediately before the instant t = t c . The magnitude of this magnetization eepends on the longitudinal relaxation time T l . When an excited layer of, for example a patient contains different types of tissues having different T l times, the local magnetizations may be either in phase or 180° out of phase. This is because, when the duration T = (t b . t c ) is smaller than T l .1n₂, so that a local magnetization is still directed along the negative z-axis prior to the influencing of the magne­ tization by a high-frequency electromagnetic 90° pulse, the ultimate magnetization will be 180° out of phase with respect to the local magnetization generated by nucle in, for example a tissue for which the duration t exceeds T₁.1n₂. The local magnetization will then be directed along the positive z-axis at the instant t = t c . The previously mentioned magnetic field gradients caused by eddy currents can disturb the phase or frequency com­ ponents in resonance signals during inversion recovery measurements. Figure 2 shows a spin-echo pulse sequence. During interval 1 there is applied a dummy gradient G in order to obtain the same G z gradient variation as during the inversion recovery sequence, so that at least substan­ tially equal eddy currents will arise. Therefore, the phase errors caused by eddy currents will be at least substantially equal in both sequences. The high-frequency electromagnetic 90° pulse generated during the interval 3, bounded by t = t c and t = t d , rotates the magnetiza­ tion formed by the spins locally present in the layer through an angle of 90° with respect to the z-axis. The high-frequency electromagnetic 180° pulse generated during the interval 5, bounded by the instants t = t e and t = t f , creates an echo resonance signal during the interval (t f , t h ). The functions of the magnetic field gradients G, G z , G and G are the same as stated in the description of the inversion recovery pulse sequence of Figure 1. Even though, for example, different types of tissue having different T₁ relaxation times may again be present in an excited layer, the local magnetizations ultimately obtained will ideally be in phase. This is because magnetizations which are directed along the negative z-axis will not occur during a measurement. Due to the described eddy currents, disturbances will also occur in an image obtained during spin-echo measurements. As will be demonstrated hereinafter, the phase errors in an inversion recovery image which could result in pixels having an incorrect intensity can be reduced by com­ bining an inversion recovery measurement and a spin-echo measurement. After double phase-sensitive detection of the resonance signals of the inversion recovery measure­ ment, the two groups of signals derived therefrom are sampled after filtering, after which a complex image is reconstructed from the two groups of sampling values by way of 2D Fourier transformation. This operation is also performed on the resonance signals of the spin-echo measurement. The value of an inversion recovery (IR) pixel and of a spin-echo (SE) pixel can be expressed in a com­ plex notation: S IR = mod (S IR )e i⊘S IR (1) S SE = mos (S SE ) e i⊘S SE (2) wherein S IR is the value of a complex IR pixel, mod (S IR ) is the modulus, and ⊘ S is the phase of the complex IR pixel. A corresponding notation is used for the spin-­ echo resonance pixel. Ideally, that is to say without the phase errors caused by eddy currents and without other phase errors, for example as caused by inhomogeneities of the static, uniform magnetic field and by imperfect high-frequency electromagnetic 90° pulses or 180° pulses, the following would hold good: ⊘ S = 0 so that S SE = mod (S SE ) for ⊘ S = 0 (4) When the phase errors caused by eddy currents and the other phase errors are taken into account in a term Δ ⊘, the expressions (1) and (2) will be: The term Δ ⊘ S can now be simply determined for each pixel according to: wherein ReS SE and Im S SE are the real and the imaginary part, respectively, of the complex pixel (2). Assuming that the phase errors of the complex pixels S SE and S IR exhibit the same behaviour, so Δ ⊘ S ≈ Δ ⊘ S (8) the phase errors in inversion recovery pixels can be reduced. This condition can indeed be substantially satisfied by choosing the experimental situation to be such that the spins during an inversion recovery measurement or, for example during a spin-echo measurement, are subject to the same effect on the phase error by the gradients. The measurements as shown in the Figures 1 and 2 are relevant examples. the modulus image of the IR pixel in accordance with formula (1) is the same as the modulus image of the IR pixel in accordance with formule (5). For correction it will merely be necessary to determine the + or - sign of themodulus in accordance with formule (3). It has been found that the modulus can be suitably estimated to obtain a + sign if - π /2 < ⊘ S + Δ ⊘ S - Δ ⊘ S < π /2 and that it obtains a - sign if π /2 < ⊘ S + Δ ⊘ S - Δ ⊘ S < 3π /2 (9) Therein, ⊘ S + Δ ⊘ S is the measured phase per pixel of the IR measurement and Δ ⊘ S is the measured phase per pixel of the SE measurement. Other estimation inter­ vals may also be chosen. In addition to providing the modulus with a + or - sign, for example the modulus may also be reduced in the vicinity of the boundaries of the estimation intervals, which implies that the sign function will have an absolute value smaller than one in the vicinity of the boundaries of the estimation intervals. Another possibility consists in the use of more advanced statistical methods for determining the + or - sign in the vicinity of the boundaries. The invention can also be used for determining the longitudinal relaxation time T₁ per pixel. From Journal of Magn. Res. in Med. 2, 136-158 (1985) it is known to derive the relacation time T₁ from the ratio of the real portions Re (S IR ) and Re (S SE ) of the pixels S IR and S SE , respectively: As appears from equation (10), this determination of the relaxation time T₁ is also subject to phase errors. The T₁ relaxation time can be determined without phase errors when in accordance with the invention the time T₁ is derived by means of: wherein the sign function generates a + sign or a - sign, respectively, in accordance with the criterion of formule (9) or another criterion. The longitudinal relaxation time T₁ per pixel can subsequently be determined by using a value R obtained by means of the equation (11) for looking up an associated relaxation time T₁ in com­ puter look-up tables in which the theoretical ratio R is stated together with the associated time T₁. Using the R obtained from equation (11), T₁ can also be solved from the theoretical ratio R by means of numerical zero­ point determination methods. The theoretical dependency of T₁ on the ratio R can be found by solution of Bloch's equations. It is to be noted that after the interval 2 in both sequences any other method can be used in order to derive the re­ sonance signals or to generate echo signals, including manipulation of gradients such that image reconstruction is also possible using, for example a projection recon­ struction method.

Claims: 1. A method of determining a nuclear magnetization distribution in a region of a body which is situated in a steady, uniform magnetic field, spin resonance signals being generated in the region of the body in two sequences by high-frequency excitation pulses and being measured in the presence of magnetic field gradients superposed on the uniform magnetic field, the excitation pulse in one sequence being preceded by a high-frequency inver­ sion pulse which inverts the nuclear magnetization, using a transformation a first image of complex image values being formed from said one sequence whilst a second image of complex image values is formed from the other sequence, each pixel in the first and the second image being represented by an image amplitude and an image phase, characterized in that for each corresponding pixel and value of a sign function is determined from the phase difference between the pixels of the first image and the second image. 2. A method as claimed in Claim 1, characterized in that the sign function has the value one if - π / 2 < the phase difference < π /2, and otherwise has the value minus one. 3. A method as claimed in Claim 1 or 2, where the variation of the magnetic field gradients is identical in the two sequences for each pair of resonance signals, characterized in that in the other sequence a dummy gradient is applied at the location which corresponds to that at which the high-frequency inversion pulse occurs in said one sequence. 4. A method as claimed in Claim 1, 2 or 3, charac­ terized in that the other sequenee is generated with a different resolution and that by interpolation and trans­ formation an image is obtained therefrom having the same resolution as an image from said one sequence. 5. A method as claimed in any one of the preceding Claims, where a third image is formed from the first and the second complex image, characterized in that each pixel of the third image is derived from the product of the image amplitude of the corresponding pixel of the first image and the associated value of the sign function. 6. A method as claimed in any one of the Claims 1 to 4, where for each pixel there is determined a ratio which depends on the longitudinal relaxation time T₁, characterized in that the ratio is the quotient of the image amplitude of the corresponding pixels of the first and the second image, multiplied by the value of the sign function for these pixels. 7. A device for determining a nuclear magnetization distribution in a region of a body, including: means for generating a steady, uniform magnetic field, means for generating magnetic field gradients which are superposed on the uniform magnetic field, means for measuring spin resonance signals, processing means for processing spin resonance signals measured, control means for controlling the above means, said control means being suitable for the measurement of different sequences, said processing means being suitable for combining mea­ surement results from two different sequences and for determining images from the sequences, characterized in that for corresponding pixels the processing means determine the value of a sign function from the phase difference between the pixels of the first and the second image.

Inventors: GROEN, JOHANNES PETRUS, IN DEN KLEEF, JOSEPHUS JACOBUS EGIDIUS
Applicants: N.V. PHILIPS' GLOEILAMPENFABRIEKEN
Application year: 1986
Main classification: G01N24
IPC codes: A61B5, G01R33, A61B10