US12429541B2 - Measurement of a gradient field in an MRT system - Google Patents
Measurement of a gradient field in an MRT systemInfo
- Publication number
- US12429541B2 US12429541B2 US18/127,644 US202318127644A US12429541B2 US 12429541 B2 US12429541 B2 US 12429541B2 US 202318127644 A US202318127644 A US 202318127644A US 12429541 B2 US12429541 B2 US 12429541B2
- Authority
- US
- United States
- Prior art keywords
- slice
- gradient
- exciting
- disruption
- variable
- 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.)
- Active, expires
Links
Images
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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/243—Spatial mapping of the polarizing magnetic field
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
-
- 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/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
-
- 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/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
Definitions
- the present embodiments relate to methods for measurement of a gradient field in an imaging region of a magnetic resonance tomography (MRT) system, and to a corresponding MRT system.
- MRT magnetic resonance tomography
- Apparatuses for Magnetic Resonance Tomography are imaging apparatuses that use a strong external magnetic field in order to align the nuclear spins of an object to be examined and, by application of an RF excitation pulse, excite the nuclear spins to precess around the corresponding alignment.
- the precession or the transition of the spins from this excited state into a state with lower energy generates an electromagnetic alternating field as a response, which may be detected via receive antennas as an MR signal.
- a position encoding may be impressed onto the signals, which subsequently allows the signal received to be assigned to a volume element of the examination object.
- the received signal may then be evaluated in order, for example, to provide a graphical representation of the examination object.
- Deviations of the temporal course of the gradient fields may be caused by eddy currents, non-linearities, controller deviations of a gradient amplifier, or other transmission effects (e.g., filter characteristics of the hardware used).
- eddy currents which may lead to a location-dependent influencing of the gradient fields generated. Eddy currents with relatively short time constants, roughly in the region of a few ms and or less, are especially relevant.
- the spatial encoding during MR imaging may be influenced by the deviations or imperfections of the temporal course of the gradient fields, for example.
- the falsifications of the readout gradients may, for example, lead to artifacts, since the intended gradient trajectory is not realized precisely, and the recorded data is thereby sorted incorrectly into the k-space. Depending on the sequence technique used, various artifacts may arise thereby.
- Typical resulting artifacts are Nyquist ghosting or N/2 ghosting in EPI measurements. This is attributable to different polarities of the readout gradients being used in EPI for different k-space rows. The aforementioned deviations of the actual gradient shape from the intended shape lead to the rows in the k-space that are recorded with positive gradient polarity now experiencing a different disruption (e.g., by shifting of the echoes or accumulated phases) to the rows measured with negative polarity in the k-space experience.
- the problem is especially relevant for measurements with data acquisition during ramp sampling.
- the data that is recorded during the ramp sampling and thus corresponds to a non-equidistant k-space sampling may be initially converted into equidistant points in the k-space.
- This process is also referred to as gridding or re-gridding.
- the knowledge of the temporal course of the output readout gradients is important.
- the deviations mentioned may lead to an incorrect gridding, which may cause artifacts, such as “ringing” in the readout direction and “ghosting” during EPI measurement with ramp sampling.
- the latter relates above all to a type of ghosting, in which edges of the object appear as an especially strong ghost.
- Eddy currents may bring about different spatial distributions of field disruptions.
- the eddy currents may, for example, be divided into B0-like eddy currents that create a spatially constant field disruption, linear or gradient-like eddy currents that create a gradient-like field disruption, which essentially has the same spatial course as the gradient that created it, cross terms, and eddy currents of higher spatial order.
- Cross terms that, however, have a course that runs orthogonally to the gradient that creates the cross terms (e.g., when a gradient in the X direction creates an eddy current that leads to a field that runs linearly in the Y direction) create a gradient-like field disruption.
- Eddy currents of higher spatial order create field disruptions of which the course depends on more than one spatial coordinate or on a power of a coordinate.
- Document US 2016/0091582 A1 describes a method in which the temporal fluctuations of the magnetic field in an MRT volume are determined based on the rise in a phase difference between the spin responses in a number of slices. These are stored and used for more exact re-gridding of the acquired k-space data before the image reconstruction.
- a disadvantage in this case is that, through inherent averagings, field disruptions of a higher spatial order are not mapped or are not mapped completely.
- the present embodiments may obviate one or more of the drawbacks or limitations in the related art.
- space disruption influences such as disruption fields that influence a temporal course of gradient fields, are determined more precisely.
- disruptions of a higher spatial order may be determined more precisely.
- the present embodiments are based on the idea of exciting two slices different from one another and intersecting each other in a test object, switching a readout gradient thereafter, and acquiring a magnetic resonance (MR) signal emitted from a corresponding intersection region of the slices.
- MR magnetic resonance
- a disruption variable is computed. The disruption variable determines a deviation of a temporal course of an amplitude of the readout gradient from a predetermined required course.
- a method for measuring a gradient field in an imaging region of a magnetic resonance tomography (MRT) system is specified.
- a first slice in a test object which is located in the imaging region, is excited (e.g., by a radio frequency (RF) transmit coil arrangement and a gradient coil arrangement of the MRT system, which are controlled by at least one control unit of the MRT system in order to excite the first slice).
- RF radio frequency
- a second slice in the test object is excited (e.g., by the gradient coil arrangement and the RF transmit coil arrangement activated by the at least one control unit).
- the second slice intersects with the first slice within the test object in an intersection region.
- a readout gradient that may also be referred to as a test gradient is switched (e.g., by the gradient coil arrangement activated by the at least one control unit), and an MR signal emitted from the intersection region is acquired (e.g., by a signal acquisition apparatus of the MRT system).
- a disruption variable is computed (e.g., by a computing unit of the MRT system). The disruption variable determines a deviation of a temporal course of an amplitude of the readout gradient from a predetermined required course for the readout gradient (e.g., for the amplitude of the readout gradient).
- a gradient e.g., of a readout gradient
- a gradient field e.g., of the readout gradient
- a gradient field e.g., of the readout gradient
- the change in the amplitude of the magnetic field as a function of the spatial coordinates then corresponds to a gradient or gradient field in the meaning used.
- a temporal course of the amplitude of the gradient or of a number of amplitudes of gradients in different directions may also be referred to as the gradient trajectory.
- the test object may involve an object to be imaged (e.g., a patient) who is to be imaged during MRT imaging, which, for example, may follow the method for measuring of the gradient field of the present embodiments.
- the test object may, however, also involve a phantom object that is employed, for example, during calibration or during tune-up of an MRT system.
- the switching of a gradient may be time-dependent activation and deactivation or control of the gradient coil arrangement in order to create a corresponding temporal course of the gradient field.
- a gradient in this case may also be a gradient pulse with one or more parts, which are also referred to as lobes.
- the first slice and the second slice are different from one another and not parallel to one another, so that the first slice and the second slice intersect with one another within the test object.
- Each slice involves an approximately three-dimensional spatial region with a constant thickness, which is also referred to as the slice thickness. Accordingly, a volume region in which the first slice and the second slice intersect, apart from deviations through the shape of the test object, is approximately given by a three-dimensional volume in the form of a column. If the first slice and the second slice are at right angles to one another, then the column has an approximately rectangular or square cross section.
- the intersection region may correspond to the entire column or to a part thereof.
- the intersection region corresponds to the entire column (e.g., the entire volume region in which the first slice and the second slice intersect).
- the intersection region may involve a part of the column.
- intersection region thus always involves a volume region that is contained in all excited slices involved. Accordingly, the MR signal is also emitted from the intersection region.
- what is achieved by the excitation of the first slice and of the second slice and also optionally the third slice is that, on switching of the readout gradient and/or after the switching of the readout gradient, the MR signal that is emitted from the intersection region is acquired.
- the MR signal involves a signal as a result of a resonant spin response of the excited nuclear spins in the intersection region. Depending on the form of embodiment of the method, this may involve a spin echo or a so-called stimulated echo.
- the MR signal involves a spin echo when just the first slice and the second slice are excited, but not the third slice, whereas the MR signal involves a stimulated echo when the first slice, the second slice, and the third slice are excited.
- the acquisition of the MR signal includes the creation of corresponding MR data or MR test data dependent on the MR signal.
- the MR signal may be detected, for example, by one or more antennas or receive coils of the signal acquisition apparatus and be converted by receive electronics of the signal acquisition apparatus into the MR test data.
- the receive electronics for example, include an Analog-to-Digital Converter (ADC).
- ADC Analog-to-Digital Converter
- the acquisition electronics e.g., the ADC
- the acquisition electronics may be activated, for example, while the readout gradient is being switched and/or after the readout gradient has been switched.
- the required course of the amplitude of the readout gradient corresponds to the temporal course of the amplitude of the readout gradient theoretically to be expected under ideal conditions (e.g., when there is no disruption), such as eddy currents or the like, present.
- the required course may, for example, be given by a temporal course of the current and/or a control signal for activation of the gradient coil apparatus for switching of the readout gradient or vice versa.
- the readout gradient is switched in accordance with the predetermined required course or depending on the predetermined required course.
- the disruption variable is, for example, a function of the time and may also be referred to as the disruption function.
- the disruption variable may directly correspond to the deviation of the temporal course of the amplitude from the predetermined required course.
- the deviation of the temporal course of the amplitude from the required course may be computed explicitly from the disruption variable and where necessary from further known information.
- the disruption variable is, for example, independent of the readout gradient and the required course. This has the advantage that the disruption variable, once determined or computed, may also be used for computing a deviation of a further readout gradient, slice selection gradient, or further gradient, from a corresponding further required course.
- the determination of the disruption variable may thus be understood as a kind of disruption calibration.
- the disruption variable in this case quantifies the influence of disruption effects, such as eddy currents, on the creation of the readout gradient and its actual gradient trajectory.
- the disruption variable however, also quantifies the influence of the disruption effects on other gradient trajectories.
- the disruption variable may therefore be determined once and then be stored for further use.
- both the first slice and also the second slice are excited before the readout gradient is switched and the MR signal is accordingly emitted not from the entire first slice and not from the entire second slice, but just from the intersection region, provides that the disruption variable represents the disruption influences in the intersection region and not the disruption influences averaged over the first or second slice.
- This thus allows a spatially-resolved measurement of the gradient field and thus of the disruption variable to be achieved, which is not just restricted to a resolution in one slice selection direction, but also in a further direction.
- the disruption variable is thus able to be quantified with an increased spatial resolution and accordingly more accurately.
- the disruption variable is computed as a variable independent of the required course.
- the disruption variable is thus, for example, computed as a time-dependent function e(t), where e(t) does not depend either on the temporal course of the readout gradients or on the required course, but merely on the disruption influences themselves. This may be achieved, for example, by a corresponding model description being chosen for the disruption influences and being reconciled with the measurements (e.g., the MR test data) that are computed based on the acquired MR signal. Since the physical interrelationships such as an eddy current or another disruption that may influence the readout gradients are known per se, with the knowledge of the required course from the measured MR test data, the independent disruption variable e(t) may thus be computed.
- the disruption variable may then be stored (e.g., on a memory medium of the MRT system) and be used for further uses (e.g., in a method for magnetic resonance imaging following on from the method for measuring the gradient field).
- the disruption variable may in this case, for example, be taken into account within the framework of the eddy current compensation mentioned at the outset and/or during gridding in order to avoid artifacts.
- the first slice is excited by a first RF pulse being emitted into the imaging region and by a first slice selection gradient being switched at least partly at the same time as the first RF pulse.
- RF in this case stands for radio frequency and is used here and below with the meaning usually used in the context of MR imaging.
- the frequency of RF pulses lies in the order of magnitude of the Larmor frequency of the nuclear spins to be excited (e.g., depending on the size of the basic magnetic field, in the region of 10 MHz through to several 100 MHz).
- the second slice is excited by a second RF pulse being emitted in the imaging region and a second slice selection gradient being switched at least partly at the same time as the second RF pulse.
- the first slice and the second slice selection gradient point in this case in different directions, so that the first slice and the second slice correspondingly intersect in the test object.
- a first slice in the test object is excited (e.g., by a further first RF pulse being emitted in the imaging region and a further first slice selection gradient being switched at least partly at the same time as the further first RF pulse).
- a further second slice is excited in the test object (e.g., by a further second RF pulse being emitted in the imaging region and a further second slice selection gradient being switched at least partly at the same time as the further second RF pulse).
- the further second slice intersects with the further first slice in this case within the test object in a further intersection region.
- a further readout gradient is switched, and a further MR signal emitted from the further intersection region is acquired.
- the disruption variable is computed depending on the MR signal and the further MR signal (e.g., by the computing unit).
- the predetermined required course is both a required course for the readout gradients and also for the further readout gradients.
- the shape of the readout gradients is the same as the shape of the further readout gradients.
- the further readout gradient is switched in accordance with predetermined required course or depending on the predetermined required course.
- the first slice is the same as the first further slice, and the second slice is parallel to the second further slice. In one embodiment, the first slice is at right angles to the second slice in this case.
- the disruption variable is computed depending on a phase difference between the MR signal and the further MR signal.
- d corresponds to the spatial distance between the second slice and the second further slice.
- d corresponds to the spatial distance between the first slice and the first further slice.
- an approach is selected for the disruption variable in accordance with which the disruption variable is defined by a plurality of parameters independent of one another, and respective values are computed for the plurality of independent parameters by a fit method.
- a i refers to the amplitudes
- ⁇ i refers to the (e.g., positive) time constants of individual eddy current amounts to the disruption variable e(t).
- N is further a predetermined natural number (e.g., greater than or equal to 1).
- the plurality of the parameters independent of one another then corresponds to A i and ⁇ i .
- optimal values for A i and ⁇ i are determined (e.g., by the measured phase difference, the predetermined required course, and the approach e(t) being set in the above relationship of the equation (1)).
- Such forms of embodiment are, for example, advantageous when the third slice is not excited and, for example, the further third slice is not excited; the MR signal or the further MR signal is thus attributable to a spin echo.
- a further required course for the gradient pulse may be corrected depending on the disruption variable, and the at least one gradient pulse may be switched as a function of the corrected required course.
- an MRT system is specified.
- the MRT system has an imaging region, an RF transmit coil arrangement, and a gradient coil arrangement.
- the imaging region is surrounded by the RF transmit coil arrangement and the gradient coil arrangement.
- the MRT system has a signal acquisition apparatus and at least one control unit that is configured to activate the RF transmit coil arrangement and the gradient coil arrangement in order to excite a first slice in a test object that is located in the imaging region.
- the at least one control unit is, for example, configured to activate the RF transmit coil arrangement to emit a first RF pulse into the imaging region, and to activate the gradient coil arrangement to switch, at least partly at the same time as the first RF pulse, first slice selection gradients in order to excite the first slice.
- the at least one control unit is configured to activate the RF transmit coil arrangement and the gradient coil arrangement in order, after the excitation of the first slice, to excite a second slice in the test object.
- the second slice intersects with the first slice (e.g., within the test object) in an intersection region.
- the at least one control unit is, for example, configured to activate the RF transmit coil arrangement to emit a second RF pulse into the imaging region and to activate the gradient coil arrangement to switch a second slice selection gradient, at least partly at the same time as the second RF pulse, in order to excite the second slice.
- the at least one control unit is configured to activate the gradient coil arrangement, after the excitation of the second slice, to switch a readout gradient and to activate the signal acquisition apparatus to acquire an MR signal emitted from the intersection region.
- the MRT system has a computing unit that is configured, depending on the MR signal, to compute a disruption variable that determines a deviation of a temporal course of an amplitude of the readout gradient from a predetermined required course.
- the MRT system directly follow the various embodiments of the method and vice versa.
- individual features and corresponding explanations with regard to the different forms of embodiment for the method may be transferred by analogy to the corresponding forms of embodiment of the MRT system.
- the MRT system is configured or programmed to carry out a method of the present embodiments.
- the MRT system of the present embodiments carries out the method of the present embodiments.
- a computer program with commands is specified.
- the commands When the commands are executed by an MRT system of the present embodiments, the commands cause the MRT system to carry out a method of the present embodiments for measuring a gradient field and/or a method for MRT of the present embodiments.
- a computer-readable memory medium e.g., a non-transitory computer-readable storage medium
- the non-transitory computer-readable storage medium stores a computer program of the present embodiments.
- a computing unit may, for example, be understood as a data processing device that contains a processing circuit.
- the computing unit may thus, for example, process data for carrying out computing operations. Where necessary, this also includes operations for carrying out indexed accesses to a data structure (e.g., a Look-Up Table (LUT)).
- a control unit may therefore represent or contain a computing unit or be part of a computing unit.
- the computing unit may, for example, contain one or more computers, one or more microcontrollers, and/or one or more integrated circuits (e.g., one or more Application-Specific Integrated Circuits (ASIC), one or more Field Programmable Gate Arrays (FPGA), or one or more systems on a ship).
- the computing unit may also contain one or more processors (e.g., one or more microprocessors, one or more Central Processing Units (CPU), one or more Graphics Processing Units (GPU), and/or one or more signal processors, such as one or more digital signal processors (DSPs)).
- the computing unit may also contain a physical or a virtual network of computers or other of the units.
- the computing unit includes one or more hardware and/or software interfaces and/or one or more memory units.
- a memory unit may be embodied as a volatile data memory (e.g., as a Dynamic Random Access Memory (DRAM) or Static Random Access Memory (SRAM) or as a non-volatile data memory, such as a Read-Only Memory (ROM), as a Programmable Read-Only Memory (PROM), as an Erasable Programmable Read-Only Memory (EPROM), as an Electrically Erasable Programmable Read-Only Memory (EEPROM), as a flash memory or flash EEPROM, as Ferroelectric Random Access Memory (FRAM), as a Magnetoresistive Random Access Memory (MRAM), or as Phase-Change Random Access Memory (PCRAM)).
- DRAM Dynamic Random Access Memory
- SRAM Static Random Access Memory
- ROM Read-Only Memory
- PROM Programmable Read-Only Memory
- EPROM Erasable Programmable Read-Only Memory
- EEPROM Electrically Erasable Programmable Read-Only Memory
- FRAM Ferroelectric Random Access Memory
- MRAM Magnetoresist
- a component of the MRT system especially the computing unit, the at least one control unit or the signal acquisition apparatus of the MRT system, being configured, embodied, designed, or the like to carry out or realize a specific function, achieve a particular effect, or serve a particular purpose
- this may be understood such that the component, above and beyond the principle or theoretical usability or suitability of the component for this function, effect, or this purpose, through a corresponding adaptation, programming, physical embodiment, and so forth, is in a position in concrete terms and actually to carry out or to realize the function, to achieve the effect, or to serve the purpose.
- FIG. 1 shows a schematic diagram of an embodiment of a magnetic resonance tomography (MRT) system
- FIG. 4 shows schematics of excited slices in a test object in accordance with an embodiment of a method for measuring a gradient field
- FIG. 5 shows schematics of excited slices in a test object in accordance with a an embodiment of a method for measuring a gradient field
- FIG. 6 shows schematics of excited slices in a test object in accordance with an embodiment of a method for measuring a gradient field
- FIG. 7 shows schematically an MR pulse diagram in accordance with a further embodiment of a method for measuring a gradient field.
- the MRT system 1 includes a magnet unit with a field magnet 3 that creates a static magnetic field for alignment of nuclear spins of an object 8 (e.g., of a patient) in an imaging region.
- the imaging region is characterized by an extremely homogeneous static magnetic field, where the homogeneity relates, for example, to the magnetic field strength or to its amplitude.
- the imaging region is located in a patient tunnel 2 that extends in a longitudinal direction Z through the magnet unit.
- the field magnet 3 may, for example, be a superconducting magnet that may create magnetic fields with a magnetic flux density of up to 3 T or more. However, for lower field strengths, permanent magnets or electromagnets with normally conducting coils may also be used.
- a patient table 7 may be movable within the patient tunnel 2 .
- the magnet unit further includes a gradient coil arrangement 5 with a number of gradient coils that serve to superimpose gradient fields (e.g., location-dependent magnetic fields) on the static magnetic field in the three spatial directions for spatial differentiation of the sampled image areas in the imaging region.
- the gradient coils of the gradient coil arrangement 5 may be configured, for example, as coils of normally conducting wires that, for example, may create fields or field gradients orthogonal to one another in the imaging region.
- the MRT system 1 also includes a control and computing system.
- the control and computing system may include a receive control unit 10 that may be seen as part of the signal acquisition apparatus and is linked to the body coil 4 and/or the local coil 6 .
- the receive control unit 10 which may include an Analog-to-Digital Converter (ADC), may create corresponding MR data (e.g., in the k-space).
- a computing unit 9 of the control and computing system may evaluate the MR data and carry out an image reconstruction, for example.
- the control and computing system includes a transmit control unit 11 that is linked to the body coil 4 and controls the coil for creation of RF pulses, such as excitation pulses and/or refocusing pulses.
- the control and computing system includes a gradient control unit 12 that is linked to the gradient coil arrangement 5 and controls the unit in order to switch slice selection gradients, gradients for frequency, and/or phase encoding and/or readout gradients.
- the MRT system 1 is, for example, capable of carrying out a method for MRT in accordance with the present embodiments, which, for example, includes a method of the present embodiments for measuring a gradient field in the imaging region.
- the gradient control unit 12 and the transmit control unit 11 activate the gradient coil arrangement 5 or the body coil 4 to first excite a first slice 17 a and then a second slice 18 a in the object, where the second slice 18 a intersects with the first slice 17 a in an intersection region 19 a , as shown in FIG. 4 .
- the first slice is at right angles to the X direction and thereby parallel to the Y-Z plane; the second slice 18 a is at right angles to the Y direction and thereby parallel to the X-Z plane.
- Other orientations are possible, however.
- embodiments in which the first slice 17 a and the second slice 18 a are not at right angles to one another may be provided.
- first slice 17 a and the second slice 18 a are, however, at right angles to one another, then the two slices 17 a , 18 a intersect with one another in a volume region with the shape of a square or a column with a rectangular or quadratic cross-section, which extends in parallel to the Z direction.
- the intersection region 19 a may correspond to this column or, in forms of embodiment in which a further third slice is excited, which, for example, is at right angles to the first slice 17 a and to the second slice 18 a , a part of the column.
- the transmit control unit 11 may, for example, activate the body coil 4 to emit a first RF pulse 13 a (e.g., an excitation pulse with a flip angle of 90°), and the gradient control unit 12 may activate the gradient coil arrangement 5 , essentially at the same time, to switch a first slice selection gradient 14 a , as is shown schematically in the pulse diagram of FIG. 3 .
- the first slice selection gradient 14 a is switched in the X direction.
- the gradient control unit 12 activates the gradient coil arrangement 5 to switch a readout gradient 15 (e.g., in the X direction) in order to detect an MR signal of a spin echo from the first intersection region 19 a .
- the ADC of the receive control unit 10 is activated in a readout window 16 .
- the readout gradient 15 may, for example, lie within the readout window 16 or come before the window.
- the readout gradient 15 shown may also be only partly output during the readout window ADC or lie before this, in order to investigate the disruption effects existing after a gradient pulse.
- curve G Shown schematically in FIG. 2 as curve G is the temporal course of the readout gradient 15 and the corresponding required course as curve G s .
- G and G s may be expressed as follows
- G ⁇ ( t ) G s ( t ) - dG s ( t ) dt * e ⁇ ( t ) , ( 3 )
- e(t) refers to the disruption variable
- the operator * represents a convolution.
- a i refers to the amplitudes
- ⁇ i refers to the time constant of the individual disruption amounts (e.g., eddy current amounts). Different methods may be used for determination of the amplitudes and time constants.
- a measurement may be carried out at at least one further position, as is shown in FIG. 5 by way of example.
- a further first slice 17 b is excited, which, for example, is parallel to the first slice 17 a .
- the first slice 17 a may be centered at an X position X1
- the further first slice 17 b may be centered at an X position X2.
- the first slice 17 a is excited once again, and thereafter, a further second slice 18 b , which, for example, is parallel to the second slice 18 a , is excited.
- the second slice 18 a may be centered at a Y position Y1
- the further second slice 18 b may be centered at a Y position Y2.
- the first slice 17 a and the further second slice 18 b may then intersect in the second intersection region 19 b .
- the readout gradient 15 may be switched once again in order to measure a further MR signal from the second intersection region 19 b.
- a phase difference 4 ( t ) may be determined, and from a first temporal derivation of the MR signal and the further MR signal, the gradient trajectory G(t) may be determined
- G ⁇ ( t ) 1 2 ⁇ ⁇ ⁇ d ⁇ ⁇ ⁇ d ⁇ ⁇ dt , ( 5 )
- ⁇ is the gyromagnetic ratio
- d in the example of FIG. 5 corresponds to the distance X2-X1 of the first slice 17 a from the further first slice 17 b
- Y2-Y1 of the second slice 18 a from the further second slice 18 b corresponds to the distance Y2-Y1 of the second slice 18 a from the further second slice 18 b.
- the amplitudes and time constants of the disruption variable e(t) may be determined.
- the temporal course of gradient fields and field disruptions of higher spatial order may be measured, which, for example, are caused by eddy currents.
- too short-term effects e.g., eddy currents with very short time constants
- a particular application purpose is the correction of gradient trajectories during gridding.
- disruption of higher spatial order may be approximated locally by gradient-like disruption.
- additional reference measurements may be carried out, in which measurements are made at the respective slice positions without output of the readout gradients 15 .
- This data from the reference measurements may, for example, be subtracted from the data of the measurement with the readout gradients 15 (e.g., sample by sample) in order not to suppress effects in focus such as drift effects of the B0-field or influences of concomitant fields.
- a specific number of repetitions of the otherwise same measurements may be output. After this, the corresponding data is then averaged.
- FIG. 4 , FIG. 5 , and FIG. 6 are examples; other possible combinations are likewise possible. It is likewise possible to use yet further gradient pulses (e.g., for spoiling of unwanted signals, such as for suppression of signals resulting from Free Induction Decay (FID)).
- FID Free Induction Decay
- FIG. 7 For a development of the method, a pulse diagram is shown in FIG. 7 .
- at least one third RF pulse 13 c with a simultaneous third slice selection gradient 14 c is used, where, for example, all slice selection gradients 14 a , 14 b , 14 c used point in different directions, (e.g., orthogonal to one another).
- the method with two RF pulses 13 a , 13 b utilizes the spin echo arising; the method with three RF pulses 13 a , 13 b , 13 c utilizes the stimulated echo arising.
- the flip angle of the second RF pulse 13 b in the first case is equal to 180° and those of the first RF pulses 13 a are equal to 90°, in order to create a maximum spin echo
- the flip angle of all three RF pulses 13 a , 13 b , 13 c is, for example, equal to 90°, in order to create a maximum stimulated echo.
- the created MR signal originates from their intersection region 19 a , 19 b .
- the signal originates from a column that is not restricted in the third coordinate.
- the edge lengths of the restricted sides of the column are given by the slice thicknesses that result from the two RF pulses 13 a , 13 b .
- the signal originates from a voxel, of which the edge lengths are given by the thicknesses of the slices excited by the three RF pulses 13 a , 13 b , 13 c.
- the spin echo signal may, if necessary, also be suppressed by spoiling.
- the choice of slice position is determined, for example, by the frequency of the corresponding RF pulses.
- measurements are carried out at two positions.
- the intersection regions 19 a , 19 b are offset in the respective direction (e.g., in the X direction), as shown in FIG. 5 .
- the present embodiments thus offer the option of measuring the temporal course of a gradient in one direction with a fixed other coordinate (e.g., in the case of FIG. 5 with Y0). This is, for example, advantageous when field disruptions of higher order that may be approximated locally by a gradient-like disruption are present.
- forms of embodiment offer the option of determining the dynamics of the created fields in a direction other than that of the gradients creating the fields, as shown in FIG. 7 .
- the cross term of the X-gradient may be determined in the Y direction, for example (e.g., the dynamic field course in the Y direction, which is created by the X gradient).
- the previously described spatial limitation for examination of a gradient may be limited in a further coordinate.
- the method with three RF pulses 13 a , 13 b , 13 c is also suitable for determining dynamic field distributions of a higher spatial order.
- measurements may be carried out at different voxels, and from the measured phase courses, the spatial distribution of the field effects of a gradient output on an axis may be measured. Thereby, three-dimensional field distributions, and thus also eddy currents of a higher spatial order may be determined.
- An analysis of the field distributions of higher spatial orders may be undertaken, for example, based on the development of the measured fields in spherical harmonics, which represent a set of orthonormal basic functions. It is thus possible, for example, from measurements at nine different positions, to determine field distributions of a higher spatial order up to the second order. With measurements at 16 positions, it is possible to determine distributions up to the third order and so forth.
- the method of the present embodiments may be carried out in different forms of embodiment in the calibration of the MRT system 1 (e.g., as part of a tune-up process).
- the information and data determined e.g., the disruption variable e(t)
- further parameters are also known, such as the positioning of the slices to be measured. Therefore, the calibration measurement may then be aligned so that, for example, the gradient course of a readout gradient may be measured directly at the position of a slice to be measured.
- the measured gradient course may be used to correct the measured data in the k-space.
- the gradient course may be measured at positions that are restricted in one or two coordinates. This enables the gridding correction to be carried out depending on the location (e.g., aligned to the measured slice in each case).
- transversal slices e.g., with slice selection in the Z direction
- readout gradients in the X direction the actual X gradient may be measured as function of the Z position. This then enables the gridding correction to be carried out as a function of the slice.
- a gridding correction may also be carried out in directions other than the readout direction (e.g., in the phase encoding or slice selection direction).
- forms of embodiment may be utilized for determination of eddy currents of a higher order and of cross terms, where shortest time constants may also be measured.
- the cross terms determined may then also be applied by pre-emphasis, for example.
- the data may also be used for further image corrections (e.g., for phase corrections in EPI measurements).
- Different forms of embodiment enable measurement time to be saved where necessary, since no complete segmentation is needed. For example, only measurements at the relevant positions may be carried out (e.g., at 16 voxels for a three-dimensional determination of the field courses in third spatial order).
- the measurements are also able to be integrated in different forms of embodiment as navigator or pre-scan measurements into imaging sequences. Further, forms of embodiment of do not require any additional sensor systems, such as field probes or the like.
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General 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)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
{dot over (ϕ)}(t)/(2πdγ)=G s(t)−∫0 t e(τ)Ġ s(τ)dτ, (1)
where t refers to the time, Gs refers to the required course, Ġs refers to the first temporal derivation of the required course, {dot over (ϕ)} refers to the first temporal derivation of the phase difference, γ refers to the gyromagnetic ratio, and d refers to a spatial distance between the intersection region and the further intersection region.
e(t)=Σi=1 N A i e −t/τ
where e(t) refers to the disruption variable, and the operator * represents a convolution. The disruption variable itself may be expressed as follows
e(t)=Σi=1 N A i e −t/τ
where γ is the gyromagnetic ratio, and d in the example of
Claims (17)
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP22164891 | 2022-03-29 | ||
| EP22164891.8 | 2022-03-29 | ||
| EP22164891.8A EP4253976A1 (en) | 2022-03-29 | 2022-03-29 | Measurement of a gradient field in an mri system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| US20230314539A1 US20230314539A1 (en) | 2023-10-05 |
| US12429541B2 true US12429541B2 (en) | 2025-09-30 |
Family
ID=80979000
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US18/127,644 Active 2043-12-13 US12429541B2 (en) | 2022-03-29 | 2023-03-28 | Measurement of a gradient field in an MRT system |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US12429541B2 (en) |
| EP (1) | EP4253976A1 (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0184840A1 (en) | 1984-12-14 | 1986-06-18 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Arrangement for the spatially resolved inspection of a sample by means of magnetic resonance of spin moments |
| EP0386822A2 (en) | 1989-03-04 | 1990-09-12 | Philips Patentverwaltung GmbH | Nuclear magnetic resonance investigation process and apparatus for utilizing this process |
| US5451877A (en) | 1993-04-23 | 1995-09-19 | Siemens Aktiengesellschaft | Method for the compensation of eddy currents caused by gradients in a nuclear magnetic resonance apparatus |
| US6335620B1 (en) * | 1998-12-22 | 2002-01-01 | Siemens Aktiengesellschaft | Method of acquiring eddy currents that are caused by switched magnetic field gradients of a magnetic resonance apparatus and that contain cross-terms |
| US7002343B2 (en) | 2003-02-13 | 2006-02-21 | Siemens Aktiengesellschaft | Method for determining a compensation setting for an eddy current field |
| WO2013002232A1 (en) | 2011-06-30 | 2013-01-03 | 株式会社 日立メディコ | Magnetic resonance imaging device and method for measuring gradient magnetic field output waveform thereof |
| US9046590B2 (en) * | 2011-03-31 | 2015-06-02 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients |
| US20160091582A1 (en) | 2014-09-30 | 2016-03-31 | Kabushiki Kaisha Toshiba | Mri gradient trajectory mapping |
| US10042019B2 (en) * | 2014-06-24 | 2018-08-07 | Siemens Aktiengesellschaft | Frequency monitoring of gradient pulses in magnetic resonance imaging |
-
2022
- 2022-03-29 EP EP22164891.8A patent/EP4253976A1/en not_active Withdrawn
-
2023
- 2023-03-28 US US18/127,644 patent/US12429541B2/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0184840A1 (en) | 1984-12-14 | 1986-06-18 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Arrangement for the spatially resolved inspection of a sample by means of magnetic resonance of spin moments |
| US4748409A (en) | 1984-12-14 | 1988-05-31 | Max-Planck-Gesellschaft Zur Foerderung Der Wissenschaften E.V. | Imaging magnetic spin resonance method and apparatus |
| EP0386822A2 (en) | 1989-03-04 | 1990-09-12 | Philips Patentverwaltung GmbH | Nuclear magnetic resonance investigation process and apparatus for utilizing this process |
| US5451877A (en) | 1993-04-23 | 1995-09-19 | Siemens Aktiengesellschaft | Method for the compensation of eddy currents caused by gradients in a nuclear magnetic resonance apparatus |
| US6335620B1 (en) * | 1998-12-22 | 2002-01-01 | Siemens Aktiengesellschaft | Method of acquiring eddy currents that are caused by switched magnetic field gradients of a magnetic resonance apparatus and that contain cross-terms |
| US7002343B2 (en) | 2003-02-13 | 2006-02-21 | Siemens Aktiengesellschaft | Method for determining a compensation setting for an eddy current field |
| US9046590B2 (en) * | 2011-03-31 | 2015-06-02 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients |
| WO2013002232A1 (en) | 2011-06-30 | 2013-01-03 | 株式会社 日立メディコ | Magnetic resonance imaging device and method for measuring gradient magnetic field output waveform thereof |
| US10042019B2 (en) * | 2014-06-24 | 2018-08-07 | Siemens Aktiengesellschaft | Frequency monitoring of gradient pulses in magnetic resonance imaging |
| US20160091582A1 (en) | 2014-09-30 | 2016-03-31 | Kabushiki Kaisha Toshiba | Mri gradient trajectory mapping |
Non-Patent Citations (7)
| Title |
|---|
| Atkinson, I. C. et al; "Characterization and Correction of System Delays and Eddy Currents for MRI Imaging with Ultrashort Echo-Time and Time-Varying Gradients"; Magnetic Resonance in Medicine; vol. 62; pp. 532-537; 2009. |
| Barmet, C. et al.; "A third-order field camera with microsecond resolution for MR system diagnostics", 2009, Proc. Intl. Soc. Mag. Reson. Med., vol. 17, pp. 781. |
| Barmet, Christoph et al.; "Spatiotemporal Magnetic Field Monitoring for MR"; Magnetic Resonance in Medicine, vol. 60, pp. 187-197, 2008. |
| Bernstein, Matt A. et al. "Handbook of MRI Pulse Sequences", Elsevier, 2004; ISBN-13:978-0-12-092861-3, Chapter 13.2.; 16 pp. |
| Ma, C. et al.; "A new eddy-current compensation method in MRI."; PIERS online vol. 3, No. 6; 2007, pp. 874-878. |
| Machine Translation of WO2013-002232A1 (Year: 2013). * |
| Zhao X. et al.; "Impact of gradient imperfections on bone water quantification with UTE MRI." Magnetic resonance in medicine vol. 84(4); 2020; pp. 2034-2047. |
Also Published As
| Publication number | Publication date |
|---|---|
| US20230314539A1 (en) | 2023-10-05 |
| EP4253976A1 (en) | 2023-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10012713B2 (en) | Method and device for determination of a magnetic resonance control sequence | |
| JP6422700B2 (en) | Method for obtaining T1 time of water and T1 time of fat and magnetic resonance apparatus | |
| US8698495B2 (en) | Flip angle imaging with improved B1 mapping for multi-RF transmit systems | |
| US9097779B2 (en) | Magnetic field insensitive CEST imaging | |
| CN103323802B (en) | Magnetic resonance system controls the determination of sequence | |
| US8558547B2 (en) | System and method for magnetic resonance radio-frequency field mapping | |
| US7999545B2 (en) | Method for obtaining amplitude and phase profiles of RF pulses for spatially selective excitation | |
| CN104101853B (en) | Method and magnetic resonance equipment for determining resonant frequency deviation | |
| JP6258050B2 (en) | MR data acquisition method, B1 magnetic field determination method and correspondingly configured magnetic resonance apparatus | |
| US9279872B2 (en) | Method for determining a set of B1 field maps | |
| US8890527B1 (en) | Methods of radio frequency magnetic field mapping | |
| US10551466B2 (en) | Correction of a magnetic resonance transmission signal | |
| US20040164737A1 (en) | Method for determining the B1 field strength in MR measurements | |
| US10416258B2 (en) | Controlling magnetic resonance systems | |
| CN103969608B (en) | Acquisition of magnetic resonance data and the method and corresponding magnetic resonance equipment for determining B1 magnetic fields | |
| CN103454606B (en) | Method to determine subject-specific b1 distribution in magnetic resonance apparatus | |
| US10126383B2 (en) | RF pulse alignment | |
| US12429541B2 (en) | Measurement of a gradient field in an MRT system | |
| US9933497B2 (en) | Method and magnetic resonance apparatus for determination of excitation profiles from excitation pulses | |
| US9146292B2 (en) | Method and device for determining a magnetic resonance system control sequence | |
| US7279894B2 (en) | Method and magnetic resonance tomography apparatus for location coding with a non-linear gradient field | |
| JP2016131847A (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
| AS | Assignment |
Owner name: SIEMENS HEALTHINEERS AG, GERMANY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SIEMENS HEALTHCARE GMBH;REEL/FRAME:066267/0346 Effective date: 20231219 Owner name: SIEMENS HEALTHINEERS AG, GERMANY Free format text: ASSIGNMENT OF ASSIGNOR'S INTEREST;ASSIGNOR:SIEMENS HEALTHCARE GMBH;REEL/FRAME:066267/0346 Effective date: 20231219 |
|
| AS | Assignment |
Owner name: SIEMENS HEALTHCARE GMBH, GERMANY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:KOEHLER, MICHAEL;REEL/FRAME:069377/0605 Effective date: 20231026 |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
| STCF | Information on status: patent grant |
Free format text: PATENTED CASE |