Method and apparatus recording two magnetic resonance images

11067655 · 2021-07-20

Assignee

Inventors

Cpc classification

International classification

Abstract

In a method for creating a first and a second image dataset of an examination object, a train of RF refocusing pulses are radiated into the examination object after the radiation of an RF excitation pulse to generate a spin echo signal after each radiated RF refocusing pulse, phase encoding gradients are activated for encoding the phases of the spin echo signals generated, and readout gradients are activated in each case in a readout window to read out the generated spin echo signals as measurement data. The readout windows alternately include a first time point at which the phases of the different spin species in the spin echo signal are the same, and a second time point at which the phases of the different spin species in the spin echo signal are not the same.

Claims

1. A method for creating a first image dataset and a second image dataset of an examination object in a measurement volume of a magnetic resonance system, said object comprising at least two different spin species, said method comprising: radiating an RF excitation pulse into the examination object; radiating a train of at least two RF refocusing pulses into the examination object after the radiation of the RF excitation pulse in order to generate a spin echo signal after each radiated RF refocusing pulse; activating phase encoding gradients for encoding the phases of the spin echo signals generated; activating readout gradients in each case in a readout window in order to read out the generated spin echo signals, as measurement data, with the readout windows alternately comprising a first time point at which the phases of the different spin species in the spin echo signal are the same, and a second time point at which the phases of the different spin species in the spin echo signal are not the same; storing the measurement data during the readout windows; reconstructing a first image from measurement data recorded in readout windows comprising the first time point; and reconstructing a second image from measurement data recorded in readout windows comprising the second time point.

2. The method as claimed in claim 1, wherein the phases of the different spin species in the spin echo signal are opposed at the second time point.

3. The method as claimed in claim 1, wherein the readout gradients are monopolar.

4. The method as claimed in claim 1, comprising generating the readout windows so the first and the second time point each lie in the middle of the respective readout windows.

5. The method as claimed in claim 1, comprising storing the measurement data that were read out during the readout windows that comprise the second time point as asymmetric measurement data.

6. The method as claimed in claim 1, comprising reading out and storing substantially all the measurement data for the reconstruction of the first image or the second image after a single RF excitation pulse.

7. The method as claimed in claim 1, comprising reading out and storing the measurement data for the reconstruction of the first image or the second image after a plurality of RF excitation pulses.

8. The method as claimed in claim 1, comprising encoding two successively generated spin echo signals with the same phase encoding gradients.

9. The method as claimed in claim 1, comprising encoding two successively generated spin echo signals with different phase encoding gradients.

10. The method as claimed in claim 1, wherein the measurement data stored for the reconstruction of the first image or the second image are incomplete with respect to the Nyquist criterion, and supplementing the incomplete data by supplementary procedures before said first and second images are reconstructed.

11. The method as claimed in claim 10, comprising acquiring calibration data for the supplementary procedures directly after a last RF refocusing pulse of a train of RF refocusing pulses.

12. The method as claimed in claim 10, comprising producing calibration data for said supplementary procedures using said measurement data.

13. The method as claimed in claim 1, comprising creating at least one specific image, which shows only one of the different spin species, from the first and second reconstructed images.

14. The method as claimed in claim 1, further comprising performing a supplementary procedure to supplement the measurement data before the first and second images are reconstructed.

15. The method as claimed in claim 14, comprising acquiring calibration data for the supplementary procedure after a last RF refocusing pulse of a train of RF refocusing pulses.

16. The method as claimed in claim 15, wherein the calibration data is acquired directly after the last RF refocusing pulse.

17. The method as claimed in claim 14, further comprising producing calibration data for the supplementary procedure based on the measurement data.

18. A magnetic resonance apparatus comprising: a magnetic resonance data acquisition scanner comprising a radio-frequency (RF) transmitting/receiving arrangement and a gradient coil arrangement; and a computer configured to: operate the magnetic resonance data acquisition unit in order to radiate an RF excitation pulse into the examination object with said RF transmitting/receiving arrangement; operate the magnetic resonance data acquisition unit in order to radiate a train of at least two RF refocusing pulses into the examination object with said RF transmitting/receiving arrangement, after the radiation of the RF excitation pulse in order to generate a spin echo signal after each radiated RF refocusing pulse; operate the magnetic resonance data acquisition unit in order to activate phase encoding gradients with said gradient coil arrangement for encoding the phases of the spin echo signals generated; operate the magnetic resonance data acquisition unit in order to activate readout gradients with said gradient coil arrangement in each case in a readout window in order to read out the generated spin echo signals, as measurement data, with the readout windows alternately comprising a first time point at which the phases of the different spin species in the spin echo signal are the same, and a second time point at which the phases of the different spin species in the spin echo signal are not the same; operate the magnetic resonance data acquisition unit in order to store the measurement data during the readout windows; reconstruct a first image from measurement data recorded in readout windows comprising the first time point; and reconstruct a second image from measurement data recorded in readout windows comprising the second time point.

19. A non-transitory, computer-readable data storage medium encoded with programming instructions, said storage medium being loaded into a computer of a magnetic resonance apparatus, and said programming instructions causing said computer to: operate the magnetic resonance apparatus in order to radiate an RF excitation pulse into the examination object; operate the magnetic resonance apparatus in order to radiate a train of at least two RF refocusing pulses into the examination object after the radiation of the RF excitation pulse in order to generate a spin echo signal after each radiated RF refocusing pulse; operate the magnetic resonance apparatus in order to activate phase encoding gradients for encoding the phases of the spin echo signals generated; operate the magnetic resonance apparatus in order to activate readout gradients in each case in a readout window in order to read out the generated spin echo signals, as measurement data, with the readout windows alternately comprising a first time point at which the phases of the different spin species in the spin echo signal are the same, and a second time point at which the phases of the different spin species in the spin echo signal are not the same; operate the magnetic resonance apparatus in order to store the measurement data during the readout windows; reconstruct a first image from measurement data recorded in readout windows comprising the first time point; and reconstruct a second image from measurement data recorded in readout windows comprising the second time point.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic depiction of a HASTE sequence.

(2) FIG. 2 shows a schematic sequence diagram of the method according to the invention.

(3) FIG. 3 shows an example of a sequence pattern.

(4) FIG. 4 shows an example of a possible sampling pattern of k-space with the sequence pattern according to FIG. 2.

(5) FIG. 5 shows a further example of a sequence pattern.

(6) FIG. 6 shows an example of a possible sampling pattern of k-space with the sequence pattern according to FIG. 4.

(7) FIG. 7-9 are schematic depictions of possible time coordinations of RF refocusing pulses and readout windows,

(8) FIG. 10 is a schematic depiction of a magnetic resonance apparatus according to the invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

(9) FIG. 2 is a schematic sequence diagram of a method according to the invention for creating a first and a second image dataset of an examination object arranged in a measurement volume of a magnetic resonance system, which comprises at least two different spin species.

(10) An RF excitation pulse is radiated into an examination object U (block 201). The RF excitation pulse manipulates the magnetization of spins in a measurement volume FOV in the examination object P.

(11) After the radiation 201 of the RF excitation pulse RF1, a train Z of at least two RF refocusing pulses RF2 is radiated into the examination object U, wherein a spin echo signal is generated after each radiated RF refocusing pulse RF2.

(12) Herein, an RF refocusing pulse RF2 can be radiated (block 203) that generates a spin echo signal E1, which is recorded in a readout window AF as measurement data MDS1 (block 205). Following this, a further RF refocusing pulse RF2 is radiated (block 207), which generates a spin echo signal E2, which in is recorded in a readout window AF as measurement data MDS2 (block 209).

(13) The readout windows belonging to successive RF refocusing pulses RF2 are arranged temporally such that the readout windows alternately comprise a first time point at which the phases of the different spin species in the spin echo signal Ei are the same and a second time point at which the phases of the different spin species in the spin echo signal Eo are not same, but, for example, are opposed. Hence, when a read-out spin echo signal E1 is a spin echo signal Ei with which the phases of the different spin species are the same, the following read-out spin echo signal E2 is a spin echo signal Eo with which the phases of the different spin species not are the same; and vice versa, when a read-out spin echo signal E1 is a spin echo signal Eo with which the phases of the different spin species not are the same, the following read-out spin echo signal E2 is a spin echo signal Ei with which the phases of the different spin species are the same.

(14) Before, during and after the radiation of the RF pulses and the recording of the measurement data in the readout windows (blocks 201-209), gradients for the spatial encoding, in particular phase encoding gradients and readout gradients, and possibly also slice selection gradients, can be activated (block 211).

(15) Examples of sequence patterns applicable in this context are depicted in FIGS. 3 and 5. These depict the respective temporal courses of radio-frequency pulses and signals (RF), activated gradients in the readout direction (GR), activated gradients in the phase encoding direction (GP) and the data recording (ADC) in the relationship to one another. Slice selection gradients are not depicted, but is also possible for slice selection gradients to be activated, in particular simultaneously with the RF excitation pulses RF1 used and/or simultaneously with the RF refocusing pulses RF2 used, for example for two-dimensional measuring methods with which, for example, measurement data are only to be recorded within the slices defined by the slice selection gradients and respective RF pulses. For 3D measuring methods, such as, for example, SPACE, slice selection gradients can be dispensed with.

(16) After an RF excitation pulse RF1, a train of RF refocusing pulses RF2 is radiated at regular time spacings. In the example depicted, the spin echo signal Ei radiated by a first radiated RF refocusing pulse RF2 is read out in a readout window AF if the spins of the different spin species are in phase.

(17) After the reading out of the spin echo signal Ei generated, a further RF refocusing pulse RF2 is radiated. The spin echo signal Eo generated thereby is now shifted, not according to the time duration ESi expected for an “in-phase” echo Ei between two “in-phase” echoes (“echo spacing”), but by a time difference Δt if the spins of the different spin species are read out in opposed phase as an “opposed-phase” echo during an associated readout window AF. The time difference is obtained from the chemical shift between the different spin species.

(18) The spin echo signal Ei generated by the following RF refocusing pulse RF2 is read out again during a readout window AF if the spins of the different spin species are in phase so that a spin echo signal Ei (“in-phase”) and a spin echo signal Eo (“opposed-phase”) are always recorded in alternation. In the examples depicted, the readout windows AF are arranged such that the first and second time point in each case lie in the middle of the respective readout windows, i.e. that the “in-phase” or “opposed-phase” condition are in each case satisfied in the temporal middle of the readout window. In this way, the measurement data is acquired symmetrically in k-space.

(19) To keep the sequence as short as possible, the regular time spacing of the RF refocusing pulses radiated after a common RF excitation pulse RF1 can be selected such that it corresponds to the short “in-phase” echo spacing ESi.

(20) During the reading out of the spin echoes, a monopolar, constant readout pulse can be activated in each case. This enables chemical shift artifacts to be avoided in a simple manner. Before and after the readout windows, further gradients can be activated in the readout direction in order, for example, to destroy unwanted signals (“crusher” gradients).

(21) For the spatial encoding of the spin echo signals Ei, Eo recorded in the readout windows as measurement data MDS1, MDS2, phase encoding gradients are also activated in the phase encoding direction GP.

(22) As in the example depicted in FIG. 3, in each case two successive spin echo signals Ei and Eo can be recorded with the same phase encoding gradients in the phase encoding direction GR, between the respective RF refocusing pulse RF2 and the following spin echo signal Ei or Eo. This enables each phase encoding step to be recorded in each case for both spin echo signals to be recorded, Ei and Eo. The fact that “in-phase” spin echo signals Ei and “opposed-phase” spin echo signals Eo are hence recorded with only a small time offset with the same spatial encoding means that the images BDS1 and BDS2 reconstructed later from the measurement data MDS1 and MDS2 can be processed particularly well with one another to form images that depict only one of the spin species without it being necessary to take account of artifacts caused by a movement of the examination object.

(23) The example of a sequence pattern shown in FIG. 5 only differs from that shown in FIG. 3 only in a different type of phase encoding in that two successively generated spin echo signals Ei and Eo are in each case encoded with different phase encoding gradients in the phase encoding direction GR. For each spin echo signal to be read out of an RF refocusing pulse RF2 of a train of RF refocusing pulses RF2 after a common RF excitation pulse RF1, in each case another phase encoding gradient is activated, from a series of phase encoding gradients required to satisfy the Nyquist condition. The “another” phase encoding gradient can be the next gradient in a series ranked according to amplitude.

(24) If, overall, more than one RF excitation pulse RF1 followed by a train of RF refocusing pulses RF2 is radiated for each measurement volume to be recorded, for example for each slice of the examination object to be recorded, the phase encoding can take place such that, during the recording of the spin echo signals of every second train of RF refocusing pulses RF2, the “in-phase” spin echo signals Ei are in each case encoded with a phase encoding gradient with which an “opposed-phase” spin echo signal Eo was encoded in the preceding train of RF refocusing pulses RF2. Also, the “opposed-phase” spin echo signals Eo are in each case encoded with a phase encoding gradient with which an “in-phase” spin echo signal Ei was encoded in the preceding train of RF refocusing pulses RF2. In this way, overall all “in-phase” spin echo signals Ei and all “opposed-phase” spin echo signals Eo can be encoded with all the phase encoding gradients required by the Nyquist theorem.

(25) The train of RF refocusing pulses RF2 includes as many RF refocusing pulses RF2 as spin echo signals E1 and E2 are to be generated and read out. This can, for example, be checked by a query 213 that ensures that further RF refocusing pulses RF2 can be radiated if not all the desired spin echo signals E1 and E2 have been recorded and stored as measurement data MDS1 and MDS2 (query 213, “n”). If desired, it is also possible to initiate the radiation of a new RF excitation pulse RF1 followed by a further train of RF refocusing pulses RF2.

(26) If all the desired/to-be-recorded measurement data MDS1 and MDS2 has been read out (query 213, “y”), a first image BDS1 can be reconstructed from measurement data MDS1 recorded in readout windows comprising the first time point recorded and a second image BDS2 reconstructed from measurement data MDS2 recorded in readout windows comprising the second time point (block 215).

(27) At least one specific image SB showing only one of the different spin species can be reconstructed from the first and second reconstructed images BDS1 and BDS2. For example, a fat and/or a water image can be created in a manner usual with Dixon methods from the first image BDS1 and the second image BDS2.

(28) Herein, the entity of k-space can be recorded (filled) for the measurement data MDS1 and MDS2 in each case. For example, to this end, it is possible with a sequence pattern according to FIG. 3, for each phase encoding gradient required to satisfy the Nyquist condition to be performed for each type of spin echo signal to be recorded. A complete recording of the measurement data MDS1 and MDS2 of this kind can improve the signal-to-noise ratio.

(29) It is also conceivable, so as to shorten the necessary length of the train of RF refocusing pulses RF2, and hence the measuring time, to read out only a part of the measurement data required by the Nyquist theorem, i.e. to record the measurement data MDS1 and MDS2 only incompletely (undersampled). The “missing” measurement data can be supplemented by means of supplementary procedures (block 217), before images BDS1, BDS2 are reconstructed. To this end, it may be necessary to determine calibration data KD (block 219).

(30) For an incomplete reading out of the measurement data MDS1 and MDS2 of this kind, it is possible, for example with a sequence pattern according to FIG. 3, to in each case fill every second k-space line to be recorded for the measurement data MDS1 and for the measurement data MDS2. As a result, the sampling patterns depicted in FIG. 4 are obtained for the measurement data MDS1 and MDS2. The recorded measurement data in each case are depicted as continuous k-space lines and omitted, non-recorded measurement data depicted as dashed k-space lines. Hence, the sampling patterns depicted correspond to conventional single-slice GRAPPA sampling in which every second k-space line in one spatial direction (here: k.sub.y direction) is omitted and hence only half k-space points are measured, corresponding to an acceleration factor of 2.

(31) It is also possible, for example with a sequence pattern according to FIG. 5, for in each case every second k-space line to be recorded for the measurement data MDS1 and the measurement data MDS2, such that sampling patterns for the measurement data MDS1 and MDS2 such as depicted by way of example in FIG. 6 are obtained. There, once again, in each case recorded measurement data is depicted as continuous k-space lines and omitted, non-recorded measurement data as dashed k-space lines. In contrast to the sampling patterns shown in FIG. 4, here k-space lines recorded in the measurement data MDS1 are those that were not recorded in the measurement data MDS2, and vice versa.

(32) The two sets of measurement data MDS1 and MDS2 can in each case be supplemented with supplementary methods, such as afore-mentioned parallel acquisition methods, for example GRAPPA.

(33) If a supplementary method is used that requires calibration data, the calibration data can be recorded directly after a last RF refocusing pulse of a train of RF refocusing pulses and/or created based on the recorded measurement data recorded.

(34) If calibration data of this kind are recorded with a different sequence type than the measurement data to be completed of the actual measurement, due to possible possibly different contrasts and/or different sensitivities to different disruptions, this can result can in artifacts in the supplementation of the measurement data. Therefore, such data are preferably also recorded with a (T)SE sequence.

(35) The time point at which the calibration data are recorded compared to the measurement data, can also result in artifacts during the supplementation, for example due to (unwanted) movements of the examination object. Therefore, the recording of the calibration data can be performed following the recording of measurement data in order to reduce sensitivity to physiological movements, as described, for example, in DE102016207641.

(36) The two sets of measurement data MDS1 and MDS2 can also be advantageously combined with one another for a supplementary method of this kind, as described in DE102016219052B3, for example. In this way, the combination of the two sets of measurement data MDS1 and MDS2 enables a higher SNR in the images BDS1 and BDS2 reconstructed from measurement data recorded after supplementation by means of the calibration data and higher acceleration factors in the context of the parallel acquisition technique with which the measurement data MDS1 and MDS2 are recorded, even if the technique is only implemented with only a few reception channels (receiving coils). Independent sets of “in-phase” and “opposed-phase” calibration data are to be created for each set of measurement data MDS1 and MDS2 in each case. Due to the offset of the read-out k-space lines in the two sets of measurement data MDS1 and MDS2, a sampling pattern according to FIG. 6 enables particularly high acceleration factors, since the two sets of measurement data can almost be considered to be different slices as in a slice multiplexing method. To be precise, the depicted offset of k-space lines has the result that k-space data contained in MDS1 have a different image space phase than k-space data contained in MDS2. Hence, sampling patterns are obtained that resemble the sampling patterns in, for example, the slice multiplexing technique CAIPIRINHA (see, for example, Breuer et al. “Controlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA) for Multi-Slice Imaging”, Magnetic Resonance in Medicine 53, pp 684-691 (2005)), as a result, the techniques used there can also be used here, wherein the different phase positions of the spins (“in-phase” and “opposed-phase”) replace the different layers of the slice multiplexing technique.

(37) With a single-shot method, such as HASTE, substantially all measurement data MDS1 and MDS2 for the reconstruction of a first and second image after a single RF excitation pulse can be generated and recorded in an advantageously short time. Hence, the measuring time can be reduced so that in particular even a breath-holding phase of a patient to be examined as an examination object U can be sufficient to record all the measurement data. In addition, the short measuring time minimizes the sensitivity of the measurement data to motion artifacts.

(38) FIGS. 7 to 9 are schematic depictions comparing possible temporal coordinations of RF refocusing pulses and readout windows for reading out “opposed-phase” spin echo signals Eo, wherein, for clarity, the readout windows AF are also depicted on the readout gradient line GR.

(39) FIG. 7 again depicts a situation such as that already shown in FIGS. 3 and 5: after an RF refocusing pulse, a spin echo signal Eo is read out that, as described above, is shifted by the aforementioned time difference Δt compared to an “in-phase” spin echo signal Ei, while a readout gradient is activated in the readout direction GR. This is followed by a further RF refocusing pulse RF2. Herein, the readout window AF is arranged such that the second time point (=when the excited different spin species have an opposed phase) lies in the middle of the associated readout window AF.

(40) FIGS. 8 and 9 illustrate strategies for reducing the time spacing of successive RF refocusing pulses by a time period Δτ. The reduction Δτ achieved in each case in the examples shown in FIGS. 8 and 9 is only shown as the same length in order to simplify the depiction. The actual reductions Δτ achieved can quite feasibly differ and the strategies described below can also be combined. The reduction of the time spacing of successive RF refocusing pulses enables the measuring time to be further reduced.

(41) The reduction strategy shown in FIG. 8 includes reading out the spin echo signals not exactly at the respective first or second time point, but in particular before the respective time points are reached. To this end, FIG. 8 shows an example of how a “opposed-phase” spin echo signal Eo′ is read out in a readout window AF with a time difference Δ′, which is shorter than the time difference Δt, after the first time point (the “in-phase” condition is achieved when the different excited spin species are in phase). The readout window AF is again arranged such that the second time point lies in the middle of the associated readout window AF. Dixon methods are known which can also be performed with measured values of this kind that are not recorded exactly at a first or second time point.

(42) The reduction strategy depicted in FIG. 9 uses an asymmetric sampling pattern. The spin echo signal is only partially read out. In the example depicted, the “opposed-phase” spin echo signal Eo is read out such that the measured values are only read out on a shortened k-space line with an incomplete right half. The recorded measurement data are stored as asymmetric measurement data. Herein, the readout window AF is hence arranged such that the second time point (fulfilment of the “opposed-phase” condition) lies in the second half of the readout window and hence more measurement data is acquired before the fulfilment of the “opposed-phase” condition than after the fulfilment of the “opposed-phase” condition. In order, nevertheless, to maintain the Carr-Purcell-Meiboom-Gill condition for the pulse sequence, after the readout window a gradient with a higher amplitude can be activated, which can also be combined with a subsequent crusher gradient.

(43) FIG. 10 is a schematic depiction of a magnetic resonance apparatus 1 according to the invention. This has a scanner 3 with a basic field magnet that generates the basic magnetic field, a gradient unit 5 that generates the gradient fields, a radio-frequency unit 7 that radiates and receives radio-frequency signals, and a control computer 9 designed (programmed or configured) to perform the method according to the invention. FIG. 10 is only a basic schematic depiction of these subunits of the magnetic resonance apparatus 1. In particular, the radio-frequency unit 7 can have a number of subunits, in particular at least two coils, such as the schematically depicted coils 7.1 and 7.2, which can be designed either only to transmit radio-frequency signals or only to receive the triggered radio-frequency signals, or for both.

(44) For the examination of an examination object U, for example a patient or a phantom, this object U can be introduced into the measurement volume FOV of the scanner 3 on a bed L.

(45) The control computer 9 controls the magnetic resonance apparatus 1 and can in particular control the gradient unit 5 by a gradient controller 5′ and the radio-frequency unit 7 by means of a radio-frequency transmitting/receiving controller 7′. The radio-frequency unit 7 can have multiple channels in which signals can be transmitted or received.

(46) The radio-frequency unit 7 together with its radio-frequency transmitting/receiving controller 7′ is responsible for generating and radiating (transmitting) a radio-frequency alternating field for manipulating the spins in a region to be manipulated, for example a slice, of the examination object U. The center-frequency of the radio-frequency alternating field, also called the B1 field, should lie close to the resonance frequency of the spins to be manipulated of the desired spin species. To generate the B1 field, currents controlled by the radio-frequency transmitting/receiving controller 7′ are applied to the RF coils in the radio-frequency unit 7.

(47) The control computer 9 also has a coordination processor 15 that coordinates readout windows with the phase positions of the different spin species in the examination object U.

(48) A control processor 13 of the control computer 9 is designed to carry out all the computing operations required for the necessary measurements and determinations. Interim results required to this end or determined herein, can be stored in a memory S of the control computer 9. The units depicted are not necessarily physically separate units—they merely represent a subdivision into coherent units, which can also be implemented, for example, in just a few or even only a single physical unit.

(49) An operator user can use an input/output device E/A of the magnetic resonance apparatus 1 to enter control commands into the magnetic resonance apparatus 1 and/or to view results of the control computer 9 such as image data.

(50) A non-transitory, electronically readable data carrier 26 can be provided with electronically readable control information (program code) stored therein. The program code cause the control computer 9 to carry out the described method when the data carrier 26 is loaded into the control computer 9.

(51) 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.