衍射增强成像提取多种信息的简便方法研究

衍射增强相位衬度成像中的信息分离研究一直以来都是相关研究人员重点关注的方向之一. 本文在利用余弦函数拟合衍射增强成像中的摇摆曲线的基础上, 推导出了峰位像、左腰像和右腰像的余弦函数表达式以及吸收像、折射角像和散射角方差像的解析表达式, 形成了一种基于衍射增强成像的简便信息分离方法. 该方法只需利用摇摆曲线左腰、右腰和峰位三幅图像, 就能提取出样品的吸收、折射和散射信息, 具有方法简便、样品所受辐射剂量低等优点. 对模型样品和真实样品的实验结果表明, 本文所提方法可成功对样品进行信息分离, 并且所获得的信息分离结果可以和目前常用的多图统计方法(至少需要7张图像)获得的实验结果相比较. 关键词： 衍射增强成像 相位衬度成像 样品信息提取     

X射线衍射增强成像中的折射衬度

根据X射线在界面处的折射过程,推导出X射线在物体中折射偏转角度的表达式;从X射线双晶衍射的原理出发,讨论了衍射增强成像(DEI)技术的空间分辨率和电子密度的衬度分辨率,明确了折射衬度在DEI技术中的表现形式,并从实验上进行了验证.利用不同能量的X射线获得了人肝病理样品中血管的DEI清晰图像,讨论了折射衬度在其中所起的作用.     

X射线衍射增强成像中吸收、折射以及散射衬度的计算层析

X射线相位衬度成像是一种新型的X射线成像技术,通过记录射线穿过物体后相位的改变对物体进行成像,可以提供比传统的X射线吸收成像更高的图像衬度以及空间分辨力。衍射增强成像方法(Diffraction enhancedimaging,DEI)是X射线相位衬度成像方法之一,利用一块放置在物体和探测器之间的分析晶体提取物体的吸收、折射以及散射信息并进行成像。将衍射增强成像方法与计算机断层成像技术(Computerized Tomography)进行结合,利用吸收、散射以及折射信息,分别采用滤波反投影以及雷登(Radon)变换,对昆虫样品——蜜蜂进行计算层析重建,获得了好于X射线吸收计算层析的重建结果,验证了衍射增强成像信息分离计算层析的优越性。     

PITRE: software for phase‐sensitive X‐ray image processing and tomography reconstruction

Synchrotron‐radiation computed tomography has been applied in many research fields. Here, PITRE (Phase‐sensitive X‐ray Image processing and Tomography REconstruction) and PITRE_BM (PITRE Batch Manager) are presented. PITRE supports phase retrieval for propagation‐based phase‐contrast imaging/tomography (PPCI/PPCT), extracts apparent absorption, refractive and scattering information of diffraction enhanced imaging (DEI), and allows parallel‐beam tomography reconstruction for conventional absorption CT data and for PPCT phase retrieved and DEI‐CT extracted information. PITRE_BM is a batch processing manager for PITRE: it executes a series of tasks, created viaPITRE, without manual intervention. Both PITRE and PITRE_BM are coded in Interactive Data Language (IDL), and have a user‐friendly graphical user interface. They are freeware and can run on Microsoft Windows systems via IDL Virtual Machine, which can be downloaded for free and does not require a license. The data‐processing principle and some examples of application will be presented.     

DEIReconstructor:a software for diffraction enhanced imaging processing and tomography reconstruction

Diffraction enhanced imaging (DEI) has been widely applied in many fields, especially when imaging low-Z samples or when the difference in the attenuation coefficient between different regions in the sample is too small to be detected. Recent developments of this technique have presented a need for a new software package for data analysis. Here, the Diffraction Enhanced Image Reconstructor (DEIReconstructor), developed in Matlab, is presented. DEIReconstructor has a user-friendly graphical user interface and runs under any of the 32-bit or 64-bit Microsoft Windows operating systems including XP and Win7. Many of its features are integrated to support imaging preprocessing, extract absorption, refractive and scattering information of diffraction enhanced imaging and allow for parallel-beam tomography reconstruction for DEI-CT. Furthermore, many other useful functions are also implemented in order to simplify the data analysis and the presentation of results. The compiled software package is freely available.     

同步辐射硬X射线衍射增强峰位成像CT研究

对于主要由轻元素组成的生物、医学样品,利用衍射增强成像技术可以观察到常规吸收成像无法观察到的内部微观结构,因而衍射增强成像具有较高的衬度和空间分辨率.选用苍蝇作为实验样品,在北京同步辐射装置上首次开展了衍射增强峰位成像CT实验.重建出的样品断层像显示其分辨率达到了几十μm水平.对于将衍射增强成像技术应用于生物和医学等领域具有重要意义.     

X射线衍射增强成像中的定量测量

提出了一种基于X射线衍射增强成像(DEI)断层计算机X射线层析术(CT)图像的物体尺寸精确测量方法.X射线衍射增强成像是一种基于相位衬度的成像技术.通过建立DEI的简化模型,研究衍射成像过程中品体转角、投影图像谷点位置、成像系统等效模糊等因素之间关系,由此具体探讨了系统模糊效应对圆物体边界成像带来的位置偏移,并以圆形被测样品为例.提出可精确测定直径的简单有效算法.通过理论仿真模型数据和北京同步辐射装置上的实测数据验证了该算法的精度.该方法实现了利用DEI图像对被测物体几何尺寸的精确测量,可用于对牛物组织样品等物体内部微小结构的尺寸的精确测量.     

衍射增强成像方法在计算机断层成像中的应用

衍射增强成像是X射线相位衬度成像方法之一，这种方法具有较高的衬度和空间分辨率. 利用它对由轻元素组成的生物、医学样品成像可以观察到常规吸收成像无法观察到的内部微细结构. 这种方法在生物、医学、材料科学等无损检测领域中的应用研究，已成为当前国际上X射线成像领域中的研究热点. 讨论衍射增强成像方法和该方法在计算机断层成像中的应用. 实验结果表明，使用衍射增强成像方法获得的数据源能够重建出微米级的生物组织结构，这些组织结构信息在常规X射线断层成像中是难以得到的. 关键词： 衍射增强成像 CT重建 同步辐射 微细结构     

Preliminary attempt on maximum likelihood tomosynthesis reconstruction of DEI data

Tomosynthesis is a three-dimension reconstruction method that can remove the effect of superimposition with limited angle projections. It is especially promising in mammography where radiation dose is concerned. In this paper, we propose a maximum likelihood tomosynthesis reconstruction algorithm (ML-TS) on the apparent absorption data of diffraction enhanced imaging (DEI). The motivation of this contribution is to develop a tomosynthesis algorithm in low-dose or noisy circumstances and make DEI get closer to clinic application. The theoretical statistical models of DEI data in physics are analyzed and the proposed algorithm is validated with the experimental data at the Beijing Synchrotron Radiation Facility (BSRF). The results of ML-TS have better contrast compared with the well known `shift-and-add' algorithm and FBP algorithm.

     

Improvement and error analysis of quantitative information extraction in diffraction-enhanced imaging

Diffraction-enhanced imaging （DEI） is a powerful phase-sensitive technique that provides higher spatial resolution and supercontrast of weakly absorbing objects than conventional radiography. It derives contrast from the X-ray absorption, refraction, and ultra-small-angle X-ray scattering （USAXS） properties of an object. The separation of different-contrast contributions from images is an important issue for the potential application of DEI. In this paper, an improved DEI （IDEI） method is proposed based on the Gaussian curve fitting of the rocking curve （RC）. Utilizing only three input images, the IDEI method can accurately separate the absorption, refraction, and USAXS contrasts produced by the object. The IDEI method can therefore be viewed as an improvement to the extended DEI （EDEI） method. In contrast, the IDEI method can circumvent the limitations of the EDEI method well since it does not impose a Taylor approximation on the RC. Additionally, analysis of the IDEI model errors is performed to further investigate the factors that lead to the image artifacts, and finally validation studies are conducted using computer simulation and synchrotron experimental data.     

X-ray diffraction enhanced imaging study of intraocular tumors in human beings

Diffraction enhanced imaging(DEI) with edge enhancement is suitable for the observation of weakly absorbing objects.The potential ability of the DEI was explored for displaying the microanatomy and pathology of human eyeball in this work.The images of surgical specimens from malignant intraocular tumor of hospitalized patients were taken using the hard X-rays from the topography station of Beamline 4W1A at Beijing Synchrotron Radiation Facility(BSRF).The obtained radiographic images were analyzed in correlation with those of pathology.The results show that the anatomic and pathologic details of intraocular tumors in human beings can be observed clearly by DEI for the first time,with good visualization of the microscopic details of eyeball ring such as sclera,choroids and other details of intraocular organelles.And the best resolution of DEI images reaches up to the magnitude of several tens of μm.The results suggest that it is capable of exhibiting clearly the details of intraocular tumor using DEI method.     

Optimization-based image reconstruction in x-ray computed tomography by sparsity exploitation of local continuity and nonlocal spatial self-similarity

The additional sparse prior of images has been the subject of much research in problems of sparse-view computed tomography(CT) reconstruction. A method employing the image gradient sparsity is often used to reduce the sampling rate and is shown to remove the unwanted artifacts while preserve sharp edges, but may cause blocky or patchy artifacts.To eliminate this drawback, we propose a novel sparsity exploitation-based model for CT image reconstruction. In the presented model, the sparse representation and sparsity exploitation of both gradient and nonlocal gradient are investigated.The new model is shown to offer the potential for better results by introducing a similarity prior information of the image structure. Then, an effective alternating direction minimization algorithm is developed to optimize the objective function with a robust convergence result. Qualitative and quantitative evaluations have been carried out both on the simulation and real data in terms of accuracy and resolution properties. The results indicate that the proposed method can be applied for achieving better image-quality potential with the theoretically expected detailed feature preservation.     

小入射角条件下气动光学成像偏移

针对超声速飞行器在大气中飞行时由于气动光学效应产生的成像偏移,推导了光线偏折角和流场折射率梯度之间的关系,给出了小入射角情况下,离散折射率场光线偏折角的一种计算方法;利用不同流场的计算流体力学数据计算出对应的折射率场;使用光线追迹法得到光线通过不同流场后的偏折角和传播路径上的折射率梯度分布。仿真结果表明,新的偏折角计算方法所得结果在入射角小于30°时与Runge-Kutta和Snell光线追迹法具有很好的一致性;同一高度下小入射角时,光线的偏折角并不随着飞行速度的增加而增加,而是基本保持稳定;光线通过流场的偏折是正负折射率梯度区域共同作用的结果,负的折射率梯度区域可以减少光线偏折,起到部分"校正"作用。     

衍射增强成像信息分离方法研究

衍射增强成像（DEI）是一种功能强大的相位衬度成像技术,对于软组织,它能获得比基于吸收的传统X射线成像技术更高的衬度.研究了一种基于DEI的信息分离方法,它能同时获得三种参数图像,即吸收像、折射像和小角散射像.信息分离利用解析方程求解,获得的参数图像质量取决于在三个不同的分析晶体角度处获得的DEI图像.利用矩阵条件数对方程求解的稳定性进行了讨论,推导出理论上的折射角和小角散射宽度的最大值.最后,利用豚鼠耳蜗DEI图像对研究结果进行了验证. 关键词： 信息分离 衍射增强成像 矩阵条件数 豚鼠耳蜗     

基于先验组分的多谱CT序列DCM融合算法研究

多谱CT成像是通过不同谱段的CT图像表征检测对象中的不同组分。为了便于在同一视图中显示所有组分的信息,需要研究多谱CT序列的融合方法;但是常用融合方法如加权平均法、小波变换融合法等都是针对图像细节信息的优化,不能表达组分的物理特性,从而导致融合图像的灰度不具有物理表征性,影响CT的定量检测。为此,结合具有物理表征特征的数据约束模型(DCM),开展了基于先验组分的多谱CT序列DCM融合算法研究。首先通过能谱滤波分离的成像方法获得多个能谱范围内的多能投影数据,采用TV-OSEM算法重建不同能谱段的CT序列;其次,利用传统DCM模型和改进DCM模型分别对多谱CT序列进行融合,传统DCM模型是严格单能的,由于滤波后能谱的非严格单能特性,其融合结果不能表征出对象序列中的全部组分。针对此问题提出了改进DCM模型。改进DCM模型选择了新的体元定义,并且在多谱CT序列融合中引入先验组分作为参照,通过先验物质对融合结果中其他物质进行校准,实现检测对象中各组分位置的准确分布。仿真实验表明,该方法可从物理表征正确性的角度,实现多谱CT序列融合,在满足CT序列中不同组分区分的同时,其融合图像的灰度具有物理可参照性,有利于后续的CT定量检测。     

A new method for information retrieval in two-dimensional grating-based X-ray phase contrast imaging

Grating-based X-ray phase contrast imaging has been demonstrated to be an extremely powerful phase-sensitive imaging technique.By using two-dimensional(2D) gratings,the observable contrast is extended to two refraction directions.Recently,we have developed a novel reverse-projection(RP) method,which is capable of retrieving the object information efficiently with one-dimensional(1D) grating-based phase contrast imaging.In this contribution,we present its extension to the 2D grating-based X-ray phase contrast imaging,named the two-dimensional reverseprojection(2D-RP) method,for information retrieval.The method takes into account the nonlinear contributions of two refraction directions and allows the retrieval of the absorption,the horizontal and the vertical refraction images.The obtained information can be used for the reconstruction of the three-dimensional phase gradient field,and for an improved phase map retrieval and reconstruction.Numerical experiments are carried out,and the results confirm the validity of the 2D-RP method.     

Novel Fourier-based iterative reconstruction for sparse fan projection using alternating direction total variation minimization

Sparse-view x-ray computed tomography(CT) imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in reconstruction speed and memory usage. A novel Fourier-based iterative reconstruction technique that utilizes non-uniform fast Fourier transform(NUFFT) is presented in this work along with advanced total variation(TV) regularization for a fan sparse-view CT. The proposition of a selective matrix contributes to improve reconstruction quality. The new method employs the NUFFT and its adjoin to iterate back and forth between the Fourier and image space. The performance of the proposed algorithm is demonstrated through a series of digital simulations and experimental phantom studies. Results of the proposed algorithm are compared with those of existing TV-regularized techniques based on compressed sensing method, as well as basic algebraic reconstruction technique. Compared with the existing TV-regularized techniques, the proposed Fourier-based technique significantly improves convergence rate and reduces memory allocation, respectively.     

Quantitative quasi-local tomography using absorption and phase contrast

Possibilities are investigated for quantitative computed tomography (CT) reconstruction of the spatial distribution of refractive index in a region of interest (ROI) inside an object which is larger than the field of view. The analysis is presented for the case of conventional (absorption) CT, as well as diffraction (phase-contrast) CT. It is shown that in both cases an accurate reconstruction can be achieved using the projection data corresponding to rays passing through a sufficiently wide vicinity of a ROI. Our analysis also indicates that X-ray phase-contrast CT can typically be localized to smaller regions compared to absorption CT. In particular, quasi-local hard X-ray micro-CT of regions of the order of 100 μm or even smaller in size, appears feasible with the use of propagation-based phase contrast. Numerical and experimental tests confirm the possibility of accurate quantitative CT reconstruction from truncated projections.     

A compressed sensing approach for enhancing infrared imaging resolution

This paper presents a novel approach for improving infrared imaging resolution by the use of Compressed Sensing (CS). Instead of sensing raw pixel data, the image sensor measures the compressed samples of the observed image through a coded aperture mask placed on the focal plane of the optical system, and then the image reconstruction can be conducted from these samples using an optimal algorithm. The resolution is determined by the size of the coded aperture mask other than that of the focal plane array (FPA). The attainable quality of the reconstructed image strongly depends on the choice of the coded aperture mode. Based on the framework of CS, we carefully design an optimum mask pattern and use a multiplexing scheme to achieve multiple samples. The gradient projection for sparse reconstruction (GPSR) algorithm is employed to recover the image. The mask radiation effect is discussed by theoretical analyses and numerical simulations. Experimental results are presented to show that the proposed method enhances infrared imaging resolution significantly and ensures imaging quality.     

结合计算层析与光学相干层析测量折射率分布

通过将计算层析技术与光学相干层析技术相结合,测量散射介质非均匀折射率分布。该方法测量散射样品折射率分布时,通过光学相干层析技术测量散射样品折射率分布在多个方向上的投影,采用计算层析技术实现对样品折射率分布的层析重建,克服了传统折射率光学测量方法如单纯的基于光学相干层析原理的焦点追迹法和光程匹配法,不能测量散射介质非均匀折射率分布的缺点。采用该方法在实验中对梯度折射率棒的径向相对折射率分布进行层析重建,取得了较好的实验结果。