基于压缩传感和代数重建法的CT图像重建

代数重建法(ART)是一种重要的CT图像重建方法,适合于不完全投影数据的图像重建,其缺点是重建速度慢。为提高图像重建的质量和速度,利用压缩传感理论提出了一种基于ART的高质量图像重建算法。该算法将CT图像的梯度稀疏性结合到ART图像重建中,在每次迭代中的投影操作结束后用梯度下降法调整全变差,减小图像梯度的l1范数。实验结果验证了该算法的有效性。     

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

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

A computed tomography reconstruction algorithm based on multipurpose optimal criterion and simulated annealing theory

Although emission spectral tomography (EST) combines emission spectral measurement with optical computed tomography (OCT), it is difficult to gain transient emission data from a large number of views,therefore, high precision OCT algorithms with few views ought to be studied for EST application. To improve the reconstruction precision in the case of few views, a new computed tomography reconstruction algorithm based on multipurpose optimal criterion and simulated annealing theory (multi-criterion simulated annealing reconstruction technique, MCSART) is proposed. This algorithm can suffice criterion of least squares, criterion of most uniformity, and criterion of most smoothness synchronously. We can get global optimal solution by MCSART algorithm with simulated annealing theory. The simulating experiment result shows that this algorithm is superior to the traditional algorithms under various noises.     

Photoacoustic image reconstruction based on Bayesian compressive sensing algorithm

The photoacoustic tomography (PAT) method, based on compressive sensing (CS) theory, requires that, for the CS reconstruction, the desired image should have a sparse representation in a known transform domain. However, the sparsity of photoacoustic signals is destroyed because noises always exist. Therefore, the original sparse signal cannot be effectively recovered using the general reconstruction algorithm. In this study, Bayesian compressive sensing (BCS) is employed to obtain highly sparse representations of photoacoustic images based on a set of noisy CS measurements. Results of simulation demonstrate that the BCS-reconstructed image can achieve superior performance than other state-of-the-art CS-reconstruction algorithms.     

抗混叠轮廓波域采用压缩感知的云图融合方法

提出一种采用压缩感知的云图融合方法.该方法针对传统轮廓波存在频谱混叠的缺点,结合抗混叠塔式滤波器组和方向滤波器组,构造出一种抗混叠的轮廓波变换,并将其引入压缩感知中的稀疏表示环节,将云图分解成稠密和稀疏两部分|对稠密成份采用传统方法进行融合,而对稀疏成份,则在压缩感知框架下,通过少数线性测量的融合,并采用二步迭代收缩的图像重构算法,在迭代时利用前面两个估计值更新当前值,得到融合结果.实验表明,该方法的融合结果无论在视觉质量及定量指标上都明显优于传统方法,有利于揭示全面的天气信息.     

基于稀疏优化的计算机断层成像图像不完全角度重建综述

计算机断层成像(computed tomography,CT)技术在医学和工业无损检测中都具有非常广泛的应用,CT重建算法是其中的核心,而不完全角度重建问题则是实际应用中重建算法研究领域的一个热点和难点问题.近年来,随着稀疏优化理论与算法的飞速发展,基于稀疏优化的重建算法已经在不完全角度重建问题中得到了较广泛的应用,且表现出了良好的精度与速度性能.本文首先对稀疏优化的基本理论结论与常用算法进行了介绍;而后对稀疏优化理论在CT图像不完全角度重建中的应用进行归纳,分类介绍了其主要研究成果及稀疏优化所发挥的作用;最后对基于稀疏优化的不完全角度重建研究进行了展望.     

基于三投影方向的层析重建分析

针对传统的层析方法中,投影方向多、系统复杂等特点,提出了基于三个投影方向的层析重建方法.论文选取适用于少投影方向的代数迭代重建算法,采用参量修正方法,能较好的提高重建准确度和抑制椒盐噪音.同时对影响重建图像质量和速度的迭代次数、不同原图像结构对重建图像质量的影响进行了分析,给出了相应分析模拟结果.结果表明,迭代300次后,重建结果基本稳定.同时对于具有轴对称结构物体,三个方向投影数据能够较好地实现三维层析重建.分析结果为后期开展的基于单幅层析全息图的三维层析重建实验研究提供了理论依据.     

基于聚类和小波变换的多光谱图像压缩算法

针对多光谱图像压缩算法现存的时空复杂度高、光谱特性利用不充分等问题,研究了多光谱图像的谱间稀疏等价表示及其聚类实现途径,进而设计了一种基于谱间自适应聚类和小波变换的多光谱图像压缩算法。算法利用吸引力传播聚类产生多光谱图像的谱间稀疏等价表示、在低复杂度下去除图像的谱间冗余,使用二维小波变换去除稀疏表示成分的空间冗余,采用分层树集合分割排序算法(SPIHT)进行压缩编码,并通过误差补偿机制提高多光谱图像重建质量。实验表明,该算法在保证较低时间和空间复杂度的基础上,较SPIHT等同类经典压缩算法,在相同的压缩比下,明显提高了重建图像的峰值信噪比,是一种通用有效的多光谱图像压缩算法。     

总变差约束的数据分离最小图像重建模型及其Chambolle-Pock求解算法

基于优化的迭代法,可以结合压缩感知和低秩矩阵等稀疏优化技术高精度地重建图像.其中,总变差最小（total variation minimization,TV）模型是一种简单有效的优化模型.传统的约束TV模型,使用数据保真项为约束项,TV正则项为目标函数.本文研究TV约束的、数据分离最小（TV constrained,data divergence minimization,TVcDM）新型TV模型及其求解算法.详细推导了TVcDM模型的Chambolle-Pock（CP）算法,验证了模型及算法的正确性;分析了算法的收敛行为;评估了模型的稀疏重建能力;分析了模型参数的选择对重建的影响及算法参数对收敛速率的影响.研究表明,TVcDM模型有高精度稀疏重建能力;TVcDM-CP算法确保收敛,但迭代过程中有振荡现象;TV限对重建有重要影响,参数值过大会引入噪声而过小会模糊图像细节;算法参数的不同选取会导致不同的收敛速率.     

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.     

基于L1范数的微分相位衬度CT稀疏角度重建算法

微分相位衬度成像及其计算层析(CT)技术是近年出现的无损检测新方法。但是,相位衬度CT往往需要对样品进行多次扫描,这必将导致非常长的辐射时间和巨大的辐射剂量。稀疏角度重建在降低辐射剂量方面有着非常明显的优势,因此,研究针对相位衬度CT的稀疏角度重建算法就显得尤为重要。在分析了相位衬度CT的特点之后,将压缩感知理论引入相位衬度CT重建中,并在该理论框架下将L1约束融入代数迭代重建(ART)算法中,提出了一种微分相位衬度CT重建算法。数值模拟和实际实验表明,该方法可以根据少量投影数据给出较好的重建结果。     

代数迭代重建算法在折射衬度CT中的应用

X射线折射衬度CT是一种基于相位衬度的断层成像技术,特别适合对由轻元素组成的生物、医学样品进行成像,可以观察到常规吸收衬度CT无法观察到的软组织内部微细结构,是一种具有巨大发展潜力的新成像方法.迭代重建算法和解析重建算法是计算机断层成像技术中并行发展的两种算法,虽然已经提出了几种X射线折射衬度CT的解析重建算法,可是还未见X射线折射衬度CT迭代重建算法的报道.研究了代数迭代重建算法在X射线折射衬度CT中的应用,比较分析了不同的投影数据排列方式对于折射衬度CT重建图像的影响,并对实验数据进行了图像重建,获得了满意的CT图像.研究结果表明,在相位衬度CT中,迭代重建算法相对于解析重建算法而言,能减少投影次数,降低曝光剂量,减少对生物样品的辐射损伤,在生物样品成像和投影数据不完整的情况下具有明显的优势. 关键词： 衍射增强成像 代数迭代算法 CT重建 同步辐射     

稀疏阈值的超分辨率图像重建

为了解决基于字典学习的超分辨重构算法耗时过长的问题,提出了基于稀疏阈值模型的图像超分辨率重建方法。首先,将联合字典理论与图像块稀疏阈值方法相结合,训练得到高、低分辨率过完备图像字典对。接着,通过稀疏阈值OMP算法对图像特征块进行稀疏表示。然后,通过高分辨率字典重构出初始的超分辨图像。最后,通过改进迭代反投影算法对初始的超分辨图像进行全局优化,从而进一步提高图像重构质量。实验结果表明,超分辨图像重构平均峰值信噪比(PSNR)为30.1 d B,平均结构自相似度(SSIM)为0.937 9,平均计算时间为10.2 s。有效提高了超分辨重构的速度,改善了重构高分辨图像的质量。     

Imaging resolution analysis in limited-view Laser Radar reflective tomography

In many applications of Laser radar reflective tomography, anisotropic resolution of images reconstructed from an incomplete set of reflective projections will degrade imaging quality. In this paper, we present the theoretical point spread function (PSF) of imaging system over full views for reflective tomography. The corresponding imaging resolution of the PSF was determined by Rayleigh's criterion, here the most common algorithm of filtered backprojection was used for image reconstruction. The theoretical imaging resolution was derived approximately from the range resolution resolved by the laser pulse. The simulated horizontal imaging resolutions of an ideal diffused single point target in different limited-views were compared with the theoretical value. Experimental reconstructed images from limit-view projections of letter “E” illustrated the effects of anisotropic resolution.     

Regularized maximum likelihood algorithm for PET image reconstruction using a detail and edges preserving anisotropic diffusion

The traditional iterative reconstruction algorithms of positron emission tomography cannot effectively suppress the noise in low SNR case. Recently anisotropic diffusion (AD) is introduced into tomography reconstruction, which can improve the reconstructed image. Although AD reconstruction algorithm can suppress noise, it does not perverse the detail edge information accurately, especially the thin edges. In order to solve the problem, we introduce a new anisotropic diffusion term, which can preserve the detail edges effectively, into the maximum likelihood algorithm, and combine with median filter, forming the regularized maximum likelihood algorithm in PET image reconstruction (PML_NewAD). Results of computer simulated demonstrate that compared with the other classical reconstruction algorithms, PML_NewAD not only availably suppress the noise and produce a higher quality image, but also preserve the structure of image's edge excellently.     

Convex relaxation based sparse algorithm for hyperspectral target detection

Target detection in hyperspectral images is an important task. In this paper, we propose a sparsity based algorithm for target detection in hyperspectral images. In sparsity model, each hyperspectral pixel is represented by a linear combination of a few samples from an overcomplete dictionary, and the weighted vector for such reconstruction is sparse. This model has been applied in hyperspectral target detection and solved with several greedy algorithms. As conventional greedy algorithms may be trapped into a local optimum, we consider an alternative way to regularize the model and find a more accurate solution to the model. The proposed method is based on convex relaxation technique. The original sparse representation problem is regularized with a properly designed weighted ?₁ minimization and effectively solved with existing solver. The experiments on synthetic and real hyperspectral data suggest that the proposed algorithm outperforms the classical sparsity-based detection algorithms, such as Simultaneous Orthogonal Matching Pursuit (SOMP) and Simultaneous Subspace Pursuit (SSP) and conventional ?₁ minimization.     

Convergence rate calculation of simultaneous iterative reconstruction technique algorithm for diffuse optical tomography image reconstruction: A feasibility study

Diffuse optical tomography (DOT) is an emerging technique for functional biological imaging. The imaging quality of DOT depends on the imaging reconstruction algorithm. The SIRT has been widely used for DOT image reconstruction but there is no criterion to truncate based on any kind of residual parameter. The iteration loops will always be decided by experimental rule. This work presents the CR calculation that can be great help for SIRT optimization. In this paper, four inhomogeneities with various shapes of absorption distributions are simulated as imaging targets. The images are reconstructed and analyzed based on the simultaneous iterative reconstruction technique (SIRT) method. For optimization between time consumption and imaging accuracy in reconstruction process, the numbers of iteration loop needed to be optimized with a criterion in algorithm, that is, the root mean square error (RMSE) should be minimized in limited iterations. For clinical applications of DOT, the RMSE cannot be obtained because the measured targets are unknown. Thus, the correlations between the RMSE and the convergence rate (CR) in SIRT algorithm are analyzed in this paper. From the simulation results, the parameter CR reveals the related RMSE value of reconstructed images. The CR calculation offers an optimized criterion of iteration process in SIRT algorithm for DOT imaging. Based on the result, the SIRT can be modified with CR calculation for self-optimization. CR reveals an indicator of SIRT image reconstruction in clinical DOT measurement. Based on the comparison result between RMSE and CR, a threshold value of CR (CR_T) can offer an optimized number of iteration steps for DOT image reconstruction. This paper shows the feasibility study by utilizing CR criterion for SIRT in simulation and the clinical application of DOT measurement relies on further investigation.     

Tomographic reconstruction of damage images in hollow cylinders using Lamb waves

Lamb wave tomography (LWT) is a potential and efficient technique for non-destructive tomographic reconstruction of damage images in structural components or materials. A two-stage inverse algorithm proposed by the authors for quickly reconstructing the damage images was applied to hollow cylinders. An aluminum hollow cylinder with an internal surface pit and a Carbon Fiber Reinforced Plastic (CFRP) laminated hollow cylinder with an artificial internal surface damage were used to validate the proposed method. The results show that the present method is capable of successfully reconstructing the images of the above damages in a larger inspection area with much less experimental data compared to some conventional ultrasonic tomography techniques.     

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.