Magnetic resonance method and device for acquiring short echo time images

Abstract: A magnetic resonance imaging method for high quality short echo time imaging is suggested using FID-acquired-echoes. Gradient recalled echoes are excited in two short echo time acquisition sequences having opposite reading gradients. Half-echoes for corresponding phase encoding steps are concatenated to reconstructed echoes. An image is reconstructed from the reconstructed echoes using a Fourier transform reconstuction method. With the method artefact free high quality images with echo times below 4 msec can be obtained, even under unfavorable conditions, such as fast blood flow.

Description: The invention relates to a magnetic resonance imaging method, in which an object is situated in a stationary homogeneous magnetic field and is subjected to an acquisition sequence for acquiring a magnetic resonance signal from a region of the object, the acquisition sequence comprising an rf-pulse for exciting atomic spins in the object, a phase encoding gradient, and a reading gradient which is reversed for acquiring the magnetic resonance signal, the acquisition sequence being repeated a number of times with a different value of the time integral of the phase encoding gradient. The invention also relates to a magnetic resonance device comprising means for generating a stationary homogeneous magnetic field, means for exciting rf-pulses, means for generating magnetic field gradients, means for detecting and sampling magnetic resonance signals obtained with the aforementioned means in an object which is situated in the homogeneous field, control means for controlling said means for exciting rf-pulses, said means for generating magnetic field gradients and said means for detecting and sampling, image reconstruction means, and programmed means programmed such that an acquisition sequence for acquiring a gradient recalled echo with a reading gradient is applied to the object by means of said means for exciting rf-pulses, said means for generating magnetic field gradients, said means for detecting and sampling and said control means. A method and device of this kind are known from the Handbook, "Practical NMR Imaging", M.A. Foster and J.M.S. Hutchinson, 1987, IRL Press, pp. 15-22, and pp. 28-38. A Fourier transform sequence is given at page 19 and a blockdiagram of an NMR imaging device is shown at page 29 of said Handbook. In the method the body is irradiated with an rf-pulse, e.g. a 90 o pulse, in the presence of a slice selection gradient, thus flipping atomic spins, such as nuclear spins or with adaptations electron spins, of a slice of the body in a rotating frame rotating with the Larmorfrequency by an angle of 90°, resulting in an FID-signal (free induction decay). Then the FID is succesively dephased by a reading gradient having a first polarity and rephased by reversing the reading gradient. During rephasing a socalled gradient recalled echo occurs. Twodimensional coding is achieved by applying a variable amplitude phase encoding gradient and by sampling the echo signal. Fourier reconstruction gives a spatially resolved spin distribution of the slice that can be displayed on a monitor as a greyscale image e.g. Extension to 3D-imaging is wellknown and straightforward, so is a modification of the sequence to a socalled projection reconstruction sequence. In the sequence the maximum of the echo occurs at a point of time t e from the centre of the pulse at t o . The time difference t e - t o is the socalled echo time, which may not be chosen too small in case a complete echo is to be sampled. A too small echo time causes echo distortion and eventually image artefacts. In case short echo times are desired, e.g. for Sodium-imaging, for less flowartefacts, or for getting a different image contast, a distorted echo could only be sampled right to the echo maximum or so (asymmetric sampling). Ideally the image would be real, and displaying the real part of the image would suffice. But due to inter alia B₀-inhomogeneities, susceptibility, eddy currents and flow, the image is not real. Thus, direct real-part reconstruction or, equivalently, data mirroring in k-space often leads to inacceptable image artefacts. Correction methods are known for correcting some of these problems e.g. from the European Patent Application No. 0 250 050, giving an iterative correction algorithm. This algorithm needs more than half of the k-space (e.g. 5/8-th of an echo is sampled), thus restricting echo time reduction. Such correction methods tend to be artefact-sensitive if 50% of the data in k-space is missing, or if the phase varies rapidly over the image in the case of blood flow e.g. Then it would be expedient to acquire data over the whole k-space. It is an object of the invention to provide an NMR-imaging method having short echo times and not leading to unacceptable artefacts due to rapid phase variations e.g. To this end a method according to the invention is characterized in that the object is subjected to further acquisition sequences differing from said acquisition sequences in having a reading gradient of opposite polarity to said reading gradient, in that the magnetic resonance signals from the acquisition sequences and the further acquisition sequences are sampled starting sampling at least approximately from maxima of the resonance signals giving a first and a second set of half-echo samples, in that the half-echo samples from the first and second set corresponding to the same value of said time integral are combined to form a reconstructed full echo, and in that an image is reconstructed from the reconstructed full echoes. Very short echo times without disturbing flowartefacts have been obtained, e.g. in a 1.5 T system between 3 and 4.5 msec with usual imaging parameters. Samples in the centre of the k-space are measured twice and an average is formed when combining "halfechoes" to full "pseudo-echoes". As the two readout gradients have opposite signs, so do the flow-induced phase shifts caused by movements in the readout direction. For the k-space-centre samples, artefacts are therefore cancelled ruing reconstruction. This effect is enhanced by the low sensitivity to flow-induced phase shifts of a short T E sequence. An embodiment of a method according to the invention is characterized in that before forming a reconstructed full echo a first and a second sampling shift are applied to the half-echo samples of the first and second set respectively, the first and second sampling shift being obtained from samples from resonance signals acquired with a given value of said time integral by determining shifts of respective first samples to signal maxima. It is thus achieved that unacceptable misalignment of half-echoes is corrected before combining the half-echoes to a reconstructed, complete, echo. The sampling shifts are determined for a zero phase encoding e.g. and may appear to be two or three sampling intervals in practice. For other phase encoding steps the determined sampling shifts are applied before combining half-echoes. The features and advantages of the method and device according to the invention will be more clearly apparent from the following detailed description, in conjunction with the accompanying drawing, wherein Fig. 1 shows, diagrammatically, a magnetic resonance device according to the invention, Fig. 2 shows a pulse sequence according to the invention, Fig. 3A, 3B and 3C show combination of half echoes to a reconstructed echo in detail, and Fig. 4 shows phase shifts induced by flowing spins in readout direction. In Fig. 1 there is shown, diagrammatically, a magnetic resonance device 1 according to the invention comprising means 2 for generating a stationary homogeneous magnetic field B₀, means 3 for exciting rf-pulses to be irradiated to an object 4, modulating means 5 for modulating the rf-pulses, detecting means 6 for detecting excited magnetic resonance signals, sampling means 7 for sampling the magnetic resonance signals and control means 8 for controlling said means 3, 5, 6 and 7. The device 1 further comprises Fourier transform (or filtered back-projection) means 9 for reconstructing an image from the sampled magnetic resonance signals, means 10 for generating magnetic field gradients which are superimposed upon the stationary magnetic field. In further detail the means 3 comprise the means 5, a power amplifier 11, a directional coupler 12 and a coil 13. The means 5 comprise a modulator 14 and an oscillator 15. The means 3 will excite nuclear spins in the object 4 when transmitting rf-pulses and the means 6 will receive magnetic resonance signals in response thereupon. The demodulating means comprise a quadrature demodulator and the sampling means 7 two A/D-converters giving quadrature detected samples of respective sin and cos components of the magnetic resonance signal. The means 10 comprise a power supply 16 to be controlled by the control means 8, for separately controlling gradient magnet coils 17, 18 and 19 for respective generation of magnetic field gradients G x , G y and G z . The gradient magnet coils may be arranged such that the field direction of the magnetic field gradient coincides with the direction of the stationary homogeneous magnetic field and that the directions of the gradients are mutually perpendicular, indicated with mutually perpendicular axes x, y and z. If the device 1 is put into operation and the object 4 is placed inside the magnet coils 2 a small excess of nuclear spins (of nuclei having a magnetic moment) will be directed into the direction of the stationary homogeneous magnetic field, which can be considered macroscopically as a magnetization M, an equilibrium magnetization. For a more detailed description of an MRI device and a general description of the principles of MRI referred is to said Handbook, "Practical NMR Imaging". In Fig. 2 there is shown a pulse sequence according to the invention as a function of time t. At the instant t=t₀ an rf-pulse p is excited with said means 3 and irradiated to the object 4. For selecting a part of the object 4 (a slice perpendicular to the z-axis e.g.) a selection gradient G z is applied during the rf-pulse p by means of the gradient magnet coil 19. The rf-pulse p selectively rotates a magnetization (90 o or less). A magnetization in the xy-plane (transverse magnetization) gives an FID-signal (not shown). A phase encoding gradient G y is applied after the rf-pulse p by means of the gradient magnet coil 18 for phase encoding nuclear spins in the part of the object. During the gradient G y also a reading gradient G x is applied by means of the gradient coil 17 for dephasing the spins. Thereafter the reading gradient G x is reversed giving a recalled FID at the instant t=t e , an echo e. The echo e is quadrature detected and sampled by means of two A/D converters, so that real and imaginary signal components are acquired. The acquisition sequence is repeated for different values of the time integral of the phase encoding gradient G y . The samples are stored, building up a first set of samples. A further acquisition sequence is applied being identical to the described sequence but for the gradient G x . A gradient G x ' therein is opposite to the gradient G x . The samples acquired with the further acquisition sequence are stored and build up a second set of samples. Samples from the first and second set are combined to form a reconstructed echo. By means of the Fourier transform means 9 a two dimensional Fourier transform is carried out on the combined sets for getting a spin-density distribution. The spin-density distribution is displayed as a greyscale image e.g. on a monitor (not shown). In Fig. 3A, 3B and 3C combination of half-echoes to a reconstructed echo is shown in detail. Fig. 3A shows real Re and imaginary Im parts of half-echoes e1 and e2 acquired with said acquisition sequences as a function of time t. Before combining the half-echoes e1 and e2 to a reconstructed echo a first and a second sampling shift ssh1 and ssh2 are determined with respect to actual first samples ts1 and ts2. The sampling shifts ssh1 and ssh2 are determined from maximum moduli of samples around the signal maxima M1 and M2 respectively. Typically the sampling shifts ssh1 and ssh2 are in the order of 2-3 sampling intervals si. The sampling shifts ssh1 and ssh2, determined for a zero phase encoding e.g., are applied for all other phase encoding steps. The first order phase correction prior to concatenating the half-echoes (FID-acquired-echoes) to a reconstructed echo ensures that symmetric reconstructed echoes are obtained, otherwise artefacts for lower spatial frequencies would occur. In fig. 3B wherein k-space (k x k y ) is shown it is achieved that k-space is scanned symmetrically, shown for G y =0 for a right scan rs and a left scan ls. In practical situations omitting phase correction would cause asymmetrical k-space scanning. In Fig. 3C a simular situation is shown when applying projection reconstruction sequences. With the method according to the invention very short echo times can be achieved. It should be noted that were there no phase errors (which does not hold in practice, generally) short echo times could also be achieved from right half k-space data e.g., but the modulus reconstruction could not be applied. Phase errors would lead to unacceptable image artefacts in practical situations when using only right half k-space data. The method according to the invention has been implemented on a whole body 1.5 T system. For a field of view of 250 mm, a 256x256 pixel matrix (giving a 1 mm image resolution), and a slice thickness of 5 mm, in a 2D-FT or 3DFT fast imaging sequence echo times between 3.0 and 4.5 msec were obtained. For small sequence repetition times (e.g. <100 msec) an opposite G y gradient and for spurious FID suppression a constant dephasing gradient G x (repetition times <200 msec) were applied after signal sampling. E.g. images have been made of a human head using headcoils or heart images using a whole body coil. In a multiple slice version together with gating and a sequence repetition time of 30 msec a heart movie was made with an image frequency of 32 images per second for the same slice. Flow images can be made when displaying a phase image instead of a modulus image. It should be noted that the order of acquisition of the signals for the first and second set can be chosen freely but it is advantageous to acquire signals with opposite reading gradients for the same value of the time integral of the phase encoding gradient consecutively. It should further be noted that the acquisition sequences may be modified, the invention being applicable to many variants of the described sequence, so the sequences may constitute a multiple-echo sequence e.g. Fig. 4 shows phase shifts induced by flowing spins in readout direction for conventional sequences such as fast imaging (FE) and according to present invention (FAcE) as a function of time t with respect to M=0 of the first moment M of the readout gradient during the sampling period and with respect to the echo maximum at t=t₀. As in present invention the two readout gradients have opposite signs, so do the flow-induced phase shifts caused by movements in the readout direction. For the k-space-centre samples, artefacts are therefore cancelled during reconstruction. Experiments with flow compensated fast image sequences and sequences according to present invention show superior artefact suppression of present invention. Due to this high flow-artefact suppression, in vivo images display very fine branching vessels in the brain, and also lung vessels.

Claims: Magnetic resonance imaging method, in which an object (4) is situated in a stationary homogeneous magnetic field and is subjected to an acquisition sequence for acquiring a magnetic resonance signal (e) from a region of the object, the acquisition sequence comprising an rf-pulse (p) for exciting atomic spins in the object, a phase encoding gradient (G y ), and a reading gradient (G x ) which is reversed for acquiring the magnetic resonance signal (e), the acquisition sequence being repeated a number of times with a different value of the time integral of the phase encoding gradient (G y ), characterized in that the object is subjected to further acquisition sequences differing from said acquisition sequences in having a reading gradient (G x ') of opposite polarity to said reading gradient (G x ), in that the magnetic resonance signals (e) from the acquisition sequences and the further acquisition sequences are sampled starting sampling at least approximately from maxima of the resonance signals giving a first and a second set of half-echo samples (e1,e2), in that the half-echo samples (e1,e2) of the first and second set corresponding to the same value of said time integral are combined to form a reconstructed full echo, and in that an image is reconstructed from the reconstructed full echoes. Magnetic resonance imaging method according to claim 1, characterized in that before forming a reconstructed full echo a first and a second sampling shift (ssh1,ssh2) are applied to the half-echo samples (e1,e2) of the first and the second set respectively, the first and the second sampling shift (ssh1,ssh2) being obtained from samples from resonance signals acquired with a given value of said time integral by determining shifts of respective first samples to signal maxima. Magnetic resonance imaging method according to claim 2, characterized in that the shifts (ssh1, ssh2) to maxima are determined from moduli of samples around the maxima. Magnetic resonance imaging method according to claim 1,2 or 3, characterized in that the image is reconstructed as a modulus image. Magnetic resonance imaging method according to Claim 1, 2, 3 or 4, characterized in that said acquisition sequences are modified to projection reconstruction sequences. Magnetic resonance imaging method according to one of the preceding Claims used for flow imaging, especially for flow in the direction of the reading gradient. Magnetic resonance device (1) comprising means (2) for generating a stationary homogeneous magnetic field, means (3) for exciting rf-pulses, means (10) for generating magnetic field gradients, means (6, 7) for detecting and sampling magnetic resonance signals obtained with the aforementioned means in an object (4) which is situated in the homogeneous field, control means (8) for controlling said means (3) for exciting rf-pulses, said means (10) for generating magnetic field gradients and said means (6, 7) for detecting and sampling, image reconstruction means (9), and programmed means programmed such that an acquisition sequence for acquiring a gradient recalled echo with a reading gradient is applied to the object (4) by means of said means (3) for exciting rf-pulses, said means (10) for generating magnetic field gradients, said means (6, 7) for detecting and sampling and said control means (8), characterized in that the device comprises further programmed means programmed such that a further acquisition sequence for acquiring a further gradient recalled echo with a reading gradient of opposite polarity to said reading gradient is applied to the object, that half-echo samples (e1,e2) at least approximately starting from maxima of the resonance signals of the respective sequences are obtained, that sampling shifts (ssh1,ssh2) are applied to the half-echo samples (e1,e2), that the half-echo samples (e1,e2) are combined to a reconstructed full echo, and that an image is reconstructed from the reconstructed full echoes.

Inventors: BOESIGER PETER, SCHEIDEGGER MARKUS, BOESIGER, PETER, SCHEIDEGGER, MARKUS, Bösiger, Peter
Applicants: PHILIPS ELECTRONICS NV, PHILIPS ELECTRONICS N.V.
Application year: 1990
Main classification: G01R33
IPC codes: G01R33