A Magnetic Resonance (MR) scanner generates images of patient anatomy using timed sequences of RF pulses. MR imaging is useful in scenarios requiring high contrast between different soft tissues. Cardiac MR (CMR) imaging is increasingly used to non-invasively evaluate myocardial structure and function without requiring ionizing radiation as in other imaging modalities.
CMR imaging often uses a segmented inversion recovery (IR)-prepared sequence which follows the injection of a T1-shortening Gadolinium-based contrast agent. Such a sequence, known as a “delayed enhancement” (DE) or “late Gadolinium enhancement” (LGE) sequence, produces T1-weighted images in which viable (i.e., healthy, living) myocardium appears dark and infarcted (i.e., dead) myocardium appears bright. The sequence also exhibits high spatial resolution and excellent sensitivity to infarcted myocardium.
Generation of a high-resolution LGE image requires more data than can be acquired in a single data readout (i.e., shot), especially for images of moving structures such as the heart. For example, an image produced from the data of a single shot would exhibit poor temporal resolution, and intricate cardiac features would appear blurred. Therefore, each image is reconstructed from raw data acquired during multiple shots wherein each shot is taken with precise temporal resolution (i.e., in the same cardiac phase but in a different heartbeat). Each of the multiple shots acquires a different subset of all the lines of the raw data space known as “k-space”.
More specifically, the data acquired by a single shot is not image data but is frequency-encoded data in k-space. k-space is defined by two orthogonal directions, the phase encoding (PE) direction and the frequency encoding (FE) direction. Different regions of k-space contain different image properties. For example, the center of k-space contains image brightness and contrast information, and the edges of k-space contain sharpness and detail information. A two-dimensional Fourier Transform (FT) may be used to calculate an image from acquired k-space data.
The relationship between the temporal order in which PE lines are acquired and their specific location in k-space is known as “reordering”. In the context of IR-prepared imaging, reordering affects the image contrast, specifically T1 contrast, and the robustness of the acquisition to motion and flow. Conventional types of reordering include continuous reordering and interleaved reordering.
According to contiguous reordering, each shot acquires a set or segment of m contiguous PE lines, such that two segments from consecutive shots share one border in k-space.
Using contiguous reordering, the different modulation corresponding to each PE line within a segment creates a periodic pattern of the modulation function in the PE direction as shown in
As illustrated in
As illustrated in
However, interleaved reordering is extremely sensitive to motion and flow of the objects to be imaged. Ghosting artifacts may result from even the slightest positional shift between shots, for example due to an imperfect breath hold. Ghosting may also occur in the presence of arrythmia or poor ECG-triggering, which may cause one or more shots of the sequence to be acquired during a different cardiac phase than the other shots. By way of explanation, positional shift modulates the amplitude and phase of the raw data and, if the positional shift is different for each shot and the data acquired by each shot is combined by interleaving, artificial frequencies are introduced into k-space. These artificial frequencies are reflected as replication of structures in image space, i.e., ghosting.
Correction of such interleaved-reordered data in k-space or image space has not been successfully performed and may be impossible, partially because the motion causing the positional shifts is non-rigid. Whereas non-rigid motion correction exists for single shot images, non-rigid motion correction is incompatible with images reconstructed from interleaved data. Each set of interleaved k-space data reconstructs into an aliased image wherein different sections of the imaged field of view in the PE direction are superimposed on one another, known as aliasing. Since each region has experienced a different non-rigid positional shift, no single motion correction of an aliased image can correct all sections. Moreover, due to the sensitivity of interleaved PE lines to amplitude and phase changes, motion-corrected interleaved data typically produces ghosting in image space.
Conversely, contiguous reordering is quite motion-robust. However, as described above, the combination of IR-preparation with contiguous reordering causes artifacts unrelated to motion. Although interleaved reordering may be employed to reduce these artifacts, this approach results in motion-sensitivity, as also described above.
Systems to efficiently generate high-quality, motion-robust images using an IR-prepared segmented pulse sequence are desired.
The following description is provided to enable any person in the art to make and use the described embodiments. Various modifications will remain apparent to those in the art.
Embodiments provide a reordering and inversion time (TI) scheme that addresses the motion- and flow-sensitivity of conventional techniques, while being compatible with IR-preparation sequences.
It is traditionally counterintuitive to use contiguous reordering for a segmented IR-prepared sequence due to the above-described incompatibility of contiguous reordering with IR-preparation. Some embodiments overcome this incompatibility by using a different TI for each shot and by acquiring, in each shot, a different number of “dummy” segments before acquisition of a shot-specific segment to be used in image reconstruction (i.e., a k-space segment). The TI of each shot refers to the time between the peak of an IR preparation pulse and a time at a midpoint of the acquisition of the k-space segment of the shot. Use of a different TI for each shot is also counterintuitive to the common belief that the same TI should be used for all shots of a segmented sequence.
These modifications to a standard contiguous reordering scheme prevent both the periodic modulation pattern and the discontinuities in k-space resulting from standard IR-prepared contiguous reordering and shown in
Embodiments described herein, in essence, reproduce the magnetization behavior of a single shot acquisition (i.e., an acquisition in which all k-space lines desired for imaging are acquired in a single segment across k-space). However, a true single shot acquisition with the same high spatial resolution as the embodiments described herein would exhibit very poor temporal resolution, taking so much time as to include the cardiac motion of a significant portion of a cardiac cycle, for example. The heart would therefore be blurry or unrecognizable in the resulting image. It would not be practical to motion-correct such a poor-quality single shot image because the large number of time points represented in the image would not be separable in post-processing.
Acquisition of a DCS or DUS according to some embodiments may utilize any suitable type of readout pulses, including but not limited to gradient echo, non-interleaved epi, steady-state free precession, and turbo spin-echo.
According to the
For each shot (#1-#5) in the pulse sequence, the time TD from the peak of the IR pulse to acquisition of the first segment (image or dummy) is fixed. Shot #1 acquires no DUS before its DCS, shot #2 acquires 1 DUS before its DCS, and shot #5 acquires 4 DUS before its DCS. Considering the generic case of n shots resulting in n k-space segments, shot #n acquires (n−1) DUS before its DCS. The time between the IR pulse of a shot to acquisition of the DCS of the shot is therefore different for each shot. Accordingly, the time delay between the ECG's R-wave and the IR pulse is different for each shot, unlike in known sequences using interleaved or continuous reordering.
The IR pulse in some embodiments may comprise a spatially non-selective IR pulse, a spatially selective IR pulse, a wideband IR pulse, a double-IR dark-blood preparation pulse, or any other suitable IR pulse. In some embodiments, a magnetization transfer preparation or a T2-preparation is applied before or after an IR pulse.
The T1-recovery curves of
In the first R-R interval of
In IR-prepared imaging, TI is defined as the time from the IR pulse to the acquisition of the center line of k-space. A TI may be selected so as to achieve a desired T1-contrast, since the center of k-space contains image brightness and contrast information. In the example of
According to some embodiments, the shot which acquires the center line of k-space determines the time delay TD between the IR pulse and the first segment of each shot. For example, a TI is selected based on a desired T1-contrast. With reference to shot #3 of
The timing of all other shots is derived therefrom by keeping the acquisition of the DCS of each shot at the same point of the R-R interval (at diastole), inserting acquisition of the required number of preceding dummy segments, and placing the IR pulse at a time TD prior to acquisition of the first segment of the shot. Shot #n is the longest in duration and determines the minimum needed acquisition window. All shots may ideally fit into one R-R interval, but imaging over two or more R-R intervals is possible, which would result in a longer shot and thus higher spatial resolution.
Process 500 and all other processes mentioned herein may be embodied in executable program code read from one or more of non-transitory computer-readable media, such as a disk-based or solid-state hard drive, a DVD-ROM, a Flash drive, and a magnetic tape, and then stored in a compressed, uncompiled and/or encrypted format. In some embodiments, hard-wired circuitry may be used in place of, or in combination with, program code for implementation of processes according to some embodiments. Embodiments are therefore not limited to any specific combination of hardware and software.
At S510, the shot number n is initialized to 1, where the total number of shots to execute is N. Flow pauses at S520 until an R peak (i.e., the start of a heartbeat interval) is detected. At S530, an IR pulse is applied as is known in the art at a timing specific to the first (n=1) shot. The timing is determined based on the predetermined timing of the DCS of the shot (i.e., at diastole), on the number of DUSs preceding the DCS and their overall acquisition duration, and on the time TD calculated as described above. Accordingly, each IR pulse of the
Next, process 500 waits at S540 for time TD. Once time TD expires, contiguous k-space lines of the first n−1 segments of k-space are acquired. These n−1 segments are referred to as dummy segments above. Since n currently equals 1 in the present example, no dummy segments are acquired at S550, as shown in
Contiguous k-space lines of the nth segment of k-space are acquired at S570. Continuing the present example, contiguous k-space lines of segment 1 of k-space are acquired during shot #1, as also shown in
The contiguous k-space lines acquired at S560 are stored at S570 and it is then determined at S580 whether N shots have been executed. If not, n is incremented at S585 and flow returns to S520 to wait for a next R peak. According to some embodiments, flow pauses at S520 until two R peaks are detected in order to allow a magnetization recovery heartbeat interval to elapse between consecutive shots as illustrated in
Once a suitable R peak is detected, an inversion recovery pulse is applied at a timing specific to shot #2. This timing must account for acquisition of the one DUS of shot #2 and therefore the inversion pulse is applied earlier in the R-R interval than in shot #1. Flow pauses at S540 for the same time TD as used in shot #1.
One (n−1=1) DUS is acquired at S550, followed by acquisition at S560 of contiguous k-space lines of the second segment of k-space. The contiguous k-space lines of the second segment of k-space are stored at S570.
Flow continues in this manner until five shots have been executed, each with the same TD and different TIs. Each shot acquires and stores contiguous k-space lines of a respective segment of k-space. Accordingly, at S590, an image is reconstructed based on the stored contiguous k-space lines as is known in the art.
Returning to
The concept of improving image quality using a different TI during acquisition of the DCS of each shot can also be applied to interleaved reordering. This approach may reduce or remove the above-mentioned small steps in the modulation function and thus prevent the slight ghosting which is inherent to interleaved reordering.
Center shot #3 includes a delay TDC from the peak of its inversion pulse to the beginning of acquisition of its DCS. TDC is equal to a pre-selected TI minus half the number of lines acquired per shot multiplied by the echo spacing (ES). The time delay TD between the IR pulse and acquisition of the DCS of shot #2 is TDC−1ΔTI, where ΔTI is equal to ES/N, and where N is the number of shots in the sequence.
Continuing this logic, the time delay TD between the IR pulse and acquisition of the DCS of shot #1 is TDC−2ΔTI. Conversely, the time delay TD for shot #4 is increased by ΔTI to TDC+1ΔTI, and for shot #5 by 2ΔTI to TDC+2ΔTI. The modulation function of interleaved reordering using the foregoing timings is shown in
Process 800 of
Process 800 waits at S840 for a shot-specific delay time before applying its DCS at S850. Accordingly, unlike the
Next, at S850, a segment is applied to acquire interleaved k-space lines of an nth segment of k-space as described above. The interleaved k-space lines are stored at S860, and flow continues in this manner for each successive shot until it is determined at S870 that M shots have been executed. As described above with respect to process 500, flow may pause at S820 until two R peaks are detected in order to allow a magnetization recovery heartbeat interval to elapse between consecutive shots as illustrated in
The order in which segments of k-space lines are acquired can be varied from that shown in
Advantageously, after the first shot, a full field of view intermediate image can be reconstructed based on the acquired center segment of k-space lines. This image is of low spatial resolution in the PE direction, but without fold-in artifacts. After acquisition of segments 2 and 4 by shots #2 and #3, segments 2, 3 and 4 may be used to reconstruct a second intermediate image with triple the PE-resolution of the first image. A full PE-resolution image can be reconstructed after executing shots #4 and #5 to acquire segments 1 and 5. The intermediate images may be displayed during acquisition in real-time as they are reconstructed.
The sequence of
The motion-robust embodiments described herein can be combined with the acquisition of a respiratory navigator, a reference line, or a low-resolution reference image. This allows free-breathing acquisitions without ghosting, which would be observed if using interleaved reordering.
Although the sequences described herein use an ECG signal as a triggering source, embodiments are compatible with any other means of triggering, such as pulse oximetry, a respiratory belt, an inductively coupled device, or a Pilot Tone device. Embodiments may be implemented without triggering by using a constant effective time-to-repeat (effective TR) in place of the above-described R-R interval.
According to MR techniques, a substance (e.g., human tissue) is subjected to a main polarizing magnetic field (i.e., B0), causing the individual magnetic moments of the nuclear spins in the substance to process about the polarizing field in random order at their characteristic Larmor frequency, in an attempt to align with the field. A net magnetic moment Mz is produced in the direction of the polarizing field, and the randomly-oriented magnetic components in the perpendicular plane (the x-y plane) cancel out one another.
The substance is then subjected to an excitation field (i.e., B1) created by emission of a radiofrequency (RF) pulse, which is in the x-y plane and near the Larmor frequency, causing the net aligned magnetic moment Mz to rotate into the x-y plane so as to produce a net transverse magnetic moment Mt, which is rotating, or spinning, in the x-y plane at the Larmor frequency. The excitation field is terminated, and signals are emitted by the excited spins as they return to their pre-excitation field state. The emitted signals are detected, digitized and processed to reconstruct an image or a spectrum using one of many well-known MR techniques.
Gradient coils 6 produce magnetic field gradients Gx, Gy, and Gz which are used for position-encoding NMR signals. The magnetic field gradients Gx, Gy, and Gz distort the main magnetic field in a predictable way so that the Larmor frequency of nuclei within the main magnetic field varies as a function of position. Accordingly, an excitation field B1 which is near a particular Larmor frequency will tip the net aligned moment Mz of those nuclei located at field positions which correspond to the particular Larmor frequency, and signals will be emitted only by those nuclei after the excitation field B1 is terminated.
Gradient coils 6 may consist of three windings, for example, each of which is supplied with current by an amplifier 8a-8c in order to generate a linear gradient field in its respective Cartesian direction (i.e., x, y, or z). Each amplifier 8a-8c includes a digital-analog converter 9a-9c which is controlled by a sequence controller 10 to generate desired gradient pulses at prescribed times.
Sequence controller 10 also controls the generation of RF pulses by RF system 11 and RF power amplifier 12. RF system 11 and RF power amplifier 12 are responsive to a scan prescription and direction from sequence controller 10 to produce RF pulses of the desired frequency, phase, and pulse amplitude waveform. The generated RF pulses may be applied to the whole of RF coil 7 or to one or more local coils or coil arrays. RF coil 7 converts the RF pulses emitted by RF power amplifier 12, via multiplexer 13, into a magnetic alternating field in order to excite the nuclei and align the nuclear spins of the object to be examined or the region of the object to be examined. As mentioned above, RF pulses may be emitted in a magnetization preparation step in order to enhance or suppress certain signals.
The RF pulses are represented digitally as complex numbers. Sequence controller 10 supplies these numbers in real and imaginary parts to digital-analog converters 14a-14b in RF system 11 to create corresponding analog pulse sequences. Transmission channel 15 modulates the pulse sequences with a radio-frequency carrier signal having a base frequency corresponding to the resonance frequency of the nuclear spins in the volume to be imaged.
RF coil 7 both emits radio-frequency pulses as described above and scans the alternating field which is produced as a result of processing nuclear spins, i.e., the nuclear spin echo signals. The received signals are received by multiplexer 13, amplified by RF amplifier 16 and demodulated in receiving channel 17 of RF system 11 in a phase-sensitive manner. Analog-digital converters 18a and 18b convert the demodulated signals into digitized real and imaginary components.
Electrocardiograph (“ECG”) monitor 19 acquires ECG signals from electrodes placed on patient 4. Such physiological signals may be used by sequence controller 10 to synchronize, or “gate”, transmitted RF pulses of a spectroscopy pulse sequence based on the heartbeat of patient 4 as described herein.
Computing system 30 receives the digitized real and imaginary components from analog-digital converters 18a and 18b and may process the components according to known techniques. Such processing may, for example, include reconstructing two-dimensional or three-dimensional images by performing a Fourier transformation of raw k-space data, performing other image reconstruction techniques such as iterative or back-projection reconstruction techniques, applying filters to raw k-space data or to reconstructed images, generating functional magnetic resonance images, calculating motion or flow images, and generating a chemical shift vs. magnitude spectrum.
System 30 may comprise any general-purpose or dedicated computing system. Accordingly, system 30 includes one or more processing units 31 (e.g., processors, processor cores, execution threads, etc.) configured to execute processor-executable program code to cause system 30 to operate as described herein, and storage device 32 for storing the program code. Storage device 32 may comprise one or more fixed disks, solid-state random access memory, and/or removable media (e.g., a thumb drive) mounted in a corresponding interface (e.g., a USB port).
One or more processing units 31 may execute program code of control program 33 to provide instructions to sequence controller 10 via MR system interface 34. For example, sequence controller 10 may be instructed to initiate a desired pulse sequence of pulse sequences 37. In particular, sequence controller 10 may be instructed to control the switching of magnetic field gradients via amplifiers 8a-8c at appropriate times, the transmission of radio-frequency pulses having a specified phase and amplitude at specified times via RF system 11 and RF amplifier 12, and the readout of the resulting MR signals. The timing of the various pulses of a pulse sequence may be based on physiological data received by ECG monitor interface 36.
Storage device 32 stores MR images 38 generated as described herein. Such images may be provided to terminal 40 via terminal interface 35 of system 30. Terminal interface 35 may also receive input from terminal 40, which may be used to provide commands to control program 33 in order to initiate single-voxel spectroscopy as described herein. Terminal 40 may comprise a display device and an input device coupled to system 30. In some embodiments, terminal 40 is a separate computing device such as, but not limited to, a desktop computer, a laptop computer, a tablet computer, and a smartphone.
Each element of system 1 may include other elements which are necessary for the operation thereof, as well as additional elements for providing functions other than those described herein. Storage device 32 may also store data and other program code for providing additional functionality and/or which are necessary for operation of system 30, such as device drivers, operating system files, etc.
Executable program code according to the above description may be stored on a form of non-transitory computer-readable media. Computer-readable media includes volatile and nonvolatile, removable, and non-removable media implemented in any method or technology for storage of information such as program code, data structures, program modules or other data. Computer-readable media includes, but is not limited to, random access memory (RAM), read-only memory (ROM), electrically erasable programmable ROM (EEPROM), flash memory or other memory technology, compact disk ROM (CD-ROM), digital volatile disks (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired instructions and which may be accessed by a system (e.g., a computer), including by internet or other computer network form of access.
The foregoing diagrams represent logical architectures for describing processes according to some embodiments, and actual implementations may include more or different components arranged in other manners. Other topologies may be used in conjunction with other embodiments. Moreover, each component or device described herein may be implemented by any number of devices in communication via any number of other public and/or private networks. Two or more of such computing devices may be located remote from one another and may communicate with one another via any known manner of network(s) and/or a dedicated connection. Each component or device may comprise any number of hardware and/or software elements suitable to provide the functions described herein as well as any other functions. For example, any computing device used in an implementation of a system according to some embodiments may include a processor to execute program code such that the computing device operates as described herein.
Embodiments described herein are solely for the purpose of illustration. Those in the art will recognize other embodiments may be practiced with modifications and alterations to that described above.
This application claims priority to U.S. Provisional Patent Application No. 63/271,832, filed Oct. 26, 2021, the disclosure of which is incorporated herein by reference for all purposes.
Number | Name | Date | Kind |
---|---|---|---|
10591568 | Rehwald et al. | Mar 2020 | B2 |
20050245809 | Wolff | Nov 2005 | A1 |
20100219829 | Rehwald | Sep 2010 | A1 |
20110074417 | Kitane | Mar 2011 | A1 |
20160132746 | Saranathan | May 2016 | A1 |
20170364252 | Deshpande | Dec 2017 | A1 |
Number | Date | Country |
---|---|---|
102005027963 | Dec 2006 | DE |
2008068089 | Mar 2008 | JP |
Entry |
---|
Kim, Raymond J., et al. “Relationship of MRI delayed contrast enhancement to irreversible injury, infarct age, and contractile function.” Circulation 100.19 (1999): 1992-2002. |
Simonetti, Orlando P., et al. “An improved MR imaging technique for the visualization of myocardial infarction.” Radiology 218.1 (2001): 215-223. |
Kellman, Peter, et al. “Motion-corrected free-breathing delayed enhancement imaging of myocardial infarction.” Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine 53.1 (2005): 194-200. |
Serai, Suraj D., et al. “Newly developed methods for reducing motion artifacts in pediatric abdominal MRI: tips and pearls.” American Journal of Roentgenology 214.5 (2020): 1042-1053. |
Zaitsev, Maxim, Julian Maclaren, and Michael Herbst. “Motion artifacts in MRI: A complex problem with many partial solutions.” Journal of Magnetic Resonance Imaging 42.4 (2015): 887-901. |
Wildgruber, M., et al. “Inversion-recovery single-shot cardiac MRI for the assessment of myocardial infarction at 1.5 T with a dedicated cardiac coil.” The British journal of radiology 85.1017 (2012): e709-e715. |
Xue, Hui, et al. “Motion correction for myocardial T1 mapping using image registration with synthetic image estimation.” Magnetic resonance in medicine 67.6 (2012): 1644-1655. |
Kellman, Peter, et al. “Phase-sensitive inversion recovery for detecting myocardial infarction using gadolinium-delayed hyperenhancement.” Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine 47.2 (2002): 372-383. |
Sievers, Burkhard, et al. “Rapid detection of myocardial infarction by subsecond, free-breathing delayed contrast-enhancement cardiovascular magnetic resonance.” Circulation 115.2 (2007): 236-244. |
Huber, Armin, et al. “Single-shot inversion recovery TrueFISP for assessment of myocardial infarction.” American Journal of Roentgenology 186.3 (2006): 627-633. |
Number | Date | Country | |
---|---|---|---|
20230132314 A1 | Apr 2023 | US |
Number | Date | Country | |
---|---|---|---|
63271832 | Oct 2021 | US |