Magnetic field and eddy current measuring method

Abstract: A sample (80) is disposed in an examination region (10) through which a strong, uniform magnetic field is generated. A gradient pulse (24, 82) is created across the magnetic field. After the gradient field has been turned off, a rapid succession of radio frequency pulses (52, 84) are applied. Each radio frequency pulse is followed by a free induction decay signal (54, 86). For each free induction decay signal, the evolution of phase vs. time is determined (60) because the slope of the changing phase vs. changing time curve is proportional to resonance frequency (66) and the resonance frequency is proportional to magnetic field strength (68). For each free induction decay signal, the magnetic field strength and the time since the determination of the gradient field pulse are determined and stored (70). From the plurality of free induction decay signals, a decaying curve (28) representing eddy current induced gradient field strength vs. time is produced.

Description: This invention relates to a magnetic field and eddy current measuring method. More especially the invention relates to such methods for use in magnetic resonance apparatus. It finds particular application in conjunction with measuring and compensating for eddy currents in magnetic resonance imagers and will be described with particular reference thereto. It is to be appreciated, however, that the invention may be utilized to measure and compensate for eddy current induced fields in other applications such as diffusion or flow studies or magnetic resonance spectroscopy. In magnetic resonance imaging, a strong substantially uniform magnetic field is generated longitudinally through an examination region. The magnetization vector of dipoles of a subject disposed in the examination region preferentially aligns with the uniform field. Radio frequency excitation pulses are supplied to cause the magnetization vectors to precess about the uniform field. Additional radio frequency pulses and magnetic field gradient pulses are applied to manipulate the precessing magnetization vector to create magnetic resonance signals, such as echo signals. Electrical current pulses are applied to the windings of gradient field magnet coils adjacent the examination region to create the gradient field pulses. A profile, or particular temporal dependence, is selected for the current pulse in accordance with the profile of the gradient magnetic field to be applied. Commonly, the current pulses strive to approximate a square wave, trapezoid, or other ideal gradient pulse profile. Inherently, the profile of the gradient magnetic field pulse does not match the profile of the electrical current pulse. A changing magnetic field induces eddy currents in adjacent conductive structures. The rising field at the leading edge of each gradient field pulse induces eddy currents that superimpose eddy current magnetic field components on the gradient pulse. The falling field at the trailing end of each pulse induces like, opposite polarity eddy currents that cause analogous eddy magnetic field components after the pulse. Thus, the eddy currents add unwanted eddy components to the gradient magnetic field pulse. The effect of the eddy current varies with the amount and conductivity of the material in which the eddy current is induced and the proximity of the material to the gradient coil. In order to improve the image quality, the shape of the electrical current pulse is commonly altered such that the magnetic field produced by the sum of the current pulse and the eddy currents approximates the desired gradient magnetic field pulse profile. Commonly, a current pulse correction or pre-emphasis circuit includes a plurality of filters whose characteristic frequencies are adjustable and an amplifier with an adjustable gain associated with each filter. The filter frequency and amplification factors are adjusted to add current components whose frequencies and amplitudes effectively cancel the induced eddy current fields. Asymmetric eddy currents are commonly dealt with by physical adjustments, such as gradient tube centering. Such physical adjustments compromise the pre-emphasis correction. The accuracy of the eddy current compensation is affected by the accuracy with which the eddy currents or eddy fields can be measured. Various measurement techniques have been developed. One eddy current measurement and compensation technique which is described in U.S. patent application serial number 118,865, utilizes a search coil and an integrator. Any change in the magnetic field strength perpendicular to a plane of the coil induces a voltage in the coil which is proportional to the change in magnetic field per unit time. The voltage wave form is integrated to yield a gradient wave form which is digitized for analysis. Although the search coil technique works well, it has some drawbacks. First, the search coil is sensitive to all magnetic field changes orthogonal to its plane, not just those changes along the magnet axes or z direction which are important in imaging. The coil is sensitive to motion, such as vibration within the main magnetic field. In order to reduce the vibration induced voltages in high field magnets, the main, uniform magnetic field may be ramped down prior to using the search coil for medium and long time constant analysis. Ramping down or turning off the main field has been found to alter reversibly the eddy current characteristics of a system. Second, the integrator circuitry tends to drift particularly during long time constant measurements. This drift or instability adversely affects the resultant measurement. Third, the search coil measures the average gradient over its area. With the relatively large search coils currently utilized, the average gradient field over the corresponding relatively large area is measured. Finer control over localization would be desirable, particularly for accurately plotting asymmetric gradient magnetic fields. The resonant frequency of a given sample is directly proportional to the magnetic field strength. There are various NMR methods that in some way use the samples' resonant frequency to measure time-dependent changes in the magnetic field strength. The free induction decay signal (FID) from a well-defined cylinder can be used to measure a constant gradient. (J.S. Murday, "Measurement of Magnetic Field Gradient by its Effect on the NMR Free Induction Decay" J. Mag. Res., Vol. 10, pp. 111-120, 1973). The Murday method has been extended to measure short time constant eddy currents by observing the phase response of a free induction decay signal from a small, uniform sample (E. Yamamoto and H. Khono, "Gradient Time-Shape Measurement by NMR" J. Phys. E: Sci. Instrum. , Vol. 19, pp. 708-711, 1986). However, spin-spin relaxation renders the Yamamoto and Khono method not applicable to medium and long time constant measurements. U.S. Patent No. 4,698,591 to G.H. Glover and N.J. Pelc extend this methodology to measure medium and long time constant eddy current effects. One FID signal is acquired from each scan sequence and its phase evolution is calculated. The absolute phase behavior of an entire series of FIDs acquired at successive delays after the gradient is forced to be continuous in order to construct a gradient response curve. Like Yamamoto and Khono, the Clover and Pelc technique is also limited by spin-spin relaxation such that a continuous series of overlapping FIDs must be acquired to map medium or long term effects. A further weakness of the Glover and Pelc technique is that any absolute phase errors accumulate in each successive FID. According to the present invention there is provided a method of measuring a magnetic field caused by eddy currents comprising disposing a sample in a main magnetic field in an examination region; causing a preselected gradient across the magnetic field in the examination region; removing the gradient; immediately following removal of the gradient, inducing a series of free induction decay signals from the sample; monitoring each free induction decay signal and determining its frequency (fr) which frequency (fr) is proportional to magnetic field strength (H), whereby the magnetic field strength (H) at a sequence of times after removal of the gradient is determined and wherein the monitoring step includes measuring signal phase () of the free induction decay signal vs. time (t) and wherein the frequency (fr) determining step includes determining a slope (d/dt) of the measured phase () vs. time (t). One advantage of the present invention is that it measures only z magnetization. It is relatively insensitive to vibration. Another advantage of the present invention is that it is particularly adept at measuring long time constants. Each point of a derived curve of the frequency vs. time is derived independently from an FID measurement relative to an RF synthesizer that has inherently greater stability than an integrator. Another advantage of the present invention is that it achieves good localization. Small spatial variations in magnetic field can be accurately measured. One apparatus and various methods in accordance with the invention will now be described, by way of example, with reference to the accompanying drawings in which:- Figure 1 is a diagrammatic illustration of the apparatus; Figure 2 illustrates a gradient pulse and free induction decay collection sequence including a typical eddy current induced gradient field response and correction; Figure 3 illustrates an alternate embodiment to the sequence of Figure 2 in which the RF pulse spaces are non-constant; Figure 4 illustrates a resonance frequency determining scheme; Figure 5 illustrates a technique for measuring spatial linearity of gradient magnetic fields in accordance with the present invention; Figure 6 illustrates an alternate embodiment for measuring temporal linearity; and, Figure 7 illustrates an alternate embodiment for measuring response linearity. With reference to Figure 1, the apparatus includes a main magnetic field means A for generating a substantially uniform, strong, magnetic field longitudinally through an examination or image region 10. The main magnetic field means includes a plurality of superconducting, permanent, or resistive main field magnets 12 which are operated under the control of a magnetic field control means and power supply 14 . Various electrical and mechanical shimming techniques are provided, as are conventional in the art, for rendering the resultant main magnetic field linear and uniform through the image region after correction for distortion caused by surrounding imager hardware. A gradient field means B selectively creates gradient magnetic fields across the main magnetic field in the examination region 10 . More specific to the preferred embodiment, linear gradients along an x, y, and z axis are selectively applied. The gradient field means includes a gradient field control means 20 which selectively provides current pulses to gradient coils 22 . The gradient field control means creates current pulses to cause a resultant gradient magnetic field pulse 24 with a preselected profile. However, due to eddy current magnetic fields, the profile of the current pulse and the gradient magnetic field pulse normally do not match. Rather, eddy currents tend to retard or round off a leading corner 26 of a gradient field pulse and cause a complementary lingering or residual magnetic field tail 28 after the end of the gradient field pulse. A pre-emphasis circuit 30 adds a current boost to the beginning portion of the current pulse to compensate for eddy currents. Typically, the current pre-emphasis means includes an adjustable frequency selecting means and an adjustable amplitude selecting means. Commonly, these are embodied in a series connected frequency filter or resonant circuit and amplifier. A plurality of these series connected filter and amplifier combinations are connected in parallel. The frequency selecting means is adjusted such that the pre-emphasis circuit adds a current component with the same frequency as one of the eddy current components. The amplitude adjusting means is adjusted to match the added current component amplitude to the eddy current component amplitude. A resonance excitation means C includes a radio frequency transmitter 40 for generating magnetic resonance excitation pulses and other magnetic resonance manipulation pulses. A radio frequency antenna 42 is disposed adjacent the examination region for transmitting radio frequency pulses into a calibration sample 44 disposed at the geometric center of the examination region. A receiver coil or antenna 46 receives radio frequency resonance signals, particularly free induction decay signals emanating from the sample. A radio frequency receiver 48 demodulates the resonance signals. In the preferred embodiment, the calibration sample is a water and copper sulfate solution at 50 grams per liter. This sample composition is selected for its very short relaxation times which allows the FID to be measured at very short intervals or repeat times. A timing and control means 50 controls the gradient field means 20 , the radio frequency transmitter 40 , and other circuit components to implement selected gradient and radio frequency pulse sequences, such as the sequence of FIGURE 2. That is, the gradient field control means 20 and the pre-emphasis circuit 24 generate current pulses to create the magnetic field gradient 24 , e.g. a z-gradient, across the examination region. Immediately following the end of the gradient 24 , the radio frequency transmitter 40 and radio frequency antenna 42 transmit a series of radio frequency magnetic resonance excitation pulses 52 in rapid succession. Typically, the radio frequency pulses are at 10 milliseconds intervals or longer and the series extends over 3.5 seconds or longer. However, greater time resolution is achieved by repetition of the gradient - [FID . . .] sequence with modified timing parameters, i.e. different delay times from the gradient ramp. Immediately following each radio frequency pulse, a resultant free induction decay signal 54 emanates from the sample, is received by the receiving means 48 , and digitized by an analog-to-digital converter 56 . Optionally, as shown in FIGURE 3, the intervals between radio frequency pulses 52 can be increased toward the end of the eddy current gradient field tail in accordance with the exponential decay characteristics. With reference to FIGURE 4, a plotting means 60 plots signal phase of the FID vs. time with phase unwrapping. That is, each time the phase crosses 2π, 2π is added to the phase such that the phase continues to grow ( 62 in the resultant plot) rather than cycling between zero and 2π ( 64 ). The slope of this plot, i.e. d/dt represents the resonant frequency. More precisely, this slope represents the median resonant frequency during the FID sampling interval. The resonant frequency fr is proportional to the magnetic field H, f r = 2πγH where γ is the gyromagnetic ratio of the sample. A slope or frequency determining means 66 measures the slope d/dt, hence, the resonance frequency. A magnetic field strength determining means 68 converts the frequency to an indication of the magnetic field strength using this relationship. A magnetic field strength vs. time memory 70 stores each determined magnetic field strength along with the time since the end of the gradient field pulse. A display means 72 , such as a CRT or a video monitor, displays the magnetic field strength vs. time curve. The eddy magnetic field 28 is roughly a decaying exponential curve. The resultant measured response is the sum of the exponential fields of each of the plurality of eddy currents. If the gradient is on for a time substantially longer than the longest time constant, the decaying exponential responses 28 to the pulse falling edge is substantially the inverse of the leading edge response 26 at the beginning of the gradient field pulse. To square the leading edge, one adds current components which compensate for the displayed eddy current magnetic field. This can be done by computing the current components that duplicate the inverse of the normalized eddy current field as it decays from zero to one. Preferably, an operator manually views the decay curve and determines a frequency or time constant and amplitude which offsets or cancels the trailing end of the curve. This is subtracted from the curve leaving a shorter exponentially decaying curve. The operator determines the time constant or frequency and amplitude for the next current component to remove the trailing edge of the remaining curve. This process is repeated until the curve is substantially removed or removed to within acceptable tolerances. The pre-emphasis circuit is adjusted in accordance with the selected frequencies and amplitudes. Alternately, because the data stored in the magnetic field strength vs. time memory 70 has already been digitized, appropriate computer implemented analysis and generation of the appropriate eddy current time constant and amplitudes may be undertaken. An appropriate computer based analysis means may be provided to perform this analysis and automatically adjust the pre-emphasis current means accordingly. The computer may use an iterative approach analogous to the above described manual approach or may calculate a multi-exponential fit to the data and determine all pre-emphasis time constant and amplitude settings in one operation. With reference to FIGURE 5, the present invention may also be utilized to monitor the spatial linearity, strength, and other characteristics of each generated gradient magnetic field. To this end, the sample is placed off center along one of the axes, e.g. at position 80 along the z axis. A gradient magnetic field 82 is applied and one or more radio frequency pulses 84 are applied to excite magnetic resonance at preselected times into the gradient magnetic field. The time into the gradient magnetic field is selected to be sufficiently long that eddy current effects have dissipated. The free-induction decay signal 86 following each radio frequency pulse is received and the resonant frequency determined. For a given physical location, the field should be constant during the pulse. Accordingly, the FID signals may be averaged. The sample 44 is then shifted to another point along the z axis, such as point 88 which is symmetric relative to the geometric center of the system to point 80 . The sequence is repeated, the resultant free induction decay measured, and resultant gradient field strength determined. If the gradient magnetic field is linear, then the magnetic field strength at these two points describes the gradient magnetic field along the z axis. Optionally, additional points may be taken along the z axis, as may be appropriate, to measure any non-linearities in the magnetic field strength. Analogous measurements may be made along the x and y or other axes. With reference to FIGURE 6, the same technique may also be utilized to measure temporal stability of the gradient field 82 . To this end, the series of radio frequency pulses 84 is applied after eddy currents die out and the FID signals 86 measured. All FIDs should be the same until the gradient current is altered. Rather than averaging, the frequency of each is recorded and analyzed for variation from other FID signals acquired under the same conditions. With reference to FIGURE 7, another use is to measure the linearity of the gradient response to an input signal. For a linear response, gradient strength, in the absence of eddy currents, should be directly proportional to the input signal to the gradient amplifiers. The response is monitored by repeating the sequence of Figure 5 at the same location but with different input levels or ramp heights 82, 82′.

Claims: A method of measuring a magnetic field caused by eddy currents comprising disposing a sample (44) in a main magnetic field in an examination region (10); causing a preselected gradient (24) across the magnetic field in the examination region (10); removing the gradient (24); immediately following removal of the gradient (24), inducing a series of free induction decay signals (54) from the sample (44); monitoring each free induction decay signal (54) and determining its frequency which frequency is proportional to magnetic field strength (H), whereby the magnetic field strength (H) at a sequence of times after removal of the gradient (24) is determined and wherein the monitoring step includes measuring signal phase () of the free induction decay signal (54) vs. time (t) and wherein the frequency determining step includes determining a slope (d/dt) of the measured phase () vs. time (t). A method according to Claim 1, wherein the step of inducing free induction decay signals (54) includes rapidly applying a series of radio frequency magnetic resonance excitation pulses (52), each of which induces magnetic resonance in the sample (44), and is followed by one of the free induction decay signals (54). A method according to Claim 2, wherein the radio frequency magnetic resonance excitation pulses (52) are applied at non-constant intervals. A method according to any one of the preceding claims further including again causing and removing the preselected gradient (24) and following removal of the gradient (24) inducing and monitoring a second series of free induction decay signals (54) with different delay times with respect to the removal of the gradient (24). A method according to any one of the preceding claims, wherein the sample (44) includes a substance that has a short relaxation time to facilitate inducing the free induction decay signals (54) in rapid succession. A method according to any one of the preceding claims, wherein the sample (44) includes a copper sulfate and water solution.

Inventors: HINKS RICHARD SCOTT, HINKS, RICHARD SCOTT
Applicants: PICKER INT INC, PICKER INTERNATIONAL, INC.
Application year: 1990
Main classification: G01R33
IPC codes: G01R33, A61B5