Avoidance of Artifacts in Measurement Data Captured Using a Magnetic Resonance System
20230213603 · 2023-07-06
Assignee
Inventors
- Adam Kettinger (Erlangen, DE)
- Flavio Carinci (Würzburg, DE)
- Miriam Van De Stadt-Lagemaat (Erlangen, DE)
- Mario Zeller (Erlangen, DE)
- Dominik Paul (Bubenreuth, DE)
Cpc classification
G01R33/3852
PHYSICS
G01R33/56554
PHYSICS
G01R33/543
PHYSICS
G01R33/56518
PHYSICS
International classification
G01R33/561
PHYSICS
G01R33/54
PHYSICS
Abstract
A method for avoiding artifacts in measurement data captured using a magnetic resonance system which has a gradient unit. The method includes loading data which characterizes the gradient unit of the magnetic resonance system; loading a measurement protocol to be used for capturing the measurement data, wherein the measurement protocol includes gradients to be switched and RF excitation pulses and RF refocusing pulses to be irradiated, wherein, after irradiation of an RF excitation pulse, a train of at least two RF refocusing pulses is irradiated and measurement data is captured after each RF refocusing pulse; determining compensation gradients which, after the capture of the measurement data, are to be switched after a final RF refocusing pulse of the train of RF refocusing pulses associated with the RF excitation pulse and before a following RF excitation pulse as a function of the loaded measurement protocol and of the data which characterizes the gradient unit; and carrying out the measurement protocol using the determined compensation gradients.
Claims
1. A method for avoiding artifacts in measurement data captured using a magnetic resonance system which includes a gradient unit, the method comprising: loading data which characterizes the gradient unit of the magnetic resonance system; loading a measurement protocol to be used for capturing the measurement data, wherein the measurement protocol includes gradients to be switched and radio frequency (RF) excitation pulses and RF refocusing pulses to be irradiated, wherein, after irradiation of an RF excitation pulse, a train of at least two RF refocusing pulses is irradiated, and the measurement data is captured after each RF refocusing pulse; and determining compensation gradients which, after the capture of the measurement data, are to be switched after a final RF refocusing pulse of the train of RF refocusing pulses associated with the RF excitation pulse and before a following RF excitation pulse as a function of the measurement protocol and of the data which characterizes the gradient unit; and carrying out the measurement protocol using the compensation gradients.
2. The method as claimed in claim 1, wherein the gradient unit comprises at least three gradient coils, and the data which characterizes the gradient unit is loaded for each of the gradient coils included in the gradient unit.
3. The method as claimed in claim 1, wherein the determination of the compensation gradients comprises a transformation of the data which characterizes the gradient unit from a physical coordinate system of the gradient coils of the gradient unit into a logical coordinate system of the measurement protocol.
4. The method as claimed in claim 1, wherein the determination of the compensation gradients comprises convolution of gradients which are to be switched and are specified by the measurement protocol in the train of RF refocusing pulses preceding the compensation gradient with the data which characterizes the gradient unit.
5. The method as claimed in claim 1, wherein the data which characterizes the gradient unit characterizes an eddy current characteristic of the gradient unit.
6. The method as claimed in claim 5, wherein the determination of the compensation gradients takes into account at least one point in time at which contributions from eddy currents are to be avoided or at least reduced, wherein the at least one point in time is selected from a group consisting of: i) point in time of irradiation of the RF excitation pulse following the compensation gradients, ii) points in time of capture of central measurement data, and iii), if the measurement protocol comprises RF saturation pulses, point in time of irradiation of an RF saturation pulse before the RF excitation pulse following the compensation gradients.
7. The method as claimed in claim 6, wherein the determination of the compensation gradient comprises minimizing a measure which measures an influence of eddy current effects for the at least one point in time.
8. The method as claimed in claim 1, wherein the measurement protocol specifies a time frame for compensation gradients to be determined.
9. The method as claimed in claim 1, wherein the data which characterizes the gradient unit is measured once on the magnetic resonance system.
10. The method as claimed in claim 1, wherein the data which characterizes the gradient unit is based on measured gradient impulse response functions.
11. The method as claimed in claim 10, wherein the measured gradient impulse response functions have been measured on the magnetic resonance system for a series of low frequencies located within a half-value width of a spectral density-frequency diagram, for frequencies of less than 10 kHz.
12. The method as claimed in claim 1, wherein the data which characterizes the gradient unit is determined on the basis of a model which comprises measured decay profiles of the gradient unit.
13. A magnetic resonance system comprising a magnet unit, a gradient unit, a radio-frequency unit and a control device with a radio-frequency transmit/receive controller and having a compensation gradient determination unit, wherein the control device is configured to carry out the method as claimed in claim 1 on the magnetic resonance system.
14. A non-transitory electronically readable data storage medium having electronically readable control information stored thereon, which comprises at least one computer program and is designed such that, when the data storage medium is used in a control device of a magnetic resonance system, the information carries out the method as claimed in claim 1.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0032] Further advantages and details of the present disclosure are revealed by the exemplary aspects described below and with reference to the drawings. The stated examples do not limit the disclosure. In the figures:
[0033]
[0034]
[0035]
DETAILED DESCRIPTION
[0036]
[0037] Data GRD which characterizes the gradient unit of the magnetic resonance system is loaded (block 101).
[0038] If the gradient unit comprises at least two, in particular three, gradient coils, the data which describe the gradient unit may be loaded for each of the coils which the gradient unit comprises, wherein each gradient coil of the gradient unit acts in one axis of the physical coordinate system (X,Y,Z).
[0039] In particular, the data which characterizes the gradient unit characterizes an eddy current characteristic of the gradient unit.
[0040] The data which characterizes the gradient unit may, as has already been mentioned above, be based on measured gradient impulse response functions (GIRFs).
[0041] The gradient impulse response functions on which the data which characterizes the gradient unit is based may have been measured on the magnetic resonance system for a series of frequencies, in particular for low frequencies, for example frequencies located within the half-value width of the spectral density-frequency diagram, for example for frequencies of less than 10 kHz. Such a measurement of the GIRFs also for long time constants with a correspondingly low frequency means that long-term eddy currents can also still be effectively described or predicted by convolution according to (1).
[0042] The data which characterizes the gradient unit may, as has already been mentioned above, alternatively be determined on the basis of a model which comprises measured decay profiles of the gradient unit.
[0043] A measurement protocol MP to be used for capturing the measurement data MD, which protocol comprises gradients to be switched and RF excitation pulses and RF refocusing pulses to be irradiated, is loaded (block 103). Spatial encoding during capture of the measurement data is carried out in a logical coordinate system (x,y,z), wherein for example the z direction corresponds to the slice selection direction, the y direction to the phase-encoding direction and the x direction to the read-out direction. The logical coordinate system is here selected in accordance with the location and position of a test report of an object under examination from which measurement data is to be captured. The logical coordinate system (x,y,z) therefore generally does not correspond to the physical coordinate system (X,Y,Z) in which the magnetic resonance system is oriented. A selected logical coordinate system (x,y,z) also depends on the location and position of the object under examination in the magnetic resonance system. The logical coordinate system (x,y,z) may for example be selected by a user or also automatically as a function of the examination region to be examined and its position and location in the magnetic resonance system, wherein it is conceivable for the logical coordinate system to be adapted over the course of a measurement, for example after a movement of the object under examination.
[0044] The RF excitation pulses, RF refocusing pulses and gradients to be switched which the measurement protocol MP comprises provide that, after irradiation of an RF excitation pulse, a train of at least two RF refocusing pulses is irradiated and measurement data is captured after each RF refocusing pulse. The measurement protocol MP therefore specifies capture of measurement data MD by means of a TSE sequence, or at least of a TSE-like sequence.
[0045] One example of a possible sequence specified by the measurement protocol MP will be shown later with the assistance of part of a pulse sequence scheme shown diagrammatically in
[0046] As a function of the loaded measurement protocol MP and of the data GRD which characterizes the gradient unit, compensation gradients KG are determined which are to be switched after capture of measurement data MD after a final RF refocusing pulse of a train of RF refocusing pulses associated with an RF excitation pulse and before a following RF excitation pulse (block 105).
[0047] Determination of a compensation gradient may comprise a transformation of the data which characterizes the gradient unit from the physical coordinate system (X,Y,Z) of the gradient coils of the gradient unit into a logical coordinate system (x,y,z) of the measurement protocol.
[0048] Data which characterizes the gradient unit present in the physical coordinate system (X,Y,Z) may here be transformed into data present in logical coordinates (x,y,z) in accordance with the same rotation matrix R which has already been used for conversion between the coordinate systems (X,Y,Z) and (x,y,z) during the generation of all the gradients, for which the following for example applies:
[0049] Thus, the following is analogously obtained for example in the case of GIRFs as the data which characterizes the gradient unit and was measured for the physical coordinate system (X,Y,Z):
[0050] Once, through such a transformation of the coordinate system, the data which characterizes the gradient unit is available in the logical coordinate system corresponding to the image to be captured, all further steps may proceed in the logical coordinate system.
[0051] Such a transformation of the coordinate system may proceed with a constantly up-to-date rotation matrix R on each determination of compensation gradients during the measurement time which is required for capturing all the desired measurement data. In this way, dynamic adaptation of compensation gradients is achieved in the event of possible changes in the location of an object under examination or of the desired location of the region under examination, and thus a possible change in the logical coordinate system and the rotation matrix R.
[0052] Determination of a compensation gradient may, for example according to one of equations (1) or (3), comprise convolution of (nominal) gradients which are to be switched and are specified by the loaded measurement protocol in a repeat time TR, in which the compensation gradient to be determined is located, with the loaded data which characterizes the gradient unit. It is conceivable for (nominal) gradients which are to be switched also to be taken into account in determining the compensation gradients before the repeat time TR in which the compensation gradient to be determined is located.
[0053] The determination of compensation gradients may be limited to the slice selection direction (z direction) and the read-out direction (x direction) in the logical coordinate system if perturbations in the phase-encoding direction (y direction) are considered negligible.
[0054] The measurement protocol may for example specify a nominal starting compensation gradient, the time profile and in particular the amplitude of which may also already be selected, for example according to a known method as a function of the nominal gradients to be switched in the measurement protocol in (and optionally before) the repeat time TR in which the compensation gradient to be determined is located. Such starting compensation gradients (G.sup.nominal) may, as described above, be transformed into a currently valid logical coordinate system (x,y,z) and the actual gradients (G.sup.actual) thereof ascertained by convolution with the data which characterizes the gradient unit, e.g. GIRFs, and has likewise been transformed into the logical coordinate system (x,y,z). For such nominal starting compensation gradients in the x direction and in the z direction, this gives rise, for example for GIRFs, to:
[0055] It should once again be mentioned that, when determining the actual gradients (G.sup.actual), account should only be taken of nominal gradients which are to be switched before an RF excitation pulse, before which compensation gradients to be determined, although the compensation gradients only take effect timewise at the earliest from this subsequent RF excitation pulse.
[0056] Determination of a compensation gradient may take account of at least one desired point in time at which artifacts, in particular due to contributions from eddy currents should be avoided or at least reduced. Such desired points in time may be selected, for example depending on the nature of the measurement to be made, from the group of: [0057] i) point in time of irradiation of the RF excitation pulse following the compensation gradient, [0058] ii) points in time of capture of central measurement data, and [0059] iii) if the measurement protocol comprises RF preparation pulses, such as for example RF saturation pulses, in particular RF fat saturation pulses, point in time of irradiation of an RF saturation pulse before the RF excitation pulse following the compensation gradient.
[0060] By taking account of such desired points in time, compensation gradients can be determined in a targeted manner in such a way that compensation is optimum at the desired points in time, such that an influence of eddy current effects is minimized in particular at the desired points in time.
[0061] To this end, a suitable measure which measures the influence of eddy current effects may for example be minimized for the desired points in time. For example the 2 norm, also known as the standard norm, of the actual gradients can be used as such a measure, which can be minimized iteratively, and, for example, separately for all logical axes x, y, and z concerned.
[0062] The desired points in time can furthermore be provided with a weighting. Such weighting may be achieved by way of weighting factors assigned to the desired points in time, which factors are for example selected such that desired points in time which are classed as more important for image quality are more heavily weighted than desired points in time which are classed as less important. Such weighting can on the one hand accelerate minimization, wherein it is on the other hand ensured that perturbing artifacts are avoided as efficiently as possible.
[0063] For the example used above of compensation gradients in the x direction and in the z direction and in the notation used above, such an optimization may be written:
with t.sub.i as the desired points in time, and w.sub.i as their weighting factors.
[0064] Once the described (optionally weighted) 2 norms have been determined, the profile of the associated compensation gradients, in particular the duration, amplitude and ramps thereof, can be modified according to a known minimization method, for example according to a conjugate gradient (CG) method, and the 2 standards redetermined until a minimum is reached. Once the minimum of the measure which measures an influence of eddy current effects has been identified, the compensation gradients are determined.
[0065] In order to accelerate minimization for all the desired points in time, a compensation gradient identified as the minimum for a desired point in time can here be used as a starting compensation gradient for further desired points in time. In particular, since convolution is a linear operation and providing no other changes are made to the pulse sequence used, it is possible to adjust such a minimum compensation gradient identified for a desired point in time only by time shifts corresponding to the time distances of the further desired points in time from the desired point in time for which a compensation gradient has already been identified, in order to obtain the compensation gradients associated with the further desired points in time.
[0066] If parameters, in particular the gradients to be switched, of the measurement protocol are changed during ongoing measurement, for example due to a change in the logical coordinate system (x,y,z), the compensation gradients must be redetermined taking the change into account.
[0067] It should be noted that, for the stated data which characterizes the gradient unit, actual gradients can in each case be determined according to the equations (1) and (3) by way of convolution with the respective data which characterizes the gradient unit, wherein convolution is convolved once directly with the nominal gradients (GIRF), and once with the first time derivative of the nominal gradients (decay profiles).
[0068] While convolution of the nominal gradients with GIRFs directly leads (1) to actual gradients (which include eddy current effects), convolution of the first time derivative of the nominal gradients with the stated decay profiles yields (3) the pure perturbation produced by eddy current effects.
[0069] As stated, however, the compensation gradients should only take effect timewise at the earliest from an RF excitation pulse following the compensation gradient. At the point in time of irradiation of this subsequent RF excitation pulse, gradients to be switched in the repeat time TR before the RF excitation pulse are already switched off, whereby nominal gradients to be switched during or before a repeat time TR have the value zero from the point in time of a subsequent RF excitation pulse. In this way, by using GIRFs to determine the data which characterizes the gradient unit, equation (1) also yields the pure perturbation produced by eddy current effects for points in time when the compensation gradients should take effect, for example at desired points in time t.sub.i.
[0070] Therefore, when the stated decay profiles are used for determining the data which characterizes the gradient unit, it is possible to proceed analogously to the above GIRF-based data which characterizes the gradient unit.
[0071] In particular, decay profiles present in the physical coordinate system (X,Y,Z)
may be transformed into the logical coordinate system (x,y,z) analogously to equation (5) by way of the rotation matrix:
[0072] Furthermore (analogously to equations 6a and 6b in conjunction with equation 3), the perturbations in the axes of the logical coordinate system (x,y,z) caused by eddy current effects may be ascertained by convolving the first time derivative of the nominal gradients with the decay profiles present in the logical coordinate system:
[0073] By taking account of above-described desired points in time, compensation gradients can likewise be determined in a targeted manner in such a way that compensation is optimum at the desired points in time, such that an influence of eddy current effects is minimized in particular at the desired points in time.
[0074] To this end, a suitable measure which measures the influence of eddy current effects may be minimized for the desired points in time, wherein the desired points in time may in turn also be provided with a weighting (cf. equations 7a and 7b). The following is analogously obtained:
[0075] Determination of a compensation gradient may thus comprise minimizing a measure which measures the influence of eddy current effects for desired points in time. Determination of a compensation gradient may furthermore comprise weighting of the desired points in time.
[0076] The loaded measurement protocol MP is carried out using the specific compensation gradients KG (block 107), whereby measurement data MD is captured and wherein changes in the parameters, in particular of the gradients to be switched and/or of the logical coordinate system (x,y,z), of the measurement protocol may optionally be taken into account on determination of the compensation gradients.
[0077] On the basis of the captured measurement data MD, image data BD can be reconstructed (block 109), in which artifacts are reduced by using the compensation gradients according to the disclosure.
[0078] By determining compensation gradients according to the disclosure as a function of both a loaded measurement protocol MP to be used and of system-specific data which characterizes the gradient unit, the compensation gradients determined in this way are optimized not only for the measurement protocol but also for the magnetic resonance system with its gradient unit which has been used. In this way, troublesome artifacts, in particular perturbations produced by system-specific eddy current effects, can simultaneously be optimally avoided not only in a protocol-specific manner but also in a system-specific manner.
[0079]
[0080] The top line (RF) shows the time profile of RF-pulses to be irradiated, the second line (x) shows associated gradients to be switched in the read-out direction x, the third line (y) shows associated gradients to be switched in the phase-encoding direction y and the fourth line (z) shows associated gradients to be switched in the slice selection direction.
[0081] For the purposes of illustration and in order to focus on essentials, only four RF refocusing pulses RF2 are shown in a train of RF refocusing pulses RF2 after an RF excitation pulse RF1, and thus within a repeat time TR. In fact, distinctly more RF refocusing pulses RF2 may be used in a repeat time TR after an RF excitation pulse RF1. A repeat time TR is the time span between two successive RF excitation pulses RF1, in which a train of at least two, in this case four, RF refocusing pulses RF2 is irradiated, after which measurement data is in each case captured.
[0082] The illustrated pulse sequence scheme of the loaded measurement protocol MP, in each case after a final read-out process in which measurement data is captured, specifies, in a repeat time TR, a time frame KGT for compensation gradients KG to be determined, which may in principle be located on all the logical axes, thus all the directions x, y, and z. Compensation gradients KG according to the disclosure, as shown in dashed lines in the lefthand time frame KGT, were, for example, determined. In the next repeat time TR (which can be seen in its entirety in
[0083] Possible desired points in time at which artifacts, in particular due to contributions from eddy currents, are to be avoided or at least reduced, are indicated by the vertical lines t1, t2, t3, t4 and t5. In the example of
[0084]
[0085]
[0086] In order to examine an object under examination U, for example a patient or also a phantom, this can be introduced on a couch L into the measuring volume of the magnetic resonance system 1. The slice or slab S.sub.i is an exemplary target volume of the object under examination from which echo signals are to be captured and acquired as measurement data.
[0087] The control device 9 serves to control the magnetic resonance system 1 and can in particular control the gradient unit 5 by means of a gradient controller 5′ and the radio-frequency unit 7 by means of a radio-frequency transmit/receive controller 7′. The radio-frequency unit 7 may here comprise a plurality of channels on which signals can be transmitted or received.
[0088] The radio-frequency unit 7, together with its radio-frequency transmit/receive controller 7′, is responsible for generating and irradiating (transmitting) a radio-frequency alternating magnetic field for manipulating spins in a region to be manipulated (for example in slices S to be measured) of the object under examination U. The center frequency of the radio-frequency alternating magnetic field, which is also known as the B1 field, is usually set such that it is as close as possible to the resonant frequency of the spins to be manipulated. Deviations of the center frequency from the resonant frequency are known as off-resonance. The radio-frequency transmit/receive controller 7′ in the radio-frequency unit 7 applies controlled currents to the RF coils in order to generate the B1 field.
[0089] The control device 9 furthermore comprises a compensation gradient determination unit 15 with which compensation gradients according to the disclosure can be determined for avoiding artifacts, in particular for compensating eddy currents, which gradients can be implemented by the gradient controller 5′. The control device 9 is configured overall to carry out a method according to the disclosure.
[0090] A computing unit 13 included in the control device 9 is configured to carry out all the computing operations necessary for the required measurements and determinations. Intermediate and final results required or ascertained for this purpose can be stored in a memory unit S of the control device 9. The illustrated units should not necessarily be understood as physically separate units, but merely represent a subdivision into units of meaning which may, however, also be embodied in fewer or even just one single physical unit.
[0091] Control commands can be sent, for example by a user, to the magnetic resonance system and/or results from the control device 9, for example in the form of image data, can be displayed via an input/output device E/A of the magnetic resonance system 1.
[0092] A method described herein may also take the form of a computer program product which comprises a program and implements the described method on a control device 9 when it is executed on the control device 9. Likewise, an electronically readable data storage medium 26 with electronically readable control information stored thereon may be present, which control information comprises at least one such computer program product as has just been described and is designed to carry out the described method in a control device 9 of a magnetic resonance system 1 using the data storage medium 26.