锥束X射线CT图像重建的新型滤波函数

第13卷第9期2008年9月中国图象图形学报Jour na l of I m ag e and G raphi csV01．13，N o．9S ep．，2008一种基于G r engeat公式的锥束C T重建算法研究张东平张定华张丰收(西北工业大学现代设计与集成制造技术教育部重点实验室，西安710072)摘要基于单圆轨迹的Fel dkam p(FD K)重建算法只有在小锥角条件下才能取得较好重建效果。

随着锥角增大，图像质量迅速下降。

因此在实际应用中受到一定限制。

在分析了物体R adon数据与单圆扫描获得的R adon数据间的关系后，结合空间可变滤波反投影重建算法(SV—FB P)，提出了一种基于G r engeat公式的叠加算法。

该重建算．法由两部分重建结果叠加而成，第1部分结果首先采用FD K算法通过对单圆扫描的投影数据进行重建来获得；然后采用外推方法获得缺失的数据。

并利用SV—FB P进行重建得到第2部分结果；最后将两部分结果进行叠加。

实验结果表明，该算法不仅有效地抑制了FD K算法重建的伪影，而且使锥角的使用范围比FD K算法提高了3—4倍。

这种新的叠加重建算法在大长物体的重建中，具有重要的理论和应用价值。

关键词单圆轨迹FD K算法空间可变滤波反投影重建算法叠加算法中图法分类号：T P391．41文献标识码：A文章编号：1006-8961(2008)09—1649-06A C one—beam C T R e cons t r uct i on A l gor i t hmB ased on G r angea t’S For m ul aZ H A N G D on g—p i ng，Z H A N G D i ng—hua，ZH A N G Feng shou(T h e K ey Labor at ory of C ont em por ar y D e si gn and如t egr at ed M an uf act u r i ng，M i ni s tr y of Ed uc at i on，N o ahw est e m Po l yt e c hn i c al U n i ver s i t y，X i’A n710072)A b st r ac t FD K al gor i t h m ba sed on s i ng l e ci rcular orbi t ca n onl y be use d w hen t he con e angl e i s sm al l．T he i m age sr econs t r uct ed by Fel dk a m p al gor i t hm bec om e bl ur r ed a nd di s tor t ed w hen t he co ne angl e i n cr ea ses．A s a r es ult t her e i s a l i m i t at ion f or pr act i cal appl i cat i on．A f t er t he r el at i onshi p be t w een t he R a d on t r ans f or m of t he object a nd t he R adon dat a acqui r ed by t he sca nni ng of t he s ingle ci r cu l ar orbi t i s anal yzed i n t he pa per a nd t he sup er po si ng al gor i t h m i s pr e sent ed ba sed on G r e ngea t’8f or m ul a acc or di ng t o t he s pace-var i a nt f il te r ed bac kpr oj ec t i on(SV—F B P)．The fi rst s ub-r es ul t i s cal cul at ed by appl yi ng t he FD K al gor i t hm t o t he pr oj ect i on dat a acqu i r ed by a s i ng l e cir cul ar orbi t s ca n．The se con d s ub—r es ult i s a cor r ect i on r ecei ved by appl y i ng SV—F B P al gor i t h m t o t he m i s si ng R adon dat a acqu i r ed by ext r a pol at i on．T he f inal r es ult i s t he s am of t he t w o sub-r esul t s．T he r es ult dem onst r a t es t he al gor i t h m i s a val i d m et hod t o r ed uce a rt if a ct s of t hei m a ges r eco nst r uct ed by FD K al go r i t hm，and t he appl i cabl e con e angl e r ange c an be3t o4t i m e s l ar g er t han t hat of t he FD Kal g or i t hm．T hi s new sup er po si ng m e t ho d w i l l posses s i m p or t antt heo r et i cal a nd appl i cat i on val ue i n t he f ie l d s of bi g a nd l ong obj ect r econst r ucti on．K eyw or ds s i n gl e cir cul ar or bi t，F D K al gor i t hm，s pace-va r i ant f il t er ed backpr oj ect i on al g or i t hm，sup er posi ng al gor i t h m1引言在计算机层析成像(C T)技术13趋成熟的背景下，关于3维物体重建方法的研究取得了突破性的进展，3维图像重建的理论体系也逐步确立起来。

三维锥束重建中滤波器的设计①高 洁1,2,3，杨 波1,2,3，孔 斌1,31(中国科学院合肥智能机械研究所仿生感知与控制研究中心，安徽合肥 230027)2(中国科学技术大学自动化系，安徽合肥 230031)3(安徽省仿生感知与先进机器人技术重点实验室，安徽合肥 230031)摘 要：设计了几种用于CT三维锥束图像重建的滤波器，并将其对图像质量的影响进行了分析。

CT是一种先进的成像技术，现已被应用于多个领域。

CT图像重建算法多采用滤波反投影法，滤波是此方法的核心部分。

本文介绍了三维锥束近似重建的基本原理，分析并设计了几种用于CT锥束图像重建的滤波器，并且通过投影数据的图像重建结果分析了不同滤波函数对断层图像重建质量的影响。

实验结果证明：所设计的几种滤波器能较好的满足实际图像重建的要求。

关键词：CT；图像重建；滤波反投影；滤波函数Design of Filter in 3D Cone Beam ReconstructionGAO Jie1,2,3, YANG Bo1,2,3, KONG Bin1,31(Center for Biomimetic Sensing and Control Research, Institute of Intelligent Machines, Chinese Academy of Sciences, Hefei 230031, China)2(Department of Automation, University of Science and Technology of China, Hefei 230027, China)3(The Key Laboratory of Biomimetic Sensing and Advanced Robot Technology, Hefei 230031, China)Abstract:This article analysises designes for several filters using in the approximate three-dimensional cone beam reconstruction principle, and analysises the effect of CT’s image quality. CT is an advanced imaging technology, and have been used in many fields. The current CT image reconstruction alrorithms mostly use filtered back projection method.Filtering is the core of this method. This article analysises and designes for several filters using in the approximate three-dimensional cone beam reconstruction principle. It analysises the different filter functions affecting the quality of CT image reconstruction according to the results of the projection data’s image reconstruction. These filters can be better to meet the requirements of the actual image reconstruction.Keywords: CT; image reconstruction; filter back projection; filter functionComputerized Tomography简称为CT，即计算机断层成像技术，是由对物体进行不同角度的扫描投影重建获得物体截面信息的成像技术[3]。

锥束X射线CT图像重建的新型滤波函数洪振厚;周彬;郭金川【摘要】The conventional filter function of the FDK algorithm in X-ray computed tomography (CT) can result in the obvious Gibbs phenomenon, which has great effect on the CT image.A NEW-MS-L filter function is deduced based on the rationale of mixed filter theory and weighted average theory.First, the S-L filter function is weighted and averaged to form M3S-L filter function.Then, as a superposition of the M3S-L filter function and the NEW filter function, the NEW-MS-L filter function is pared with the other filter functions such as NEW, R-L-S-L, R-L-NEW, and R-L-MS-L, the NEW-MS-L can suppress the Gibbs phenomenon while maintaining high resolution.The results show that the NEW-MS-L filter function can provide X-ray CT with better reconstruction images.%针对X射线计算机断层成像(computed tomography,CT)图像重建FDK算法中,采用通常的滤波函数会导致明显的Gibbs现象,影响重建图像的质量.基于混合滤波和加权平均理论,设计了一种新型的NEW-MS-L滤波函数.先将S-L 滤波函数加权平均为M3S-L滤波函数,再与NEW滤波函数叠加,得到NEW-MS-L 滤波函数.通过数值仿真,分别对比NEW、R-L-S-L、R-L-NEW、R-L-MS-L和NEW-MS-L滤波函数的重建图像效果,结果表明,NEW-MS-L滤波函数能够在保持较高图像分辨率的情况下,更有效地抑制Gibbs现象,使重建图像效果更佳.【期刊名称】《深圳大学学报（理工版）》【年(卷),期】2017(034)003【总页数】6页(P284-289)【关键词】光学工程;计算机断层成像;FDK算法;滤波函数;Gibbs现象;混合滤波;加权平均【作者】洪振厚;周彬;郭金川【作者单位】深圳大学光电工程学院,光电子器件与系统教育部/广东省重点实验室,广东深圳518060;深圳大学光电工程学院,光电子器件与系统教育部/广东省重点实验室,广东深圳518060;深圳大学光电工程学院,光电子器件与系统教育部/广东省重点实验室,广东深圳518060【正文语种】中文【中图分类】TP391.41随着计算机断层成像(computed tomography, CT)在医疗诊断、工业无损探伤和食品安全检测等领域广泛的应用[1-5]，其对图片质量的要求也越来越高，进而对软硬件特别是算法的要求也越来越高．在众多算法中，FDK算法[6]应用最为普遍．FDK算法是由Feldkamp、Davis和Kress提出的一种基于圆轨道扫描的近似重建算法，其本质是滤波反投影(filtered back projection，FBP)，它对图像质量的影响非常明显．传统的滤波函数能保持较高的空间分辨率，但同时引起明显的Gibbs现象(即有明显的振荡效应)，致使重建效果较差．因此，在保持高图像分辨率的情况下设计新的滤波函数，降低Gibbs现象就尤为重要．为此，研究人员相继提出NEW滤波函数[7]、R-L-S-L、R-L-NEW和R-L-MS-L等混合滤波函数[8-10]来消除Gibbs现象．虽然这些滤波函数都具有保持高空间分辨率的同时减小Gibbs现象的作用，但无法完全消除Gibbs现象，主要原因是与之结合的R-L滤波函数的近旁瓣突出，远旁瓣的幅度和宽度较大[9-10]，而NEW滤波函数和R-L滤波函数一样，由理想滤波器推导出，所以其近旁瓣也较为突出[7].本研究基于FDK算法，设计了一种新型滤波函数进行图像重建，在保持较高空间分辨率的情况下大幅减小Gibbs现象，提高了密度分辨率，进而改进了重建图像质量．FDK算法因其高效、简便以及易于图像处理器(graphics processing unit，GPU)加速，至今仍被大量使用．FDK算法本质上是FBP，把得到的锥束投影数据进行滤波，然后利用扇形束近似重建而得．FDK算法的计算公式[11]如式(1)和式(2)．对投影数据P的加权滤波为其中， R为光源到物体中心的距离；β为源绕中心旋转轴z轴的旋转角度； a为旋转中心的横坐标； b为旋转中心的纵坐标；P(β, a, b)为投影数据； G(a)为滤波函数；符号*表示卷积运算．反投影重建结果为β其中，U(x, y, β)=R+xcos β+ysin β为加权函数，这里x和y为探测器坐标；P′(β, a, b)为加权滤波后的投影数据；其他变量定义如式(1)．影响重建结果好坏的关键在于滤波函数的选择、重建过程中插值的方式，以及锥角的大小等．从式(2)可直观地看出，滤波函数能直接影响到反投影重建的结果．所以，合适的滤波器是实现高质量重建图像的关键因素之一．2.1 选择要混合和加权平均的滤波函数要选择合适的滤波函数进行混合和加权平均，必须要先判断滤波函数的优劣特性．判断一个滤波函数对重建结果的影响主要有：主瓣、近旁瓣和远旁瓣．主瓣高而窄，说明空间分辨率好；远旁瓣的幅度和幅值越小，说明Gibbs现象越小、密度分辨率越好[12-13]．范惠荣等[14]研究表明，NEW滤波函数能够保持空间分辨率的同时减小Gibbs现象，而S-L滤波函数[11]可通过这3点加权平均使远旁瓣的幅度和幅值大幅减小，其近旁瓣收敛相比R-L[11]、S-L和NEW滤波函数更快，更能有效地抑制Gibbs现象，但因其主瓣变矮，空间分辨率会变得很差，如图1．其中， n为采样点； h[n]为采样点n所对应的滤波函数值．图2中除R-L滤波函数外，剩下3种滤波函数的远旁瓣几乎重叠，表明M3S-L滤波函数[10]的远旁瓣的幅度和幅值非常小，与R-L滤波函数相比无明显振荡，说明Gibbs现象非常小．因此，可将NEW滤波函数和加权平均后的M3S-L滤波函数混合叠加，得到新的NEW-MS-L滤波函数，大幅减小Gibbs现象，且能保持与单独使用NEW 滤波函数后相当的空间分辨性能．2.2 构建新型滤波函数NEW-MS-L新型滤波函数是根据加权平均和混合滤波的思想构建的．首先需将S-L滤波函数进行加权平均．S-L滤波函数的离散形式[11]为hS-L[n]=其中，n=0,±1,±2,…, 为采样点(后面各式含义相同)；τ为探测器单元的大小，一般设为1．基于文献[15]的研究结果，本研究对S-L滤波函数进行3点加权平均，记为M3S-L，该函数的离散形式[10]为hM3S-L[n]= 0.6hS-L[n]+0.2hS-L[n-1]+0.2hS-L[n+1]将M3S-L滤波函数和NEW滤波函数进行混合滤波，得到新的滤波函数NEW-MS-L．NEW滤波函数[7]和NEW-MS-L滤波函数的离散形式分别为hNEW[n]=-hNEW-MS-L[n]= k1×hM3S-L [n]+k2×hNEW[n]其中， k1和k2为权重， k1+k2=1．本研究取k1=0.5； k2=0.5．NEW-MS-L、NEW和M3SL滤波函数在空域的主瓣分布对比如图3．从图3可见，新型滤波函数的hNEW-MS-L[n]主瓣高度与NEW滤波函数的相当，近旁瓣幅度比NEW滤波函数的小，且收敛更快．远旁瓣的幅度和幅值很小，表明NEW-MS-L滤波函数可实现与NEW滤波函数相当的空间分辨率，且能更有效抑制Gibbs现象，进一步提高密度分辨率，获得更优的图像质量．基于FDK算法，分别用NEW、R-L-S-L、R-L-N、R-L-MS-L和NEW-MS-L滤波函数对尺寸为256×256×256像素的Shepp-Logan三维头模型[16]进行重建，然后观察中心面处的截面及其纵向第128行的灰度曲线图，结果如图4．图5(a)至图5(f)分别为采用NEW、R-L-S-L、R-L-NEW、NEW-MS-L和R-L-MS-L滤波函数重建的纵向第128行的灰度曲线图．灰度曲线图的波动大小反应与原始模型的差异，波动越大说明差异大，表现为图像越粗糙．由图5可见，NEW-MS-L较其他滤波函数能够更好的还原原始模型．仅通过灰度曲线图的对比来说明NEW-MS-L滤波函数的优越性是不够的，还需通过计算NEW、R-L-S-L、R-L-N、R-L-MS-L和NEW-MS-L滤波函数的归一化均方误差d和归一化平均绝对误差r，如式(7)和式(8)，进行定量对比.r=其中， i和j分别为图像的横坐标和纵坐标； N为横纵坐标的点数； fijrec为重建结果的灰度值； fij为原始模型重建结果的灰度值；为原始模型重建结果的灰度值的平均值．d反映少数点的大误差， d=0表示重建后图像完全再现了原始模型图像； d越小表示两者的偏差越小． r反映多数点的小误差， r=0表示重建图像与头模型原始图像无误差； r越小说明误差越小．表1为分别采用NEW、R-L-S-L、R-L-NEW、R-L-MS-L和NEW-MS-L滤波函数重建图像后的d值和r值对比．由表1可见，采用NEW-MS-L滤波函数后的d值与采用NEW滤波函数后的d值相当，且所得r值较其他滤波函数的都小，说明采用NEW-MS-L滤波函数所得的重建效果更佳．上述结果还可通过调整k1和k2值来进一步优化[10]．d和r随k1的变化如表2，由表2可知，当k1=0.8时， d值最小；但当k1=0.7、 0.8和0.9时， d值虽然相差无几，但r值随k1的增加而变大，且变化较明显．所以，本研究权衡d和r的值，选取k1=0.7．当k1=0.7， k2=0.3时，可得到d=0.315 5， r=0.737 3．虽然r值增大了，但d值进一步减小．将NEW-MS-L滤波函数的d值和r值与表1的其他滤波函数的d值和r值进行对比，结果反映新滤波函数性能比其他滤波函数要好．基于混合滤波和加权平均理论，构建新型滤波函数NEW-MS-L．其空域主瓣分布图具有主瓣高而窄、近旁瓣小、远旁瓣衰减快等特征，能很好地抑制Gibbs现象，保持较高的空间分辨率．仿真结果表明，通过对比NEW、R-L-S-L、R-L-NEW、R-L-MS-L滤波函数重建结果与灰度曲线图，并分析归一化均方误差和归一化平均绝对误差，发现NEW-MS-L滤波函数可在保持较高空间分辨率的情况下更有效地抑制Gibbs现象，密度分辨率更高，重建结果也更平滑．此外，根据实际情况适当调整k1和k2，能使图像效果达到符合预期的要求．NEW-MS-L滤波函数在本次仿真实验中主要针对传统基于吸收效应的CT，今后可进一步探讨NEW-MS-L滤波函数在相衬CT[17-22]中的应用．【相关文献】[1] Egan C K,Jacques S D M,Wilson M D, et al.3D chemical imaging in the laboratory by hyperspectral X-ray computed tomography[J].Scientific Reports,2015,5: 15979.[2] Bieberle M,Barthel bined phase distribution and particle velocity measurementin spout fluidized beds by ultrafast X-ray computed tomography[J].Chemical Engineering Journal,2016,285:218-227.[3] Edlund R,Skorpil M,Lapidus G,et al.Cone-beam CT in diagnosis of scaphoid fractures[J].Skeletal Radiology,2016,45(2):197-204.[4] Senyshyn A,Mühlbauer M J,Dolotko O,et al.Homo geneity of lithium distribution in cylinder-type Li-ion batteries[J].Scientific Reports,2015,5:18380.[5] Plessis A D, Olawuyi B J, Boshoff W P, et al. Simple and fast porosity analysis of concrete using X-ray computed tomography[J]. Materials and Structures, 2016, 49(1): 553-562.[6] Feldkamp L A, Davis L C, Kress J W. Practical cone-beam algorithm[J]. Journal of the Optical Society of America A, 1984, 1(6): 612-619.[7] Wei Yuchuan, Wang Ge, Hsieh J. An intuitive discussion on the ideal ramp filter in computed tomography (I)[J]. Computers & Mathematics with Applications, 2005, 49(5/6): 731-740.[8] 谢强．计算机断层成像技术:原理、设计、伪像和进展[M]．张朝宗，郭志平，王宗贤，译．北京：科学出版社，2006． Hsieh puted tomography: principle, design, artifacts and recent advances[M]. Zhang Chaozong,Guo Zhiping, Wang Zongxian, trans. Beijing: Science Press,2006.(in Chinese)[9] 胡永梅，全红，徐利明，等．新型混合滤波器在FDK算法中的应用[J]．中国医学影像技术，2011，27(2)：397-400． Hu Yongmei, Quan Hong, Xu Liming, et al. Application of a new mixed filter in FDK algorithm [J]. Chinese Journal of Medical Imaging Technology, 2011, 27(2): 397-400.(in Chinese)[10] 王晓鹏，王明泉，侯慧玲．基于R-L-MS-L滤波函数的CT图像重建[J]．电视技术，2014，38(7)：26-28． Wang Xiaopeng,Wang Mingquan,Hou Huilin. CT image reconstruction based on R-L-MS-L filter function[J].Video Engineering, 2014,38(7): 26-28.(in Chinese) [11] 闫镔，李磊．CT图像重建算法[M]．北京：科学出版社，2014． Yan Bin,Li Lei.CT image reconstructinalgorithm[M].Beijing:Science Press,2014.(in Chinese)[12] 郑君里．信号与系统[M]．北京：高等教育出版社，2002． Zhen Junli. Signals and systems[M].Beijing:Higher Education Press,2002.(in Chinese)[13] 杨位钦，顾岚．时间序列分析与动态数据建模[M]．北京：北京工业学院出版社，1986． Yang Weiqin,Gu Lan. Time series analysis and dynamic data modeling[M].Beijing:Beijing Institute of Tchnology Press,1986.(in Chinese)[14] 范惠荣，徐茂林，邱钧，等．关于CBP算法的一种新型滤波函数和它的性质[J]．电子学报，2004，32(2)：232-235． Fan Huirong,Xu Maolin,Qiu Jun,et al. A new filter function in CBP and its properties[J].Acta Electronica Sinica,2004,32(2):232-235.(in Chinese)[15] 刘晓，杨朝文．RL滤波函数的改进对卷积反投影图像重建的影响[J]．四川大学学报自然科学版， 2004，41(1)：112-117． Liu Xiao,Yang Chaowen. The influence of modified-RL filterfunction on the quality of image reconstruction using convolution backprojection algorithm [J].Journal of Sichuan University Natural Science Edition,2004,41(1):112-117.(in Chinese)[16] 刘泽，孙丰荣，李艳玲，等．基于3D Shepp-Logan头部模型的三维医学图像重建仿真[J]．生物医学工程学杂志，2006，23(5)：938-943． Liu Ze,Sun Fengrong,Li Yanling, etal.The three-dimension medical image reconstruction simulation on 3D Shepp-Logan head phantom [J].Journal of Biomedical Engineering,2006,23(5):938-943.(in Chinese) [17] Lei Yaohu,Huang Jianheng,Liu Xin,et al.Improvement of visibility of moir fringe in X-ray differential phase-contrast imaging[J].Journal of Shenzhen University Science and Engineering,2016, 33(5):506-510.[18] 刘鑫，郭金川．部分相干光源微分相衬双图像相位恢复[J]．深圳大学学报理工版，2014，31(2)：169-173． Liu Xin,Guo Jinchuan.Dual-images phase retrieval in differential phase-contrast imaging with partial coherence source[J].Journal of Shenzhen University Science and Engineering,2014,31(2):169-173.(in Chinese)[19] Jerjen I, Revol V, Kottler C, et al. Phase contrast cone beam tomography with an X-ray grating interferometer [C]// International Conference on Advanced Phase Measurement Methods in Optics & Imaging.Ascona, Switzerland: American Institute of Physics, 2010, 1236(1):227-231.[20] Huang Zhifeng, Kang Kejun, Li Zheng, et al. Direct computed tomographic reconstruction for directional-derivative projections of computed tomography of diffraction enhanced imaging[J]. Applied Physics Letters, 2006, 89(4): 041124.[21] McCann M T, Nilchian M, Stampanoni M, et al. Fast 3D reconstruction method for differential phase contrast X-ray CT[J]. Optics Express, 2016, 24(13): 14564-14581. [22] Fu Jian, Hu Xinhua, Velroyen A, et al. 3D algebraic iterative reconstruction for cone-beam X-ray differential phase-contrast computed tomography[J]. PloS One, 2015, 10(3): e0117502.。

第16卷第4期 CT理论与应用研究V ol. 16 No. 42007年12月（31-37） CTTheoryandApplications Dec.,2007 文章编号：1004-4140（2007）04-0031-07锥束CT图像重建算法的快速实现吴胜利，潘瑞谊，文斌（西安交通大学 电子物理与器件教育部重点实验室，西安710049）摘要：本文基于锥束CT滤波反投影重建的FDK算法，通过两种算法改进并结合基于共享内存的Open MP并行技术和代码优化，实现了锥束CT图像的快速重建。

基于锥束CT实际投影数据的重建结果表明，图像重建速度得到了较大的提高，断层图像重建质量与FDK原型算法相当。

关键词：锥束CT；图像重建；FDK滤波反投影法；Open MP并行技术中图分类号：TP391.41 文献标识码：A锥束CT利用锥状X光束和面阵探测器[1－3]可以一次扫描获得被检物体全周投影数据，具有扫描速度快、分辨率高、重建图像质量好等优点，是CT领域中活跃的研究课题之一。

锥束CT实用化的关键在于提高图像重建速度和图像质量，这方面主要取决于图像重建算法的进步，因此，重建算法的研究一直都是CT研究非常重要的一个方面。

锥束CT图像重建算法分为解析法和迭代法两大类[4]。

迭代法直接重建三维图像，是一种高分辨率重建法，但重建处理速度慢[5]，需要更多的计算资源，实际中为了保证高速重建，一般采用解析法。

解析法分为近似算法和精确算法，近似算法由于数学上简单、容易实现，且当锥角较小时能够取得较好的重建效果,因而在实际中得到广泛的应用[6]。

目前应用最多最广的算法是FDK滤波反投影重建法，其后由该算法改进而发展出的一些算法如P-FDK、T-FDK[7]、S-FDK[8]也都属于近似算法，锥束CT精确重建算法发展于2002年以后[9－10]。

对FDK算法改进并结合数据并行处理可以快速重建三维图像[11]，本文通过对FDK滤波反投影法的算法改进并基于具有双核处理器的计算机上的代码优化，在保证重建图像质量的同时提高了重建速度、节约了资源。

CT图像重建中基于指数形式的滤波函数优化马思汉;张催;陈章谷;张益海;潘小东;商宏杰;李公平【期刊名称】《无损检测》【年(卷),期】2017(039)012【摘要】计算机断层成像(CT)技术广泛应用于工业、医疗、国防、航空航天等领域.随着科学技术的不断发展,对使用CT来获得更加精确的材料内部结构信息提出了更高的要求,希望能在重建图像中更加完整地表达组织结构信息.为了获取更好的图像质量,提出一种基于指数形式的滤波函数(EBFF),在MATLAB平台下采用EBFF 对图像进行重建,并通过相似系数、归一化平均值以及归一化均方根距离等评价手段与几种传统的滤波函数重建图像结果进行对比,同时在锥束CT实验平台下采用EBFF进行图像重建.模拟与试验均表明EBFF对重建图像质量有一定改善,其较传统滤波函数能给出更好的重建结果,有助于图像重建结果的优化.【总页数】6页(P1-6)【作者】马思汉;张催;陈章谷;张益海;潘小东;商宏杰;李公平【作者单位】兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000;兰州大学核科学与技术学院,兰州 730000【正文语种】中文【中图分类】TG115.28;TL99【相关文献】1.基于能谱滤波和图像残差重投影的CT图像金属伪影校正方法 [J], 李磊;韩玉;席晓琦;王敬雨;闫镔2.粒子滤波算法在ECT图像重建中的应用 [J], 吴新杰;黄国兴;王静文3.CT图像重建滤波反投影算法中指数滤波器的研究 [J], 骆岩红4.Kalman滤波在ECT图像重建中的应用 [J], 霍启正;5.CT图像重建中的一种新型滤波器 [J], 杨坪坪;冯汉升;许继伟;宋云涛因版权原因，仅展示原文概要，查看原文内容请购买。

重排的半覆盖螺旋锥束CT的反投影滤波重建邹晓兵;曾理【期刊名称】《光学精密工程》【年(卷),期】2010(018)009【摘要】由于传统的螺旋锥束CT系统采用的扫描方法使视场范围受到限制,本文深入分析了螺旋锥束CT的半覆盖扫描方式,提出了螺旋锥束半覆盖扫描的重建算法.该扫描方式只要求射线束覆盖物体横截面的一半以上,探测器的宽度和投影数据量大约是传统螺旋锥束CT的一半,从而节约了探测器的成本.结合扇形束投影的特点和二维反投影滤波(BPF)来实现这种横向截断投影的重建,研究了重排的加权BPF重建方法.该方法首先根据重建的层将半覆盖的螺旋锥束CT投影数据重排为半覆盖的扇形束CT投影,然后用加权的扇形束BPF算法进行重建.同传统视场区域横截面全覆盖的FDK算法相比,该方法即使在大螺距的情况下其重建结果也更接近于真实值,而且均方误差减少了一半,重建时间节约了3/5.【总页数】9页(P2077-2085)【作者】邹晓兵;曾理【作者单位】重庆大学,光电技术及系统教育部重点实验室,ICT研究中心,重庆,400044;重庆教育学院,计算机科学系,重庆,400067;重庆大学,光电技术及系统教育部重点实验室,ICT研究中心,重庆,400044;重庆大学,数理学院,重庆,400044【正文语种】中文【中图分类】TP391;R814.42【相关文献】1.改进的大锥角锥束CT圆轨迹反投影滤波算法 [J], 李增光;韩玉;李磊;闫镔2.非PI线的半覆盖螺旋锥束CT的Katsevich重建 [J], 邹晓兵;李宗剑3.基于混合滤波的避投影数据微分螺旋锥束CT重建 [J], 马建华;陈凌剑;颜刚;陈武凡4.螺旋锥束计算机断层成像倾斜扇束反投影\r滤波局部重建算法 [J], 席晓琦;韩玉;李磊;闫镔5.一种用于小体积偏置探测器锥束CT系统的反投影滤波重建算法 [J], 李亮; 陈志强; 张丽; 邢宇翔; 康克军因版权原因，仅展示原文概要，查看原文内容请购买。

螺旋锥束CT重建的近似逆算法胡红莉;张建州【摘要】三维螺旋锥束CT以扫描速度快、成像分辨率高等诸多优点成为现代CT 技术的一个重要发展方向.Katsevich精确FBP算法的提出,使得三维锥束CT研究获得了突破性进展.由于该算法的复杂性,应用中受到了限制.研究了Katsevich算法在检测板上沿滤波线展开的形式,其滤波运算由Hilbert核函数构成,利用近似逆的思想提出了融合的CT重建算法.该算法将Katsevich 公式改写成近似逆的形式,得到了重建核的具体形式.%The 3-dimensional helical cone-beam CT has many advantages of fast scanning, high resolution imaging, and becomes an important development direction of modern CT technology.Katsevich proposes exact FBP algorithm, which promotes access to a breakthrough in 3-dimensional cone-beam CT studies.With the complexity of the algorithm, the application has been limitedAn integration CT reconstruction algorithm based on Katsevich algorithm and approximate inverse algorithm is proposed.The filtering operation is constituted by the Hilbert kernel function.The Katsevich algorithm is rewritten into the form of approximate inverse,and the specific forms of reconstructed kernel are given.【期刊名称】《计算机工程与应用》【年(卷),期】2011(047)021【总页数】3页(P199-201)【关键词】近似逆;重建核;锥束螺旋CT重建;Katsevich算法【作者】胡红莉;张建州【作者单位】四川大学计算机学院,成都610051;四川大学计算机学院,成都610051【正文语种】中文【中图分类】TP1811 引言近年来采用滑环技术和先进的重建技术使计算机层析成像技术（CT）克服了许多实际问题，并在此基础上设计出了螺旋CT。