MR imaging with spiral acquisition

Abstract

The present invention relates to a method for MR imaging of an object (10) placed in an examination volume of an MR device (1). The aim of the invention is to enable efficient spiral MR imaging without blur artifacts, even in the case of strong B0 inhomogeneities. The method of the present invention comprises the steps of subjecting a subject 10 to an imaging sequence including at least one RF excitation pulse and modulated magnetic field gradients, acquiring MR signals along two or more planar spiral k-space trajectories 31, 32, 33, where the radial k-space velocity, i.e., the rate of variation of the radial distance from the spiral origin, is essentially constant along each planar spiral k-space trajectory, and the two or more k-space trajectories 31, 32, 33 are offset in-plane from one another, and reconstructing an MR image from the acquired MR signals. Furthermore, the present invention relates to an MR apparatus computer program.

Description

[Technical Field]

[0001] The present invention relates to the field of magnetic resonance (MR) imaging. The present invention relates to a method for MR imaging of an object. The present invention further relates to an MR device and to a computer program running on an MR device. [Background technology]

[0002] MR imaging methods, which utilize the interaction between magnetic fields and atomic nuclear spins to form two- or three-dimensional images, are widely used today, especially in the field of medical diagnostics, because they are superior in many ways to other imaging methods for imaging soft tissues, do not require ionizing radiation, and are usually non-invasive.

[0003] According to a typical MR method, the subject, e.g., the body of a patient to be examined, is placed in a strong, uniform magnetic field, the direction of which defines the axis (usually the z-axis) of the coordinate system on which the measurements are based. The magnetic field generates different energy levels for individual nuclear spins, depending on the field strength (spin resonance), which can be excited by the application of an alternating electromagnetic field (RF field) of a defined frequency (the so-called Larmor frequency or MR frequency). From a macroscopic perspective, the distribution of individual nuclear spins generates an overall magnetization, which can be deflected out of equilibrium by the application of an electromagnetic pulse (RF pulse) of an appropriate frequency, resulting in a precession of the magnetization around the z-axis. The precession describes the surface of a cone whose opening angle is called the flip angle. The magnitude of the flip angle depends on the strength and duration of the applied electromagnetic pulse. In the case of a so-called 90° pulse, the spins are deflected in a plane transverse to the z-axis (flip angle 90°).

[0004] After the RF pulse ends, the magnetization returns to its original equilibrium state, where the magnetization in the z direction reestablishes with a first time constant T1 (spin-lattice or longitudinal relaxation time), while the magnetization in the direction perpendicular to the z direction relaxes with a second time constant T2 (spin-spin or transverse relaxation time). The magnetization fluctuations can be detected by a receive RF coil positioned within the examination volume of the MR system and oriented so that the magnetization fluctuations are measured in a direction perpendicular to the z axis. The decay of the transverse magnetization, for example, after the application of a 90° pulse, involves the transition of nuclear spins (caused by inhomogeneities in the local magnetic field) from an ordered state with the same phase to a state in which all phase angles are uniformly distributed (dephasing). The dephasing can be compensated for, for example, by a refocusing pulse (e.g., a 180° pulse). This generates an echo signal (spin echo) in the receive coil.

[0005] To achieve spatial resolution within the body, magnetic field gradients extending along the three principal axes are superimposed on a uniform magnetic field, resulting in a linear spatial dependence of the spin resonance frequency. The signals picked up in the receive coils then contain different frequency components that can be associated with different locations within the body. The signal data obtained through the receive coils corresponds to the spatial frequency domain and is called k-space data. The k-space data set is converted into an MR image by an image reconstruction algorithm. Summary of the Invention [Problem to be solved by the invention]

[0006] Spiral imaging is a high-speed MR imaging technique that benefits from efficient k-space coverage and low sensitivity to motion and flow artifacts. Spiral k-space trajectories enable efficient and time-flexible sampling of k-space because shorter paths are required to cover the desired k-space region, and signal acquisition can begin at the center of k-space. However, spiral imaging techniques are vulnerable to amplitude inhomogeneities in the main magnetic field B0, which cause blurring and degrade image quality. For example, images obtained by spiral MR imaging of the brain are commonly contaminated by off-resonance signal contributions from the sagittal sinus and nasal cavity. Due to strong field inhomogeneities, the shape of a true spiral k-space trajectory deviates significantly from the theoretical spiral shape, and therefore, in practice, usable signal data may not be sampled from certain regions of k-space. Magnetic field inhomogeneities induced by patient anatomy can produce local field gradients of up to 0.5 μT / mm. In such cases, conventional single-shot spiral imaging (sampling the complete k-space region of interest along one spiral k-space trajectory after a single RF excitation) with a readout time of approximately 25 ms results in images exhibiting artifacts that cannot be resolved by known deblurring methods that apply post-processing phase correction. It is known that shortening the readout time, in combination with interleaved multi-shot trajectories, can reduce blur artifacts to some extent if necessary to obtain sufficient k-space coverage. However, this significantly reduces the efficiency of the method.

[0007] US Patent Application Publication No. 2014 / 0218028 discloses a spin-echo implementation with a single Archimedean spiral sampling orbit in k-space with a constant radial velocity.

[0008] From the above, it can be readily appreciated that there is a need for improved MR imaging techniques. It is an object of the present invention to address the above-mentioned limitations and enable efficient spiral MR imaging with reduced artifacts, even in the case of strong B0 inhomogeneities. [Means for solving the problem]

[0009] According to the present invention, a method for MR imaging of an object positioned within an examination volume of an MR device is disclosed. The method includes the steps of subjecting the object to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals along two or more planar spiral k-space trajectories, where the radial k-space velocity, i.e., the rate of variation of the radial distance from the spiral origin, is essentially constant along each planar spiral k-space trajectory, and the two or more k-space trajectories are offset in-plane from one another; and reconstructing an MR image from the acquired MR signals. This method is based on the insight that artifacts caused by local B0 inhomogeneities are related to the radial k-space velocity of the planar spiral k-space acquisition. This artifact occurs in the presence of magnetic field inhomogeneities combined with too slow a progression of sampling along the planar spiral k-space trajectory from the k-space center outward to its periphery. The sampling rate along the planar spiral k-space trajectory is limited toward the k-space periphery (due to limitations of the gradient hardware of the MR device used). Therefore, conventional spiral sampling strategies (typically using an Archimedean spiral trajectory) require the spiral's radial velocity to decelerate with increasing distance from the k-space origin in order to obtain sufficiently dense k-space coverage. The present invention instead proposes applying a substantially constant radial k-space velocity along the entire spiral trajectory. According to the present invention, since the progression speed along the trajectory is not decelerated, the sampling density inevitably decreases toward the k-space periphery. Nevertheless, to ensure sufficient k-space density, the present invention proposes dividing the scan into two or more in-plane k-space spirals offset from each other. According to the present invention, as the k-space sampling progresses rapidly outward from the k-space center along the spiral trajectory, blurring artifacts in the presence of main magnetic field inhomogeneities are significantly reduced without reducing the efficiency of the method.

[0010] The present invention relates to a method for MR imaging of an object, including sampling k-space along multiple spiral trajectories. In an outward sampling implementation, k-space sampling is performed along a current spiral trajectory that begins at a predetermined initial k-space position and continues outward to a current final k-space position where the radial k-space sampling density is below a current predetermined threshold. Then, starting from the next initial k-space position, a next spiral trajectory begins near the current final k-space position and continues outward to a current final k-space position where the radial k-space sampling density is below the next threshold. Thus, over the multiple concatenated spiral trajectories, the radial k-space velocity remains above the current and next predetermined thresholds, limited by the capabilities of the gradient system. In the succession of spiral trajectories, the radial k-space sampling density does not fall below the current and next predetermined thresholds, thereby remaining within a range limited by nominal values ​​related to the thresholds and the capabilities of the gradient system. By selecting the current and next thresholds within a narrow range below the nominal radial k-space velocity related to the gradient system's capabilities in terms of gradient strength and slew rate, the radial k-space sampling density remains essentially constant over the concatenated k-space trajectories. The first of the spiral trajectories can start at (near) the origin of k-space. In an actual implementation, the next initial k-space position either coincides with the current final k-space position or is separated from it by a sampling distance. The current and next predetermined thresholds may be equal so that the residual variation in radial k-space velocity is the same from one spiral trajectory to the next.

[0011] In an inward sampling implementation, the current spiral trajectory starts from a start k-space position and continues inward to a current end k-space position. At the start k-space position, the radial k-space sampling density is at least a current threshold value, and the radial k-space sampling density increases to a current ceiling value at the current end k-space position. Then, starting from a next start k-space position near the current end k-space position, the next spiral trajectory begins toward the next end k-space position. At the next start k-space position, the radial k-space sampling density is at least the next threshold value and also increases to the next ceiling value. In an actual implementation, the next start k-space position coincides with the current end k-space position or is a sampling distance away from the current end k-space position. The current predetermined ceiling value and the next predetermined ceiling value can be equal so that the residual variation in the radial k-space sampling density is the same from the current spiral trajectory to the next spiral trajectory. These ceiling values ​​can actually be related to the capabilities of the gradient system.

[0012] The k-space sampling along the spiral trajectory can be achieved iteratively, such that after the next spiral trajectory, a further next spiral trajectory is traversed, etc. The more different spiral trajectories are used, the smaller the variation in radial k-space sampling density over the succession of spiral trajectories, or the larger the region of k-space that is sampled with a given (narrow) variation range of radial k-space sampling density.

[0013] A magnetic resonance image can be reconstructed from the acquired magnetic resonance signals. The reconstruction can be implemented as a reconstruction module in software installed on a host computer of the magnetic resonance examination system. Alternatively, the reconstruction can be performed remotely, for example in a cloud at a location and time different from the magnetic resonance examination system and the acquisition of the MR signals. That is, the present invention can include configuring a magnetic resonance image to be reconstructed from magnetic resonance signals acquired by sampling k-space along a spiral trajectory.

[0014] In a preferred embodiment, MR signals are acquired along two or more planar spiral k-space trajectories after a single RF excitation pulse. This corresponds to a single-shot implementation of the method of the present invention. Two or more planar spiral k-space trajectories are acquired one after the other after a single RF excitation. Therefore, in terms of scan efficiency, the method of the present invention is identical to conventional single-shot spiral imaging techniques.

[0015] In another preferred embodiment, the starting points of two or more planar spiral k-space trajectories (i.e., their respective starting points in k-space) are offset in-plane from each other. Different origins of the two or more spiral trajectories can also be used to obtain the desired sufficiently dense k-space coverage. The offset origins in-plane result in the two or more spiral trajectories being interleaved in k-space, such that sampling gaps of one trajectory are filled by another trajectory.

[0016] In yet another possible embodiment, two or more planar spiral k-space trajectories are counter-wound relative to each other about their spiral axes, which also allows for an interleaved set of planar spiral k-space trajectories with sufficiently dense k-space coverage.

[0017] According to a further embodiment, at least one of the planar spiral k-space trajectories is offset in-plane from the k-space origin. It is sufficient that only one of the two or more k-space trajectories starts at the k-space origin. The starting points of the other spiral trajectories can be freely chosen to optimize the sampling density around the k-space.

[0018] In one embodiment of the method of the present invention, the (radial) distance between turns of the planar spiral k-space trajectory increases with increasing distance from the spiral origin. This is the result of a combination of a constant radial velocity of the spiral k-sampling and the inevitable gradient limitation. Therefore, the spiral trajectories used in accordance with the present invention are generally not of the Archimedean type.

[0019] According to yet another preferred embodiment, a B0 map is derived by comparing MR signals acquired along various planar spiral k-space trajectories. An MR image can be reconstructed from the MR signals of each planar spiral k-space trajectory (ideally limited to a central k-space region with sufficient sampling density). The B0 map can be derived by comparing different phases of the reconstructed images. Then, a final MR image can be reconstructed with B0 inhomogeneity corrected based on the derived B0 map. Similarly, the effect of T2* can be addressed. Magnetic resonance signals along different spiral trajectories are affected differently by T2* relaxation. T2* can be determined by comparing the MR signals of the different spiral trajectories and then used to correct the MR signals accordingly. In the sense of the present invention, a planar spiral k-space trajectory is defined as a curve that winds around the origin at continuously increasing or decreasing distances. This definition covers (among other things) spirals in three-dimensional k-space, such as a conical spiral that winds around an axis, and spirals in two-dimensional k-space, which are curves in a plane that wind around a fixed central point at continuously increasing (outward spiral) or decreasing (inward spiral) distances from the central point, i.e., the central point constitutes the origin, or starting point, of the spiral.

[0020] The MR signals can be acquired, for example, as free induction decay (FID) immediately after the emission of each RF excitation pulse. The MR signals can also be acquired as spin echo or gradient echo signals.

[0021] The inventive method described above can be performed by an MR apparatus having at least one main magnet coil for generating a homogeneous static magnetic field in an examination volume, multiple gradient coils for generating switched magnetic field gradients in different spatial directions in the examination volume, at least one RF coil for generating RF pulses in the examination volume and / or receiving MR signals from a subject arranged 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 inventive method can be realized, for example, by a corresponding program in the reconstruction unit and / or the control unit of the MR apparatus.

[0022] The method of the present invention can be advantageously implemented in most MR devices currently in clinical use. For this purpose, it is only necessary to utilize a computer program that controls the MR device to perform the above-described method steps of the present invention. The computer program may be present on a data carrier or in a data network, so as to be downloaded for installation in the control unit of the MR device.

[0023] The accompanying drawings disclose preferred embodiments of the present invention. It is to be understood, however, that the drawings are designed for purposes of illustration only and not as a definition of the limits of the invention. [Brief explanation of the drawings]

[0024] [Figure 1] 1 is a diagram showing an MR apparatus for carrying out the method of the present invention. [Figure 2] 1 is a diagram of k-space illustrating the spiral sampling technique of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0025] Referring to Figure 1, there is shown a schematic representation of an MR system 1. The system has superconducting or resistive main magnet coils 2 such that a substantially uniform and temporally constant main magnetic field is generated along the z-axis through an examination volume.

[0026] The magnetic resonance generation and manipulation system applies a series of RF pulses and switching magnetic field gradients to invert or excite nuclear magnetic spins, induce magnetic resonance, refocus magnetic resonance, manipulate magnetic resonance, spatially and otherwise encode magnetic resonance, saturate spins, and the like to perform MR imaging.

[0027] More specifically, gradient pulse amplifiers 3 apply current pulses to selected ones of whole-body gradient coils 4, 5, and 6 along the x-, y-, and z-axes of the examination volume. A digital RF frequency transmitter 7 transmits RF pulses or pulse packets to whole-body volume RF coil 9 via transmit / receive switch 8 to transmit RF pulses into the examination volume. A typical MR imaging sequence consists of packets having multiple short-duration RF pulse segments that, together with any applied magnetic field gradients, achieve selected manipulation of nuclear magnetic resonance. The RF pulses are used to saturate resonance, excite resonance, invert magnetization, refocus resonance, or manipulate resonance, as well as to select portions of a body 10 located within the examination volume. MR signals are also picked up by the whole-body volume RF coil 9.

[0028] To generate MR images of a limited region of the body 10, a set of local array RF coils 11, 12, 13 are positioned adjacent to the region selected for imaging. The array coils 11, 12, 13 can be used to receive MR signals caused by body coil RF transmissions.

[0029] The resulting MR signals are picked up by a whole-body volume RF coil 9 and / or array RF coils 11, 12, 13 and demodulated by a receiver 14, preferably having a preamplifier (not shown). The receiver 14 is connected to the RF coils 9, 11, 12, and 13 via a transmit / receive switch 8.

[0030] A host computer 15 controls the gradient pulse amplifier 3 and the transmitter 7 to generate multiple MR imaging sequences, such as spin-echo imaging, and acquire MR signals along spiral k-space trajectories according to the present invention. For a selected sequence, the receiver 14 receives single or multiple MR data along each k-space trajectory. A data acquisition system 16 performs analog-to-digital conversion of the received signals and converts each MR signal into a digital format suitable for further processing. In modern MR devices, the data acquisition system 16 is a separate computer dedicated to acquiring raw image data.

[0031] Ultimately, the digital raw image data is reconstructed into an image representation by a reconstruction processor 17, which applies interpolation or re-gridding of the data from the spiral acquisitions before a Fourier transform or other appropriate reconstruction algorithm. The MR image may represent a planar slice across the patient, an array of parallel planar slices, a three-dimensional volume, etc. The image is then stored in an image memory, which can be accessed to convert slices, projections, or other portions of the image representation into an appropriate format for visualization, for example, via a video monitor 18, which provides a human-readable display of the resulting MR image.

[0032] The MR apparatus 1 is configured, for example, by suitable programming of a host computer 15 and a reconstruction processor 17 to carry out the imaging method of the invention as described above and below.

[0033] With continuing reference to FIG. 1 and further reference to FIG. 2, an embodiment of the imaging approach of the present invention will now be described.

[0034] FIG. 2 illustrates the spiral acquisition strategy of the present invention. The present invention proposes the application of two or more planar spiral k-space trajectories 31, 32, 33, where the radial k-space velocity, i.e., the rate of variation of the radial distance from the spiral origin, is essentially constant along each planar spiral k-space trajectory, and the two or more k-space trajectories 31, 32, 33 are offset in-plane from one another. The first planar spiral k-space trajectory is designated by reference numeral 31. At the k-space center, which coincides with the origin of the spiral 31, the angular velocity of sampling along the trajectory is high. However, the progression speed along the spiral trajectory 31 is limited to a reduced angular velocity (due to the limited capabilities of the gradient systems 3, 4, 5, 6 in terms of strength and slew rate). If the spiral trajectory is of the Archimedean type (as in the prior art), the radial velocity is also reduced. According to the present invention, the radial velocity of sampling along the spiral trajectory 31 is kept substantially constant, so that the shape of the trajectory 31 deviates from an Archimedean spiral, resulting in a lower sampling density at the periphery of k-space. Here, the sampling density is too low to obtain the desired quality of the final reconstructed MR image. To compensate for the reduced density, the sampling is divided into two or more (in the illustrated embodiment, three) planar spiral k-space trajectories 31, 32, 33. After completion of sampling along the first trajectory 31 (which takes 10 ms in the illustrated example), sampling transitions to the second trajectory 32 (e.g., in 0.7 ms, as indicated by the dashed curve). In this example, the second trajectory 32 does not start at the origin of k-space (kx = ky = 0), but rather at a point where the sampling density of the first trajectory 31 falls below a given threshold, e.g., passing the first trajectory 31 3 ms after the start of the scan. Thus, the second trajectory requires only 7 ms to reach the outer limit of the required k-space region. The second trajectory 32 then transitions to the start of a third trajectory 33 which takes 5 ms to reach the outer limit of the required k-space region. The illustrated planar spiral k-space sampling scheme achieves a higher and substantially constant density of spiral arms 31, 32, 33 than conventionally applied single spiral k-space trajectories, while keeping the radial k-space velocity constant.In this way, the level of blur artifacts caused by B0 inhomogeneities in the final reconstructed magnetic resonance image can be reduced without significantly adversely affecting sampling efficiency.

[0035] In the presence of B0 inhomogeneity, MR signals acquired along different planar spiral k-space trajectories 31, 32, 33 accumulate different amounts of phase. Furthermore, the MR signals of the different k-space trajectories 31, 32, 33 experience different T2* effects. Intermediate MR images can be reconstructed from MR signals acquired from central k-space portions along the individual trajectories 31, 32, 33. By comparing these images, a B0 map and T2* can be obtained. Knowing these parameters, a final MR image can be reconstructed from the MR signals of all three planar spiral k-space trajectories 31, 32, 33. The following describes embodiments of the present invention. (Appendix 1) 1. A method for MR imaging of an object positioned within an examination volume of an MR device, comprising: subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing with a radial k-space velocity, wherein a current spiral trajectory is traversed outward from a predetermined current initial k-space location to a current final k-space location where the radial k-space velocity is below a predetermined current threshold, and a next spiral trajectory is traversed outward from a next initial k-space location to a current final k-space location where the radial k-space velocity is below a predetermined next threshold; A method having the following. (Appendix 2) 1. A method for MR imaging of an object positioned within an examination volume of an MR device, comprising: subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing with a radial k-space velocity, wherein a current spiral trajectory is traversed inward from a predetermined starting k-space location where the radial k-space velocity is at least a current threshold value to a current ending k-space location where the radial k-space velocity increases to a current ceiling value, and a next spiral trajectory is traversed inward from a next starting k-space location near the current ending k-space location where the radial k-space velocity is at least a next threshold value to a next ending k-space location where the radial k-space velocity reaches a next ceiling value; A method having the following. (Appendix 3) 2. The method of claim 1, wherein the MR signals are acquired along the two or more planar spiral k-space trajectories after a single RF excitation pulse. (Appendix 4) 4. The method of claim 3, wherein the origins of the two or more planar spiral k-space trajectories are offset from one another in a plane. (Appendix 5) 5. The method of any one of claims 1 to 4, wherein the two or more planar spiral k-space trajectories are counter-wound relative to each other about their spiral axes. (Appendix 6) 6. The method of any one of claims 1 to 5, wherein the two or more planar spiral k-space trajectories are interleaved with one another. (Appendix 7) 7. The method of any one of claims 1 to 6, wherein at least one of the planar spiral k-space trajectories is offset in-plane from the k-space origin. (Appendix 8) 8. The method of any one of claims 1 to 7, wherein the distance between turns of the planar spiral k-space trajectory increases with increasing distance from the spiral origin. (Appendix 9) 9. The system of any one of claims 1 to 8, wherein a B0 map is derived by comparing the MR signals acquired along the planar spiral k-space trajectory. (Appendix 10) 10. The system of claim 9, wherein the MR image is reconstructed by correcting B0 inhomogeneity based on the derived B0 map. (Appendix 11) 1. An MR apparatus comprising: at least one main magnet coil for generating a homogeneous static magnetic field in an examination volume; several gradient coils for generating switched magnetic field gradients in different spatial directions in the examination volume; at least one RF coil for generating RF pulses in the examination volume and / or receiving MR signals from a subject positioned in the examination volume; a control unit for controlling the time-sequential RF pulses and the switched magnetic field gradients; and a reconstruction unit for reconstructing MR images from the received MR signals, subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing a radial k-space velocity, wherein a current spiral trajectory is traversed outward from a current initial k-space location to a current final k-space location where the radial k-space velocity is below a predetermined current threshold, and a next spiral trajectory is traversed outward from a next initial k-space location to a next final k-space location where the radial k-space velocity is below a predetermined next threshold; or subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing with a radial k-space velocity, wherein a current spiral trajectory is traversed inward from a predetermined start k-space location where the radial k-space velocity is at least a current threshold value to a current end k-space location where the radial k-space velocity increases to a current ceiling value, and a next spiral trajectory is traversed inward from a next start k-space location near the current end k-space location where the radial k-space velocity is at least a next threshold value to a next end k-space location where the radial k-space velocity reaches a next ceiling value; The MR apparatus is configured to perform the steps of: (Appendix 12) A computer program executed on an MR device, subjecting the subject to an imaging sequence having at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing with a radial k-space velocity, wherein a current spiral trajectory is traversed outward from a current initial k-space location to a current final k-space location where the radial k-space velocity is below a predetermined current threshold, and a next spiral trajectory is traversed outward from a next initial k-space location to a next final k-space location where the radial k-space velocity is below a predetermined next threshold; or subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing with a radial k-space velocity, wherein a current spiral trajectory is traversed inward from a predetermined start k-space location where the radial k-space velocity is at least a current threshold value to a current end k-space location where the radial k-space velocity rises to a current ceiling value, and a next spiral trajectory is traversed inward from a next start k-space location near the current end k-space location where the radial k-space velocity is at least a next threshold value to a next end k-space location where the radial k-space velocity reaches a next ceiling value; 1. A computer program comprising:

Claims

1. 1. A method of MR imaging of an object positioned within an examination volume of an MR device, comprising: subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing at a constant radial k-space velocity, wherein a current spiral trajectory is traversed outward from a predetermined current initial k-space location to a current final k-space location where the radial k-space sampling density is below a predetermined current threshold, and a next spiral trajectory is traversed outward from a next initial k-space location to a next final k-space location where the radial k-space sampling density is below a next predetermined threshold, and at least one of the planar spiral k-space trajectories is offset in-plane from the k-space origin to increase the sampling density around the k-space; A method having the following.

2. 1. A method of MR imaging of an object positioned within an examination volume of an MR device, comprising: subjecting the subject to an imaging sequence including at least one RF excitation pulse and a modulated magnetic field gradient; acquiring MR signals by sampling k-space along a plurality of planar spiral k-space trajectories within k-space while traversing at a constant radial k-space velocity, wherein a current spiral trajectory is traversed inward from a predetermined starting k-space location where the radial k-space sampling density is at least a current threshold value to a current ending k-space location where the radial k-space sampling density increases to a current ceiling value, and a next spiral trajectory is traversed inward from a next starting k-space location near the current ending k-space location where the radial k-space sampling density is at least a next threshold value to a next ending k-space location where the radial k-space sampling density reaches a next ceiling value, and at least one of the planar spiral k-space trajectories is offset in-plane from the k-space origin to increase the sampling density around k-space; A method having the following.

3. The method of claim 1 , wherein the MR signals are acquired along two or more of the planar spiral k-space trajectories after a single RF excitation pulse.

4. The method of claim 3 , wherein the origins of two or more of the planar spiral k-space trajectories are offset from one another in a plane.

5. 5. The method of claim 1, wherein two or more of the planar spiral k-space trajectories are counter-wound relative to each other about their spiral axes.

6. 6. The method of claim 1, wherein two or more of the planar spiral k-space trajectories are interleaved with each other.

7. 7. The method of claim 1, wherein the distance between turns of the planar spiral k-space trajectory increases with increasing distance from the origin of the trajectory.

8. 8. The method of claim 1, wherein a B0 map is derived by comparing the MR signals acquired along the planar spiral k-space trajectory.

9. The method of claim 8 , wherein an MR image is reconstructed by correcting for B0 inhomogeneity based on the derived B0 map.

10. 3. An MR apparatus comprising: at least one main magnet coil for generating a homogeneous static magnetic field in an examination volume; several gradient coils for generating switched magnetic field gradients in different spatial directions in the examination volume; at least one RF coil for generating RF pulses in the examination volume and / or receiving MR signals from a subject positioned in the examination volume; a control unit for controlling the temporally successive RF pulses and the switched magnetic field gradients; and a reconstruction unit for reconstructing MR images from the received MR signals, wherein the MR apparatus is configured to perform the method according to claim 1 or 2.

11. A computer program running on an MR device, causing a computer of the MR device to carry out the method according to claim 1 or 2.

Images