几种CT图像重建算法的研究和比较.pdfVIP

第 1O卷 第 4期 2001年 I1月(4～9) CT理论与应用研究 CT Theory and Applications 、bI．10．NO 4 Nov一20o1 几种CT图像重建算法的研究和比1殴★ 张 朋 张兆田 北京100101：2)国家 北京 10OO85 摘要 本姗 究了滤波反投影算法cF日P]、局部重建算法(Localn)、修正的代数重建算法 (黼 )．用特别设计的 仿鹁蝴断 高 螺 ，分别对这三种募出韭行了成像数值试验。分析和比较了各算法的重建速度、空间分辨率、密度分 辨率，及其优缺点和适用范围。 关建最 计算机层析或I拳 重建算法 重建速度 分辨率 碍sN文章 号{1∞“14o(2001)O4-0OO4436 中圈分类号l TP391 41 如辨际硼玛：I A Research and Comparison oil Several CT Reconstruction Algorithms ZHANG Peng ZHANGZhaoti~ 2 1) Di,i~mofA 醐 ]vlafiael／l~cs,B卿_岫窟 _帆T 帆Rd1【l0Io hl_衄呜Ch 珏B~ijing100101 1) 删 oflafcrmaaonscien~NationalNatlaz]ScienceFotmda6on。f0咄 glO0。B5 ABSTRACT：Filter back projection O：BP)aI local recznmuctkm algorithm 4Dcal RA)and modifed algebraic reconsm~fion technique(MARI3 a researched．The 0 e e“矗 of ime~ g re。0ns ∞ are made to these algorithms by using desi ed s~nulzfion dam 龃d real XCT dam respectively． Analysis and comparison ofthese aI酬 thms On recons虮 on sPe吐 spatial resolution,densityresolution,as wel as em：lvantagesmaddisadv．,mtagesare m Keywords：Comptm~tomof*d~ ,Recomtn~on algorithm，Reconsm~on Sp∞d， 1a resolution． 1引 言 可分为： 精确算法、近1l算法和迭代算法。近似算法中，以滤波反投 影算法 (Filter back projection，abbr．FBP)最具代表性，应用最为广泛。选代算法中，代数重建 算法 (Algebmic reconstruction technique，abbr．ART)是提出最早井最为人们熟悉的算法。选 代型算法 (如代数重建算法等)具有许多优点，但由 柏事量大、重建时间长．在很长一段时间内限 制了其在医学和工业 CT领域的应用。 们提出了不少提高迭代算崩 十算速度的方法，加 E近年 — 直是 们?关注的问题。近年来人 受到人们青睐。此外，由于应用的需要，局部重建算法( a王Reconstruction hlgorith~abbr． LocalP~)也在近十年中有了较大的发展。在传统全局 cr算法中．即使重建物体断面中—个小区域的 图像，也得围绕整个断面采集投影数据。而局部重建算法， 据，即可重建感兴趣区域的图像。局部重建算法 体)的放射摄入量。 队 体(或生物 近年来，找们先后研究了滤波反投影算法、代数 彭 局部重毫算法，提出了—种基于块 选代的修正的代数重建算法 (MART)，并又多次优化改进，该算法较传统 A盯 算法重建速度太幅度提 高；设计了分别针对 FBP和LocalRA算法的滤波器：研究了LocalRA影像补 撤正 。采用—些典型的 仿真数据和实采数据，对 FBP、MA肝、LocalRA算法 亍了重建试验，就重建速度、重建图像质量、 收稿日期；2001-09—28．本研究得到国家自然科学基金和教育部高等学校骨干教师资助计划资助 4期 垄塑董! !里堡重堡蔓堕塑塑塞塑生墼 了比较、分析，指出了各算法的优缺点，以及适用范囝。本文将报告 t述 研究结果。以下各节安排如下：第2节简介FBP、 和LocMRA算法，第3节报告数值试验结果， 第 4 2算法介绍 设厂( )为物体被检测断面的线性衰减系数的分布，则 ，1( )为R 上的函数，满足 ／( )=0，当 =( f，X2)诺Q时。记0 =(OOS ，si力 ) ，0 =(一sin~,cos~) ．R 中 直线三可表示为： ·0=P，