Continual trajectory correction in magnetic resonance imaging
11585884 · 2023-02-21
Assignee
Inventors
- Manuel Stich (Parkstein, DE)
- Dominik Paul (Bubenreuth, DE)
- Mario Zeller (Erlangen, DE)
- Nadine Dispenza (Erlangen, DE)
Cpc classification
G01R33/5608
PHYSICS
International classification
G01R33/565
PHYSICS
Abstract
A correction method for reducing temperature-related deviations in a gradient response of an MR pulse sequence in MR imaging is provided. An MR pulse sequence that includes at least one nominal test gradient is run. A gradient response to the at least one nominal test gradient is repeatedly acquired by a magnetic field measurement in an examination region. A gradient system transfer function is determined based on the gradient response. A corrected MR pulse sequence is determined based on the gradient system transfer function and of the at least one nominal test gradient.
Claims
1. A correction method for reducing temperature-related deviations in a gradient response of a magnetic resonance (MR) pulse sequence in MR imaging, the correction method comprising: applying an MR pulse sequence that includes at least one nominal test gradient, wherein the MR pulse sequence includes the at least one nominal test gradient as an additional gradient pulse; repeatedly acquiring a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region during an MR imaging process; determining a gradient characterization function based on the gradient response; and determining a corrected MR pulse sequence based on the gradient characterization function and the at least one nominal test gradient.
2. The correction method of claim 1, wherein repeatedly acquiring the gradient response comprises acquiring the gradient response at least after every fifth repetition interval.
3. The correction method of claim 2, wherein the repeatedly acquiring the gradient response comprises acquiring the gradient response in each repetition interval.
4. The correction method of claim 1, wherein the at least one nominal test gradient has a shape of a rectangle or a trapezium.
5. The correction method of claim 4, wherein the MR pulse sequence includes the at least one nominal test gradient as an integral component.
6. The correction method of claim 1, wherein the MR pulse sequence does not include any rectangular- or trapezoid-shaped gradient pulses, such that the at least one nominal test gradient has a shape that differs from a rectangle or a trapezium.
7. The correction method of claim 6, wherein the at least one nominal test gradient comprises at least one chirp pulse, and wherein determining the gradient characterization function comprises acquiring a gradient response of the at least one chirp pulse.
8. The correction method of claim 1, wherein the MR pulse sequence includes at least one spiral-shaped gradient.
9. The correction method of claim 1, wherein based on the gradient characterization function and the at least one nominal test gradient in a repetition interval, a gradient in the MR pulse sequence that has already been applied is corrected.
10. The correction method of claim 1, wherein based on the gradient system transfer function and the at least one nominal test gradient of a repetition interval, a corrected gradient in the MR pulse sequence of a subsequent repetition interval is determined.
11. A method for recording magnetic resonance image data relating to an examination object, the method comprising: acquiring raw data that is generated by a gradient pulse sequence; determining a corrected gradient pulse sequence, the determining of the corrected gradient pulse sequence comprising: applying the gradient pulse sequence, which includes at least one nominal test gradient as an additional gradient pulse; repeatedly acquiring a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region during an MR imaging process; determining a gradient characterization function based on the gradient response; and determining the corrected gradient pulse sequence based on the gradient characterization function and the at least one nominal test gradient; and reconstructing magnetic resonance (MR) image data based on the acquired raw data and the corrected gradient pulse sequence.
12. A correction facility comprising: an output interface configured to run a magnetic resonance (MR) pulse sequence that includes at least one nominal test gradient, wherein the MR pulse sequence includes the at least one nominal test gradient as an additional gradient pulse; an input interface configured to repeatedly acquire a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region; a determination unit configured to determine a gradient characterization function based on the gradient response; and a correction unit configured to determine a corrected MR pulse sequence based on the gradient characterization function and the at least one nominal test gradient.
13. A magnetic resonance (MR) image reconstruction facility comprising: a correction facility comprising: an output interface configured to run a magnetic resonance (MR) pulse sequence that includes at least one nominal test gradient as an additional gradient pulse; a first input interface configured to repeatedly acquire a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region; a determination unit configured to determine a gradient characterization function based on the gradient response; and a correction unit configured to determine a corrected MR pulse sequence based on the gradient characterization function and the at least one nominal test gradient a second input interface configured to acquire MR raw data from the examination region; and an image reconstruction unit configured to reconstruct image data based on the acquired MR raw data and the corrected MR pulse sequence.
14. A magnetic resonance installation comprising: a main field magnet system; a radio frequency (RF) transmission antenna system; a gradient system; an RF receive antenna system comprising a control facility for activating the main field magnet system, the RF transmission antenna system, the gradient system, and the RF receive antenna system; and a correction facility comprising: an output interface configured to run a magnetic resonance (MR) pulse sequence that includes at least one nominal test gradient as an additional gradient pulse; an input interface configured to repeatedly acquire a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region; a determination unit configured to determine a gradient characterization function based on the gradient response; and a correction unit configured to determine a corrected MR pulse sequence based on the gradient characterization function and the at least one nominal test gradient.
15. In a non-transitory computer-readable storage medium that stores instructions executable by a computation unit to reduce temperature-related deviations in a gradient response of a magnetic resonance (MR) pulse sequence in MR imaging, the instructions comprising: applying an MR pulse sequence that includes at least one nominal test gradient, wherein the MR pulse sequence includes the at least one nominal test gradient as an additional gradient pulse; repeatedly acquiring a gradient response to the at least one nominal test gradient using a magnetic field measurement in an examination region during an MR imaging process; determining a gradient characterization function based on the gradient response; and determining a corrected MR pulse sequence based on the gradient characterization function and the at least one nominal test gradient.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
DETAILED DESCRIPTION
(14)
(15) In a lower partial diagram 12d, both gradients G.sub.nom, and G.sub.real are shown together. As shown by the partial diagram 12e, the real gradient G.sub.real is somewhat staggered over time with respect to the nominal gradient G.sub.nom.
(16) A section of the partial diagram 12d is shown enlarged in a partial diagram 12e. In the partial diagram 12e, the real gradient G.sub.real or the amplitude of the real gradient G.sub.real at t=122 s drops below the nominal gradient G.sub.nom and at around t=130 s assumes a higher value of magnitude than the nominal gradient G.sub.nom in order to coincide with the nominal gradient again on the zero line at values of t=160 s.
(17)
(18) As shown from diagram 20, peak temperatures of up to 60° C. are recorded at times, and the temperature range extends from around 20° C. to 60° C. These temperature differences lead to different response characteristics of the entire system to applied gradients. However, since the trajectory that the gradients describe in k-space should be known as precisely as possible for the reconstruction in order to obtain a good image quality, it is desirable to know the temperature-related deviations in the real gradient G.sub.real from the nominal gradient G.sub.nom or optionally even to compensate for these.
(19)
(20)
(21)
(22)
(23) The top partial diagram in
(24) The gradients G.sub.nom shown in
(25)
(26) A first partial diagram 70a shows a pulse sequence diagram illustrating a gradient sequence with a spiral-shaped trajectory for the gradients G.sub.x, G.sub.y in the x-direction and in the y-direction. The spiral-shaped gradients are not suitable for calculating a gradient system-transfer function since the spiral-shaped gradients do not have a broad frequency range. For this reason, in each repetition interval, an additional chirp gradient G.sub.CH is added to the respective spiral-shaped gradient. This chirp gradient G.sub.CH is shown in detail in a partial diagram 70b bottom left. The chirp gradient G.sub.CH includes a number of frequencies and is temporally short enough in order to fit into a repetition interval. The partial diagram 70b also shows a comparison of a nominal input gradient G.sub.CH_nom in chirp form with the actual gradient G.sub.CH_real that has been output and measured. For this purpose, a section of the chirp pulse illustrated is shown enlarged. The enlargement makes it clear that the amplitude of the nominal gradient G.sub.CH_nom is somewhat greater and also that a time lag exists between the nominal gradient G.sub.CH_nom and the real gradient G.sub.CH_real.
(27) In a third partial diagram 70c, the respective GSTF function or the magnitude thereof M is shown as a function of the frequency. The gradient system transfer function GSTF is derived by dividing the output gradient G.sub.CH_real by the input gradient G.sub.Ch_nom. The diagonal values of the gradient system transfer function GSTF or of the matrix thereof are shown.
(28) In order to obtain test gradients with a broad frequency spectrum despite this in the event of a non-rectangular- or non-trapezoid-shaped gradient shape, a chirp pulse is therefore added to the gradients G.sub.x, G.sub.y. A chirp pulse is a short frequency-sweep over a sufficient frequency range. This pulse has the advantage that the pulse is very short and therefore fits into a repetition time without any problems. A combination of triangular-shaped gradients is not suitable because the triangular-shaped would increase the measurement time considerably, which would lead to a current gradient system status not being able to be corrected soon enough.
(29) For both cases shown in
(30) In a first partial diagram 80a,
(31)
(32) In act 9.I, an MR pulse sequence is first applied (e.g., an EPI sequence), as shown in
(33)
(34)
(35)
(36)
(37) The magnetic resonance scanner 2 is equipped in the usual manner with a main field magnet system 4, a gradient system 6, and also an RF transmission antenna system 5 and an RF receive antenna system 7. In the exemplary embodiment shown, the RF transmission antenna system 5 is a whole body coil fixedly incorporated in the magnetic resonance scanner 2, whereas the RF receive antenna system 7 consists of local coils to be arranged on the patient or test subject (symbolized in
(38) The magnetic resonance tomography system shown in
(39) The magnetic resonance tomography system 1 further includes a central control facility 13 that is used to control the MR system 1. This central control facility 13 includes a sequence control unit 14. This controls the sequence of radiofrequency (RF) pulses and of gradient pulses as a function of a selected pulse sequence PS or of a succession of a plurality of pulse sequences for recording a plurality of slices in an interesting volume region of the examination object within a measurement session. Such a pulse sequence PS may be set and parameterized within a measurement or control protocol P, for example. Various control protocols P for different measurements or measurement sessions are usually stored in a memory 19 and may be selected by an operator (and optionally changed where necessary) and are then used to perform the measurement. In the present case, the control facility 13 contains pulse sequences for acquiring the raw data.
(40) To output the individual RF pulses in a pulse sequence PS, the central control facility 13 has a radiofrequency transmission facility 15 that generates the RF pulses, amplifies the RF pulses, and feeds the RF pulses via an appropriate interface (not shown in detail) into the RF transmission antenna system 5. To control the gradient coils in the gradient system 6 in order to switch the gradient pulses appropriately according to the set pulse sequence PS, the control facility 13 has a gradient system interface 16. The diffusion gradient pulses and also, for example, spoiler gradient pulses may be applied via this gradient system interface 16. The sequence control unit 14 communicates in an appropriate manner, for example, by transmitting sequence control data SD, with the radiofrequency transmission facility 15 and the gradient system interface 16 to carry out the pulse sequence PS.
(41) In addition, the control facility 13 has a radiofrequency receive facility 17 (likewise communicating in an appropriate manner with the sequence control unit 14) in order to receive magnetic resonance signals in a coordinated manner and therefore to acquire the raw data using the RF receive antenna system 7 within the readout windows ADC set in the pulse sequence PS.
(42) For example, a reconstruction unit 120a receives the acquired raw data and reconstructs magnetic resonance image data therefrom. This reconstruction may likewise ensue based on parameters that may be set in the respective measurement or control protocol P. This image data may then be stored, for example, in a memory 19.
(43) It is basically known in detail to a person skilled in the art how suitable raw data may be acquired, and MR images or parameter maps may be reconstructed therefrom by irradiating RF pulses and switching on gradient pulses; therefore, this is not explained in further detail here.
(44) The correction facility 110 is in data contact with the other units (e.g., with the gradient system interface 16 or the sequence control unit 14). Alternatively, the correction facility 110 may also be part of the sequence control unit 14. The correction facility 110 includes a plurality of units for determining or acquiring various parameters. First, appropriate test and reference pulse sequences g.sub.nom,l(t) are generated by the correction facility 110 and transmitted to the gradient system interface 16 and to the sequence control unit 14. Second, corrected gradients g.sub.pre,l(t) are generated by the correction facility 110 based on a determined gradient characterization function GSTF and are transmitted by the correction facility 110 to the reconstruction unit 120a, which carries out an image reconstruction on the basis thereof and on the basis of raw data.
(45) The central control facility 13 may be operated via a terminal 11 with an input unit 10 and a display unit 9, through which the entire magnetic resonance tomography system 1 may also be operated by an operator. The display unit 9 may also be used to display magnetic resonance tomography images, and using the input unit 10, optionally in combination with the display unit 9, measurements may be planned and started and, for example, control protocols P may be selected and optionally modified.
(46) The magnetic resonance tomography system 1 according to the present embodiments and, for example, the control facility 13 may in addition include a number of further components that are not shown here individually but are usually present in installations of this kind, such as, for example, a network interface to connect the entire system with a network and in order to be able to exchange raw data and/or image data or parameter maps and also further data, such as patient-relevant data, for example, or control protocols.
(47) It is basically known to a person skilled in the art how suitable raw data may be acquired and magnetic resonance tomography images may be reconstructed therefrom by irradiating RF pulses and generating gradient fields, and this is therefore not explained in further detail here. Likewise, a very wide range of measurement sequences, such as, for example, EPI measurement sequences or other measurement sequences for generating image data, and also diffusion-weighted image data are basically known to a person skilled in the art.
(48) The methods and structures described in detail in the aforementioned are exemplary embodiments, and the basic principle may also be varied in many ways by a person skilled in the art without departing from the scope of the invention insofar as set out in the claims. Further, the method described is likewise not restricted to medical applications. For the sake of completeness, the use of the indefinite article “a” or “an” does not preclude the relevant features from also being present in plurality. Likewise, the term “unit” does not preclude this from consisting of a plurality of components that may optionally also be spatially distributed.
(49) The elements and features recited in the appended claims may be combined in different ways to produce new claims that likewise fall within the scope of the present invention. Thus, whereas the dependent claims appended below depend from only a single independent or dependent claim, it is to be understood that these dependent claims may, alternatively, be made to depend in the alternative from any preceding or following claim, whether independent or dependent. Such new combinations are to be understood as forming a part of the present specification.
(50) While the present invention has been described above by reference to various embodiments, it should be understood that many changes and modifications can be made to the described embodiments. It is therefore intended that the foregoing description be regarded as illustrative rather than limiting, and that it be understood that all equivalents and/or combinations of embodiments are intended to be included in this description.