NUCLEAR MAGNETIC RESONANCE IMAGING METHODS

Abstract: To derive NMR imaging information an object is subjected to a static magnetic field on which is superimposed a sinusoidally time varying magnetic gradient field. A sequence of narrow 90 DEG rf pulses are applied around the instants when the gradient field has a zero value. The rf pulses have relative phase quadrature such that excited nuclei precess with accumulated phase in each sequence.

Description: NUCLEAR MAGNETIC RESONANCE IMAGING METHODS This invention relates to NMR imaging methods and has particular application in methods of obtaining NMR imaging information from solid objects. The NMR imaging of solid and other materials which have short spin-spin relaxation times has not been developed to the same level of sophistication as the NMR imaging of liquid and quasi-liquid materials. This is because of the Inherently greater Hnewidth in solids compared to liquids which forces the use of excessively large magnetic field gradients In prohibitively short times if traditional techniques are to be adopted. One kind of method that has been proposed to overcome the difficulties associated with imaging in solids uses multlpulse sequences to provide line narrowing. However, the biggest drawback with ultipulse sequencing concerns the large radio-frequency pulse power which is required. In part at least this is needed to overcome the large spectral width 1n an imaging experiment resulting from the presence of magnetic field gradients during the time of application of rf pulses. It is an object of the invention to provide a method in which the requirements of rf power are reduced. According to the invention a method of obtaining NMR imaging information from a solid object comprising subjecting the object to a static magnetic field, applying a gradient magnetic field which varies sinusoidally in amplitude about a zero value, applying a repetitive sequence of radio frequency pulses at instants of zero value of the said gradient field each of such frequency, magnitude and duration and successive pulses of the sequence being of such relative rf phase, that selected nuclei in the object precess with accumulated phase in each sequence, and measuring the resulting NMR signal from the object. Since the magnetic field gradients are at or close to zero during the time of application of rf pulses, considerable advantage is gained. Furthermore while rapid gradient switching is difficult to achieve on the tlmescale of multipulse experiments, it can nevertheless be achieved more easily when the coil windings generating the magnetic gradient fields are incorporated in a tuned circuit driven resonantly at a frequency with cycle time equal to twice the rf Interpulse period τ, thus generating a sinusoidally time varying magnetic gradient field. In carrying out the Invention the rf pulse sequence should preferably be such that the pulses of the sequence are separated by constant interpulse window lengths, equal to half the magnetic gradient field cycle time 2τ. In this arrangement the sign of the magnetic gradient field reverses in successive windows and the sequence therefore Is such as to make the spins accumulate phase in the rotating frame from an oscillating magnetic field offset as opposed to a static field offset. The spins specifically do not accumulate phase from a static resonance offset. This prevents magnet inhomogeneities and mistuning introducing errors into the phase accumulated by the spins due to the field gradients. As an added bonus It also removes chemical shift broadening. In some embodiments of the invention a sinusoidal offset magnetic field of spatially uniform amplitude 1s superimposed on the sinusoidally time-varying gradient magnetic field to shift spikes 1n the observed signal after Fourier transformation caused by rf phase errors away from areas of Interest. A sequence of 90° pulses and windows of length τ which fulfills the above requirements 1s 90_y-τ-(90_ χ -τ-90__y-τ-90_ χ -τ-90_ χ -τ-90_ y -τ-90_ χ -τ-) n (1 ) where the suffixes "-y" and "-x" refer to the rf phase of different pulses of a sequence relative to each other and n 1s an Integer. The product of the six pulses of the repeated cycle 1s unity as can be easily verified using simple matrices. For the case of an oscillating magnetic field resonance offset of mean magnitude B g the mean Ham1lton1an acting 1n odd numbered windows of the cycle 1s H n = -γf. B fl Σ l + Σ A (j .jj -3I 1 l (2) 9 j 1<j 1j z z and in even numbered windows 1 s Hr = +γf. B α Σ A + ∑ A (j . T.3 -3I 1 l (3) y j 1 <j 1 j z z In equations ( 2) and (3) the first term on the right-hand side is the Zeeman Interaction due to the offset field and the second term 1s the van-Vleck truncated dipolar Ham1lton1an. The symbols have their usual meaning. The procedure for calculation of the average Hamlltonlan acting over the complete six pulse cycle and hence the observed transverse magnetisation is well established. The Hamlltonlan "states" 1n the toggling reference frame to which the Hamiltonian in each window is transformed can be indicated as (Y,X,Z,X,7,Z). Thus to first order the average Hami-ltonian acting over the complete cycle is - j R = y 3 γ B g Σ I z (4) j It is seen that the dipolar Hamiltonian has been removed so achieving the aim of line narrowing, indeed it 1s removed over each consecutive half cycle. If the signal is observed 1n window six of each cycle quadrature magnetisations M χ . = M 0 cos (γB g t/3) (5) and M y ' = M 0 sin (γB g t/3) (6) are measured. M 0 is the equilibrium magnetisation. The main objective has been achieved. The spins precess 1n the rotating frame at a frequency proportional to Y- of the mean magnitude of the oscillating offset they experience. In the usual way the offset field 1s made spatially dependent by the Introduction of a gradient field of mean magnitude g. Including time dependence it is, for the x direction B α = ιrgx/2.sin(ιrt/τ) (7) Including the spin density p(x) the observed signal (in complex notation) is M(t) -- M χ .(t) + 1M y .(t) = M 0 p(x)exp(iγgxt/3)dx (8 ) x=-= This can be Fourier transformed with respect to time to recover the spin density profile p(x) = l/τrM 0 .γg/3.Re M(t) exp(-iγgχt/3)dt (9) t=o In the case that the field offset does not oscillate but is of constant magnitude, B c , then H 0 = He = - γt. B c Σ + Σ A (I 1 l3 -3I 1 l j z 1<j ij z z (10) The Hamiltonian states in the toggling frame are (Y,X,Z,X,7,Z), so that R = 0 (11) again to first order. This is an Important result. It means that the sequence 1s insensitVve to DC magnetic field offsets. If the mean Hamiltonian were not equal to zero- for constant fields, spins experiencing a combination of toggling and static fields would precess about an intermediate axis making recovery of spatial information difficult. This result also implies that chemical shift Information 1s lost and therefore that chemical shift broadening 1s removed. Susceptibility broadening is also removed. In order that the Invention may be more fully understand reference will now be made to the accompanying drawings 1n which: Figure 1 Illustrates diagramrnatlcally a pulse sequence embodying the Invention, Figure 2 Illustrates an image of a solid object obtained with the sequence of Figure 1, and. Figure 3 is a perspective view of the object Imaged in Figure 2. Referring now to Figure 1 there is shown therein a sinusoidal curve which represents an oscillating magnetic gradient field B g of cycle time period 2τ which 1s applied to an object In addition to a static magnetic field. This field is applied by standard gradient colls tuned by a series capacitor. In addition there are further windings in the form of a Hel holtz pair to provide a spatially uniform oscillating field. The object 1s also subjected to a radio-frequency pulse sequence which 1s shown in Figure 1 on the same tl escale. The sequence comprises 90° pulses applied at Instants when the gradient field B g has zero value. The pulses of the sequence are applied with equal interpulse spacing and the pulses of the sequence have a time period τ between successive pulses. The sequence comprises firstly a preparation pulse 90_ y followed after time interval τ by a repeated cycle of six pulses of cycle time 6τ. The relative rf phase of each pulse, is as indicated in Figure 1 and corresponds to equation (1) above. There is a time window marked S after the preparation pulse and at the end of each pulse cycle in which sampling of the observed signal takes place. The samples can be averaged and after Fourier transformation represent the spin density distribution of the excited nuclei, usually protons, in the gradient direction of the magnetic gradient field. Imaging information relating to the object can be built up by repeating the experiment with different orientations of gradient direction of the gradient field relative to the object. Figure 2 shows a two dimensional proton density cross sectional image across a rubber bung 7 mm In diameter and 8 mm in length with a 1.5 mm slot cut across a diameter along its length. The actual bung is shown 1n perspective in Figure 3. The image was obtained by back projection of 36 profiles each of 71 data points. 256 cycles of the sequence were used for each profile with a pulse gap of 20 μs, the oscillating frequency of the sinusoidal gradient field B g being 25 kHz. The sample was rotated by 5° between projection measurements. Alternatively in a system with multiple orthogonal gradient coils rotation of the gradient direction can be effected electrically. A suitable arrangement is a pair of magnetic coils positioned 1n relative quadrature to each other. The direction of the gradient magnetic field 1s determined by the ratio of the electric currents In the two coils. By changing this ratio the gradient direction can be changed. In addition a sinusoidal offset field of spatially uniform amplitude was superimposed on the gradient field. This is because radio frequency pulse phase transients and other effects give rise to spurious precession of the magnetisation of order a few hundred Hertz in multi-pulse work and consequently ,to a low frequency spike in the Fourier transformed spectrum. To overcome the problem spectra are often recorded a little off resonance so that the low frequency artefacts in the signal do not corrupt the results. The equivalent procedure 1n this case 1s to include a sinusoidal offset field so that the oscillating field (including the gradient) is non-zero everywhere over the sample. The low frequency («0HZ) spike then occurs at the side of the profile where it can be suppressed without loss of Image information. As an alternative to the six pulse sequence defined by equation (1) a 12 pulse sequence can be used by taking the six pulse sequence and following 1t by its Her itlan conjugate to give the following sequence: 90_y-τ-(90_ χ -τ-90_y-τ-90_ χ -τ-90_ χ -τ-90_ y -τ-90_ χ -τ-90 χ -τ-90y-τ-90 χ -τ-90 χ -τ-90y-τ-90 χ -τ-) n (12) In the sequence defined by equation (12) sampling takes place in the sixth and twelfth window of the pulse cycle. An advantage of the 12 pulse sequence is that it 1s less sensitive to pulse flip angle defects than the six pulse sequence. In either the six pulse or the 12 pulse sequence it is possible to sample 1n every third window in order to Increase the sampling rate. If this rate of sampling is used then 1t is necessary to interchange the quadrature magnetisations of each alternate sampling.

Claims: CLAIMS 1. A method of obtaining NMR imaging information from a solid object comprising subjecting the object to a static magnetic field, applying a gradient magnetic field which varies sinusoidally 1n amplitude about a zero value, applying a repetitive sequence of radio frequency pulses at Instants around the zero value of the said gradient field each of such frequency, magnitude and duration and successive pulses of the sequence being of such relative rf phase, that selected nuclei in the object precess with accumulated phase 1n each sequence, and measuring the resulting NMR signal from the object. 2. The method as claimed in Claim 1 in which the sequence of radio frequency pulses is such that the pulses of the sequence are separated by constant Interpulse window lengths equal to half the magnetic gradient field cycle time. 3. The method as claimed in either one of the preceding claims 1n which a sinusoidal offset magnetic field of spatially uniform amplitude is superimposed on the sinusoidally time-varying gradient magnetic field to shift spikes in the measured signal after Fourier transformation caused by rf phase errors away from areas of interest. 4. The method as claimed in any one of the preceding claims in which an initial radio frequency preparation pulse is applied before the repetitive sequence of radio frequency pulses. 5. The method as claimed in Claim 4 in which the preparation pulse and the sequence of pulses each comprise 90° rf pulses, the rf phases of which are 1n relative phase quadrature with each other and according to the following pattern 90_y-τ-(90_ χ -τ-90_y-τ-90_ χ -τ-90_ χ -τ-90_y-τ-90_ χ -τ-) n where -y and -x indicate rf relative phase quadrature, τ Is a predetermined time interval and n 1s an Integer equal to the number of repetitions of the sequence. 6. The method as claimed In any one of the preceding claims In which the measuring of the resultant nmr signal from the object takes place 1n the time Interval between the preparation pulse and the sequence of pulses and between the last two pulses of a repeated sequence. 7. The method as claimed 1n any one of the preceding claims which is repeated with different directions of gradient magnetic field relative to the sample. 8. The method as claimed in Claim 7 1n which the direction of a gradient magnetic field is determined by applying two sinusoidally varying electric currents of predetermined ratio to two magnetic colls positioned 1n relative quadrature to each other.

Inventors: MCDONALD PETER JOHN, TOKARCZUK PAWEL FRANCISZEK, MCDONALD, PETER, JOHN, TOKARCZUK, PAWEL, FRANCISZEK
Applicants: NAT RES DEV, NATIONAL RESEARCH DEVELOPMENT CORPORATION
Application year: 1990
Main classification: G01R33
IPC codes: A61B5, G01N24, G01R33