MRI simulator for instruction in pulse-sequence selection

An ordinary desk-top microcomputer was programmed to simulate MR images for specified spin-echo pulse sequences. Model pixel maps of proton density and T1 and T2 relaxation times were made from published estimated values for regions of the human head, neck, and spine. Images were generated and displayed from the model maps and user-specified pulse-sequence parameters in less than 30 sec/image. Models for various pathologic conditions, including calcification, subacute hemorrhage, porencephaly, lipoma, and multiple sclerosis, were superimposed on the images of normal anatomy to create unknown cases. Simulated images can easily demonstrate the effect of pulse-sequence selection on the contrast of normal structures and pathologic conditions. Use of simulated images is an excellent technique for gaining experience in pulse-sequence selection. Low-cost microcomputers can provide adequate image detail and reasonable image display time of synthetic MR images for teaching purposes.     

Multispectral analysis of magnetic resonance images

Magnetic resonance (MR) imaging systems produce spatial distribution estimates of proton density, relaxation time, and flow, in a two dimensional matrix form that is analogous to that of the image data obtained from multispectral imaging satellites. Advanced NASA satellite image processing offers sophisticated multispectral analysis of MR images. Spin echo and inversion recovery pulse sequence images were entered in a digital format compatible with satellite images and accurately registered pixel by pixel. Signatures of each tissue class were automatically determined using both supervised and unsupervised classification. Overall tissue classification was obtained in the form of a theme map. In MR images of the brain, for example, the classes included CSF, gray matter, white matter, subcutaneous fat, muscle, and bone. These methods provide an efficient means of identifying subtle relationships in a multi-image MR study.     

Uterine leiomyoma: correlation between signal intensity on magnetic resonance imaging and pathologic characteristics.

To correlate the signal intensity of uterine leiomyoma with its pathologic characteristics, with particular emphasis on the fibrous component, 33 magnetic resonance (MR) examinations that revealed 93 leiomyomas were prospectively studied. All patients were imaged in axial and sagittal planes with different spin-echo pulse sequences to obtain T1-, T2-weighted, and proton density images. Nondegenerative leiomyomas (n = 62) showing a homogeneous signal of low intensity, and degenerative leiomyomas (n = 31) with a heterogeneous signal of variable intensity on T2-weighted images could be correlated. Histopathological assessment of fiber constitution and degeneration, and MR intensity were interpreted by independent observers. There was excellent accord between the averages for MR intensity, T2 relaxation time and fiber content, although the intensity values in each fiber grade showed a wide range. The greater the fiber content the lower the MR intensity on T2-weighted images, and the shorter the T2 relaxation time (p < 0.0001). In addition, the manner in which fiber distribution affected MR appearance was also elucidated. These data contribute guidelines for precise tissue differentiation of myogenic tumors on MR images, and for MR imaging tissue diagnosis of any lesion with a considerable fibrous element.     

Statistical error mapping for reliable quantitative T1 imaging

We developed a statistically based error image for rapid appreciation of unreliable regions in quantitative water proton T1 images. The chi-squared error and coefficient of variation of the fitted parameter were used to estimate uncertainties in the goodness-of-fit to monoexponential T1 relaxation and the reliability of the calculated T1, respectively, for each pixel. Errors exceeding a statistical threshold based on a .1 acceptance criterion were displayed as a color-coded overlay on the T1 image. Error maps of quantitative T1 images from 31 healthy volunteers showed a characteristic error structure; few pixels within the parenchyma had excessive errors. Clinical cases with stroke and sickle cell disease showed deviations from the normal pattern in the spatial distribution and magnitude of chisquared errors. Disease states may deviate from monoexponential T1 relaxation more than normal brain does. The color-coded error map is a valuable tool for investigators using quantitative MR imaging to determine tissue relaxation parameters.     

Reproducibility of relaxation and spin-density parameters in phantoms and the human brain measured by MR imaging at 1.5 T

The reproducibility of T1, T2, and proton density, measured in phantoms and the human brain was evaluated by proton imaging techniques. The sequence used to derive T1 and density values was a multiple-saturation recovery which consists of four pairs of 90 degrees pulses, followed by a 180 degrees phase reversal pulse, generating four T1-weighted images. T2 was derived from a multiple-echo sequence, generating four T2-weighted images. The data were analyzed by fitting the pixel intensities to the respective equations by means of nonlinear multiparameter least-squares analysis. Short-term reproducibility between four consecutive scans was evaluated to be 1-4% depending on location with a phantom covering the entire span of physiologic T1 and T2 values. A second phantom containing a series of identical samples served to study the dependence of the apparent T1 and T2 on position, both radially and axially, with respect to magnet isocenter. Reproducibility across the field of view was found to be better than 7% (T1 and T2). This phantom was further used to evaluate effects of long-term reproducibility, which at each location varied from 5-14% (T1) and 2-10% (T2). Finally, interinstrument reproducibility, tested by means of the same protocol on three different instruments, all operating at the same magnetic field and using largely identical hardware for each location, was found to be 1-14% (T1) and 2-10% (T2). The positional dependence of the apparent relaxation times appears to be systematic and may be due to variations in the effective field, caused by magnet and rf inhomogeneity. Finally, brain tissue relaxation and spin-density data were determined using the same protocol in 37 scans performed on 27 normal volunteers. The tissues analyzed were putamen, thalamus, caudate nucleus, centrum semiovale, internal capsule, and corpus callosum. Excellent accordance was further obtained between left and right hemispheres.     

MR imaging using stimulated echoes (STEAM)

The introduction of STEAM (stimulated echo acquisition mode) magnetic resonance (MR) sequences provides access to a variety of MR parameters. T1-weighted and calculated T1 proton MR images of the head of healthy volunteers and a patient with an astrocytoma are presented. MR examinations were performed with a 2.0-T whole-body system. The STEAM T1 method can be used to characterize multiexponential relaxation behavior, to evaluate T1 relaxation times, and to improve the T1 contrast within MR images. Both the measuring time and the spatial resolution are the same as for a conventional image.     

Model‐based Acceleration of Parameter mapping (MAP) for saturation prepared radially acquired data

A reconstruction technique called Model‐based Acceleration of Parameter mapping (MAP) is presented allowing for quantification of longitudinal relaxation time and proton density from radial single‐shot measurements after saturation recovery magnetization preparation. Using a mono‐exponential model in image space, an iterative fitting algorithm is used to reconstruct one well resolved and consistent image for each of the projections acquired during the saturation recovery relaxation process. The functionality of the algorithm is examined in numerical simulations, phantom experiments, and in‐vivo studies. MAP reconstructions of single‐shot acquisitions feature the same image quality and resolution as fully sampled reference images in phantom and in‐vivo studies. The longitudinal relaxation times obtained from the MAP reconstructions are in very good agreement with the reference values in numerical simulations as well as phantom and in‐vivo measurements. Compared to available contrast manipulation techniques, no averaging of projections acquired at different time points of the relaxation process is required in MAP imaging. The proposed technique offers new ways of extracting quantitative information from single‐shot measurements acquired after magnetization preparation. The reconstruction simultaneously yields images with high spatiotemporal resolution fully consistent with the acquired data as well as maps of the effective longitudinal relaxation parameter and the relative proton density. Magn Reson Med 70:1524–1534, 2013. © 2013 Wiley Periodicals, Inc.     

Multiparametric display of spin-echo data from MR studies of brain

A magnetic resonance (MR) image processing technique that uses a single color image for simultaneous presentation of spin-echo information and its application to MR studies of the brain is described. Relaxation rate and proton-density maps were calculated from 160 brain MR studies performed at 1.5 and 1.0 T with standard spin-echo sequences. Maps were fused into single color images, with R1, R2. and proton density coded, respectively, by red, green, and blue. The possibility of standardizing the technique was evaluated. Comparative analysis of color and conventional MR images of white matter disease and brain tumors was performed to assess intra- and interob-server variability. Unequivocal and reproducible chromatic characterization of normal brain structures and a variety of lesions was obtained. Intra-and interobserver analysis showed that color images can be used as a diagnostic tool. The technique may provide a simplified and timesaving approach for interpretation and presentation of brain MR studies.     

Theoretical-experimental comparison of contrast in magnetic resonance

The intensity of MR signal depends on several parameters, such as proton density [N(H)], relaxation times (T1 and T2), repetition time (TR), and echo time (TE). A theoretical model describes this dependence, which is currently employed for image optimization. It allows the evaluation of image contrast once the tissue parameters are known. The above-mentioned theoretical model was tested with the use of CuSO4 samples at various concentrations for which T1 and T2 values were known from the literature. Our unit was an ESATOM MR 5000 which employed a 0.5 Tesla magnetic field. We used spin-echo sequences with TR = 500, 1000 ms and TE ranging from 50 to 150 ms. Signal intensity was measured both by direct access to the data matrix and with the use of the pixel intensity calculation program for regions of interest. The difference in the signals corresponding to the various samples were determined to evaluate the contrast. Our results are in strict agreement with those from the theoretical model. The latter can thus be employed for image optimization.     

Analysis of multiple T2 proton relaxation processes in human head and imaging on the basis of selective and assigned T2 values

Two-dimensional T2-selective proton imaging of human head has been performed at 10 MHz employing Carr-Purcell-Gill-Meiboom pulse trains with echo separation of 6 and 12 msec. Using the information of 36 spin echoes and applying a recovery time of 2 sec the magnetization decay has been traced. The multi-exponential T2 relaxation curve of each pixel of the image has been deconvoluted into up to three mono-exponential functions which are defined by 1/T2,i and alpha i. The T2,i values are represented in a T2 histogram of the slice and then selected for generating images which relate to protons with specific T2. The alpha i values indicate the relative amounts of the T2-selected protons. Imaging on the basis of alpha i values increases the contrast of the image. The multiple T2-selective imaging technique leads to head images which show selectively cerebrospinal fluid, gray matter, different types of white matter, and nonassigned fast-relaxing proton classes.     

Uptake of a superparamagnetic contrast agent imaged by MR with high spectral and spatial resolution.

Conventional MRI implicitly treats the proton signal as a single, narrow Lorentzian. However, water signals in vivo are often in homogeneously broadened and have multiple resolvable components. These components represent discrete populations of water molecules within each pixel which are affected differently by physiology and contrast agents. Accurate measurement of each component of the water resonance can improve anatomic and functional MR images and provide insight into the structure and dynamics of subpixelar microenvironments. This report describes high spectral and spatial resolution (HiSS) MR imaging of rodent prostate tumors before and after injection of a superparamagnetic contrast agent. HiSS datasets were used to synthesize images in which intensity is proportional to peak height, peak frequency, and linewidth. These images showed anatomic features which were not clearly delineated in conventional T(2) and gradient echo images. HiSS images obtained after injection of the contrast agent showed T *(2) and T(1) changes which were not seen in conventional images. These changes are associated with microvessel density and permeability. The results suggest HiSS with superparamagnetic contrast agents has the potential to improve characterization of tumors.     

A simple method to improve image nonuniformity of brain MR images at the edges of a head coil

One of the major sources of image nonuniformity in the high field MR scanners is the radiofrequency (RF) coil inhomogeneity. It degrades conspicuity of lesion(s) in the MR images of the brain and surrounding tissues and reduces accuracy of image postprocessing particularly at the edges of the coil. In this investigation, we have devised and tested a simple method to correct for nonuniformity of MR images of the brain at the edges of the RF head coil. Initially, a cylindrical oil phantom, which fit exactly in the head coil, was scanned on a 1.5 T imager. Then, a correction algorithm identified a reference pixel value in the phantom at the most homogeneous region of the RF coil. Next, every pixel inside the phantom was normalized relative to this reference value. The resulting set of coefficients or "correction matrices" was obtained for different types of MR contrast agent. Finally, brain MR images of normal subjects and multiple sclerosis patients were acquired and processed by the corresponding correction matrices obtained with different pulse sequences. Application of correction matrices to brain MR images showed a gain in pixel intensity particularly in the slices at the edge of the coil.     

Noninvasive in vivo magnetic resonance imaging of injury-induced neointima formation in the carotid artery of the apolipoprotein-E null mouse

Mice deficient in apolipoprotein-E (apoE) experience severe hypercholesterolemia, are prone to atherosclerosis, and recently have emerged as a powerful tool in the study of plaque formation. In this study, we developed magnetic resonance (MR) imaging methods to detect the progression of atherosclerosis noninvasively in a mouse model of arterial injury. Four 14-week-old apoE-deficient mice were imaged 5 weeks after beginning an atherogenic Western diet and 4 weeks after wire denudation injury of the left common carotid artery (LCCA). Information from several images was combined into high-information content images using methods previously developed. The image resolution was 47 x 47 x 750 microm(3). We acquired T1-, T2-, and proton density (PD)-weighted images (TR/TE 650/14, 2000/60, and 2000/14 msec, respectively). Each 8-bit image was placed in a separate color channel to produce a 24-bit color image (red = T1, green = PD, and blue = T2). The composite image created contrast between different tissue types that was superior to that of any single image and revealed significant luminal narrowing of the LCCA, but not the uninjured RCCA. MR images were compared with corresponding histopathology cross sections and luminal area measurements from each method correlated(r2= 0.61). Atherosclerotic luminal narrowing was successfully detected through MR imaging in a mouse model of arterial injury that is small, reproduces quickly, and lends itself to genetic analysis and manipulation.     

Correlation of MR imaging and histologic findings in mouse melanoma.

M2R melanoma tumors in male C57 black mice were used to correlate magnetic resonance (MR) images with the corresponding histologic slices and to determine if analysis of the achievable correlation can provide a basis for predicting gross histologic features with MR imaging alone. The MR imaging sections obtained at 4.7 T were each 680 microns thick, with an in-plane resolution of 195 microns. The distribution of melanin within the histologic slices correlated well with the high-signal-intensity regions on the T1-weighted images (T1WIs), while these regions had low signal intensity on the T2-weighted images (T2WIs), providing evidence that melanin or melanin-associated paramagnetic species are responsible for the observed proton relaxation rate enhancement. Viable melanoma cells typically showed intermediate signal intensity on T2WIs, T1WIs, and proton-density images. Necrosis typically had high signal intensity on T2WIs, T1WIs, and proton-density images. Quantitation of the MR imaging results, followed by statistical analysis, demonstrated statistically significant differences between melanin-rich, viable-melanoma, and necrotic regions on MR images.     

Temperature mapping with MR imaging of molecular diffusion: application to hyperthermia

Efficacy and safety considerations for hyperthermia (HT) cancer therapy require accurate temperature measurements throughout the heated volume. Noninvasive thermometry methods have been proposed, including magnetic resonance (MR) imaging based on the temperature dependence of the relaxation time T1. However, the temperature accuracy achieved to date with T1 measurements does not fulfill the HT requirements (1 degree C/cm). The authors propose to use molecular diffusion, for which temperature dependence is well known. Molecular diffusion is more sensitive than T1 and can be determined with high accuracy with MR imaging. Diffusion and derived temperature images were obtained with a 2 X 2-mm pixel size in a polyacrylamide gel phantom heated inside the head coil of a clinical 0.5-T whole-body MR imaging system by means of a modified clinical HT device made compatible with the system. Temperatures determined from these images with 0.8-cm2 regions of interest were found to be within 0.5 degrees C of those recorded with thermocouples placed inside the gel. The utility of this method in clinical hyperthermia is enhanced by its potential to also help monitor blood perfusion.     

Intracranial meningiomas: correlation between MR imaging and histology in fifty patients

The magnetic resonance (MR) findings in 50 surgically verified intracranial meningiomas were reviewed. An attempt was made to correlate their signal intensity on spin echo (SE) T1-weighted, proton density [N(H)], and T2-weighted images with the different histologic subtypes. The T1-weighted images were nonspecific in differentiating the subtypes of meningiomas. On proton density and T2-weighted images, more information was available, but there remained large group (46%) of meningiomas that were not classifiable. The average signal intensity scores on T1-weighted, proton density, and T2-weighted images in the different histologic subtypes were correlated with each other using the Student t test. Only one significant correlation (psammomatous-anaplastic) and three almost significant correlations (syncytial-transitional or psammomatous and transitional-psammomatous) were found. Different histologic subtypes may have a different MR appearance, but this does not suffice to reach a histologic diagnosis by MR imaging.     

Human malignant melanomas with varying degrees of melanin content in nude mice: MR imaging, histopathology, and electron paramagnetic resonance

The etiology of the paramagnetic relaxation enhancement seen in malignant melanoma on proton magnetic resonance (MR) images has been the subject of many recent investigations and has been ascribed to iron from associated hemorrhage or chelated metal ions, rather than directly due to melanin. The purpose of this study was to correlate proton relaxation times on MR images in malignant melanomas with histopathologic features (i.e., degree of pigmentation, iron deposition, and necrosis), water content, and electron paramagnetic resonance (EPR) spectra to elucidate the etiology of the relaxation behavior demonstrated by these neoplasms. Cultured cells derived from human malignant melanoma metastases were implanted subcutaneously into nude mice. Twelve separate lesions were evaluated in 10 mice. Magnetic resonance imaging was performed in vivo at 1.9 T using spin echo and inversion recovery acquisitions for the purposes of calculating T1, T2, and proton density [N(H)]. Histopathologic examination was performed on specimens resected immediately after imaging, using hematoxylin/eosin, Prussian blue, and Fontana stains to assess tumor necrosis, and iron and melanin content. Dry/wet weight ratios and EPR spectra were also obtained on resected specimens. Our results indicate that T1 shortening correlates with increasing melanin content and not with increasing iron deposition, EPR-active metallic cations, necrosis, or water content. In fact, a presumably unrelated statistical correlation was found between increased iron and T1 prolongation. The T2 relaxation times did not correlate with the presence of any single factor other than proton density. Although the unique relaxation behavior of nonhemorrhagic malignant melanoma in vivo cannot be traced to a single cause, our data suggest that, contrary to previous investigations, it is strongly influenced by the presence of melanin rather than iron or other naturally occurring paramagnetic ions.     

Magnetic resonance imaging and mathematical modeling of progressive formalin fixation of the human brain.

The temporal magnetic resonance (MR) appearance of human brain tissue during formalin fixation was measured and modeled using a diffusion mathematical model of formalin fixation. Coronal MR images of three human brains before formalin fixation and at multiple time points thereafter were acquired. T1 relaxation, T2 relaxation, water apparent diffusion coefficient (ADC), and proton density (PD) maps were calculated. The size of a light "formalin band" region, visible in T1 weighted images, was compared to a mathematical model of diffusive mass transfer of formalin into the brain. T1 relaxation, T2 relaxation, and PD all decreased, in both gray and white matter, as formalin fixation progressed. The ADC remained more or less constant. The location of the inner boundary of the formalin band followed a time course consistent with the steepest formalin concentration gradient in the mathematical model. Based on the diffusion model, the brain is not completely saturated in formalin until after 14.8 weeks of formalin immersion and, based on the observed changes in T1, T2, and PD, fixation is not complete until after 5.4 weeks. During fixation, the ongoing attenuation of T1 relaxation, T2 relaxation, and PD must be taken into consideration when performing postmortem MRI studies.     

Frequency encoding for simultaneous display of multimodality images.

An original method for simultaneous display of functional and anatomic images, based on frequency encoding (FE), merges color PET with T1-weighted MR brain images, and grayscale PET with multispectral color MR images. A comparison with two other methods reported in the literature for image fusion (averaging and intensity modulation techniques) was performed. METHODS: For FE, the Fourier transform of the merged image was obtained summing the low frequencies of the PET image and the high frequencies of the MR image. For image averaging, the merged image was obtained as a weighted average of the intensities of the two images to be merged. For intensity modulation, the red, green and blue components of the color image were multiplied on a pixel-by-pixel basis by the grayscale image. A comparison of the performances of the three techniques was made by three independent observers assessing the conspicuity of specific MRI and PET information in the merged images. For evaluation purposes, images from seven patients and a computer-simulated MRI/PET phantom were used. Data were compared with a chi-square test applied to ranks. RESULTS: For the depiction of MRI and PET information when merging color PET and T1-weighted MR images, FE was rated superior to intensity modulation and averaging techniques in a significant number of comparisons. For merging grayscale PET with multispectral color MR images, FE and intensity modulation were rated superior to image averaging in terms of both MRI and PET information. CONCLUSION: The data suggest that improved simultaneous evaluation of MRI and PET information can be achieved with a method based on FE.     

RLSQ: T1, T2, and rho calculations, combining ratios and least squares

A method is proposed for calculating images of spin density rho and relaxation times T1, T2 from MR imaging sequences with the following characteristics: it avoids nonlinear iterative techniques and is consequently sufficiently fast to allow calculation of rho, T1, and T2 on a pixel by pixel basis routinely; it uses all available data simultaneously in an approximately optimal fashion. The method is applied to the combined SE + IR sequence but can be used for other acquisition schemes as well.