(Radiology. 2000;216:891-899.)
© RSNA, 2000
Optimized Single-Slab Three-dimensional Spin-Echo MR Imaging of the Brain1
John P. Mugler, III, PhD,
Sumi Bao, PhD,
Robert V. Mulkern, PhD,
Charles R. G. Guttmann, MD,
Richard L. Robertson, MD,
Ferenc A. Jolesz, MD and
James R. Brookeman, PhD
1 From the Dept of Radiology, Box 800170, University of Virginia School of Medicine, Charlottesville, VA 22908 (J.P.M., J.R.B.); Dept of Nuclear Engineering, Massachusetts Institute of Technology, Cambridge (S.B.); and Depts of Radiology, Childrens Hospital (R.V.M., R.L.R.) and Brigham and Womens Hospital (R.V.M., C.R.G.G., F.A.J.), Harvard Medical School, Boston, Mass. Received Jan 4, 1999; revision requested Mar 17; final revision received Nov 26; accepted Dec 17. J.P.M., S.B., R.V.M., C.R.G.G., F.A.J., and J.R.B. supported in part by grant NS-35142 from the National Institute of Neurological Disorders and Stroke. Address correspondence to J.P.M. (e-mail: jpm7r@virginia.edu).
 |
ABSTRACT
|
|---|
The development and optimization of spin-echobased, single-slab, three-dimensional techniques for magnetic resonance imaging of the whole brain are described. T1-weighted and T2-weighted image sets with a volume resolution of 1 mm3 and fluid-attenuated inversion-recovery image sets with a volume resolution of 3 mm3 were obtained in acquisition times of less than 10 minutes per image set.
Index terms: Brain, abnormalities, 10.121411, 10.363, 10.3637, 10.871 Brain, MR, 10.121411, 10.121413, 10.121416 Magnetic resonance (MR), pulse sequences, 10.121411, 10.121413, 10.121416 Magnetic resonance (MR), rapid imaging, 10.121411, 10.121413, 10.121416 Magnetic resonance (MR), three-dimensional, 10.121411, 10.121413, 10.121416
 |
INTRODUCTION
|
|---|
Magnetic resonance (MR) imaging is well established as the modality of choice for the diagnostic exploration of many brain diseases. The majority of MR techniques used for brain imaging are based on spin-echo (SE) acquisitions because they provide a wide range of useful image contrast properties that highlight pathologic changes and are resistant to image artifacts from a variety of sources such as radio-frequency or static-field inhomogeneities. However, despite numerous advances in MR technology, a major benefit of clinical MR brain imaging envisioned more than a decade ago has not been realizednamely, to combine the desirable properties of SE-based imaging with the intrinsic capability of MR imaging for true volume acquisitions and thereby routinely obtain high-spatial-resolution images of arbitrary orientations (1).
In practice, the limitations of three-dimensional (3D) techniques have hindered their routine use for imaging the brain. The contrast of 3D gradient-recalled-echo techniques is, in many instances, unsatisfactory when compared with that of SE imaging. The 3D gradient-recalled-echo techniques are more susceptible to image artifacts from off-resonance effects, and true T2 weighting is difficult to attain. As a result, only short-repetition-time, T1-weighted (eg, fast low-angle shot) or short-repetition-time, steady-state (eg, gradient-recalled acquisition in the steady state) sequences have gained some clinical acceptance. T2-weighted imaging still depends largely on two-dimensional (2D), multisection conventional SE or, more prevalently, fast SE methods (25).
Multislab 3D fast SE techniques (6,7) offer the potential for high-spatial-resolution imaging of the whole brain but nonetheless possess several important limitations. Radio-frequency (RF) excitation along the section-select direction is nonuniform, and, as a result, 3D reformations generally are contaminated by artifacts at the junctions between slabs. In addition, because of slab profile effects, some of the outer sections in each slab typically are discarded, thus decreasing the efficiency. Power deposition is relatively high and may compromise the coverage attained per unit time, particularly at a field strength of 1.5 T and greater. Each slab undergoes unwanted magnetization-transfer effects from the large number of high-flip-angle off-resonance RF pulses applied to the other slabs during the acquisition (8). The use of additional contrast preparations, such as magnetization-transfer contrast or diffusion, requires that the repetition time be increased or the anatomic coverage compromised.
Recent advances in whole-body MR imager technology, especially higher performance gradient and receiver systems, now present the opportunity for developing single-slab 3D techniques that capitalize on SE acquisitions. In this article, we report our initial efforts to develop and optimize single-slab 3D brain imaging based on fast SE techniques (25). The major technical issues that had to be addressed are presented, followed by preliminary imaging results for 3D T1-weighted, T2-weighted, and fluid-attenuated inversion-recovery, or FLAIR, imaging.
 |
Materials and Methods
|
|---|
Theoretic Simulations
A computer-based theoretic model of MR imaging, previously described (9), was used in conjunction with experimental results to understand the relationship between various types of nonideal imager behavior (eg, gradient-waveform errors, RF-pulse phase instabilities) and potential image artifacts. Simulated images of uniform objects were generated for both ideal conditions and with specific perturbations such as an error in the area of the preread gradient waveform. Computer programs were written in FORTRAN and implemented on a workstation (Ultra-60; Sun Microsystems, Mountain View, Calif).
Pulse Sequence Design
The major technical issues that arose during the development of the single-slab 3D pulse sequences included receiver-system dynamic range limitations due to the large imaging volume and image ghost artifacts related to the fidelity of the gradient waveforms. Other important sequence design considerations included the configuration of the spoiler (crusher) gradient waveforms and the implementation of a half-Fourier acquisition mode.
Dynamic range limitations.It is well known that the signal dynamic range required for Fourier-based encoding of large volumes, such as the whole brain, can exceed that provided by the 16-bit analog-to-digital convertors typically used in the receiver systems of MR imagers (10,11,12). With insufficient dynamic range, quantization (digitization) noise can dominate random thermal noise and consequently degrade image quality, primarily by decreasing the image signal-to-noise ratio (SNR). Compared with a thin-section 2D acquisition, a thick 2D section or 3D slab requires that the gain setting of the receiver is lowered to accommodate the larger signal amplitude. As a consequence, the amplitude of the thermal noise from the coil and subject is decreased relative to the peak signal amplitude, and at sufficiently low receiver-gain settings, the amplitude of the noise will be less than the resolution of the digitized signal. For these conditions, the noise intensity is determined by the resolution of the analog-to-digital convertor, not by the inherent thermal noise, and is effectively increased relative to the peak signal amplitude, which consequently reduces the SNR.
Ultimately, whether the thermal noise is digitized adequately depends not only on the resolution of the analog-to-digital convertor but also on the sampling frequency (ie, on the degree of oversampling), the RF coil characteristics, the receiver-gain setting, and the noise performance of the other components in the receiver system. Nonetheless, experimental characterization of the noise as a function of receiver gain can be used to establish whether quantization noise is an issue. To perform such a characterization of the noise for our specific application, the relationship between receiver gain and image noise intensity was determined by acquiring image sets with the standard head RF coil over the relevant range of receiver-gain settings. The noise intensity was measured as the mean of the random noise in the image background within a region of interest that included a minimum of 5,000 pixels and excluded any visible artifacts from motion or data truncation.
Although several methods for overcoming dynamic range limitations have been described (1013), none has been adopted for widespread use, since the majority of clinical MR imaging applications currently depend on thin-section 2D acquisitions. Of the previously described methods, only a variable gain approach (14,15), wherein the receiver gain is varied during the acquisition to match the signal level, requires no modification of the pulse sequence. To address the dynamic range limitation associated with imaging the whole brain by using a single thick slab, an adaptation of the gain-changing approach was developed as illustrated in Figure 1. This approach takes advantage of the multiple receiver channels available for phased-array RF coils to simultaneously record the signal from the standard head coil at two receiver-gain settings. The low-gain data, which are contaminated by quantization noise, are used for only the central, high-signal portion of spatial-frequency space. In practice, approximately 0.03% of the total number of spatial-frequency-space components in the 3D data set are from the low receiver-gain setting. The reconstruction software is programmed to automatically scale and combine the data from the two receiver channels. This two-receiver-channel method has been integrated into one of the two imagers used in this study (Vision; Siemens Medical Systems, Iselin, NJ) and is currently being developed for the other imager (Signa; GE Medical Systems, Milwaukee, Wis).

View larger version (20K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 1. Schematic shows the acquisition methodology used to overcome receiver dynamic range limitations. Data from the head coil are acquired through the standard receiver channel at high gain and also through a second (phased-array) channel at low gain. The low-gain central spatial-frequency space data and high-gain outer spatial-frequency space data are combined automatically during image reconstruction. ADC = analog-to-digital convertor, CH = channel, PREAMP = preamplifier.
|
|
In the context of dynamic range limitations, another important consideration is the numeric representation used for the image reconstruction calculations. The sampled data from the analog-to-digital convertors, typically 16-bit integers as discussed earlier, can be converted into any arbitrary numeric representation, for example, 32-bit integer or floating-point numbers, for mathematic calculations. This is particularly relevant from a historic perspective because many early (1980s) MR systems used 16-bit integer arithmetic for image calculations. The 216 numeric values available by using 16-bit integers generally are adequate for 2D image calculations, but they can introduce additional noise from mathematic round-off errors for 3D data sets. A 32-bit numeric representation, which reduces the effect of round-off errors to a negligible level, was used for the image reconstruction calculations in all 3D protocols in this study. Thus, in contrast to the limitations of the sampling process discussed earlier, the image SNR was not (further) reduced owing to an insufficient dynamic range during the calculation process.
Gradient-related ghost artifacts.The high sensitivity of Carr-Purcell-MeiboomGill SE trains, as used in fast SE sequences, to errors in the phase of the transverse magnetization is well documented (16,17). Common sources of these errors include eddy currents and RF-pulse phase instabilities. However, any differences between 2D and 3D sequences in their sensitivity to phase errors has not, to our knowledge, been recognized. Our investigation was prompted by the finding that preliminary versions of the 3D pulse sequences demonstrated substantial ghost artifacts (peak artifact intensity of about 10% of the mean signal intensity within the object of interest), whereas comparable 2D sequences showed negligible ghost artifacts.
For a given sequence configuration, the specific source and form of gradient-induced phase errors is dependent on the details of the gradient-system and eddy-current behaviors associated with a particular imager. For an example of a potential source for gradient-induced phase errors, we chose to study an error in the area of the preread gradient waveform (ie, the gradient waveform applied along the signal readout direction between the excitation and first refocusing RF pulses). Experimental and simulated images from 2D and 3D sequences with identical timing were assessed for the presence and severity of ghost artifacts as a function of an error in the area of the preread gradient waveform that ranged between 0% and 0.2% of the area corresponding to the readout gradient waveform. Simulated images also were used to predict the fraction of image intensity present as ghost artifacts as a function of gradient-area error and refocusing RF-pulse flip angle.
Spoiling configuration.Spoiler or crusher gradients are used in fast SE pulse sequences between the refocusing RF pulse and data sampling window of each echo to dephase the free induction decay generated by the refocusing RF pulse. These gradient waveforms are often placed along the section-select direction (7). However, in 3D sequences, this configuration may substantially decrease the efficiency of the pulse sequence because data sampling does not begin until both the spoiling and 3D phase-encoding gradients are applied. A more time-efficient approach, illustrated in Figure 2, is simply to extend the duration of the readout gradient waveform a sufficient amount to displace signals originating at the refocusing RF pulses out of the data sampling window. Thus, encoding in the third dimension and spoiling are performed simultaneously, which permits a decreased echo spacing compared with that achieved when spoiling is performed along the section-select direction.

View larger version (22K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 2. Diagram shows the time-efficient gradient configuration for spoiling the free induction decay from the refocusing RF pulses. The diagram illustrates the portion of the pulse sequence between two consecutive refocusing RF pulses. The duration of the readout gradient waveform is extended beyond the data acquisition period to displace signals originating at the refocusing RF pulses out of the data sampling window. DAQ = data acquisition period; GPE = phase-encoding gradient; GRO = readout gradient; GSS = section-select gradient; RF = applied RF pulses.
|
|
Half-Fourier.To decrease imaging times, half-Fourier acquisition mode was implemented as an option in the sequences. The half-Fourier image reconstructions were obtained by using a straightforward extension of a previously described algorithm (18) to three dimensions.
The implementation of the half-Fourier acquisition differs in an important way from that used in another recently developed 3D SE-based pulse sequence (19). In our case, a full-Fourier acquisition was obtained along the phase-encoding direction in spatial-frequency space that is modulated by signal decay during the echo train (the in-plane phase-encoding direction), and a half-Fourier acquisition was obtained along the direction of spatial-frequency space that has constant signal-strength weighting (the 3D phase-encoding direction).
This is in contrast to the method used for 3D inversion-recovery half-Fourier single-shot rapid acquisition with relaxation enhancement (19), which places the half-Fourier acquisition along the phase-encoding direction of spatial-frequency space that is modulated by signal decay during the echo train. Thus, with half-Fourier single-shot rapid acquisition with relaxation enhancement, the weighting along the half-Fourier direction of spatial-frequency space has the form of a low-pass filter, which potentially results in blurring for tissues with relatively short T2 relaxation times.
MR Imaging
Imaging was performed by using two 1.5-T commercial whole-body imagers (Vision and Signa). Both MR systems were equipped with echo-planar-capable gradient systems having maximum gradient strengths and minimum rise times of 25 mT/m and 300 µsec (Vision) and 23 mT/m and 300 µsec (Signa). Initial testing and optimization of the sequences were performed by using phantoms that were supplied with the imagers by the manufacturers and that contained water doped with nickel sulfate or nickel chloride. Preliminary images obtained in the heads of healthy volunteers and patients were then acquired by using the optimized sequence configurations after obtaining informed consent. Institutional review board approval was obtained for all human imaging protocols. The pulse sequence parameters used for 3D T1-weighted, T2-weighted, and fluid-attenuated inversion-recovery imaging are summarized in the Table.
The volunteers (14 men, four women; age range, 2558 years) included subjects who were in good health and were not seeing a physician for any neurologic disorder. A small number of patients (four female patients; age range, 1245 years) were included irrespective of specific diseases to demonstrate the image contrast following administration of gadopentetate dimeglumine (Magnevist [0.1 mmol per kilogram of body weight]; Berlex Laboratories, Wayne, NJ). The clinical diagnoses in these patients were optic glioma (n = 1), medulloblastoma (n = 1), and multiple sclerosis (n = 2).
A consensus evaluation of the images from volunteers and patients was performed by four of the authors (C.R.G.G., J.P.M., R.L.R., R.V.M.). The gross contrast behaviors were assessed in comparison with those established for 2D conventional SE and fast SE imaging, and the images were evaluated for the presence and severity of physiologic (eg, motion) and nonphysiologic artifacts. In this preliminary study, no formal grading of the image characteristics was performed.
 |
Results
|
|---|
Dynamic Range Limitations
The noise intensity as a function of receiver gain for both imagers is shown in Figure 3. The image noise intensity was approximately constant for receiver gains ranging from the highest value to an intermediate value roughly midway between the gain limits. As the receiver gain was decreased below this intermediate value, the noise intensity increased, at first slowly and then more rapidly, as quantization noise secondary to insufficient analog-to-digital convertor dynamic range overcame thermal noise as the dominant noise source. The receiver-gain range over which the noise intensity was constant included the gain settings typically used for 2D clinical MR acquisitions. At the receiver-gain values appropriate for single-slab 3D imaging (<20 dB in Fig 3), these results indicated that dynamic range limitations may increase the noise intensity substantially above the lower limit determined by thermal noise, and consequently decrease the image SNR.
The images in Figure 4 from two 3D T2-weighted acquisitions of the whole brain demonstrate the improvement in SNR attained by using the two-receiver-channel method. For Figure 4a, the image was acquired by using a single receiver gain of 71 dB (see Fig 3; an absolute gain of 71 dB corresponds to a relative gain of 15 dB in the figure). This was the gain value required to accommodate the highest signal amplitude from the whole head for the selected repetition and echo times and hence also the gain value that would be selected by the standard prescan software. From Figure 3, a noise intensity of 20.4 corresponds to this gain setting (Vision MR system). The image in Figure 4b was acquired by using receiver gains of 56 and 82 dB. These absolute gains correspond to relative gains of 0 and 26 dB, respectively, in Figure 3. The low-gain value was chosen to accommodate the highest signal amplitude from the whole head for the range of the repetition and echo times used in all experiments; the high-gain value was chosen to be within the range in which the image noise intensity was constant.

View larger version (131K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4a. Transverse sections from 3D T2-weighted MR acquisitions of the whole brain acquired by using (a) a single absolute receiver gain of 71 dB and (b) the two-receiver-channel method with absolute gains of 56 and 82 dB. The SNR of b is more than double that of a. Pulse sequence parameters were as follows: 2,400/128 (repetition time msec/effective echo time msec); matrix, 256 (readout) x 140 (first phase-encoding direction) x 220 (second [3D] phase-encoding direction); field of view, 22 x 16 x 23 cm; echo spacing, 4.0 msec; echo train length, 70; half-Fourier; acquisition time, 9.4 minutes.
|
|

View larger version (126K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4b. Transverse sections from 3D T2-weighted MR acquisitions of the whole brain acquired by using (a) a single absolute receiver gain of 71 dB and (b) the two-receiver-channel method with absolute gains of 56 and 82 dB. The SNR of b is more than double that of a. Pulse sequence parameters were as follows: 2,400/128 (repetition time msec/effective echo time msec); matrix, 256 (readout) x 140 (first phase-encoding direction) x 220 (second [3D] phase-encoding direction); field of view, 22 x 16 x 23 cm; echo spacing, 4.0 msec; echo train length, 70; half-Fourier; acquisition time, 9.4 minutes.
|
|
The noise intensity values corresponding to these gain settings (Fig 3) are 91.5 and 9.0. Considering that more than 99.9% of the spatial-frequency space data for Figure 4b is from the high receiver-gain setting, the noise intensity should correspond to that for the high receiver-gain setting; thus, the noise for Figure 4b should decrease by a factor of 20.4 divided by 9.0, or 2.3, compared with that in Figure 4a. This prediction agrees well with the experimental resultthe noise for the image in Figure 4b is decreased by a factor of 2.4 compared with that for the image in Figure 4a, which resulted in an increase in the SNR by the same factor because the signal level is unaffected. The actual noise intensities in Figure 4a and 4b were scaled by using an arbitrary multiplicative factor, compared with the specific values stated earlier; however, this scaling is identical for the two images.
Gradient-related Ghost Artifacts
Images from the 3D pulse sequences demonstrated visible ghost artifacts for preread gradient-area errors of only a few hundredths of a percent, whereas images from 2D pulse sequences in which identical timing parameter values were used showed negligible ghost artifacts over the complete range of gradient-area errors. The only substantial artifact noted in the 2D images was an intensity shading along the readout direction for the larger values of the gradient-area error. Figure 5 shows a comparison of 2D and 3D images for a gradient-area error of 0.1% and a refocusing-RF-pulse flip angle of 180°. The 3D images (Fig 5b, 5d) demonstrate both diffuse and focal degradation compared with the 2D images (Fig 5a, 5c).

View larger version (115K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5a. Comparison of ghost artifacts for transverse T2-weighted (a, c) 2D and (b, d) whole-brain 3D MR acquisitions obtained by using identical pulse sequence timing and a preread gradient-area error of 0.1%. The number of signals acquired for the 2D sequence was adjusted to yield an SNR equivalent to that for the 3D sequence. Pulse sequence parameters for both sequences were as follows: 3,000/163; matrix, 256 (readout) x 140 (in-plane phase-encoding direction); field of view, 22 x 16.5 cm; echo spacing, 5.1 msec; echo train length, 70; refocusing RF pulse, 180°. For the 3D sequence, the matrix size in the second phase-encoding direction was 80, and the corresponding field of view was 20 cm. The section thicknesses for the 2D and 3D sequences were 5.0 and 2.5 mm, respectively. Images c and d show the same data as images a and b, respectively, with the display contrast adjusted to allow better depiction of the ghost artifacts that appear in d as discrete regions of high signal intensity outside of the head. These ghost artifacts represent signals from within the head that were encoded incorrectly owing to the preread gradient-area error and thus appear displaced along the in-plane phase-encoding direction (left-to-right in d).
|
|

View larger version (129K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5b. Comparison of ghost artifacts for transverse T2-weighted (a, c) 2D and (b, d) whole-brain 3D MR acquisitions obtained by using identical pulse sequence timing and a preread gradient-area error of 0.1%. The number of signals acquired for the 2D sequence was adjusted to yield an SNR equivalent to that for the 3D sequence. Pulse sequence parameters for both sequences were as follows: 3,000/163; matrix, 256 (readout) x 140 (in-plane phase-encoding direction); field of view, 22 x 16.5 cm; echo spacing, 5.1 msec; echo train length, 70; refocusing RF pulse, 180°. For the 3D sequence, the matrix size in the second phase-encoding direction was 80, and the corresponding field of view was 20 cm. The section thicknesses for the 2D and 3D sequences were 5.0 and 2.5 mm, respectively. Images c and d show the same data as images a and b, respectively, with the display contrast adjusted to allow better depiction of the ghost artifacts that appear in d as discrete regions of high signal intensity outside of the head. These ghost artifacts represent signals from within the head that were encoded incorrectly owing to the preread gradient-area error and thus appear displaced along the in-plane phase-encoding direction (left-to-right in d).
|
|

View larger version (81K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5c. Comparison of ghost artifacts for transverse T2-weighted (a, c) 2D and (b, d) whole-brain 3D MR acquisitions obtained by using identical pulse sequence timing and a preread gradient-area error of 0.1%. The number of signals acquired for the 2D sequence was adjusted to yield an SNR equivalent to that for the 3D sequence. Pulse sequence parameters for both sequences were as follows: 3,000/163; matrix, 256 (readout) x 140 (in-plane phase-encoding direction); field of view, 22 x 16.5 cm; echo spacing, 5.1 msec; echo train length, 70; refocusing RF pulse, 180°. For the 3D sequence, the matrix size in the second phase-encoding direction was 80, and the corresponding field of view was 20 cm. The section thicknesses for the 2D and 3D sequences were 5.0 and 2.5 mm, respectively. Images c and d show the same data as images a and b, respectively, with the display contrast adjusted to allow better depiction of the ghost artifacts that appear in d as discrete regions of high signal intensity outside of the head. These ghost artifacts represent signals from within the head that were encoded incorrectly owing to the preread gradient-area error and thus appear displaced along the in-plane phase-encoding direction (left-to-right in d).
|
|

View larger version (82K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5d. Comparison of ghost artifacts for transverse T2-weighted (a, c) 2D and (b, d) whole-brain 3D MR acquisitions obtained by using identical pulse sequence timing and a preread gradient-area error of 0.1%. The number of signals acquired for the 2D sequence was adjusted to yield an SNR equivalent to that for the 3D sequence. Pulse sequence parameters for both sequences were as follows: 3,000/163; matrix, 256 (readout) x 140 (in-plane phase-encoding direction); field of view, 22 x 16.5 cm; echo spacing, 5.1 msec; echo train length, 70; refocusing RF pulse, 180°. For the 3D sequence, the matrix size in the second phase-encoding direction was 80, and the corresponding field of view was 20 cm. The section thicknesses for the 2D and 3D sequences were 5.0 and 2.5 mm, respectively. Images c and d show the same data as images a and b, respectively, with the display contrast adjusted to allow better depiction of the ghost artifacts that appear in d as discrete regions of high signal intensity outside of the head. These ghost artifacts represent signals from within the head that were encoded incorrectly owing to the preread gradient-area error and thus appear displaced along the in-plane phase-encoding direction (left-to-right in d).
|
|
Figure 6 shows a theoretic calculation of the fraction of the total image intensity present as ghost artifacts as a function of both the preread gradient-area error and the flip angle of the refocusing RF pulses. These results indicate that the severity of ghost artifacts for the 3D sequence is a strong function of the flip angle of the refocusing RF pulses and is substantially worse than that for a comparable 2D sequence when relatively high refocusing RF-pulse flip angles are used. The predicted rapid decrease in the artifact fraction for the 3D sequence as the flip angle is decreased was corroborated by experimental measurements. Note that the artifact fraction for the 2D sequence is lower than that for the 3D sequence regardless of the refocusing RF-pulse flip angle for the 3D sequence.

View larger version (30K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 6. Graph shows the fraction of the total image intensity present as ghost artifacts as a function of the preread gradient-area error and the flip angle of the refocusing RF pulses. The severity of ghost artifacts for the 3D sequence is a strong function of the refocusing RF-pulse flip angle and decreases rapidly with decreasing flip angle. The predicted artifact fraction for the 2D sequence is lower than that for the 3D sequence, regardless of flip angle. Simulation parameters were as follows: 2,500/144; matrix, 256 (readout) x 140 (phase encoding); field of view, 22 x 22 cm; echo spacing, 4.0 msec; echo train length, 70; T2, 300 msec.
|
|
The results in Figure 6 suggest that an effective strategy to mitigate artifacts from gradient-area errors in the 3D sequences is to use low-flip-angle refocusing RF pulses. However, simply using constant, low-flip-angle refocusing RF pulses yields a substantial loss of signal. As an alternative, pseudosteady-state reduced flip angles, developed by Alsop (20) for application to 2D fast SE sequences, provide an improved transfer of magnetization to the transverse plane (ie, increased signal) and allow the desired contrast (ie, the contrast corresponding to a flip angle of 180°) to be maintained over a wide range of T1 and T2 values. The pseudosteady-state condition is a constant echo amplitude during the SE train for the idealized case of no relaxation (ie, infinite T1 and T2 values). Pseudosteady-state flip angles for our sequence configurations were derived for a terminal flip-angle value of 90° (20). This refocusing RF-pulse flip-angle series provided comparable artifact reduction but increased signal intensities compared with those associated with a constant flip angle of 90°.
Although 2D images acquired by using a 180° refocusing RF-pulse flip angle were predicted (Fig 6) to yield an artifact fraction only slightly lower than that for 3D images acquired by using a 90° flip angle, the severity of artifacts for the 2D case, on the basis of visual inspection of both simulated and experimental images, was substantially lower. This apparent discrepancy was explored by using simulated images to investigate the nature of the ghost artifacts as a function of the position within the 2D section profile. Evaluation of these images indicated that the low level of ghost artifacts in the 2D case arises for two reasons: (a) The section-selective refocusing RF pulse contains a range of flip angles between 0° and 180°; as the flip angle decreases with distance from the center of the section, so does the relative intensity of the ghost artifacts. (b) The location of the ghost artifacts is a strong function of the flip angle.
The relationship between the location of ghosts and the flip angle was recognized previously by Norris et al (16) in connection with fast 2D rapid acquisition with relaxation enhancement pulse sequences in which low-flip-angle refocusing RF pulses were used. Since the artifact seen in a 2D image is the result of integrating over the section, and hence over the corresponding range of flip angles, there is limited constructive interference among the artifacts arising from different portions of the section. Thus, in contrast to the 3D case, the artifacts are spread more evenly across the image and appear primarily as a low-intensity diffuse smearing that for most practical purposes is irrelevant.
Whole-Brain 3D Imaging Examples
Figures 7 and 8 present sagittal, coronal, and transverse T1- and T2-weighted images reconstructed at 1.0-mm thickness from 9.4-minute 3D acquisitions covering the whole head. For both image sets, the spatial resolution along each dimension was approximately 1 mm. These image sets demonstrate the typical image quality and contrast obtained from the T1- and T2-weighted single-slab 3D acquisitions, and they illustrate the capability to obtain high-quality, high-spatial-resolution image reformations in any plane. Figure 9 presents an additional example of sagittal, coronal, and transverse 3D T1-weighted images, reconstructed at 1.0-mm thickness, in a patient with an optic glioma. This gadolinium-enhanced image set demonstrates high contrast between the tumor and surrounding tissue, which is comparable to that for gadolinium-enhanced 2D T1-weighted conventional SE imaging. Again, with 1-mm spatial resolution along each dimension, high-spatial-resolution images can be obtained in any plane.

View larger version (167K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 7. T1-weighted sagittal (left column), coronal (middle column), and transverse (right column) MR images of the whole brain reconstructed at 1.0-mm thickness from a 9.4-minute 3D acquisition. These images demonstrate contrast comparable to that seen on 2D conventional SE images and illustrate the capability to obtain high-spatial-resolution T1-weighted images in any orientation from the single-slab 3D acquisition. Pulse sequence parameters were as follows: 480/28; matrix, 256 (readout) x 140 (first phase-encoding direction) x 220 (second [3D] phase-encoding direction); field of view, 22 x 16.5 x 22 cm; echo spacing, 4.0 msec; echo train length, 14; half-Fourier.
|
|

View larger version (151K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 8. T2-weighted sagittal (left column), coronal (middle column), and transverse (right column) MR images reconstructed at 1.0-mm thickness from a 9.4-minute 3D acquisition of the whole brain. These images demonstrate contrast comparable to that seen on 2D fast SE images and illustrate the capability to obtain high-spatial-resolution T2-weighted images in any orientation from the single-slab 3D acquisition. Pulse sequence parameters were as follows: 2,400/128; matrix, 256 (readout) x 140 (first phase-encoding direction) x 220 (second [3D] phase-encoding direction); field of view, 22 x 16.5 x 22 cm; echo spacing, 4.0 msec; echo train length, 70; half-Fourier.
|
|

View larger version (58K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 9. Gadolinium-enhanced T1-weighted sagittal (left), coronal (middle), and transverse (right) MR images with 1.0-mm thickness from a 3D acquisition of the whole brain in a patient with an enhancing lesion (straight arrows) identified as an optic glioma. Note the high contrast between the tumor and surrounding tissue, which is comparable to that seen on gadolinium-enhanced 2D T1-weighted conventional SE images. The isotropic 1-mm resolution permits high-spatial-resolution images to be obtained in any plane. The artifactual intensities (curved arrows) seen on the image obtained in the sagittal plane are an aliasing artifact secondary to a suboptimal choice of the primary imaging plane and slab thickness. Pulse sequence parameters were as follows: 600/12; matrix, 220 (readout) x 160 (first phase-encoding direction) x 140 (second [3D] phase-encoding direction); field of view, 22 x 16 x 14 cm; echo spacing, 6.5 msec; echo train length, 16; acquisition time, 14 minutes; full-Fourier.
|
|
Transverse fluid-attenuated inversion-recovery images are shown in Figure 10 from a 9.5-minute 3D acquisition covering the whole head. In addition to thin, contiguous sections, the single-slab 3D fluid-attenuated inversion-recovery technique provided the advantage over a 2D or multislab 3D (21) fluid-attenuated inversion-recovery technique of eliminating potential signal ambiguities arising from cerebrospinal fluid motion during the long inversion delay (22). Thus, as illustrated in Figure 10, the signal intensity from cerebrospinal fluid is effectively nulled throughout the brain, including regions of relatively fast cerebrospinal fluid motion such as the aqueduct and the area surrounding the brainstem.

View larger version (73K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 10. Transverse fluid-attenuated inversion-recovery MR images with 3.0-mm thickness from a 9.5-minute 3D acquisition. Thin, contiguous sections covering the whole brain were obtained with a single 3D acquisition. The single-slab nature of the pulse sequence results in uniform suppression of the signal intensity from cerebrospinal fluid throughout the brain, including the aqueduct and the area surrounding the brainstem. Pulse sequence parameters were as follows: 6,200/128; inversion delay, 2,000 msec; matrix, 256 (readout) x 140 (first phase-encoding direction) x 76 (second [3D] phase-encoding direction); field of view, 22 x 16.5 x 23 cm; echo spacing, 4.0 msec; echo train length, 70; half-Fourier.
|
|
The gross contrast behaviors of the T1-weighted, T2-weighted, and fluid-attenuated inversion-recovery image sets, as illustrated in Figures 710, were judged to be comparable to their 2D counterparts; however, further evaluation will be required to determine if there are any subtleties to the contrast properties beyond those already documented for 2D fast SE techniques (5). Early versions of the pulse sequences used for human imaging demonstrated ghost artifacts related to gradient-system performance, as described in the previous section. These artifacts were not visible in the optimized pulse sequences. Ghost artifacts from pulsatile blood flow were not observed; the absence of these commonly seen artifacts is attributed to the combination of the nonsection-selective RF pulses and a spatial-presaturation slab inferior to the brain.
The only other artifact sometimes noted was aliasing (wraparound) of signal from outside the imaging volume into the images because of a suboptimal choice of parameter values. An example of this is shown in the sagittal image in Figure 9. The combination of a transverse primary orientation and an insufficient slab thickness resulted in some aliasing artifact (eg, in the corpus callosum), even though a spatial-presaturation slab was placed inferior to the imaging volume.
 |
Discussion
|
|---|
The two-receiver-channel scheme used for overcoming the dynamic range limit imposed by a 16-bit analog-to-digital convertor provided a substantial increase in the image SNR compared with using the standard single receiver-gain setting. When pushing the limits of acquisition speed and spatial resolution, such gains in SNR are obviously crucial to obtaining the highest quality results. Similar experience has been reported recently for 3D gradient-recalled-echo techniques (15). We hope that the development of new 3D imaging methods, such as those demonstrated in this article, that depend on signal dynamic ranges greater than 16 bits will motivate the MR imager manufacturers to implement an increased dynamic range as a standard feature of their products.
In developing the 3D pulse sequences, we discovered an important difference between 2D and 3D SE-trainbased techniques related to gradient-system performance. For the specific type of gradient-waveform error studied, the 2D sequences did not demonstrate substantial shading artifacts until the gradient error was on the order of 0.1%, whereas the 3D sequences demonstrated visible ghost artifacts for a gradient error that was an order of magnitude smaller. Thus, successful implementation of the 3D SE-train sequences requires substantially higher gradient-waveform fidelity than that needed for 2D SE-train sequences. In practice, we found that a combination of precautions was effective in maintaining gradient-related ghost artifacts at a negligible level (ie, peak artifact intensity less than approximately 1% of the mean signal intensity within the object of interest) in the 3D sequences. These included (a) ensuring that eddy-current compensation was within manufacturer specifications, (b) using consistent waveform slew rates and spacing when gradient waveforms were placed immediately adjacent to short-duration RF pulses, and (c) using low-flip-angle refocusing RF pulses for the shortest echo spacings.
The 3D T1-weighted images yielded contrast comparable to that of conventional 2D SE images, but the contrast between brain white and gray matter appeared to be less than that typically achieved by using 3D gradient-recalled- echo techniques such as fast low-angle shot or magnetization-prepared rapid gradient echo. However, for contrast agentenhanced 3D imaging, the SE-based technique may be preferred, considering the reported limitations of gradient-recalled-echo techniques for depicting contrast-enhanced lesions (23,24). Compared with the current standard for contrast-enhanced imaging, conventional 2D SE, the 3D SE-based technique also provides the advantages of thin contiguous sections and fewer flow-related artifacts, which suggests the potential for the 3D technique to replace conventional 2D gadolinium-enhanced SE imaging.
Despite the promising results presented for single-slab 3D SE-based imaging in this article, the issue of which type of 3D SE-based imagingsingle-slab or multislabwill prove to be superior for imaging the whole brain in arbitrary orientations with isotropic resolution is still open to debate. With multislab 3D imaging, increased echo spacing is possible while still providing sufficient coverage to image the whole brain with an acquisition time comparable to that used with single-slab 3D techniques. Thus, the multislab approach potentially can achieve a higher sampling efficiency (the fraction of the echo spacing for which the analog-to-digital convertor is on) and therefore an increased SNR per unit time.
However, multislab implementations typically discard a number of the outer sections from each slab to decrease slab boundary artifacts, which in turn decreases sampling efficiency. In addition, the use of echo trains that combine gradient echoes and SEs with 3D SE-based imaging (25,26) decreases the theoretic sampling efficiency advantage of multislab techniques. Although substantial pulse sequence improvements have been made with regard to reducing slab-boundary artifacts (27,28), these artifacts are yet to be completely eliminated. Recently, a potential alternative, the shifted interleaved multivolume acquisition, has been introduced (29). This method eliminates slab-boundary artifacts in exchange for ghosting in an orthogonal direction. Overall, considering the results presented to date for the single-slab and multislab approaches, further development and direct comparison of both techniques seems warranted.
Although our studies primarily concentrated on obtaining image sets with approximately isotropic 1-mm spatial resolution, different configurations of the single-slab 3D techniques may prove useful for other applications, particularly when high-spatial-resolution image reformations in arbitrary planes are not needed. One example is the 3D fluid-attenuated inversion-recovery images in Figure 10 in which 1.0-mm sections were traded for an extended repetition time. With comparable or increased SNR, the current T1- and T2-weighted sequences could acquire 2.0-mm-thick sections covering the whole head in less than 5 minutes by using 1-mm in-plane resolution or acquire 3.5-mm-thick sections in less than 10 minutes by using 0.5-mm in-plane resolution.
In conclusion, MR pulse sequences for single-slab 3D brain imaging, based on fast SE techniques, have been developed. Optimization of the pulse sequence structure to overcome SNR losses induced by quantization noise, to suppress ghost artifacts secondary to gradient-system performance, and to maximize sampling efficiency permitted acquisition of 3D image sets of the whole brain in less than 10 minutes by using T1- or T2-weighted contrast and a volume resolution of 1 mm3 or fluid-attenuated inversion-recovery contrast and a volume resolution of 3 mm3. The development of these techniques revealed important pulse sequence design issues unique to the single-slab 3D approach. In our preliminary experience, these pulse sequences routinely provided high-quality image sets, free of any substantial image artifacts. These techniques have the potential to improve the state of the art both for routine MR imaging and for specialized applications such as surgical planning or lesion quantification for therapy monitoring. The ability to collect such high-spatial-resolution 3D data sets in a clinically reasonable time will affect the diagnosis, treatment, and research for disseminated diseases of the brain, such as multiple sclerosis, Alzheimer disease, and metastatic neoplasms.
 |
FOOTNOTES
|
|---|
J.P.M. and J.R.B. receive research support from Siemens Medical Systems, and R.V.M. and F.A.J. receive research support from GE Medical Systems.
Author contributions: Guarantors of integrity of entire study, J.P.M., F.A.J.; study concepts, J.P.M., R.V.M., C.R.G.G., F.A.J., J.R.B.; study design, J.P.M.; definition of intellectual content, J.P.M., R.V.M., C.R.G.G., J.R.B.; literature research, J.P.M., R.V.M.; clinical studies, R.L.R., R.V.M., S.B., C.R.G.G.; experimental studies, J.P.M., S.B., R.V.M., C.R.G.G.; data acquisition, J.P.M., S.B., R.V.M., C.R.G.G.; data analysis, J.P.M., S.B.; manuscript preparation, J.P.M.; manuscript editing, J.P.M., R.V.M., C.R.G.G., J.R.B.; manuscript review, R.L.R., F.A.J.
 |
REFERENCES
|
|---|
-
Pykett IL, Buonanno FS, Brady TJ, Kistler JP. True three-dimensional nuclear magnetic resonance neuro-imaging in ischemic stroke: correlation of NMR, X-ray CT and pathology. Stroke 1983; 14:173-177.[Abstract/Free Full Text]
-
Hennig J, Nauerth A, Friedburg H. RARE imaging: a fast imaging method for clinical MR. Magn Reson Med 1986; 3:823-833.[Medline]
-
Mulkern RV, Wong STS, Winalski C, Jolesz FA. Contrast manipulation and artifact assessment of 2D and 3D RARE sequences. Magn Reson Imaging 1990; 8:557-566.[Medline]
-
Melki PS, Jolesz FA, Mulkern RV. Partial RF echo planar imaging with the FAISE method. I. Experimental and theoretical assessment of artifact. Magn Reson Med 1992; 26:328-341.[Medline]
-
Jones KM, Mulkern RV, Schwartz RB, Oshio K, Barnes PD, Jolesz FA. Fast spin-echo MR imaging of the brain and spine: current concepts. AJR Am J Roentgenol 1992; 158:1313-1320.[Free Full Text]
-
Oshio K, Jolesz FA, Melki PS, Mulkern RV. T2-weighted thin-section imaging with the multislab three-dimensional RARE technique. J Magn Reson Imaging 1991; 1:695-700.[Medline]
-
Yuan C, Schmiedl UP, Weinberger E, Krueck WR, Rand SD. Three-dimensional fast spin-echo imaging: pulse sequence and in vivo image evaluation. J Magn Reson Imaging 1993; 3:894-899.[Medline]
-
Melki PS, Mulkern RV. Magnetization transfer effects in multislice RARE sequences. Magn Reson Med 1992; 24:189-195.[Medline]
-
Mugler JP, III, Brookeman JR. Off-resonance image artifacts in interleaved-EPI and GRASE pulse sequences. Magn Reson Med 1996; 36:306-313.[Medline]
-
Maudsley AA. Dynamic range improvement in NMR imaging using phase scrambling. J Magn Reson 1988; 76:287-305.
-
Wedeen VJ, Chao YS, Ackerman JL. Dynamic range compression in MRI by means of a nonlinear gradient pulse. Magn Reson Med 1988; 6:287-295.[Medline]
-
Oh CH, Hilal SK, Wu EX, Cho ZH. Phase-scrambled RF excitation for 3D volume-selective multislice NMR imaging. Magn Reson Med 1992; 28:290-299.[Medline]
-
Delsuc MA, Lallemand JY. Improvement of dynamic range in NMR by oversampling. J Magn Reson 1986; 69:504-507.
-
Oh CH, Hilal SK, Johnson G, Mun IK, Javaid T, Cho ZH. Effective dynamic range improvement of NMR signal detection by using analog programmable attenuators (abstr). Book of abstracts: Society of Magnetic Resonance in Medicine 1989; 181.
-
Elliott MA, Insko EK, Greenman RL, Leigh JS. Improved resolution and signal-to-noise ratio in MRI via enhanced signal digitization. J Magn Reson 1998; 130:300-304.[Medline]
-
Norris DG, Bornert P, Reese T, Leibfritz D. On the application of ultra-fast RARE experiments. Magn Reson Med 1992; 27:142-164.[Medline]
-
Norris DG, Bornert P. Coherence and interference in ultrafast RARE experiments. J Magn Reson A 1993; 105:123-127.
-
Margosian P, Schmitt F, Purdy D. Faster MR imaging: imaging with half the data. Health Care Instrum 1986; 1:195-197.
-
Wielopolski P, Achten E, de Groot JC, et al. A dual echo three-dimensional inversion recovery HASTE sequence: initial experience and potential applications for the quantification of white matter lesions and brain tissue volumes (abstr) In: Proceedings of the Fourth Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1996; 554.
-
Alsop DC. The sensitivity of low flip angleRARE imaging. Magn Reson Med 1997; 37:176-184.[Medline]
-
Barker GJB. 3D fast FLAIR: a CSF-nulled 3D fast spin-echo pulse sequence. Magn Reson Imaging 1998; 16:715-720.[Medline]
-
Rydberg JN, Hammond CA, Grimm RC, et al. Initial clinical experience in MR imaging of the brain with a fast fluid-attenuated inversion-recovery pulse sequence. Radiology 1994; 193:173-180.[Abstract/Free Full Text]
-
Brant-Zawadzki M, Gillan GD, Nitz WR. MP RAGE: a three-dimensional, T1-weighted, gradient-echo sequenceinitial experience in the brain. Radiology 1992; 182:769-775.[Abstract/Free Full Text]
-
Chappell PM, Pelc NJ, Foo TK, Glover GH, Haros SP, Enzmann DR. Comparison of lesion enhancement on spin-echo and gradient-echo images. AJNR Am J Neuroradiol 1994; 15:37-44.[Abstract]
-
Feinberg DA, Oshio K. GRASE (gradient- and spin-echo) MR imaging: a new fast clinical imaging technique. Radiology 1991; 181:597-602.[Abstract/Free Full Text]
-
Mugler JP, III. Improved three-dimensional GRASE imaging with the SORT phase-encoding strategy. J Magn Reson Imaging 1999; 9:604-612.[Medline]
-
Le Roux P, Hinks RS. Stabilization of echo amplitudes in FSE sequences. Magn Reson Med 1993; 30:183-190.[Medline]
-
Le Roux P, Sun Y. Slab boundary artefact reduction in 3DFSE (abstr) In: Proceedings of the Fourth Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1996; 112.
-
Liu K, Xu Y, Loncar M. Artifact transformation technique: shifted interleaved multi-volume acquisition (SIMVA) for 3D FSE (abstr) In: Proceedings of the Sixth Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1998; 572.
This article has been cited by other articles:

|
 |

|
 |
 
S. Roosendaal, B Moraal, P. Pouwels, H Vrenken, J. Castelijns, F Barkhof, and J. Geurts
Accumulation of cortical lesions in MS: relation with cognitive impairment
Multiple Sclerosis,
June 1, 2009;
15(6):
708 - 714.
[Abstract]
[PDF]
|
 |
|

|
 |

|
 |
 
Y. Kato, S. Higano, H. Tamura, S. Mugikura, A. Umetsu, T. Murata, and S. Takahashi
Usefulness of Contrast-Enhanced T1-Weighted Sampling Perfection with Application-Optimized Contrasts by Using Different Flip Angle Evolutions in Detection of Small Brain Metastasis at 3T MR Imaging: Comparison with Magnetization-Prepared Rapid Acquisition of Gradient Echo Imaging
AJNR Am. J. Neuroradiol.,
May 1, 2009;
30(5):
923 - 929.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. Moraal, D. S. Meier, P. A. Poppe, J. J. G. Geurts, H. Vrenken, W. M. A. Jonker, D. L. Knol, R. A. van Schijndel, P. J. W. Pouwels, C. Pohl, et al.
Subtraction MR Images in a Multiple Sclerosis Multicenter Clinical Trial Setting
Radiology,
February 1, 2009;
250(2):
506 - 514.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
K. J. Stevens, R. F. Busse, E. Han, A. C. S. Brau, P. J. Beatty, C. F. Beaulieu, and G. E. Gold
Ankle: Isotropic MR Imaging with 3D-FSE-Cube--Initial Experience in Healthy Volunteers
Radiology,
December 1, 2008;
249(3):
1026 - 1033.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. K. Kuhl, F. Traber, and H. H. Schild
Whole-Body High-Field-Strength (3.0-T) MR Imaging in Clinical Practice * Part I. Technical Considerations and Clinical Applications
Radiology,
March 1, 2008;
246(3):
675 - 696.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
Y. Duan, P.G. Hildenbrand, M.P. Sampat, D.F. Tate, I. Csapo, B. Moraal, R. Bakshi, F. Barkhof, D.S. Meier, and C.R.G. Guttmann
Segmentation of Subtraction Images for the Measurement of Lesion Change in Multiple Sclerosis
AJNR Am. J. Neuroradiol.,
February 1, 2008;
29(2):
340 - 346.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
G. E. Gold, R. F. Busse, C. Beehler, E. Han, A. C. S. Brau, P. J. Beatty, and C. F. Beaulieu
Isotropic MRI of the Knee with 3D Fast Spin-Echo Extended Echo-Train Acquisition (XETA): Initial Experience
Am. J. Roentgenol.,
May 1, 2007;
188(5):
1287 - 1293.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R J Mills, C A Young, and E T S Smith
3D MRI in multiple sclerosis: a study of three sequences at 3 T
Br. J. Radiol.,
May 1, 2007;
80(953):
307 - 320.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. J. W. Pouwels, J. P. A. Kuijer, J. P. Mugler III, C. R. G. Guttmann, and F. Barkhof
Human Gray Matter: Feasibility of Single-Slab 3D Double Inversion-Recovery High-Spatial-Resolution MR Imaging
Radiology,
December 1, 2006;
241(3):
873 - 879.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. J. G. Geurts, P. J. W. Pouwels, B. M. J. Uitdehaag, C. H. Polman, F. Barkhof, and J. A. Castelijns
Intracortical Lesions in Multiple Sclerosis: Improved Detection with 3D Double Inversion-Recovery MR Imaging
Radiology,
July 1, 2005;
236(1):
254 - 260.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
I L. Tan, R. A. van Schijndel, P. J. W. Pouwels, H. J. Ader, and F. Barkhof
Serial Isotropic Three-Dimensional Fast FLAIR Imaging: Using Image Registration and Subtraction to Reveal Active Multiple Sclerosis Lesions
Am. J. Roentgenol.,
September 1, 2002;
179(3):
777 - 782.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. F. Kallmes, F. K. Hui, and J. P. Mugler III
Suppression of Cerebrospinal Fluid and Blood Flow Artifacts in FLAIR MR Imaging with a Single-Slab Three-dimensional Pulse Sequence: Initial Experience
Radiology,
October 1, 2001;
221(1):
251 - 255.
[Abstract]
[Full Text]
[PDF]
|
 |
|