Properties of the phase-alternating phase-shift (PHAPS) multiple spin-echo protocol in MRI: a study of the effects of imperfect RF pulses

Multiexponential T2 relaxometry of benign and malignant adipocytic tumours

Background

We investigated a recently proposed multiexponential (Mexp) fitting method applied to T2 relaxometry magnetic resonance imaging (MRI) data of benign and malignant adipocytic tumours and healthy subcutaneous fat. We studied the T2 distributions of the different tissue types and calculated statistical metrics to differentiate benign and malignant tumours.

Methods

Twenty-four patients with primary benign and malignant adipocytic tumours prospectively underwent 1.5-T MRI with a single-slice T2 relaxometry (Carr-Purcell-Meiboom-Gill sequence, 25 echoes) prior to surgical excision and histopathological assessment. The proposed method adaptively chooses a monoexponential or biexponential model on a voxel basis based on the adjusted R² goodness of fit criterion. Linear regression was applied on the statistical metrics derived from the T2 distributions for the classification.

Results

Healthy subcutaneous fat and benign lipoma were better described by biexponential fitting with a monoexponential and biexponential prevalence of 0.0/100% and 0.2/99.8% respectively. Well-differentiated liposarcomas exhibit 17.6% monoexponential and 82.4% biexponential behaviour, while more aggressive liposarcomas show larger degree of monoexponential behaviour. The monoexponential/biexponential prevalence was 47.6/52.4% for myxoid tumours, 52.8/47.2% for poorly differentiated parts of dedifferentiated liposarcomas, and 24.9/75.1% pleomorphic liposarcomas. The percentage monoexponential or biexponential model prevalence per patient was the best classifier distinguishing between malignant and benign adipocytic tumours with a 0.81 sensitivity and a 1.00 specificity.

Conclusions

Healthy adipose tissue and benign lipomas showed a pure biexponential behaviour with similar T2 distributions, while decreased adipocytic cell differentiation characterising aggressive neoplasms was associated with an increased rate of monoexponential decay curves, opening a perspective adipocytic tumour classification.

Inverse Laplace transform and multiexponential fitting analysis of T2 relaxometry data: a phantom study with aqueous and fat containing samples

BACKGROUND

The inverse Laplace transform (ILT) is the most widely used method for T2 relaxometry data analysis. This study examines the qualitative agreement of ILT and a proposed multiexponential (Mexp method) regarding the number of T2 components. We performed a feasibility study for the voxelwise characterisation of heterogeneous tissue with T2 relaxometry.

METHODS

Eleven samples of aqueous, fatty and mixed composition were analysed using ILT and Mexp. The phantom was imaged using a 1.5-T system with a single slice T2 relaxometry 25-echo Carr-Purcell-Meiboom-Gill sequence in order to obtain the T2 decay curve with 25 equidistant echo times. The adjusted R² goodness of fit criterion was used to determine the number of T2 components using the Mexp method on a voxel-based analysis. Comparison of mean and standard deviation of T2 values for both methods was performed by fitting a Gaussian function to the ILT resulting vector.

RESULTS

Phantom results showed pure monoexponential decay for acetone and water and pure biexponential behaviour for corn oil, egg yolk, and 35% fat milk cream, while mixtures of egg whites and yolks as well as milk creams with 12-20% fatty composition exhibit mixed monoexponential and biexponential behaviour at different fractions. The number of T2 components by the Mexp method was compared to the ILT-derived spectrum as ground truth.

CONCLUSIONS

Mexp analysis with the adjusted R² criterion can be used for the detection of the T2 distribution of aqueous, fatty and mixed samples with the added advantage of voxelwise mapping.

T2, T2* and spin coupling ratio as biomarkers for the study of lipomatous tumors

BACKGROUND

Subcutaneous fat may have variable signal intensity on T2w images depending on the choice of imaging parameters. However, fatty components within tumors have a different degree of signal dependence on the acquisition scheme. This study examined the use of T2, T2^* relaxometry and spin coupling related signal changes (Spin Coupling ratio, SCr) on two different imaging protocols as clinically relevant descriptors of benign and malignant lipomatous tumors.

MATERIALS AND METHODS

20 patients with benign lipomas or liposarcomas of variable histologic grade were examined at an 1.5 T scanner with Multi Echo Spin Echo (MESE) different echo spacing (ESP) in order to produce bright fat T2w images (ESP: 13.4 ms, 25 equidistant echoes) and dark fat images (ESP: 26.8 ms with 10 equidistant echoes). T2^* relaxometry acquisition comprises 4 sets of in-opposed echoes (2.4-19.2 ms, ESP: 2.4 ms) Multi Echo Gradient Echo (MEGRE) sequence. All parametric maps were calculated on a pixel basis.

RESULTS

Significant differences of SCr were found for five different types of lipomatous tumors (Pairwise t-test with Bonferroni correction): lipomas, well differentiated liposarcomas, myxoid liposarcomas, pleomorphic liposarcomas and poorly differentiated liposarcomas. SCr surpassed the classification performance of T2 and T2^* relaxometry.

DATA CONCLUSION

A novel biomarker based on spin coupling related signal loss, SCr, is indicative of lipomatous tumor histological grading. We concluded that T2, T2^* and SCr can be used for the classification of fat containing tumors, which may be important for biopsy guidance in heterogeneous masses and treatment planning.

Texture Analysis in Quantitative MR Imaging

The diagnostic potential of texture analysis in quantitative tissue characterisation by MR imaging at 1.5 T was evaluated in the brain of 6 healthy volunteers and in 88 patients with intracranial tumours. Texture images were computed from calculated T1 and T2 parameter images by applying groups of common first-order and second-order grey level statistics. Tissue differentiation in the images was estimated by the presence or absence of significant differences between tissue types. A fine discrimination was obtained between white matter, cortical grey matter, and cerebrospinal fluid in the normal brain, and white matter was readily separated from the tumour lesions. Moreover, separation of solid tumour tissue and peritumoural oedema was suggested for some tumour types. Mutual comparison of all tumour types revealed extensive differences, and even specific tumour differentiation turned out to be successful in some cases of clinical importance. However, no discrimination between benign and malignant tumour growth was possible. Much texture information seems to be contained in MR images, which may prove useful for classification and image segmentation.

Compressed sensing CPMG with group-sparse reconstruction for myelin water imaging

PURPOSE

Myelin content is a marker for nervous system pathology and is quantifiable by myelin water imaging using multi-echo CPMG sequence, which is inherently slow. One way to accelerate the scan is to utilize compressed sensing. However, reconstructing the images piecemeal by standard compressed sensing methods is not the optimal solution, because it only exploits intraimage spatial redundancy. It does not recognize that the different T2 weighted images are scans of the same anatomical volume and hence correlated. The purpose of this work is to test the feasibility of compressed sensed CPMG with group-sparsity promoting optimization for myelin water imaging.

METHODS

Group-sparse reconstruction was performed at various simulated and actual undersampling factors for an electronic phantom, ex vivo rat spinal cord, and in vivo rat spinal cord. Normalized mean square error was used as the metric for comparison.

RESULTS

For both simulated undersampling and the actual undersampling, the method was found to minimally impact myelin water fraction map quality (normalized mean square error < 0.25) when acceleration factor was below two.

CONCLUSION

Compressed sensed CPMG with group-sparse reconstruction is useful for achieving a shorter scan time than traditionally possible.

On electrophoretic NMR. Exploring high conductivity samples

The performance of a new electrophoretic NMR (eNMR) method that uses a Carr-Purcell-Meiboom-Gill echo train with repeated electric field reversal is investigated. We show that this pulse sequence, with acronym CPMGER, yields strongly reduced artifacts from convective flow effects caused by the simultaneous presence of electroosmotic and thermal driving forces. We demonstrate the achieved improvements in various aqueous solutions. Ultimately, the method can be used for obtaining electrophoretic mobilities by eNMR without relying on uncharged reference molecules, otherwise a significant limitation for electrophoretic experiments performed with nuclei other than (1)H.

Sensitivity-enhanced chemical exchange saturation transfer (CEST) MRI with least squares optimization of Carr Purcell Meiboom Gill multi-echo echo planar imaging

Chemical exchange saturation transfer (CEST) imaging is a novel MRI technique that is sensitive to biomolecules, local pH and temperature, and offers considerable advantages for in vivo applications. However, the magnitude of CEST effect for dilute CEST agents undergoing slow or intermediate chemical exchange is typically small, requiring the use of signal averaging to enhance its sensitivity. Given that T2 -induced signal loss can be normalized by asymmetry analysis, the magnitude of CEST effect is independent of echo time. Therefore, CEST MRI with multi-echo echo planar imaging (EPI) readout should yield the same CEST effect as conventional single echo acquisition. Importantly, CEST multi-echo (CESTme) EPI images can be averaged to enhance CEST MRI sensitivity. The goal of this study was to validate CESTme EPI using a creatine-agarose gel CEST phantom with similar T2 as biological tissue. Using least-squares optimization, we found that the sensitivity of CESTme sequence was significantly higher than that obtained by conventional single echo CEST-EPI acquisition. Specifically, signal-to-noise ratio and contrast-to-noise ratio from the proposed CESTme EPI were approximately equivalent to that obtained by doubling the number of signal averages of the standard single echo CEST MRI sequence. In summary, our results demonstrated CESTme EPI for sensitivity-enhanced CEST imaging.

Effects of MR image noise on estimation of short T2 values from T2 -weighted image series

We examined how the noise from magnetic resonance (MR) imaging affects the calculation of T(2) in skeletal muscle, a tissue with short T(2) values. The measured pixel intensity of the MR image (: the magnitude image) was the superimposed signal which was composed of the MR signal and the noise, and we demonstrated that noise from a magnitude image matches the DC component of the T(2) decay curve. In materials with long T(2) values, the noise has no influence on the selective echo time (TE) in calculating T(2). However, in materials with short T(2) values, noise clearly influences the selective TE. In this study, we proposed a T(2) effective signal-ratio, T(2)SR, as an index for determining whether the noise of the magnitude image can be ignored in calculating T(2). When T(2)SR and the signal-to-noise ratio (SNR), an index of image quality, were compared as indices to evaluate the influence of noise in the calculation of T(2), T(2)SR was useful and SNR was not. The use of multiple spin echo (MSE) technique shortened imaging time, but required detailed understanding of the MSE. Our results indicated that T(2) can be calculated correctly for skeletal muscle and other tissues with short T(2) even when the receiver coil has a low SNR and few measurement points are available.

Optimized clinical T2 relaxometry with a standard CPMG sequence

PURPOSE

To optimize the accuracy and precision of T2 measurements using the standard Carr-Purcell-Meiboom-Gill (CPMG) sequence. T2 values obtained with this technique are normally sensitive to imperfect refocusing due to the formation of unwanted stimulated echoes.

MATERIALS AND METHODS

Modifications are made to the refocusing slice selection width and the interleaving scheme. A widened refocusing slice improves the uniformity of the refocusing flip angle across the slice. A slow spin echo acquisition provided "gold standard" T2 values. Repeated T2 measurements in phantom and human studies are used to compare the accuracy and precision of the optimized and non-optimized CPMG implementations.

RESULTS

The required slice thickness ratio between refocusing and excitation slice widths is found to be 3:1 for typical optimized radiofrequency pulses. T2 values obtained using this optimized implementation more closely correspond to "gold standard" values. Repeated T2 measurements indicate significantly improved correspondence between data and model. A reduction in the fitting error of approximately 70% is demonstrated for phantoms.

CONCLUSION

We demonstrate that a relatively simple change to the CPMG relaxometry sequence parameters from the default setup yields significant improvements in the accuracy and precision of T2 measurements.

Investigation of vacuum pumping on the dose response of the MAGAS normoxic polymer gel dosimeter

The effect of vacuum pumping on the dose response of the MAGAS polymer gel dosimeter has been investigated. A delay of several days post-manufacture before irradiation was previously necessary due to the slow oxygen scavenging of ascorbic acid. The MAGAS polymer gel dosimeter was vacuum pumped before gelation to remove dissolved oxygen. The MAGAS polymer gel dosimeter was poured into glass screw-top vials, which were irradiated at various times, post-manufacture to a range of doses. Magnetic resonance imaging techniques were used to determine the R2-dose response and R2-dose sensitivity of the MAGAS polymer gel. The results were compared with a control batch of MAGAS polymer gel that was not vacuum pumped. It was shown that vacuum pumping on the MAGAS polymer gel solution immediately prior to sealing in glass screw-top vials initially increases the R2-dose response and R2-dose sensitivity of the dosimeter. An increase in the R2-dose response and R2-dose sensitivity was observed with increasing time between manufacture and irradiation. Over the range of post-manufacture irradiation times investigated, the greatest R2-dose response and R2-dose sensitivity occurred at 96 hours.

Investigation of vacuum pumping on the dose response of the MAGAS normoxic polymer gel dosimeter

The effect of vacuum pumping on the dose response of the MAGAS polymer gel dosimeter has been investigated. A delay of several days post-manufacture before irradiation was previously necessary due to the slow oxygen scavenging of ascorbic acid. The MAGAS polymer gel dosimeter was vacuum pumped before gelation to remove dissolved oxygen. The MAGAS polymer gel dosimeter was poured into glass screw-top vials, which were irradiated at various times, post-manufacture to a range of doses. Magnetic resonance imaging techniques were used to determine the R2-dose response and R2-dose sensitivity of the MAGAS polymer gel. The results were compared with a control batch of MAGAS polymer gel that was not vacuum pumped. It was shown that vacuum pumping on the MAGAS polymer gel solution immediately prior to sealing in glass screw-top vials initially increases the R2-dose response and R2-dose sensitivity of the dosimeter. An increase in the R2-dose response and R2-dose sensitivity was observed with increasing time between manufacture and irradiation. Over the range of post-manufacture irradiation times investigated, the greatest R2-dose response and R2-dose sensitivity occurred at 96 hours.

Comparison of TSE, TGSE, and CPMG measurement techniques for MR polymer gel dosimetry

A radiation dose distribution that optimally conforms to the target volume is of major interest for stereotactic radiotherapy. For this purpose treatment plans have to be verified experimentally before transferring to the patient. The requirements regarding dose accuracy and spatial resolution can be fulfilled with tissue equivalent polymer gel dosimeters which offer the possibility to visualize 3D dose distributions. Herewith, dosimetry can be performed by the spin-spin relaxation rate R2 which varies with the absorbed dose. In this work, different MR measurement techniques were evaluated: The standard Carr-Purcell-Meiboom-Gill (CPMG) method, a modified Turbo-Spin-Echo (TSE) sequence, and a modified Turbo-Gradient-Spin-Echo (TGSE) sequence. Experiments were performed both with a homogeneous water phantom and an irradiated polymer gel. The results show that TGSE and especially TSE are suited well for MR polymer gel dosimetry: The acquisition time of both techniques can be reduced in comparison to CPMG by a factor of 5. The accuracy of dose determination for doses between 2 Gy and 13 Gy lies between 5.6% and 2.0% (TSE), 9.0% and 3.2% (TGSE), and 7.9% and 2.7% (CPMG). These investigations show that especially TSE can be handled as a substitute or at least an alternative to CPMG for the verification of treatment plans in stereotactic radiotherapy.

Double inversion black-blood fast spin-echo imaging of the human heart: a comparison between 1.5T and 3.0T

PURPOSE

To evaluate the effectiveness of blood suppression and the quality of black-blood cardiac images acquired at 3.0 Tesla using a double-inversion recovery fast spin-echo sequence by comparing data acquired at 3.0T to data acquired at 1.5T.

MATERIALS AND METHODS

Black-blood T2-weighted fast spin-echo images of the heart were acquired from five normal volunteers at 1.5T and five normal volunteers at 3.0T. Region-of-interest signal intensity measurements were performed at several locations in the suppressed blood regions of the left and right ventricles and around the left ventricle walls to assess the effectiveness and uniformity of the blood suppression, the myocardial signal-to-noise ratio (SNR), and the signal uniformity at both field strengths. B1 field maps were produced in phantoms and in subjects at both field strengths.

RESULTS

Blood suppression performance is equivalent at 1.5T and 3.0T. The improvement in SNR at 3.0T compared with 1.5T is less than has been predicted in previous studies. The signal uniformity is significantly poorer at 3.0T than at 1.5T due to dielectric effects and shorter radio frequency wavelengths (P < 0.005).

CONCLUSION

Spin-echo and spin-echo echo-train sequences that perform well at 1.5T will produce large signal variations in the chest cavity at 3.0T without modifications. B1 insensitive methods must be explored and implemented for spin-echo sequences to fully realize the advantages of using these sequences for high-field MRI.

Artefacts in multi-echo T2 imaging for high-precision gel dosimetry: III. Effects of temperature drift during scanning

In high-precision 3D gel dosimetry, long MR measurement times together with a high amount of RF energy being absorbed by the phantom are very common, and result in a spatially dependent temperature rise in the gel. As T2 of the dosimeter gel is temperature dependent, dose estimation will be affected. In this study we assess the temperature rise in the dosimeter gel by use of MR temperature mapping and computer modelling. It is shown that in conventional MR sequences. where linear k-space sampling is used, a temperature rise of 3 C results in a dose underestimation of 10% over the whole dose map. To correct for these dose errors, a compensation method involving centric k-space ordering is suggested. Computer simulations have been performed to analyse the robustness of the proposed method. Applying the compensated sequence, a temperature rise of 3 C leads to a narrow dose artefact of the order of 3% for a 'worst case' situation in which a single pixel dose gradient is assumed. Negligible deviations are found in the rest of the dose map.

High-resolution polymer gel dosimetry by parameter selective MR-microimaging on a whole body scanner at 3T

High dose variations across small spatial distances, as present in brachytherapeutic applications or radiosurgery and especially gamma-knife therapy, are difficult to quantify by standard dosimetry. We demonstrate the possibility to obtain planar spatial resolutions for dose imaging at pixel sizes below 200 microm within multislice parameter selective MR imaging on polymer gels. The sensitivity of the transversal and longitudinal relaxation time as well as diffusivity on dose is shown. High spatial resolution is achieved by parameter selective microimaging of polymer gels on a high-field (3 T) whole-body MR system equipped with a dedicated strong gradient system and a small probe head matched to the sample size. In addition to the spin-spin relaxation rate R2 = 1/T2 we investigate the sensitivity of the longitudinal relaxation rate R1 = 1/T1 and the diffusivity Dapp in acrylic polymer gels on irradiation up to dose levels of about 20 Gy. Dose images are obtained after calibration of the corresponding MR parameters by known dose levels of gamma irradiation. Also the MR-parameter T1 may be used for dose imaging. The impact of all of the three parameters T1, T2, and diffusivity on obtained signal intensities in irradiated regions has to be taken into account in nonoptimized pulse sequences. Further, very high spatial resolution imposes several restrictions on the evaluation of R2, which have to be considered for quantitative dosimetry. These restrictions are discussed in detail. We also demonstrate the importance of such a high spatial resolution in case of a set of differently sized gamma-knife stereotactic irradiation schemes. Gel dosimetry based on MR parameter selective microimaging represents a potent alternative for the detection of dose distributions characterized by steep dose gradients, typical in brachytherapeutic and radiosurgical applications.

Artefacts in multi-echo T2 imaging for high-precision gel dosimetry: I. Analysis and compensation of eddy currents

In BANG gel dosimetry, the spin-spin relaxation rate, R2 = I/T2, is related to the radiation dose that has been delivered to the gel phantom. R2 is calculated by fitting the pixel intensities of a set of differently T2-weighted base images. In gel dosimetry for radiotherapy, an accuracy of 5% in dose and 3 mm spatially, whichever is lower, is the objective. Therefore, possible sources of artefacts must be considered and dealt with. To obtain a set of base images a multiple spin-echo sequence is used. However, in a conventional MR scanner eddy currents will be provoked by switching the imaging gradients. As the eddy currents change in the course of the sequence, the net magnetization will be affected accordingly. Hence, eddy currents may have a significant influence on the quantitative R2 images themselves as well as on their slice position. In this study, we report an analysis of the eddy currents as they appear in the multiple spin-echo sequence. Eddy currents are measured using a frequency shift method resulting in eddy current field maps. The related geometrical displacements are obtained by use of a pyramidal phantom. The R2 versus dose relation is determined in the three main directions of the magnet, revealing a dependence of the measured R2 on slice orientation. The time course of eddy currents is then used in a computer simulation to estimate the effects they produce in the recorded R2 images. A compensation method for eddy current effects in multi-echo T2 mapping is described.

Artefacts in multi-echo T2 imaging for high-precision gel dosimetry: II. Analysis of B1-field inhomogeneity

In BANG gel dosimetry, the spin-spin relaxation rate, R2 = 1/T2, is related to radiation dose that has been delivered to a gel phantom. R2 is calculated by fitting the pixel intensities of a set of differently T2-weighted base images. The accuracy that is aimed for in this quantitative MR application is about 5% relative to the maximum dose. In a conventional imaging MR scanner, however, several imaging artefacts may perturb the final dose map. These deviations manifest themselves as either a deformation of the dose map or an inaccuracy of the dose pixel value. Inaccuracies in the dose maps are caused by both spatial and temporal deviations in signal intensities during scanning. This study deals with B1-field inhomogeneities as a source of dose inaccuracy. First, the influence of B1-field inhomogeneities on slice profiles is investigated using a thin-slice phantom. Secondly, a FLASH sequence is used to map the B1-field by assessing the effective flip angle in each voxel of a homogeneous phantom. In addition, both experiments and computer simulations revealed the effects of B1 field inhomogeneities on the measured R2. This work offers a method to correct R2 maps for B1 -field inhomogeneities.

High-resolution dose profile studies based on MR imaging with polymer BANG(TM) gels in stereotactic radiation techniques

High-resolution dose profiles produced by the Leksell Gamma Knife were obtained in BANG(TM) polymer gel, using a 3 T whole-body scanner upgraded by a magnetic resonance microscopy unit. The gel was contained in 22.3 mm diameter flasks that were inserted into a solid, tissue-equivalent head phantom irradiated by fields of by 8 and 14 mm collimators. Dose profiles were obtained from a linear dose-response curve (R(2) vs. Dose). Excellent agreement was obtained when the gel data were compared to film dosimetry and calculated data.

An investigation of the chemical stability of a monomer/polymer gel dosimeter

The aim of this work is to investigate the temporal stability of a polyacrylamide gelatin hydrogel used for 3D monomer/polymer gel dosimetry techniques involving different methods of analysis. Long-term instabilities for a similar gel have recently been reported, but differ markedly from those described in this work. Two kinds of long-term instabilities are described. One affects the slope of the dose-R2 plot and is related to post-irradiation polymerization of the comonomer/polymer aggregates. It is observed that post-irradiation polymerization only lasts 12 hours after irradiation. The other instability affects the intercept of the dose-R2 plot, lasts for up to 30 days and is related to the gelation process of gelatin. Further studies were performed on gelatin gels of varying compositions to obtain a better understanding of the molecular mechanism that causes the instability due to gelation. The studies included observations of the spin-spin and spin-lattice relaxation rates in combination with diffusion measurements and optical measurements. It is shown that the heating history during the manufacture of the gel affects the absolute R2 value of the gel but not its variation. The findings presented in this study may help in producing more stable and reproducible monomer/polymer gel dosimeters.

Validation of MR-based polymer gel dosimetry as a preclinical three-dimensional verification tool in conformal radiotherapy

The aim of this work was to investigate MR-based polymer gel dosimetry as a three-dimensional (3D) dosimetry technique in conformal radiotherapy. A cylindrical container filled with polymer gel was placed in a water-filled torso phantom to verify a treatment plan for the conformal irradiation of a mediastinal tumor located near the esophagus. Magnetic resonance spin-spin relaxation rate images were acquired and, after calibration, converted to absorbed dose distributions. The dose maps were compared with dose distributions measured using radiographic film. The average root-mean-square structural deviation, for the complete dose distribution, amounted to less than 3% between gel and film dose maps. It may be expected that MR gel dosimetry will become a valuable tool in the verification of 3D dose distributions. The influence of imaging artifacts arising from eddy currents, temperature drift during scanning, and B1 field inhomogeneity on the dose maps was taken into account and minimized.

Non-Fourier encoding with multiple spin echoes

The advantages and limitations of multiple spin-echo sequences for non-Fourier encoding are investigated. Complications caused by improper encoding of alternate magnetization pathways due to imperfect refocusing pulses are analyzed. It is shown that mirror image ghosts result if the encoding RF pulse matrix is real-valued. These ghosts can be avoided as long as the rows of the RF pulse matrix are conjugate symmetric, which implies that spatial profiles are real valued. Non-Fourier encoding using bases derived from wavelet, Hadamard, and other real-valued orthogonal functions does not result in a mirror ghost artifact. A RARE sequence for non-Fourier encoding has been implemented on a clinical imaging system and successfully applied for brain imaging.

T2 maximum likelihood estimation from multiple spin-echo magnitude images

An optimal maximum likelihood (ML) method is described for an unbiased estimation of monoexponential T2 from magnitude spin-echo images. The algorithm is based on a Gaussian assumption of noise distribution. The validity of this assumption was checked by a statistical chi 2 test on spin-echo and fast low-angle shot surface coil images. Monte-Carlo simulations of magnitude data showed that the ML estimate standard deviation is lower than that produced by a weighted least-squares fitting on signal logarithm. Correction schemes are proposed to reduce bias deriving from magnitude reconstruction. The variance of the ML estimate converged rapidly toward the theoretical algebraic expression of the Cramér-Rao lower bound.

Characterization of parotid gland tissue: a description of an MRI protocol set-up and results of in-vivo applications

A reliable protocol for proton T2 mapping of the parotid region was set up for future characterization of parotid gland disease. A Carr-Purcell-Meiboom-Gill sequence, phase compensated, available on our 1.5 T imager, was selected and acquisition parameters were chosen on the basis of tests performed on phantoms (agarose-doped gels with T2 in the physiological range). Some experiments were carried out to evaluate the accuracy of T2 calculations for selective and nonselective refocussing pulses, for image uniformity corrections, and for different situations of slice shift and repetition times. The chosen protocol was then applied to in vivo evaluations to check the long-term precision by means of repeated measurements performed on the same subject over a 2-month period. Two or more reference gels were positioned both in the phantom and volunteer at the edge of the field-of-view (FOV). Image postprocessing consisted of an automatic procedure, written by the authors in Fortran 77, that selected the best fit for each pixel between mono- and biexponential decay models, and prepared four parametric images (T2 and Rho slow and fast contribution, Rho being a function of proton density and of T1) that may be used for future elaborations. The phantom experiment results showed an accuracy of 2.5% if a linear correction was performed using the reference gels at the edge of the FOV. No significant differences in accuracy were found between selective and nonselective refocussing pulse, and a homogeneity correction was not demonstrated necessary. The measurements performed on four volunteers showed that the best decaying model for healthy parotid tissue was monoexponential. Evaluated T2 resulted 80.18 +/- 6.11 ms (72.96 +/- 4.97 ms for uncorrected results). Long-term reproducibility of the group of measurements from one volunteer, summarizing all the measurement errors, ranged from 0.9 to 8.5%. The two-way ANOVA that was carried out considering the two classes of volunteers and of parotid positions (right or left) showed that differences found between the two parotids were not significant, while T2 differences among individuals are significant if a probability level higher than 1.1% is accepted. As in this case, the main source of error can be attributed to the biological variations among individuals. Future statistics collected on patients for the T2 evaluations of the pathologic tissue will clarify whether the T2 relaxation is a sufficient parameter for T2 discrimination of healthy and pathologic tissue.

Dependence on T1 of the echo amplitudes from multiple spin-echo sequences with equidistant echoes: simulation studies

Transverse relaxation times were estimated from numerical simulations on spin systems using multi-echo spin-echo MRI protocols. The influence of T1 on the echo amplitudes via stimulated echo components was studied. The resulting effects on T2 estimates from the Carr-Purcell (CP), Carr-Purcell-Meiboom-Gill (CPMG), and Phase-Alternating-Phase-Shift (PHAPS; combination of CP and CPMG), multiple echo schemes were examined. Protocols with either spatially selective or nonselective refocusing pulses were studied. An intravoxel static field inhomogeneity of 0.1, 1, and 10 ppm was stimulated. The dependence on T1 of the T2 estimates was notable for T1 values below approximately 800 msec for all protocols. The PHAPS scheme provided rather accurate, but underestimated, T2 values when selective refocusing was used. With nonselective refocusing, PHAPS T2 values were overestimated and demonstrated a pronounced dependence on magnetic field inhomogeneity. In general, long T2 values were erroneous with the PHAPS protocol. The results indicate that a CPMG protocol structure provides a more robust method for T2 estimations than the PHAPS protocol.

MRI measurements of the dependence on T1 of the echo amplitudes using a multiple spin-echo scheme

Analytical calculations using the Bloch formalism were performed to assess the dependence on T1 of the echo amplitudes for the Phase-Alternating Phase-Shift (PHAPS) multiple spin-echo protocol. Measurements in a 0.5 T MR imaging unit were performed to ratify the analytical results. Especially for low T2 values, the echo amplitudes were erroneous, with an increasing contribution from stimulated echo components with increasing T1. Apart from affecting T2 estimates, stimulated echoes generated a non-monoexponential signal decay of the echo trains. The results confirmed previous simulation studies as regards the dependence on T1 of T2 estimates from PHAPS.