第 42 卷 第 5 期 电 子 科 技 大 学 学 报 Vol.42 No.5
2013年9月 Journal of University of Electronic Science and Technology of China Sep. 2013
圆弧加直线轨迹Katsevich类型的锥束CT重建算法
王 瑜
1
,陈 亮
1
,欧宗瑛
2
,董彦磊
1
,黄建龙
1
(1. 电子科技大学机械电子工程学院 成都 611731; 2. 大连理工大学精密与特种加工教育部重点实验室 辽宁 大连 116024)
【摘要】由于圆与直线轨迹的交点处存在轨迹变量的偏导数,Katsevich算法直接忽略了关于构造因子的分析,从而对后
续重建计算产生影响,导致重建结果出现伪像。通过对该问题进行研究,确定轨迹交点处的构造因子,并推导出在平板检测
器下与两段轨迹相对应的重建公式。据此引申出一个不含有关于轨迹变量偏导的等价公式,从而计算出轨迹交点处相应的重
建贡献值。同时,提出使用光滑过渡函数改善因反投影区域边界的不连续性而对重建结果产生的影响。利用Forbild头模型进
行仿真模拟,实验结果表明,改进的算法有效地抑制了重建结果的伪像。
关 键 词 计算机断层成像; 锥束重建; 滤波反投影; Katsevich重建算法
中图分类号 TP391.41; TG115.28 文献标志码 A doi:10.3969/j.issn.1001-0548.2013.05.030
Cone-Beam CT Reconstruction of Katsevich-Type Algorithm for
Orthogonal Arc-and-Line Trajectory
WANG Yu
1
, CHEN Liang
1
, OU Zong-ying
2
, DONG Yan-lei
1
, and HUANG Jian-long
1
(1. School of Mechatronics Engineering, University of Electronic Science and Technology of China Chengdu 611731;
2. Key Laboratory for Precision and Non-Traditional Machining Technology of Ministry of Education,
Dalian University of Technology Dalian Liaoning 116024)
Abstract Because Katsevich algorithm includes the partial derivative with respect to the variable of
trajectory, the analysis of structure factor and subsequent computation of reconstruction are ignored at the
intersection point of two source orbits consisting of a circle and orthogonal lines, which may induce artifacts in the
reconstructed image. This paper investigates the Katsevich algorithm for this geometry, determines the structure
factor at the intersection point, and deduces the corresponding reconstruction formula about two source trajectories
using the planar detector. As a result, an equivalence formula is achieved excluding the partial derivative with
respect to the variable of trajectory and can be used to calculate the contribution of critical plane at the intersection
point. In addition, a smooth transitional function is suggested to improve the discontinuity of contribution at the
border of projection region defined by PI-lines. Computer simulations about Forbild phantom show a reduction of
the artifacts that are found with the Katsevich algorithm and verify the validity of the proposed algorithm.
Key words computed tomography, cone-beam reconstruction, filtered backprojection, Katsevich cone-
beam
收稿日期:
2011 12 25; 修回日期:2012 09 02
基金项目:
辽宁省高等学校科学研究青年基金(2004F057);中央高校基本科研业务费专项资金(ZYGX2012J097)
作者简介:王瑜(1972 ),女,博士,主要从事CT重建理论方面的研究.
文献[1]提出具有滤波反投影结构的锥束精确重
建算法,取得CT重建理论领域内突破性的进展。由
于Katsevich重建公式是一种通用表达式,其具体执
行细节紧密地依赖于扫描轨迹与检测器的几何形状
[2]
,为此许多学者展开了大量的研究
[2-11]
。在医学领
域内,螺旋CT是最常见的断层成像设备,针对螺旋
轨迹的重建算法研究的较多
[2-6]
。近年来,因平板检
测器的成像优势
[12-13]
,利用血管照影C臂设备产生清
晰的断层图象引起了关注
[14-16]
。该设备主要可生成
单圆弧、圆弧加圆弧或直线轨迹。单圆弧轨迹不满
足锥束精确重建的充要条件,为获得精确的重建结
果,需要依靠圆弧加轨迹
[7-11]
。相对而言,直线轨迹
更容易采集单圆弧轨迹缺失的锥束投影数据。本文
对基于平板检测器对圆加直线轨迹上的Katsevich重
建算法
[10-11]
进行了研究;分析了轨迹交点位置处存
在 构 造 因 子 不 为 零 的 Radon 平 面 通 过 推 导 出
Katsevich重建公式的具体表达,解决了Katsevich算
法无法准确计算轨迹交点处关于扫描轨迹变量的偏
导数的问题。另一方面,依据与重建点相关联的PI
线段,在某些投影视角上,投影数据仅在检测器平