基于螺旋锥束CT的两种解析重建算法的研究.pdfVIP

研讨会·论文集 全陶射线数字成像与CT新技术 中国·上海 基于螺旋锥束CT的两种解析重建算法的研究 王勃，潘晋孝，孔慧华 (中北大学电子测试技术国家重点实验室，太原030051) 改进，利Jfj插值多项式函数的导数来作为所求导数的近似值，然后比较了这两种重建算法。实 验结果表明：当投影数据没有噪声的时候，FDK算法和Katsevich算法均能取得较好的效果，当 对实际物体进行重建的时候，由J：投影数据含有噪卢，而Katsevich算法需要对投影数据求导， 对投影数据的噪声较敏感，重建质鬣有所下降。 关键词：螺旋cT；图像重建；FDK算法；Katsevich算法 近年来，螺旋锥束cT在医学和工业上得到越来越广泛的应用，同时螺旋锥束图像重建算 法也在迅速发展。螺旋CT重建算法可分为：近似重建算法、精确畦建算法等。近似重建算 法以FDK算法为代表，其实现简单、成像速度快，主要缺点是X射线z轴向锥角较大时远离 有突破性进展。该算法虽然是在理论上是精确的，但现实中要实现精确却很困难，Yu Hengyong等‘对实现该算法有比较好的结果，并对实现重建中的伪影进行了分类¨1。 Katsevich重建算法在重建过程中需要对投影数据求导，但投影数据本身是离散的，不能直 接求导，通常对离散数据求导采用差商法来作为所求导数的近似值，本文考虑利用插值多 项式函数的导数来作为所求导数的近似值，通过改变对投影数据的求导方法来改进 算法，并对这两种算法进行了比较研究，最后通过仿真数据和实际采集的实验数据对这两 种重建算法进行比较。 1 算法介绍 螺旋锥束CT几何位置犬系如图l所示 物体 平板探测器 y 图l螺旋锥束cT几何位置关系 下勃等：摹于螺旋锥束CT的两种解析重建算法的研究 7l 射线源S的扫描轨迹为(Rcoss，Rsins，兰)，R为射线源到物体中心的距离，即旋转半 21t 径，h为螺距，D为射线源到探测器的距离，J为旋转角度，尸为物体上的任意一点，在s 角度下，P点在探测器上的投影的位置为Q，其在140’v嫩标系下的班标为(甜，1，)，凶此投 影书可表示为g(s，U，v)。 1．1螺旋FDK算法 螺旋FDK算法哺1是从圆轨迹FDK算法推广而来的，其重建公式如F： 舳力一』鲁∽训灿 其帕加降一鸭睾+刁 U=R+xsins-ycoss (2) ～．1， 材：2：!兰!竺!兰±上!!翌兰2：2：!兰!竺!兰±羔兰!翌兰2 (3) U R+xsinS—ycoss 】V：：一=一 v：土篓!：尘兰 (4) U R+xsinS—ycoss 喜cs，甜，V，=(：志·gcs，甜，V，]}，，c“， (5) ^(”)为斜坡滤波函数，幸表示卷积符号。 扫描视角内根据式(1)进行累加重加。 1．2螺旋Katsevich重建算法及其改进 锥束重建理论的+‘项突破性进展，根据螺旋轨迹Katsevich重建算法的丛本原理，可将其分 为对投影数据关于旋转角度求导