This application is a continuation-in-part of commonly assigned copending U.S. patent application Ser. No. 07/545,068 filed Jun. 29, 1990 entitled "Switched-Field MRI" and naming Messrs. Carlson, Crooks and Kaufman as joint inventors. The entirety of this related copending application is hereby incorporated by reference.
This invention relates to method and apparatus for performing nuclear magnetic resonance (NMR) relaxometry. In vivo magnetic resonance imaging (MRI) may be simultaneously employed so as to provide NMR relaxometry measurements for the voxels of two or three dimensional images of living tissue. In effect, NMR relaxometry is used in imaging sequences having an acceptable total imaging time.
The referenced related copending application describes MRI apparatus and method wherein the background magnetic strength Bo may be selectively switched to different magnitudes at different times in NMR pulse sequences for different effects. However, it does not explicitly discuss NMR relaxometry.
NMR relaxometry is per se known. For example, Lundbom et al in "Magnetic Field Dependence of 1/T1 of Human Brain Tumors," Investigative Radiology, November 1990, Vol. 25, pp 397-1205 describe NMR relaxometry providing plots of proton relaxation rates T1 versus proton Larmor frequency/magnetic field strength. As Lundbom et al have shown, different types of tissue may be distinguished (i.e., identified) as a function of such relaxometry measurements.
Koenig also describe "A Field-Cycling Relaxometer" in a chapter on Relaxometry of Tissue within the text of "NMR Spectroscopy of Cells & Organisms," Volume II, Raj K. Gupta, Editor, CRC Press, Boca Raton, FL (1987), pp 108-114. Koenig et al first equilibrate nuclei at either zero Gauss or very high Gauss (e.g., 1.2 T). Then the field is switched to a T1 measurement value and, after a "measurement time to allow exponential field change at this new field, the field is changed to a fixed NMR value at which a 90°-180° spin echo measurement is quickly made. Since this SE magnitude is proportional to the magnetization achieved during the selected measurement time, one can effectively plot the exponential change in nuclei magnetization (and thus discern T1) by repeating the experiment for different measurement times. Koenig et al describe 23 cycles for each T1 measurement at a given field strength. However, as noted, all RF signals are conveniently made at a fixed magnetic field value.
Traditional NMR relaxometry (e.g., such as that in Koenig) is performed on an entire sample (e.g., as contained in a test tube within laboratory NMR apparatus). Such is not conducive to in vivo measurements--such as in the context of magnetic resonance imaging (MRI).
Ramping or changing the background magnetic field strength (e.g., as might be suggested by Lundbom et al) is not possible in conventional permanent magnet MR systems and is, at best, time consuming in a resistive or superconductive solenoidal magnet MRI system (especially if the time to reach stable magnet operation is added to the ramping of the magnetic field strength itself).
Such a straight forward approach to in vivo NMR/MRI relaxometry also would require the RF transmitter/receiver frequency (ann all attendant RF circuits, coils, amplifiers, etc.) to be changed and/or retuned to the new NMR Larmor frequency associated with the new background field strength.
Koenig et al is at least a more interesting approach since all RF measurements are made at a fixed magnetic field strength--and hence at fixed RF frequency bands. However, the Koenig et al approach will cause nuclei magnetization to partly equilibrate during the NMR measurement phase. If only a small field change is involved, the change in equilibrium magnetization may also be very small. To enhance the change in magnetization, Koenig et al allows the nuclei magnetization to first come to equilibrium with either a very strong (e.g., 1.2 T) or very weak (e.g., zero) magnetic field. Furthermore, this prior "soak" field is re-applied before each of many (e.g., 23) data point measurement cycles (e.g., using many different measurement time intervals at each given measurement field value). Such a technique may be excellent for laboratory samples without imaging. However, it would take far too long for use in effective MR imaging procedures with in vivo patient anatomy.
Koenig et al use a technique of (1) applying a soak field, (2) applying a "relaxation rate field" for some amount of time and (3) doing NMR at a third field. The three steps are repeated using variable time delays in step (2).
Now, however, we have discovered a technique for using the field-switched MRI apparatus of our earlier above-referenced copending application so as to perform NMR relaxometry measurements at different effective background magnetic field strengths Bo --even in a two or three dimensional MRI context--in vivo, and still without the need to change RF operating frequencies from those employed in regular MRI procedures using a steady background magnetic field Bo. In effect, we dispense with step (1) of Koenig et al and apply step (2) with only one time delay. There are at least two exemplary ways to achieve this desired result.
1. One way is to first measure T1 at the base field of the magnet using any one of the common conventional techniques. Subsequent measurements at a switched field require solution of the Bloch equations for the evolution in a multistep, switched field. Finding T1 at the switched field involves solving an equation, an example is given below. This task can be programmed quite easily and solved numerically.
2. A second technique is to apply a 90 degree NMR RF pulse prior to the switched field. This is a distinct way of doing the measurement but may not turn out to be the most accepted technique. This is physically analogous to applying a zero "soak" field.
In our exemplary embodiment, the main magnet of a MRI system need not be changed or ramped. Rather, a switched additive/subtractive field can be raised or lowered relative to the main field and utilized only at times when RF signals are not being transmitted or received. In this way, the RF receiver/transmitter frequencies may be fixed to the usually expected frequencies and bandwidths and RF coils and other RF circuitry need not be unusually changed or manipulated.
Furthermore, in the second exemplary embodiment, instead of a long pre-"soak" period at an extremely high or extremely low magnetic field, we use a simple 90° NMR RF nutation pulse. This is roughly equivalent in results to applying a zero field for a very long time--but much easier and faster to implement.
We have also devised a technique for measuring T1 with just one application of the desired measurement field strength. While this may lower the effective signal-to-noise ratio (e.g., due to use of a reduced minimum size data set), the procedure is much more time-tolerable to an MRI patient.
In brief summary, the desired relaxometry background field strength is applied and maintained for a known duration Δt. It may be for a relatively long duration (e.g., on the order of the T1 parameter value being measured or longer). During this time no RF signals are being transmitted or received. However, in an immediately following time period--when the background field is at its usual nominal and fixed value Bo the effective T1 value for the NMR nuclei is measured. If the measurement time period is substantially shorter than the relaxometry field period, the T1 value associated with the desired relaxometry field strength value will predominate and produce the desired measurements providing T1 versus magnetic field strength data for in vivo tissue containing the NMR nuclei.
If the measurement time period is not substantially shorter than the relaxometry field period, then a penalty is attached to this multi-field strength approach. Namely, within a single relaxometry data acquisition cycle, the calculation of T1 is somewhat more complicated--since NMR relaxation actually occurs at two different field strengths during the MRI relaxometry sequence.
Furthermore, by using an appropriate NMR/MRI sequence such as explained below in one of the exemplary embodiments, the computation can be simplified. For example, in one exemplary embodiment of this invention, a 90° RF NMR nutation pulse is introduced just prior to switching "on" the switched background field. Thereafter, immediately after switching the supplemental background field "off" (and therefore returning to the same constant background field strength Bo), a 90°-180° RF NMR nutation pulse sequence is performed so as to generate a spin echo RF NMR response within a relatively short time (i.e., a short TE). In such case, the T1 relaxation at the switched "on" composite background field strength dominates the observed change in NMR signal strength. Two such NMR measurement sequences using different time intervals between the 90° RF nutation pulses of appropriate lengths (e.g., as related by integer multiples) permit the T1 relaxation parameter to be determined in a relatively straightforward manner (e.g., in accordance with the teaching of Ortendahl et al, "Analytical Tools for Magnetic Resonance Imaging", Radiology, Vol. 153, No. 2, November 1984, pp 479-488) for one field strength. Thereafter, if other imaging parameters are maintained constant, T1 can be determined in only a single further measurement cycle for each different field strength. That is, N relaxometry data points can be determined in only N+1 measurement cycles.
It is believed that this novel switched-field approach to NMR relaxometry provides a more practical way to achieve in vivo relaxometry in magnetic resonance imaging than heretofore. (For example, see an earlier in vivo MRI relaxation technique as described in Crooks et al, "Magnetic Resonance Imaging: Effects of Magnetic Field Strength," Radiology, Vol. 151, No. 1, April 1984, pp 127-133.)
Thus, compared to prior relaxometry approaches (even the non-imaging approaches), we have achieved considerable advantage in at least each of the following ways:
1. Use of an initial 90° NMR RF nutation pulse (in one of the exemplary embodiments) to more quickly initialize nuclei magnetization state;
2. acquisition of N relaxometry data points (T1,B) using only N+1 measurement cycles (in both exemplary embodiments); and
3. MRI of relaxometry data measured in such time efficient procedures (in both exemplary embodiments)--e.g., thus permitting interleaved imaging measurements so as to preserve registration of patient anatomy vis-a-vis the acquired relaxometry data for each imaged voxel.
These as well as other objects and advantages of this invention will be more completely understood and appreciated by careful study of the following detailed description of the presently preferred exemplary embodiments of this invention taken in conjunction with the accompanying drawings, of which:
FIG. 1 is a simplified general block diagram of a switched-field MRI system suitable for practicing the exemplary embodiment of this invention;
FIG. 2 is a timing diagram showing an exemplary sequence of RF nutation pulses switched magnetic fields, magnetic gradients, etc., for an exemplary embodiment of this invention;
FIG. 3 is a flowchart of an exemplary program segment that may be incorporated into the programmed MRI system control of FIG. 1 so as to practice the exemplary embodiment of this invention; and
FIG. 4 is a flowchart of another exemplary program segment that may be incorporated into the programmed MRI system control of FIG. 1 so as to practice the exemplary embodiment of this invention.
FIG. 1 is a simplified schematic depiction of a typical MRI system that has been adapted to practice this invention. For example, it includes the usual gantry system 10 in which an image volume is adapted to receive a living human subject 12. A suitable generator 14 of a substantially homogeneous background magnetic field Bo is provided. This may be a solenoidal cryogenic electromagnet (with the main background magnetic field oriented horizontally), a suitable permanent magnet structure (typically with the field oriented vertically) or the like.
An additional switched "background" magnetic field LB generator 16 is also provided as described in our above-referenced related copending parent application Ser. No. 07/545,068 so as to selectively alter the background magnetic field strength at selected times during MRI data gathering cycles. The usual RF receiving and/or transmitting coils 18 are also provided and closely RF coupled to the image volume. The usual gradient coil structures Gx, Gy, Gz, 20 are also provided.
The controllable portions of the gantry system 10 are controlled by a suitable MRI system control 30 having a system of computers, input/output drivers, RF receiver/transmitter circuits, etc., as will be known to those skilled in this art. Typically, the interface for operator control may be an operator console 40 having the usual keyboard and screen for inputting control commands for viewing MRI data, parameters, images, etc. As will be appreciated, other peripheral devices are typically provided for recording MRI data, images, etc.
The MRI system of FIG. 1 is, in one exemplary embodiment, programmed so as to effect the exemplary MRI relaxometry cycles depicted at FIG. 2. Here, immediately after a 90° NMR RF nutation pulse, the switched field is turned "on" to effect a composite background magnetic field strength of a selected magnitude. This switched field strength is maintained for a time on the order of the value of the NMR T1 parameter to be measured (e.g., 100-300 milliseconds). In other words, the NMR nuclei are preferably permitted to approach saturation magnetization in accordance with the switched field strength magnitude.
Immediately thereafter, 90° and 180° NMR RF nutation pulses are imposed in relatively quick succession so as to elicit an NMR RF spin echo signal response still having a characteristic magnitude related to the earlier switched background magnetic field strength. The time to echo or TE parameter is preferably made relatively short compared to the expected T1 value. For example, TE is, in the exemplary embodiment, maintained less than about 20 to 30 milliseconds.
The Gz, Gy and Gx gradient field pulses depicted in FIG. 2 are typical of those used for two-dimensional Fourier Transform (2DFT) MRI. In this manner, spin echo response magnitudes for each voxel of a two-dimensional image "slice" volume through the patient 12 may be obtained. As should be apparent, three-dimensional Fourier Transform MRI, projection reconstruction MRI or other types of MRI measurement cycles could also be used. Furthermore, in vivo relaxometry measurements could be made on selected "slices" as a whole or other portions within the image volume or engross for the entire image volume without any imaging per se.
The elicited spin echo (SE) RF responses are proportional (assuming all other imaging parameters to remain constant) to the following function of the time duration Δt and the T1 parameter associated with the switched background magnetic field:
SE∝1-exp(-Δt/T1) (Equation 1)
At one particular field strength, if TE is kept constant, then two different measurements (at different TR values) are needed to solve for T1 since there are two unknowns. One unknown is T1, the other is a multiplicative combination of equilibrium magnetization density, receiver efficiency and gain, and an exponential factor of T2 decay. If the field is now switched (e.g., upwards by 50%), the "other" factor (involving equilibrium magnetization density, receiver efficiency, gain and T2 etc) goes linearly with the field strength. Thus, after the first data point is measured using as few as two measurement cycles, there remains effectively only one unknown for each different field strength to be probed. Accordingly, by measuring SE signal using two different Δt values at one field strength, then measuring one point with a Δt on the order of T1 at each new field strength (TE being held constant) one can calculate the desired T1 value at each new probed magnetic field strength using but a single further measurement cycle.
Since there is no read-out actually taken during the time the switched field is present, there is less requirement for the switched field to be absolutely constant, homogeneous, etc. of course to the extent the switched field can be made constant and homogeneous, that would be desirable.
To implement the saturation recovery NMR relaxometry sequence depicted in FIG. 2, the NMR system control 30 of FIG. 1 can be programmed as generally depicted in the flowchart of FIG. 3. Here, upon entry at 300 to this program segment, initialization is effected at 302 (e.g., so as to initialize the values stored in the switched field magnitude register ΔB, the phase encoding register PE and the switched field duration register Δt).
Thereafter, a 90° RF NMR nutation pulse is transmitted at the usual Larmor frequency associated with the nominal background Bo at 304. Immediately thereafter, a switched field is switched on for the ΔAt duration at 306. Immediately, thereafter, a 90°-180° NMR sequence is performed at 308 so as to capture suitable spin echo data for the existing phase encoded Gy value. A test is then made at 310 (if magnetic resonance imaging is to be employed) so as to see if the phase encoding register has yet reached the maximum desired phase encoded value. If not, the phase encoding register PE is incremented at 312 and the cycle is repeated by returning control to block 306.
Eventually, all of the desired phase encoded spin echo responses will have been captured (e.g., 256 of them if MRI is being performed or zero if no imaging is to be performed) and the phase encoding register PE will be reset at 314 in readiness for a subsequent MRI cycle. A test is then made at 316 to see if the duration register Δt is yet at its maximum value. If not, that register is incremented at 318 (e.g., by an integer multiple) and control is again returned to block 306 for another MRI cycle at the next Δt increment. As will be understood, in the exemplary embodiment, the maximum Δt value will be reset at 318 so that after the first relaxometry data point has been captured only one At value is used for each new ΔB measurement field.
After the desired number of MRI cycles for different Δt increments (e.g., two on the first ΔB value, one thereafter), the Δt register is reset at block 320 in readiness for a future data gathering cycle and a test is made at 322 to see if the switched field magnitude register ΔB is yet at its desired maximum value. If not, the switched field is incremented at 324 and control is again returned back to block 306 for yet another MRI data gathering cycle. Finally, after data is gathered for all of the desired ΔB increments, the ΔB register is reset at 326.
At that time (or subsequently if desired) the desired relaxometry measurements of T1 versus ΔB may be determined for each MRI voxel (e.g., the first using two Δt values in accordance with the teachings of the related Ortendahl et al teachings and subsequent points using only single Δt values per point as noted above). Ultimately, the program segment of FIG. 3 is then exited at 330 and control is returned to other conventional programming within the MRI system control 30.
As should now be apparent, individual relaxometry measurement cycles (as depicted in FIG. 2) can also be interleaved between conventional MRI data acquisition cycles to better preserve data/image registration of the patient.
Another exemplary embodiment (as depicted in FIG. 4) avoids the use of a 90° RF saturation pulse prior to each change in the background field Bo. Here, use is made of the well-known relationship between NMR nuclei magnetization (and the related NMR spin echo (SE) signal magnitude) changes with respect to the well-known NMR parameters TE, TR, T1 and T2:
SE∝[exp(-TE/T2)][1+exp(-TR/T1)-2exp{-(TR-TE/2)/T1}] (Equation 2)
First, by keeping TE constant, the TE and T2 dependence can be reduced to just another constant factor of proportionality:
SE∝1+exp(-TR/T1)-2exp[-(TR-TE/2)/T1] (Equation 3)
Then, the two NMR spin echo experiments (e.g., SE=985 and 633 signal strength units, respectively) depicted at 400, 402 at the constant base background field Bo (e.g., 650 Gauss) and constant TE (e.g., 30 ms) but different TR values (e.g., 300 ms and 150 ms) permits determinations of T1 (e.g., 164 ms) at step 404 using conventional exponential curve fitting techniques (e.g., empirical trial and error using different known exponential curves).
Once T1 at the base field Bo is thus derived, the asymptote (e.g., 1216 signal strength units) is similarly ascertainable at 406 for an infinite TR spin echo experiment using the same other NMR parameters.
Thereafter, as depicted at steps 408 and 410, for each changed background field strength ΔB pulse of Δt duration, one immediately following NMR spin echo experiment (using the same constant TE at the same base field Bo) suffices to determine T1 at the changed background field value:
T1=-Δt/ln[1-(SE.sub.m -SE.sub.x)/(M·SE.sub.a -SE.sub.x)](Equation 4)
where
T1=the T1 NMR parameter value at the changed background field strength
Δt=the time duration of the changed field pulse
SEm =the spin echo signal magnitude measured immediately following the changed field pulse
SEx =the spin echo signal magnitude which would result at hypothetical TR interval which would result if Δt is reduced to zero
M=the ratio of the changed magnetic field to the background magnetic field
SEa =the asymptotic spin echo signal magnitude for the background magnetic field strength
In the example depicted at steps 408 and 410, the measured spin echo magnitude was 1186 of the same arbitrary signal strength units used in the other FIG. 4 measurements. Thus, T1 at the boosted 1300 Gauss field strength is determined to be 410 ms. Repetition of steps 408 and 410 for other field strengths will yield the respectively corresponding T1 values as depicted at 412 in FIG. 4.
To summarize, first, with a TR=300 ms, TE=30 ms, constant field sequence one measures a signal of 985 (arbitrary units). Repeating this with TR 150 ms, the same echo time TE gives a measured signal of 633 of the same arbitrary units. One can now solve the T1 (e.g., 164 ms) at the base field strength and the signal (e.g., 1216) we would expect to measure in an infinite TR spin echo experiment at the same constant TE. Now, with one more measurement using a TR=300 ms, using a field boost of 650 Gauss for 150 ms produces a measured signal of 1186. We can use this one piece of data to get the new T1. After 150 ms at the lower field, if the NMR spin echo experiment were repeated at that TR (i.e., TR=150 ms), then the SE signal and/or nuclei magnetization would be 633 units. However, in reality in this experiment, in the remaining 150 ms, the magnetization grows towards the new equilibrium asymptote value of 2432 (this is 2 times the old asymptote; the factor of two comes from the 2× boost in the field) and actually is measured as 1186. This is enough to give T1 at the higher field ##EQU1##
As will be apparent to those in the art, other cases such as for different times for high and low field strengths and multi-slice sequences can be similarly derived. It will also be apparent that the process of FIG. 4 can be viewed as a flowchart description of a computer program subroutine module.
While only two exemplary embodiments have been described in detail, those skilled in the art will realize that many variations and modifications may be made in such exemplary embodiments while yet retaining many of the novel features and advantages of this invention. Accordingly, all such variations and modifications are intended to be included within the scope of the appended claims.