ART算法中几种加权因子模型的研究.docVIP

ART算法中几种加权因子模型的研究 　　文章编号： 1005??5630(2014)02??0142??05 　　收稿日期： 2013??11??28 　　摘要： ART(algebraic reconstruction technique)算法是一种适合于投影数据采集量比较少的情况的图像重建算法。利用其进行图像重建时的主要工作是计算加权因子,该计算方法严重影响图像重建的重建质量和重建速度。讨论、研究并仿真了加权因子的三种计算模型,经比较分析最后得出了一种最优的重建模型。 　　关键词： 图像重建; 代数重建算法; 加权因子 　　中图分类号： TP 391文献标志码： Adoi： 10.3969/j.issn.1005??5630.2014.02.011 　　Research on several weight factor models of ART 　　CHEN Liang 　　(Navy Representative Office of Jiangnan Shipyad Co., Ltd., Shanghai 201913, China) 　　Abstract： ART (algebraic reconstruction technique) algorithm is fit for data acquisition from incomplete projections. The main work is to solve of weight factor in ART algorithm and the computing method has a strong impact on the quality and speed of image reconstruction. Finally, by discussing, studying and simulating these three computing methods, an optimal reconstruction model is obtained through comparative analysis. 　　Key words： image reconstruction; algebraic reconstruction technique(ART); weight factor 　　引言图像重建算法主要分为解析法和迭代法两大类。解析法中传统的滤波反投影(FBP)［1］算法仍然被广泛使用,它可以在投影数据采集完备时重建出较为准确的图像,但当投影数据采集不完备时,重建出的图像会出现严重的伪影。这时一般采用迭代法来解决,迭代重建法的经典方法是Gorden等提出的ART(algebraic reconstruction technique)算法［2］,该算法可以在投影数据比较少的情况下重建出精确的图像,其原理是将重建问题转换成用迭代方法求解线性方程组的问题。ART算法的特点是:适合于不完备投影数据的图像重建,重建质量好,图像的密度分辨率高,同时对空间分辨率也有很好的保证;其次,重建算法计算简单,不同形式的采样数据重建都适用;另外还可以结合一些已知的先验知识得到更准确的解。由于在ART算法的迭代过程中,每次投影都要计算修正值,但这些修正值并不是完全相同的,穿过同一个像素网格时,图像的误差修正将会给重建区域带来严重的噪声,且要想得到较好的重建效果,算法需要进行多次迭代,从而导致计算量大,重建时间长,因而人们一直关注的是如何提高该算法的重建效率。ART算法中花费时间最长的是加权因子的计算,本文首先研究了ART算法的3种加权因子计算模型,然后通过理论和仿真实验对用3种模型来重建图像的速度和质量进行比较分析,得出了实用的结论。图1Radon变换示意图 　　Fig.1Diagram of Radon transform1原理 　　1.1Radon变换在一个平面内对函数f(x,y)沿不同的直线做线积分,就得到了函数f(x,y)的Radon变换［3］。也就是说,平面上每个点的像函数值都有一个原始函数的线积分值与之对应,如图1所示。光学仪器第36卷 　　第2期陈亮：ART算法中几种加权因子模型的研究 　　图中每条射线的方程用xcosθ+ysinθ=t来表示,夹角为θ的射线的投影值可以表示为Pθ(t)=∫∞－∞∫f(x,y)δ(xcosθ+ysinθ－t)dxdy(1)图像重建的任务就是由各个方向不同位置的投影值P(t)恢复出原图像的分布函数f(x,y)。 　　1.2ART算法原理重建图像的方法有很多种,本文着重介绍和研究ART算法。将待重建的区域离散成一个N=n×n的正方形小网格［4］,其中每个小网格内的值fj是指所在的第j个像素的平均值