WO2012129433A2 - System and method of performing magnetic resonance spectroscopic imaging - Google Patents
System and method of performing magnetic resonance spectroscopic imaging Download PDFInfo
- Publication number
- WO2012129433A2 WO2012129433A2 PCT/US2012/030176 US2012030176W WO2012129433A2 WO 2012129433 A2 WO2012129433 A2 WO 2012129433A2 US 2012030176 W US2012030176 W US 2012030176W WO 2012129433 A2 WO2012129433 A2 WO 2012129433A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- magnetic resonance
- compartment
- interest
- resonance spectroscopy
- phase
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/485—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy based on chemical shift information [CSI] or spectroscopic imaging, e.g. to acquire the spatial distributions of metabolites
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
Definitions
- the field of the currently claimed embodiments of this invention relates to systems and methods of spatially localized magnetic resonance spectroscopy.
- MRS multi-voxel magnetic resonance spectroscopy
- CSI chemical shift imaging
- a method of performing spatially localized magnetic resonance spectroscopy includes receiving a magnetic resonance image of an object; identifying a plurality C of compartments that generate magnetic resonance spectroscopy signals in the object including at least one compartment of interest; segmenting in at least one spatial dimension the magnetic resonance image of the object into the C compartments; acquiring magnetic resonance spectroscopy signals from the compartments by applying a plurality of M' phase encodings applied in the at least one spatial dimension, wherein M' > C; calculating a spatially localized magnetic resonance chemical shift spectrum from the at least one compartment of interest; and rendering a spatially localized magnetic resonance spectrum that is substantially equal to a spatial average of magnetic resonance chemical shift spectra from the at least one compartment of interest.
- a magnetic resonance localized spectroscopy and imaging system includes a magnetic resonance imaging scanner and a data processing system configured to communicate with the magnetic resonance imaging scanner to receive magnetic resonance spectroscopy signals of an object.
- the data processing system is configured to receive a magnetic resonance image of the object; display the magnetic resonance image to permit identification of a plurality C of compartments that generate magnetic resonance spectroscopy signals in the object and that includes at least one compartment of interest; segment in at least one spatial dimension the magnetic resonance image of the object into the C compartments; receive magnetic resonance spectroscopy signals from the object corresponding to the magnetic resonance image by applying a plurality M' of phase encodings in at least one spatial dimension, wherein M' > C; calculate a spatially localized magnetic resonance chemical shift spectrum from the at least one compartment of interest; and provide a spatially localized magnetic resonance spectrum substantially equal to the spatial average of the magnetic resonance chemical shift spectra from the at least one compartment of interest.
- FIG. 1 is a schematic illustration of magnetic resonance localized spectroscopy and imaging system according to an embodiment of the current invention.
- FIG. 2 is a flow chart depicting implementation of the presently disclosed
- FIGS. 3A-3D show simulated 16-step phosphorus ( 31 P) one-dimensional (ID)
- the reconstructed SLAM chest (3B) and heart (3D) spectra are indistinguishable from the originals.
- FIGS. 4A-4E show cardiac model (4A) and Monte Carlo simulation of the effect of noise and 30% ( ⁇ 15%) inhomogeneity on the accuracy of SLAM signal reconstruction vs. CSI (4B-4E).
- the chest-to-heart signal ratio is held constant at 4 in (4B, 4C) depicted by the dark continuous curve in (4A).
- the ratio is 2.5 scaled by the experimental surface coil sensitivity profile depicted by the blue dashed curve in (4A).
- the largest errors in the heart correspond to configurations #1 (2cm chest, 2cm heart, no separation between chest and heart), #6 (2cm chest, 2cm heart, 1cm gap), #11 (3cm, 2cm, 0cm), and #16 (3cm, 2cm, 1cm).
- (5A) and (5C) are for a 4- voxel thick heart
- (5B) and (5D) are for a 3 -voxel-thick heart compartment, all with a 2- voxel thick chest compartment.
- FIGS. 6A-6D show the spatial response function for the heart compartment
- FIGS. 7A-7E show CSI and SLAM spectra reconstructed from the standard
- Philips Medical System's 31 P test phantom comprised of a H 3 P0 4 disk on the bottom (7A, 7C), and a H 3 P0 2 disk on top (7B, 7D), as shown in the image (7E).
- the CSI spectra (7 A, 7B) are the sum of the spectra from the voxels (red horizontal lines) containing the disks and were acquired with 16 phase-encoding gradients (-8... +7).
- the SLAM spectra (7 A, 7D) were acquired 4-times faster with just 4 phase-encodes (-1,-2,0,1).
- the SNR for the CSI spectra are 660 (7 A) and 638 (7B), compared to 528 (7C) and 482 (7D) for SLAM.
- the signal at ⁇ 0 ppm is a contaminant present only in the H 3 P0 2 disk.
- FIGS. 8A-8C show (8A) human leg 31 P spectrum acquired by SLAM (top) and CSI (lower) from the same 6-voxel volume in the same scan time (2.1 min).
- FIG. 8B 31 P spectra acquired from a normal human heart from the same 4-voxel volume, using ID CSI in 8.4 min, and SLAM spectra reconstructed with a subset of 4 central £-space phase- encodes and a 4-compartment model.
- FIG. 8C Spectra acquired with just two phase- encodes and a 2-compartment model (chest and heart). The effective SLAM acquisition times were l/4 th and l/8 th of CSI.
- FIGS. 9A-9D show fitting results reconstructed by SLAM from a subset of 4 of the 16 CSI phase encoding steps acquired from the 24 heart patients and control subjects, as compared to the CSI results.
- FIGS. 10A-10D show (10A) CSI, (10B) SLAM and (IOC) error-minimized fSLAM spectra, all normalized to constant noise on the same volunteer with the same total scan time and total voxel volume. Gradient encoding steps of -8 to +7 (integer) were used for standard CSI; integer steps -2, -1, 0, 1 repeated 4 times were used for SLAM; and fSLAM used non-integer steps -2.13, -0.73, +0.73, +2.13 repeated 4 times.
- FIGS. 11 is a table showing a comparison of some embodiments of the current invention with conventional approaches.
- FIGS. 12A and 12B show an example of ID cardiac 31 P MRS (same volume,
- FIGS. 13A and 13B show an example of 2D brain 1H MRS (same volume,
- FIGS. 14A and 14B show an example of a 3D phantom 31 P MRS (same volume, SLAM 100 times faster).
- SNR signal-to-noise ratio
- MRS magnetic resonance spectroscopy
- Matching voxels to anatomical compartments a priori yields better SNR than the spectra created by summing signals from constituent chemical-shift-imaging (CSI) voxels post-acquisition.
- CSI chemical-shift-imaging
- SLAM linear algebraic modeling
- SLAM CSI surface coil phosphorus MRS in phantoms, the human leg and the heart on a 3T clinical scanner. Its SNR performance, accuracy, sensitivity to registration errors and inhomogeneity, are evaluated. Compared to one-dimensional CSI, SLAM yielded quantitatively the same results 4-times faster in 24 cardiac patients and healthy subjects, and 45% higher cardiac SNR when applied pro-actively to 6 additional subjects.
- SLAM can be further extended according to an embodiment of the current invention with fractional phase- encoding gradients that optimize SNR and/or minimize both inter- and intra-compartmental contamination.
- fractional-SLAM fSLAM
- SLAM and fSLAM are simple to implement and reduce minimum scan times for CSI, which otherwise limits scan-time reductions achievable with higher SNR and field strengths.
- the broad concepts of the current invention are not limited to these particular embodiments and examples.
- Scan-time and signal-to-noise ratio are major problems for in vivo spatially localized magnetic resonance spectroscopy (MRS) of low-concentration metabolites.
- MRS magnetic resonance spectroscopy
- SNR is proportional to voxel size
- matching the voxel to the desired anatomical compartment a priori yields the best SNR for a fixed scan time [1].
- CSI chemical shift imaging
- the same principle applies in general wherever the CSI voxel size is smaller than the object of interest.
- the SNR gain factor for a fixed scan-time obtained by correctly encoding a compartment at the outset, as compared to adding signals from individual CSI voxels to form the equivalent-sized compartment post-acquisition, is: notwithstanding the effects of nonuniform sensitivity and concentration distributions, or differences in the integrated spatial response function (SRF).
- SRF integrated spatial response function
- SLIM[3], GSLIM[4] and SLOOP[5] could realize the g-fold SNR gain if the desired compartments were prescribed from scout MRI prior to acquisition, and if an appropriately SNR-optimized gradient set were then applied.
- the compartment's signal is modeled as the integral of phase-encoded signal contributions in each compartment, assumed homogeneous. The approach is prone to inter- and intra-compartmental errors when metabolite distributions are non-uniform between [6] and within each compartment, and as the number of phase-encoding gradient steps are reduced.
- GSLIM[4] and SLOOP[5] were introduced to minimize the inter-compartmental errors.
- GSLIM does this by applying non- Fourier, generalized series modeling to the SLIM result[4, 6].
- SLOOP minimizes the inter- compartmental error by optimizing the SRF for the desired compartment, ideally by specifically tailoring the phase-encoding gradient set for the acquisition[5].
- Several other proposed improvements add constraints to deal with inhomogeneity in the main (B 0 ) fleld[7- 9], registration errors[9], and multi-element receivers [10].
- SLAM linear algebraic modeling
- the SLAM pulse sequence differs in that the number of phase-encoding steps is essentially C, and they are always selected from the center of the integer-stepped k-space of CSI where SNR is highest. Other than determining the number, C, the need for image-guided gradient optimization, prescription and implementation at the scanner-side prior to acquisition, is thus avoided.
- SLAM g-fold SNR gains of 30-200% SNR are demonstrated in 3T phosphorus ( P) studies of the human leg and heart in vivo, compared to conventional [ 17- 22] one-dimensional (ID) CSI spectra from the same net volume and scan-time.
- fSLAM phase-encoding gradients
- the phase-encoding gradients are pro- actively optimized at the scanner-side to maximize SNR and/or minimize both the inter- compartmental leakage as well as the intra-compartmental errors produced by nonuniform signal distributions. Intra-compartmental errors have not been addressed in prior methods [3- 5].
- An example of fSLAM according to an embodiment of the current invention is demonstrated in pro-active human cardiac 31 P studies.
- FIG. 1 is a schematic illustration of a magnetic resonance spectroscopy and imaging (MRSI) system 100 according to an embodiment of the current invention.
- the MRSI system 100 includes a magnetic resonance scanner 101, a data storage unit 108, and a data processing unit 109.
- Magnetic resonance scanner 101 has a main magnet 105 providing a substantially uniform main magnetic field B 0 for a subject (or object) 102 under observation on scanner bed 103, a gradient system 106 providing a perturbation of the main magnetic field B 0 to encode spatial information of the constituent molecules of subject 102 under observation, and a radio-frequency (RF) coil system 107 to transmit electromagnetic waves and to receive magnetic resonance signals from subject 102.
- RF radio-frequency
- Data storage unit 108 may be, for example, a hard disk drive, a network area storage (NAS) device, a redundant array of independent disks (RAID), a flash drive, an optical disk, a magnetic tape, a magneto-optical disk, etc.
- NAS network area storage
- RAID redundant array of independent disks
- flash drive an optical disk
- magnetic tape a magnetic tape
- magneto-optical disk etc.
- the data storage unit 108 is not limited to these particular examples. It can include other existing or future developed data storage devices without departing from the scope of the current invention.
- the data processing system 109 is in communication with magnetic resonance scanner 101 to receive magnetic resonance signals for forming magnetic resonance images of subject 102.
- Data processing system 109 may be partially or totally incorporated within a structure housing magnetic resonance scanner 101.
- Data processing system 109 may be partially or totally incorporated in a workstation that is structurally separate from and in communication with magnetic resonance scanner 101.
- Data processing system 109 may be incorporated in a workstation that is structurally separate from and in communication with magnetic resonance scanner 101.
- An operator 113 may interact with the MRSI system 100 with input/output device 112.
- the data processing system 109 is configured to receive a magnetic resonance image of the object; display the magnetic resonance image to permit identification of a plurality C of compartments that generate magnetic resonance spectroscopy signals in the object and that includes at least one compartment of interest; segment in at least one spatial dimension the magnetic resonance image of the object into the C compartments; receive magnetic resonance spectroscopy signals from the object corresponding to the magnetic resonance image by applying a plurality M' of phase encodings in at least one spatial dimension, where M' > C; calculate a spatially localized magnetic resonance chemical shift spectrum from at least one compartment of interest; and provide a spatially localized magnetic resonance spectrum substantially equal to the spatial average of magnetic resonance chemical shift spectra from the at least one compartment of interest.
- the calculating of the spatially localized magnetic resonance chemical shift spectrum from the at least one compartment of interest can use a linear algebraic method.
- the magnetic resonance spectroscopy and imaging scanner can be further configured to permit identification and segmentation of the C compartments prior to receiving the magnetic resonance spectroscopy signals from the object.
- the magnetic resonance spectroscopy and imaging scanner can be further configured to optimize the M' phase encodings for at least one of the signal-to-noise ratio or the spatial selection in the at least one compartment of interest.
- the plurality M' of phase encodings can be provided from a central portion of a k-space corresponding to the at least one spatial dimension of the magnetic resonance image of the object.
- the at least one of the plurality M' of phase encodings can be an integer multiple of a smallest non-zero phase encoding.
- the at least one of the plurality M' of phase encodings can be a non-integer multiple of a smallest non-zero phase encoding.
- the plurality M' of phase encodings can be determined from the magnetic resonance image after segmentation, the data processing system can be further configured to perform at least one of (1) optimization of the signal-to-noise ratio in the at least one compartment of interest which includes the selection of phase encodings from a central portion of the k-space of the magnetic resonance image of the object, or (2) optimization of the spatial selection of the at least one compartment of interest by minimization of at least one of the magnetic resonance spectroscopy signal arising from outside of the compartment of interest, or of erroneous signals arising from non-uniform magnetic resonance spectroscopy signal distributions arising within the compartment of interest.
- the plurality M' of phase encodings can be provided by including a metric that optimizes both the signal-to-noise ratio and the spatial selection in the at least one compartment of interest.
- the at least one spatial dimension is one of two spatial dimensions or three spatial dimensions
- the plurality M' of phase encodings includes two or three subsets of phase encodings that are applied in the two or the three spatial dimensions, respectively.
- the number of phase encodings provided in each of the subset of phase encodings is greater than or equal to the number of compartments generating magnetic resonance spectroscopy signals that are segmented in the corresponding spatial dimensions of the object.
- the at least one compartment of interest can be a plurality of compartments of interest.
- Each row of the known signal matrix, M*N , on the left side of the equation is an N-point array, where N is the number of time-domain data points.
- the first matrix on the right side is the phase-encoding FT operator (PE), and each term of the unknown spectral matrix, p, is also an N-point array.
- PE phase-encoding FT operator
- the goal of the CSI experiment is to reconstruct the M unknown spectra in matrix p of Eq. (3), from the M known signals (S) acquired with M different phase-encodes.
- S M known signals
- p has just C ⁇ M MRS compartments of interest, as well as the spatial position of each compartment. Theoretically, only C measurements with C phase-encoding steps are needed to unambiguously solve p and reconstruct the C spectra.
- Eq. (5) is now over-determined and can be solved with three phase-encoding rows.
- the minimum number of required phase-encoding steps is reduced from 4 to 3.
- ⁇ C N [ s a submatrix of MxM X ⁇ ⁇ ⁇ retaining the C non-eliminated rows; is a submatrix of M xM MxM that retains the C columns corresponding to the C non-
- M xN is a submatrix of MxN signals acquired from the sample using a subset of M' «M phase-encoding steps.
- Solution of Eq. (8) results in a set of spectra, each of which closely approximates the average spectrum of each ID CSI compartment.
- Steps 1-5 are performed with Steps 1-5 as follows:
- the choice of the M' phase-encoding steps need not be limited to the original basis set of M CSI steps corresponding to integer k's in Eq. (3).
- the M' phase-encoding gradients can be chosen to optimize desired properties of the reconstruction. For example, as we now show, the gradients can be optimized to maximize the SNR, and/or minimize the inter-compartmental signal contamination, and/or mimmize the intra-compartmental error due to nonuniform signal sources. This effectively involves allowing for fractional k's in the CSI Eq. (3), with all other experimental parameters left unchanged. Unlike SLAM, this fractional SLAM method, denoted fSLAM, does require scanner-side gradient optimization and prescription. fSLAM with maximum SNR
- Eq. (9) is modified to include noise terms the time-domain signal:
- ⁇ CxN + ⁇ CxN PE CxM' X (S M , xN + f M' xN (10) where CxN is the noise in the reconstructed spectra.
- each row in ⁇ MxN is an average spectrum of its compartment and each row in MxN is the deviation of the true spectrum from its compartmental average.
- each row in MxN is the deviation of the true spectrum from its compartmental average.
- PE CxM x S M' *N ⁇ Cx g N +PE ⁇ xM . x PE M . xM x K" (15) PE + x PE . x jO inhom )
- PE' PE + x PE
- This coefficient set has a smaller sum-of-the-squares and is not affected by differences in the mean coefficient of each compartment.
- ⁇ c (16) [0056]
- vt ⁇ . can reflect, for example, intrinsic differences in metabolite concentrations between compartments.
- the fSLAM experiment is performed using the same Steps 1-5 as the SLAM protocol ( Figure 2) except that the phase-encoding gradients in Step 2 are obtained by minimizing either the SNR cost-function in Eq. (12) or the error cost-function in Eq. (18).
- the different optimizations will result in different sets of phase- encoding gradients. If a gradient set optimized for both SNR and minimum error is being sought, minimization of the sum of the cost-functions in Eqs. (12) and (18) cannot be used because their scales differ. Instead, minimization of a weighted sum of the ratio of cost functions for fSLAM to those for SLAM can suffice. The choice of the weighting will depend on the application and error tolerance.
- the phase-encoding gradients in Step 2 are typically fractional.
- the first integral in Eq. (21) is the chest to heart leakage, ⁇ . f (x) ⁇ Af (x)
- a random inhomogeneity of ⁇ 15% (30% total) in the resultant signal was then simulated for both scenarios.
- the mean signal was determined for each compartment by adding signals from the corresponding voxels of the full CSI set to serve as a reference. Then, white noise was added such that the SNR in the heart compartment was 20.
- the mean error and the SD of the error were determined after 1000 Monte Carlo simulation runs.
- SLAM with SLIM[3], with respect to registration errors A ID cardiac 31 P model with chest from -60mm to -30mm, heart from -30mm to 10mm, and a chest-to-heart signal ratio of 4 was assumed as in scenario- 1, above ( Figure 4A).
- a random segmentation error between -2 mm and +2 mm was introduced at the edges of either compartment: (i) with the chest and heart stationary (no partial volume error); and (ii) with the chest and heart also moved by ⁇ 2mm (partial volume error). The chest was constrained never to overlap the heart.
- Both SLAM and SLIM were simulated with four CSI phase-encodes from central k-space.
- SLIM reconstruction was performed as prescribed[3], by integrating the phase-encoding coefficients over the 3 -compartment model of heart, chest and Other' and generating a 4x3 'G'-matrix[3].
- the mean ( ⁇ SD) % error between the reconstructed signal and the true or the CSI result was calculated for 1000 runs.
- SRF h was calculated from Eq. (19) for 4-step SLAM, 4-step fSLAM, 16-step
- MRI/MRS system on phantoms, the human leg, and the human heart.
- the phantom studies were done with a 14-cm diameter single loop transmit/receive coil, and the human studies used a 17-cm/l 1-cm diameter dual loop transmit and a 8-cm diameter single loop receive P coil set described previously [22]. All human studies were approved by the Johns Hopkins Medicine Institutional Review Boards and all participants provided informed consent.
- the individual CSI spectra from all of the volume elements constituting each compartment were co-added post-acquisition for all comparisons of spectra from the equivalent volumes reconstructed using SLAM and fSLAM.
- Phantom studies were performed on two standard Philips 31 P test disks 15-cm in diameter and 2.5-cm thick. One contained 300 mM H 3 P0 2 , the other had 300 mM H 3 P0 .
- FSC frequency-sweep-cycled
- AHP adiabatic half passage
- Spectra were fit by the circle-fit method[25] to provide a quantitative comparison of PCr and ⁇ - ⁇ peak areas measured by SLAM with those from conventional CSI (the localization and spectral analysis method are independent). Spectra were exponential-filtered (15-Hz line-broadening) and zero-filled 4 times to 2048 points.
- SLAM used the same 4 middle k-space phase-encoding steps for each exam, repeated four times for the same total scan-time as CSI.
- fSLAM phase-encoding employed 4, typically-fractional gradient steps, specifically optimized for minimum compartmental leakage in the heart compartment for each volunteer, after manual segmentation of the scout MRI using the scanner's cursor function. As in the simulations, optimization was performed using Matlab on a lap-top computer at the scanner-side, with weighting factors set to unity. The four gradient values were manually entered as experimental parameters in the fSLAM pulse sequence on the scanner. The four steps were repeated four times for the same total scan- time as the CSI.
- Figures 3A-3D show that SLAM spectra of the chest and heart, reconstructed using only the three middle (k-space) phase-encoding steps of the original 16, are indistinguishable from the original simulated spectra in the absence of inhomogeneity or noise.
- the effect of adding noise and inhomogeneity on SLAM spectra reconstructed for a range of different chest-muscle and heart compartment distributions, is illustrated by the Monte Carlo simulations for both models of concentration and sensitivity variations in Figures 4A-4E. These show that the accuracy of the reconstruction, as indexed by the mean of the error ⁇ 10% for all chest/heart anatomical combinations. As might be expected, the higher the concentration or larger the compartment size, the smaller the error SD. For the heart, the simulations predict highest errors when the effective extent of the cardiac compartment is smallest.
- Table 1 Monte Carlo analysis of the effect of ⁇ 2mm misregistration on accuracy of cardiac PCr measurements for a 30mm chest/40mm heart model.
- compartment error factor, for SLAM and fSLAM, as compared with 16- voxel ID CSI of the heart, are shown in Figure 5 with 3- and 4- voxel thick cardiac compartments.
- the maximum SNR results from choosing the phase-encoding steps closest to central k-space. Because SLAM is confined to the CSFs integer phase-encoding set, its SNR advantage fades as more-and-more of the high k-space phase-encodes are used. Optimum SNR performance for SLAM occurs when the number of phase-encodes approximates the number of compartments, wherein its performance approximates that of fSLAM.
- the best strategy is to choose the M' ⁇ C non-equal CSI phase-encoding steps at or closest to the center of k-space, and repeat or average the acquisitions up to the allotted scan time, rather than add any higher k-space phase-encodes.
- fSLAM always achieves 1.5- 1.8 times the SNR of standard CSI independent of the number of phase-encoding steps that are allowed. This reflects the fact that fSLAM is free to choose an array of fractional phase- encodes that all fall close to central k-space.
- the additiona -encodes offer the added benefit of reduced signal bleed ( Figures 5C, 5D).
- the errors, for fSLAM decay faster than SLAM as phase-encodes are added, indicating better error suppression, with larger compartments generating less error than smaller ones.
- Figures 6A-6D plots SRF h for 16- and 4-step CSI, 4-step SLAM and 4-step fSLAM. It is important to recognize that the signal derives from the integral of the curve over each compartment, resulting in cancellation of signal outside the heart.
- the cancellation is essentially perfect in the case of SLAM and fSLAM but not CSI (Table 2, first row).
- the signal in the chest compartment varies by up to 30% peak-to-peak, the upper bound for contamination of the heart compartment rises to 12-14% for SLAM and fSLAM. This compares to 9% for 16-step CSI, while the 4-step CSI is basically unusable (Table 2, second row).
- H 3 P0 4 has a single 31 P peak at 2.9 ppm, while the H 3 P0 2 resonance is a triplet centered about 13.5 ppm (coupling constant, 545 Hz), due to heteronuclear coupling with hydrogen.
- the SLAM spectra from the two disks are very similar to the summed CSI spectrum from the same compartment volumes, with negligible leakage consistent with the simulations (Figure 4).
- the SLAM spectrum has 2.1 times better SNR than CSI, and shows negligible signal contamination or bleed from the H 3 PO4 phantom positioned above the leg.
- Figure 8B shows 3 P heart spectra from a 16-step ID CSI (averaging 4 voxels), 4-step SLAM and 2-step SLAM from the same volume.
- the baseline roll is due to the acquisition delay for the phase- encoding gradient.
- CSI spectra from 8 healthy subjects and 16 patients are presented in Figure 9 .
- the ratio of PCr signal in chest to that in heart compartments was at or below ⁇ 5.
- the PCr and ⁇ - ⁇ peak areas from the SLAM reconstruction agree with those from CSI reconstruction ( Figures 9A, 9B).
- the myocardial PCr/ ATP ratio for the pooled patients and healthy subjects was the same (1.94 ⁇ 0.60 in CSI vs. 1.90 ⁇ 0.67 in SLAM), consistent with negligible contamination from chest muscle with its much higher PCr/ ATP ratio of ⁇ 4[19].
- the data also show that the total of the PCr in the chest plus that in the heart compartments measured by SLAM, is equal to the total measured by CSI ( Figure 9C).
- the fraction of cardiac PCr to the total PCr measured by SLAM is also equal to that measured by CSI ( Figure 9D). This means that the contamination of heart spectra from chest muscle in SLAM is not distinguishable from that in CSI. This result is consistent with Table 2. Importantly, all these SLAM results correspond to acquisitions effectively taking 1/4 th of the scan-time of CSI.
- Figures 10A-10D compare CSI, SLAM and fSLAM 31 P cardiac spectra proactively acquired in the same total scan-time from the same volume size in the same healthy volunteer.
- the time taken to implement fSLAM scanner-side was l-2min to manually segment the scout MRI, plus several seconds to optimize the gradient set on the lap-top computer.
- the SLAM and fSLAM spectra both have higher SNR than CSI from the same volume, while a possible bleed signal from the coil marker in the SLAM spectrum is absent in the fSLAM spectrum.
- Table 3 lists the SNR of human cardiac PCr in same-sized voxels for CSI, SLAM and fSLAM in 6 volunteers acquired in the same scan time.
- the mean SNR improvement for SLAM vs CSI for the six studies is 1.42 ⁇ 0.23.
- the mean SNR improvement for fSLAM vs CSI for the six studies is 1.34 ⁇ 0.19.
- this SNR gain would be consistent with a cardiac compartment equivalent to two of the 1-cm CSI voxels even though the reconstruction assumed a 4-voxel cardiac compartment. This likely reflects the combined effect of the decline in surface coil sensitivity with depth, and the 1-2 cm thickness of the anterior myocardial wall.
- the PCr/ ATP ratio for SLAM and fSLAM was not significantly different from that measured in either the first or the repeated last CSI scans, and the absolute metabolite signal levels do not change (Figure 10D).
- the highest SNR efficiency of CSI is only realized when the spatial resolution imposed at the time of acquisition, matches the desired compartment size[l].
- CSI's spatial resolution is usually set not by the size of the desired compartment, but by the geometry of the tissue that it must be distinguished from (eg, the chest in heart or liver studies, the scalp in brain studies).
- Alternative approaches that localize spectra to pre-selected compartments based on anatomical MRI information, are not new.
- the SLIM, GSLIM and SLOOP methods were originally proposed some 20 years ago[3-5], but see little use today compared to CSI or even PRESS, STEAM or ISIS.
- SLAM differs from SLIM, GSLIM and SLOOP in both the pulse sequence that is applied, and in MRS reconstruction.
- the SLAM pulse sequence is based on a CSI sequence from which essentially all of the high-order gradient phase-encoding steps are eliminated except for the C phase-encoding steps closest to central k-space.
- SLAM SLAM aims to generate spectra that are at best equal to the compartmental average CSI spectra
- SLIM, GSLIM and SLOOP use MRI-based constrained reconstruction or SRF optimization to obtain optimally-localized compartment spectra. Because of the relatively coarse resolution of CSI, this renders SLAM relatively insensitive to registration errors in segmenting the compartments-compared to SLIM for example (Table 1), where problems were noted previously[l 1, 29].
- SLAM can be performed retroactively on raw data sets that are accompanied by a scout MRI, simply by applying the algorithm to a subset of frames in each CSI data set. The result can be compared with the summed CSI from the same-sized compartments analogous to Figure 9.
- fSLAM extends SLAM by removing the limitation that the phase-encodes be selected from the set of integer-stepped CSI gradients. Instead, they are adjusted to minimize leakage or errors due to inhomogeneity and/or maximize SNR. We observed that maximizing SNR alone can produce unacceptable error if the clustering of phase-encodes at the center of k-space is unchecked (Figure 5).
- the SRF is not global but is specific to the cardiac model. Inter- and intra-compartmental leakage occurs only when the integral over the entire compartment is non-zero or in the presence of significant heterogeneity. Compartmental segmentation in SLAM ensures that the integral of the SRF vanishes over other compartments, while fSLAM minimizes the effect of heterogeneity within the compartment of interest as well. Ultimately however, the spatial responses for SLAM and fSLAM and their compartmental contamination are fully characterized by determining the accuracy of the solutions and leakage errors, for which CSI is used as the standard in the current work ( Figures 5, 7-10, Tables 1, 2).
- the SLAM and fSLAM methods yield spectra comparable to the average of same-sized CSI compartments but with large scan-time reductions, SNR gains, and manageable, if not insignificant, bleed artifacts.
- the SNR gains predicted by Eq. (1) will be moderated in practice by the depth-dependence of the surface coil sensitivity, as well as the actual metabolite distribution (in our case, the myocardial wall thickness).
- SLAM and fSLAM Independent of the SNR gain, SLAM and fSLAM reduce the minimum scan-time required for localization from M acquisitions in CSI, to C or M' «M.
- SLAM is implemented with the steps: (i) Acquire MRI; (ii) Segment MRI into C compartments and overlay on the CSI grid; (iii) Apply M' central k-space phase encodes; and (iv) Reconstruct the spectra using SLAM.
- 2D- and 3D SLAM experiments were done in a 3T Philips MRI system on the human brain and a phosphate phantom ( 31 P ), respectively.
- the compartments were: scalp, brain, lateral ventricle and background (2D); and H 3 P0 4 , H 3 P0 2 disk phantoms plus background (3D).
- An additional 31 P heart study with chest, heart and background (ID) compartments is also shown.
- SLAM spectra were reconstructed with central 4 (ID), 7x7 (2D) and 2x4x2 (3D) phase-encodes, and compared with compartmental average CSI spectra obtained from the whole datasets with 16 (ID), 32x25 (2D), and 10x20x8 (3D) phase encodes.
- ( ⁇ C) must be chosen for each of the two or three spatial dimensions.
- subsets of M' x , M' y , and M' z gradients must be chosen to apply phase-encoding in the X-, y-, and z- Cartesian directions respectively, as provided by the MRI/MRS scanner's spatially encoding gradient system.
- M' x + M' y , + M' Z M'.
- M' x , M' y precede or M' z is omitted depending on which dimensions are being encoded.
- the number of phase encodes used in any one direction, M' x , M'ylit or M' z , should at least not be less than the number of signal-generating compartments C x , C y , and C Z; that can be segmented in that dimension, ie, M' x >C X , M' y necessarily >Cy ; and M' z >C Z , etc.
- the particular strategies for choosing the individual phase-encoding gradients in the subsets for each dimension, M' x , M' y , and M' 2 are as described above applied to each of the dimensions.
- the gradients are chosen from central k-space (for each dimension, k x , k y , and k z ), and can be integer multiples for SLAM, or fractional multiples for fSLAM, while not repeating the zero phase-encode (other than for the purpose of signal averaging).
- gradient optimization for fSLAM can be treated as a separate application of the optimization algorithm for each M' x , M' y , and M' z in each dimension.
- Figures 12A and 12B show 3I P spectra for the same-sized chest (12A) and heart (12B) compartments reconstructed from ID CSI and SLAM.
- Figures 13A and 13B show 2D SLAM and CSI spectra from the same brain (13A) and lateral ventricle (13B) volume.
- Figures 14A and 14B show H 3 Po 2 (14A) and H 3 Po 4 (14B) phantom spectra reconstructed from 3D CSI and SLAM.
- the speedup for the ID, 2D and 3D SLAM compared with CSI are 4-, 16- and 100-fold, respectively: the SNR cost is 14%, ⁇ 30% and 50%.
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Health & Medical Sciences (AREA)
- Optics & Photonics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014501253A JP6070963B2 (en) | 2011-03-22 | 2012-03-22 | System and method for performing magnetic resonance spectroscopy imaging |
US14/006,069 US10209330B2 (en) | 2011-03-22 | 2012-03-22 | System and method of performing magnetic resonance spectroscopic imaging |
EP12760705.9A EP2688471B1 (en) | 2011-03-22 | 2012-03-22 | System and method of performing magnetic resonance spectroscopic imaging |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161466188P | 2011-03-22 | 2011-03-22 | |
US61/466,188 | 2011-03-22 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2012129433A2 true WO2012129433A2 (en) | 2012-09-27 |
WO2012129433A3 WO2012129433A3 (en) | 2012-12-27 |
Family
ID=46880047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2012/030176 WO2012129433A2 (en) | 2011-03-22 | 2012-03-22 | System and method of performing magnetic resonance spectroscopic imaging |
Country Status (4)
Country | Link |
---|---|
US (1) | US10209330B2 (en) |
EP (1) | EP2688471B1 (en) |
JP (1) | JP6070963B2 (en) |
WO (1) | WO2012129433A2 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9791489B2 (en) | 2011-03-22 | 2017-10-17 | The Johns Hopkins University | High dynamic range RF power monitor |
WO2012129433A2 (en) | 2011-03-22 | 2012-09-27 | The Johns Hopkins University | System and method of performing magnetic resonance spectroscopic imaging |
WO2015126499A2 (en) * | 2013-12-02 | 2015-08-27 | Andrew Irish | Systems and methods for gnss snr probabilistic localization and 3-d mapping |
EP2983193B1 (en) * | 2014-08-05 | 2021-10-20 | Aselta Nanographics | Method for determining the parameters of an ic manufacturing process model |
US10120049B2 (en) * | 2015-05-18 | 2018-11-06 | The Johns Hopkins University | System and method of obtaining spatially-encoded NMR parameters from arbitrarily-shaped compartments and linear algebraic modeling |
SG10201912062UA (en) | 2015-07-17 | 2020-02-27 | Univ California | System and method for localization and tracking using gnss location estimates, satellite snr data and 3d maps |
US10386493B2 (en) | 2015-10-01 | 2019-08-20 | The Regents Of The University Of California | System and method for localization and tracking |
DE102018208118A1 (en) * | 2018-05-23 | 2019-11-28 | Robert Bosch Gmbh | Method and apparatus for authenticating a message transmitted over a bus |
Family Cites Families (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3620901A1 (en) | 1986-06-21 | 1988-01-14 | Krupp Gmbh | CUTTING TOOL |
US5214372A (en) | 1989-05-11 | 1993-05-25 | Nokia Mobile Phones Ltd. | Linearizing circuit for dectection of a RF-power sensor |
JP3104985B2 (en) * | 1990-10-23 | 2000-10-30 | 株式会社東芝 | Magnetic resonance diagnostic equipment |
US5128629A (en) | 1991-04-22 | 1992-07-07 | Hughes Aircraft Company | Method for controlling the output power of digital cellular telephones |
JP3538249B2 (en) * | 1995-01-19 | 2004-06-14 | ジーイー横河メディカルシステム株式会社 | MRIS device |
US5982165A (en) | 1996-11-29 | 1999-11-09 | Wavetek Corporation | RF circuit for use in a combined leakage detector and signal level monitor |
US6310646B1 (en) | 1996-11-29 | 2001-10-30 | Wavetek Corporation | Method and apparatus for measuring a radio frequency signal having program information and control information |
US5903149A (en) | 1997-04-11 | 1999-05-11 | The Trustees Of The University Of Pennsylvania | Three-dimensional localized proton NMR spectroscopy using a hybrid of one-dimensional hadamard with two-dimensional chemical shift imaging |
US6255635B1 (en) | 1998-07-10 | 2001-07-03 | Ameritherm, Inc. | System and method for providing RF power to a load |
JP4416221B2 (en) * | 1999-09-28 | 2010-02-17 | 株式会社日立メディコ | Magnetic resonance imaging system |
US6426623B1 (en) | 2000-09-27 | 2002-07-30 | Mayo Foundation For Medical Education And Research | MRI RF power monitor |
US6853227B2 (en) | 2001-04-17 | 2005-02-08 | K-Tek Corporation | Controller for generating a periodic signal with an adjustable duty cycle |
GB2385674B (en) | 2002-02-20 | 2005-04-06 | Racal Instruments Ltd | Radio frequency power generation and power measurement |
ATE470998T1 (en) | 2002-03-07 | 2010-06-15 | Nokia Corp | POWER CONTROL DEVICE AND METHOD FOR CALIBRATING THE POWER OF A TRANSMITTER OR RECEIVER IN A MOBILE COMMUNICATIONS NETWORK |
US6803818B2 (en) | 2002-11-26 | 2004-10-12 | Agere Systems Inc. | Method and apparatus for improved output power level control in an amplifier circuit |
US7460837B2 (en) | 2004-03-25 | 2008-12-02 | Cisco Technology, Inc. | User interface and time-shifted presentation of data in a system that monitors activity in a shared radio frequency band |
US20070242715A1 (en) | 2006-04-18 | 2007-10-18 | Johan Gustavsson | Mode and polarization control in vcsels using sub-wavelength structure |
US20070249928A1 (en) * | 2006-04-19 | 2007-10-25 | General Electric Company | Method and system for precise repositioning of regions of interest in longitudinal magnetic resonance imaging and spectroscopy exams |
US7683614B2 (en) * | 2006-04-27 | 2010-03-23 | Stefan Posse | Magnetic resonance spectroscopy with sparse spectral sampling and interleaved dynamic shimming |
WO2007131187A2 (en) | 2006-05-05 | 2007-11-15 | Viasat, Inc. | Rf power sensor with chopping amplifier |
US7899432B2 (en) | 2006-06-19 | 2011-03-01 | California Institute Of Technology | Submillimeter wave heterodyne receiver |
US7902991B2 (en) | 2006-09-21 | 2011-03-08 | Applied Materials, Inc. | Frequency monitoring to detect plasma process abnormality |
US20080177163A1 (en) | 2007-01-19 | 2008-07-24 | O2 Medtech, Inc. | Volumetric image formation from optical scans of biological tissue with multiple applications including deep brain oxygenation level monitoring |
KR101727290B1 (en) | 2007-07-27 | 2017-04-27 | 페어차일드코리아반도체 주식회사 | Converter and the driving method thereof |
WO2009032336A1 (en) | 2007-09-07 | 2009-03-12 | The Johns Hopkins University | Sar dosimeter for rf power deposition in mri and methods and systems related thereto |
WO2009070695A1 (en) * | 2007-11-26 | 2009-06-04 | Mayo Foundation For Medical Education And Research | System and method for phase offset and time delay correction in magnetic resonance spectroscopy data |
US7822565B2 (en) | 2007-12-31 | 2010-10-26 | Advanced Energy Industries, Inc. | System, method, and apparatus for monitoring characteristics of RF power |
EP2338563B1 (en) | 2009-12-22 | 2013-11-20 | BIOTRONIK CRM Patent AG | MRI signaling |
WO2011153521A2 (en) | 2010-06-04 | 2011-12-08 | Mclean Hospital Corporation | Multi-modal imaging of blood flow |
US9791489B2 (en) | 2011-03-22 | 2017-10-17 | The Johns Hopkins University | High dynamic range RF power monitor |
WO2012129433A2 (en) | 2011-03-22 | 2012-09-27 | The Johns Hopkins University | System and method of performing magnetic resonance spectroscopic imaging |
JP6076677B2 (en) | 2011-11-25 | 2017-02-08 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging system |
-
2012
- 2012-03-22 WO PCT/US2012/030176 patent/WO2012129433A2/en active Application Filing
- 2012-03-22 EP EP12760705.9A patent/EP2688471B1/en active Active
- 2012-03-22 JP JP2014501253A patent/JP6070963B2/en active Active
- 2012-03-22 US US14/006,069 patent/US10209330B2/en active Active
Non-Patent Citations (37)
Title |
---|
A. BASHIRD.A. YABLONSKIY: "Natural linewidth chemical shift imaging (NL-CSI", MAGNETIC RESONANCE IN MEDICINE, vol. 56, 2006, pages 7 - 18 |
A.-M. EL-SHARKAWYM. SCHARR. OUWERKERKR.G. WEISSP.A. BOTTOMLEY: "Quantitative cardiac 31P spectroscopy at 3 Tesla using adiabatic pulses", MAGNETIC RESONANCE IN MEDICINE, vol. 61, 2009, pages 785 - 795 |
C.S. SMITHP.A. BOTTOMLEYS.P. SCHULMANG. GERSTENBLITHR.G. WEISS: "Altered Creatine Kinase Adenosine Triphosphate Kinetics in Failing Hypertrophied Human Myocardium", CIRCULATION, vol. 114, 2006, pages 1151 - 1158 |
DONG Z ET AL., MRM, vol. 55, 2006, pages 1447 - 1453 |
H.R. BROOKERT.H. MARECIJ. MAO: "Selective Fourier transform localization", MAGNETIC RESONANCE IN MEDICINE, vol. 5, 1987, pages 417 - 433, XP002633857 |
HU X ET AL., MRM, vol. 8, 1988, pages 314 - 322 |
I. KHALIDOVD. VAN DE VILLEM. JACOBF. LAZEYRASM. UNSER: "BSLIM: Spectral Localization by Imaging With Explicit B0 Field Inhomogeneity Compensation", MEDICAL IMAGING, IEEE TRANSACTIONS, vol. 26, 2007, pages 990 - 1000, XP011186646, DOI: 10.1109/TMI.2007.897385 |
J. FRAHMH. BRUHNM.L. GYNGELLK.D. MERBOLDTW. HANICKER. SAUTER: "Localized high-resolution proton NMR spectroscopy using stimulated echoes: Initial applications to human brain in vivo", MAGNETIC RESONANCE IN MEDICINE, vol. 9, 1989, pages 79 - 93, XP000004945 |
J.A. KMIECIKC.D. GREGORYZ.P. LIANGP.C. LAUTERBURM.J. DAWSON: "Lactate quantitation in a gerbil brain stroke model by GSLIM of multiple-quantum-filtered signals", J MAGN RESON IMAGING, vol. 9, 1999, pages 539 - 543 |
L. ANS. WARACHJ. SHEN: "Spectral localization by imaging using multielement receiver coils", MAGNETIC RESONANCE IN MEDICINE, 2011 |
L.P. PANYCHL. ZHAOR.V. MULKERN: "PSF-choice: A novel MRI method for shaping point-spread functions in phase-encoding dimensions", MAGNETIC RESONANCE IN MEDICINE, vol. 54, 2005, pages 159 - 168 |
LIANG ZP ET AL., IEEE TMI, vol. 10, 1991, pages 132 - 137 |
LOFFLER R ET AL., JMR, vol. 134, 1998, pages 287 - 299 |
M. BEERT. SEYFARTHJ. SANDSTEDEW. LANDSCHUTZC. LIPKEH. KOSTLERM. VON KIENLINK. HARRED. HAHNS. NEUBAUER: "Absolute concentrations of high-energy phosphate metabolites in normal, hypertrophied, and failing human myocardium measured noninvasively with 3 IP-SLOOP magnetic resonance spectroscopy", J AM COLL CARDIOL, vol. 40, 2002, pages 1267 - 1274 |
M. JACOBX. ZHUA. EBELN. SCHUFFZ.P. LIANG: "Improved model-based magnetic resonance spectroscopic imaging", MEDICAL IMAGING, IEEE TRANSACTIONS, vol. 26, 2007, pages 1305 - 1318, XP011193207, DOI: 10.1109/TMI.2007.898583 |
M. MEININGERW. LANDSCHUTZM. BEERT. SEYFARTHM. HORNT. PABSTA. HAASED. HAHNS. NEUBAUERM. VON KIENLIN: "Concentrations of human cardiac phosphorus metabolites determined by SLOOP 31P NMR spectroscopy", MAGNETIC RESONANCE IN MEDICINE, vol. 41, 1999, pages 657 - 663 |
M. VON KIENLINM. BEERA. GREISERD. HAHNK. HARREH. KOSTLERW. LANDSCHUTZT. PABSTJ. SANDSTEDES. NEUBAUER: "Advances in human cardiac 31P-MR spectroscopy: SLOOP and clinical applications", J MAGN RESON IMAGING, vol. 13, 2001, pages 521 - 527 |
M. VON KIENLINR. MEJIA: "Spectral localization with optimal pointspread function", JOURNAL OF MAGNETIC RESONANCE, vol. 94, 1991, pages 268 - 287, XP023957643, DOI: 10.1016/0022-2364(91)90106-4 |
M.A. CONWAYP.A. BOTTOMLEYR. OUWERKERKG.K. RADDAB. RAJAGOPALAN: "Mitral regurgitation: Impaired systolic function, eccentric hypertrophy, and increased severity are linked to lower phosphocreatine/ATP ratios in humans", CIRCULATION, vol. 97, 1998, pages 1716 - 1723 |
P.A. BOTTOMLEY: "Encyclopedia of Magnetic Resonance", 2009, JOHN WILEY, article "NMR Spectroscopy of the Human Heart" |
P.A. BOTTOMLEY: "Spatial localization in NMR spectroscopy in vivo", ANNAL NY ACAD SCI, vol. 508, 1987, pages 333 - 348 |
P.A. BOTTOMLEYC.J. HARDY: "Strategies and Protocols for Clinical 31P Research in the Heart and Brain", PHIL. TRANS. R. SOC. LOND. A, vol. 333, 1990, pages 531 - 544, XP009181094, DOI: 10.1098/rsta.1990.0181 |
P.A. BOTTOMLEYC.J. HARDYP.B. ROEMERR.G. WEISS: "Problems and expediencies in human 31P spectroscopy. The definition of localized volumes, dealing with saturation and the technique-dependence of quantification", NMR IN BIOMEDICINE, vol. 2, 1989, pages 284 - 289 |
R. LOFFLERR. SAUTERH. KOLEMA. HAASEM. VON KIENLIN: "Localized spectroscopy from anatomically matched compartments: improved sensitivity and localization for cardiac 31P MRS in humans", JOURNAL OF MAGNETIC RESONANCE, vol. 134, 1998, pages 287 - 299, XP004407377, DOI: 10.1006/jmre.1998.1497 |
R. ORDIDGEA. CONNELLYJ. LOHMAN: "Image-selected in vivo spectroscopy (ISIS). A new technique for spatially selective NMR spectroscopy", JOURNAL OF MAGNETIC RESONANCE, vol. 66, 1969, pages 283 - 294, XP023957448, DOI: 10.1016/0022-2364(86)90031-4 |
R.A. HORNC.R. JOHNSON: "Matrix Analysis", 1990, CAMBRIDGE UNIVERSITY PRESS |
R.E. GABRR. OUWERKERKP.A. BOTTOMELY: "Quantifying in vivo MR spectra with circles", JOURNAL OF MAGNETIC RESONANCE, vol. 179, 2006, pages 152 - 163 |
R.G. WEISSG. GERSTENBLITHP.A. BOTTOMLEY: "ATP flux through creatine kinase in the normal, stressed, and failing human heart", PROC NATL ACAD SCI USA, vol. 102, 2005, pages 808 - 813 |
R.G. WEISSP.A. BOTTOMLEYC.J. HARDYG. GERSTENBLITH: "Regional Myocardial Metabolism of High-Energy Phosphates during Isometric Exercise in Patients with Coronary Artery Disease", N ENGL J MED, vol. 323, 1990, pages 1593 - 1600 |
See also references of EP2688471A4 |
T.R. BROWNB.M. KINCAIDK. UGUIBIL: "NMR chemical shift imaging in three dimensions", PROC. NATL ACAD SCI USA, vol. 79, 1982, pages 3523 - 3526, XP000653160, DOI: 10.1073/pnas.79.11.3523 |
VON KIENLIN M ET AL., JMR, vol. 94, 1991, pages 268 - 287 |
X. HUD.N. LEVINP.C. LAUTERBURT. SPRAGGINS: "SLIM: spectral localization by imaging", MAGNETIC RESONANCE IN MEDICINE, vol. 8, 1988, pages 314 - 322 |
Y. ZHANGR.E. GABRM. SCHARH. ZHUP. BARKERR.G. WEISSP.A. BOTTOMLEY: "Dramatic speedup in 1D-, 2D- and 3D- MRS scan times with linear algebraic modeling (SLAM", PROCEEDINGS OF THE INTERNATIONAL SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE, IN PRESS, 2012 |
Z. DONGJ.H. HWANG: "Lipid signal extraction by SLIM: Application to 1H MR spectroscopic imaging of human calf muscles", MAGNETIC RESONANCE IN MEDICINE, vol. 55, 2006, pages 1447 - 1453 |
Z.P. LIANGF. BOADAR. CONSTABLEE. HAACKEP. LAUTERBURM. SMITH: "Constrained reconstruction methods in MR imaging", REV MAGN RESON MED, vol. 4, 1992, pages 67 - 185, XP009035614 |
Z.P. LIANGP.C. LAUTERBUR: "A generalized series approach to MR spectrocopic imaging", IEEE TRANSACTIONS MEDICAL IMAGING, vol. 10, 1991, pages 132 - 137, XP000234747, DOI: 10.1109/42.79470 |
Also Published As
Publication number | Publication date |
---|---|
JP2014511727A (en) | 2014-05-19 |
US20140015529A1 (en) | 2014-01-16 |
JP6070963B2 (en) | 2017-02-01 |
EP2688471B1 (en) | 2023-06-21 |
US10209330B2 (en) | 2019-02-19 |
WO2012129433A3 (en) | 2012-12-27 |
EP2688471A2 (en) | 2014-01-29 |
EP2688471A4 (en) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Magnetic resonance Spectroscopy with Linear Algebraic Modeling (SLAM) for higher speed and sensitivity | |
EP2688471B1 (en) | System and method of performing magnetic resonance spectroscopic imaging | |
Nelson | Analysis of volume MRI and MR spectroscopic imaging data for the evaluation of patients with brain tumors | |
Lee et al. | MR‐based conductivity imaging using multiple receiver coils | |
US7924003B2 (en) | Method for water-fat separation and T2* estimation in an MRI system employing multiple fat spectral peaks and fat spectrum self-calibration | |
Cruz et al. | 3D free‐breathing cardiac magnetic resonance fingerprinting | |
Zhu et al. | Coil combination methods for multi-channel hyperpolarized 13C imaging data from human studies | |
Kühn et al. | T1 bias in chemical shift‐encoded liver fat‐fraction: role of the flip angle | |
Löffler et al. | Localized spectroscopy from anatomically matched compartments: improved sensitivity and localization for Cardiac31P MRS in humans | |
Zhang et al. | Highly-accelerated quantitative 2D and 3D localized spectroscopy with linear algebraic modeling (SLAM) and sensitivity encoding | |
An et al. | Spectral localization by imaging using multielement receiver coils | |
Maril et al. | Strategies for shimming the breast | |
Cao et al. | Optimized multi‐axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole‐brain high‐isotropic‐resolution quantitative imaging | |
Yang et al. | Fat fraction bias correction using T1 estimates and flip angle mapping | |
Vareth et al. | A comparison of coil combination strategies in 3D multi‐channel MRSI reconstruction for patients with brain tumors | |
Moran et al. | High‐resolution 3D radial bSSFP with IDEAL | |
Rowland et al. | Whole brain P MRSI at 7T with a dual‐tuned receive array | |
Vemuri et al. | Coil sensitivity estimation for optimal SNR reconstruction and intensity inhomogeneity correction in phased array MR imaging | |
Gast et al. | Assessing muscle‐specific potassium concentrations in human lower leg using potassium magnetic resonance imaging | |
US20080146951A1 (en) | Regional cerebral volume flow using quantitative magnetic resonance angiography | |
Mirkes et al. | Triple‐quantum‐filtered sodium imaging at 9.4 Tesla | |
Zhao et al. | Quantitative mapping of total choline in healthy human breast using proton echo planar spectroscopic imaging (PEPSI) at 3 Tesla | |
Edden et al. | Proton MR spectroscopic imaging of the medulla and cervical spinal cord | |
Goelman et al. | Reducing voxel bleed in Hadamard‐encoded MRI and MRS | |
Sung et al. | B compensation in 3T cardiac imaging using short 2DRF pulses |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12760705 Country of ref document: EP Kind code of ref document: A2 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14006069 Country of ref document: US |
|
ENP | Entry into the national phase |
Ref document number: 2014501253 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012760705 Country of ref document: EP |