EP0182873A1 - Prise d'image de fourier nmr a partir d'echos multiples - Google Patents

Prise d'image de fourier nmr a partir d'echos multiples

Info

Publication number
EP0182873A1
EP0182873A1 EP19850902856 EP85902856A EP0182873A1 EP 0182873 A1 EP0182873 A1 EP 0182873A1 EP 19850902856 EP19850902856 EP 19850902856 EP 85902856 A EP85902856 A EP 85902856A EP 0182873 A1 EP0182873 A1 EP 0182873A1
Authority
EP
European Patent Office
Prior art keywords
echoes
applying
pulses
data set
pulse
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.)
Withdrawn
Application number
EP19850902856
Other languages
German (de)
English (en)
Inventor
Andrew A. Maudsley
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.)
Advanced NMR Systems Inc
Original Assignee
Advanced NMR Systems Inc
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
Application filed by Advanced NMR Systems Inc filed Critical Advanced NMR Systems Inc
Publication of EP0182873A1 publication Critical patent/EP0182873A1/fr
Withdrawn legal-status Critical Current

Links

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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0

Definitions

  • U.S. Patent No. 4,290,019, issued September 15, 1981 to J.M.S. Hutchison et al. discloses another technique for NMR Fourier imaging.
  • This technique includes the application of a. sequence of gradient fields and high frequency pulses to obtain a spin echo of the FID. After an echo is obtained, the atoms being imaged are permitted to relax before the sequence is repeated to obtain another echo. The echoes are processed to obtain an image, and the processing includes Fourier transformation.
  • the invention provides NMR Fourier imaging techniques with more sensitive data acquisition and increased image contrast between different types of tissue.
  • the invention is based on the discovery that sensitive data acquisition and increased image contrast can be obtained in NMR Fourier imaging by applying a pulse and gradient switching sequence which produces a series of echoes of each FID. These echoes can be combined into a data set which is then transformed into the image.
  • the method of the invention thus includes applying a sequence of stimulating signals, the pulse and gradient switching sequence, to the object to be imaged to produce a series of echoes.
  • the echoes are then combined to obtain a data set, and the data set is then transformed into an NMR image of the object.
  • the apparatus of the invention includes means for applying stimulating signals to an object, the applying means being operable for applying the pulse and gradient switching sequence which produces a plurality of echoes.
  • the apparatus includes combining means for combining the echoes to obtain a combined data set and transforming means for transforming the combined data set into an NMR image of the object.
  • the sequence of stimulating signals which is applied preferably includes an excitation pulse such as a 90° radio frequency (rf) pulse, which rotates the resonant nuclei in the object through 90°.
  • an excitation pulse such as a 90° radio frequency (rf) pulse
  • rf radio frequency
  • This pulse is then followed by a series of refocusing pulses which will be 180° rf pulses, each of which rotates the nuclei through 180°.
  • the phase error introduced by the pulses is compensated if the excitation and refocusing pulses follow a specific sequence.
  • the specific pulse and gradient switching sequence may be described in terms of the axes along which the pulses and gradients are applied.
  • the axis along which the constant magnetic field is applied is generally designated the Z axis, and each plane perpendicular to that axis will have a corresponding X axis and Y axis.
  • the X and Y axes are riot stationary, however, but are rotating in a frame of reference known as the "rotating frame", which rotates at the nuclear spin resonant frequency. Therefore, the pulses and spin magnetizations which are applied along the the X and Y axes in the following discussion are actually applied in the rotating frame of reference.
  • the envelope of the given signal is at the nuclear spin resonant frequency, with its phase determined according to the axis along which it is applied.
  • a signal applied along the Y axis will accordingly be 90o out of phase with a signal applied along the X axis.
  • the compensating sequence of refocusing pulses according to the invention can be a series of 180° rf pulses applied in groups of 8 pulses along the positive and negative X and Y axes in the following order: +X, +Y, +X, +Y, -X, -Y, -X, -Y.
  • the pulses are applied alternately along the X and Y axes, and the first four are along the positive axes, while the second four are along the negative axes. If a sequence of stimulating signals like that above is applied, a series of groups of four echoes will be produced.
  • each echo may be analyzed into an in phase magnetization component which has the same phase as the initial FID resulting from the excitation pulse and an out-of-phase component which is 90o out of phase from the initial FID.
  • the presence of the two components results from the application of phase encoding gradients, so that both components must be preserved in order to preserve the phase information. Therefore, the coherent addition of the echoes requires that both the in-phase components and the out-of-phase components may be added without cancellation to produce combined in-phase and out-of-phase components.
  • the process of coherent addition may be illustrated in terms of the pulse sequence described above.
  • the second of the echoes being added has an in-phase component which refocuses with inverted phase from that of the first echo, while its out-of-phase component refocuses with the same phase. Therefore, the first and second echoes may be coherently added after reversing the time axis of the second echo, multiplying it by -1 and taking its complex conjugate. As a result, each component of the second echo will be of the same sign as the corresponding component of the first, and the two echoes can be added together without cancelling.
  • alternate echoes may be time reversed and all of the even echoes may be added to form one data set while all of the odd echoes may be added to form another. In this case, the data set formed from the even echoes will contain information about spatial frequency from the positive phase encoding gradient, while the data set formed from the odd echoes will contain information about spatial frequency from the negative phase encoding gradient.
  • the echoes will be transformed, before or after coherent adding, to obtain a data set encoded according to spatial frequency.
  • This transform may be a Fourier transform or other similar transform, such as the FFT.
  • the terms "Fourier transform” and “transform” refer not only to the Fourier transform in its continuous and discrete forms, but also other similar transforms which depend on phase information, including digital approximations to the Fourier transform such as the FFT and the Hadamard transform. More generally, such a transform will be referred to herein as a "phase-dependent transform".
  • the data set encoded according to spatial frequency When the data set encoded according to spatial frequency has been completed, it may then be transformed to provide a spatial image of the object. This transformation may be in two or three dimensions, for example. Because the echoes are combined to form the data set, the resulting image has substantially better contrast than with the prior art techniques.
  • the method and apparatus of the present invention may reduce data acquisition time for NMR Fourier imaging. Rather than taking several separate FID measurements for signal averaging, as described above in relation to the prior art, the invention makes it possible to obtain several echoes from a single excitation. This is possible because the real spin-spin relaxation time, T 2, following an excitation is typically much longer than the effective relaxation time, T 2 *, of each echo. As a result, several echoes can be measured after each excitation and combined to obtain significantly greater sensitivity.
  • Fig. 1 is a cross sectional view of an arrangement of NMR gradient and rf coils around a subject.
  • Fig. 2A is a functional block diagram showing the basic operations of NMR Fourier imaging according to the invention.
  • Fig. 2B is a flow chart showing the basic steps of the pulse and gradient switching sequence of the invention.
  • Fig. 2C is a detailed flow chart showing one embodiment of coherent adding of echoes according to the invention.
  • Fig. 5A is a graphical illustration of the real and imaginary parts of an odd echo.
  • Fig. 5B is a graphical illustration like that of
  • Fig. 5C is a graphical illustration of the time reversal of the waveform of Fig. 5B.
  • Fig. 5D is a graphical illustration of the complex conjugate of the waveform of Fig. 5C.
  • Fig. 7 is a simplified representation of an image which would result from the transformation of the data set shown in Fig. 6.
  • Fig. 9 is a graph showing the improvement of the relative signal-to-noise ratio as the number of echoes sampled increases.
  • Fig. 10 is a sequence of graphs showing the phases of the components of a series of echoes.
  • rf coil 17 is tuned to the resonant frequency of the substance being imaged. Therefore, rf coil 17 may operate both to provide rf pulses and to detect rf signals from body 18. These signals may then be processed to provide information about body 18, such as an image.
  • Fig. 2A is a functional block diagram showing the basic operations of NMR Fourier imaging 20 according to the invention. Block 30 illustrates that the first operation is to apply the pulse and gradient switching sequence. The gradients are applied using gradient coils 14, 15, 16 as shown in Fig. 1, while the pulses are applied using rf coil 17. As shown in block 40, if the sequence of pulses and gradients is appropriate, rf coil 17 will receive echoes containing both phase and amplitude information.
  • the first of the echoes is produced by the first refocusing pulse, shown in block 34a in Fig. 2B.
  • the mechanism by which the refocusing pulses produce echoes will be discussed in greater detail below in relation to Figs. 10 and 11.
  • the excitation pulse is a 90° pulse along the positive Y axis
  • the first refocusing pulse of pulse and gradient switching sequence 30 will be a 180° rf pulse along the X axis. This pulse serves to reflect the spin magnetization vectors around the X axis, so that as the T 2 decay continues, the magnetization is refocused to form a spin echo decaying according to T 2 *.
  • the excitation pulse is a 90° pulse along the positive Y axis
  • the first refocusing pulse of pulse and gradient switching sequence 30 will be a 180° rf pulse along the X axis.
  • This pulse serves to reflect the spin magnetization vectors around the X axis, so that as the T 2 decay continues, the magnetization is refocused
  • Pulse and gradient switching sequence 30 continues with a third, refocusing pulse, as shown in block 34c in Fig. 2B.
  • This third refocusing pulse is a 180o rf pulse along the X axis.
  • a fourth refocusing pulse is applied as shown in block 34d in the form of a 180o rf pulse along the Y axis.
  • the pulse and gradient switching sequence 30 of the invention also includes a series of refocusing pulses along the negative X and Y axes, as shown in Fig. 2B.
  • the fifth refocusing pulse, shown in block 36a is a 180° rf pulse along the negative X axis
  • the sixth refocusing pulse, shown in block 36b is a 180o rf pulse along the Y axis.
  • the seventh and eighth refocusing pulses shown in blocks 36c and 36d are also 180o rf pulses along the negative X and Y axes. Therefore, the fifth through eighth refocusing pulses compensate pulse phase error in the same manner as the first through fourth refocusing pulses shown in blocks 34a-34d.
  • each refocusing pulse rather than being a full 180o pulse, is an effective pulse of 170o, and therefore introduces an error of 10o per pulse. This error may result from inaccurate setting or other causes.
  • both the in-phase and out-of-phase components of the echoes are fully compensated for the error after each group of eight spin echoes.
  • both components have values close to the full relative magnetization even for those echoes which are not fully compensated, with the maximum loss of signal less than 8%.
  • the in-phase components of the echoes resulting from the CP and CPMG pulse sequences are compensated, but the out-of- phase components deteriorate rapidly.
  • the in-phase CP and CPMG components are very well compensated, with a loss of 1.5% on alternate refocusing pulses.
  • the out-of-phase CP and CPMG components show no substantial compensation, so that the amplitude varies widely. After nine refocusing pulses, each in error by 10°, the amplitude goes to zero, indicating that the full 90° of error is uncompensated. As a result, the phase information is lost with the CP and CPMG pulse sequences.
  • the refocusing pulses of the invention are applied along X and Y axes in the rotating frame of reference.
  • gradient fields along the X, Y and Z axes in the laboratory frame of reference may also be applied for spatial encoding of the NMR data.
  • Figs. 3A and 3B show in greater detail the timing of the pulses discussed above.
  • gradients 82 and 84 are applied after the excitation pulse 72.
  • Gradient 82 may include both G x and G z components. These components of gradient 82 are applied only during the evolution period between initial excitation pulse 72 and the first refocusing pulse 76.
  • the spin axes of the atoms begin in a plane perpendicular to the Y axis and move toward the Y axes at a rate determined by gradient 84 which has a component G y .
  • Gradient 82 performs a phase encoding function, which encodes the phase of the magnetization in the X-Z plane according to the spatial location along the X and Z axes in the laboratory frame of reference. The amount of phase encoding depends on the strength of the components of the field gradient 82 and the location of the spins in the
  • the G x and G z components of gradient 82 may be applied successively.
  • the G y component of gradient 84 is necessary because refocusing gradients 86 and 86' also include G y components.
  • gradient 84 and 86 must satisfy the condition expressed in the following equation:
  • t represents the time from excitation pulse 72
  • t 1 represents the time from the peak of the first echo 78
  • t 1 represents the time from the peak of the first echo 78
  • t 2 represents the time between the excitation pulse 72 and the first refocusing pulse 76.
  • ⁇ is chosen to provide an optimal evolution, and the time between refocusing pulses will be 2 .
  • the first half of each refocusing gradient 86, 86' is a phase encoding gradient over time, while the second half serves to re-encode the phase of the magnetization along the Y axis. Because of this use of G y for two purposes, gradient 84 is necessary for encoding before the echoes begin.
  • a single excitation pulse 72 is followed by refocusing pulses, including first refocusing pulse 76 and second refocusing pulse 76'.
  • refocusing pulses including first refocusing pulse 76 and second refocusing pulse 76'.
  • a series of observations must be made, each beginning with an excitation pulse 72 followed by gradient 82 or gradient 88.
  • the phase encoding which may be measured in degrees, depends on the strength of these gradients. Therefore, a series of measurements may be taken using gradients of different strengths to obtain data sets with different amounts of phase encoding. These data sets may then be ordered according to phase encoding to provide a combined data set which is then transformed.
  • the pulse and gradient switching sequence 30 of the invention provides a series of echoes having specific characteristics.
  • the application of a series of pulse and gradient switching sequences 30 with different degrees of phase encoding will produce a set of echo trains, each of which may be combined to form one or more data sets according to the invention.
  • a series of pulse and gradient switching sequences 30 with different degrees of phase encoding will produce a set of echo trains, each of which may be combined to form one or more data sets according to the invention.
  • Equation (2) the signal is expressed as S, M 0 represents the equilibrium magnetization at the point x 0 , y 0 , z0 , and w is equal to the product of the gyromagnetic ratio ⁇ with the coordinate z 0 and the gradient G z .
  • ⁇ is the phase encoding term, given by:
  • each signal may be expressed as the sum of a real term and a complex term.
  • Each of these terms will be dependent on the phase encoding term ⁇ , resulting from the successive applications of G x and G z , as they affect the location x 0 ,
  • the resulting spin packet at the location x 0 , y 0 , z 0 will be the sum of all spin magnetization arising in a voxel at that position, and the phase of this spin packet may vary over a full 360° or 2 ⁇ radians.
  • Signals S A and S B are orthogonal, and together are equal to the entire signal S.
  • S A and S B may be obtained from a phase sensitive detector. Therefore, the signals obtained at the first echo for each of these components will be:
  • Fig. 10 shows schematically the magnetization resulting from some pulses of pulse and gradient switching sequence 30.
  • excitation pulse 72 a 90° rf pulse along the Y axis in the rotating frame of reference, causes the magnetization to rotate through 90° from the Z axis to the X axis, resulting in initial magnetization 110.
  • initial magnetization 110 evolves into a form having in-phase component A and out-of-phase component B.
  • first refocusing pulse 76a a 180° rf pulse along the X axis
  • first echo magnetizations 112 and 122 result along the X and Y axes, respectively.
  • the components A and B of the fourth echo will have the same phase relation to the components of the third echo as the components of the second echo have to the first echo.
  • the echo from each even-numbered refocusing pulse of pulse and gradient switching sequence 30 will have this same relation to the preceding odd-numbered echo, so that we may generally treat equations 4A and 4B as describing the odd echoes and equations 5A and 5B as describing the even echoes.
  • each echo received as a result of the pulse and gradient switching sequence 30 contains both amplitude and phase information which can be represented by the real and imaginary components of a complex function.
  • Figs. 2C and 2D are flow charts showing two ways in which these echoes may be coherently added according to the invention. In both cases, the coherent addition results in a data set in which both phase and amplitude information are present.
  • the process begins with block 42, in which the echoes are received and converted into stored digital signals. This may be done, for example, by a phase-sensitive detector followed by an analog-to-digital converter.
  • the result will preferably be a group of measurements for the in-phase and out-of-phase components similar to those in equations 4 and 5, above.
  • the out-of-phase components represented in equations 4B and 5B
  • this can be done by time reversing the alternate echoes, as shown in block 52. If we treat the second echo as the echo to be reversed in time, we may rewrite equations 5A and 5B by replacing t 2 with -t 2 as follows:
  • the next step, shown in block 54, is to take the complex conjugate of the alternate echoes, which yields:
  • equations 7A and 7B As can be seen by comparing equations 7A and 7B with equations 4A and 4B, the signs of the terms are now the same for the first echo and the second echo, and the two data sets can be added together, term by term. Furthermore, subsequent echoes can be treated in the same manner such that all odd numbered echoes are treated alike and all even numbered echoes are also treated alike. Therefore, the final step of this embodiment is to add the echoes as shown in blcck 56 of Fig. 2C. Since the time reversal of block 52 and the complex conjugate of block 54 do not affect the sign of the real components, the addition of the echoes in block 56 also results in a cancellation of any existing zero frequency or baseline components.
  • the real parts of the complex data set will be summed with the real parts of the complex data set will be summed with the signs (+, -, -, +) for each group of four echoes. Since the complex conjugate operation of block 54 changes the sign of imaginary zero frequency components, they will be summed with the signs (+, +, -, -), which will also result in cancellation of the baseline offsets. Therefore, the use of the pulse and gradient switching sequence 30 with coherent addition 50 results in cancellation of baseline offsets.
  • Figs. 5A-5D illustrate graphically the process of coherent addition according to Fig. 2C.
  • the in-phase signal S A is shown in solid line, while the out-of-phase signal S B is shown in dashed line.
  • Fig. 5A shows real components 42a and imaginary components 42b of an odd echo, each with an arbitrary illustrative waveshape.
  • Fig. 5B shows real components 42c and imaginary components 42d of an even echo, which will bear a characteristic relation to the components 42a and 42b of the odd echo.
  • Fig. 5C shows the even echo of Fig.
  • the second echo is also modified and becomes:
  • Equation 9A and 9B Comparing equations 9A and 9B with equations 4A and 4B, it is apparent that the second echo obtained from the negative application of the spatial encoding gradient is identical, except for relaxation effects, to the first echo from the positive application of the same spatial encoding gradient. Therefore, these two data sets may be added together without any additional data processing. Similarly, comparing equations 8A and 8B with equations 5A and 5B, it is apparent that the first echo from the negative spatial encoding gradient may be added to the second echo from the positive encoding gradient without additional data processing.
  • Fig. 2D shows an alternative method of coherently adding the echoes to form data sets which follows from the above discussion. In this method, the echoes are received and converted to stored digital signals in block 42, in the same manner as in Fig.
  • computer 90 is represented as including a central processing unit (CPU) 92, program memory 94, and a data memory 96.
  • the operation of the CPU 92 includes sending and receiving appropriate signals from the gradient coil control circuitry 102a, the rf pulse control circuitry 102b and the receiving circuitry 102c.
  • CPU 92 is connected to an appropriate operator interface 104, which may be a standard operator terminal.
  • CPU 92 is connected to a display device 106, which may be a CRT display or any other appropriate display having the necessary degree of resolution.
  • CPU 92 will be initiated by appropriate signals from the operator interface 104. In response to those signals, CPU 92 will initiate the pulse and gradient switching sequence 30 in response to the switching sequence program 94a stored in program memory 94. Execution of this program will cause CPU 92 to send signals to gradient coil control circuitry 102a and rf pulse control circuitry 102b to produce the pulse and gradient switching sequence 30. As a result, echoes will be received by receiving circuitry 102c and changed into digital form, as described above. As a part of the switching sequence program 94a, CPU 92 also operates to store the digitized echoes in the echo data memory 96a.
  • CPU 92 executes transform algorithm 94c in order to obtain the basic image data which is stored in display data memory 96c.
  • the display data may then be read out by CPU 92 to the display device 106 for display to the operator.
  • the switching sequence program 94a and coherent addition program 94b will perform the steps described above in order to preserve the phase information necessary to perform Fourier imaging.
  • these programs will preferably conform to the additional design constraints set forth below.
  • T 2 * used above to represent the apparent relaxation time, may be more narrowly defined as the relaxation time which takes into account the additional contribution from local field inhomogeneity over the resolution element or voxel in question.
  • T 2 ** can then be defined as the relaxation time which takes into account the field inhomogeneity over the whole object within the imaging region.
  • t s should be as short as possible, and should therefore have a value of twice T 2 **. This will result in an image which is relatively free of spatial distortions due to field inhomogeneity.
  • T 2 * places a lower limit on the strength of the field gradient necessary to produce an acceptable image.
  • the image will, however, be spatially distorted due to the effect of field inhomogeneity. This distortion may be acceptable for the particular imaging requirements or, alternatively, the distortion may be corrected by suitable software.
  • the S/N ratio may be improved by reducing the band width over which the NMR signal is received, because the noise power at each point is proportional to the square root of the frequency band at that point.
  • the frequency band is the reciprocal of the sampling time t s .
  • the signal S(M) which can be obtained by summing echoes, keeping all other parameters constant, is given by:
  • the ratio of T 2 /T 2 ** is much larger at higher fields than at lower fields, and the value of R(M) in equation 14 is increased.
  • T 2 /T 2 ** 3
  • the fourth echo has an intensity of 37% of the first echo and a sensitivity improvement of 30% is achieved by sampling four echoes.
  • T 2 /T 2 ** 30
  • a sensitivity improvement of 355% is achievable by sampling 38 echoes.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (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)

Abstract

Des images de Fourier basées sur des échos NMR (résonance magnétique nucléaire) multiples peuvent être obtenues en utilisant le procédé et l'appareil ci-décrit. Une séquence d'impulsions et de gradients sont appliqués à l'objet à prendre en image (30). Les impulsions commencent avec une impulsion d'excitation (72) qui est suivie d'une série d'impulsions de refocalisation (76, 76'). De manière à compenser l'erreur de phase dans les impulsions de refocalisation, celles-ci sont appliquées le long des axes alternatifs X et Y dans le cadre rotatif de référence. Les échos (78, 78') sont ensuite additionnés de manière cohérente pour conserver les informations de phase (56, 56'). Des échos alternés peuvent être inversés (52) dans le temps et leur conjugué complexe peut être pris (54) pour permettre une addition cohérente (56). Alternativement, les échos pairs (54') peuvent être ajoutés séparément par rapport aux échos impairs (56') pour produire des ensembles de données séparés. La séquence de commutation d'impulsions de gradients (30) fournit des inversions de signes de sorte que toute composante de ligne de base est annulée par l'addition cohérente. Avant ou après une addition cohérente, les échos sont transformés, et après formation d'un ensemble de données combinées à partir d'une série de séquences d'échos (50), l'ensemble de données combinées est de nouveau transformé pour produire une image de l'objet (60).
EP19850902856 1984-06-01 1985-05-21 Prise d'image de fourier nmr a partir d'echos multiples Withdrawn EP0182873A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US61628384A 1984-06-01 1984-06-01
US616283 1984-06-01

Publications (1)

Publication Number Publication Date
EP0182873A1 true EP0182873A1 (fr) 1986-06-04

Family

ID=24468776

Family Applications (1)

Application Number Title Priority Date Filing Date
EP19850902856 Withdrawn EP0182873A1 (fr) 1984-06-01 1985-05-21 Prise d'image de fourier nmr a partir d'echos multiples

Country Status (2)

Country Link
EP (1) EP0182873A1 (fr)
WO (1) WO1985005693A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115951281A (zh) * 2022-12-08 2023-04-11 无锡鸣石峻致医疗科技有限公司 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2574551B1 (fr) * 1984-12-12 1986-12-26 Commissariat Energie Atomique Procede de generation et de traitement de signaux pour l'obtention par resonance magnetique nucleaire d'une image exempte de distorsions a partir d'un champ de polarisation inhomogene
JPH0811112B2 (ja) * 1985-03-11 1996-02-07 株式会社日立製作所 核磁気共鳴を用いた検査装置
JPS62231647A (ja) * 1986-03-31 1987-10-12 株式会社東芝 磁気共鳴映像装置
IL82878A0 (en) * 1986-06-18 1987-12-20 Philips Nv Method of and device for reconstructing a nuclear magnetization distribution from a partial magnetic resonance measurement
DE3855288T2 (de) * 1987-06-24 1996-09-19 Picker Int Inc Verfahren und Apparat zur Bilderzeugung mittels magnetischer Resonanz
GB2217847B (en) * 1988-04-27 1993-02-03 Nat Res Dev Improvements in or relating to nmr imaging systems
JPH02215440A (ja) * 1989-02-17 1990-08-28 Toshiba Corp 磁気共鳴映像装置
GB9217718D0 (en) * 1992-08-20 1992-09-30 British Tech Group Method of and apparatus for nmr testing
US5304930A (en) * 1993-02-01 1994-04-19 Panacea Medical Laboratories Remotely positioned MRI system

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4297637A (en) * 1978-07-20 1981-10-27 The Regents Of The University Of California Method and apparatus for mapping lines of nuclear density within an object using nuclear magnetic resonance
DE2951537A1 (de) * 1979-12-20 1981-07-02 Siemens AG, 1000 Berlin und 8000 München Zeugmatografieverfahren
US4431968A (en) * 1982-04-05 1984-02-14 General Electric Company Method of three-dimensional NMR imaging using selective excitation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO8505693A1 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115951281A (zh) * 2022-12-08 2023-04-11 无锡鸣石峻致医疗科技有限公司 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质
CN115951281B (zh) * 2022-12-08 2023-10-20 无锡鸣石峻致医疗科技有限公司 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质

Also Published As

Publication number Publication date
WO1985005693A1 (fr) 1985-12-19

Similar Documents

Publication Publication Date Title
US4689567A (en) NMR Fourier imaging from multiple echoes
US5942897A (en) Magnetic resonance imaging apparatus
EP2233940B1 (fr) Excitation et acquisition simultanée en résonance magnétique
CN100370948C (zh) 磁共振成像方法和磁共振成像系统
US4993075A (en) Image reconstruction method in NMR imaging
US4748410A (en) Rapid NMR imaging system
JP3529446B2 (ja) Ｅｐｉ及びｇｒａｓｅ ｍｒｉにおける読み出し傾斜磁界極性の補正方法
JP3734086B2 (ja) 核磁気共鳴イメージング装置
US5652516A (en) Spectroscopic magnetic resonance imaging using spiral trajectories
US4628262A (en) Multiple echo chemical shift imaging
JPH0224135B2 (fr)
JPH0213571B2 (fr)
Talagala et al. Introduction to magnetic resonance imaging
EP0145276B1 (fr) Représentation du quotient complexe en résonance magnétique nucléaire
US5025216A (en) Magnetic resonance imaging of short T2 species
US4985677A (en) Magnetic resonance imaging and spectroscopy using an excitation pulse for multiple-dimensional selectivity
US4868501A (en) Method and means for magnetic resonance spin-echo imaging using an adiabatic three pi pulse sequence
US4689566A (en) NMR phase encoding using phase varying rf pulses
EP0182873A1 (fr) Prise d'image de fourier nmr a partir d'echos multiples
US4706027A (en) Method for correcting phase errors in magnetic resonance imaging data
US4697149A (en) NMR flow imaging using a composite excitation field and magnetic field gradient sequence
US4853635A (en) Method of reconstructing a nuclear magnetization distribution from a partial magnetic resonance measurement
CN110730912A (zh) 并行多切片mr成像
JP3173612B2 (ja) 磁気共鳴分光法および映像法における、またはそれに関連する改良
US20210356547A1 (en) Magnetic resonance imaging using motion-compensated image reconstruction

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): AT BE CH DE FR GB IT LI LU NL SE

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

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 19860506

RIN1 Information on inventor provided before grant (corrected)

Inventor name: MAUDSLEY, ANDREW, A.