针对X射线锥束成像模拟计算量大、速度慢的问题，提出了一种基于图形处理器(GPU)的快速成像模拟算法。该算法沿着每条射线累加所经过体素对投影值的贡献量，实现了对X射线成像的模拟。在计算射线与体素的交线长时，采用分类处理交线的方法，减少了增量Siddon算法的动态分支计算。为了提高投影图像质量，该算法还用GPU硬件线性插值采样取代Siddon算法的邻近插值采样。对三维Shepp-logan模型的测试结果表明，该算法的速度比基于GPU的增量Siddon算法平均提高了44%，而且图像质量明显提高。最后，用实测数据进一步验证了算法的有效性。

所属栏目

科研成果与学术交流国家自然科学基金资助项目(60972140，60971131);北京市自然科学基金资助项目(3102009);北京市教委科技发展计划资助项目(KM201010028001)

收稿日期

2011/1/18

作者单位

杨涛:首都师范大学 检测成像工程中心，北京100048
赵星:首都师范大学 检测成像工程中心，北京100048

备注

杨涛(1982-)，男，博士研究生，主要从事计算机图形学、科学可视化等方面的研究。

引用该论文:

YANG Tao,ZHAO Xing.A GPU-Based Algorithm for the Simulation of X-Ray Cone-Beam Imaging[J].Nondestructive Testing,2011,33(12):6～13
杨涛,赵星.基于图形处理器的X射线锥束成像模拟算法[J].无损检测,2011,33(12):6～13

参考文献

【1】

Duvauchelle P, Freud N, Kaftandjian V, et al. A computer code to simulate X-ray imaging techniques[J]. Nuclear Instruments and Methods in Physics Research Section B,2000,170(1):245-258.

【2】

Freud N, Duvauchelle P, Letang J M, et al. Fast and robust ray casting algorithms for virtual X-ray imaging[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms,2006,248(1):175-180.

【3】

邵军明，路宏年，蔡慧.X射线成像的一种点扩展函数模型[J].光学学报,2005,25(8):1148-1152.

【4】

孔凡琴，路宏年.数字X射线图像散射的自适应抑制算法[J].光电工程,2005,32(6):86-88.

【5】

邹晶，孙艳勤，张朋.由少量投影数据快速重建图像的迭代算法[J].光学学报,2009,29 (5):1198-2004.

【6】

赵星，胡晶晶，潘晓川，等.一种新的基于GPU实现的锥束CT正投影算法[J].电子学报,2009,37(6):1165-1169.

【7】

Xu F. Fast implementation of iterative reconstruction with exact ray-driven projector on GPUs[J]. Tsinghua Science and Technology,2010,15(1):30-35.

【8】

Siddon R L. Fast calculation of the exact radiological path for a three dimensional CT Array[J]. Medical Physics,1985,12(2):252-255.

【9】

Jacobs F, Sundermann E, Sutter B D, et al. A fast algorithm to calculate the exact radiological path through a pixel or voxel space[J]. Journal of Computing and Information Technology,1998,6(1):89-94.

【10】

Han G, Liang Z, You J. A fast ray tracing technique for TCT and ECT studies[J]. IEEE Transactions on Medical Imaging,2000,19(5):556-564.

【11】

Zhao H, Reader A J. Fast ray-tracing technique to calculate line integral paths in voxel arrays[C]. // Proceedings of IEEE Nuclear Science Symposium. Portland, USA,2003:211-218.

【12】

Greef M, Crezee J, Eijk J C, et al. Accelerated ray tracing for radiotherapy dose calculations on a GPU[J]. Medical Physics,2009,36(9):4095-4102.