Embodiments of the present invention relate generally to magnetic resonance imaging (“MRI”) systems, and more specifically, to a system and method for performing wave-encoded magnetic resonance imaging of an object.
MRI is a widely accepted and commercially available technique for obtaining digitized visual images representing the internal structure of objects having substantial populations of atomic nuclei that are susceptible to nuclear magnetic resonance (“NMR”). Many MRI systems use superconductive magnets to scan a subject/patient via imposing a strong main magnetic field on the nuclei in the subject to be imaged. The nuclei are excited by a radio frequency (“RF”) signal/pulse transmitted by a RF coil at characteristics NMR (Larmor) frequencies. By spatially disturbing localized magnetic fields surrounding the subject and analyzing the resulting RF responses from the nuclei as the excited protons relax back to their lower energy normal state, a map or image of these nuclei responses as a function of their spatial location is generated and displayed. An image of the nuclei responses provides a non-invasive view of a subject's internal structure.
Wave-encoded MRI is a form of MRI in which oscillating magnetic gradients are applied during the readout of each k-space line on the phase and slice-encoding axes. The resulting sampling strategy spreads the aliasing in all spatial directions and allows highly accelerated imaging with minimal residual artifacts due to efficient use of the three-dimensional (“3D”) distribution of RF coil sensitivities.
Many traditional MR acquisition techniques for wave-encoded MRI, hereinafter also referred to simply as “wave-encoded MR acquisition techniques”, such as Wave-CAIPI and Wave-CS, however, often require a separate calibration scan to estimate a wave-encoded point-spread function (“PSF”). The addition of a separate calibration scan, however, increases the complexity and total time of a wave-encoded MRI scan, hereinafter also referred to simply as “MRI scan” and/or “scan”. While some wave-encoded MR acquisition techniques utilize theoretically estimated wave PSF, such wave-encoded MR acquisition techniques often produce image artifacts, e.g., distortion and blurring. Additionally, while some MRI systems utilize auto-calibrated wave-CAIPI reconstruction for estimating wave PSF with CAPIRINHA-type uniform under-sampling, such MRI systems typically have excessively long time requirements, to the point of being impractical for many sampling patterns, e.g., variable-density and pseudo-random sampling.
Further, many traditional wave-encoded MR acquisition techniques usually require additional prephasers and rephasers to refocus the MR signal to central k-space. The use of prephasers and rephasers, however, often increases the total scan time and/or echo spacing requirements for fast spin echo based MR acquisitions.
Further still, many traditional wave-encoded MR acquisition techniques require an additional calibration scan for estimating RF coil sensitivity maps, which adds time and complexity to an MRI scan.
Yet further still, the k-space modulation introduced by wave gradients often makes applying known motion-correction methods to traditional wave-encoded MR acquisition techniques difficult.
What is needed, therefore, is an improved system and method for performing wave-encoded magnetic resonance imaging of an object.
In an embodiment, a method for performing wave-encoded magnetic resonance imaging of an object is provided. The method includes applying one or more wave-encoded magnetic gradients to the object, and acquiring MR signals from the object. The method further includes calibrating a wave point-spread function and reconstructing an image from the MR signals based at least in part on the calibrated wave point-spread function. Calibration of the wave point-spread function is based at least in part on one or more intermediate images generated from the MR signals.
In another embodiment, an MRI system for performing wave-encoded magnetic resonance imaging of an object is provided. The MRI system includes an MRI controller operative to: apply one or more wave-encoded magnetic gradients to the object; acquire MR signals from the object; calibrate a wave point-spread function; and reconstruct an image from the MR signals based at least in part on the calibrated wave point-spread function. The MRI controller calibrates the wave point-spread function based at least in part on one or more intermediate images generated from the MR signals.
In yet another embodiment, a non-transitory computer readable medium storing instructions is provided. The stored instructions are configured to adapt an MRI controller to: apply one or more wave-encoded magnetic gradients to an object; acquire MR signals from the object; calibrate a wave point-spread function; and reconstruct an image from the MR signals based at least in part on the calibrated wave point-spread function. Calibration of the wave point-spread function is based at least in part on one or more intermediate images generated from the MR signals.
The present invention will be better understood from reading the following description of non-limiting embodiments, with reference to the attached drawings, wherein below:
Reference will be made below in detail to exemplary embodiments of the invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference characters used throughout the drawings refer to the same or like parts, without duplicative description.
As used herein, the terms “substantially,” “generally,” and “about” indicate conditions within reasonably achievable manufacturing and assembly tolerances, relative to ideal desired conditions suitable for achieving the functional purpose of a component or assembly. As used herein, “electrically coupled,” “electrically connected,” and “electrical communication” mean that the referenced elements are directly or indirectly connected such that an electrical current may flow from one to the other. The connection may include a direct conductive connection, i.e., without an intervening capacitive, inductive or active element, an inductive connection, a capacitive connection, and/or any other suitable electrical connection. Intervening components may be present.
Further, while the embodiments disclosed herein are described with respect to an MRI system, it is to be understood that embodiments of the present invention may be applicable to other imaging systems which utilize wave-encoded magnetic gradients. Further still, as will be appreciated, embodiments of the present invention related imaging systems may be used to analyze tissue generally and are not limited to human tissue.
Referring now to
The MRI system control 32 includes a set of modules connected together by a backplane 38. These include a CPU module 40 and a pulse generator module 42, which connects to the operator console 12 through a serial link 44. It is through link 44 that the system control 32 receives commands from the operator to indicate the scan sequence that is to be performed. The pulse generator module 42 operates the system components to execute the desired scan sequence and produces data which indicates the timing, strength and shape of the RF pulses produced, and the timing and length of the data acquisition window. The pulse generator module 42 connects to a set of gradient amplifiers 46, to indicate the timing and shape of the gradient pulses that are produced during the scan. The pulse generator module 42 can also receive patient data from a physiological acquisition controller 48 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes attached to the patient. And finally, the pulse generator module 42 connects to a scan room interface circuit 50, which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 50 that a patient positioning system 52 receives commands to move the patient to the desired position for the scan.
The pulse generator module 42 operates the gradient amplifiers 46 to achieve desired timing and shape of the gradient pulses that are produced during the scan. The gradient waveforms produced by the pulse generator module 42 are applied to the gradient amplifier system 46 having Gx, Gy, and Gz amplifiers. Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly, generally designated 54, to produce the magnetic field gradients used for spatially encoding acquired signals. The gradient coil assembly 54 forms part of a magnet assembly 56, which also includes a polarizing magnet 58 (which in operation, provides a homogenous longitudinal magnetic field B0 throughout a target volume 60 that is enclosed by the magnet assembly 56) and a whole-body (transmit and receive) RF coil 62 (which, in operation, provides a transverse magnetic field B1 that is generally perpendicular to B0 throughout the target volume 60). As will be appreciated, however, embodiments of the invention may use head coils and/or other localized RF coils.
The resulting signals emitted by the excited nuclei in the patient may be sensed by the same RF coil 62 and coupled through the transmit/receive switch 64 to a preamplifier 66. The amplifier MR signals are demodulated, filtered, and digitized in the receiver section of a transceiver 68. The transmit/receive switch 64 is controlled by a signal from the pulse generator module 42 to electrically connect an RF amplifier 70 to the RF coil 62 during the transmit mode and to connect the preamplifier 66 to the RF coil 62 during the receive mode. The transmit/receive switch 64 can also enable a separate RF coil (for example, a surface coil) to be used in either transmit or receive mode.
The MR signals picked up by the RF coil 62 are digitized by the transceiver module 68 and transferred to a memory module 72 in the system control 32. A scan is complete when an array of raw k-space 74 (
As illustrated in
Turning now to
Moving now to
As will be appreciated, in certain aspects of the invention, the wave-encoded magnetic gradients 90 may be self-refocusing with respect to a center 104 (also shown in
Additionally, the wave-encoded magnetic gradients 90 may also be configured to be second-order differential between neighboring/adjacent lobes, e.g., lobes 120 and 123, and lobes 122 and 125, which may mitigate and/or avoid rapid changes in the gradient slew rate. While
Turning to
Referring now to
PSF[kx,y,z]=exp(−iγ∫0t[k
where gy and gz are the wave magnetic gradients in the phase- and slice-encoding directions, respectively, and y and z are the position in image space. For 2D imaging applications, the above equation may be simplified to the following:
PSF[kx,y]=exp(−iγ∫0t[k
The image m can then be reconstructed 140 from the zero-filled wave-encoded k-space s using the theoretical PSF:
m=
x
−1·PSF·y−1z−1s
Next, as represented by steps 142, 144, 146, and/or 148, the following metric:
G
obj(m)=Σi,j(|gi,j(m)|/Σi,j|gi,j(m)|)2
where gi,j is the image gradient at coordinate i,j, may be maximized iteratively through an optimization method, e.g., the Nelder-Mead simplex method, to yield the calibrated wave PSF. As will be appreciated, other appropriate metrics may be used. Thus, the optimization problem can be formulated as:
argminvariables{−Gobj(x−1·PSF(variables)·y−1z−1s)}
where the optimization variables are the amplitude and time delay of the wave-encoded magnetic gradients 90, and the center of the prescribed field of view (“FOV”) in the phase-encoding direction.
Moving to
S
wave
[k
x
,y,z]=PSF[kx,y,z]·SCartesian(kx,y,z)
which can be expressed as:
where C1 and C2 are constants associated with the gy and gz gradients, y and z are positions in image space, Δy and Δz are the corresponding spatial resolutions, and iy and iz are spatial indices. Accordingly, the wave-encoded central k-space 74 can be treated as an independent low-resolution k-space modulated by the same calibrated wave PSF. Thus, by replacing Δy and Δz with
respectively, where Ny and Nz define the size of the acquisition matrix, and Ny,calibration and Nz,calibration represent the size of the calibration region, i.e., the size of the sampled central k-space 74, a low resolution wave PSF 152 can be obtained, which in turn can be demodulated from the wave-encoded MR signals to subsequently reconstruct the Cartesian central k-space using an inverse Fourier transformation. As will be appreciated, the resulting Cartesian k-space center can be used to estimate/calibrate the coil sensitivity map. Accordingly, the k-space 74 in
Thus, in embodiments, reconstruction 96 (
argminm|Dyz·PSF·xEm−s∥22+λ∥Ψm∥1
where D is a k-space sampling operator, PSF is the wave PSF, E contains the RF coil 62 sensitivities, m is the desired image, s is the under-sampled wave-encoded k-space, W is a wavelet transform operator such that the term ∥Ψm∥1 promotes sparsity of the solution, and A is the regularization parameter related to the L1-norm of the wavelet coefficients m.
As stated above, embodiments of the method 86 (
S
m
[n]=S
0·exp(i·2π·(kx[n]·dx+ky[n]·dy+kz[n]·dz))·exp(−C1·dy−Cz·dz)
where dx, dy, and dz are the estimated distances of linear translation in three orthogonal directions and n denotes a particular sampling time point. As shown below, starting from a reference image m0[x, y, z], the previous relation can be derived according to Fourier theory and the expression of the wave-encoding PSF as follows:
Therefore, any motion-correction technique that uses linear translation estimates to correct motion corrupted MR signals can be applied to wave-encoded data by simply demodulating the following quantity: exp (i2π(kxdx+kydy+kzdz))·exp(−C1dy−C2 dz). As will be appreciated, the exp (−C1dy−C2 dz) term compensates for phase shifts due to linear translation of the wave-encoding PSF. Accordingly, other motion-correction methods may also be applied with and/or in place of the above described linear phase correction. For example, embodiments of the invention may apply advanced motion-correction methods based on higher-order phase correction, e.g., auto-focusing, or methods that do not rely on phase correction, e.g., soft-gating and/or locally low rank.
Finally, it is also to be understood that the system 10 may include the necessary electronics, software, memory, storage, databases, firmware, logic/state machines, microprocessors, communication links, displays or other visual or audio user interfaces, printing devices, and any other input/output interfaces to perform the functions described herein and/or to achieve the results described herein. For example, as previously mentioned, the system may include at least one processor and system memory/data storage structures, which may include random access memory (RAM) and read-only memory (ROM). The at least one processor of the system 10 may include one or more conventional microprocessors and one or more supplementary co-processors such as math co-processors or the like. The data storage structures discussed herein may include an appropriate combination of magnetic, optical and/or semiconductor memory, and may include, for example, RAM, ROM, flash drive, an optical disc such as a compact disc and/or a hard disk or drive.
Additionally, a software application that adapts the controller to perform the methods disclosed herein may be read into a main memory of the at least one processor from a computer-readable medium. The term “computer-readable medium,” as used herein, refers to any medium that provides or participates in providing instructions to the at least one processor of the system 10 (or any other processor of a device described herein) for execution. Such a medium may take many forms, including but not limited to, non-volatile media and volatile media. Non-volatile media include, for example, optical, magnetic, or opto-magnetic disks, such as memory. Volatile media include dynamic random access memory (DRAM), which typically constitutes the main memory. Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, a RAM, a PROM, an EPROM or EEPROM (electronically erasable programmable read-only memory), a FLASH-EEPROM, any other memory chip or cartridge, or any other medium from which a computer can read.
While in embodiments, the execution of sequences of instructions in the software application causes at least one processor to perform the methods/processes described herein, hard-wired circuitry may be used in place of, or in combination with, software instructions for implementation of the methods/processes of the present invention. Therefore, embodiments of the present invention are not limited to any specific combination of hardware and/or software.
It is further to be understood that the above description is intended to be illustrative, and not restrictive. For example, the above-described embodiments (and/or aspects thereof) may be used in combination with each other. Additionally, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from its scope.
For example, in an embodiment, a method for performing wave-encoded magnetic resonance imaging of an object is provided. The method includes applying one or more wave-encoded magnetic gradients to the object, and acquiring MR signals from the object. The method further includes calibrating a wave point-spread function and reconstructing an image from the MR signals based at least in part on the calibrated wave point-spread function. Calibration of the wave point-spread function is based at least in part on one or more intermediate images generated from the MR signals. In certain embodiments, the one or more wave-encoded magnetic gradients are self-refocusing with respect to a center of a corresponding k-space. In certain embodiments, the one or more wave-encoded magnetic gradients are symmetrical. In certain embodiments, the one or more wave-encoded magnetic gradients are at least one of sinusoidal, triangular, square, and trapezoidal. In certain embodiments, a first integral of a first lobe and a second integral of a final lobe of at least one of the one or more wave-encoded magnetic gradients are less than or equal to about half of a third integral of a middle lobe of the same wave-encoded magnetic gradient. In certain embodiments, the one or more wave-encoded magnetic gradients are second-order differentiable between neighboring lobes. In certain embodiments, acquiring MR signals from the object via the MRI system includes receiving the wave-encoded MR signals via two or more RF coils. In such embodiments, the method further includes calibrating a coil sensitivity map of the two or more RF coils based at least in part on the MR signals and the calibrated wave point-spread function, and reconstructing an image from the MR signals is further based at least in part on the calibrated coil sensitivity map. In certain embodiments, calibrating a coil sensitivity map based at least in part on the MR signals and the calibrated wave point-spread function includes generating a low resolution wave point-spread function from the calibrated wave point-spread function. In such embodiments, the coil sensitivity map is further based at least in part on the low resolution wave point-spread function and central k-space calibration data within the acquired MR signals. In certain embodiments, the method further includes motion-correcting the acquired MR signals. In certain embodiments, motion-correcting the acquired MR signals is based at least in part on localized translation motion estimates.
Other embodiments provide for an MRI system for performing wave-encoded magnetic resonance imaging of an object. The MRI system includes an MRI controller operative to: apply one or more wave-encoded magnetic gradients to the object; acquire MR signals from the object; calibrate a wave point-spread function; and reconstruct an image from the MR signals based at least in part on the calibrated wave point-spread function. The MRI controller calibrates the wave point-spread function based at least in part on one or more intermediate images generated from the MR signals. In certain embodiments, the one or more wave-encoded magnetic gradients are self-refocusing with respect to a center of a corresponding k-space. In certain embodiments, the one or more wave-encoded magnetic gradients are symmetrical. In certain embodiments, the one or more wave-encoded magnetic gradients are at least one of sinusoidal, triangular, square, and trapezoidal. In certain embodiments, a first integral of a first lobe and a second integral of a final lobe of at least one of the one or more wave-encoded magnetic gradients are less than or equal to about half of a third integral of a middle lobe of the same wave-encoded magnetic gradient. In certain embodiments, the one or more wave-encoded magnetic gradients are second-order differentiable between neighboring lobes. In certain embodiments, the MRI system acquires the MR signals via two or more RF coils. In such embodiments, the MRI controller is further operative to calibrate a coil sensitivity map of the two or more RF coils based at least in part on the MR signals and the calibrated wave point-spread function. In certain embodiments, reconstruction of the image via the MRI controller is further based at least in part on the calibrated coil sensitivity map. In certain embodiments, the MRI controller is further operative to generate a low resolution wave point-spread function from the calibrated wave point-spread function. In such embodiments, calibration of the coil sensitivity map via the MRI controller is further based at least in part on the low resolution wave point-spread function and central k-space calibration data within the acquired MR signals. In certain embodiments, the MRI controller is further operative to motion correct the acquired MR signals.
Yet still other embodiments provide for a non-transitory computer readable medium storing instructions. The stored instructions are configured to adapt an MRI controller to: apply one or more wave-encoded magnetic gradients to an object; acquire MR signals from the object; calibrate a wave point-spread function; and reconstruct an image from the MR signals based at least in part on the calibrated wave point-spread function. Calibration of the wave point-spread function is based at least in part on one or more intermediate images generated from the MR signals.
Accordingly, by incorporating self-refocusing wave-encoded magnetic gradient waveforms, some embodiments of the present invention may provide for an Wave-encoded MR acquisition technique that does not require prephasers and rephasers. By incorporating an auto-calibrating technique for estimating wave PSF, some embodiments of the present invention may eliminate the need for external calibrations. As will be appreciated, by excluding prephasers and rephasers and/or external calibrations, some embodiments of the present invention may have greatly improved scan efficiencies over traditional wave-encoded MR acquisition techniques. Such embodiments may also provide for image artifacts due to imperfections in the magnetic gradients to be modeled and corrected simultaneously.
Further, by utilizing a data-driven coil sensitivity map calibration, some embodiments of the present invention may provide for significant improvements in the accuracy of the coil sensitivity map in the presence of motion, which in turn may further provide a more simplified workflow over traditional wave-encoded MR acquisition techniques.
Yet further still, by utilizing motion-correction based on localized translational motion estimates applied to wave-encoded k-space, some embodiments of the present invention may provide for the correction of artifacts due to intra-scan motion by the imaged subject/object.
Additionally, some embodiments of the present invention provide for a wave-encoded MR acquisition technique that incorporates parallel imaging having voxel spreading in the FE direction. Thus, such embodiments make better/more efficient use of RF coil sensitivities so as to achieve lower g-factors than traditional Cartesian acquisitions. As will be appreciated, such embodiments may also provide for higher acceleration factors than traditional Cartesian acquisitions without degrading image quality, and/or having a more uniform distribution of motion artifacts, which in turn may further provide for improved motion robustness over traditional Cartesian acquisitions.
Moreover, by utilizing self-refocusing wave-encoding magnetic gradients, some embodiments of the present invention provide for a wave-encoded MR acquisition technique having a readout duration identical to that of a traditional Cartesian acquisition while maintaining contrast within the generated images. As will also be appreciated, such embodiments provide for a wave-encoded MR acquisition that has the same k-space modulation due to T2 decay as conventional acquisitions for fast spin echo imaging.
Also, by calibrating the wave-encoded magnetic gradients in the image domain, some embodiments of the present invention provide for time improvements over other SENSE-based reconstruction methods, and can additionally be applied to variable-density sampling patterns. Thus, some embodiments of the present invention may provide for simultaneous correction of eddy currents and gradient delays by including their effect in the gradient waveform model.
In summary, some embodiments of the present invention provide for an efficient and motion-robust wave-encoded MR acquisition technique that performs highly accelerated wave-encoded MRI for various clinical applications without increasing the overall scan duration and with minimal impact on workflow.
Additionally, while the dimensions and types of materials described herein are intended to define the parameters of the invention, they are by no means limiting and are exemplary embodiments. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. In the appended claims, the terms “including” and “in which” are used as the plain-English equivalents of the respective terms “comprising” and “wherein.” Moreover, in the following claims, terms such as “first,” “second,” “third,” “upper,” “lower,” “bottom,” “top,” etc. are used merely as labels, and are not intended to impose numerical or positional requirements on their objects. Further, the limitations of the following claims are not written in means-plus-function format are not intended to be interpreted as such, unless and until such claim limitations expressly use the phrase “means for” followed by a statement of function void of further structure.
This written description uses examples to disclose several embodiments of the invention, including the best mode, and also to enable one of ordinary skill in the art to practice the embodiments of invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to one of ordinary skill in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.
As used herein, an element or step recited in the singular and proceeded with the word “a” or “an” should be understood as not excluding plural of said elements or steps, unless such exclusion is explicitly stated. Furthermore, references to “one embodiment” of the present invention are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. Moreover, unless explicitly stated to the contrary, embodiments “comprising,” “including,” or “having” an element or a plurality of elements having a particular property may include additional such elements not having that property.
Since certain changes may be made in the above-described invention, without departing from the spirit and scope of the invention herein involved, it is intended that all of the subject matter of the above description shown in the accompanying drawings shall be interpreted merely as examples illustrating the inventive concept herein and shall not be construed as limiting the invention.
This application claims the benefit of U.S. Provisional Patent Application No. 62/425,966 filed on Nov. 23, 2016, which is herein incorporated by reference in its entirety.
This invention was made with Government support under contracts EB009690, EB019241, and EB015891 awarded by the National Institutes of Health. The Government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
62425966 | Nov 2016 | US |