US10775467B2 - Magnetic resonance method and apparatus for echo-planar imaging with a zigzag-type trajectory in raw-data space - Google Patents

Magnetic resonance method and apparatus for echo-planar imaging with a zigzag-type trajectory in raw-data space Download PDF

Info

Publication number
US10775467B2
US10775467B2 US15/830,329 US201715830329A US10775467B2 US 10775467 B2 US10775467 B2 US 10775467B2 US 201715830329 A US201715830329 A US 201715830329A US 10775467 B2 US10775467 B2 US 10775467B2
Authority
US
United States
Prior art keywords
raw
data space
data
space
zigzag
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active, expires
Application number
US15/830,329
Other versions
US20180156884A1 (en
Inventor
Robin Heidemann
Patrick Liebig
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthineers AG
Original Assignee
Siemens Healthcare GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Healthcare GmbH filed Critical Siemens Healthcare GmbH
Assigned to SIEMENS HEALTHCARE GMBH reassignment SIEMENS HEALTHCARE GMBH ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: HEIDEMANN, ROBIN, LIEBIG, PATRICK, MR.
Publication of US20180156884A1 publication Critical patent/US20180156884A1/en
Application granted granted Critical
Publication of US10775467B2 publication Critical patent/US10775467B2/en
Assigned to Siemens Healthineers Ag reassignment Siemens Healthineers Ag ASSIGNMENT OF ASSIGNOR'S INTEREST Assignors: SIEMENS HEALTHCARE GMBH
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/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/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • 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
    • 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/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
    • 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
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • G06T12/10—
    • G06T12/20—
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4084—Scaling of whole images or parts thereof, e.g. expanding or contracting in the transform domain, e.g. fast Fourier transform [FFT] domain scaling
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/52—Parallel processing
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography

Definitions

  • the present invention concerns a method for echo-planar MR imaging with which MR signals are acquired with a zigzag-type trajectory in the raw-data space (k-space).
  • the invention further relates to a magnetic resonance apparatus and an electrically readable data storage medium for implementing such a method.
  • EPI echo-planar imaging
  • MRT magnetic resonance tomography
  • One EPI recording technique uses the so-called zigzag trajectory with which continuous phase encoding is used with simultaneous sinusoidal readout gradients. This technique has the advantage of producing less noise when recording MR signals, and achieves a shorter echo time as a result of the more effective, continuous data recording.
  • the data points are not arranged on a Cartesian grid in k-space, which is a basic prerequisite for reconstructing image data from the k-space data using a reconstruction algorithm of the type employed in Cartesian parallel acquisition imaging, such as GRAPPA (generalized autocalibrating partially parallel acquisition).
  • GRAPPA generalized autocalibrating partially parallel acquisition
  • the data are already undersampled fourfold.
  • the higher reconstruction factor in the GRAPPA reconstruction results in an increase in the g factor, which, together with the acceleration factor for the acquisition, describes the signal-to-noise loss caused by the parallel imaging reconstruction. If additional data are reconstructed from undersampled data by the use of a parallel imaging reconstruction algorithm, there is now a region with twofold undersampling in the middle of the k-space lines. This means that it is not possible to combine the GRAPPA reconstruction with an acceleration factor and to use iFT.
  • An object of the present invention is to use an accelerated recording with the use of parallel imaging techniques while simultaneously not increasing the noise in the recorded MR image to an excessive degree.
  • MR signals are entered in raw-data space along a zigzag-type trajectory.
  • MR raw-data points are recorded with a sequence of readout gradients and phase-encoding gradients that cause a zigzag-type undersampled trajectory in the acquired raw-data space.
  • the acquisition is performed such that the signal echoes are shifted by up to a quarter of the acquired raw-data space in the direction of the readout axis. Furthermore, an additional zigzag-type trajectory in the acquired raw-data space is reconstructed by a parallel-imaging method, and using the acquired zigzag-type trajectory. Then, at least one MR image is reconstructed using interlaced Fourier transform (iFT).
  • iFT interlaced Fourier transform
  • Shifting the signal echoes by up to a quarter of the acquired raw-data space in the direction of the readout axis causes the signal echoes to occur in a region of the raw-data space in which the spacing of the trajectories in the phase-encoding direction correspond to one ⁇ ky.
  • ⁇ ky is proportional to the integral over the gradient along the phase-encoding axis, multiplied by the time in which the gradient is present.
  • the difference between two trajectories in the direction of the phase-encoding axis is now less than one ⁇ ky, which means that the Nyquist theorem is satisfied here. If the echo is shifted by less than a quarter, the raw-data center is present in a region in which the spacing between two points in the phase-encoding direction is alternately smaller and greater than the minimum spacing required by the Nyquist theorem for reconstruction without aliasing. This means that the average spacing between adjacent points in the phase-encoding direction corresponds to the Nyquist criterion and the data can be reconstructed successfully.
  • the raw data in the raw-data space are extended such that an extension segment is formed, this extension taking place with the application of a point symmetry in the raw-data space with respect to the raw-data center, so as to create an extended raw-data space encompassing the acquired raw-data space and the extension segment.
  • the raw data in the extension segment are calculated using point reflection at the raw-data center.
  • the raw data in the extension segment can be determined using the so-called zero-filling technique or with the partial Fourier technique.
  • the signal echoes can be shifted in the direction of a readout axis by shortening or lengthening the dephasing pre-gradient before the signal readout.
  • the signal echoes are shifted in a first direction of the readout axis, for example by a quarter.
  • the extension segment then also lies in the first direction of the readout axis relative to the acquired raw-data space.
  • One iFT can, for example, be applied to the acquired raw-data space to create a hybrid space and one iFT can be applied to the extension segment to create a hybrid extension segment, wherein the hybrid space and the hybrid extension segment are merged to represent a total hybrid space.
  • This total hybrid space can be subjected to a Fourier transform along this axis to create the at least one MR image.
  • the use of the extension segment provides two possibilities.
  • the recording time is reduced while the resolution remains constant, with the other possibility, the recording time is left the same and the resolution is increased in the direction of the readout axis.
  • this is achieved in that the at least one MR image that is reconstructed has a desired resolution, wherein the recording time required for the desired resolution of the at least one MR-image is reduced in that the raw-data space required for the desired resolution is not completely recorded.
  • the raw-data space required for the desired resolution corresponds to the extended raw-data space with the extension segment. Since only a partial area was acquired in the extended raw-data space, the recording time is shorter than with a recording with which the entire extended raw-data space with MR signals was acquired for the desired resolution.
  • the MR image is recorded with a defined recording time, wherein the raw data recorded in the defined recording time without using the extension segment can be used to reconstruct an MR reference image with a defined resolution.
  • the MR image reconstructed with the iFT using the extension segment then has higher resolution than the defined resolution, namely higher resolution in the direction of the readout axis.
  • the reconstruction of the additional zigzag-type trajectory is preferably performed using the GRAPPA method.
  • the additional zigzag-type trajectory can be reconstructed such that an incline of the last segment of the acquired zigzag-type trajectory in the direction of the positive phase-encoding direction has an opposite preliminary sign than an incline of the last segment of the reconstructed zigzag-type trajectory in the direction of the positive phase-encoding direction.
  • the invention encompasses a magnetic resonance apparatus with an image-sequence controller that operates an MR data acquisition scanner of the apparatus so as to produce the zigzag-type trajectory, and a computer for the reconstruction of the at least one MR image as described in detail above and below.
  • the present invention also encompasses a non-transitory, computer-readable data storage medium encoded with programming instructions that, when the storage medium is loaded into a computer of a magnetic resonance apparatus, cause the computer to operate the magnetic resonance apparatus according to any or all of the embodiments of the method according to the invention, as described above.
  • FIG. 1 schematically illustrates an MR apparatus with which MR images can be calculated according to the invention.
  • FIG. 2 schematically illustrates raw-data space and the relationship between zigzag-type trajectories and the spacing between those trajectories.
  • FIG. 3 schematically illustrates an example of raw-data space with a recorded and reconstructed trajectory, wherein image data cannot be reconstructed from the raw-data space shown in FIG. 3 without artifacts.
  • FIG. 4 shows an example of raw-data space, which is filled with zigzag-type trajectories by extension in the readout direction, such that artifact-free reconstruction of an MR image is possible with undersampling.
  • FIG. 5 shows an example of raw-data space, which was completely recorded, but with which artifact-free reconstruction is not possible with undersampling.
  • FIG. 6 shows the recording of raw-data space with which, with a given recording time, the resolution can be increased by extension in the readout direction.
  • FIG. 7 shows the recording of raw-data space with which the recording time can be reduced while the resolution remains constant.
  • FIG. 8 is a flowchart for artifact-free reconstruction of MR images by extension of the raw-data space in the readout direction in accordance with the invention.
  • an MR apparatus 10 is explained with which, as explained below, MR images are recorded with echo-planar imaging and undersampling, wherein, despite the undersampling, artifact-free images are generated.
  • the MR apparatus 10 has a scanner 11 that has a basic field magnet that generates a polarization field B 0 .
  • An examination subject 13 on a bed 12 is moved into the scanner 11 in order to record spatially encoded magnetic resonance signals from the examination subject 13 .
  • the coils used for signal recording such as a whole-body coil or local coils are not depicted, for clarity.
  • the invention can be applied with so-called parallel imaging with which the MR signals are recorded simultaneously with multiple local coils, i.e. a coil array or local coils.
  • the radiation of radio-frequency pulses and switching of magnetic field gradients can cause the magnetization generated by the polarization field B 0 to be deflected out of the equilibrium position and spatially encoded and the resulting magnetization is detected by the receiver coils.
  • Methods for generating MR images by radiating the RF pulses and switching magnetic field gradients in different combinations and sequences of MR images are in principle known to those skilled in the art, and so need not be explained in further detail herein.
  • the MR apparatus furthermore has a control computer 20 that controls the MR apparatus 10 .
  • the control computer 20 has a gradient controller 15 for controlling and switching the necessary magnetic field gradients.
  • An RF controller 14 is provided to control and generate the RF pulses for deflecting the magnetization.
  • An image-sequence controller 16 controls the sequence of magnetic field gradients and RF pulses and hence indirectly the gradient controller 15 and the RF controller 14 .
  • An operator can control the MR apparatus 10 via an input unit 17 and MR images and other information required for the control can be displayed on a display monitor 18 .
  • a computing unit 19 with at least one processor (not shown) is provided to control the different units in the control computer 20 .
  • a memory 21 is provided in which program modules or programs can be stored, which, when they are executed by the computing unit 19 or its processor, control the operation of the MR apparatus 10 .
  • the computing unit 19 is designed to calculate, according to the invention, MR images from the acquired MR signals such that undersampling is possible with a zigzag trajectory, wherein simultaneously a substantially artifact-free reconstruction with a satisfactory signal-noise ratio of the MR images is achieved.
  • a combination of interlaced Fourier transform (iFT) and parallel-imaging requires the accelerated data to be recorded such that, following a parallel image reconstruction, the condition for iFT are satisfied.
  • This condition is that the spacing of the raw data around the center of raw-data space, i.e. k-space, causes the Nyquist theorem to be satisfied, or in other words is smaller or equal to ⁇ ky is.
  • the first possibility is to record only 2 ⁇ 3 in the readout direction of raw-data space or k-space. This shortens the recording time and, following a partial Fourier reconstruction in the readout direction of a readout axis, produces the same resolution as that originally desired for the k-space.
  • the k-space center is then shifted so as to be located in a region in which the Nyquist criterion is satisfied.
  • Another possibility for achieving this condition is to keep the recording time the same and to shift the recorded k-space by changing the pre-gradient in the direction of the readout axis. Since in this embodiment the position of the k-space center does not change, k-space is ultimately sampled asymmetrically. This again has the result that the k-space center is located on a region in which the Nyquist criterion is satisfied.
  • a raw-data space 30 in the kx and ky directions, wherein kx is the readout axis and ky the phase-encoding axis, wherein a recorded trajectory 31 fills this raw-data space.
  • the maximum spacing in the center of the raw-data space is 1 ⁇ ky, while at the edge it is 2 ⁇ ky.
  • This raw dataset can be reconstructed with the iFT such that the signal-noise ratio is satisfactory and no substantial artifacts occur; this will be described as artifact-free below.
  • a dashed line indicates a raw data-zigzag trajectory 32 with which the recording was made with an acceleration factor of 2, i.e. simple undersampling.
  • the individual segments of the trajectory 32 have a spacing of 2 ⁇ ky in the middle and a maximum spacing of 4 ⁇ ky at the edge. This dataset 32 cannot be reconstructed to form an image without artifacts with the iFT.
  • FIG. 3 depicts raw-data space 40 with a recorded zigzag trajectory 41 .
  • the zigzag-type trajectory takes place via a periodically switched readout gradient and via a phase-encoding gradient corresponding to the absolute value of the readout gradient with a correspondingly lower amplitude.
  • the trajectory was recorded with an acceleration factor of 2, i.e. only 50% of the raw data was recorded, while the trajectory 42 was reconstructed with parallel imaging techniques, for example GRAPPA.
  • This raw dataset depicted in FIG. 3 cannot be reconstructed to form an image without artifacts with the iFT since the maximum spacing is 2 ⁇ ky in the center and 2 ⁇ ky at the edge.
  • a further problem is that the points of intersection in the central raw-data region lead to artifacts with the iFT. Artifact-free reconstruction would require a further GRAPPA reconstruction.
  • FIG. 4 depicts raw-data space with which the reconstruction takes place such that artifact-free reconstruction of MR images is possible even with undersampling with a factor of 2.
  • the raw-data space has an acquired raw-data space 51 in which a zigzag-type trajectory 52 was recorded as depicted. Reducing the dephasing pre-gradient relative to the gradient during the readout time causes the signal echo to occur earlier and not in the geometric center of the acquired raw-data space, but in the raw-data center 60 .
  • an additional zigzag-type trajectory 53 is calculated, for example using GRAPPA reconstruction.
  • raw-data space is now read out asymmetrically and it is possible for an extended raw-data space to be formed, two thirds of which are the acquired raw-data space 51 and one third is an extension segment 55 .
  • the extended raw-data space 50 now includes the acquired raw-data space 51 and the extension segment 55 .
  • the raw-data points in the extension segment 55 can be determined by so-called zero filling or with partial Fourier. This uses the point symmetry of the extended raw-data space.
  • the raw-data center 60 is used for this.
  • the raw-data points in the extension segment 55 are calculated by point reflection at the raw-data center 60 . If the extended raw-data space 50 is now considered, it is identifiable that the spacing is only 1 ⁇ ky in the raw-data center while it is maximum 2 ⁇ ky at the edge. Furthermore, the points of intersection of the individual zigzag-type trajectories are not in the raw-data center 60 . This enables artifact-free reconstruction with the extended raw-data space 50 with the iFT. Hence, it is possible to record a two-fold accelerated MR image with echo-planar imaging using a GRAPPA reconstruction.
  • a raw-data space 61 depicted in FIG. 5 has the same extension in the direction of the readout axis kx as the acquired raw-data space 51 in FIG. 4 .
  • the maximum spacing of the segments of a trajectory is 2 ⁇ ky so that here once again artifact-free reconstruction is not possible. Artifact-free reconstruction is not possible without shifting the signal echoes.
  • FIG. 6 depicts an embodiment similar to that in FIG. 4 , wherein once again raw-data space is divided into a recorded raw-data space 71 and an extension segment 75 in the readout direction.
  • the length of the acquired raw-data space 71 with the acquired trajectory 73 and the reconstructed trajectory 72 once again corresponds to the length of the raw-data space 61 in FIG. 5 .
  • the asymmetrical signal readout and the extension with point reflection about the center 76 in the kx direction with filling in the extension segment 75 again causes the raw-data space to be changed such that there is a maximum spacing of 1 ⁇ ky in the raw-data center, while the maximum spacing at the edge is 2 ⁇ ky.
  • the extended raw-data space 70 can be reconstructed artifact-free.
  • the recording time is the same as with the recording of the raw-data space in FIG. 5 , but the resolution is higher in the kx direction than it is in an MR image calculated with the raw-data space in FIG. 5 .
  • the present invention enables the recording of raw-data points with undersampling using continuous phase encoding.
  • the raw-data points do not lie on a Cartesian grid, but this is a perquisite for the use of Cartesian parallel imaging such as GRAPPA.
  • the k-space center is no longer located in the middle of the acquired raw-data space as a result of which, first, the maximum spacing of the raw-data points is reduced and, furthermore, the points of intersection of the individual segments of the trajectory are not located in the raw-data center. It is now possible, for example, only to perform an iFT in the segment with the k-space center in order to obtain an image.
  • extension segment primarily containing radio-frequency k-space information.
  • This image could be combined with the first image in order to increase the resolution.
  • a further possibility consists in not increasing the resolution, but recording less data in the kx direction. In this case, the recording time would be shorter than with a completely sampled raw-data space.
  • the missing segment, the extension segment can be replaced either by zero filling or by a partial Fourier reconstruction from the data of the two recorded segments, since the data is point-symmetric with respect to the k-space center. This embodiment is described in FIG. 7 .
  • the entire raw-data space 90 depicted in FIG. 7 has the same extension in the kx direction as the acquired raw-data space 71 .
  • the method is summarized in FIG. 8 .
  • the method starts in step S 80 and in step S 81 zigzag-type segments of a trajectory are recorded, wherein following signal recording, the raw dataset is divided into even and odd echoes.
  • the signal recording is performed such that the maximum signal intensity, the center of the signal echoes, is shifted, preferably by a quarter of the recorded raw-data space, in the direction of the readout axis.
  • step S 82 an additional trajectory is recorded for the even and odd recorded raw-data points in each case using a parallel recording technique.
  • step S 83 the raw-data space is extended in the readout direction, as was described in connection with FIGS. 4, 6 and 7 .
  • step S 84 an iFT is performed each for each segment of the extended raw-data space along the phase-encoding direction.
  • step S 84 and the iFT there is now a part-transformed raw-data space and a part-transformed extension segment.
  • Part-transformed means that a Fourier transform was performed only along one direction, namely along the phase-encoding direction.
  • step S 85 the part-transformed raw-data space and the part-transformed extension segment are now used to apply an inverse Fourier transform back into the raw-data space (inverse Fourier transform in the phase-encoding direction). After step S 85 , this results in a hybrid space and a hybrid extension segment.
  • step S 86 the elements generated in step S 85 , the hybrid space and the hybrid extension segment are merged to form a total hybrid space.
  • step S 87 the Fourier transform can, as known, be performed to create the at least one MR image. The method ends in step S 88 .
  • the above-described method makes it possible to apply a method such as is known from parallel imaging with EPI data recorded on a zigzag trajectory.
  • a method such as is known from parallel imaging with EPI data recorded on a zigzag trajectory.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)

Abstract

In a method and magnetic resonance (MR) apparatus for echo-planar MR imaging with which MR signals are entered in raw-data space with a zigzag-type trajectory, a sequence of readout gradients and phase-encoding gradients is applied such that a zigzag-type undersampled trajectory in raw-data space is filled, such that with the signal echoes being shifted by up to a quarter of the acquired raw-data space in the direction of the readout axis. Image data are reconstructed using a parallel-imaging reconstruction algorithm that operates on the raw data acquired in the zigzag-type trajectory, based on an interlaced Fourier transform.

Description

BACKGROUND OF THE INVENTION Field of the Invention
The present invention concerns a method for echo-planar MR imaging with which MR signals are acquired with a zigzag-type trajectory in the raw-data space (k-space). The invention further relates to a magnetic resonance apparatus and an electrically readable data storage medium for implementing such a method.
Description of the Prior Art
High-speed recording methods using the EPI (echo-planar imaging) principle are known in magnetic resonance tomography (MRT). One EPI recording technique uses the so-called zigzag trajectory with which continuous phase encoding is used with simultaneous sinusoidal readout gradients. This technique has the advantage of producing less noise when recording MR signals, and achieves a shorter echo time as a result of the more effective, continuous data recording. However, with such continuous phase encoding, the data points are not arranged on a Cartesian grid in k-space, which is a basic prerequisite for reconstructing image data from the k-space data using a reconstruction algorithm of the type employed in Cartesian parallel acquisition imaging, such as GRAPPA (generalized autocalibrating partially parallel acquisition). Known from DE 10 2013 100 349 A1 is the use of a phase encoding that makes the data lie on a mathematically correct zigzag trajectory. It is also possible to select other forms; they only need to be an exactly scaled amount of the readout gradient.
It is possible to reconstruct image data from the, on-average, non-undersampled MR data by the use of an interlaced Fourier transform (iFT). With undersampled, i.e. accelerated, data acquisition, it is no longer possible to reconstruct the data with iFT. In such a case, it is necessary to use a GRAPPA reconstruction taking account of the fact that, with an assumed acceleration factor of 2 in the acquisition, the GRAPPA reconstruction has to be used with a reconstruction factor of 4. This is due to the fact that the GRAPPA reconstruction is applied separately in the even or odd k-space lines in the raw-data space. These separated datasets are now two-fold undersampled for a non-accelerated recording. For a two-fold accelerated recording, the data are already undersampled fourfold. The higher reconstruction factor in the GRAPPA reconstruction results in an increase in the g factor, which, together with the acceleration factor for the acquisition, describes the signal-to-noise loss caused by the parallel imaging reconstruction. If additional data are reconstructed from undersampled data by the use of a parallel imaging reconstruction algorithm, there is now a region with twofold undersampling in the middle of the k-space lines. This means that it is not possible to combine the GRAPPA reconstruction with an acceleration factor and to use iFT.
SUMMARY OF THE INVENTION
An object of the present invention is to use an accelerated recording with the use of parallel imaging techniques while simultaneously not increasing the noise in the recorded MR image to an excessive degree.
In the method for echo-planar MR imaging according to the invention, MR signals are entered in raw-data space along a zigzag-type trajectory. In accordance with the invention, MR raw-data points are recorded with a sequence of readout gradients and phase-encoding gradients that cause a zigzag-type undersampled trajectory in the acquired raw-data space.
The acquisition is performed such that the signal echoes are shifted by up to a quarter of the acquired raw-data space in the direction of the readout axis. Furthermore, an additional zigzag-type trajectory in the acquired raw-data space is reconstructed by a parallel-imaging method, and using the acquired zigzag-type trajectory. Then, at least one MR image is reconstructed using interlaced Fourier transform (iFT).
Shifting the signal echoes by up to a quarter of the acquired raw-data space in the direction of the readout axis causes the signal echoes to occur in a region of the raw-data space in which the spacing of the trajectories in the phase-encoding direction correspond to one Δky. This means that a new raw-data center is present at a location where the sampling has taken place densely enough for the performance of a reconstruction using iFT. As is known, Δky is proportional to the integral over the gradient along the phase-encoding axis, multiplied by the time in which the gradient is present. The difference between two trajectories in the direction of the phase-encoding axis is now less than one Δky, which means that the Nyquist theorem is satisfied here. If the echo is shifted by less than a quarter, the raw-data center is present in a region in which the spacing between two points in the phase-encoding direction is alternately smaller and greater than the minimum spacing required by the Nyquist theorem for reconstruction without aliasing. This means that the average spacing between adjacent points in the phase-encoding direction corresponds to the Nyquist criterion and the data can be reconstructed successfully.
Preferably, before the reconstruction of the at least one MR image, the raw data in the raw-data space are extended such that an extension segment is formed, this extension taking place with the application of a point symmetry in the raw-data space with respect to the raw-data center, so as to create an extended raw-data space encompassing the acquired raw-data space and the extension segment. The raw data in the extension segment are calculated using point reflection at the raw-data center.
The raw data in the extension segment can be determined using the so-called zero-filling technique or with the partial Fourier technique.
The signal echoes can be shifted in the direction of a readout axis by shortening or lengthening the dephasing pre-gradient before the signal readout. In the case of a reduced dephasing pre-gradient, signals occur earlier than in the case where the signal echo occurs in the middle of a readout window. Here, the signal echoes are shifted in a first direction of the readout axis, for example by a quarter. The extension segment then also lies in the first direction of the readout axis relative to the acquired raw-data space.
One iFT can, for example, be applied to the acquired raw-data space to create a hybrid space and one iFT can be applied to the extension segment to create a hybrid extension segment, wherein the hybrid space and the hybrid extension segment are merged to represent a total hybrid space. This total hybrid space can be subjected to a Fourier transform along this axis to create the at least one MR image.
The use of the extension segment provides two possibilities. In a first possibility, the recording time is reduced while the resolution remains constant, with the other possibility, the recording time is left the same and the resolution is increased in the direction of the readout axis. With the first possibility, this is achieved in that the at least one MR image that is reconstructed has a desired resolution, wherein the recording time required for the desired resolution of the at least one MR-image is reduced in that the raw-data space required for the desired resolution is not completely recorded. The raw-data space required for the desired resolution corresponds to the extended raw-data space with the extension segment. Since only a partial area was acquired in the extended raw-data space, the recording time is shorter than with a recording with which the entire extended raw-data space with MR signals was acquired for the desired resolution.
In the second option with the increased resolution, the MR image is recorded with a defined recording time, wherein the raw data recorded in the defined recording time without using the extension segment can be used to reconstruct an MR reference image with a defined resolution. The MR image reconstructed with the iFT using the extension segment then has higher resolution than the defined resolution, namely higher resolution in the direction of the readout axis.
The reconstruction of the additional zigzag-type trajectory is preferably performed using the GRAPPA method. However, it is possible to use other methods. It is generally possible to use Cartesian parallel-imaging methods, for example SENSE (SENSitivity Encoding).
The additional zigzag-type trajectory can be reconstructed such that an incline of the last segment of the acquired zigzag-type trajectory in the direction of the positive phase-encoding direction has an opposite preliminary sign than an incline of the last segment of the reconstructed zigzag-type trajectory in the direction of the positive phase-encoding direction. In other words, this means that the acquired trajectory and the reconstructed trajectory end or start at opposite corners of the raw-data space in the phase-encoding direction.
Furthermore, the invention encompasses a magnetic resonance apparatus with an image-sequence controller that operates an MR data acquisition scanner of the apparatus so as to produce the zigzag-type trajectory, and a computer for the reconstruction of the at least one MR image as described in detail above and below.
The present invention also encompasses a non-transitory, computer-readable data storage medium encoded with programming instructions that, when the storage medium is loaded into a computer of a magnetic resonance apparatus, cause the computer to operate the magnetic resonance apparatus according to any or all of the embodiments of the method according to the invention, as described above.
The above-described features and features described below can be used not only in the corresponding explicitly described combinations, but also in further combinations, unless expressly stated otherwise, or individually, without departing from the scope of protection of the present invention.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 schematically illustrates an MR apparatus with which MR images can be calculated according to the invention.
FIG. 2 schematically illustrates raw-data space and the relationship between zigzag-type trajectories and the spacing between those trajectories.
FIG. 3 schematically illustrates an example of raw-data space with a recorded and reconstructed trajectory, wherein image data cannot be reconstructed from the raw-data space shown in FIG. 3 without artifacts.
FIG. 4 shows an example of raw-data space, which is filled with zigzag-type trajectories by extension in the readout direction, such that artifact-free reconstruction of an MR image is possible with undersampling.
FIG. 5 shows an example of raw-data space, which was completely recorded, but with which artifact-free reconstruction is not possible with undersampling.
FIG. 6 shows the recording of raw-data space with which, with a given recording time, the resolution can be increased by extension in the readout direction.
FIG. 7 shows the recording of raw-data space with which the recording time can be reduced while the resolution remains constant.
FIG. 8 is a flowchart for artifact-free reconstruction of MR images by extension of the raw-data space in the readout direction in accordance with the invention.
 DESCRIPTION OF THE PREFERRED EMBODIMENTS
The following describes the present invention with reference to preferred embodiments and also with reference to the drawings. In the figures, the same reference characters identify the same or similar elements. Furthermore, the figures are schematic representations of different embodiments of the invention. Elements depicted in the figures are not necessarily depicted true-to-scale. Instead, the elements shown in the figures are reproduced such that their function and purpose are understandable to the person skilled in the art. The links between functional units or other elements depicted in the figures can also be implemented as indirect links, wherein a link can be wireless or hardwired. Functional units can be implemented as hardware, software or a combination of hardware and software.
With reference to FIG. 1, an MR apparatus 10 is explained with which, as explained below, MR images are recorded with echo-planar imaging and undersampling, wherein, despite the undersampling, artifact-free images are generated. The MR apparatus 10 has a scanner 11 that has a basic field magnet that generates a polarization field B0. An examination subject 13 on a bed 12 is moved into the scanner 11 in order to record spatially encoded magnetic resonance signals from the examination subject 13. The coils used for signal recording, such as a whole-body coil or local coils are not depicted, for clarity. The invention can be applied with so-called parallel imaging with which the MR signals are recorded simultaneously with multiple local coils, i.e. a coil array or local coils. The radiation of radio-frequency pulses and switching of magnetic field gradients can cause the magnetization generated by the polarization field B0 to be deflected out of the equilibrium position and spatially encoded and the resulting magnetization is detected by the receiver coils. Methods for generating MR images by radiating the RF pulses and switching magnetic field gradients in different combinations and sequences of MR images are in principle known to those skilled in the art, and so need not be explained in further detail herein.
The MR apparatus furthermore has a control computer 20 that controls the MR apparatus 10. The control computer 20 has a gradient controller 15 for controlling and switching the necessary magnetic field gradients. An RF controller 14 is provided to control and generate the RF pulses for deflecting the magnetization. An image-sequence controller 16 controls the sequence of magnetic field gradients and RF pulses and hence indirectly the gradient controller 15 and the RF controller 14. An operator can control the MR apparatus 10 via an input unit 17 and MR images and other information required for the control can be displayed on a display monitor 18. A computing unit 19 with at least one processor (not shown) is provided to control the different units in the control computer 20. Furthermore, a memory 21 is provided in which program modules or programs can be stored, which, when they are executed by the computing unit 19 or its processor, control the operation of the MR apparatus 10. As explained below, the computing unit 19 is designed to calculate, according to the invention, MR images from the acquired MR signals such that undersampling is possible with a zigzag trajectory, wherein simultaneously a substantially artifact-free reconstruction with a satisfactory signal-noise ratio of the MR images is achieved.
As explained in detail below, a combination of interlaced Fourier transform (iFT) and parallel-imaging requires the accelerated data to be recorded such that, following a parallel image reconstruction, the condition for iFT are satisfied. This condition is that the spacing of the raw data around the center of raw-data space, i.e. k-space, causes the Nyquist theorem to be satisfied, or in other words is smaller or equal to Δky is. There are two ways to achieve this. The first possibility is to record only ⅔ in the readout direction of raw-data space or k-space. This shortens the recording time and, following a partial Fourier reconstruction in the readout direction of a readout axis, produces the same resolution as that originally desired for the k-space. The k-space center is then shifted so as to be located in a region in which the Nyquist criterion is satisfied. Another possibility for achieving this condition is to keep the recording time the same and to shift the recorded k-space by changing the pre-gradient in the direction of the readout axis. Since in this embodiment the position of the k-space center does not change, k-space is ultimately sampled asymmetrically. This again has the result that the k-space center is located on a region in which the Nyquist criterion is satisfied.
In connection with FIGS. 2 and 3, an explanation is given of the conditions and trajectories with which artifact-free reconstruction of an MR image using the iFT is possible or not possible. Herein, it is necessary to take account of the fact that, with artifact-free reconstruction in particular account is taken of the loudness of the noise in the generated MR image. In the following discussion, it is assumed that, with a spacing of raw-data segments of the zigzag trajectory, the Nyquist theorem is satisfied when the spacing in the direction of the phase-encoding axis is 1 Δky (one times Δky). 2 Δky (two times Δky) means simple undersampling, etc. FIG. 2 depicts a raw-data space 30 in the kx and ky directions, wherein kx is the readout axis and ky the phase-encoding axis, wherein a recorded trajectory 31 fills this raw-data space. The maximum spacing in the center of the raw-data space is 1 Δky, while at the edge it is 2 Δky. This raw dataset can be reconstructed with the iFT such that the signal-noise ratio is satisfactory and no substantial artifacts occur; this will be described as artifact-free below. Furthermore, a dashed line indicates a raw data-zigzag trajectory 32 with which the recording was made with an acceleration factor of 2, i.e. simple undersampling. As can be identified, the individual segments of the trajectory 32 have a spacing of 2 Δky in the middle and a maximum spacing of 4 Δky at the edge. This dataset 32 cannot be reconstructed to form an image without artifacts with the iFT.
FIG. 3 depicts raw-data space 40 with a recorded zigzag trajectory 41. As known to the person skilled in the art, the zigzag-type trajectory takes place via a periodically switched readout gradient and via a phase-encoding gradient corresponding to the absolute value of the readout gradient with a correspondingly lower amplitude. The trajectory was recorded with an acceleration factor of 2, i.e. only 50% of the raw data was recorded, while the trajectory 42 was reconstructed with parallel imaging techniques, for example GRAPPA. This raw dataset depicted in FIG. 3 cannot be reconstructed to form an image without artifacts with the iFT since the maximum spacing is 2 Δky in the center and 2 Δky at the edge. A further problem is that the points of intersection in the central raw-data region lead to artifacts with the iFT. Artifact-free reconstruction would require a further GRAPPA reconstruction.
FIG. 4 depicts raw-data space with which the reconstruction takes place such that artifact-free reconstruction of MR images is possible even with undersampling with a factor of 2. The raw-data space has an acquired raw-data space 51 in which a zigzag-type trajectory 52 was recorded as depicted. Reducing the dephasing pre-gradient relative to the gradient during the readout time causes the signal echo to occur earlier and not in the geometric center of the acquired raw-data space, but in the raw-data center 60. Furthermore, in the example depicted, similarly to in FIG. 3, an additional zigzag-type trajectory 53 is calculated, for example using GRAPPA reconstruction. With this acquired raw-data space, the maximum spacing is 2 Δky in the center and at the edge so that, as explained in connection with FIGS. 2 and 3, artifact-free reconstruction is not possible. According to the invention, raw-data space is now read out asymmetrically and it is possible for an extended raw-data space to be formed, two thirds of which are the acquired raw-data space 51 and one third is an extension segment 55. The extended raw-data space 50 now includes the acquired raw-data space 51 and the extension segment 55. The raw-data points in the extension segment 55 can be determined by so-called zero filling or with partial Fourier. This uses the point symmetry of the extended raw-data space. The raw-data center 60 is used for this. The raw-data points in the extension segment 55 are calculated by point reflection at the raw-data center 60. If the extended raw-data space 50 is now considered, it is identifiable that the spacing is only 1 Δky in the raw-data center while it is maximum 2 Δky at the edge. Furthermore, the points of intersection of the individual zigzag-type trajectories are not in the raw-data center 60. This enables artifact-free reconstruction with the extended raw-data space 50 with the iFT. Hence, it is possible to record a two-fold accelerated MR image with echo-planar imaging using a GRAPPA reconstruction.
A raw-data space 61 depicted in FIG. 5 has the same extension in the direction of the readout axis kx as the acquired raw-data space 51 in FIG. 4. However, it is identifiable that, both in the raw-data-space center or k-space center and also at the edge, the maximum spacing of the segments of a trajectory is 2 Δky so that here once again artifact-free reconstruction is not possible. Artifact-free reconstruction is not possible without shifting the signal echoes.
FIG. 6 depicts an embodiment similar to that in FIG. 4, wherein once again raw-data space is divided into a recorded raw-data space 71 and an extension segment 75 in the readout direction. The length of the acquired raw-data space 71 with the acquired trajectory 73 and the reconstructed trajectory 72 once again corresponds to the length of the raw-data space 61 in FIG. 5. The asymmetrical signal readout and the extension with point reflection about the center 76 in the kx direction with filling in the extension segment 75 again causes the raw-data space to be changed such that there is a maximum spacing of 1 Δky in the raw-data center, while the maximum spacing at the edge is 2 Δky. Hence, the extended raw-data space 70 can be reconstructed artifact-free. The recording time is the same as with the recording of the raw-data space in FIG. 5, but the resolution is higher in the kx direction than it is in an MR image calculated with the raw-data space in FIG. 5.
The present invention enables the recording of raw-data points with undersampling using continuous phase encoding. With this continuous phase encoding, the raw-data points do not lie on a Cartesian grid, but this is a perquisite for the use of Cartesian parallel imaging such as GRAPPA. As a result of the extension in the kx direction and asymmetric acquisition, the k-space center is no longer located in the middle of the acquired raw-data space as a result of which, first, the maximum spacing of the raw-data points is reduced and, furthermore, the points of intersection of the individual segments of the trajectory are not located in the raw-data center. It is now possible, for example, only to perform an iFT in the segment with the k-space center in order to obtain an image. It is additionally possible to generate an image by using iFT in the extension segment primarily containing radio-frequency k-space information. This image could be combined with the first image in order to increase the resolution. A further possibility consists in not increasing the resolution, but recording less data in the kx direction. In this case, the recording time would be shorter than with a completely sampled raw-data space. The missing segment, the extension segment, can be replaced either by zero filling or by a partial Fourier reconstruction from the data of the two recorded segments, since the data is point-symmetric with respect to the k-space center. This embodiment is described in FIG. 7. The entire raw-data space 90 depicted in FIG. 7 has the same extension in the kx direction as the acquired raw-data space 71. However, unlike the case with the acquired raw-data space 71, only the raw-data space 91 is acquired and the raw data in the extension segment 92 is reconstructed using the point symmetry and point reflection at the raw-data center 93. Unlike the case with the recorded raw-data space 71, now the same resolution is achieved in the kx direction as with an MR image that was only achieved using the acquired raw-data space 71, but the recording time is shorter in this case. In FIG. 7, a shorter signal readout time is used with each signal echo than with the embodiment in FIG. 6 and thus the overall recording time is reduced.
The method is summarized in FIG. 8. The method starts in step S80 and in step S81 zigzag-type segments of a trajectory are recorded, wherein following signal recording, the raw dataset is divided into even and odd echoes. The signal recording is performed such that the maximum signal intensity, the center of the signal echoes, is shifted, preferably by a quarter of the recorded raw-data space, in the direction of the readout axis. In step S82, an additional trajectory is recorded for the even and odd recorded raw-data points in each case using a parallel recording technique. In step S83, the raw-data space is extended in the readout direction, as was described in connection with FIGS. 4, 6 and 7. Then, in step S84, an iFT is performed each for each segment of the extended raw-data space along the phase-encoding direction. This means that an iFT is performed for example for the segment 71 in FIG. 6 and separately for the extension segment 75 or for the segment 51 in FIG. 4 and separately therefrom for the segment 55. Following this step S84 and the iFT, there is now a part-transformed raw-data space and a part-transformed extension segment. Part-transformed means that a Fourier transform was performed only along one direction, namely along the phase-encoding direction. In step S85, the part-transformed raw-data space and the part-transformed extension segment are now used to apply an inverse Fourier transform back into the raw-data space (inverse Fourier transform in the phase-encoding direction). After step S85, this results in a hybrid space and a hybrid extension segment. In step S86, the elements generated in step S85, the hybrid space and the hybrid extension segment are merged to form a total hybrid space. Following the merging of the data in step S86, in step S87 the Fourier transform can, as known, be performed to create the at least one MR image. The method ends in step S88.
The above-described method makes it possible to apply a method such as is known from parallel imaging with EPI data recorded on a zigzag trajectory. Herein, it is possible to use the iFT with which the even and odd echoes are phase-corrected individually so that they lie on Cartesian coordinates. This enables the achievement of higher acceleration factors and better image quality. It is also possible to perform partial-Fourier acquisition in the kx direction, which further reduces the acquisition time.
Although modifications and changes may be suggested by those skilled in the art, it is the intention of the Applicant to embody within the patent warranted hereon all changes and modifications as reasonably and properly come within the scope of the Applicant's contribution to the art.

Claims (15)

The invention claimed is:
1. A method for operating a magnetic resonance (MR) apparatus, comprising:
operating an MR data acquisition scanner to execute an echo-planar MR data acquisition sequence in which signal echoes are acquired by applying readout gradients along a readout axis, kx, and phase-encoding gradients along a phase-encoding axis, ky, with said signal echoes being entered asymmetrically into a memory organized as raw-data space yielding a zigzag-type undersampled trajectory such that said signal echoes are shifted in the raw-data space to cause a raw-data center associated with the raw-data space to be offset from a geometric center of the raw-data space with respect to the kx and ky axes, the shift in the raw data space being limited to a maximum of a quarter of the raw-data space in a direction of said readout axis;
in a computer, reconstructing an additional zigzag-type trajectory in the raw-data space by executing a parallel imaging reconstruction algorithm using the zigzag-type trajectory of the acquired signal echoes; and
in said computer, reconstructing at least one MR image from the acquired signal echoes in the raw-data space, using an interlaced Fourier transform and said additional zigzag-type trajectory.
2. A method as claimed in claim 1 comprising, before reconstructing said at least one MR image, extending raw-data points represented by said signal echoes in the raw-data space to form an extension segment, using a point symmetry of the raw-data space with respect to the center of raw-data space, and thereby producing an extended raw-data space comprising said raw-data space and said extension segment, with raw data in said extension segment being calculated using point reflection with respect to said center of the raw-data space.
3. A method as claimed in claim 2 comprising shifting said signal echoes in said direction of said readout axis by a quarter of the raw-data space,
wherein said extension segment is in said direction of the readout axis.
4. A method as claimed in claim 2, further comprising:
applying said interlaced Fourier transform to said signal echoes in the raw-data space to produce a hybrid space;
applying said interlaced Fourier transform to said extension segment to produce a hybrid extension segment;
merging said hybrid space and said hybrid extension segment to produce a total hybrid space; and
subjecting said total hybrid space to a Fourier transform along said readout axis to produce said at least one MR image.
5. A method as claimed in claim 2 wherein said at least one MR image has a resolution, and wherein a recording time required for acquiring said at least one MR image with said resolution is reduced by not filling the raw-data space completely with raw data that would be necessary for obtaining said resolution, but which corresponds to said extended raw-data space with said extension segment.
6. A method as claimed in claim 2 comprising recording said at least one MR image with a defined recording time, with raw data recorded in said defined recording time, without using said extension segment, being used to reconstruct an MR reference image with a defined resolution, said at least one MR image being reconstructed with said interlaced Fourier transformation using said extension segment so as to give said at least one MR image a higher resolution than said defined resolution.
7. A method as claimed in claim 1 comprising shifting said signal echoes along said readout axis only by an amount that makes a spacing of said zigzag-type undersampled trajectory, relative to the reconstructed additional zigzag-type trajectory, in a direction of the phase-encoding axis of said phase-encoding gradients, smaller than one sampling step in the ky-direction (Δky).
8. A method as claimed in claim 1 comprising reconstructing said additional zigzag-type trajectory using a Generalized Autocalibrating Partially Parallel Acquisition (GRAPPA) reconstruction algorithm.
9. A method as claimed in claim 1 comprising reconstructing said additional zigzag-type trajectory so that an incline of a last segment of said reconstructed zigzag-type trajectory is opposite, in a positive direction of said phase-encoding gradients, to an incline of a last segment of the acquired zigzag-type trajectory.
10. A method as claimed in claim 1, further comprising:
prior to reconstructing said at least one MR image, extending raw-data points represented by said signal echoes in the raw-data space to form an extension segment; and
shifting said signal echoes in said direction of said readout axis by a quarter of the raw-data space,
wherein said extension segment is in said direction of the readout axis.
11. A method as claimed in claim 1, further comprising:
prior to reconstructing said at least one MR image, extending raw-data points represented by said signal echoes in the raw-data space to form an extension segment;
applying said interlaced Fourier transform to said signal echoes in the raw-data space to produce a hybrid space; and
applying said interlaced Fourier transform to said extension segment to produce a hybrid extension segment.
12. A method as claimed in claim 11, further comprising:
merging said hybrid space and said hybrid extension segment to produce a total hybrid space; and
subjecting said total hybrid space to a Fourier transform along said readout axis to produce said at least one MR image.
13. A method as claimed in claim 1, wherein an average spacing between adjacent points of the raw-data space in a phase-encoding direction satisfies the Nyquist criterion.
14. A magnetic resonance (MR) apparatus comprising:
an MR data acquisition scanner;
a computer configured to operate said MR data acquisition scanner to execute an echo-planar MR data acquisition sequence in which signal echoes are acquired by applying readout gradients along a readout axis, kx, and phase-encoding gradients along a phase-encoding axis, ky, with said signal echoes being entered asymmetrically into a memory organized as raw-data space yielding a zigzag-type undersampled trajectory such that said signal echoes are shifted in the raw-data space to cause a raw-data center associated with the raw-data space to be offset from a geometric center of the raw-data space with respect to the kx and ky axes, the shift in the raw data space being limited to a maximum of a quarter of the raw-data space in a direction of said readout axis;
said computer being configured to reconstruct an additional zigzag-type trajectory in the raw-data space by executing a parallel imaging reconstruction algorithm using the zigzag-type trajectory of the acquired signal echoes; and
said computer being configured to reconstruct at least one MR image from the acquired signal echoes in the raw-data space, using an interlaced Fourier transform and said additional zigzag-type trajectory.
15. A non-transitory, computer-readable data storage medium encoded with programming instructions, said storage medium being loaded into a computer of a magnetic resonance (MR) apparatus comprising an MR data acquisition scanner, and said programming instructions causing said computer to:
operate said MR data acquisition scanner execute an echo-planar MR data acquisition sequence in which signal echoes are acquired by applying readout gradients along a readout axis, kx, and phase-encoding gradients along a phase-encoding axis, ky, with said signal echoes being entered asymmetrically into a memory organized as raw-data space yielding a zigzag-type undersampled trajectory such that said signal echoes are shifted in the raw-data space to cause a raw-data center associated with the raw-data space to be offset from a geometric center of the raw-data space with respect to the kx and ky axes, the shift in the raw data space being limited to a maximum of a quarter of the raw-data space in a direction of said readout axis;
reconstruct an additional zigzag-type trajectory in the raw-data space by executing a parallel imaging reconstruction algorithm using the zigzag-type trajectory of the acquired signal echoes; and
reconstruct at least one MR image from the acquired signal echoes in the raw-data space, using an interlaced Fourier transform and said additional zigzag-type trajectory.
US15/830,329 2016-12-06 2017-12-04 Magnetic resonance method and apparatus for echo-planar imaging with a zigzag-type trajectory in raw-data space Active 2038-06-09 US10775467B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE102016224175.3A DE102016224175B4 (en) 2016-12-06 2016-12-06 EPI with zigzag trajectory and parallel imaging
DE102016224175 2016-12-06
DE102016224175.3 2016-12-06

Publications (2)

Publication Number Publication Date
US20180156884A1 US20180156884A1 (en) 2018-06-07
US10775467B2 true US10775467B2 (en) 2020-09-15

Family

ID=62163663

Family Applications (1)

Application Number Title Priority Date Filing Date
US15/830,329 Active 2038-06-09 US10775467B2 (en) 2016-12-06 2017-12-04 Magnetic resonance method and apparatus for echo-planar imaging with a zigzag-type trajectory in raw-data space

Country Status (2)

Country Link
US (1) US10775467B2 (en)
DE (1) DE102016224175B4 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12108993B2 (en) 2021-08-05 2024-10-08 GE Precision Healthcare LLC Methods and system for guided device insertion during medical imaging

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017203936A1 (en) 2017-03-09 2018-09-13 Siemens Healthcare Gmbh Silent Echo Planar Imaging
US12299889B2 (en) 2018-06-14 2025-05-13 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for magnetic resonance image reconstruction
US10657646B2 (en) * 2018-06-14 2020-05-19 Uih America, Inc. Systems and methods for magnetic resonance image reconstruction
CN112557978B (en) * 2020-11-26 2022-03-25 上海交通大学 Multi-echo sampling method based on three-dimensional non-Cartesian trajectory
CN115598575B (en) * 2021-06-28 2025-10-28 上海联影医疗科技股份有限公司 Magnetic resonance imaging method, apparatus, computer equipment and storage medium
JP7623260B2 (en) * 2021-09-24 2025-01-28 富士フイルム株式会社 Magnetic resonance imaging apparatus and image processing method

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5672969A (en) * 1996-04-25 1997-09-30 General Electric Company Reduction of Nyquist ghost artifacts in oblique echo planar imaging
US6560353B1 (en) * 2000-03-21 2003-05-06 Magnetic Resonance Innovations, Inc. Method of MRI image reconstruction from partially acquired data in two or more dimensions using a multidimensional inverse transform technique
US20040204643A1 (en) * 2003-04-09 2004-10-14 Andrzej Jesmanowicz Perfusion magnetic resonance imaging using encoded RF tagging pulses
US20080024132A1 (en) * 2006-03-30 2008-01-31 Brau Anja C Method and appraratus of multi-coil mr imaging with hybrid space calibration
US20080068016A1 (en) * 2006-09-15 2008-03-20 Ajeetkumar Gaddipati System and method of accelerated mr propeller imaging
US20090072826A1 (en) * 2007-09-13 2009-03-19 The Board Of Trustees Of The Leland Stanford Junior University Magnetic resonance imaging with bipolar multi-echo sequences
US7592808B1 (en) * 2008-05-06 2009-09-22 General Electric Company System and method for reducing MR scan time using partial fourier acquisition and compressed sensing
US7941204B1 (en) * 2004-11-16 2011-05-10 Yi Wang Magnetic resonance imaging concepts
US20110148410A1 (en) * 2009-12-23 2011-06-23 Universitaetsklinikum Freiburg Method for data acquisition acceleration in magnetic resonance imaging (MRI) using receiver coil arrays and non-linear phase distributions
US8089278B1 (en) * 2008-04-11 2012-01-03 The Regents Of The University Of California Time-resolved contrast-enhanced magnetic resonance (MR) angiography
US8570034B2 (en) * 2009-03-23 2013-10-29 Siemens Aktiengesellschaft Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
US20140197834A1 (en) * 2013-01-14 2014-07-17 David Andrew Porter Method and magnetic resonance apparatus for echo planar imaging with data entry into k-space along a zigzag trajectory
US20140266191A1 (en) * 2013-03-12 2014-09-18 The Board Of Trustees Of The Leland Stanford Junior University Distortion-free magnetic resonance imaging near metallic implants
US20150253408A1 (en) 2014-03-04 2015-09-10 Siemens Aktiengesellschaft Method and apparatus for magnetic resonance imaging
US9234953B2 (en) * 2010-01-18 2016-01-12 Max-Planck-Gesellschaft Zur Foerderung Der Wissenschaften E.V. Method and device for magnetic resonance spectroscopic imaging
US10209331B2 (en) * 2015-05-11 2019-02-19 Universitätsspital Basel Magnetic resonance imaging method with asymmetric radial acquisition of K-space data

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5672969A (en) * 1996-04-25 1997-09-30 General Electric Company Reduction of Nyquist ghost artifacts in oblique echo planar imaging
US6560353B1 (en) * 2000-03-21 2003-05-06 Magnetic Resonance Innovations, Inc. Method of MRI image reconstruction from partially acquired data in two or more dimensions using a multidimensional inverse transform technique
US20040204643A1 (en) * 2003-04-09 2004-10-14 Andrzej Jesmanowicz Perfusion magnetic resonance imaging using encoded RF tagging pulses
US7941204B1 (en) * 2004-11-16 2011-05-10 Yi Wang Magnetic resonance imaging concepts
US20080024132A1 (en) * 2006-03-30 2008-01-31 Brau Anja C Method and appraratus of multi-coil mr imaging with hybrid space calibration
US20080068016A1 (en) * 2006-09-15 2008-03-20 Ajeetkumar Gaddipati System and method of accelerated mr propeller imaging
US20090072826A1 (en) * 2007-09-13 2009-03-19 The Board Of Trustees Of The Leland Stanford Junior University Magnetic resonance imaging with bipolar multi-echo sequences
US8089278B1 (en) * 2008-04-11 2012-01-03 The Regents Of The University Of California Time-resolved contrast-enhanced magnetic resonance (MR) angiography
US7592808B1 (en) * 2008-05-06 2009-09-22 General Electric Company System and method for reducing MR scan time using partial fourier acquisition and compressed sensing
US8570034B2 (en) * 2009-03-23 2013-10-29 Siemens Aktiengesellschaft Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
US20110148410A1 (en) * 2009-12-23 2011-06-23 Universitaetsklinikum Freiburg Method for data acquisition acceleration in magnetic resonance imaging (MRI) using receiver coil arrays and non-linear phase distributions
US9234953B2 (en) * 2010-01-18 2016-01-12 Max-Planck-Gesellschaft Zur Foerderung Der Wissenschaften E.V. Method and device for magnetic resonance spectroscopic imaging
US20140197834A1 (en) * 2013-01-14 2014-07-17 David Andrew Porter Method and magnetic resonance apparatus for echo planar imaging with data entry into k-space along a zigzag trajectory
US20140266191A1 (en) * 2013-03-12 2014-09-18 The Board Of Trustees Of The Leland Stanford Junior University Distortion-free magnetic resonance imaging near metallic implants
US20150253408A1 (en) 2014-03-04 2015-09-10 Siemens Aktiengesellschaft Method and apparatus for magnetic resonance imaging
US10185014B2 (en) * 2014-03-04 2019-01-22 Siemens Aktiengesellschaft Method and apparatus for EPI magnetic resonance with slew rate controlled and kspace entry optimized
US10209331B2 (en) * 2015-05-11 2019-02-19 Universitätsspital Basel Magnetic resonance imaging method with asymmetric radial acquisition of K-space data

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Scott B. Reeder et al., "Iterative Decomposition of Water and Fat With Echo Asymmetry and Least-Squares Estimation (IDEAL): Application With Fast Spin-Echo Imaging", Magnetic Resonance in Medicine, 54, 2005, pp. 636-644. (Year: 2005). *
Sekihara et al: "New Reconstruction Technique for Echo-Planar Imaging to Allow Combined Use of Odd and Even Numbered Echoes", Magnetic Resonance in Medicine 5, pp. 485-491; (1987).

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12108993B2 (en) 2021-08-05 2024-10-08 GE Precision Healthcare LLC Methods and system for guided device insertion during medical imaging

Also Published As

Publication number Publication date
DE102016224175A1 (en) 2018-06-07
US20180156884A1 (en) 2018-06-07
DE102016224175B4 (en) 2025-02-27

Similar Documents

Publication Publication Date Title
US10775467B2 (en) Magnetic resonance method and apparatus for echo-planar imaging with a zigzag-type trajectory in raw-data space
US9778336B2 (en) System and method for rapid, multi-shot segmented magnetic resonance imaging
US10823806B2 (en) Magnetic resonance imaging using dataset undersampling
US11860257B2 (en) Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset
EP2933651B1 (en) MRI Method of Hybrid Acquisition in 3D TSE
US9753110B2 (en) Method and magnetic resonance system for acquiring magnetic resonance data
US10203394B2 (en) Metal resistant MR imaging
US20130279786A1 (en) Rapid parallel reconstruction for arbitrary k-space trajectories
US10241184B2 (en) EPI ghost correction involving sense
Berman et al. Ultra‐high spatial resolution BOLD fMRI in humans using combined segmented‐accelerated VFA‐FLEET with a recursive RF pulse design
EP3988957B1 (en) Method for acquiring an mr-image dataset of at least two slices by means of simultaneous multislice excitation
US10502801B2 (en) Method and magnetic resonance apparatus for generating a weighting matrix for reducing artifacts with parallel imaging
US12523728B2 (en) Methods for acquiring a magnetic resonance image dataset and for generating a motion-corrected image dataset
US11543482B2 (en) Magnetic resonance imaging using motion-compensated image reconstruction
EP3683595B1 (en) Segment-specific undersampling in readout-segmented magnetic resonance imaging
US20180356488A1 (en) Method for calibration in a magnetic resonance imaging procedure
JP2019535435A (en) Propeller MR imaging
US10670681B2 (en) Magnetic resonance apparatus and method for parallel imaging with a reference data set for determining the weighting matrix
US11474178B2 (en) Method for generating a magnetic resonance image
US10663544B2 (en) Method and magnetic resonance apparatus for quiet echo-planar imaging
US11016155B2 (en) Method, apparatus and storage medium for recording a magnetic resonance dataset
US11029382B2 (en) Echo sharing in imaging sequences with multiple delays and multiple spin echoes
Tian et al. Overlap‐Kernel EPI: Estimating MRI Shot‐to‐Shot Phase Variations by Shifted‐Kernel Extraction From Overlap Regions at Arbitrary k‐Space Locations
US20240206819A1 (en) Generating a motion-corrected magnetic resonance image dataset
US20250093449A1 (en) Method for generating a series of magnetic resonance images with cross-frame iterative reconstruction

Legal Events

Date Code Title Description
FEPP Fee payment procedure

Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

AS Assignment

Owner name: SIEMENS HEALTHCARE GMBH, GERMANY

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:HEIDEMANN, ROBIN;LIEBIG, PATRICK, MR.;REEL/FRAME:045318/0032

Effective date: 20180110

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS

STPP Information on status: patent application and granting procedure in general

Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED

STCF Information on status: patent grant

Free format text: PATENTED CASE

AS Assignment

Owner name: SIEMENS HEALTHINEERS AG, GERMANY

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SIEMENS HEALTHCARE GMBH;REEL/FRAME:066267/0346

Effective date: 20231219

Owner name: SIEMENS HEALTHINEERS AG, GERMANY

Free format text: ASSIGNMENT OF ASSIGNOR'S INTEREST;ASSIGNOR:SIEMENS HEALTHCARE GMBH;REEL/FRAME:066267/0346

Effective date: 20231219

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 4