Susceptibility-matched envelope for the correction of EPI artifacts

Fast gradient echo sequences, such as echo planer imaging (EPI) and spiral imaging, are vulnerable to artifacts resulting from B(0) inhomogeneities. A major contribution to these artifacts is the susceptibility variation across the head, which is most severe in regions adjacent to air-tissue interfaces, such as the mouth, nasal sinuses, ears and the cortex. Susceptibility artifacts can cause geometrical distortions in the image as well as loss of signal due to T(2)* dephasing. The extent of these artifacts increases with the main field, thus compromising the signal-to-noise ratio (SNR) benefit gained in higher fields. In the current work, inhomogeneity caused by susceptibility variations at the external boundary of the human body has been corrected by surrounding the organs with a liquid without hydrogen atoms and whose susceptibility is similar to that of the imaged organ. EPI experiments were conducted on head-sized phantom, human brain, hand and legs. This method causes minimal patient inconvenience and no interference with any function of the scanner, thus yielding a simple and efficient solution for the correction of B(0) variation.     

Concurrent correction of geometric distortion and motion using the map-slice-to-volume method in echo-planar imaging

The accuracy of measuring voxel intensity changes between stimulus and rest images in fMRI echo-planar imaging (EPI) data is severely degraded in the presence of head motion. In addition, EPI is sensitive to susceptibility-induced geometric distortions. Head motion causes image shifts and associated field map changes that induce different geometric distortion at different time points. Conventionally, geometric distortion is "corrected" with a static field map independently of image registration. That approach ignores all field map changes induced by head motion. This work evaluates the improved motion correction capability of mapping slice to volume with concurrent iterative field corrected reconstruction using updated field maps derived from an initial static field map that has been spatially transformed and resampled. It accounts for motion-induced field map changes for translational and in-plane rotation motion. The results from simulated EPI time series data, in which motion, image intensity and activation ground truths are available, show improved accuracy in image registration, field corrected image reconstruction and activation detection.     

Referenceless one-dimensional Nyquist ghost correction in multicoil single-shot spatiotemporally encoded MRI

Single-shot spatiotemporally encoded (SPEN) MRI is a novel fast imaging method capable of retaining the time efficiency of single-shot echo planar imaging (EPI) but with distortion artifacts significantly reduced. Akin to EPI, the phase inconsistencies between mismatched even and odd echoes also result in the so-called Nyquist ghosts. However, the characteristic of the SPEN signals provides the possibility of obtaining ghost-free images directly from even and odd echoes respectively, without acquiring additional reference scans. In this paper, a theoretical analysis of the Nyquist ghosts manifested in single-shot SPEN MRI is presented, a one-dimensional correction scheme is put forward capable of maintaining definition of image features without blurring when the phase inconsistency along SPEN encoding direction is negligible, and a technique is introduced for convenient and robust correction of data from multi-channel receiver coils. The effectiveness of the proposed processing pipeline is validated by a series of experiments conducted on simulation data, in vivo rats and healthy human brains. The robustness of the method is further verified by implementing distortion correction on ghost corrected data.     

小动物高分辨扩散加权成像在临床 MRI 上的实现

在临床用MRI系统上对小动物扩散加权成像一般采用回波平面成像序列,但是回波平面成像易受偏共振效应的影响,得到的图像伪影大、几何变形严重、图像分辨率低,无法探究微小的生物组织结构. 该文报道了在临床用3 T MRI系统上采用自旋回波序列实现了高分辨扩散加权成像. 为减少运动伪影,序列中整合了导航回波矫正技术. 对脑缺血模型大鼠脑部的扫描结果显示,自旋回波扩散加权序列获得的图像基本没有发生形变,并且具有较高的分辨率和较好的信噪比.     

Improved B0 field map estimation for high field EPI

Echo planar imaging (EPI) is an ultrafast magnetic resonance imaging (MRI) technique that allows one to acquire a 2D image in about 100 ms. Unfortunately, the standard EPI images suffer from substantial geometric distortions, mainly originating from susceptibility differences in adjacent tissues. To reduce EPI distortions, correction methods based on a field map, which is a map of the off-resonance frequencies, have been developed. In this work, a nonlinear least squares estimator is used to optimize the estimation of the field map of the B₀ field. The model of the EPI and reference data includes parameters for the phase evolution, the complex magnitude, the relaxation of the MRI signal and the EPI-specific phase difference between odd and even echoes, and from these parameters, additional corrections might be computed. The reference data required to estimate the field map can be acquired with a modified EPI-sequence. The proposed method is tested on simulated as well as experimental data and proves to be significantly more robust against noise, compared to the previously suggested method.     

Accurate and sensitive measurements of magnetic susceptibility using echo planar imaging

Susceptibility differences are common causes for artifacts in magnetic resonance (MR); therefore, it is important to choose phantom materials in a way that these artifacts are kept at a minimum. In this study, a previously proposed MR imaging (MRI) method [Beuf O, Briguet A, Lissac M, Davis R. Magnetic resonance imaging for the determination of magnetic susceptibility of materials. J Magn Reson 1996; Series B(112):111-118] was improved to facilitate sensitive in-house measurements of different phantom materials so that such artifacts can more easily be minimized. Using standard MRI protocols and distilled water as reference, we measured magnetic volume susceptibility differences with a clinical MR system. Two imaging techniques, echo planar imaging (EPI) and spin echo, were compared using liquid samples whose susceptibilities were verified by MR spectroscopy. The EPI sequence has a very narrow bandwidth in the phase-encoding direction, which gives an increased sensitivity to magnetic field inhomogeneities. All MRI measurements were evaluated in two ways: (1) manual image analysis and (2) model fitting. The narrow bandwidth of the EPI made it possible to detect very small susceptibility differences (equivalent susceptibility difference, Deltachi(e)> or =0.02 ppm), and even plastics could be measured. Model fitting yielded high accuracy and high sensitivity and was less sensitive to other image artifacts as compared with manual image analysis.     

High-resolution diffusion imaging using phase-corrected segmented echo-planar imaging

Diffusion magnetic resonance imaging (MRI) was performed with a high-resolution segmented echo-planar imaging technique, which provided images with substantially less susceptibility artifacts than images obtained with single-shot echo-planar imaging (EPI). Diffusion imaging performed with any multishot pulse sequence is inherently sensitive to motion artifacts and in order to reduce motion artifacts, the presented method utilizes navigator echo phase corrections, performed after a one-dimensional Fourier transform along the frequency-encoding direction. Navigator echo phases were fitted to a straight line prior to phase correction to avoid errors from internal motion. In vivo imaging was performed using electro cardiographic (ECG) triggering. Apparent diffusion coefficient (ADC) maps were calculated on a pixel-by-pixel basis using up to seven diffusion sensitivities, ranging from b = 0 to 1129 x 10(6) s/m(2).     

Off-resonance correction using an estimated linear time map

Images acquired in the presence of magnetic field deviations and reconstructed without taking into account the off-resonance, are distorted and corrupted with artifacts. Several post-processing algorithms have been developed for correcting the distortion when it is not possible to fix the field inhomogeneities. These off-resonance correction methods are, in general, slow and computing intensive. To make them faster they are usually adapted to a particular situation or approximated. One of these approximations is to assume that the field map is linear. Although this assumption makes the algorithm fast and robust it is not well suited for arbitrary field maps. On the other hand, there are k-space trajectories with an almost linear time map (time at which each k-space value is acquired), such as 2DFT and EPI. This paper presents an algorithm for off-resonance correction based on a linear time map approximation. This approximation allows a fast algorithm that takes advantage of the almost linearity of the time map and uses the whole field map to correct the images. The proposed correction algorithm reduces the off-resonance induced artifacts while being fast. The linear approximation of the time map needs to be done only once for each trajectory because it does not depend on the acquired image or field map data. The method can also be extended to a multi-plane approximation for sequences with more complex time maps.     

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.     

Reduced distortion artifact whole brain CBF mapping using blip-reversed non-segmented 3D echo planar imaging with pseudo-continuous arterial spin labeling

PurposeTo implement and evaluate interleaved blip-up, blip-down, non-segmented 3D echo planar imaging (EPI) with pseudo-continuous arterial spin labeling (pCASL) and post-processing for reduced susceptibility artifact cerebral blood flow (CBF) maps.Materials and methods3D EPI non-segmented acquisition with a pCASL labeling sequence was modified to include alternating k-space coverage along phase encoding direction (referred to as “blip-reversed”) for alternating dynamic acquisitions of control and label pairs. Eight volunteers were imaged on a 3T scanner. Images were corrected for distortion using spatial shifting transformation of the underlying field map. CBF maps were calculated and compared with maps obtained without blip reversal using matching gray matter (GM) images from a high resolution 3D scan. Additional benefit of using the correction for alternating blip-up and blip-down acquisitions was assessed by comparing to corrected blip-up only and corrected blip-down only CBF maps. Matched Student t-test of overlapping voxels for the eight volunteers was done to ascertain statistical improvement in distortion.ResultsMean CBF value in GM for the eight volunteers from distortion corrected CBF maps was 50.8 ± 9.9 ml/min/100 gm tissue. Corrected CBF maps had 6.3% and 4.1% more voxels in GM when compared with uncorrected blip up (BU) and blip down (BD) images, respectively. Student t-test showed significant reduction in distortion when compared with blip-up images and blip-down images (p < 0.001). When compared with corrected BU and corrected BD only CBF maps, BU and BD corrected maps had 2.3% and 1% more voxels (p = 0.006 and 0.04, respectively).ConclusionPseudo-continuous arterial spin labeling with non-segmented 3D EPI acquisition using alternating blip-reversed k-space traversal and distortion correction provided significantly better matching GM CBF maps. In addition, employing alternating blip-reversed acquisitions during pCASL acquisition resulted in statistically significant improvement over corrected blip-up and blip-down CBF maps.     

基于极坐标变换去除计算机层析图像环形伪影

探测器像元响应的不一致、光源能不稳定等因素使得计算机层析(CT)图像中含有较多的环形伪影,严重降低了图像质量,影响了图像的三维重建和量化分析,为此提出了基于极坐标变换与傅里叶变换后低通滤波的算法去除环形伪影。通过极坐标变换将直角坐标下的环形伪影转化为极坐标下的线性伪影,然后对线性伪影图像进行傅里叶变换获得频谱图像,进而设计二维低通滤波器进行滤波处理,最后通过傅里叶逆变换与坐标逆变换获得校正后的图像。利用Matlab软件,编写程序对算法进行验证,结果表明,该算法能够有效地去除环形伪影,使图像内部细节清晰可见,并且保护了图像边缘信息,提高了图像的信噪比;另外,使用该方法处理100张切片图像只需3.5min,可满足批量处理的需求。     

Joint B0 and image estimation integrated with model based reconstruction for field map update and distortion correction in prostate diffusion MRI

In prostate Diffusion Weighted MRI, differences in susceptibility values exist at the interface between the prostate and rectal-air. This can result in off-resonance magnetic field leading to geometric distortions including signal stretching and signal pile-up in the reconstructed images. Using a set of EPI data acquired with blip-up and blip-down phase encoding gradient directions, model based reconstruction has recently been proposed that can correct these distortions by using a B₀ field estimated from a separate B₀ scan. However, change in the size of the rectal air region across time can occur that can result in a mismatch of the B₀ field to the EPI scan. Also, the measured B₀ field itself can be erroneous in regions of low Signal to Noise ratio around the prostate rectal air interface. In this work, using a set of single shot EPI data acquired with blip-up and blip-down phase encoding gradient directions, a novel joint model based reconstruction is proposed that can account for changes in the off resonance effects between the B₀ and EPI scans. For ten prostate patients, using a measured B₀ field as an initial B₀ estimate, on a 5-point scale (1–5) image quality scores evaluated by an experienced radiologist, the proposed framework achieved scores of 3.50 ± 0.85 and 3.40 ± 0.51 for b-values of 0 and 500 s/mm², respectively compared to 3.40 ± 0.70 and 3.30 ± 0.67 for model based reconstruction. The proposed framework is also capable of estimating a distortion corrected EPI image even without an initial B₀ field estimate in situations where a separate B₀ scan cannot be obtained due to time constraint.     

Parallel EPI artifact correction (PEAC) for N/2 ghost suppression in neuroimaging applications

Echo Planar Imaging (EPI) is a neuroimaging tool for clinical practice and research investigation. Due to odd-even echo phase inconsistencies, however, EPI suffers from Nyquist N/2 ghost artifacts. In standard neuroimaging protocols, EPI artifacts are suppressed using phase correction techniques that require reference data collected from a reference scan. Because reference-scan based techniques are sensitive to subject motion, EPI performance is sub-optimal in neuroimaging applications. In this technical note, we present a novel EPI data processing technique which we call Parallel EPI Artifact Correction (PEAC). By introducing an implicit data constraint associated with multi-coil sensitivity in parallel imaging, PEAC converts phase correction into a constrained problem that can be resolved using an iterative algorithm. This enables “reference-less” EPI that can improve neuroimaging performance. In the presented work, PEAC is investigated using a standard functional magnetic resonance imaging (fMRI) protocol with multi-slice 2D EPI. It is demonstrated that PEAC can suppress ghost artifacts as effectively as the standard reference-scan based phase correction technique used on a clinical MRI system. We also found that PEAC can achieve dynamic phase correction when motion occurs.     

Color Correction of Pathological Images Based on Dye Amount Quantification

Pathological images are color images of stained tissue slides, the color of which varies depending on staining conditions. For reliable diagnosis, the color variation must be corrected in these images. This paper proposes a color correction method for hematoxylin and eosin (H&E) stained pathological images in which the amounts of H&E dyes are estimated based on multispectral imaging technique and Beer Lambert law, and the color image is generated corresponding to the adjusted amount of dyes. This enables us to correct an image to an arbitrary or specified optimal staining-condition image. Through experiments using H&E stained human liver slide images, the effectiveness of the proposed method was confirmed.     

低能X射线工业CT图像杯状伪影校正

为了去除X射线工业CT图像中的杯状伪影,提高CT图像的识别能力和量化分析精度,提出一种基于分度投影和权函数的射束硬化校正方法。首先分析得出杯状伪影主要是由X射线连续谱穿过被测物体过程中出现的射束硬化所导致。然后扫描阶梯模型,采集不同厚度下的投影数据并求出线衰减系数,通过拟合曲线,得到硬化模型函数和权函数校正模型函数,并确定权函数。接着,扫描被测圆柱形工件,采集不同分度下的投影数据。最后,针对每一个分度投影数据,采用权函数与当前分度投影数据乘积的方法进行硬化校正。对含有杯状伪影的实际CT图像进行了校正实验,结果表明,与多项式拟合法相比,该方法校正后的灰度图像没有放大噪声,且信噪比提高3.29%,有效地消除了杯状伪影,同时较好地保留了图像边界细节。     

Susceptibility, field inhomogeneity, and chemical shift-corrected NMR microscopy: Application to the human finger in vivo

Spectroscopic proton image data recorded with the aid of a gradient-echo spectroscopic imaging pulse sequence are reported. A postdetection processing method is suggested which permits correction of artifacts due to inhomogeneity, susceptibility, and chemical-shift resonance offsets. That is, apart from the spectral information available in this way, better spatial resolutions can be achieved. The method is demonstrated by resonance-offset corrected images of the human finger in vivo. Moreover, resonance-line selective and spectroscopically resolved diffusion-weighted images and diffusivity maps rendered with the aid of the same postdetection procedure are shown.     

A dual image approach for bias field correction in magnetic resonance imaging

In this paper, we propose a dual image approach to correcting intensity inhomogeneities for MR images acquired using surface coils. Previous methods are usually not satisfactory due to restricted application domains, considerable human interactions, or some undesirable artifacts. The proposed algorithm provides nice correction results for a variety of surface-coil MR images. It is accomplished by using an additional body-coil MR image of a smaller size captured at the same position as that of the surface-coil image to facilitate the estimation of the bias field function. The correction algorithm consists of aligning the surface-coil image with the body-coil image and fitting a spline surface from a sparse set of data points for the associated bias field function. Experiments on some real images show satisfactory correction results by using the proposed algorithm.     

颅脑梯度回波成像:呼吸伪影和导航回波矫正

梯度回波序列是磁共振成像中常用的脉冲序列,然而梯度回波对主磁场波动非常敏感,呼吸等生理运动引起的信号波动会导致图像伪影.该文报道了采用导航回波技术获取呼吸运动导致的局部磁场波动,用以矫正图像回波中随时间变化的相位波动,并将该技术应用于三维多回波梯度回波成像和T2*定量图研究.研究结果显示:矫正前,相位波动幅度随回波时间增长而增大,模图和T2*定量图在相位编码方向有明显伪影,并且男女呼吸伪影水平有显著性差异;矫正后,相位波动幅度大幅下降,图像伪影水平有显著性下降.     

Correction of concomitant gradient artifacts in experimental microtesla MRI

Magnetic resonance imaging (MRI) suffers from artifacts caused by concomitant gradients when the product of the magnetic field gradient and the dimension of the sample becomes comparable to the static magnetic field. To investigate and correct for these artifacts at very low magnetic fields, we have acquired MR images of a 165-mm phantom in a 66-microT field using gradients up to 350 microT/m. We prepolarize the protons in a field of about 100 mT, apply a spin-echo pulse sequence, and detect the precessing spins using a superconducting gradiometer coupled to a superconducting quantum interference device (SQUID). Distortion and blurring are readily apparent at the edges of the images; by comparing the experimental images to computer simulations, we show that concomitant gradients cause these artifacts. We develop a non-perturbative, post-acquisition phase correction algorithm that eliminates the effects of concomitant gradients in both the simulated and the experimental images. This algorithm assumes that the switching time of the phase-encoding gradient is long compared to the spin precession period. In a second technique, we demonstrate that raising the precession field during phase encoding can also eliminate blurring caused by concomitant phase-encoding gradients; this technique enables one to correct concomitant gradient artifacts even when the detector has a restricted bandwidth that sets an upper limit on the precession frequency. In particular, the combination of phase correction and precession field cycling should allow one to add MRI capabilities to existing 300-channel SQUID systems used to detect neuronal currents in the brain because frequency encoding could be performed within the 1-2 kHz bandwidth of the readout system.     

自导航快速自旋回波在低场MRI上的实现

低场磁共振成像仪一般需采用数据累加的办法来提高图像信噪比,这样会延长扫描时间,因此更易受运动伪影的影响. 为了解决运动伪影问题,本文在低场磁共振成像仪上实现了自导航快速自旋回波去运动伪影成像技术,并且与常规快速自旋回波序列进行了临床对比实验. 结果表明,与常规快速自旋回波序列相比,采用自导航快速自旋回波技术后,由于病人运动导致的伪影得到明显地抑制. 