PubMed Health⌕ Search

PubMed · 10873196

An efficient geometric image distortion correction method for a biplanar planar gradient coil.

Abstract

Since the spatial field non-linearity of gradient coils translates into image geometric distortion in MRI, in many applications, such as cardiac function analysis and interventional MR-based device tracking/guidance, where the precise geometric information is needed, the presence of geometric image distortion can not be simply ignored. To address the concern for geometric image distortion, we have developed and validated a general and efficient numerical technique for parameterizing the global image distortion for a bi-planar gradient coil as well as accomplishing image restoration as a post-imaging processing. This image correction methodology is based on a global distortion coordinate mapping function which can be systematically defined directly from the gradient field non-linearity in 3-dimension (3D) of a given gradient coil. The image correction was carried out in two steps: (1) map each pixel of the corrected image representation onto its distorted image according to the distortion mapping; (2) interpolate the pixel intensity in the distorted image using its neighboring points via a bi-linear interpolation procedure. The results showed clearly that the distortion correction method was robust in term of the capability of reducing image geometric distortion dramatically. Also it is shown that the magnetic field non-linearity or the image distortion of a typical bi-planar gradient coil can be adequately parameterized using a finite Taylor series expansion based on its design parameters. Furthermore, this image distortion correction method is very efficient in practice for performing 3D correction for any image orientation since a compact parameterized field expression contains non-zero terms.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

H Liu. 2000. An efficient geometric image distortion correction method for a biplanar planar gradient coil.. https://doi.org/10.1007/bf02601841

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Searching for the best model: ambiguity of inverse solutions and application to fetal magnetoencephalography.

Fetal brain signals produce weak magnetic fields at the maternal abdominal surface. In the presence of much stronger interference these weak fetal fields are often nearly indistinguishable from noise. Our initial objective was to validate these weak fetal brain fields by demonstrating that they agree with the electromagnetic model of the fetal brain. The fetal brain model is often not known and we have attempted to fit the data to not only the brain source position, orientation and magnitude, but also to the brain model position. Simulation tests of this extended model search on fetal MEG recordings using dipole fit and beamformers revealed a region of ambiguity. The region of ambiguity consists of a family of models which are not distinguishable in the presence of noise, and which exhibit large and comparable SNR when beamformers are used. Unlike the uncertainty of a dipole fit with known model plus noise, this extended ambiguity region yields nearly identical forward solutions, and is only weakly dependent on noise. The ambiguity region is located in a plane defined by the source position, orientation, and the true model centre, and will have a diameter approximately 0.67 of the modelled fetal head diameter. Existence of the ambiguity region allows us to only state that the fetal brain fields do not contradict the electromagnetic model; we can associate them with a family of models belonging to the ambiguity region, but not with any specific model. In addition to providing a level of confidence in the fetal brain signals, the ambiguity region knowledge in combination with beamformers allows detection of undistorted temporal waveforms with improved signal-to-noise ratio, even though the source position cannot be uniquely determined.

Biophysical Phenomena↗

Recovery of total I-131 activity within focal volumes using SPECT and 3D OSEM.

We experimentally investigated the SPECT recovery of I-131 activity in multiple spheres located simultaneously at different locations within a cylindrical phantom that had an elliptical cross section. The sphere volumes ranged from 209 cc down to 4.2 cc. A Prism 3000 camera and two types of parallel-hexagonal-hole collimation were employed: high energy (HE) and ultra high energy (UHE). Using appropriately-different 3D models of the point source response function for the two types of collimation, approximately the same recovery of activity could be achieved with either collimation by 3D OSEM reconstruction. The recovery coefficient was greater with no background activity in the phantom by 0.10, on average, compared to that with background. In the HE collimation case, the activity recovery was considerably better for all volumes using 3D OSEM reconstruction than it had been in the past using 1D SAGE reconstruction. Recovery-coefficient-based correction in a simulated patient case involving spherical tumours moderately improved the activity estimates (average error reduced from 14% to 9% for UHE collimation, and from 15% to 11% for HE collimation). For a test case with HE collimation, increasing the projection-image sampling density while decreasing the image voxel size increased the recovery coefficient by 0.075 on average, and, if used in a full set of calibration measurements of recovery coefficient versus volume, might lead to further improvement in accuracy for the patient case.

Biophysical Phenomena↗

Improvement of spatial resolution in the longitudinal direction for isotropic imaging in helical CT.

Experiments were conducted to confirm the isotropic spatial resolution of multislice CT with a 0.5 mm slice thickness. Isotropic spatial resolution means that the spatial resolution in the transaxial plane (X-Y plane) and that in the longitudinal direction (Z direction) are equivalent. To obtain point spread function (PSF) values in the X-Y-Z directions, three-dimensional voxel data were obtained by helical scanning of a bead phantom. The modulation transfer function (MTF) values were then obtained by three-dimensional Fourier transform of the PSF. Evaluation of the spatial resolution in the X-Y-Z directions by the MTF values showed that the spatial resolution in the Z direction does not depend on the reconstruction kernel used. It was also found that the spatial resolution in the Z direction, as compared with that in the X-Y plane, is superior with the standard kernel for the abdomen and is inferior with the high-definition kernel for the ears/bones. By performing sharpening filter processing in the Z direction with a high-definition kernel, comparable spatial resolution could be obtained in the X-Y-Z directions. It was confirmed that adjusting the spatial resolution in the Z direction with the reconstruction kernel used is an effective method for isotropic imaging.

Biophysical Phenomena↗