Method, apparatus, and article for frequency-overlapped 3-D multispectral magnetic resonance images
10061007 ยท 2018-08-28
Assignee
Inventors
Cpc classification
International classification
Abstract
A method for acquiring 3D multispectral MRI of a target includes scanning a spectrum of spectral windows with an MRI scanner, wherein each spectral window of the spectrum defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart by substantially uniform frequency offsets such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass.
Claims
1. A method for acquiring 3D multispectral MRI of a target, comprising: scanning a spectrum of spectral windows with an MRI scanner, wherein each spectral window of the spectrum defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart by substantially uniform frequency offsets such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass.
2. The method as claimed in claim 1, wherein adjacent spectral windows overlap by between about 0.25 standard deviation and about 0.75 standard deviation.
3. The method as claimed in claim 1, wherein adjacent spectral windows overlap by between about 0.4 standard deviation and about 0.6 standard deviation.
4. The method as claimed in claim 1, wherein adjacent spectral windows overlap by about 0.5 standard deviation.
5. The method as claimed in claim 1, wherein adjacent scan frequencies are spaced apart by about 1 standard deviation of their respective spectral windows.
6. The method as claimed in claim 1, further comprising: superimposing images produced by scanning a plurality of the spectrum of spectral windows.
7. The method as claimed in claim 1, wherein each spectral window of each pass has a scan frequency adjacent to the scan frequency of the corresponding spectral window within the next pass, and within each pass, consecutive spectral windows have maximally separated scan frequencies.
8. The method as claimed in claim 7, wherein the number of passes is chosen as a divisor of the total number of scan frequencies, according to a minimum echo time required for scanning each spectral window.
9. The method as claimed in claim 1, wherein a controller of the MRI scanner is configured to scan the spectral windows in a sequence of passes according to
f(b=1,p)f.sub.maxf.sub.sep(i1),
and
f(b=2:Nb,p)=f(b1,p)2*f.sub.sepN.sub.int for frequency offsets of passes indexed from p=1:.sub.int, and spectral windows indexed from b=1:N.sub.b in each pass, where f.sub.max is a maximum frequency offset within the acquisition spectrum, f.sub.sep is a frequency offset between adjacent scan frequencies, and i is the sequence index of the pass p.
10. An apparatus configured to perform 3-D multispectral MRI, comprising: a magnet assembly; and an MRI controller that is configured to operate the magnet assembly for scanning a spectrum of spectral windows, wherein each spectral window defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass.
11. The apparatus as claimed in claim 10, wherein adjacent spectral windows overlap by between about 0.25 standard deviation and about 0.75 standard deviation.
12. The apparatus as claimed in claim 10, wherein adjacent spectral windows overlap by between about 0.4 standard deviation and about 0.6 standard deviation.
13. The apparatus as claimed in claim 10, wherein adjacent spectral windows overlap by about 0.5 standard deviation.
14. The apparatus as claimed in claim 10, wherein adjacent scan frequencies are spaced apart by about 1 standard deviation of their respective spectral windows.
15. The apparatus as claimed in claim 10, wherein the MRI controller is further configured to superimpose images produced by a plurality of the spectrum of spectral windows.
16. An article of computer readable media encoded with a magnetic resonance image by a process comprising: scanning a spectrum of spectral windows with an MRI scanner, wherein each spectral window defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass.
17. The article as claimed in claim 16, wherein the process further includes superimposing images produced by a plurality of the spectrum of spectral windows.
18. The article as claimed in claim 16, wherein spectral windows are scanned in a sequence of passes according to
f(b=1,p)f.sub.maxf.sub.sep(i1),
and
f(b=2:Nb,p)=f(b1,p)2*f.sub.sepN.sub.int for frequency offsets of passes indexed from p=1: N.sub.int, and spectral windows indexed from b=1:N.sub.b in each pass, where f.sub.max is a maximum frequency offset within the acquisition spectrum, a f.sub.sep is a frequency offset between adjacent spectral windows, and i is the sequence index of the pass p.
19. The article as claimed in claim 18, wherein N.sub.int is chosen as a divisor of the total number of spectral windows, according to minimum echo time t.sub.min required for scanning each spectral window.
20. The article as claimed in claim 16, wherein each spectral window of each pass has a scan frequency adjacent to the scan frequency of the corresponding spectral window within the next pass, and within each pass, consecutive spectral windows have maximally separated scan frequencies.
Description
DRAWINGS
(1) The present invention will be better understood from reading the following description of non-limiting embodiments, with reference to the attached drawings, wherein below:
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION
(8) Reference will be made below in detail to exemplary embodiments of the invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference characters used throughout the drawings refer to the same or like parts, without duplicative description. Exemplary embodiments of the present invention are described with respect to 3D multispectral magnetic resonance imaging (3D MS MRI).
(9)
(10) The system control 32 includes a set of modules connected together by a backplane 32a. These include a CPU module 36 as well as a pulse generator module 38, which connects to the operator console 12 through a serial link 40. It is through link 40 that the system control 32 receives commands from the operator to indicate the scan sequence that is to be performed. The pulse generator module 38 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 38 connects to several components that are operated by the MRI controller 33, including a set of gradient amplifiers 42, a physiological acquisition controller 44, and a scan room interface circuit 46.
(11) In addition to the above, the pulse generator module 38 operates the gradient amplifiers 42 to achieve desired timing and shape of the gradient pulses that are produced during the scan. The module 38 also receives patient data from the physiological acquisition controller 44, which, in turn, receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes. The pulse generator module 38 receives signals from various sensors associated with the condition of the patient and the magnet system via the scan room interface circuit 46. The MRI controller 33 utilizes the scan room interface circuit 46 to command a patient positioning system 48 to move the patient or client C to a desired position for the scan.
(12) The gradient waveforms produced by the pulse generator module 38 are applied to the gradient amplifier system 42 having Gx, Gy, and Gz amplifiers. Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly, generally designated 50, to produce the magnetic field gradients used for spatially encoding acquired signals. The gradient coil assembly 50 forms part of a magnet assembly 52, which also includes a polarizing magnet 54 and a whole-body RF coil 56. In an embodiment of the invention, RF coil 56 is a multi-channel coil. A transceiver module 58 in the system control 32 produces pulses that are amplified by an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive switch 62. The resulting signals emitted by the excited nuclei in the patient may be sensed by the same RF coil 56 and coupled through the transmit/receive switch 62 to a preamplifier 64. The amplified MR signals are demodulated, filtered, and digitized in the receiver section of the transceiver 58. The transmit/receive switch 62 is controlled by a signal from the pulse generator module 32 to electrically connect the RF amplifier 60 to the coil 56 during the transmit mode and to connect the preamplifier 64 to the coil 56 during the receive mode. The transmit/receive switch 62 can also enable a separate RF coil (for example, a surface coil) to be used in either transmit mode or receive mode.
(13) The MR signals picked up by the multi-channel RF coil 56 are digitized by the transceiver module 58 and are transferred to a memory module 66, or other computer readable media, via the system control 32. Computer readable media may include, for example, structures configured so that electrical, optical, or magnetic states may be fixed in a manner perceptible and reproducible by a conventional computer: e.g., text or images printed to paper or displayed on a screen, optical discs, or other optical storage media; flash memory, EEPROM, SDRAM, or other electrical storage media; floppy or other magnetic discs, magnetic tape, or other magnetic storage media.
(14) A scan is complete when an array of raw k-space data has been acquired in the computer readable media 66. This raw k-space data is rearranged into separate k-space data arrays for each image to be reconstructed, and each of these is input to an array processor 68 which operates to Fourier transform the data into an array of image data. This image data is conveyed through the serial link 34 to the computer system 20 where it is stored in memory. In response to commands received from the operator console 12, this image data may be archived in long-term storage or it may be further processed by the image processor 22 and conveyed to the operator console 12 and presented on the display 16.
(15) When scanning a target (client C) who has a paramagnetic inclusion, e.g., a bone screw or joint replacement, the implant distorts the gradient waveforms Gx, Gy, Gz in such a way that the frequency response at any given point within the target may be offset from an expected zero frequency, as shown in
(16) Accordingly, in an exemplary aspect of the invention, as shown in
(17) Referring to
(18) At a step 212, the MRI controller 33 receives from the operator console 12 a user-selectable value TR that prescribes the interval duration desired in order to obtain useful T1/T2 contrasts. At subsequent steps, the MRI controller 33 then establishes a sequence for scanning the spectrum of spectral windows in a plurality of scan intervals p.
(19) In order to appropriately schedule the scan intervals p, at step 214 the MRI controller 33 implements PSD calculations, familiar to the skilled worker, utilizing systemic hardware and specific-absorption-rate (SAR) limits, to compute for each spectral window 300 a minimum echo train listening time t.sub.min.
(20) At step 216 the MRI controller 33 calculates a local variable
mMBPP=min(N.sub.T/2,TR/t.sub.min).
Then at step 218, the MRI controller 33 computes a preliminary number N.sub.int of scan intervals p, based on divisors of N.sub.T, e.g., for N.sub.T=24, N.sub.int=f(mMBPP): mMBPP<6, 6; 6mMBPP<8, 4; 8mMBPP<12, 3; 12mMBPP, 2
and also computes a number N.sub.b of bins to be included in each pass p, such that N.sub.b is the largest number that evenly divides the total number of spectral windows N.sub.T among N.sub.int passes or intervals.
(21) The MRI controller 33 then at step 220 determines across all intervals p a minimum interval duration TR.sub.min in light of all system parameters, i.e.
TR.sub.min=max(p=1 . . . N.sub.int).sub.p(1,N.sub.b)t.sub.min(i).
(22) At step 222 the MRI controller 33 compares TR to the PSD calculated TR.sub.min. In case TR.sub.min>TR, then at step 224 the MRI controller 33 modifies the bin construction to reduce N.sub.b, i.e. increase N.sub.int for more intervals with fewer bins per interval so as to reduce TR.sub.min. Steps 220-224 can iterate up to N.sub.int(max)=GCD(N.sub.T), e.g., as many as eight intervals for an exemplary embodiment with N.sub.T=24.
(23) Bin sorting algorithms, according to embodiments of the invention, schedule adjacent scan frequencies to be consecutively scanned at time intervals according to their respective repetition times TR, and schedule non-adjacent scan frequencies to be scanned according to their respective acquisition times t.sub.min, so that the nearest neighbor scan frequencies within each interval are scanned at a maximum temporal spacing within the interval. Generally, an aspect of the invention is that frequency-adjacent bins are scanned at TR periodicity in order to build-up T1 signal by aliasing, but within TR intervals, frequency overlap is avoided or minimized among time-adjacent bins. For example, as shown in
(24) The above algorithm 228 also can be described via a recursion, which provides a deterministic bin arrangement mechanism for frequency offsets of passes or intervals indexed from p=1:N.sub.int, and spectral windows indexed from b=1:N.sub.b in each interval. The initial bin offsets in each pass (p) are given by
f(1,p)=f.sub.maxf.sub.sep(i1),
where f.sub.max is a maximum frequency offset within the acquisition spectrum, f.sub.sep is a frequency offset between adjacent spectral windows, and i is the sequence index of the pass p. For example, for bins in a range b=2:Nb, frequency offsets then are computed according to
f(b,p)=f(b1,p)2*f.sub.sepN.sub.int.
(25) This recursion relation must also enforce periodic boundary conditions, which can be expressed through the condition:
if [f(b,p)<(f.sub.maxf.sub.sepN.sub.T)] then [f(b,p)=f(1,p)N.sub.int].
(26) It should be noted that the exemplary bin-sorting algorithm inherently restricts frequency of bin set aliasing to only a single alias per interval. For longer TR periods needed for PD and STIR contrasts, the additional time for longitudinal relaxation recovery accommodates overlap between n+2 bins that have maximum temporal separation from one another within a single interval, mitigating crosstalk artifacts that would otherwise emerge. This approximation allows for shorter ETL values to be used for PD and STIR contrasts, which moderates blurring from short T2 tissues.
(27) To minimize residual spatial distortion artifacts, it is helpful to select spectral windows that have continuously-differentiable frequency distributions (in some embodiments, Gaussian distributions) and that overlap substantially uniformly, for example, by between about 0.25 and about 0.75, as shown in
(28) For longer TR periods, which utilize more spectral windows in each interval, only n and n+1 bins are separated into separate intervals. For shorter TR periods needed for enhanced T1 contrast, the n, n+1, and n+2 spectral windows all should be acquired in separate intervals (i.e. N.sub.int3). This constraint can easily be implemented by setting the minimum ETL for the sequence to be roughly 8. When short TR periods are prescribed under such a constraint, the presented algorithm automatically sorts the bins into a least 3 intervals.
(29) In summary, the presented algorithm allows for a prospective automated determination of optimal interleaved bin arrangements that minimize crosstalk between overlapping bins for any prescribed TR period. In particular, this feature enables flexibility of TR prescription, which adds T1 contrast capabilities to implementations of heavily overlapping and interleaved 3D MS MRI.
(30) Having established an appropriate sequence of intervals, at step 230 the MRI controller 33 operates the patient positioning system 48 and the magnet assembly 52 in order to perform MRI of the target, according to the established sequence. The MRI controller 33 then additively superimposes the plurality of image data arrays generated from the plurality of spectral windows, thereby generating the 3-D multispectral magnetic resonance image 200 having enhanced T1 contrast relative to conventional modes of frequency-overlapped 3-D multispectral MRI.
(31) In certain embodiments, the invention provides a method for acquiring 3D multispectral MRI of a target. The inventive method includes scanning a spectrum of spectral windows with an MRI scanner, wherein each spectral window of the spectrum defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart by substantially uniform frequency offsets such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass. For example, adjacent spectral windows may overlap by between about 0.25 and about 0.75; by between about 0.4 and about 0.6; or by about 0.5. Adjacent scan frequencies may be spaced apart by about 1 of their respective spectral windows. The method may also include superimposing images produced by scanning a plurality of the spectrum of spectral windows. For example, each spectral window of each pass may have a scan frequency adjacent to the scan frequency of the corresponding spectral window within the next pass, and within each pass, consecutive spectral windows may have maximally separated scan frequencies. The number of passes may be chosen as a divisor of the total number of spectral windows, according to a minimum echo time required for scanning each spectral window. Particular embodiments may further include configuring a controller of the MRI scanner to scan the spectral windows in a sequence of passes according to
f(b=1,p)=fmaxfsep(i1),
and
f(b=2:Nb,p)=f(b1,p)2*fsepNint
for frequency offsets of passes indexed from p=1:Nint, and spectral windows indexed from b=1:Nb in each pass, where fmax is a maximum frequency offset within the acquisition spectrum, fsep is a frequency offset between adjacent scan frequencies, and i is the sequence index of the pass p.
(32) Other embodiments of the invention provide an apparatus configured to perform 3-D multispectral MRI. The apparatus includes a magnet assembly, and an MRI controller that is configured to operate the magnet assembly for scanning a spectrum of spectral windows, wherein each spectral window defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass. For example, adjacent spectral windows may overlap by between about 0.25 and about 0.75; by between about 0.4 and about 0.6; or by about 0.5. Adjacent scan frequencies may be spaced apart by about 1 of their respective spectral windows. The MRI controller may also be configured to superimpose images produced by scanning a plurality of the spectrum of spectral windows.
(33) In other embodiments, an article of computer readable media is encoded with a magnetic resonance image by a process that includes scanning a spectrum of spectral windows with an MRI scanner, wherein each spectral window defines a continuously-differentiable distribution of frequencies around a scan frequency and adjacent scan frequencies are spaced apart such that adjacent spectral windows substantially uniformly overlap, wherein selected adjacent spectral windows are scanned in consecutive passes, and nearest neighbor spectral windows within each pass are scanned at a maximum temporal spacing within the pass. The process may also include superimposing images produced by scanning a plurality of the spectrum of spectral windows. For example, the spectral windows may be scanned in a sequence of passes according to
f(b=1,p)=fmaxfsep(i1),
and
f(b=2:Nb,p)=f(b1,p)2*fsepNint
for frequency offsets of passes indexed from p=1:Nint, and spectral windows indexed from b=1:Nb in each pass, where fmax is a maximum frequency offset within the acquisition spectrum, fsep is a frequency offset between adjacent spectral windows, and i is the sequence index of the pass p. N.sub.int may be chosen as a divisor of the total number of spectral windows, according to minimum echo time t.sub.min required for scanning each spectral window. In certain embodiments, each spectral window of each pass may have a scan frequency adjacent to the scan frequency of the corresponding spectral window within the next pass, and within each pass, consecutive spectral windows may have maximally separated scan frequencies.
(34) It is to be understood that the above description is intended to be illustrative, and not restrictive. For example, the above-described embodiments (and/or aspects thereof) may be used in combination with each other. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from its scope. While the dimensions and types of materials described herein are intended to define the parameters of the invention, they are by no means limiting and are exemplary embodiments. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. In the appended claims, the terms including and in which are used as the plain-English equivalents of the respective terms comprising and wherein. Moreover, in the following claims, terms such as first, second, third, upper, lower, bottom, top, etc. are used merely as labels, and are not intended to impose numerical or positional requirements on their objects. Further, the limitations of the following claims are not written in means-plus-function format and are not intended to be interpreted based on 35 U.S.C. 112, sixth paragraph, unless and until such claim limitations expressly use the phrase means for followed by a statement of function void of further structure.
(35) This written description uses examples to disclose several embodiments of the invention, including the best mode, and also to enable one of ordinary skill in the art to practice embodiments of the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to one of ordinary skill in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.
(36) As used herein, an element or step recited in the singular and proceeded with the word a or an should be understood as not excluding plural of the elements or steps, unless such exclusion is explicitly stated. Furthermore, references to one embodiment of the present invention are not to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. Moreover, unless explicitly stated to the contrary, embodiments comprising, including, or having an element or a plurality of elements having a particular property may include additional such elements not having that property.
(37) Since certain changes may be made in the above-described method, apparatus, or article, without departing from the spirit and scope of the invention herein involved, it is intended that all of the subject matter of the above description or shown in the accompanying drawings shall be interpreted merely as examples illustrating the inventive concept herein and shall not be construed as limiting the invention.