Patent Application
20160202336
The present disclosure generally relates to magnetic resonance spectroscopic imaging.
This disclosure relates generally to diagnostic imaging and, more particularly, to an apparatus and method of high resolution magnetic resonance spectroscopic imaging.
When a substance such as human tissue is subjected to a uniform magnetic field (such as a polarizing field B0), the individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but precess about it in random order at their characteristic Larmor frequency. If the substance, or tissue, is subject to a magnetic field (excitation field B1) in the x-y plane and which is near the Larmor frequency, the net aligned moment, or “longitudinal magnetization”, Mz, may be rotated or “tipped” into the x-y plane to produce a net transverse magnetic moment Mt. A signal is emitted by the excited spins after the excitation signal B1 is terminated and this signal may be received and processed to form an image.
When utilizing these signals to produce images, magnetic field gradients (Gx, Gy, and Gz) are employed. Typically, the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method used. The resulting set of received signals are digitized and processed to reconstruct the image using one of many known reconstruction techniques.
Magnetic resonance spectroscopic imaging (MRSI) is a noninvasive imaging method that provides spectral information in addition to the structural information that is generated by magnetic resonance imaging (MRI) alone. Traditional MRI generates a black-and-white image in which brightness is determined primarily by the water molecule concentrations, and the T1 and/or T2 relaxation times of the tissue being imaged; the spectral information obtained in MRSI provides additional information about metabolic activity. MRSI can be performed on a standard MRI scanner, and has broad applications in medicine, oncology, and general physiological studies. And, when hydrogen is the target element, MRSI is also called 1H-nuclear magnetic resonance spectroscopic imaging and proton magnetic resonance spectroscopic imaging. Similarly, we have 31P magnetic resonance spectroscopic imaging and 13C magnetic resonance spectroscopic imaging.
MRSI has been recognized as a powerful tool for noninvasive metabolic studies, but clinical and research applications of this technology have been developing more slowly than expected. Reasons for the slow development of the technology include but are not limited to long data acquisition time, poor spatial resolution, and low signal-to-noise ratio (SNR), as examples.
Significant efforts have been made to address the above issues, resulting in a large number of new data acquisition and reconstruction methods for spectroscopic imaging. For fast data acquisition, one approach is to incorporate echo-planar-type of data acquisition schemes with spectroscopic imaging. Many methods (and pulse sequences) have been proposed to implement this data acquisition strategy. And, although echo-planar spectroscopic imaging (EPSI) methods can significantly reduce the data acquisition time for spatiospectral encoding, it is at the expense of SNR. Another approach to accelerated spectroscopic imaging is to use parallel imaging in which parallel data acquisition occurs using phased array coils.
Advanced MRSI reconstruction has focused on using prior information to compensate for the lack of sufficient measurements or SNR. To this end, a number of reconstruction models have been proposed, but reconstruction methods alone may not provide adequate levels of improvements in spatial resolution, data acquisition speed, and SNR needed to have a major impact on in vivo MR spectroscopic imaging.
As such, there is a need for improved MR spectroscopic imaging.
Embodiments are directed toward a method and apparatus to enable high-speed, high-resolution MR spectroscopic imaging.
According to one aspect, a device to acquire magnetic resonance spectroscopic imaging (MRSI) data includes a plurality of gradient coils positioned about a bore of a magnet, an RF transceiver controlled by a pulse module to transmit RF signals to an RF coil assembly and receive imaging data, and a computer. The computer can be programmed to control the data acquisition system of an MRI scanner to acquire MRSI data with special sampling patterns in (k,t)-space (corresponding to special spatiospectral encoding strategies). One unique feature of the disclosed sampling patterns is variable sampling densities in (k,t)-space to enable high-speed MRSI data acquisition with sparse sampling in k-space and/or time; another unique feature of the disclosed sampling patterns is that it enables the use of subspace (or low-rank matrix/tensor) models to recover high-quality, high-resolution spatiospectral functions from sparsely sampled (k,t)-space MRSI data. The MRSI data collected using such sampling schemes can be rearranged in multiple ways. One way is to organize the acquired MRSI data into two data sets: one data set (called “training” data) with relatively high SNR, and high temporal resolution (often covering central k-space region), and another data set (called “imaging” data) with extended k-space coverage to achieve high spatial resolution. The “imaging” data can sample k-space and/or time sparsely with or without signal averaging to achieve high-speed data acquisition.
According to another aspect, a method of magnetic resonance spectroscopic imaging (MRSI) that builds on subspace (or low-rank matrix/tensor) models and uses the models for data acquisition, processing and image reconstruction, wherein a corresponding (k,t)-space is sparsely sampled in variable density and variable SNR that enable high-speed data acquisition while providing sufficient data for determination of subspace structures and for recovery of high-quality, high-resolution spatiospectral functions using low-rank matrix/tensor models. Various other features and advantages will be made apparent from the following detailed description and the drawings.
An apparatus and method is provided that acquires multiple 2D or 3D MRSI data sets using special spatiospectral encoding schemes that cover (k,t)-space sparsely. These data enable the recovery of high-quality, high-resolution spatiospectral functions using subspace (or low-rank matrix/tensor) models. The spatiospectral information thus obtained provides biochemical information of the object being imaged.
The following definitions can be utilized by one or more of the exemplary embodiments.
A device and/or method is provided for acquiring spatiospectral distributions from an object. A magnetic field generator can generate a strong static magnetic field where the object is arranged in. A plurality of gradient coils can be positioned about a bore of a magnet. An RF transceiver can be controlled by a pulse module to transmit RF signals to an RF coil assembly and to receive imaging data. A computer can control the data acquisition system to acquire sparsely sampled (k,t)-space data and can reconstruct spatiospectral functions from the sparsely sampled (k,t)-space data.
One or more of the exemplary embodiments can utilize modeling of the spatiospectral function of interest using a Partial Separability (PS) model as
where x1, x2, . . . x{circumflex over (d)} are {circumflex over (d)} separable groups of the complete set of variables x1, x2, . . . , xd for the image domain, such that xi=(xi
For {circumflex over (d)}>2, the model in Eq. (I) can be generalized as
which is equivalent to the form in Eq. (I) for {circumflex over (d)}≦2. Note that both Eq. (I) and Eq. (II) represent a subspace model that is the basis of the data acquisition and reconstruction schemes associated with the exemplary embodiments.
The model in Eqs. (I) and (II) can be expressed in different forms for magnetic resonance spectroscopic imaging (MRSI) in different contexts, specifically
ρ(r,t)=ρ(x,y,t)=Σl=1Lcl(r)ψl(t) or ρ(r,f)=ρ(x,y,f)=Σl=1Lcl(r)ψl(f),
Partial separability of the underlying spatiospectral function implies that the tensor will be rank-deficient. In other words, although the tensor, in general, is in a very high-dimensional space; those tensors representing meaningful spatiospectral functions reside in a very low-dimensional subspace. This property is effectively utilized in the disclosed method to achieve high-speed, high-resolution MRSI.
The data generation process in general magnetic resonance imaging experiments can be modeled as
where kx
a)s(k,t)=s(kx,ky,t)=∫∫ρ(x,y,t)e−i2π(k
s(k,t)=s(kx,ky,t)=∫∫(∫ρ(x,y,f)e−i2πftdf)e−i2π(k
b)s(k,t)=s(kx,ky,kx,t)=∫∫∫ρ(x,y,z,t)e−i2π(k
s(k,t)=s(kx,ky,kx,t)=∫∫∫ρ(x,y,z,f)e−i2πftdf)e−i2π(k
Referring to
System control 128 includes a set of modules connected together by a backplane (not shown). These include a CPU module 132 and a pulse generator module 134 which connects to the operator console 102 through a serial link 136. It can be through link 136 that the system control 128 receives commands from the operator to indicate a scan sequence that is to be performed. The pulse generator module 134 operates the system components to carry out 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 134 connects to a set of gradient amplifiers 138, to indicate the timing and shape of the gradient pulses that are produced during the scan. The pulse generator module 134 can also receive patient data from a physiological acquisition controller 140 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes attached to the patient. The pulse generator module 134 connects to a scan room interface circuit 142 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 142 that a patient positioning system 144 receives commands to move the patient to the desired position for the scan.
The gradient waveforms produced by the pulse generator module 134 are applied to the gradient amplifier system 138 having Gx, Gy, and Gz amplifiers. Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly generally designated 144 to produce the magnetic field gradients used for spatially encoding acquired signals. The gradient coil assembly 144 forms part of a magnet assembly 146 which includes a polarizing magnet 148 and a whole-body RF coil 150. In one embodiment, RF coil 150 can be a multi-channel coil. A transceiver module 152 in the system control 128 produces pulses which are amplified by an RF amplifier 154 and coupled to RF coil 150 by a transmit/receive switch 156. The resulting signals emitted by the excited nuclei in the patient may be sensed by the same RF coil 150 and coupled through the transmit/receive switch 156 to a preamplifier 158. The amplified MR signals are demodulated, filtered, and digitized in the receiver section of the transceiver 152. The transmit/receive switch 156 can be controlled by a signal from the pulse generator module 134 to electrically connect the RF amplifier 154 to the RF coil 150 during the transmit mode and to connect the preamplifier 158 to the RF coil 150 during the receive mode. The transmit/receive switch 156 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 multi-channel RF coil 150 are digitized by the transceiver module 152 and transferred to a memory module 160 in the system control 128. A scan can be complete when an array of raw k-space data has been acquired in memory module 160. This raw k-space data can be rearranged into separate k-space data arrays for each image to be reconstructed, and each of these can be input to an array processor 162 which operates to Fourier transform the data into an array of image data. This image data can be conveyed through serial link 130 to computer system 120 where it can be stored in memory 126. In response to commands received from operator console 102, image data may be archived in long term storage or it may be further processed by image processor 122 and conveyed to operator console 102 and presented on display 106.
Magnetic resonance spectroscopy (MRS) is used to measure the levels and changes of different metabolites in body tissues. The MR signal produces a spectrum of resonances that corresponds to different molecular arrangements of the isotope being “excited”. This signature is used to diagnose certain metabolic disorders, such as those affecting the brain, and to provide information on tumor metabolism.
Magnetic resonance spectroscopic imaging (MRSI) combines both spectroscopic and imaging methods to produce spatially localized spectra from within the sample or patient. The spatial resolution may be lower (limited by the available SNR), but the spectra in each voxel contains information about many metabolites. Because the available signal is used to encode spatial and spectral information, MRSI is typically deemed as a slow imaging modality and prefers higher field strengths, such as 3T and above, to achieve desirable SNR.
The disclosed method accelerates spectroscopic imaging using sparse sampling of (k,t)-space and subspace (or low-rank) modeling to enable high-resolution metabolic imaging with good SNR by exploiting a unique property known as partial separability (PS) of spectroscopic signals. This property indicates that high-dimensional spectroscopic signals reside in a very low-dimensional subspace and enables special data acquisition and image reconstruction strategies to be used to obtain high-resolution spatiospectral distributions from the imaging object with good SNR. More specifically, a variable density sampling pattern is disclosed for sparse sampling of (k,t)-space, which can be implemented in practical experiments using a 2D or 3D hybrid chemical shift imaging (CSI) and echo-planar spectroscopic imaging (EPSI) pulse sequence or a variable density 2D or 3D EPSI sequence; and a PS/low-rank model-based algorithm is disclosed for subspace estimation and image reconstruction from sparse data with the capability to incorporate prior information and field inhomogeneity correction. In one example, an EPSI pulse sequence is used having a range of readout bandwidth that allows for a customized EPSI sequence such that an echo spacing does not satisfy a spectral Nyquist criterion. In another example, a set of navigator data is included for correcting bipolar effects from the EPSI readouts.
In practical 2D/3D spectroscopic imaging experiments, the measured data can be further modeled as (based on the data acquisition equations described above):
where ρ(r,f) is the desired 2D/3D spatiospectral function, Sc (r) the sensitivity profile of the cth acquisition coil (with c being the coil index), V the excited volume of interest, Ωf the spectral bandwidth, γ the gyromagnetic ratio, ξ(k,t) the measurement noise (often assumed to be white Gaussian), and ΔB(r) the B0 field inhomogeneity (referred to as field inhomogeneity hereafter).
The disclosed method exploits an important property of ρ(r,t), known as spatiotemporal partial separability (PS). More specifically, ρ(r,t) may be expressed as:
where {φl(t)}l=1L can be viewed as a set of temporal basis functions, {cl(r)}l=1L are the corresponding spatial coefficients (which can also be viewed as spatial basis), and L is the model order (also called the order of separability). In spectroscopic imaging, the PS model may be used because the object imaged has a finite (L) number of tissue types, each of which has a distinct spectral signatures. It has been shown that PS implies matrix or tensor low-rankness. For example, a Casorati matrix formed from ρ(r,t) over any point set {(rp,tq)}p,q=1P,Q:
has a rank upper-bounded by L. This implies that ρ(r,t) resides in a low-dimensional space (or has a small number of degrees of freedom). More specifically, assume that ρ(r,t) can be represented by {ρ(rm,tn)}m,n=1M,N in a conventional pixel representation. Then, when viewed as a vector (formed from p(rm,tn)), it resides in a M N-dimensional space. However, the PS property implies that {ρ(rm,tn)}m,n=1M,N actually resides in a very low-dimensional subspace (spanned by {φl(t)}l=1L, for example). This property enables special data acquisition and image reconstruction strategies to be used to obtain high-resolution spatiospectral distributions with good SNR.
Data Acquisition: Extended (k,t) Coverage with Variable Density Sparse Sampling
To exploit the PS property for high-speed, high-resolution spectroscopic imaging, the disclosed method uses a special data acquisition (or spatiospectral encoding) strategy to achieve extended (k,t) coverage with sparse sampling and variable density. While there are many practical ways to implement the disclosed PS-based data acquisition (spatiospectral encoding) strategy, one unique feature of the disclosed sampling patterns is variable sampling densities in (k,t)-space to enable high-speed MRSI data acquisition with sparse sampling in k-space and/or time; another unique feature of the disclosed sampling patterns is that it enables the use of subspace (or low-rank matrix or tensor) models to recover high-quality, high-resolution spatiospectral functions from sparsely sampled (k,t)-space MRSI data. The MRSI data collected using such sampling schemes can be rearranged in multiple ways. One way is to organize the acquired MRSI data into two data sets: one data set (called “training” data, denoted as D1) with relatively high SNR, and high temporal resolution (often covering central k-space region), and another data set (called “imaging” data, denoted as D2) with extended k-space coverage to achieve high spatial resolution. The “imaging” data can sample k-space and/or time sparsely with or without signal averaging to achieve high-speed data acquisition. One exemplary practical implementation of the disclosed data acquisition strategy is illustrated in
Referring to
Chemical shift is the phenomenon that is seen when an isotope possessing a nuclear magnetic dipole moment resonates at a spectrum of resonance frequencies in a given magnetic field. These resonance frequencies, or chemical shifts, depend on the chemical environments of particular nuclei. Mapping the spatial distribution of nuclei associated with a particular chemical shift (e.g., hydrogen nuclei associated with water molecules or with lipid groups) is called chemical shift imaging.
Chemical shift imaging (CSI) is an extension of MR spectroscopy, allowing metabolite information to be measured in an extended region and to add the chemical analysis of body tissues to the potential clinical utility of magnetic resonance. The spatial location is phase encoded and a spectrum is recorded at each phase encoding step to allow the spectral acquisition in a number of volumes covering the whole sample. CSI provides mapping of chemical shifts, analog to individual spectral lines or groups of lines. Spatial resolution can be in one, two or three dimensions, but with long acquisition times of full 3D CSI. Commonly a slice-selected 2D acquisition is used. The chemical composition of each voxel is represented by spectra, or as an image in which the signal intensity depends on the concentration of an individual metabolite. Alternatively frequency-selective pulses excite only a single spectral component.
There are several known methods of performing CSI, e.g. the inversion recovery method, chemical shift selective imaging sequence, chemical shift insensitive slice selective RF pulse, the saturation method, spatial and chemical shift encoded excitation, and quantitative chemical shift imaging, as examples.
To reduce the relatively long acquisition time required for spatial encoding, fast spectroscopic imaging approaches were introduced. Among them was an echo planar-based encoding technique. In echo planar spectroscopic imaging (EPSI), an alternating current, which simultaneously encodes space and chemical shift (time), is applied along one of the spatial (readout) directions, thereby reducing the phase-encoding dimension by one and increasing speed. Recent advances in gradient hardware, which have increased both the speed of gradient switching and gradient strength achievable and improved eddy current performance of gradients, led to successful implementations of EPSI. EPSI allows for a larger coverage of the k-space, thus permitting a relatively high spatial resolution within practically acceptable acquisition time. When signal-to-noise ratio (SNR) is not a limiting factor, EPSI can be used to obtain spectroscopic images with resolutions comparable to those of anatomical images in practical scan time. Although for most spectroscopic applications the improved spatial resolution afforded by EPSI is not utilized due to SNR limitations, it can be used for imaging applications in which spectroscopic information may be beneficial.
The disclosed data acquisition scheme can be implemented using a hybrid CSI/EPSI pulse sequence, an example of which is shown in
One notable difference occurs, illustrated in
Thus,
Reconstruction from Sparse Data
Assume that D1={s1(km,tn)}m,n=1M
Determination of the temporal basis (or temporal subspace), {φ1(tn)}l=1L for n=1, 2, . . . , N1, from D1 is straightforward if the effect of field inhomogeneity ΔB(r) on s1(km,tn) is negligible. In practice, ΔB(r) can be significant, and its effect on s1(km, tn) should be reduced or removed. Assuming that ΔB(r) is available (which can be acquired during the MRSI experiment), field inhomogeneity correction on limited k-space spectroscopic data can be done using known methods, such as a penalized maximum-likelihood method, which enables the use of anatomical constraints for improving the correction performance with only limited k-space data. The corrected data, denoted as ŝ1(km,tn), is then used to form an M1×N1 Casorati matrix C C({ŝ(km,tn)}m,n=1M
After {φl(tn)} are determined, the spatial coefficients, {cl(rm)}l=1L, can be determined from the sparse data in D2 by solving a least-squares problem. For notation convenience, {ρ(rm,tn)}m,n=1M,N={Σl=1Lcl(rm)φl(tn)}m,n=1M,N is rewritten as CΦ (with CεM×Lφε
L×N such that Cml=cl(rm) and Φln=φl(tn), where M=M2 and N=N1 for the proposed data acquisition scheme), and use s2εCM
where the first term measures the data consistency of a reconstruction, and Ψ(•) is a regularization functional with regularization parameter λ. There are many choices for Ψ(•) to incorporate prior information about ρ(r,t) or ρ(r,f) (including both quadratic and sparsity-promoting penalties). Thus, the disclosed method includes demonstrating the concept and potential using Ψ(CΨ)=∥WDCΨ∥F2, where D is a finite difference operator and W contains edge weights derived from a high-resolution anatomical image. Integrating this regularization functional into Eq. (V) yields a weighted-l2 regularized least-squares problem that can be solved efficiently. As such, image reconstruction may be done by solving a regularized least-squares optimization problem.
Thus, referring to
Dataset D1 is corrected for B0 field inhomogeneity, at block 522, and the corrected field inhomogeneity data is generated at block 524. As described and as one example, such field inhomogeneity correction methods may include a penalized maximum-likelihood method, which enables the use of anatomical and/or low-rank constraints for improving the correction performance with only limited k-space data. At block 524, corrected field inhomogeneity data, denoted as ŝ1(km,tn), is generated. In one example, the field inhomogeneity effect is removed using a regularized reconstruction performed on a higher resolution grid.
At block 526 an spatiospectral function is reconstructed using the first and second datasets D1 and D2. That is, datasets D1 and D2 as modified (removal of water-fat of both datasets, and field inhomogeneity correction of D1), are used to reconstruct the spatiospectral (or spatiotemporal) function. Image reconstruction block 526 includes, as described above, a two-step procedure that includes a determination of the temporal basis at block 528, and a determination of the spatial coefficients at block 530. In one example, reconstruction or spatiospectral function generation includes (a) representing a desired spatiotemporal function using a low-rank matrix multiplication model (i.e., the spatial coefficient matrix multiplied with the temporal basis matrix), (b) estimating the spatial coefficients, and (c) forming the final spatiotemporal function by combining the spatial coefficients and the temporal bases. In one example, the spatial coefficients are estimated by solving a regularized least-squares optimization problem with either a quadratic or a non-quadratic regularization function incorporating anatomical information. The process ends at step 532.
In one example, a CSI pulse sequence is used to generate dataset D1, and an EPSI pulse sequence is used to generate dataset D2. However, in one implementation the CSI component used to acquire D1 in the disclosed sequence can be replaced by a traditional EPSI scheme if further improvement in data acquisition speed is desired. And, in some implementations, water-fat separation 504, water and fat mask generation steps 506 and 508, and water-fat removal steps 514 and 516 is not included. That is, water-fat removal is a specific problem for 1H-nuclear MRSI. However, water-fat removal may not be necessary when imaging other nuclei. However, in one example, water and nuisance lipid signals are removed using specialized methods as the first processing step for 1H magnetic resonance spectroscopic imaging experiments.
In one example, the disclosed method included dataset D1 containing 8×8 CSI encodings, each with 512 FID samples, D2 containing 64 echo shifts each having 128×128 spatial encodings and L=10. A CSI reconstruction included 12×12 spatial encodings while the EPSI reconstruction was obtained from 128×128 spatial encodings, both with 512 FID samples, resulting in improved resolution and SNR. The disclosed method achieves, in one example, both high spatial resolution (2.5 mm in plane) and improved SNR with comparable acquisition time. It is noted that the disclosed method includes a reconstruction having comparable resolution to a high-resolution CSI, having in one example a factor of 10 acceleration in data acquisition. In one observation, if a full (instead of sparsely sampled) (k,t)-space EPSI data set is collected for determination of the spatial coefficients, then the disclosed is equivalent to subspace denoising.
A technical contribution for the disclosed method and apparatus is that it provides for a computer-implemented apparatus and method for improving MR spectroscopic imaging.
In accordance with one embodiment, a device to acquire magnetic spectroscopic resonance imaging (MRSI) images includes a plurality of gradient coils positioned about a bore of a magnet, an RF transceiver controlled by a pulse module to transmit RF signals to an RF coil assembly and receive imaging data, and a computer. The computer can be programmed to acquire MRSI data with sparse sampling of (k,t)-space in variable density and SNR. The data can be rearranged to give one set with relatively high SNR and temporal resolution and another set with extended k-space coverage to provide the desired spatial resolution, wherein the second data set may have sparse sampling in both k-space and time to reduce data acquisition time depending on the practical experimental requirements, and reconstruct an spatiospectral function from the acquired MRSI data.
In accordance with another embodiment, a method of magnetic resonance spectroscopic imaging (MRSI) includes acquiring MRSI data covering a (k,t)-space sparsely with variable density and SNR, and generating a spatiospectral function using PS/low-rank models from the acquired data.
An implementation of system 100 in an example comprises a plurality of components such as one or more of electronic components, hardware components, and/or computer software components. A number of such components can be combined or divided in an implementation of the system 100. An exemplary component of an implementation of the system 100 employs and/or comprises a set and/or series of computer instructions written in or implemented with any of a number of programming languages, as will be appreciated by those skilled in the art. An implementation of system 100 in an example comprises any (e.g., horizontal, oblique, or vertical) orientation, with the description and figures herein illustrating an exemplary orientation of an implementation of the system 100, for explanatory purposes.
An implementation of system 100 in an example employs one or more computer readable signal bearing media. A computer-readable signal-bearing medium in an example stores software, firmware and/or assembly language for performing one or more portions of one or more implementations. An example of a computer-readable signal-bearing medium for an implementation of the system 100 comprises the recordable data storage medium of the image processor 122, and/or CPU 124. A computer-readable signal-bearing medium for an implementation of the system 100 in an example comprises one or more of a magnetic, electrical, optical, biological, and/or atomic data storage medium. For example, an implementation of the computer-readable signal-bearing medium comprises floppy disks, magnetic tapes, CD-ROMs, DVD-ROMs, hard disk drives, and/or electronic memory. In another example, an implementation of the computer-readable signal-bearing medium comprises a modulated carrier signal transmitted over a network comprising or coupled with an implementation of the system 100, for instance, one or more of a telephone network, a local area network (“LAN”), a wide area network (“WAN”), the Internet, and/or a wireless network.
When introducing elements of various embodiments of the present invention, the articles “a,” “an,” “the,” and “said” are intended to mean that there are one or more of the elements. The terms “comprising,” “including,” and “having” are intended to be inclusive and mean that there may be additional elements other than the listed elements. Furthermore, any numerical examples in the following discussion are intended to be non-limiting, and thus additional numerical values, ranges, and percentages are within the scope of the disclosed embodiments.
Dedicated hardware implementations including, but not limited to, application specific integrated circuits, programmable logic arrays and other hardware devices can likewise be constructed to implement the methods described herein. Application specific integrated circuits and programmable logic array can use downloadable instructions for executing state machines and/or circuit configurations to implement embodiments of the subject disclosure. Applications that may include the apparatus and systems of various embodiments broadly include a variety of electronic and computer systems. Some embodiments implement functions in two or more specific interconnected hardware modules or devices with related control and data signals communicated between and through the modules, or as portions of an application-specific integrated circuit. Thus, the example system is applicable to software, firmware, and hardware implementations.
In accordance with various embodiments of the subject disclosure, the operations or methods described herein are intended for operation as software programs or instructions running on or executed by a computer processor or other computing device, and which may include other forms of instructions manifested as a state machine implemented with logic components in an application specific integrated circuit or field programmable gate array. Furthermore, software implementations (e.g., software programs, instructions, etc.) including, but not limited to, distributed processing or component/object distributed processing, parallel processing, or virtual machine processing can also be constructed to implement the methods described herein. It is further noted that a computing device such as a processor, a controller, a state machine or other suitable device for executing instructions to perform operations or methods may perform such operations directly or indirectly by way of one or more intermediate devices directed by the computing device.
While the tangible computer-readable storage medium is an example embodiment as a single medium, the term “tangible computer-readable storage medium” should be taken to include a single medium or multiple media (e.g., a centralized or distributed database, and/or associated caches and servers) that store the one or more sets of instructions. The term “tangible computer-readable storage medium” shall also be taken to include any non-transitory medium that is capable of storing or encoding a set of instructions for execution by the machine and that cause the machine to perform any one or more of the methods of the subject disclosure. The term “non-transitory” as in a non-transitory computer-readable storage includes without limitation memories, drives, devices and anything tangible but not a signal per se.
The term “tangible computer-readable storage medium” shall accordingly be taken to include, but not be limited to: solid-state memories such as a memory card or other package that houses one or more read-only (non-volatile) memories, random access memories, or other re-writable (volatile) memories, a magneto-optical or optical medium such as a disk or tape, or other tangible media which can be used to store information. Accordingly, the disclosure is considered to include any one or more of a tangible computer-readable storage medium, as listed herein and including art-recognized equivalents and successor media, in which the software implementations herein are stored.
Although the present specification describes components and functions implemented in the embodiments with reference to particular standards and protocols, the disclosure is not limited to such standards and protocols. Each of the standards for Internet and other packet switched network transmission (e.g., TCP/IP, UDP/IP, HTML, HTTP, CoAP, MQTT) represent examples of the state of the art. Such standards are from time-to-time superseded by faster or more efficient equivalents having essentially the same functions. Wireless standards for device detection (e.g., RFID), short-range communications (e.g., Bluetooth®, WiFi, Zigbee®), and long-range communications (e.g., WiMAX, GSM, CDMA, LTE) can be used by a computer system to implement one or more functions described herein.
The illustrations of embodiments described herein are intended to provide a general understanding of the structure of various embodiments, and they are not intended to serve as a complete description of all the elements and features of apparatus and systems that might make use of the structures described herein. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The exemplary embodiments can include combinations of features and/or steps from multiple embodiments. Other embodiments may be utilized and derived therefrom, such that structural and logical substitutions and changes may be made without departing from the scope of this disclosure. Figures are also merely representational and may not be drawn to scale. Certain proportions thereof may be exaggerated, while others may be minimized. Accordingly, the specification and drawings are to be regarded in an illustrative rather than a restrictive sense.
Although specific embodiments have been illustrated and described herein, it should be appreciated that any arrangement which achieves the same or similar purpose may be substituted for the embodiments described or shown by the subject disclosure. The subject disclosure is intended to cover any and all adaptations or variations of various embodiments. Combinations of the above embodiments, and other embodiments not specifically described herein, can be used in the subject disclosure. For instance, one or more features from one or more embodiments can be combined with one or more features of one or more other embodiments. In one or more embodiments, features that are positively recited can also be negatively recited and excluded from the embodiment with or without replacement by another structural and/or functional feature. The steps or functions described with respect to the embodiments of the subject disclosure can be performed in any order. The steps or functions described with respect to the embodiments of the subject disclosure can be performed alone or in combination with other steps or functions of the subject disclosure, as well as from other embodiments or from other steps that have not been described in the subject disclosure. Further, more than or less than all of the features described with respect to an embodiment can also be utilized.
Less than all of the steps or functions described with respect to the exemplary processes or methods can also be performed in one or more of the exemplary embodiments. Further, the use of numerical terms to describe a device, component, step or function, such as first, second, third, and so forth, is not intended to describe an order or function unless expressly stated so. The use of the terms first, second, third and so forth, is generally to distinguish between devices, components, steps or functions unless expressly stated otherwise. Additionally, one or more devices or components described with respect to the exemplary embodiments can facilitate one or more functions, where the facilitating (e.g., facilitating access or facilitating establishing a connection) can include less than every step needed to perform the function or can include all of the steps needed to perform the function.
In one or more embodiments, a processor (which can include a controller or circuit) has been described that performs various functions. It should be understood that the processor can be multiple processors, which can include distributed processors or parallel processors in a single machine or multiple machines. The processor can be used in supporting a virtual processing environment. The virtual processing environment may support one or more virtual machines representing computers, servers, or other computing devices. In such virtual machines, components such as microprocessors and storage devices may be virtualized or logically represented. The processor can include a state machine, application specific integrated circuit, and/or programmable gate array including a Field PGA. In one or more embodiments, when a processor executes instructions to perform “operations”, this can include the processor performing the operations directly and/or facilitating, directing, or cooperating with another device or component to perform the operations.
The Abstract of the Disclosure is provided with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. In addition, in the foregoing Detailed Description, it can be seen that various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separately claimed subject matter.
While the disclosure has been described in detail in connection with a limited number of embodiments, it should be readily understood that the invention is not limited to such disclosed embodiments. Rather, that disclosed can be modified to incorporate any number of variations, alterations, substitutions or equivalent arrangements not heretofore described, but which are commensurate with the spirit and scope of the invention. Furthermore, while single energy and dual-energy techniques are discussed above, that disclosed encompasses approaches with more than two energies. Additionally, while various embodiments of the invention have been described, it is to be understood that disclosed aspects may include only some of the described embodiments. Accordingly, that disclosed is not to be seen as limited by the foregoing description, but is only limited by the scope of the appended claims.
One or more of the exemplary embodiments can include a device to acquire MRSI data. The device can include a plurality of gradient coils positioned about a bore of a magnet; an RF transceiver controlled by a pulse module to transmit RF signals to an RF coil assembly and receive imaging data; and a computer. The computer can be programmed to acquire MRSI data with sparse sampling of (k,t)-space in variable density and SNR. The data can be rearranged to give one set with relatively high SNR and temporal resolution and another set with extended k-space coverage (may be with sparse sampling in both k-space and time to reduce data acquisition time depending on the practical experimental requirements). The computer can reconstruct the desired spatiospectral function from the acquired MRSI data to provide biochemical information of the object being imaged.
In one embodiment, the computer can be further programmed to reconstruct the spatiospectral function by first determining a temporal basis and then determining spatial coefficients. Joint determination of the temporal (spectral) basis and the spatial coefficients is also possible using a variation of the disclosed method. If the temporal (spectral) basis is known (or predetermined in a different experiment, or based on prior knowledge), it is also possible to eliminate the step of determining the basis in the disclosed method. In one embodiment, the computer can reduce the effect of field inhomogeneity by mapping the field in homogeneity and then determine the temporal (or spectral) basis based on the mapped field inhomogeneity. In one embodiment, the computer can reconstruct the spatio spectral function by estimating a set of temporal basis functions from field inhomogeneity-corrected data and using the estimated temporal basis functions to reconstruct the spatiospectral function. In one embodiment, the computer can use a chemical shift imaging pulse sequence to generate the first data set, wherein the pulse sequence is one of a 2D or a 3D pulse sequence. In one embodiment, the computer can use an echo-planar spectroscopic imaging pulse sequence for rapid acquisition of spatiospectral encodings, where the pulse sequence is one of a 2D or a 3D pulse sequence. In one embodiment, the computer can generate fat-water images; separate water and fat; generate a water mask and a fat mask; and remove water and fat from the acquired MRSI data using the water mask and the fat mask. In one embodiment, the computer can reconstruct the spatiospectral function using a regularized least-squares optimization. In one embodiment, the computer can reconstruct the spatiospectral function by representing a desired function using a low-rank matrix multiplication model; estimating spatial coefficients; and forming the final spatiospectral function by combining the spatial coefficients and temporal bases.
One or more of the exemplary embodiments can be a method of MRSI. The method can include acquiring MRSI data with relatively high SNR and temporal resolution suitable for determination of temporal (or spectral basis) used in a subspace (or low-rank matrix/tensor) model; acquiring MRSI data with extended k-space coverage to give high spatial resolution.
In one embodiment, the method can include determining a temporal basis based on the MRSI data with relatively high SNR and temporal resolution; and determining spatial coefficients based on MRSI data with extended k-space coverage. In one embodiment, the method can include mapping the field inhomogeneity to reduce the effect of field inhomogeneity; and determining the temporal basis based on the mapped field inhomogeneity. In one embodiment, the method can include estimating a set of temporal basis functions from field inhomogeneity-corrected data; and using the estimated temporal basis functions to generate the desired spatiospectral function. In one embodiment, the method can include acquiring MRSI data with relatively high SNR and temporal resolution using a chemical shift imaging pulse sequence, where the pulse sequence is one of a 2D or a 3D CSI pulse sequence. In one embodiment, the method can include acquiring MRSI data with extended k-space coverage using an echo-planar spectroscopic imaging pulse sequence, where the pulse sequence is one of a 2D or a 3D pulse sequence. In one embodiment, the method can include generating fat-water echo images; separating water and fat; generating a water mask and a fat mask; and removing water and fat from the acquired MRSI data. In one embodiment, the method can include reconstructing the spatiospectral function using a regularized least-squares optimization. In one embodiment, the method can include reconstructing the spatiospectral function by representing a desired spatiotemporal function using a low-rank matrix (or tensor) multiplication model; estimating spatial coefficients; and forming the final spatiospectral function by combining the spatial coefficients and temporal bases.
In one embodiment, the regularization terms are either quadratic or nonquadratic. In one embodiment, the spatial support information is incorporated into the optimization problems. In one embodiment, the optimization problems can be solved by solving a system of linear equations if the regularization term is quadratic. In one embodiment, the computer is further programmed to generate: i) a set of anatomical images; and ii) a set of field inhomogeneity maps to be used in processing and reconstruction schemes. In one embodiment, the MRSI data are acquired at different echo times (in a time-interleaved data acquisition strategy). In one embodiment, the above described methods are used to acquire and process MRSI data from multiple coils. In one embodiment, acquisition and reconstruction is performed on systems with different magnetic field strengths. In one embodiment, spatiospectral functions are generated for chemical species from one or more of the following nuclei: 1H, 31P, 19F, 13C, 7Li, 15N, 17O, 19F, 23Na, and 129Xe.
One unique feature of the disclosed data acquisition strategy is variable sampling densities in (k,t)-space to enable high-speed MRSI data acquisition with sparse sampling in k-space and/or time; another unique feature of the disclosed sampling patterns is that it enables the use of subspace (or low-rank matrix or tensor) models to recover high-quality, high-resolution spatiospectral functions from sparsely sampled (k,t)-space MRSI data. The MRSI data collected using such sampling schemes can be rearranged in multiple ways. One way is to organize the acquired MRSI data into two data sets: one data set (called “training” data, denoted as D1) with relatively high SNR, and high temporal resolution (often covering central k-space region), and another data set (called “imaging” data, denoted as D2) with extended k-space coverage to achieve high spatial resolution. The “imaging” data can sample k-space and/or time sparsely with or without signal averaging to achieve high-speed data acquisition.
One or more of the exemplary embodiments can utilize one or more features of the imaging systems and/or imaging methods described in Appendix A of U.S. Provisional Application 62/102,315, filed Jan. 12, 2015, the disclosure of which is hereby incorporated by reference herein in its entirety. The one or more features of the imaging systems and/or imaging methods described in Appendix A can be utilized in place of and/or in addition to one or more features described herein with respect to the exemplary embodiments.
This application claims the benefit of U.S. Provisional Application Ser. No. 62/102,315, filed Jan. 12, 2015, which is incorporated by reference herein in its entirety including Appendices.
This invention was made with government support under contract numbers EB015904, EB001977, and EB013695 awarded by the National Institutes of Health. The government has certain rights in the invention.
| Number | Date | Country | |
|---|---|---|---|
| 62102315 | Jan 2015 | US |
