MULTIPLEXED MAGNETIC RESONANCE IMAGING OF VOLUMETRIC REGIONS

Abstract: In a magnetic resonance imaging system, a multiplexing technique is used to image volumetric regions by performing a number of scans within the recovery time of the spin systems. The multiple scans acquire signal information from sub-regional volumes by phase-encoding in the direction of the excitation gradient.

Description: This invention relates to magnetic resonance imaging (MRI) method and, in particular, to techniques for simultaneously gathering nuclear magnetic resonance (NMR) signal information for a plurality of regions of a volume being imaged. In magnetic resonance imaging it is possible to create several different types of images of a subject. One common image type is the two-dimensional (2-D) image of a planar region or slice through the subject. In planar imaging a thin slice of the subject is selectively excited to the exclusion of the surrounding volume. The two major dimensions of the selected slice are then spatially encoded by the application of orthogonal magnetic field gradients. In the Fourier technique, one is a phase encoding gradient and the other a frequency encoding gradient. The gathered NMR signal information is processed by two-dimensional Fourier transformation to reconstruct an image of the desired plane. In two-dimensional imaging the time required to gather the necessary signal information from a live subject, the scan time, is generally relatively short. In this short imaging time the signal information will be relatively uncontaminated by artifacts due to breathing, movement of the subject, the other motion-related factors. The resultant image, with its low content of motion artifacts, will be a relatively sharp and clearly defined picture of the internal structure or anatomy of the subject. A second image type is the three-dimensional image. The three-dimensional image advantageously presents structural and anatomical detail in all three dimensions. Such an image may be formed by selectively exciting a predetermined volume of the subject. The volume is once again spatially encoded by magnetic field gradients and is scanned through to gather NMR signal information from points throughout the volume. The gathered NMR signal information is processed to reconstruct an image of the selected volume. However, the time required to scan through the entire volume is considerably longer than the time needed to scan a single plane. Hence, the three-dimensional image is intrinsically more sensitive to motion of the subject. The high degree of motion artifact contamination thus usually results in a decision to perform multiple 2-D scans rather than the three-dimensional scan. By forming images of a plurality of parallel slices of the volume, the same region can be viewed by sequentially stepping through the sequence of 2-D images instead of viewing the three-dimensional volume image. However, the taking of separate 2-D images is also a lengthy process. The scan time can be reduced, though, by acquiring signal information in rapid sequence from small groups of the full number of slices. After one slice has been interrogated, there is normally a waiting period, or recovery time, before the spin system of the slice has relaxed and the slice can be interrogated again. During this recovery time others of the desired slices can be interrogated. The acquisition of NMR information from a number of slices during the waiting period is commonly known as multiple slice NMR imaging. However, it is generally not possible to interrogate all slices, but only a limited number of slices, during a single waiting period. Moreover, the multiple slice technique suffers from a poorer signal to noise ratio than the three-dimensional imaging technique under the conditions that the desired waiting period before the sequence can be repeated is too short to obtain all of the NMR information from all of the slices. In medical diagnostic MRI the investigator tries to optimize the choices of the recovery time, termed T R , and the time interval between the excitation of a slice and the subsequent reception of spin echo signals from the slice, termed T E , or echo time. The choices are made in order to distinguish the different properties of tissue as revealed by their NMR characteristics. In general, a long T E interval will yield better T₂ contrast between different tissue, and a short T R period will yield better T₁ contrast. But a long T E and a short T R will both result in reducing the number of slices that can be interrogated during a single T R period. The optimum choices of radio frequency (RF) pulse timings for good tissue contrast are thus in opposition to the choices that optimize the speed of multiple slice data acquisition. Accordingly, it would be desirable to provide greater flexibility to the investigator who is trading off these choices. Furthermore, it would be desirable to perform magnetic resonance imaging in a manner which preserves all the information of the three-dimensional image but with shorter scan times that are less sensitive to motion. Also, the desired technique should provide a signal to noise characteristic at least equal to that of the multiple slice technique. EP-A-0 091 008 (General Electric Company) discloses a method for imaging a slab in a subject using a magnetic resonance imaging system, the method comprising: (i) selectively exciting the slab with a radio-frequency (RF) pulse; (ii) applying to the slab a field gradient in a first direction, in time coincidence with said RF pulse; (iii) applying to the slab a frequency-encoding gradient in a second direction, the frequency-encoding gradient having a direction of rotation in a first plane orthogonal to said first direction and a direction of rotation in a second place which includes the first direction; (iv) acquiring NMR signal information; (v) repeating steps (a) to (d) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient in one of the first and second planes; (vi) repeating steps (a) to (e) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient in the other of the first and second planes; and (vii) reconstructing a image using the acquired NMR signal information. In one mode of operating the method disclosed in EP-A-0 091 008, steps (i) to (iv) are repeated rapidly on the same slab so that the magnetization of the spins does not have time to return to its equilibrium between sequences. This enables information concerning the spatial distribution of T₁ values (the spin lattice relaxation time) in the slab to be collected. According to the present invention, there is provided a method for imaging a region in a subject, the region being conceptually divided into a plurality of slabs and each slab being conceptually divided into a plurality of slices using a magnetic resonance imaging system, the method comprising: (a) exciting a first one of said slabs or a first set of said slices, the set comprising one slice from each slab, with a radio-frequency (RF) pulse having a first frequency characteristic; (b) applying to the region a slab or slice selective field gradient in a first direction, in time coincidence with said RF pulse; (c) applying to the region a frequency-encoding gradient in a second direction orthogonal to said first direction, the frequency-encoding gradient having a direction of rotation in a first plane orthogonal to the first direction and a direction of rotation in a second plane which includes the first direction; (d) acquiring NMR signal information; (e) repeating steps (a) to (d) a plurality of times, each time using an RF excitation pulse having a different frequency characteristic; (f) repeating steps (a) to (e) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient in one of the first and second planes; (g) repeating steps (a) to (f) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient in the other one of the first and second planes; and (h) reconstructing an image of the region using the acquired NMR signal information, wherein steps (a) to (e) are carried out within a predetermined recovery time for the spin system being imaged. Preferred features of the method of the present invention are recited in the dependent claims. In one embodiment, the present invention provides a magnetic imaging technique in which spatially encoded image information is gathered in a multiplexing manner. During a single predetermined recovery time interval of the spin system, a plurality of volumetric region or "slabs" are sequentially selectively excited. As each region is excited, its volume is phase encoded by magnetic field gradients in two directions, one of which is in the direction of the selective excitation gradient. After the last region has been scanned, the sequence repeats, but with a change in one of the phase encoding gradients. Periodically, the other phase encoding gradient is also varied. Thus, during the recovery time of a single region, a plurality of regions are scanned, with multiplexed variation of the phase encoding gradients over a number of sequences. In this way, volumetric information from a plurality of regions is gathered quickly and with good sensitivity. The principles of the present invention are applicable to Fourier transform imaging and other imaging geometries such as 3-D back projection as well. In the drawings: FIGURE 1 illustrates a block of material with a volumetric region which is to be imaged; FIGURE 2 illustrates multiple slice imaging of the volumetric region of FIGURE 1; FIGURE 3 illustrates selective excitation of the multiple slices of FIGURE 2; FIGURE 4 illustrates the division of the volumetric region of FIGURE 1 into slabs for multiplexed imaging in accordance with the present invention; FIGURES 5a-5d illustrate multiplexed imaging of the slabs of FIGURE 4 ; FIGURE 6 illustrates selective excitation of the slabs of FIGURE 4; FIGURE 7 illustrates a second division of the volumetric region of FIGURE 1 into slabs for multiplexed imaging in accordance with the present invention; FIGURE 8a-8d illustrate multiplexed imaging of the slabs of FIGURE 7; FIGURE 9 illustrates selective excitation of the slabs of FIGURE 7; FIGURE 10 illustrates signal waveforms for the acquisition of NMR information in accordance with the principles of the present invention; and FIGURE 11 illustrates a technique of ordering contiguous slices of a multiple number of slabs. Referring to FIGURE 1, a rectangular block 100 of a material which is to be imaged is shown. As an example, suppose that volumetric region 102 of the block is to be imaged. FIGURE 1 is shown with reference to the axes drawn adjacent the FIGURE. That is, the longest dimension of the block 100 is aligned with the z axis, and the rectangular face of the block is in the x-y plane Suppose then that a simple selective three dimensional image is formed of the region 102 using the conventional spin echo signal technique. This technique is described in the article "Operating Principles of NMR Imaging", by B. C. Hill. et al ., published in Computed Tomography of the Whole Body , Vol. 2, ed. by J. R. Haagg et al. (1983), pp. 1038 et seq . at page 1051. If the block is a volume of a living subject with a typical T₁ for animal anatomy, an NMR spin echo signal might typically be acquired 30 milliseconds after excitation of a given region by a radio frequency (RF) excitation pulse, ie. T E =30 ms. However, the spin systems must be allowed to recover their orientation with the static magnetic field to a certain degree before the spin systems can be reexcited. Typically, this recovery time T R may be 200 milliseconds. Two further assumptions are necessary to complete the three dimensional imaging example, which are the Fourier encoding of the y and z dimensions of the region. It will be assumed that the y dimension comprises 256 data lines and that the z dimension comprises sixteen phase-encoded gradient levels for the Fourier transform processing of image information. With these assumptions, it is seen that a single NMR signal can be acquired during every successive 200 millisecond interval. To acquire the necessary signals for the y-direction dimension, a sequence of 256 such acquisitions or scans would be performed, each with a different y-direction gradient value. Since the volume comprises 16 phase-encoded levels in the z-direction, sixteen such sequences of 256 scans are performed, each of the sixteen with a different z-direction gradient value. To reconstruct an image of the volume, three Fourier transforms are performed on the scan data, a so-called 3D Fourier transform. The time required to complete the scanning is thus (200 ms) (256 y-gradients) (16 z-gradients) = 13.6 minutes. In addition, reconstruction using a 3D Fourier transform of this amount of data requires significant processing time. However, the three dimensional image will exhibit a good signal to noise characteristic, since the entire volume is excited during each scan. For comparison purposes, a signal-to-noise factor of 100 will be designated for the three dimensional image. The same volume may be imaged by acquiring the NMR signal information using a two-dimensional multiple slice technique. In this technique, image information from several planar slices is acquired sequentially during a single recovery interval T R . Referring to FIGURE 2, for multiple slice imaging the volume 102 is conceptually divided into sixteen parallel planar slices which are to be imaged. The slices are not imaged in their numerical sequence, but in a staggered order. This is to prevent the selection of immediately adjacent slices simultaneously, or in immediately consecutive order. Since the RF excitation pulses normally will not sharply define the planar boundaries of a slice, there is in practice a small region of overlap at the boundary region when adjacent slices are excited simultaneously or in immediate succession. Signals from the overlap area will be saturated due to multiple excitations, resulting in the loss of proper image contrast in these areas. To prevent this undesirable condition, the slices are selected as represented by FIGURE 3. Referring to FIGURE 3, a multiple slice selection sequence is illustrated. The block 100 is shown in relation to the variation in frequency ω z which is provided by the slice selecting z gradient. The frequency characteristics of the RF slice selecting pulses are illustrated below the block 100. The first slice is selected as shown by the S1 characteristic, which selects slice 1 with an ω₁ frequency characteristic. Using the same timing relationship as in the first example, an NMR spin echo signal may be acquired 30 milliseconds after the slice-selecting RF pulse is applied. Approximately 20 milliseconds are allowed for the completion of NMR signal acquisition and for the effects of gradient switching to subside, whereupon slice 5 is selected by an RF pulse exhibiting an ω₅ characteristic as shown by the S2 characteristic. Spin echo signals from slice 5 are acquired 30 milliseconds later and 20 milliseconds thereafter slices 9 and 13 are excited in a similar manner, as indicated by the S3 and S4 characteristics. Thus, a different slice is selected and NMR signal information acquired every 50 milliseconds. Within the 200 millisecond recovery time T R , four slices are interrogated. The four slice sequence is then repeated for slices 2, 6, 10 and 14, followed by two additional four slice sequences for the remaining slices. The pulse and gradient timings for the multiple slice sequence are shown in EP-A-219206 Using the same criteria as in the three dimensional example, the total scanning time is calculated as follows. Since four slices are interrogated every 200 milliseconds, it takes 800 milliseconds to interrogate all sixteen slices. To acquire the necessary information to discriminate in the y-direction, the sixteen slices must be interrogated 256 times, each with a different y-gradient level. The minimum total scanning time is thus 4 scans times (200 msec) (256 y-gradients) = 3.4 minutes. This is an appreciable time saving as compared to the three dimensional technique, since each of the four scans takes less than a minute to complete. The signal to noise characteristic is based upon the number of independent measurements that are combined in each image, which is one-sixteenth that of the 3D example. The square root of sixteen is four, indicating a factor of four reduction in signal to noise. Using the same basis as the factor of 100 for the three dimensional technique, the multiple slice technique will have a signal to noise factor of 25. In accordance with the present invention NMR signal information is acquired from the volumetric region 102 of the block 100 by conceptually dividing region 102 into a number of volumetric slabs which are spatially encoded by phase encoding in the direction of the slice or slab-selecting gradient. Referring to FIGURE 4, the region 102 is divided into four slabs, A, B, C and D. The slabs are respectively phase-encoded in the slab-selecting z-gradient direction as indicated by phases ₁, ₂, ₃, and ₄. The shadings of the slabs correspond to this phase-encoding. To provide the same resolution as in the multiple slice technique, individual slices of the slabs are selected as illustrated in FIGURES 5a-5d. In a first selective excitation, slices A1, B1, C1 and D1 are selected as indicated in FIGURE 5a. Since these four slices are excited simultaneously, NMR signal information will be produced by the four slices at the same time. The signals from the four slices are discriminated by the phase encoding in the z direction. The signals from the respective slices will exhibit respectively different phases. After the first four slices are selected and their signal information gathered, slices A2, B2 C2 and D2 are similarly interrogated as indicated in FIGURE 5b. Likewise, the remaining eight slices are interrogated as indicated in FIGURES 5c and 5d. The charateristics of the RF excitation pulses used to select the slices in FIGURES 5a-5d are shown in FIGURE 6. The first RF excitation pulse is represented by the ω₁ frequency characteristic, and contains frequency bands ω A1 ω B1 , ω C1 , and ω D1 to select the first four slices. The frequency characteristics of the second, third and fourth RF pulses are shown as characteristics ω₂, ω₃ and ω₄. Returning to FIGURES 5a-5d, the sequence there represented is repeated 256 times for the 256 levels of the phase-encoding gradient in the y-direction. Then, 256 additional sequences are performed with z-gradient variation, this time with different phases being assigned to the four slabs according to their z-position. In practice, a new level is used for the phase-encoding gradient pulse in the z-direction to accomplish this. Two more sequences of 256 y-gradient lines are then performed, each with a different z-gradient level, for a total of four iterations. This corresponds to the selection of four slices by each RF pulse and provides information needed to discriminate the four slices of each excitation during Fourier transform reconstruction. Using the same criteria as in the previous examples, the scanning time is calculated as follows. The interrogation of each group of four slices in FIGURES 5a-5d is performed in 30 milliseconds. Allowing another 20 milliseconds between each scan for gradient effects to dissipate, it is seen that information from each group of four slices may be gathered every 50 milliseconds. Thus, all sixteen slices are interrogated every 200 milliseconds. Multiplying this period by the number of y-gradient levels gives (200 msec) (256 y-gradients) = 51.2 seconds. Repeating this scanning period by the four iterations necessary for z-gradient variation gives (51.2 sec) (4 z-gradients) = 3.4 minutes, the same time saving as in the multiple slice example. The comparable signal to noise factor for this embodiment of the present invention is 50, since each separate slice location is selected only one-fourth as often as the 3D example, but four times more than in the 2D multiple slice technique. A major difference as compared with the multiple slice exmaple is that phase-encoding in the z-gradient direction is used to discriminate the signals received simultaneously from four slices, a condition not present in the multiple slice technique. In accordance with a further aspect of the present invention, another imaging technique is provided which retains the shorter scanning times of the 2D multiple slice technique but with an improved signal-to-noise figure. In this second technique, the region 102 is again conceptually divided into four slabs A, B, C and D as shown in FIGURE 7. NMR signal information is then sequentially acquired from each slab in its entirety. Within each slab, z-direction phase encoded slices are designated as shown in FIGURE 7. Phase-encoded signal information is thus gathered from every slice in the slab simultaneously and is discriminated within each slab by Fourier transformation of the phase-encoded information. The selective excitation of each slab in sequence is performed by application of a z directed gradient together with RF pulses of frequency content illustrated in FIGURE 9. As in FIGURES 3 and 6, line ω z indicates the direction of variation of the z gradient. Slab A is selected by an RF pulse having a frequency content ω A as illustrated by characteristic ω₁; slab B by characteristic ω₂; slab C by characteristic ω₃; and slab D by characteristic ω₄. The slab-selecting RF pulses are easier to tailor in terms of their frequency content than the RF pulses used to select multiple thin slices in the 2D multiple slice technique discussed above, and generally can be applied in less time than the RF pulses needed to select thin slices. The sequence of signal acquisition for the second embodiment is illustrated in FIGURES 8a-8d. In FIGURE 8a, the first slab A is interrogated, with phase encoding in the z-direction of the slab as indicated by ₁ - ₄. Thereafter, slab B is similarly interrogated as shown by FIGURE 8b, followed by slabs C and D as shown by FIGURES 8c and 8d. In each case, the same phase-encoding in the z-direction of the slab may be employed. A pulse sequence useful for performing the multislab techniques of the present invention is shown in FIGURE 10. A 90° RF pulse 20 is applied to the block 100. This RF pulse exhibits a frequency content ω₁, corresponding to the spectra of ω₁ in either FIGURE 6 or FIGURE 9. The RF pulse 20 is coincident in time with a slice or slab-selecting z-directed gradient G zs , shown as pulse 21 in FIGURE 10. Following slice or slab selection three additional gradients are applied. The G x gradient 60 preconditions the spin system for the frequency encoding gradient in the x-direction. One level of a level-variable G y gradient 64 is applied. The slab is phase encoded in the z direction by applying one level of a z-directed phase encoding G zp gradient shown at 66. In the present examples of FIGURES 4 and 7, with four phases in each excitation, the G zp gradient would exhibit one of four possible levels. After the three spatially encoding gradients have been applied a 180° RF pulse 22 is applied to the slab or slices in the presence of a G zs gradient 23. This insures that only the spin systems of the selected slab or slices are inverted for the generation of a spin echo signal, leaving the spin systems of other slabs or slices unaffected while they recover to their orientation with the static magnetic field. Thereafter spin echo signals 24 develop in the selected slab or slices, and are sampled during a sampling period 26 in the presence of a G x gradient 62. When the effects of the switched gradients have subsided a scan is performed on slab B or the second group of slices. A 90° RF pulse 30 is applied to the block 100 in the presence of slice or slab selecting gradient G z , shown at 31. The frequency content ω₂ of the RF pulse 30 is that of the ω₂ waveforms in FIGURE 6 or FIGURE 9. Spatially encoding gradients 70, 74 and 76 are applied as before, and the spin systems are inverted by a frequency selective 180° RF pulse 32 in the presence of G zs gradient 33. Spin echo signals 34 develop from slab B or slices A2, B2, C2 and D2, which are sampled during an interval 36 in the presence of a G x gradient 72. The breaks in the lines of FIGURE 10 indicate that numerous scans of various slabs or slice groups may be performed within the selected recovery time T R of the spin systems. The last scan of the sequence is shown to the right of the breaks, where a 90° RF pulse 50 selects the last slab or slice group n, which in the example of four slabs or slice groups is slab D or slices A4, B4, C4 and D4. The RF pulses 50 and 52, the slice selecting G zs gradients 51 and 53, the spatially encoding gradients 90-96, and the sampling interval 56 and spin echo signals 54 are similar to those of the preceding scans, but are specific to the last slab or slice group n, which are slab D and slices A4, B4, C4 and D4 in the examples. After the last spin echo signals 54 have been acquired and gradient effects have subsided, the full sequence of FIGURE 10 is repeated, this time using the same value as before for one of the phase-encoding gradients and a different value for the other phase-encoding gradient. In a constructed embodiment, the G zp gradient 66 is held constant while the G y gradient 64 is changed to the next line value. The term "line" as used here does not refer to a single physical line through a slab, since phase-encoded information is gathered from the entire volumetric slab or several slice volumes during each scan. Rather, the term line means a particular value of the G y gradient which is used in Fourier transformation of the acquired data. The new value for the G y gradient line is used for gradients 64, 74 and 94 during the second sequence. The sequence is thereafter repeated, each time with a new G y line value, until the G y gradient has been sequenced through all its line values, which is 256 in the present examples. Thereafter, the sequence of 256 y-line variable iterations is repeated once more for each value of the G zp gradient. For the four phase variations of each slab or slices shown in FIGURES 4 and 7, the scanning is completed when the G zp gradient has been sequenced through all four of its possible levels. In the present examples one slab or slice group is scanned once in about 50 milliseconds and all four slabs or groups are therefore scanned in 200 milliseconds. The complete scanning time for the block 100 thus takes (200 milliseconds) (256 y-gradients) ( 4 z-gradients) = 3.4 minutes. The signal to noise factor for the technique illustrated by FIGURES 7 and 8a-8d is the same as that of the technique illustrated by FIGURES 4 and 5a-5d. When a subject is imaged using the multislab technique one of the slab volumes is generally, although not necessarily, centered with respect to the center, or null point of the slab selection gradient field. As shown in FIGURE 11, slab F is centered with respect to the center Z o of the G zp gradient field. On one side of the Z o point the phase rotation due to the G zp gradient is relatively positive, and on the other side of the Z o point the G zp gradient phase rotation is relatively negative. Since the Z o point is located at the interface of phase encoded slices 2 and 3 of slab F, the interface will have a known phase with respect to the reference 300 used during phase sensitive detection of the NMR signals. Accordingly, the other slice boundaries of slab F will also have known phases. In the simplified example of FIGURE 11, the Z o reference phase is shown as ₁, and the sequence of phase values of slab F is ₃, ₄, ₁, ₂, and ₃, which phase sequence is represented by the slice shadings. The phase sequence properly identifies the sequence of the slices in the slab during image reconstruction. While slab F of FIGURE 9 is centered with respect to the gradient null point Z o , the same is not true of the other slabs E and G. These slabs will be in regions of the G zp gradient field where the gradient phase rotation is either entirely positive or entirely negative with respect to the field at Z o . While a uniform bias gradient could be used to shift the effective G zp gradient null point to each slab center during slab selection and spatial encoding as described in EP-A-0 134 701 it may be desirable to use the same G zp gradient for each slab. In such case, those slabs remote from the center slab will be encoded with phase values dependent upon their spatial location, causing a reordering of the phases at the respective slice boundaries. For instance, the interface between slices 2 and 3 of slab E in FIGURE 11 is at a point Z E in the G zp gradient field. This causes the boundary of slices 2 and 3 of slab E to be at a phase of ₂ with respect to the phase reference 300, instead of ₁ as was the case of slab F. Similarly, the center of slab G is at a phase of ₄ with respect to the phase reference 300. This reordering of the sequence of phases within each slab would cause the slices to be organized in an incorrect sequence in the different slabs unless this resequencing is taken into consideration during image reconstruction. To reorder the slices of each slab to their proper physical sequence, a correction is made which is a function of the G zp gradient field and the spatial locations of the respective slabs. The G zp gradient is known to have a characteristic variation, expressed as a change in Gauss (10 -4 T) per centimeter. For slab E, the physical distance between point Z E and the null point Z o is also a known quantity. A correction factor is thus determined by multiplying G zp (Gauss/cm) by slab distance (Z E - Z O ) and multiplying this intermediate term by the gyromagnetic ratio γ for the magnetic field. This correction factor is then used as a product term applied to the slices of slabs remote from the center slab. The correction reorders the phases of slices of the remote slabs to the same sequence as the slices of the center slab for presentation of all slabs and slices in their correct physical sequence. The present inventive techniques are also applicable to other reconstruction techniques, such as three dimensional back projection reconstruction. In 3D back projection, the excitations and signal readout will be multiplexed as described above, with variation of the direction of rotation of the frequency-encoding gradient in a plane orthogonal to that of the slab selecting gradient. This multiplexing is then repeated with variation in a plane which includes that of the slab selecting gradient. Image reconstruction is then performed on the acquired information using 3D back projection. Moreover, the back projection and Fourier techniques can be combined, as by performing rotation perpendicular to the slab-selecting gradient and phase-encoding in parallel with the slab-selecting gradient.

Claims: A method for imaging a region (102) in a subject (100), the region (102) being conceptually divided into a plurality of slabs (A, B, C, D) and each slab (A, B, C, D) being conceptually divided into a plurality of slices (A1-4, B1-4, C1-4, D1-4), using a magnetic resonance imaging system, the method comprising: (a) exciting a first one (A) of said slabs or a first set (A1, B1, C1, D1) of said slices, the set comprising one slice from each slab, with a radio-frequency (RF) pulse (20) having a first frequency characteristic (ω₁); (b) applying to the region (102) a slab or slice selective field gradient (G zs ) in a first direction, in time coincidence with said RF pulse (20); (c) applying to the region (102) a frequency-encoding gradient (G x ) in a second direction orthogonal to said first direction, the frequency-encoding gradient (G x ) having a direction of rotation in a first plane orthogonal to the first direction and a direction of rotation in a second plane which includes the first direction; (d) acquiring NMR signal information (54); (e) repeating steps (a) to (d) a plurality of times, each time using an RF excitation pulse having a different frequency characteristic (ω₂, ω₃, ω₄); (f) repeating steps (a) to (e) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient (G x ) in one of the first and second planes; (g) repeating steps (a) to (f) a plurality of times, each time using a different direction of rotation of the frequency-encoding gradient (G x ) in the other one of the first and second planes; and (h) reconstructing an image of the region using the acquired signal information, wherein steps (a) to (e) are carried out within a predetermined recovery time for the spin system being imaged. The method of claim 1, wherein the directions of rotation of the frequency-encoding gradient (G x ) in the first and second planes are determined by applying to the region (102): (a) a level-variable, phase-encoding gradient (G y ) in a third direction orthogonal to said first and second directions, at a first level; and (b) a level-variable, phase-encoding gradient (G zp ) in said first direction, at a first level, and wherein the different directions of rotation in steps (f) and (g) are determined by changes in the levels of the level-variable gradients (G y , G zp ). The method of claim 1 or claim 2, wherein the RF pulse (20) excites a slab (A, B, C, or D). The method of claim 1 or claim 2, wherein the RF pulse (20) excites a set of slices (A1, B1, C1, D1). The method of claim 4, wherein the slices in each set are contiguous. The method of claim 4, wherein the slices in each set are not contiguous. The method of any one of claims 4 to 6, wherein the slices in each set are similarly phase-encoded. The method of any one of claims 4 to 6, wherein each slice in each set are differently phase-encoded. The method of any one of claims 1 to 8, wherein in step (h) the image is reconstructed by three-dimensional Fourier transformation of the acquired NMR-signal information (54). The method of any one of claims 1 to 8, wherein in step (h) the image is reconstructed by 3-D back projection.

Inventors: COMPTON, RUSSELL ALLEN, KRAMER, DAVID MOSES, YEUNG, HONG-NING
Applicants: SIEMENS AKTIENGESELLSCHAFT
Application year: 1986
Main classification: G01R33
IPC codes: A61B5, A61B10, G01R33