Classifications

• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
• • 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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/443—Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
• • 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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
  • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
  • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
  • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
• • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
  • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
  • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
  • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0

Definitions

• the invention relates to the field of magnetic resonance (MR) imaging. It concerns a method of MR imaging of an object placed in a main magnetic field within an examination volume of an MR system.
• the invention also relates to an MR system and to a computer program to be run on an MR system.
• Image-forming MR methods which utilize the interaction between magnetic fields and nuclear spins in order to form two-dimensional or three-dimensional images are widely used nowadays, notably in the field of medical diagnostics, because for the imaging of soft tissue they are superior to other imaging methods in many respects, do not require ionizing radiation and are usually not invasive.
• Geometric distortion of clinical MR images needs to be assessed and corrected if necessary.
• System-level distortions are sufficiently corrected in state-of-the-art MR systems, but the subject-induced distortions due to susceptibility and chemical shift effects still need to be accounted for and corrected.
• a common approach for quantifying and correcting subject-induced geometric distortions is to acquire a map of the main magnetic field (B o map) in a separate pre-scan, i.e. prior to the acquisition of the actual MR image.
• the B o map is used to assess the inhomogeneity of the main magnetic field from the phase difference of two gradient echo images acquired at different echo times.
• the Bo map can be used to correct for main field inhomogeneity during processing of the acquired MR signals so as to reduce or eliminate geometrical distortions in the reconstructed MR images. This applies in particular to imaging techniques such as EPI (echo planar imaging) and to techniques that use non-Cartesian k-space trajectories, such as spiral k-space scans.
• imaging techniques such as EPI (echo planar imaging) and to techniques that use non-Cartesian k-space trajectories, such as spiral k-space scans.
• Bo shimming uses electric currents applied to dedicated coils that generate corrections to the main magnetic field in the examination volume of the MR system, making the field inside that volume as uniform as possible. Those corrections are generated in the form of field gradients of different geometries.
• the process of adjusting the currents passing through the shim coils is referred to as shimming and can be performed any time.
• the determined B o map can be used to derive appropriate shimming currents to compensate for any detected main field inhomogeneity. J.P. Stockmann et al.
• B o map estimates of the B o map faces several challenges, some of which are common to all body sites, such as poor signal-to-noise ratio and phase-wrapping at locations with large susceptibility effects, and others that are more relevant for specific body portions, such as chemical shift. Due to signal contributions from fat and water protons, the B o map can be corrupted by chemical shift artefacts, particularly for body sites with substantial fat content.
• More advanced methods for B o determination estimate the B o map as well as fat and water contributions, thereby reducing or essentially eliminating chemical shift artefacts.
• the relative contributions from fat and water to the MR signal can be calculated if information from two or more corresponding echoes, acquired at different echo times, is combined. This can be considered as chemical shift encoding, in which an additional dimension, the chemical shift dimension, is defined and encoded by acquiring two or more MR images at slightly different echo times.
• these types of experiments are often referred to as Dixon-type of measurements.
• Dixon MR imaging or Dixon water/fat MR imaging a water/fat separation is achieved by calculating contributions of water and fat from two or more corresponding echoes, acquired at different echo times. In general, this separation is possible because there is a known precessional frequency difference of hydrogen in fat and water.
• a method of MR imaging of an object placed in a main magnetic field within an examination volume of an MR system comprises the steps of: subjecting the object to a first imaging sequence to generate a first set of MR signals at a first echo time, acquiring the first set of MR signals from the object, subjecting the object to a second imaging sequence to generate a second set of MR signals at a second echo time which is different from the first echo time, wherein the first and second imaging sequences further differ with respect to at least one of: field of view, image resolution, repetition time, k-space sampling pattern, relaxation time weighting, acquiring the second set of MR signals from the object, and deriving a water map, a fat map, and a B o map from the first and second sets of MR signals.
• the known Dixon approach is based on acquiring MR signal data at least at two different echo times and using computation to separate water and fat signal contributions and determine the B o map.
• the invention generalizes this concept with no additional acquisition step, namely by utilizing MR signal data acquired at two or more different echo times, but also using two (or more) different (first and second) imaging sequences.
• the first and second imaging sequences correspond to different (gradient echo) acquisitions. They may differ with respect to the field of view, the image resolution, the repetition time, the k-space sampling pattern, the relaxation time weighting etc.
• the invention thus enables using distinct gradient echo scans for water/fat separation and Bo mapping.
• the method of the invention saves scan time by removing the need for an explicit B o acquisition in a separate pre-scan.
• this concept is implemented in that the method of MR-imaging involves a preparatory phase and a diagnostic phase.
• the preparatory phase includes the spatially low-resolution survey scan than may be employed to generate a geometry planning and the spatially low-resolution calibration scan to determine the coil sensitivity profiles of the RF receiver array’s multiple coil elements.
• the diagnostic phase includes the diagnostic imaging scan in which the k- space data are acquired at spatially high-resolution encoding with the use of the geometry planning.
• the spatially low-resolution survey scan and the spatially low-resolution calibration scan form the first and second imaging scans at respective different echo times from which using the Dixon method, the water map, fat map and the Bo-map can be derived.
• the B 0 -map may be employed for accurate Bo- shimming and compensation for Bo-inhomogeneities during the diagnostic phase.
• the water and/or fat images may be employed for water-fat separation or fat suppression during the diagnostic phase.
• the first and second imaging sequences are associated with the survey sequence and a calibration sequence, respectively, wherein the survey sequence is used to provide an overview image of the object, and wherein the calibration sequence is used to determine spatial sensitivity profdes of an RF receiver array of multiple receiving antenna (coil) elements used in acquiring the MR signals (e.g. for parallel imaging).
• Such survey and calibration scans are typically executed anyway at the beginning of each imaging session, so they do not involve any additional scan time. They use different imaging sequences (in terms of field of view and image resolution).
• the invention proposes to use the MR signal data acquired in these pre-scans for deriving the B o map (and a water map and a fat map so as to avoid chemical shift artefacts). Since the Bo map smoothly varies in space, the coarse resolution of the survey and calibration scans is no obstacle.
• the B o mapping may simplify to a mere phase subtraction technique. While this is conventionally done using a dedicated acquisition (Bo pre-scan), the invention eliminates that need, and instead uses already acquired data (from the survey and calibration scans) to compute the B o map. Further, since the imaging sequences used for the survey and calibration scans are typically factory-set protocols, they can be tailored to improve the B o mapping.
• the computation involved for the Bo mapping corresponds essentially to a conventional (two-point) Dixon method (see Eggers et al., J. Magn. Reson. Imaging 2014; 40:251-268).
• a signal model is used which includes the signal contributions from water and fat at each voxel position, the phase offset between water and fat induced by the difference in chemical shift (resulting from the chosen spectral model of fat and the two different echo time values of the first and second imaging sequences) and a (main field inhomogeneity-induced) phase error.
• two-point Dixon methods typically perform the separation in three steps: First, potential solutions for the phase error are calculated. This is independently done for each voxel. Then, one solution is selected for each voxel by maximizing the spatial smoothness of the phase error and, finally, the water and fat signal contributions are estimated given the phase error.
• the water map, the fat map, and the B o map are derived for an overlap area of the fields of view of the first and second imaging sequences.
• the object is subjected to a third imaging sequence to generate a third set of MR signals, wherein the image resolution of the third imaging sequence is higher than the image resolution of the first and/or second imaging sequences, wherein the third set of MR signals is acquired and reconstructed into an MR image, the reconstruction involving compensation for B o inhomogeneity based on the B o map derived from the first and second sets of MR signals.
• the first and second imaging sequences may represent pre-scans, such as the mentioned survey and calibration scans executed at a lower image resolution, wherein the third imaging sequence represents the main scan executed for the actual image acquisition at a higher image resolution.
• the B o map derived from the MR signal data acquired using the first and second imaging sequences can be used for correcting or suppressing artefacts related to main field inhomogeneity and chemical shift in the MR image reconstructed from the MR signal data acquired using the third imaging sequence.
• B o inhomogeneity can be compensated for during the application of the third imaging sequence by shimming of the main magnetic field based on the B o map derived according to the invention.
• the field of view of the third imaging sequence (or at least that part of the field of view in which a Bo inhomogeneity compensation is desired) should be contained in the overlap area of the fields of view of the first and second imaging sequences in order to have reliable B o data available throughout the relevant field of view of the third imaging sequence.
• the method of the invention described thus far can be carried out by means of an MR system including at least one main magnet coil for generating an essentially uniform, static main magnetic field B o within an examination volume, a number of gradient coils for generating switched magnetic field gradients in different spatial directions within the examination volume, at least one body RF coil for generating RF pulses within the examination volume and/or for receiving MR signals from a body of a patient positioned in the examination volume, a control unit for controlling the temporal succession of RF pulses and switched magnetic field gradients, and a reconstruction unit for reconstructing MR images from the received MR signals.
• the method of the invention can be implemented by a corresponding programming of the reconstruction unit and/or the control unit of the MR system.
• the method of the invention can be advantageously carried out on most MR systems in clinical use at present. To this end it is merely necessary to utilize a computer program by which the MR system is controlled such that it performs the above-explained method steps of the invention.
• the computer program may be present either on a data carrier or be present in a data network so as to be downloaded for installation in the control unit of the MR system.
• Fig. 1 shows an MR system for carrying out the method of the invention
• Fig. 2 shows an embodiment of the method of the invention as a flow chart
• Fig. 3 shows a Bo map derived from survey and calibration scan data according to the invention.
• an MR system 1 is shown as a block diagram.
• the device comprises superconducting or resistive main magnet coils 2 such that a substantially uniform, temporally constant main magnetic field B o is created along a z-axis through an examination volume.
• the device further comprises a set of (1 st , 2 nd , and - where applicable - 3 rd order) shimming coils 2’, wherein the current flow through the individual shimming coils of the set 2’ is controllable for the purpose of minimizing Bo deviations within the examination volume.
• a magnetic resonance generation and manipulation system applies a series of RF pulses and switched magnetic field gradients to invert or excite nuclear magnetic spins, induce magnetic resonance, refocus magnetic resonance, manipulate magnetic resonance, spatially and otherwise encode the magnetic resonance, saturate spins, and the like to perform MR imaging.
• a gradient pulse amplifier 3 applies current pulses to selected ones of whole-body gradient coils 4, 5 and 6 along x, y and z-axes of the examination volume.
• a digital RF frequency transmitter 7 transmits RF pulses or pulse packets, via a send/receive switch 8 and a body RF coil 9, into the examination volume.
• a typical MR imaging sequence is composed of a packet of RF pulse segments of short duration which, together with any applied magnetic field gradients, achieve a selected manipulation of nuclear magnetic resonance.
• the RF pulses are used to saturate resonance, excite resonance, invert magnetization, refocus resonance, or manipulate resonance and select a portion of a body 10 positioned in the examination volume.
• the MR signals are also picked up by the body RF coil 9.
• a set of local array RF coils 11, 12, 13 are placed contiguous to the region selected for imaging.
• the array coils 11, 12, 13 can be used to receive MR signals induced by RF transmissions via the body RF coil.
• the resultant MR signals are picked up by the body RF coil 9 and/or by the array RF coils 11, 12, 13 and demodulated by a receiver 14 preferably including a preamplifier (not shown).
• the receiver 14 is connected to the RF coils 9, 11, 12 and 13 via the send/receive switch 8.
• a host computer 15 controls the shimming coils 2’ as well as the gradient pulse amplifier 3 and the transmitter 7 to generate the imaging sequences of the invention.
• the receiver 14 receives a single or a plurality of MR data lines in rapid succession following each RF excitation pulse.
• a data acquisition system 16 performs analog-to-digital conversion of the received signals and converts each MR data line to a digital format suitable for further processing. In modem MR systems the data acquisition system 16 is a separate computer which is specialized in acquisition of raw image data.
• the digital raw image data are reconstructed into an image representation by a reconstruction processor 17 which applies a Fourier transform or other appropriate reconstruction algorithms, such as SENSE.
• the MR image may represent a planar slice through the patient, an array of parallel planar slices, a three-dimensional volume, or the like.
• the image is then stored in an image memory where it may be accessed for converting slices, projections, or other portions of the image representation into appropriate format for visualization, for example via a video monitor 18 which provides a man-readable display of the resultant MR image.
• the host computer 15 and the reconstruction processor 17 are programmed to execute the method of the invention as described above and in the following.
• a conventional survey scan is performed using a first imaging sequence.
• the survey scan provides a first set of MR signals at a first echo time, wherein an overview image of the body 10 is reconstructed from the first set of MR signals.
• a conventional calibration scan is performed to determine the spatial sensitivity profiles of the RF coils 11, 12, 13 used for acquiring the MR signals.
• the second set of MR signals is acquired at a second echo time which is different from the first echo time, wherein the first and second imaging sequences of the survey scan and the calibration scan, respectively, further differ with respect to the image resolution and the field of view, they apply different k-space sampling patterns, they are different, independent gradient echo sequences.
• the (first and second sets of) MR signal data acquired in these two pre-scans are used in step 23 for deriving a Bo map, a water map and a fat map using a conventional (generalized) Dixon method.
• the deriving of the B o map involves re-sampling of the complex MR signal data of at least one of the two (first and second) acquisitions to an identical, common k-space grid.
• a magnitude scaling factor is applied to at least one of the (re-sampled) first and second MR signal sets to take different relaxation time weightings or other amplitude effects associated with the two different imaging sequences into account.
• Such scaling factors can be estimated by simulations or by calibration.
• the Bo map, the water map and the fat map are computed by a conventional (two-point) Dixon method (see above).
• a third imaging sequence is applied to generate a third set of MR signals, wherein the image resolution of the third imaging sequence is higher than the image resolution of the first and second imaging sequences in steps 21, 22.
• the third imaging sequence represents the main scan executed for the actual image acquisition at a higher image resolution.
• the Bo map derived in step 23 is used for correcting or suppressing artefacts related to main field inhomogeneity and chemical shift in the MR image reconstructed from the MR signal data acquired in the main scan in step 24.
• the B o map can also be used for B o shimming during acquisition in step 24 (to reduce distortions and/or improve the SNR).
• Conventional Bo eddy current correction, concomitant field correction etc. may be applied in all three acquisition steps 21, 22, 24.
• Fig. 3 shows a B o map (in three planes) derived according to the invention from a survey scan and a calibration scan in a brain imaging session.
• MR signal data were acquired at 3T.
• the echo time of the calibration scan was 0.78 ms while the echo time of the survey scan was 1.4 ms.
• the calibration and survey scans had different fields of view and different image resolutions.
• the MR signal data of the two pre-scans were brought to the same field of view and resolution during processing, i.e. prior to the application of the Dixon algorithm.
• the determined B o values vary in a range of approx. 600 Hz and follow the pattern expected in the brain.

Landscapes

• Physics & Mathematics (AREA)
• Condensed Matter Physics & Semiconductors (AREA)
• General Physics & Mathematics (AREA)
• High Energy & Nuclear Physics (AREA)
• Engineering & Computer Science (AREA)
• Signal Processing (AREA)
• Health & Medical Sciences (AREA)
• General Health & Medical Sciences (AREA)
• Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
• Radiology & Medical Imaging (AREA)
• Magnetic Resonance Imaging Apparatus (AREA)

Applications Claiming Priority (2)

Publications (1)

Family

ID=83438323

Family Applications (2)

Family Applications Before (1)

Country Status (4)

Family Cites Families (1)

• 2022
  • 2022-09-23 EP EP22197279.7A patent/EP4343356A1/de not_active Withdrawn
• 2023
  • 2023-09-19 JP JP2025514739A patent/JP2025530297A/ja active Pending
  • 2023-09-19 WO PCT/EP2023/075702 patent/WO2024061852A1/en not_active Ceased
  • 2023-09-19 CN CN202380067468.XA patent/CN119895280A/zh active Pending
  • 2023-09-19 EP EP23772275.6A patent/EP4591084A1/de active Pending

Also Published As

Similar Documents

Legal Events