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 I) 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 normal 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 approximately 0.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."     

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.     

Characteristics of geometric distortion correction with increasing field-of-view in open-configuration MRI

Open-configuration magnetic resonance imaging (MRI) systems are becoming increasingly desirable for musculoskeletal imaging and image-guided radiotherapy because of their non-claustrophobic configuration. However, geometric image distortion in large fields-of-view (FOV) due to field inhomogeneity and gradient nonlinearity hinders the practical applications of open-type MRI. We demonstrated the use of geometric distortion correction for increasing FOV in open MRI. Geometric distortion was modeled and corrected as a global polynomial function. The appropriate polynomial order was identified as the minimum difference between the coordinates of control points in the distorted MR image space and those predicted by polynomial modeling. The sixth order polynomial function was found to give the optimal value for geometric distortion correction. The area of maximum distortion was < 1 pixel with an FOV of 285 mm. The correction performance error was increased at most 1.2% and 2.9% for FOVs of 340 mm and ~ 400 mm compared with the FOV of 285 mm. In particular, unresolved distortion was generated by local deformation near the gradient coil center.     

Correction of spatial distortion in magnetic resonance angiography for radiosurgical treatment planning of cerebral arteriovenous malformations.

A treatment planning system based on magnetic resonance (MR) angiographic imaging data for the radiosurgery of inoperable cerebral arteriovenous malformations is reported. MR angiography was performed using a three-dimensional (3D) velocity-compensated fast imaging with steady-state precession (FISP) sequence. Depending on the individual MR system, inhomogeneities and nonlinearities induced by eddy currents during the pulse sequence can distort the images and produce spurious displacements of the stereotactic coordinates in both the x-y plane and the z axis. If necessary, these errors in position can be assessed by means of two phantoms placed within the stereotactic guidance system--a "2D-phantom" displaying "pincushion" distortion in the image, and a "3D-phantom" displaying displacement, warp, and tilt of the image plane itself. The pincushion distortion can be "corrected" (reducing displacements from 2-3 mm to 1 mm) by calculations based on modeling the distortion as a fourth order 2D polynomial. Displacement, warp, and tilt of the image plane may be corrected by adjustment of the gradient shimming currents. After correction, the accuracy of the geometric information is limited only by the pixel resolution of the image (= 1 mm). Precise definition of the target volume could be performed by the therapist either directly in the MR images or in calculated projection MR angiograms obtained by a maximum intensity projection algorithm. MR angiography provides a sensitive, noninvasive 3D method for defining target volume and critical structures, and for calculating precise dose distributions for radiosurgery of cerebral arteriovenous malformations.     

Field-map correction in read-out segmented echo planar imaging for reduced spatial distortion in prostate DWI for MRI-guided radiotherapy applications

Diffusion-weighted echo planar imaging (DW-EPI) suffers from geometric distortion due to low phase-encoding bandwidth. Read-out segmented echo planar imaging (RS-EPI) reduces distortion but residual distortion remains in extreme cases. Additional corrections need to be applied, especially for radiotherapy applications where a high degree of accuracy is needed. In this study the use of magnetic field map corrections are assessed in DW-EPI and RS-EPI, to reduce geometric uncertainty for MRI-guided radiotherapy applications. Magnetic field maps were calculated from gradient echo images and distortion corrections were applied to RS-EPI images. Distortions were assessed in a prostate phantom by comparing to the known geometry, and in vivo using a modified Hausdorff distance metric using a T2-weighted spin echo as ground truth. Across 10 patients, field map-corrected RS-EPI reduced maximum distortion by 5 mm on average compared to DW-EPI (σ = 1.9 mm). Geometric distortions were also reduced significantly using field mapping with RS-EPI, compared to RS-EPI alone (p ≤ 0.05). The increased geometric accuracy of these techniques can potentially allow diffusion-weighted images to be fused with other MR or CT images for radiotherapy treatment purposes.     

Analysis of machine-dependent and object-induced geometric distortion in 2DFT MR imaging.

Geometric distortion in MR imaging predominantly arises from the inhomogeneity of the static field and the nonlinearity of the gradients. It is the purpose of this paper to analyse the object and machine related contributions to geometric distortion in order to determine which corrections are necessary for attaining a specified precision. System related imperfections were measured by systematic variation of the strength, direction, and polarity of the read-out gradient in imaging experiments on a grid of cylindrical sample tubes. For the 1.5-T system used in this study, static field related errors up to 7 mm and gradient related errors up to 4 mm were observed (midcoronal plane, FOV 400-mm, G-read between 0.5 and 3.0 mT/m). Field related errors were shown to be inversely proportional to gradient strength, whereas gradient related errors turned out to be virtually independent of gradient strength. It therefore seems recommendable to always apply the strongest available selection and read-out gradients when geometric fidelity is given preference to signal-to-noise considerations. Correction of system related geometric distortions in MR images can readily be performed by table lookup. Object-induced distortions of the gradient fields were studied by experiments on a grid of sample tubes immersed into a cylindrical water bath of variable saline concentration. These experiments revealed a negligible influence of the object on the gradient error distribution, and lead to the conclusion that correction for the nonlinearity of the gradients only requires the application of system dependent correction factors. Object-related distortions of B0 were studied by conventional SE and fat-suppressed IR experiments on phantoms and human subjects. In these experiments the polarity of the read-out gradient was reversed. Subtraction images showed significant object-induced inhomogeneities of the static field at tissue-air interfaces and in the immediate vicinity of the object being imaged. A first attempt to correct for object related B0 inhomogeneities was made by contour analysis of the source images. At present this correction still has to be done manually.     

A hybrid strategy for correcting geometric distortion in echo-planar images

A hybrid strategy for geometric distortion correction of echo-planar images is demonstrated. This procedure utilizes standard field mapping for signal displacement correction and the so-called reverse gradient acquisition for signal intensity correction. (The term reverse gradient refers to an acquisition of two sets of echo-planar images with phase encoding gradients of opposite polarity.) The hybrid strategy is applied to human brain echo-planar images acquired with and without diffusion-weighting. A comparison of the hybrid distortion corrected images to those corrected with standard field mapping only demonstrates much better performance of the hybrid method. A variant of the hybrid method is also demonstrated which requires the acquisition of only one pair of opposite polarity images within a set of images.     

Field camera versus phantom-based measurement of the gradient system transfer function (GSTF) with dwell time compensation

PurposeThe gradient system transfer function (GSTF) can be used to describe the dynamic gradient system and applied for trajectory correction in non-Cartesian MRI. This study compares the field camera and the phantom-based methods to measure the GSTF and implements a compensation for the difference in measurement dwell time.MethodsThe self-term GSTFs of a MR system were determined with two approaches: 1) using a dynamic field camera and 2) using a spherical phantom-based measurement with standard MR hardware. The phantom-based GSTF was convolved with a box function to compensate for the dwell time dependence of the measurement. The field camera and phantom-based GSTFs were used for trajectory prediction during retrospective image reconstruction of 3D wave-CAIPI phantom images.ResultsDifferences in the GSTF magnitude response were observed between the two measurement methods. For the wave-CAIPI sequence, this led to deviations in the GSTF predicted trajectories of 4% compared to measured trajectories, and residual distortions in the reconstructed phantom images generated with the phantom-based GSTF. Following dwell-time compensation, deviations in the GSTF magnitudes, GSTF-predicted trajectories, and resulting image artifacts were eliminated (< 0.5% deviation in trajectories).ConclusionWith dwell time compensation, both the field camera and the phantom-based GSTF self-terms show negligible deviations and lead to strong artifact reduction when they are used for trajectory correction in image reconstruction.     

Artifacts caused by transcranial magnetic stimulation coils and EEG electrodes in T(2)*-weighted echo-planar imaging

We investigated the effects of transcranial magnetic stimulation (TMS) coils and electroencephalographic (EEG) electrodes on T(2)*-weighted echo-planar images (EPI) at 2.0 T (gradient-echo EPI, mean TE = 53 ms, 2x2x4 mm(3)). In comparison with anatomic gradient-echo images (3D FLASH, TE = 4 ms, 1x1x1 mm(3)), T(2)*-weighted EPI acquisitions of a water-filled spherical phantom revealed severe signal losses and geometric distortions in the vicinity of TMS coils. Even remote effects were observed for image orientations perpendicular to the coil plane. EEG electrodes and the fixation gel caused milder localized distortions. In humans, complications were avoided by the large distance between the TMS coil and the cortical surface and when using an EPI orientation parallel to the plane of the coil. It is concluded that T(2)*-weighted EPI studies of human brain function may be performed without distortions caused by TMS coils and EEG electrodes.     

Computer-generated fMRI phantoms with motion-distortion interaction

As an extension of previous work on computer-generated phantoms, more accurate, realistic phantoms are generated by integrating image distortion and signal loss caused by susceptibility variations. With the addition of real motions and activations determined from actual functional MRI studies, these phantoms can be used by the fMRI community to assess with higher fidelity pre-processing algorithms such as motion correction, distortion correction and signal-loss compensation. These phantoms were validated by comparison to real echo-planar images. Specifically, studies have shown the effects of motion–distortion interactions on fMRI. We performed motion correction and activation analysis on these phantoms based on a block paradigm design using SPM2, and the results demonstrate that interactions between motion and distortion affect both motion correction and activation detection and thus represent a critical component of phantom generation.     

Mapping eddy current induced fields for the correction of diffusion-weighted echo planar images

Diffusion-weighted echo-planar magnetic resonance imaging is potentially of great importance as a diagnostic imaging tool; however, the technique currently suffers a number of limitations, including the image distortion caused by the eddy current induced fields when the diffusion-weighting magnetic field gradient pulses are applied. The distortions cause mis-registration between images with different diffusion-weighting, that then results in artifacts in quantitative diffusion images. A method is presented to measure the magnetic fields generated from the eddy currents for each of three orthogonal gradient pulse vectors, and then to use these to ascertain the image distortion that occurs in subsequent diffusion-weighted images with arbitrary gradient pulse vector amplitude and direction, and image plane orientation. The image distortion can then be reversed. Both temporal and spatial dependence of the residual eddy current induced fields are included in the analysis. Image distortion was substantially reduced by the correction scheme, for arbitrary slice position and angulation. This method of correction is unaffected by the changes in image contrast that occur due to diffusion weighting, and does not need any additional scanning time during the patient scan. It is particularly suitable for use with single-shot echo planar imaging.     

MRI phase offset correction method impacts quantitative susceptibility mapping

Individual channel ultra-high field (7T) phase images have to be phase offset corrected prior to the mapping of magnetic susceptibility of tissue. Whilst numerous methods have been proposed for gradient recalled echo MRI phase offset correction, it remains unclear how they affect quantitative magnetic susceptibility values derived from phase images. Methods already proposed either employ a single or multiple echo time MRI data. In terms of the latter, offsets can be derived using an ultra-short echo time acquisition, or by estimating the offset based on two echo points with the assumption of linear phase evolution with echo time. Our evaluation involved 32 channel multi-echo time 7T GRE (Gradient Recalled Echo) and ultra-short echo time PETRA (Pointwise Encoding Time Reduction with Radial Acquisition) MRI data collected for a susceptibility phantom and three human brains. The combined phase images generated using four established offset correction methods (two single and two multiple echo time) were analysed, followed by an assessment of quantitative susceptibility values obtained for a phantom and human brains. The effectiveness of each method in removing the offsets was shown to reduce with increased echo time, decreased signal intensity and reduced overlap in coil sensitivity profiles. Quantitative susceptibility values and how they change with echo time were found to be method specific. Phase offset correction methods based on single echo time data have a tendency to produce more accurate and less noisy quantitative susceptibility maps in comparison with methods employing multiple echo time data.     

Correction of B0 susceptibility induced distortion in diffusion-weighted images using large-deformation diffeomorphic metric mapping

Geometric distortion caused by B0 inhomogeneity is one of the most important problems for diffusion-weighted images (DWI) using single-shot, echo planar imaging (SS-EPI). In this study, large-deformation, diffeomorphic metric mapping (LDDMM) algorithm has been tested for the correction of geometric distortion in diffusion tensor images (DTI). Based on data from nine normal subjects, the amount of distortion caused by B0 susceptibility in the 3-T magnet was characterized. The distortion quality was validated by manually placing landmarks in the target and DTI images before and after distortion correction. The distortion was found to be up to 15 mm in the population-averaged map and could be more than 20 mm in individual images. Both qualitative demonstration and quantitative statistical results suggest that the highly elastic geometric distortion caused by spatial inhomogeneity of the B0 field in DTI using SS-EPI can be effectively corrected by LDDMM. This postprocessing method is especially useful for correcting existent DTI data without phase maps.     

Experimental and computational analyses of the effects of slice distortion from a metallic sphere in an MRI phantom

Susceptibility artifacts due to metallic prostheses are a major problem in clinical magnetic resonance imaging. We theoretically and experimentally analyze slice distortion arising from susceptibility differences in a phantom consisting of a stainless steel ball bearing embedded in agarose gel. To relate the observed image artifacts to slice distortion, we simulate images produced by 2D and 3D spin-echo (SE) and a view angle tilting (VAT) sequence. Two-dimensional SE sequences suffer from extreme slice distortion when a metal prosthesis is present, unlike 3D SE sequences for which--since slices are phase-encoded--distortion of the slice profile is minimized, provided the selected slab is larger than the region of interest. In a VAT sequence, artifacts are reduced by the application of a gradient along the slice direction during readout. However, VAT does not correct for the excitation slice profile, which results in the excitation of spins outside the desired slice location and can lead to incorrect anatomical information in MR images. We propose that the best sequences for imaging in the presence of a metal prosthesis utilize 3D acquisition, with phase encoding replacing slice selection to minimize slice distortion, combined with excitation and readout gradient strengths at their maximum values.     

CSRm闭轨畸变及其校正的模拟研究

在给定的磁铁安装误差和磁场加工误差的条件下,对兰州重离子加速器冷却储存环主环的闭轨畸变及其校正进行了计算机模拟研究,在典型的误差分布下,校正前的水平方向及垂直方向的最大闭轨畸变分别为3.08mm和2.73mm,校正模拟的结果显示CSRm的闭轨畸变可以控制在足够小的范围内.     

In vivo isotropic 3D diffusion tensor mapping of the rat brain using diffusion-weighted 3D MP-RAGE MRI

The purpose of this study was to examine the potential of diffusion-weighted (DW) three-dimensional (3D) MP-RAGE MRI for diffusion-tensor mapping of the rat brain in vivo. A DW-3D-MP-RAGE (3D-DWI) sequence was implemented at 2.0 T using six gradient orientations and a b value of 1000 s/mm2. In this sequence, the preparation sequence with a "90 degrees RF-motion proving gradient (MPG): MPG-180 degrees RF-MPG-90 degrees RF" pulse train (DW driven equilibrium Fourier transform) was used to sensitize the magnetization to diffusion. A centric k-space acquisition order was necessary to minimize saturation effects (T1 contamination) from tissues with short relaxation time. The image matrix was 128x128x128 (interpolated from 64x64x64 acquisitions), which resulted in small isotropic DW image data (voxel size: 0.273x0.273x0.273 mm3). Moreover, 3D-DWI-derived maps of the fractional anisotropy (FA), relative anisotropy (RA) and main-diffusion direction were completely free of susceptibility-induced signal losses and geometric distortions. Two well-known commissural fibers, the corpus callosum and anterior commissure, were indicated and shown to be in agreement with the locations of these known stereotaxic atlases. The experiment took 45 min, and shorter times should be possible in clinical application. The 3D-DWI sequence allows for in vivo 3D diffusion-tensor mapping of the rat brain without motion artifacts and susceptibility to distortion.     

Deconvolution of vibroacoustic images using a simulation model based on a three dimensional point spread function

Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7–0.8 mm. Current VA systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric mean filters. To assess the filters’ performance on the restored images, we use an image quality index that accounts for correlation loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener filter is 0.9 (when the index is 1.0 this indicates perfect restoration). This filter yielded the best result in comparison with the other ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted in the difference-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient information to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information to spot the wires in the images. It is demonstrated that the phase of the theoretical and the experimental PSFs are dissimilar. This fact prevents VA image restoration with the current theoretical PSF. This study is a preliminary step towards understanding the restoration of VA images through the application of deconvolution filters.     

Thermodynamics of black plane solution

We obtain a new phantom black plane solution in $4$ 4 D of the Einstein–Maxwell theory coupled with a cosmological constant. We analyse their basic properties, as well as its causal structure, and obtain the extensive and intensive thermodynamic variables, as well as the specific heat and the first law. Through the specific heat and the so-called geometric methods, we analyse in detail their thermodynamic properties, the extreme and phase transition limits, as well as the local and global stabilities of the system. The normal case is shown with an extreme limit and the phantom one with a phase transition only for null mass, which is physically inaccessible. The systems present local and global stabilities for certain values of the entropy density with respect to the electric charge, for the canonical and grand canonical ensembles.     

平面物体在曲面状态下扫描仪图像的校正理论

平面物体在曲面状态下经扫描仪扫描后,其图像将发生复杂的畸变。提出将其分类为灰度畸变、投影畸变和成像畸变。通过理论分析,提出了在二元曲面模型下对投影畸变和成像畸变进行数字校正的方法,给出了对灰度畸变进行数字校正的实用方法。     

Automatic correction of in-plane bulk motion artifacts in self-navigated radial MRI

Radial MRI is typically used for scans that are sensitive to unavoidable motion. While the translational motion artifact can be easily removed from the radial trajectory data by phase correction, correction of rotational motion still remains a challenge in radial MRI. We present a novel method to refocus the image corrupted by view-to-view motion in the view-interleaved radial MRI data. In this method, the error in rotational view angles was modeled as a polynomial function of the view order and the model parameters were estimated by minimizing the self-navigator image metrics such as image entropy, gradient entropy, normalized gradient squared and mean square difference. Translational motion correction was conducted by aligning the projection profiles. Simulation studies were conducted to demonstrate the robustness of both translational and rotational motion correction methods in different noise levels. The proposed method was successfully applied to correct for motion of healthy subjects. Substantial motion correction with relative error of less than 5% was achieved by using either first- or second-order model with the image metrics. This study demonstrates the potential of the method for motion-sensitive applications.