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.     

High resolution diffusion weighted magnetic resonance imaging of the pancreas using reduced field of view single-shot echo-planar imaging at 3 T

Diffusion weighted magnetic resonance imaging (DWI) has been mostly acquired using single-shot echo-planar imaging (ss EPI) to minimize motion induced artifacts. The spatial resolution, however, is inherently limited in ss EPI especially for abdominal imaging, even with the advances in parallel imaging. A novel method of reduced Field of View ss EPI (rFOV ss EPI) has achieved high resolution DWI in human carotid artery, spinal cord with reduced blurring and higher spatial resolution than conventional ss EPI, but it has not been used to pancreas imaging. In the work, comparisons between the full FOV ss-DW EPI and rFOV ss-DW EPI in image qualities and ADC values of pancreatic tumors and normal pancreatic tissues were performed to demonstrate the feasibility of pancreatic high resolution rFOV DWI. There were no significant differences in the mean ADC values between full FOV DWI and rFOV DWI for the 17 subjects using b = 600 s/mm² (P = 0.962). However, subjective scores of image quality was significantly higher at rFOV ss DWI (P = 0.008 and 0.000 for b-value = 0 s/mm² and 600 s/mm² respectively). The spatial resolution of DWI for pancreas was increased by a factor of over 2.0 (from almost 3.0 mm/pixel to 1.25 mm/pixel) using rFOV ss EPI technique. Reduced FOV ss EPI can provide good DW images and is promising to benefit applications for pancreatic diseases.     

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 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.     

300万像素手机镜头设计

设计一种300万像素的手机镜头。该镜头由3片塑料非球面透镜和1个红外滤光片组成,采用三星公司的一款300万像素0.635cm(1/4inch)CMOS作为该镜头的图像传感器,像素颗粒大小为1.75μm,其分辨率极限为285lp/mm,即为奈奎斯特频率。镜头的光圈值F为2.85,视场2ω为62°。该镜头有较好的成像质量,在奈奎斯特频率1/2处绝大部分视场的MTF值大于0.5,波前均方差(RMS wavefront error)小于0.14λ(λ为波长),最大畸变为-0.2%。     

Adaptive field-of-view imaging for efficient receive beamforming in medical ultrasound imaging systems

Quadrature demodulation-based phase rotation beamforming (QD-PRBF) is commonly used to support dynamic receive focusing in medical ultrasound systems. However, it is computationally demanding since it requires two demodulation filters for each receive channel. To reduce the computational requirements of QD-PRBF, we have previously developed two-stage demodulation (TSD), which reduces the number of lowpass filters by performing demodulation filtering on summation signals. However, it suffers from image quality degradation due to aliasing at lower beamforming frequencies. To improve the performance of TSD-PRBF with reduced number of beamforming points, we propose a new adaptive field-of-view (AFOV) imaging method. In AFOV imaging, the beamforming frequency is adjusted depending on displayed FOV size and the center frequency of received signals. To study its impact on image quality, simulation was conducted using Field II, phantom data were acquired from a commercial ultrasound machine, and the image quality was quantified using spatial (i.e., axial and lateral) and contrast resolution. The developed beamformer (i.e., TSD-AFOV-PRBF) with 1024 beamforming points provided comparable image resolution to QD-PRBF for typical FOV sizes (e.g., 4.6% and 1.3% degradation in contrast resolution for 160 mm and 112 mm, respectively for a 3.5 MHz transducer). Furthermore, it reduced the number of operations by 86.8% compared to QD-PRBF. These results indicate that the developed TSD-AFOV-PRBF can lower the computational requirement for receive beamforming without significant image quality degradation.     

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.     

Optical design of a full-pupil field,non-contact and high resolution corneal curved objective lens

It has been a challenge to overcome the corneal curvature radius to design a full-pupil field, non-contact and high resolution corneal curved objective lens, which covers the cornea full-pupil field and has the ability to resolve corneal cells. In this paper, we report an optical design of a full-pupil field, non-contact corneal curved objective lens for high resolution cornea imaging. The advantages of this lens are that it has a wide field of view (FOV) with the corneal curved image surface, maintains the beam normal incidence, as well as non-contact lens imaging, and offers a cell-level lateral resolution of cornea structure. The analysis of optimization shows that the system achieves diffraction limit in a circular FOV of 4 mm diameter covering the full-pupil zone. The theoretical lateral resolution is about 2.5 μm with an image space NA of 0.16, which is sufficient to resolve corneal cells of 7 μm diameter, and the working distance is larger than 15 mm which is enough for a non-contact objective lens. So the optical design is effectively and efficiently meeting the demand of specifications.     

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."     

Distortion-free single point imaging of multi-layered composite sandwich panel structures

The results of a magnetic resonance imaging (MRI) investigation concerning the effects of an aluminum honeycomb sandwich panel on the B1 and B0 fields and on subsequent image quality are presented. Although the sandwich panel structure, representative of an aircraft composite material, distorts B0 and attenuates B1, distortion-free imaging is possible using single point (constant time) imaging techniques. A new expression is derived for the error caused by gradient field distortion due to the heterogeneous magnetic susceptibility within a sample and this error is shown not to cause geometric distortion in the image. The origin of the B0 distortion in the sample under investigation was also examined. The graphite-epoxy 'skin' of the panel is the principal source of the B0 distortion. Successful imaging of these structures sets the stage for the development of methods for detecting moisture ingress and degradation within composite sandwich structures.     

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.     

500万像素手机镜头设计

 利用ZEMAX光学工程设计软件,设计了一款500万像素的手机镜头。该镜头由4片塑料非球面透镜和1片红外滤光片组成,其光圈值F为2.85,视场2ω为60°,采用Aptina公司的一款500万像素7.94mm(1/3.2)英寸CMOS作为该镜头的图像传感器,该图像传感器的像素颗粒大小为1.75μm,截止频率为285lp/mm,即为奈奎斯特频率。设计结果显示,该镜头在奈奎斯特频率处,0.7视场以内的MTF值大于0.3,在奈奎斯特频率1/2处视场的MTF值均大于0.5,波前均方差(RMS wavefront error)小于0.1,畸变小于1%。     

小视场长焦距镜头畸变高精度测量研究

利用由精密测角仪、显微摄像测量系统、微型双光栅平面干涉仪、平行光管以及星点组成畸变测量系统,对小视场长焦距的镜头进行畸变测量。在计算镜头畸变中,利用中心视场区域内畸变设计无穷小,采用三次多项式拟和的方法,计算镜头理论焦距;在边视场采用像高高次方和视场角正弦高次方加权平均的方法对测量偏差角进行修正,得到各视场的相对和绝对畸变。通过实际测量和计算验证,镜头全视场畸变测量精度可达到0.02%。     

Functional magnetic resonance imaging in a stereotactic setup

The localization of critical structures within the brain is important for the planning of therapeutic strategies. Functional MRI is capable to assess functional response of cortical structures to certain stimuli. The authors present two techniques for functional MRI (fMRI) in a stereotactic set-up. The skull of the patients has been immobilized for stereotactic treatment planning either with a self developed stereotactic ceramic frame and bony fixation or with an individual precision mask system made of light cast. It has been shown that this frame does not produce any image distortion. fMRI was performed using a modified FLASH sequence on a conventional 1.5 T MRI scanner with a specially developed linear polarized head coil. The imaging technique used was an optimized conventional 2D and 3D, first order flow rephased, gradient echo sequence (FLASH) with fat-suppression and reduce bandwidth (16–28 Hz/pixel) and TR = 80–120 ms, TE = 60 ms, flip angle = 40°, matrix = 128 × 128, FOV = 150–250 mm, slice-thickness = 2–5 mm, NEX = 1, and a total single scan time for one image of about 7 sec. The motor cortex stimulation was achieved by touching each finger to thumb in a sequential, self-paced, and repetitive manner. Statistical parametric maps based on student's test were calculated. Pixels with a highly significant signal increase (p < 0.001) are overlaid on T1w SE images. The primary motor and sensory cortex could be visualized with this method in all 10 patients that were imaged in this study. Due to tight fixation of the patient's skull there have been no motion artifacts. These results show that functional MRI is feasible in an stereotactic set-up with an standard 1.5 T scanner. This is a prerequisite for the exact pre therapeutic assessment of the function of cortical centers.     

基于共心球透镜的多尺度广域高分辨率计算成像系统设计

针对实时广域高分辨率成像需求,充分利用具有对称结构的多层共心球透镜视场大且各轴外视场成像效果一致性好的特点,设计基于共心球透镜的多尺度广域高分辨率计算成像系统.该系统基于计算成像原理,通过构建像差优化函数获得光学系统设计参数,结合球形分布的次级相机阵列进行全局性优化,提高系统性能的同时有效简化光学设计过程、降低系统设计难度.系统稳定性测试结果表明,该成像系统的MTF(modulation transmission function)值在截止频率处接近衍射极限,弥散斑均方根恒小于探测器像元尺寸,整机实景实时成像效果良好,无视觉可见畸变.该系统不仅有效解决了传统成像中广域和高分辨率成像矛盾的问题,而且为计算光学成像系统设计奠定了一定研究基础.     

Chemical exchange saturation transfer MR imaging of articular cartilage glycosaminoglycans at 3 T: Accuracy of B0 Field Inhomogeneity corrections with gradient echo method

Glycosaminoglycan Chemical Exchange Saturation Transfer (gagCEST) is an important molecular MRI methodology developed to assess changes in cartilage GAG concentrations. The correction for B₀ field inhomogeneity is technically crucial in gagCEST imaging. This study evaluates the accuracy of the B₀ estimation determined by the dual gradient echo method and the effect on gagCEST measurements. The results were compared with those from the commonly used z-spectrum method. Eleven knee patients and three healthy volunteers were scanned. Dual gradient echo B₀ maps with different ?TE values (1, 2, 4, 8, and 10 ms) were acquired. The asymmetry of the magnetization transfer ratio at 1 ppm offset referred to the bulk water frequency, MTR_asym(1 ppm), was used to quantify cartilage GAG levels. The B₀ shifts for all knee patients using the z-spectrum and dual gradient echo methods are strongly correlated for all ?TE values used (r = 0.997 to 0.786, corresponding to ?TE = 10 to 1 ms). The corrected MTR_asym(1 ppm) values using the z-spectrum method (1.34% ± 0.74%) highly agree only with those using the dual gradient echo methods with ?TE = 10 ms (1.72% ± 0.80%; r = 0.924) and 8 ms (1.50% ± 0.82%; r = 0.712). The dual gradient echo method with longer ?TE values (more than 8 ms) has an excellent correlation with the z-spectrum method for gagCEST imaging at 3 T.     

用于监控系统的鱼眼镜头光学设计

为了满足监控系统中单镜头可实现全景监控的需求,运用仿生学原理,在原有鱼眼镜头的基础上结合监控系统的需求,运用ZEMAX软件设计了一款视场为180°,相对孔径为1/1.6的鱼眼镜头.该镜头由9片透镜组成,总长度为68.5mm,采用1/3英寸CCD作为图像接收器件.在120lp/mm时,其MTF曲线在轴上大于0.4,在80...     

机载双视场中波红外光学系统优化设计

为了满足机载红外搜索与跟踪系统的实嘎使用要求,根据变焦系统的基本理论及中波红外系统的特点,设计了320pixel×256pixel的中波制冷型焦平面阵列探测器的双视场中波红外光学系统。系统采用了二次成像结构,并在系统的第一像面位置安装光阑,以此减小杂散光对系统的影响。设计结果表明：系统具有100％冷光栏效率,在仅移动一片透镜的情况下可实现在800和400mln的两档变焦,系统F数为4且恒定不变,像面保持稳定,系统场曲〈0．04mm,畸变〈2．5％,在探测器的Nyquist频率16lp／mm处光学传递函数的峰值〉0．5,表明光学系统的像质满足使用要求。     

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.     

星载均匀像面低畸变广角气溶胶探测仪的研制

为满足空间遥感的迫切需求,设计并研制了一个星载均匀像面低畸变广角气溶胶探测仪样机.该仪器通过利用光阑像差产生的有效像差渐晕提高像面照度的均匀性,解决了广角系统中像面照度不均匀的问题.合理选择结构型式校正了畸变,并且采用全球面光学系统,易于加工和检测.广角气溶胶探测仪的中心波长为670 nm,带宽20 nm,全视场72°,相对孔径1/3.6,焦距20 mm.实验结果表明:研制的星载广角气溶胶探测仪镜头其入瞳大小5.6 mm,边缘视场的相对照度达到95.6%,在36 lp/mm处,轴上视场的调制传递函数值大于0.61,轴外视场的调制传递函数值高于0.58,最大畸变量为-1.95%,完全满足设计指标要求,体积小,适合空间遥感应用.