Magnetic resonance imaging apparatus and method

Abstract

To accelerate magnetic resonance imaging.SOLUTION: A magnetic resonance imaging apparatus comprises a generation unit, an acquisition unit, and a reconstitution unit. The generation unit generates a pulse sequence to collect magnetic resonance imaging (MRI) data by applying a high frequency pulse and a slice selection gradient magnetic field in the same period every repeated time, the pulse sequence being provided with time deviation between the center of the high frequency pulse and the center of the waveform of the slice selection gradient magnetic field so as to generate an offset of a space moment in a slice encoding direction, every repeated time mentioned above. The acquisition unit collects the MRI data by using the pulse sequence. The reconstitution unit reconstitutes an MRI image from the MRI data.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] SUMMARY OF THE INVENTION The embodiments disclosed herein and in the drawings relate to magnetic resonance imaging apparatus and methods. [Background technology]

[0002] The "Background" set forth herein is provided for the purpose of generally providing a context for the present disclosure. Work by the inventors named herein, to the extent described in this Background section, together with aspects of this specification that are not otherwise considered prior art as of the filing date, are not admitted expressly or impliedly as prior art to the present disclosure.

[0003] Magnetic resonance imaging (MRI) is an imaging scanning method in which the nuclear spins of a subject placed in a static magnetic field are magnetically excited by radio frequency (RF) pulses at the Larmor frequency, and an image is generated from the magnetic resonance signal data generated by this excitation.

[0004] MRI can be accelerated using the Controlled Aliasing In Parallel Imaging (CAIPI) technique. For example, the number of samples can be reduced by thinning and zero-padding the k-space sampling points. However, thinning the k-space sampling points can cause aliasing in the spatial domain, which can be difficult to separate. By introducing the CAIPI technique to stagger samples in k-space, aliasing in accelerated data can be more easily eliminated. Current methods achieve staggered sampling in k-space by modulating RF pulses (e.g., phase modulation) or by modulating magnetic field gradients (e.g., amplitude modulation of phase-encoding gradients).

[0005] Other methods use k-space-based methods for reconstruction of CAIPI acquisitions, such as Sensitivity Encoding (SENSE), Generalized Autocalibrating Partial Parallel Acquisition (GRAPPA), slice GRAPPA, and autocalibrated CAIPI.

[0006] It is therefore desirable to improve the method for shifting the moments of samples in k-space, which can improve the simplicity and robustness of the CAIPI technique. [Prior art documents] [Patent documents]

[0007] [Patent Document 1] U.S. Patent No. 7,002,344 [Patent Document 2] U.S. Patent No. 8,405,395 [Patent Document 3] US Patent Application Publication No. 2019 / 0120922 [Patent Document 4] U.S. Patent No. 9,989,610 [Patent Document 5] US Patent Application Publication No. 2015 / 0177353 [Patent Document 6] US Patent Application Publication No. 2018 / 0164395 [Patent Document 7] European Patent Application Publication No. 3413074 [Patent Document 8] Chinese Patent Application Publication No. 108957375 [Non-patent literature]

[0008] [Non-Patent Document 1] Felix A. Breuer, Martin Blaimer, Robin M. Heidemann, Matthias F. Mueller, Mark A. Griswold, Peter M. Jakob, “Controlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA) for Multi-Slice Imaging”, Magnetic Resonance in Medicine, 53:684-691 (2005) [Non-patent document 2] Kawin Setsompop, Borjan A. Gagoski, Jonathan R. Polimeni, Thomas Witzel, Van J. Wedeen, Lawrence L. Wald, “Blipped-Controlled Aliasing in Parallel Imaging for Simultaneous Multislice Echo Planar Imaging With Reduced g-Factor Penalty”, Magnetic Resonance in Medicine.67:1210-1224(2012) Summary of the Invention [Problem to be solved by the invention]

[0009] One of the problems to be solved by the embodiments disclosed in this specification and the drawings is to accelerate magnetic resonance imaging. However, the problems to be solved by the embodiments disclosed in this specification and the drawings are not limited to the above problem. Problems corresponding to the effects of each configuration shown in the embodiments described below can also be positioned as other problems. [Means for solving the problem]

[0010] A magnetic resonance imaging apparatus according to an embodiment includes a generator, an acquirer, and a reconstructor. The generator generates a pulse sequence for acquiring MRI (Magnetic Resonance Imaging) data by applying a radio frequency pulse and a slice-selective gradient magnetic field for the same period for each repetition time, with a time shift provided between the center of the radio frequency pulse and the center of the waveform of the slice-selective gradient magnetic field so as to generate a spatial moment offset in the slice encoding direction for each repetition time. The acquirer acquires the MRI data using the pulse sequence. The reconstructor reconstructs an MRI image from the MRI data. [Brief explanation of the drawings]

[0011] The present disclosure and many of the attendant advantages thereof will become more fully understood as the same becomes better understood by reference to the following detailed description, when considered in conjunction with the accompanying drawings. [Figure 1] FIG. 1 is a flowchart of a method for generating CAIPI MRI data by time-shifting waveforms in an MRI pulse sequence and reconstructing an MRI image from the MRI data, according to one embodiment. [Figure 2] FIG. 2 is a diagram illustrating a first example of an MRI apparatus according to one embodiment. [Figure 3A] FIG. 3A illustrates a full sampling pattern in the ky direction (phase encoding direction) and kz direction (slice encoding direction) (i.e., a sampling pattern of Ry=1 and Rz=1 with no acceleration in either the ky direction or the kz direction) according to one embodiment. [Figure 3B] FIG. 3B illustrates an image of a cone that is scanned to generate an MRI image, according to one embodiment. [Figure 3C] FIG. 3C illustrates sagittal and axial images produced when a Fourier transform is applied to MRI data acquired according to the sampling pattern of FIG. 3A, according to one embodiment. [Figure 3D]FIG. 3D illustrates a sampling pattern accelerated in the ky direction (i.e., Ry=2 and Rz=1) according to one embodiment. [Figure 3E] FIG. 3E illustrates sagittal and axial images produced when a Fourier transform is applied to MRI data acquired according to the sampling pattern of FIG. 3D, according to one embodiment. [Figure 3F] FIG. 3F illustrates a kz-accelerated (i.e., Ry=1 and Rz=2) sampling pattern according to one embodiment. [Figure 3G] FIG. 3G illustrates sagittal and axial images produced when a Fourier transform is applied to MRI data acquired according to the sampling pattern of FIG. 3F, according to one embodiment. [Figure 3H] FIG. 3H illustrates a CAIPI sampling pattern accelerated by 2x, according to one embodiment. [Figure 3I] FIG. 3I illustrates sagittal and axial images produced when Fourier transforming MRI data acquired according to the sampling pattern of FIG. 3H, according to one embodiment. [Figure 3J] FIG. 3J illustrates a process for generating the CAIPI sampling pattern of FIG. 3H from the sampling pattern of FIG. 3F by adding a moment Δkz in the kz direction every other repetition time (TR), according to one embodiment. [Figure 4A] FIG. 4A illustrates a timing sequence of a pulse sequence that generates field echoes without time shifting to generate a CAIPI sampling pattern, according to one embodiment. [Figure 4B] FIG. 4B shows a timing diagram of a pulse sequence in which the RF pattern is time-shifted to generate a CAIPI sampling pattern, according to one embodiment. [Figure 4C]FIG. 4C shows a timing diagram of a pulse sequence in which the slice-selective gradient waveforms are time-shifted to generate a CAIPI sampling pattern, according to one embodiment. [Figure 5] FIG. 5 illustrates MRI data acquired with the CAIPI technique using a reduced field of view (FOV), and the SENSE method used to eliminate aliasing in the reconstructed image, according to one embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0012] The following description, set forth in connection with the accompanying drawings, is intended to illustrate various embodiments of the disclosed subject matter and is not necessarily intended to represent the only one or more embodiments. In some cases, the description includes specific details for the purpose of providing an understanding of the disclosed subject matter. However, it will be apparent to those skilled in the art that embodiments may be practiced without these specific details. In some instances, well-known structures and components are shown in block diagram form to avoid obscuring the concepts of the disclosed subject matter.

[0013] References throughout this specification to "one embodiment" or "embodiment" mean that a particular feature, structure, characteristic, scan, or function described in connection with an embodiment is included in at least one embodiment of the disclosed subject matter. Thus, any appearances of "in one embodiment" or "in an embodiment" herein do not necessarily refer to the same embodiment. Moreover, the particular features, structures, characteristics, scans, or functions may be combined in any suitable manner in one or more embodiments. Furthermore, embodiments of the disclosed subject matter can and are intended to cover variations and modifications of the described embodiments.

[0014] While related methods achieve CAIPI by modulating either the phase of an excitation RF pulse or the amplitude of a gradient magnetic field waveform to achieve different sampling patterns in the phase-encode (PE) and slice-encode (SE) directions, the method described herein advantageously uses time delays to achieve a desired k-space sampling pattern without modulating the RF pulse or the gradient magnetic field waveform. This CAIPI technique can accelerate MRI data acquisition without significantly degrading image quality by minimizing residual signals resulting from aliasing. Note that, in this specification, a time lag is also referred to as a "time shift," and a shift in time is also referred to as "time shifting" or "shifting in time."

[0015] The method described herein is not limited to using the SENSE method. For example, other parallel imaging methods such as the GRAPPA method or the ARC (Autocalibrating Reconstruction for Cartesian Imaging) method may also be used. Below, this method is illustrated and described using a non-limiting example in which separate images of slices are reconstructed using the SENSE method.

[0016] Referring now to the drawings, where like reference numerals indicate the same or corresponding parts in the several views, FIG. 1 shows a flowchart of a method 50 for generating CAIPI MRI data using a time shift to adjust the k-space sampling pattern of signals acquired in the PE and SE directions and reconstructing an image from the CAIPI MRI data. In steps 10 and 20, a field echo (FE) pulse sequence (hereinafter, FE sequence) 25 is generated and modified to sample two phase-encoding dimensions with an acceleration factor R using a CAIPI sampling pattern. Here, the CAIPI sampling pattern is achieved using a time shift between RF pulses and slice-selective gradient waveforms, as described in detail below. In step 30, MRI data 35 is acquired using the generated FE sequence 25. In step 40, an MRI image 45 is reconstructed from the MRI data 35.

[0017] FIG. 2 shows a non-limiting example of an MRI apparatus 100. The MRI apparatus 100 shown in FIG. 2 includes a gantry 101 (shown in schematic cross section) and various associated MRI system components 103 that operate in conjunction with the gantry 101. Typically, at least the gantry 101 is placed in a shielded room. The structure of the MRI apparatus shown in FIG. 2 includes a static magnetic field B0 magnet 111 arranged in a generally coaxial cylindrical shape, a set of Gx, Gy, and Gz gradient coils 113, and a large RF whole-body coil (WBC) 115. An imaging volume 117 is shown along the horizontal axis of this cylindrically arranged element, substantially surrounding the head of a patient 119 supported by a patient couch 120.

[0018] Within the imaging volume 117, one or more small RF coils, known as array coils (ACs) 121, may be mounted closer to the patient's head (referred to herein, for example, as the "subject being scanned" or "subject"). As will be appreciated by those skilled in the art, relatively small coils and / or arrays, such as surface coils, compared to WBCs are often customized for specific body parts (e.g., arms, shoulders, elbows, wrists, knees, legs, chest, spine, etc.). Such small RF coils are referred to herein as ACs or PACs (Phased-Array Coils). These include at least one coil configured to transmit RF signals into the imaging volume and multiple receiver coils configured to receive RF signals from the subject, such as the patient's head, in the imaging volume.

[0019] The MRI apparatus 100 also includes an MRI system controller 130 having input / output ports connected to a display 124, a keyboard 126, and a printer 128. As will be appreciated, the display 124 may be a touch screen to also provide control input. Additionally, other I / O devices, such as a mouse, may be provided.

[0020] The MRI system controller 130 interfaces with an MRI sequence controller 140, which in turn controls the Gx, Gy, and Gz gradient coil drivers 132, the RF transmitter 134, and the transmit / receive switch 136 (if the same RF coil is used for both transmit and receive). The MRI sequence controller 140 includes a suitable program code structure 138 for implementing MRI (also known as nuclear magnetic resonance imaging or NMR imaging) techniques, including parallel imaging. The MRI sequence controller 140 can be configured to implement MRI with or without parallel imaging. Furthermore, the MRI sequence controller 140 can execute one or more preparatory scan (prescan) sequences and a scan sequence for acquiring a main scan MR image (referred to as a diagnostic image). For example, the MR data obtained in the prescan can be used to calculate sensitivity maps (also called coil sensitivity maps or spatial sensitivity maps) of the WBC 115 and / or AC 121, as well as unfolding maps for parallel imaging.

[0021] The MRI system component 103 includes an RF receiver 141 that provides input to a data processor 142 to generate processed image data that is sent to the display 124. The MRI data processor 142 is further configured to access a memory 146 containing previously generated MRI data, MR images, and / or maps, such as coil sensitivity maps, unfolding maps for parallel imaging, distortion maps, and / or system configuration parameters, and MRI image reconstruction program code structures 144 and 150.

[0022] In one embodiment, the MRI data processor 142 includes processing circuitry, which may include application-specific integrated circuits (ASICs), configurable logic devices (e.g., simple programmable logic devices (SPLDs), complex programmable logic devices (CPLDs), field programmable gate arrays (FPGAs), and other circuit components arranged to perform the functions described in this disclosure. Here, the MRI data processor 142 is an example of a generator, acquirer, and reconstructor.

[0023] The MRI data processor 142 executes one or more sequences of one or more instructions contained in the program code structures 144 and 150. Alternatively, the instructions may be read from another computer-readable medium, such as a hard disk or removable media. Alternatively, one or more processors in a multi-processing arrangement may be used to execute the sequences of instructions contained in the program code structures 144 and 150. In other embodiments, hardwired circuitry may be used in place of or in combination with software instructions. Thus, the disclosed embodiments are not limited to any particular combination of hardware circuitry and software. For example, the program code structure 150 may store instructions that are executed to implement the method 50.

[0024] Additionally, the term "computer-readable medium" as used herein refers to non-transitory storage media used to provide instructions to processor 142 for execution. Computer-readable media may take many forms, including but not limited to, non-volatile or volatile media. Non-volatile media include, for example, optical disks, magnetic disks, magneto-optical disks, or removable media drives. Volatile media include dynamic memory.

[0025] 2 further illustrates a generalized diagram of the MRI system program storage (memory) 150, in which program code structures, such as instructions for executing method 50, are stored on a computer-readable medium accessible to the various data processing components of MRI apparatus 100. As will be appreciated by those skilled in the art, program storage 150 may be partitioned and, at least in part, directly connected to the different processing computers included in MRI system components 103 that most immediately require the program code structures in their normal operation (i.e., rather than being normally stored and directly connected to MRI system controller 130).

[0026] Furthermore, the MRI device 100 shown in Figure 2 may be utilized to implement the exemplary embodiments described herein below, where the MRI system components 103 are divided into different logical collections of "boxes" and may generally include multiple Digital Signal Processors (DSPs), microprocessors, and special purpose processing circuits (e.g., high speed A / D conversion, fast Fourier transform, array processing, etc.). Each processor is typically a clocked "state machine" in which the occurrence of one clock cycle (or a predetermined number of clock cycles) causes the physical data processing circuitry to transition from one physical state to another.

[0027] Furthermore, not only do the physical states of processing circuitry (e.g., CPUs, registers, buffers, arithmetic units, etc.) change incrementally from one clock cycle to another, but the physical states of associated data storage media (e.g., bit storage sites of a magnetic storage medium) are transformed from one state to another during operation of such systems. For example, at the end of an image reconstruction process and / or sometimes the process of generating image reconstruction maps (e.g., coil sensitivity maps, unfolding maps, ghost maps, distortion maps, etc.), arrays of computer-readable and accessible data value storage areas in the physical storage media are transformed from one prior state to a new state, and the physical states of the physical sites of such arrays change between minimum and maximum values to represent real-world physical events and conditions. As will be appreciated by those skilled in the art, such arrays of stored data values represent and constitute physical structures, as well as the specific structures of computer control program code that, when sequentially loaded into instruction registers and executed by one or more CPUs of the MRI system 100, generates a particular sequence of operating states to transition within the MRI system 100.

[0028] The MRI apparatus 100 configured to perform the method 50 has several advantages. First, the slice-selective gradient field (hereinafter referred to as G SS By shifting the RF pulse in time relative to the waveform of the RF pulse, the complexity of modulating the waveform of the RF pulse or the waveform of the gradient magnetic field can be avoided. For example, SS A time-shifted pattern of 1000 s may be generated so that the waveforms and amplitudes of all gradient fields may remain the same amplitude and shape throughout the imaging process, even as the time position within the sequence changes.

[0029] Furthermore, G SS By using a pulse sequence in which the time shift of the RF pulse or G SSThis can eliminate the need to modulate the waveform of the RF pulse. For example, modulating the RF pulse can create complications if the RF pulse is distorted by nonlinearities in downstream electronics (e.g., amplifiers). Furthermore, G SS The use of a time-shifted pulse sequence helps to maintain the gradient amplifiers, gradient power supply, duty cycle, and acoustic noise over time. In other words, rather than adding anything extra to the gradient field, the G SS By shifting the center of the waveform in time, the MRI system 100 can be made more efficient and prevent reaching the limits of practical hardware.

[0030] Returning to FIG. 1, in step 10 of method 50, an initial FE sequence 15 is generated. For example, FIG. 4A shows an example of an FE sequence implementation. In this example, an excitation RF pulse excites the transverse magnetization of a given slice. That is, the RF pulse excites a flip angle α in each spin contained in that slice. Also, a G SS is applied, where G SS The waveform of FIG. 4A has a flat top portion (hereinafter referred to as a flat top), and the period of the flat top portion is longer than the period of the RF pulse so that the RF pulse can be shifted in time within the period. The FE sequence 25 also includes the application of one or more phase encoding gradient magnetic fields and readout gradient magnetic fields (e.g., frequency encoding gradient magnetic fields). In FIG. 4A, the first RF pulse occurs at time zero. In one embodiment, the time zero is defined to be the center of each of the first and second RF pulses. The time between the time zero of the RF pulse and the readout signal is the echo time (TE), and the time between the time zero of the first RF pulse and the time zero of the second RF pulse is the repetition time (TR).

[0031] 4A , other FE sequences may be used. For example, the FE sequence 15 may include spoiler gradient pulses and / or rewind gradient pulses. Furthermore, as will be understood by those skilled in the art, the FE sequence 15 may be a steady-state free precession (SSFP) sequence, a balanced SSFP sequence, or the like. For example, the pulse sequence may refocus magnetization by rewinding the moment of the gradient magnetic field.

[0032] In step 20 of the method 50, a time shift to implement the CAIPI technique is applied to the FE sequence 25. For example, in the SE direction (k z A time shift is applied to one of the RF pulses to induce a moment shift in the direction of the sine wave (G). For example, FIG. 4B shows the second RF pulse shifted back in time by an amount Δt. Also, FIG. 4C shows the second G SS is shifted forward in time by an amount Δt.

[0033] Method 50 illustrates one non-limiting example in which FE sequence 25 is generated using two steps: first, generating an initial FE sequence 15 without a CAIPI k-space sample offset, and then obtaining the CAIPI FE sequence 25 by modifying the initial FE sequence 15 with a time shift that results in a CAIPI k-space sample offset. Note that the methods described herein also include non-limiting examples in which the CAIPI FE sequence 25 is generated directly in one step.

[0034] For example, as shown in FIG. 4C, k z Sample shift Δk in k-space in the SE direction z A time shift is selected to derive (e.g., a sample shift at one location in k-space is Δk z =1 / FOV z3A-3I show the various effects of different sampling patterns (including CAIPI) in the SE direction of k-space. y Directional acceleration coefficient R y and k z Directional acceleration coefficient R z 3B shows the sampling of each point in k-space without using a perfect sampling pattern. If the cone shape shown in Figure 3B is imaged using this perfect sampling pattern, the sagittal and axial images shown in Figure 3C are produced, i.e., no aliasing occurs.

[0035] Figure 3D shows the acceleration coefficient R y = 2 and R z = 1, i.e., k y The sampling points are thinned out every other row in the z direction. By zero-padding the thinned sampling points and then performing a Fourier transform, the sagittal and axial images shown in Figure 3E are generated. In this case, there is no aliasing in the z direction, but there is aliasing in the y direction.

[0036] Figure 3F shows the acceleration coefficient R y = 1 and R z = 2, i.e., k z The sampling points are thinned out every other row in the y-direction. By zero-padding the thinned sampling points and then performing a Fourier transform, the sagittal and axial images shown in Figure 3G are generated. In this case, there is no aliasing in the y-direction, but there is aliasing in the z-direction.

[0037] Figure 3H shows sampling using a CAIPI pattern with a total acceleration coefficient of R = 4. Zero-padding the thinned sampling points before performing a Fourier transform produces the sagittal and axial images shown in Figure 3I. In this case, the aliased image is shifted in both the y and z directions, creating a large gap between the original and aliased images. This is advantageous for enabling separation of these images. Furthermore, for the aliased image, the coil sensitivity is shifted by a larger amount than in either of the cases shown in Figures 3E and 3G. Therefore, using this CAIPI pattern with SENSE processing provides better separation than performing SENSE processing with either of the patterns shown in Figures 3D and 3F.

[0038] The CAIPI sampling pattern may start with the pattern of FIG. 3F, and then proceed to k, as shown in FIG. 3J. y For every other row in the direction, k z This can be achieved by inducing a phase addition of Δφ in the direction of the FE sequence shown in Figure 4A. For example, if the FE sequence shown in Figure 4A provides the sampling pattern of Figure 3F, a time shift of Δt can be added to each of the uniform RF excitation pulses to achieve the CAIPI sampling pattern.

[0039] The CAIPI sampling pattern shown in FIG. 3H is one of many CAIPI sampling patterns. As will be appreciated by those of ordinary skill in the art, other CAIPI sampling patterns may be used, for example, for R=3, 4, 5, ..., 16, .... For example, k z The acceleration coefficient in the direction is R z For R = 2, the initial FE sequence shown in Figure 3F is generated. Then, as shown in Figure 3J, k z The initial value of k y For each sampling point along the direction, SSA CAIPI sampling pattern is generated from the initial FE sequence by changing the relative timing of the RF pulses with respect to the waveform of R. These relative time shifts are calculated to apply a phase offset of Δφ = θ + mod(m(n-1)2π / R,2π) to each sampling point of the non-CAIPI sampling pattern, where mod(x,y) is the modulus operation of x with respect to y, θ is an arbitrary phase where the same phase is applied to all sampling points, m is an integer chosen from the set {1,…,R-1}, and n is the number of k samples. y is the index of the sampling point in the direction. For example, Figure 3J shows an example where R = 2, m = 1, and θ = 0, and k y In every other row in the direction, the sampling points where n is odd are k z In this example, the n-th sampling point is shifted by Δφ=π in the direction of the arrow, and the n-th sampling point is shifted by Δφ=0.

[0040] These other CAIPI sampling patterns consist of RF excitation pulses and G SS The position of the sampling point is determined by the relative timing with the waveform of k. z These shifts are generated by k z The excitation RF pulse and the G SS The time shift with respect to the waveform of can be in either direction (i.e., positive or negative). SS The flat-top portion of the waveform can be extended in either the positive or negative direction by a sufficient distance to accommodate either positive or negative time shifts. SS If the flat-top portion of the waveform extends further in the positive time direction than in the negative time direction relative to time zero, then a positive time shift is more desirable than a negative time shift.

[0041] A related method to the method described herein is to modulate the phase of the RF pulse or modulate the gradient magnetic field,z Unlike related methods, the method described herein uses RF pulses and G SS The relative time shift of the waveform of k z Realize directional displacement.

[0042] In step 30 of method 50 , MRI data 35 is acquired using the FE sequence 25 generated in step 10 and then time-shifted in step 20 .

[0043] In step 40 of method 50, an MRI image 45 is reconstructed from MRI data 35 using a CAIPI sampling pattern, using techniques such as the SENSE and GRAPPA methods. Even when a non-CAIPI sampling pattern, such as that shown in FIG. 3F, is used, there is still some displacement between the original image and the aliased image, as shown in FIG. 3G. This allows parallel imaging techniques (e.g., the SENSE and GRAPPA methods) to perform image de-aliasing. On the other hand, when MRI data 35 is acquired using a CAIPI sampling pattern, the displacement between the original image and the aliased image is increased, as shown in FIG. 3I, thereby improving image de-aliasing performed using parallel imaging techniques. That is, k to realize the CAIPI sampling pattern z The shift in direction minimizes aliasing overlap, thereby reducing the g-factor penalty of exploiting in-plane variations in coil sensitivity. In the SENSE method, receive coil (Rx) sensitivity spatial profiles (called Rx maps) are used to resolve individual slice images.

[0044] Figure 5 shows an addition to the technique described above in connection with CAIPI imaging, in which a shifted aliased image 256 is produced by imaging a volume with a reduced FOV and CAIPI sampling pattern.

[0045] Consider an image 260 acquired with in-plane accelerated imaging at double acceleration. With an acceleration factor of 2, lines of k-space are acquired with a step size twice that of a full FOV acquisition, resulting in a half-FOV image 260 from each coil. This half-FOV image 260 contains aliasing, including the outer corners of the circular center that are positioned around or folded back on the circular center. SENSE reconstruction separates the image from the aliasing in the image by using the coil sensitivities of the multiple coils and the aliased coil sensitivities as two separate coils. This results in an alias-free image 258.

[0046] In some embodiments, the processing logic of the SENSE reconstruction module 252 may be the same as the SENSE reconstruction logic of the reconstruction module 202. The coil sensitivity maps 254 provided as input to the SENSE reconstruction module 252 include at least a multi-coil receive sensitivity map.

[0047] In one embodiment, the magnitude of the time shift is approximately an inverse function of the bandwidth of the RF pulse, G SS The size of the RF pulse depends on the bandwidth and slice dimension (FOV Z ) and the ratio of the FOV thickness to the FOV Z is approximately equal to the slice thickness multiplied by the number of slices, and the spatial moment shift is Z and the time shift is G SS is the ratio of the spatial moment shift to the magnitude of

[0048] In one embodiment, as will be appreciated by those of ordinary skill in the art, an FE sequence may include two phase encoding gradients at spatial shifts in the spatial dimensions corresponding to two phase encoding directions (y and z), and k-space (k y and k z ) may be provided by a pattern of moment shifts in each of both phase encoding directions.

[0049] To generate a given k-space moment for each individual repetition of data sampling, G SS A time shift is used between the waveform and the RF pulse. The resulting k-space moment shift is calculated by the time shift and G SS This time shift is proportional to the product of the magnitude of G while keeping the duration of the RF pulse constant. SS By shifting the time of the waveform of SS A third option is to use an RF pulse and a G SS The waveforms may each be partially shifted so that their relative time difference provides the desired time shift.

[0050] Regarding whether the time shift of the RF pulse adversely affects image quality and chemical species separation, for example, in an FE sequence such as VIBE (Volumetric Interpolated Breath-hold Examination), the time shift is approximately 70 microseconds, which is 3% of the TE corresponding to 10° of water-fat separation. Therefore, the effect of the time shift is small, and the effect on image quality is negligible.

[0051] In a preferred embodiment, G SS The waveforms of the RF pulses are time-shifted, while the timing of the RF pulses is kept constant. This embodiment can avoid the shift of TE as described above. This embodiment can also avoid the TR from repeatedly changing, thereby keeping the timing of the sequence constant.

[0052] Furthermore, as described in U.S. Patent Application No. 15 / 791,898, which is incorporated by reference in its entirety, the readout gradient and phase encoding gradient can be fixed, and G SS The waveform of G can maintain the same amplitude and waveform. SSThe magnetic field waveform can be time shifted, and the extra SS moment in each echo can cause a slice-position dependent shift in the SE direction.

[0053] In one embodiment, G SS The "flat top" period of the waveform can be extended beyond the minimum allowable period, allowing for G SS The flat top is the G SS is the region of uniform amplitude in the center of the waveform.

[0054] In some embodiments, images can be reconstructed using parallel imaging techniques (e.g., SENSE, GRAPPA, etc.), as described in Deshmane, A. et al., "Parallel MR Imaging," J Magn Reson Imaging, vol. 36, pp. 55-72 (2012) and Breuer, F. A. et al., "High Acceleration Achieved by Managing Aliasing in Parallel Imaging for Multislice Imaging (CAIPIRINHA)," Magn Reson Med, vol. 53, pp. 684-691 (2005), both of which are incorporated by reference in their entireties.

[0055] When the acceleration factor is high, parallel imaging methods can suffer from a loss of SNR due to coil geometry constraints. This reduction in SNR due to coil geometry is called the geometry factor penalty, or more commonly, the g-factor penalty. The g-factor can be thought of as a noise amplification factor. For a given receiver coil geometry, regions of the image reconstructed with high g-factor map values are expected to have high noise.

[0056] To better understand how CAIPI techniques improve the g-factor, it is helpful to consider that aliases can overlap with each other. If there is not enough variation in the sensitivity of the receive coil along the undersampled dimension, it can be difficult to separate signals from different aliases. Therefore, the SNR of the reconstructed slice images will be lower due to the increased g-factor penalty. The SNR can be improved by shifting the aliases from each other in the PE and SE directions during acquisition using methods such as CAIPI techniques. CAIPI techniques take full advantage of the sensitivity variations in conventional receive coils. z The directional shift causes voxels in one slice to alias with voxels in one or more other slices that have sufficiently orthogonal receive sensitivity values, thereby reducing the relatively high g-factor penalty associated with separating closely spaced or overlapping aliases, as described in U.S. Patent Application No. 7,002,344, incorporated herein by reference in its entirety.

[0057] While certain embodiments have been described, these embodiments are presented by way of example only and are not intended to limit the technology of the present disclosure. Indeed, the novel methods, apparatus, and systems described herein may be embodied in a variety of other forms. Furthermore, various omissions, substitutions, and changes may be made in the form of the methods, apparatus, and systems described herein without departing from the spirit of the present disclosure.

Claims

1. a pulse sequence for collecting MRI (Magnetic Resonance Imaging) data by applying a radio frequency pulse and a slice selective gradient magnetic field for the same period for each repetition time, wherein a pulse sequence is generated in which a time shift is provided between the center of the radio frequency pulse and the center of the waveform of the slice selective gradient magnetic field for every other repetition time so that a k-space sampling pattern is generated in which the positions of sampling points are shifted in the slice encoding direction for every other row in the phase encoding direction; an acquisition unit that acquires the MRI data using the pulse sequence; a reconstruction unit that reconstructs an MRI image from the MRI data; A magnetic resonance imaging apparatus comprising:

2. the generation unit generates the pulse sequence so that k-space sampling is performed in a sampling pattern similar to a CAIPI (Controlled Aliasing In Parallel Imaging) sampling pattern.

2. The magnetic resonance imaging apparatus according to claim 1.

3. the generator generates the pulse sequence so that the radio frequency pulse excites transverse magnetization of spins in an imaging field of view in a slice encoding direction; the reconstruction unit reconstructs the MRI image using a spatial offset from an image in which aliasing occurs so that aliasing of the MRI image is removed; 3. A magnetic resonance imaging apparatus according to claim 1.

4. the reconstruction unit reconstructs the MRI image using a SENSE (Sensitivity Encoding) method or a GRAPPA (Generalized Autocalibrating Partial Parallel Acquisition) method; 4. The magnetic resonance imaging apparatus according to claim 1.

5. the generating unit generates the pulse sequence so as to shift a center time of the radio frequency pulse and maintain a constant center time of the waveform of the slice selection gradient magnetic field in each period of the repetition time.

5. The magnetic resonance imaging apparatus according to claim 1.

6. the generator generates the pulse sequence so as to further include a waveform of a readout gradient magnetic field, a first waveform of a first phase encoding gradient magnetic field corresponding to a phase encoding direction, and a second waveform of a second phase encoding gradient magnetic field corresponding to another phase encoding direction.

6. The magnetic resonance imaging apparatus according to claim 1.

7. the generating unit generates the pulse sequence so that the waveform of the readout gradient magnetic field is shifted by the same amount as the center of the radio frequency pulse during each period of the repetition time.

7. The magnetic resonance imaging apparatus according to claim 6.

8. the generating unit generates the pulse sequence so as to shift a center time of the waveform of the slice selection gradient magnetic field and maintain a center time of the radio frequency pulse constant during each period of the repetition time.

5. The magnetic resonance imaging apparatus according to claim 1.

9. the generating unit generates the pulse sequence such that, in each period of the repetition time, the amount of the time lag is based on a ratio between a phase offset in each period and a magnitude of the slice selection gradient magnetic field.

9. The magnetic resonance imaging apparatus according to claim 1.

10. the generation unit acquires an initial pulse sequence having a non-CAIPI sampling pattern in which an acceleration coefficient in the slice encoding direction is R, and generates the pulse sequence by applying a spatial moment offset, expressed as Δkz=m / FOVz, where m is an integer selected from the set of {0, 1, 2, ..., R-1}, to each sampling point of the non-CAIPI sampling pattern, so that the acceleration coefficient in the phase encoding direction becomes R.

8. A magnetic resonance imaging apparatus according to claim 6 or 7.

11. the generation unit generates the pulse sequence such that, when an acceleration coefficient in the phase encoding direction is R and m is an integer selected from the set of {0, 1, 2, ..., R-1}, a phase offset is zero and the time lag is zero for an R+m-th sampling point in the other phase encoding direction.

8. A magnetic resonance imaging apparatus according to claim 6 or 7.

12. a radio frequency transmitter that generates a first radio frequency pulse and a second radio frequency pulse; a radio frequency coil for receiving FE to generate the MRI data; a gradient magnetic field coil for generating the slice-selective gradient magnetic field; 12. The magnetic resonance imaging apparatus according to claim 1, further comprising:

13. a step of generating a pulse sequence for collecting MRI (Magnetic Resonance Imaging) data by applying a radio frequency pulse and a slice selective gradient magnetic field for the same period for each repetition time, wherein a time shift is provided between the center of the radio frequency pulse and the center of the waveform of the slice selective gradient magnetic field at every other repetition time so that a k-space sampling pattern is generated in which the positions of sampling points are shifted in the slice encoding direction for every other row in the phase encoding direction; acquiring the MRI data using the pulse sequence; reconstructing an MRI image from the MRI data; A magnetic resonance imaging method comprising:

Images