A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Recently, a 3-dimensional phantom that can provide a comprehensive, accurate and complete measurement of the geometric distortion in MRI has been developed. In this paper, a scheme for characterizing the measured geometric distortion using the 3-D phantom is described. In the proposed scheme, a number of quantitative measures are developed and used to characterize the geometric distortion. These measures encompass the overall and spatial aspects of the geometric distortion. Two specific types of volume of interest, rectangular parallelepipeds (including cubes) and spheres are considered in the proposed scheme. As an illustration, characterization of the geometric distortion in a Siemens 1.5T Sonata MRI system using the proposed scheme is presented. As shown, the proposed scheme provides a comprehensive assessment of the geometric distortion. The scheme can be potentially used as a standard procedure for the assessment of geometric distortion in MRI. (C) 2004 American Association of Physicists in Medicine.

Geometric distortion in clinical MRI systems - Part II: correction using a 3D phantom

In this paper, we present the correction of the geometric distortion measured in the clinical magnetic resonance imaging (MRI) systems reported in the preceding paper (Part 1) using a 3D method based on the phantom-mapped geometric distortion data. This method allows the correction to be made on phantom images acquired without or with the vendor correction applied. With the vendor's 2D correction applied, the method corrects for both the residual geometric distortion still present in the plane in which the correction method was applied (the axial plane) and the uncorrected geometric distortion along the axis non-nal to the plane. The evaluation of the effectiveness of the correction using this new method was carried out through analyzing the residual geometric distortion in the corrected phantom images. The results show that the new method can restore the distorted images in 3D nearly to perfection. For all the MRI systems investigated, the mean absolute deviations in the positions of the control points (along x-, y- and z-axes) measured on the corrected phantom images were all less than 0.2 mm. The maximum absolute deviations were all below similar to0.8 mm. As expected, the correction of the phantom images acquired with the vendor's correction applied in the axial plane performed equally well. Both the geometric distortion still present in the axial plane after applying the vendor's correction and the uncorrected distortion along the z-axis have all been restored. (C) 2004 Elsevier Inc. All rights reserved.

Geometric distortion in clinical MRI systems - Part I: evaluation using a 3D phantom

Recently, a 3D phantom that can provide a comprehensive and accurate measurement of the geometric distortion in MRI has been developed. Using this phantom, a full assessment of the geometric distortion in a number of clinical MRI systems (GE and Siemens) has been carried out and detailed results are presented in this paper. As expected, the main source of geometric distortion in modern superconducting MRI systems arises from the gradient field nonlinearity. Significantly large distortions with maximum absolute geometric errors ranged between 10 and 25 mm within a volume of 240 x 240 x 240 mm(3) were observed when imaging with the new generation of gradient systems that employs shorter coils. By comparison, the geometric distortion was much less in the older-generation gradient systems. With the vendor's correction method, the geometric distortion measured was significantly reduced but only within the plane in which these 2D correction methods were applied. Distortion along the axis normal to the plane was, as expected, virtually unchanged. Two-dimensional correction methods are a convenient approach and in principle they are the only methods that can be applied to correct geometric distortion in a single slice or in multiple noncontiguous slices. However, these methods only provide an incomplete solution to the problem and their value can be significantly reduced if the distortion along the normal of the correction plane is not small. (C) 2004 Elsevier Inc. All rights reserved.

Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

A pain neuromatrix approach to patients with chronic pain

This paper presents an approach to rehabilitation of pain patients. The fundamental principles of the approach are (i) pain is an output of the brain that is produced whenever the brain concludes that body tissue is in danger and action is required, and (ii) pain is a multisystem output that is produced when an individual-specific cortical pain neuromatrix is activated. When pain becomes chronic, the efficacy of the pain neuromatrix is strengthened via nociceptive and non-nociceptive mechanisms, which means that less input, both nociceptive and non-nociceptive, is required to produce pain. The clinical approach focuses on decreasing all inputs that imply that body tissue is in danger and then on activating components of the pain neuromatrix without activating its output. Rehabilitation progresses to increase exposure to threatening input across sensory and non-sensory domains. (C) 2003 Elsevier Ltd. All rights reserved.

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Noninvasive localization of petroleum-derived spray oil in plants with chemical shift selective magnetic resonance imaging

Nuclear magnetic resonance (NMR) spectroscopy and magnetic resonance imaging (MRI) were used to detect petroleum-derived spray oils (PDSOs) in citrus seedlings and trees. The NMR spectrum of the phantom containing 10% (v/v) of a nC24 agricultural mineral oil (AMO) showed the resonance of the water protons at delta = 5 ppm, while the resonance of the oil protons at delta = 1.3 to 1.7 ppm. The peak resolution and the chemical shift difference of more than 3.3 ppm between water and oil protons effectively differentiated water and the oil. Chemical shift selective imaging (CSSI) was performed to localize the AMO within the stems of Citrus trifoliata L. seedlings after the application of a 4% (v/v) spray. The chemical shift selective images of the oil were acquired by excitation at delta = 1.5 ppm by averaging over 400 transients in each phase-encoding step. Oil was mainly detected in the outer cortex of stems within 10 d of spray application; some oil was also observed in the inner vascular bundle and pith of the stems at this point. CSSI was also applied to investigate the persistence of oil deposits in sprayed mature Washington navel orange (Citrus x aurantium L.) trees in an orchard. The trees were treated with either fourteen 0.25%, fourteen 0.5%, four 1.75%, or single 7% sprays of a nC23 horticultural mineral oil (HMO) 12 to 16 months before examination of plant tissues by CSSI, and were still showing symptoms of chronic phytotoxicity largely manifested as reduced yield. The oil deposits were detected in stems of sprayed flushes and unsprayed flushes produced 4 to 5 months after the last spray was applied, suggesting a potential movement of the oil via phloem and a correlation of the persistence of oil deposit in plants and the phytotoxicity. The results demonstrate that MRI is an effective method to probe the uptake and localization of PDSOs and other xenobiotics in vivo in plants noninvasively and nondestructively.

Method for a detailed measurement of image intensity nonuniformity in magnetic resonance imaging

In magnetic resonance imaging (MRI), the MR signal intensity can vary spatially and this spatial variation is usually referred to as MR intensity nonuniformity. Although the main source of intensity nonuniformity arises from B, inhomogeneity of the coil acting as a receiver and/or transmitter, geometric distortion also alters the MR signal intensity. It is useful on some occasions to have these two different sources be separately measured and analyzed. In this paper, we present a practical method for a detailed measurement of the MR intensity nonuniformity. This method is based on the same three-dimensional geometric phantom that was recently developed for a complete measurement of the geometric distortion in MR systems. In this paper, the contribution to the intensity nonuniformity from the geometric distortion can be estimated and thus, it provides a mechanism for estimation of the intensity nonuniformity that reflects solely the spatial characteristics arising from B-1. Additionally, a comprehensive scheme for characterization of the intensity nonuniformity based on the new measurement method is proposed. To demonstrate the method, the intensity nonuniformity in a 1.5 T Sonata MR system was measured and is used to illustrate the main features of the method. (c) 2005 American Association of Physicists in Medicine.