Image reconstruction algorithm for diffraction enhanced imaging-based computed tomography

A number of reconstruction algorithms for differential phase-contrast computed tomography reconstruct the refractive index decrement using the refraction angles of X-rays. However, these approaches require that the refraction angle at each view be extracted from several raw images captured at the same view, leading to unacceptably long exposure times and huge X-ray doses. We develop an iterative reconstruction algorithm to reconstruct the complex refractive index from raw images gathered only once at each view angle. Using the Tikhonov regularization method as basis, we simultaneously obtain the refractive index decrement and attenuation coefficient by minimizing the cost function. The experimental results show that the proposed algorithm effectively yields reconstructed images of low noise level.     

Cosine fitting radiography and computed tomography

A new method in diffraction-enhanced imaging computed tomography(DEI-CT) that follows the idea developed by Chapman et al.[Chapman D,Thomlinson W,Johnston R E,Washburn D,Pisano E,Gmur N,Zhong Z,Menk R,Arfelli F and Sayers D 1997 Phys.Med.Biol.42 2015] in 1997 is proposed in this paper.Merged with a "reverse projections"algorithm,only two sets of projection datasets at two defined orientations of the analyzer crystal are needed to reconstruct the linear absorption coefficient,the decrement of the real part of the refractive index and the linear scattering coefficient of the sample.Not only does this method reduce the delivered dose to the sample without degrading the image quality,but,compared with the existing DEI-CT approaches,it simplifies data-acquisition procedures.Experimental results confirm the reliability of this new method for DEI-CT applications.     

Cell refractive index tomography by digital holographic microscopy

For what we believe to be the first time, digital holographic microscopy is applied to perform optical diffraction tomography of a pollen grain. Transmission phase images with nanometric axial accuracy are numerically reconstructed from holograms acquired for different orientations of the rotating sample; then the three-dimensional refractive index spatial distribution is computed by inverse radon transform. A precision of 0.01 for the refractive index estimation and a spatial resolution in the micrometer range are demonstrated.     

基于单元件干涉仪的计算机断层扫描重建光纤三维折射率分布

提出了一种基于单元件干涉仪的计算机断层扫描方案,用于测量光纤三维折射率分布。该单元件干涉仪基于显微成像的原理,提高了系统的横向分辨率。利用快速傅里叶变换提取相位,采用滤波反投影算法重建了光纤的折射率分布。搭建了实验测量装置,实际测量了单模和多模光纤的折射率分布。结果表明,提供的方法可以简单快速地得到全光场数据,同时具有无损和非接触的优点,并且光路结构紧凑稳定,可为计算机断层扫描设备的小型化提供一种新的思路。     

Single-shot grating-based x-ray differential phase contrast imaging with a modified analyzer grating

X-ray grating interferometer has attracted widely attention in the past years due to its capability in achieving x-ray phase contrast imaging with low brilliance source. However, the widely used phase stepping information extraction method reduces system stability and prolongs data acquisition time by several times compared with conventional x-ray absorptionbased imaging. The mechanical stepping can be avoided by using a staggered grating, but at the cost of low vertical spatial resolution. In this paper, employing a modified staggered grating and the angular signal radiography, we proposed a single-shot grating-based x-ray differential phase contrast imaging with decent vertical spatial resolution. The theoretical framework was deduced and proved by numerical experiments. Absorption, phase, and scattering computed tomography can be performed without phase stepping. Therefore, we believe this fast and highly stable imaging method with decent resolution would be widely applied in x-ray grating-based phase contrast imaging.     

基于成像差分吸收光谱技术和压缩感知理论的烟囱烟羽断层重建研究

烟羽断层重建质量受两方面条件限制：其中一个限制条件是遥感设备的时间分辨率。以往的研究多使用多轴差分吸收光谱仪（MAX-DOAS）进行CT重建,受采集数据速度的限制,重建图像的时间分辨率较低。另一个限制条件是,采集到的数据量有限,是典型的不完全角度重建。过去多使用代数迭代重建算法或统计迭代重建算法,重建图像受测量误差的影响比较大,分辨率较低且伪影较多。构造了基于成像差分吸收光谱技术（IDOAS）的光谱数据采集系统,与多轴差分吸收光谱仪构造的系统相比,数据采集的时间分辨率提高了160多倍,基本解决了时间分辨率的问题。提出了一种基于压缩感知理论和低三阶导数模型的烟羽断层重建算法--投影凸函数集低三阶导数法,简称为POCS-LTD。在投影的过程中,使用代数重建算法使重建图像符合投影方程;在全变分迭代的过程中使用了优化算法,将低三阶导数模型的全变分归一化值作为优化算法的迭代方向,前次迭代运算结果与本次投影运算的差值的模作为迭代步长。对重建算法进行了数值模拟,并以重建图像的接近度和一致性相关因子为指标,对重建结果进行了分析。数值模拟表明,算法具有良好的抗误差能力,与传统的低三阶导数法相比,本文提出的算法将重建接近度减小了80%以上。使用烟羽数据采集系统进行了外场实验,用POCS-LTD算法对外场实验的数据进行了烟羽重建,重建图像显示烟羽图像清晰,伪影得到了较好的抑制。介绍的烟羽断层数据采集系统和烟羽断层重建算法,提高了烟羽断层重建图像的时间分辨率,减少了重建图像的伪影,扩大了光谱测量技术的应用范围。     

衍射增强计算机断层技术研究

对衍射增强计算机断层技术进行了分析研究，提出了一种新方法.应用该方法只需在摇摆曲线一点处对物体进行360°范围内旋转投影成像，即可获得物体的折射率梯度分布以及吸收系数与散射系数和分布.该方法较以前的方法，大大简化了实验步骤. 关键词： 同步辐射 衍射增强 计算机断层     

高分辨率快速数字化光声CT乳腺肿瘤成像

提出了一种基于聚焦线性阵列探测器的快速光声计算机断层成像技术（光声CT）.在光声二维图像重建中,根据阵列探测器机械扫描和电子扫描相结合的组合扫描模式,提出了改进的有限场滤波反投影重建算法.一方面该算法适合多元探测器旋转扫描模式,另一方面探测器的指向性函数作为反投影的权重因子提高了系统的横向分辨率.同时,该成像系统还利用柱面声透镜实现Z轴方向上的聚焦扫描以实现三维层析成像.实验中,这套成像系统空间分辨率达到0.2mm,Z轴方向分辨率为1.5mm,扫描一幅二维图像仅需150s,得到 关键词： 光声CT 有限场滤波反投影算法 声透镜聚焦 乳腺肿瘤检测     

Projected index computed tomography

Projected index computed tomography (PICT) is a new imaging technique that provides a computed reconstruction of the index of refraction of a sample. PICT makes use of data from standard optical coherence tomography images taken from several view angles to determine a mapping of the refractive indices of the sample. A rectilinear propagation model is assumed, so the data are understood to be related to the line integral of the refractive index in the beam paths. These data thus provide a set of angular projections of the sample. The spatial distribution of the index of the object may then be reconstructed by use of standard filtered backprojection techniques. The resultant PICT images are free of the spatial distortion that is inherent in standard optical cross-sectional images and correspond well to the manufactured dimensions of specific samples.     

代数迭代法用于傅里叶变换计算层析成像光谱仪数据重建

主要介绍代数迭代法在傅里叶变换计算层析成像光谱仪数据处理中的应用。根据4种不同投影数进行了光谱切片重建,并对光谱相似性进行了比较。针对ART算法可应用少数投影重建的特点,专门在45个投影数时,采用不同迭代系数进行了重建,并进行了比较。最后,对实际场景进行了光谱切片重建。     

Feasibility of photoacoustic tomography for ophthalmology

For the eyeball composed of membrane and liquid, the contrast of ultrasound imaging is not high due to its small variance in acoustic impedance. As a new imaging modality, photoacoustic tomography combines the advantages of pure optical and ultrasonic imaging together and can provide high resolution, high contrast images. In this paper, the feasibility of photoacoustic tomography for ophthalmology is studied experimentally. A Q-switched Nd:YAG pulsed laser with 7-ns pulse width is used to generate photoacoustic signal of a porcine eyeball in vitro. The two-dimensional (2D) optical absorption image of the entire eyeball is reconstructed by time-domain spherical back projection algorithm. The imaging results agree well with the histological structure of the eyeball and show a high imaging contrast.     

Optical Computed Tomography Imaging of Absorbers in a Highly Scattering Medium by Detection of Photons Preserving Initial Polarization

A method utilizing light polarization was applied to take the projection data of absorbers in a highly scattering medium. The floor level of the projection data due to diffused light was significantly suppressed by detection of photons which preserved initial polarization. An optical computed tomgraphy image was reconstructed from the projection data taken by this method and exhibited higher spatial resolution than the image reconstructed from the projection data by conventional method.     

针孔单光子发射计算机断层成像的空间分辨率研究

针孔单光子发射计算机断层(SPECT)成像的空间分辨率通常是根据Anger经验公式来进行估算,与实际测量存在较大偏差.本文通过对针孔成像的物理过程进行分析,提出了一个近似度更高的计算公式.利用精确的蒙特卡罗方法模拟针孔SPECT成像,采用OSEM(ordered subsets expectation maximization)算法对投影数据进行图像重建,并与模具实验进行比较,验证了理论公式的适用性.同时还讨论了体素尺寸、几何映射获取投影矩阵以及探测器尺寸与成像物体尺寸比值对断层图像空间分辨率的影响.实验结果显示,该理论公式所估算的空间分辨率比实验值平均偏小约10%,而Anger经验公式所估算的空间分辨率比实验值平均偏大约60%.因此,该理论公式能更好地估算针孔SPECT成像的空间分辨率,可为针孔SPECT系统的设计和使用提供有价值的参考.     

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.     

Elemental x-ray imaging using Zernike phase contras

We develop an element-specific x-ray microscopy method by using Zernike phase contrast imaging near absorption edges, where a real part of refractive index changes abruptly. In this method two phase contrast images are subtracted to obtain the target element: one is at the absorption edge of the target element and the other is near the absorption edge. The x-ray exposure required by this method is expected to be significantly lower than that of conventional absorption-based x-ray elemental imaging methods. Numerical calculations confirm the advantages of this highly efficient imaging method.     

Elemental x-ray imaging using Zernike phase contrast

We develop an element-specific x-ray microscopy method by using Zernike phase contrast imaging near absorption edges, where a real part of refractive index changes abruptly. In this method two phase contrast images are subtracted to obtain the target element: one is at the absorption edge of the target element and the other is near the absorption edge. The x-ray exposure required by this method is expected to be significantly lower than that of conventional absorption-based x-ray elemental imaging methods. Numerical calculations confirm the advantages of this highly efficient imaging method.     

Quantitative in-line phase-contrast imaging with multienergy X rays.

We present a new method for quantitative nondestructive characterization of objects by x-ray phase-contrast imaging. Spatial distributions of the projected values of the complex refractive index in the sample are reconstructed by processing near-field images collected at a fixed sample-to-detector distance using a polychromatic incident beam and an energy-sensitive area detector, such as a CCD used in the photon-counting spectroscopy mode. The method has the potential advantages of decreased radiation dose and increased accuracy compared to conventional techniques of x-ray imaging.     

Polynomial approximation method for tomographic reconstruction of three-dimensional refractive index fields with limited data

We propose a polynomial approximation method (PAM) for reconstruction of three-dimensional refractive index fields by interferometric tomography using limited data. Based on the assumption that the fields to be reconstructed are usually smooth and can be decomposed into a finite order of (orthogonal) polynomials, a set of linear equations can be constructed using both the measured projection data and the Radon transform of the basis functions. By solving these equations, the least-squares solutions of expansion coefficients can be obtained and then substituted back to yield the desired fields. Numerical results have demonstrated that the proposed method is fast, robust to noise and can achieve satisfactory results for refractive index fields with limited projection views and large opaque objects.     

Sheet-beam geometry for in vivo fluorescent x-ray computed tomography: proof-of-concept experiment in molecular imaging

We propose a fluorescent x-ray computed tomography method using an array of detectors with an incident sheet beam, aimed at providing molecular imaging with high sensitivity and good spatial resolution. In this study, we prove the feasibility of this concept and investigate its imaging properties, including spatial and contrast resolutions and quantitativeness, by imaging an acrylic phantom and a normal mouse brain using a preliminary imaging system with monochromatic synchrotron x rays.     

Algorithm Study on Reconstruction of Refractive Angles in Fan Beam Diffraction Enhanced Computed Tomography

Based on the 360° computing method of refractive angle for parallel beam diffraction enhanced imaging computed tomography （DE-CT） technique, a new algorithm used to calculate the refractive angle for fan-beam DE-CT technique is developed. The refractive index gradient can be obtained by using the new algorithm with projection data of an object in the range of 0 - 360°, and the new algorithm only needs to set the analyser at half slope position of the rocking curve.