一种圆轨迹锥束CT中截断投影数据的高效重建算法

汪先超; 闫镔; 刘宏奎; 李磊; 魏星; 胡国恩

doi:10.7498/aps.62.098702

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名

邮箱

手机号码

标题

留言内容

验证码

摘要

本文基于数据重排方法, 提出了T-BPF (Tent-BPF)算法, 该算法先将锥束投影数据重排成平行投影数据, 然后使用一种推导的BPF型算法重建重排后的平行投影数据. T-BPF算法将原BPF算法反投影中变化的角度积分限变成固定的, 反投影中各层循环之间没有了相关性, 这意味着T-BPF算法较原BPF算法具有更好的可并行性. 实验结果显示: 使用GPU对2563的Shepp-Logan体模的图像重建进行并行加速, T-BPF算法在保证重建质量的前提下, 加速比达到了1036, 较原BPF算法有很大提升. T-BPF算法为截断投影数据的3D图像快速重建提供了方法.

关键词:

X射线光学 /
CT /
图像重建 /
GPU

作者及机构信息

1.
国家数字交换系统工程技术研究中心, 郑州 450002;

2.
大连理工大学建设工程学院, 大连 116000

基金项目: 国家高技术研究发展计划(批准号:2012AA011603)资助的课题.

参考文献

[1]	Zou Y, Pan X C 2004 Phys. Med. Biol. 49 941
[2]	Zou Y, Pan X C 2004 Phys. Med. Biol. 49 383
[3]	Zou Y, Pan X C 2004 Phys, Med. Biol. 49 2717
[4]	Clackdoyle R, Noo F, Guo J, Roberts J 2004 IEEE Trans Nucl Sci 51 2570
[5]	Luo Z Y, Yang X Q, Meng Y Z, Deng Y 2010 Acta Phys. Sin. 59 8237 (in China) [罗召洋, 杨孝全, 孟远征, 邓勇 2010 物理学报 59 8237]
[6]	Jia P X, Zhang F, Yan B, Bao S L 2010 Chin. Phys. B 19 087802
[7]	Li L, Kang K J, Chen Z Q, Zhang L, Xing Y X 2008 IEEE Nuclear Science Symposium and Medical Imaging Conference Paper No M10-346
[8]	Wang X C, Yan B, Li L, Hu G E 2012 Chin. Phys. B 21 118702
[9]	Pan X C, Xia D, Zou Y, Yu L F 2004 Phys. Med. Biol. 49 4349
[10]	Yu L F, Zou Y, Sidky E Y, Pelizzari C A, Munro P, Pan X C 2006 IEEE Trans Med Imag 25 869
[11]	Yu L F, Xia D, Zou Y, Pan X C, Pelizzari C A, Munro P 2005 IEEE Nuclear Science Symposium Conference Record 412
[12]	Yu L F, Pelizzari C A, Pan X C 2006 IEEE Nuclear Science Symposium Conference Record 2889
[13]	Wang X C, Li L, Yu C Q, Yan B, Bao S L 2012 J. of X-ray Sci. and Tech. 20 69
[14]	Grass M, Köhler Th, Proksa R 2000 Phys. Med. Biol. 45 329
[15]	Han Y, Yan B, Li L, Yu C Q, Li J X, Bao S L 2012 Chin. Phys. B 21 068701
[16]	Noo F, Clackdoyle R, Pack J D 2004 Phys. Med. Biol. 49 3903
[17]	Shen L, Xing Y X 2010 Nuclear Techniques 33 857(in Chinese) [沈乐, 邢宇翔 2010 核技术 33 857]
[18]	Wang Y J, Hu H F, Xing Y X 2009 The 10th international meeting on fully 3D image reconstruction in radiology and nuclear medicine 57
[19]	Sidky E Y, Zou Y, Pan X C 2005 Phys. Med. Biol. 50 1643
[20]	Zou Y, Pan X C, Sidky E Y 2005 J. Opt. Soc. Am. 22 2372
[21]	Tricomi F C 1951 Quarterly Journal of Mathematics 2 199
[22]	Sidky E Y, Pan X C 2005 IEEE Signal Process Lett. 12 97
[23]	Feldkamp L A, Davis L C, Kress J W 1984 J. Opt. Soc. Amer. A 1 612
[24]	Noo F, Defrise M, Clackdoyle R, Kudo H 2002 Phys. Med. Biol. 47 2525
[25]	Dennerlein F, Noo F, Hornegger J, Lauritsch G 2007 Phys. Med. Biol. 52 3227
[26]	Shepp L A, Logan B F1974 IEEE Trans Nucl. Sci. 21 21
[27]	Han X, Bian J G, Eaker D R, Kline T L, Sidky E Y, Ritman E L, Pan X C 2011 IEEE Trans Med. Imag. 30 606

施引文献

[1]	Zou Y, Pan X C 2004 Phys. Med. Biol. 49 941
[2]	Zou Y, Pan X C 2004 Phys. Med. Biol. 49 383
[3]	Zou Y, Pan X C 2004 Phys, Med. Biol. 49 2717
[4]	Clackdoyle R, Noo F, Guo J, Roberts J 2004 IEEE Trans Nucl Sci 51 2570
[5]	Luo Z Y, Yang X Q, Meng Y Z, Deng Y 2010 Acta Phys. Sin. 59 8237 (in China) [罗召洋, 杨孝全, 孟远征, 邓勇 2010 物理学报 59 8237]
[6]	Jia P X, Zhang F, Yan B, Bao S L 2010 Chin. Phys. B 19 087802
[7]	Li L, Kang K J, Chen Z Q, Zhang L, Xing Y X 2008 IEEE Nuclear Science Symposium and Medical Imaging Conference Paper No M10-346
[8]	Wang X C, Yan B, Li L, Hu G E 2012 Chin. Phys. B 21 118702
[9]	Pan X C, Xia D, Zou Y, Yu L F 2004 Phys. Med. Biol. 49 4349
[10]	Yu L F, Zou Y, Sidky E Y, Pelizzari C A, Munro P, Pan X C 2006 IEEE Trans Med Imag 25 869
[11]	Yu L F, Xia D, Zou Y, Pan X C, Pelizzari C A, Munro P 2005 IEEE Nuclear Science Symposium Conference Record 412
[12]	Yu L F, Pelizzari C A, Pan X C 2006 IEEE Nuclear Science Symposium Conference Record 2889
[13]	Wang X C, Li L, Yu C Q, Yan B, Bao S L 2012 J. of X-ray Sci. and Tech. 20 69
[14]	Grass M, Köhler Th, Proksa R 2000 Phys. Med. Biol. 45 329
[15]	Han Y, Yan B, Li L, Yu C Q, Li J X, Bao S L 2012 Chin. Phys. B 21 068701
[16]	Noo F, Clackdoyle R, Pack J D 2004 Phys. Med. Biol. 49 3903
[17]	Shen L, Xing Y X 2010 Nuclear Techniques 33 857(in Chinese) [沈乐, 邢宇翔 2010 核技术 33 857]
[18]	Wang Y J, Hu H F, Xing Y X 2009 The 10th international meeting on fully 3D image reconstruction in radiology and nuclear medicine 57
[19]	Sidky E Y, Zou Y, Pan X C 2005 Phys. Med. Biol. 50 1643
[20]	Zou Y, Pan X C, Sidky E Y 2005 J. Opt. Soc. Am. 22 2372
[21]	Tricomi F C 1951 Quarterly Journal of Mathematics 2 199
[22]	Sidky E Y, Pan X C 2005 IEEE Signal Process Lett. 12 97
[23]	Feldkamp L A, Davis L C, Kress J W 1984 J. Opt. Soc. Amer. A 1 612
[24]	Noo F, Defrise M, Clackdoyle R, Kudo H 2002 Phys. Med. Biol. 47 2525
[25]	Dennerlein F, Noo F, Hornegger J, Lauritsch G 2007 Phys. Med. Biol. 52 3227
[26]	Shepp L A, Logan B F1974 IEEE Trans Nucl. Sci. 21 21
[27]	Han X, Bian J G, Eaker D R, Kline T L, Sidky E Y, Ritman E L, Pan X C 2011 IEEE Trans Med. Imag. 30 606

引用本文:

Citation:

计量

文章访问数: 1529
PDF下载量: 434
被引次数: 0

一种圆轨迹锥束CT中截断投影数据的高效重建算法

1. 国家数字交换系统工程技术研究中心, 郑州 450002;

2. 大连理工大学建设工程学院, 大连 116000

基金项目:

国家高技术研究发展计划(批准号:2012AA011603)资助的课题.

收稿日期: 2012-10-15

修回日期: 2013-01-04

刊出日期: 2013-05-05

摘要: 本文基于数据重排方法, 提出了T-BPF (Tent-BPF)算法, 该算法先将锥束投影数据重排成平行投影数据, 然后使用一种推导的BPF型算法重建重排后的平行投影数据. T-BPF算法将原BPF算法反投影中变化的角度积分限变成固定的, 反投影中各层循环之间没有了相关性, 这意味着T-BPF算法较原BPF算法具有更好的可并行性. 实验结果显示: 使用GPU对2563的Shepp-Logan体模的图像重建进行并行加速, T-BPF算法在保证重建质量的前提下, 加速比达到了1036, 较原BPF算法有很大提升. T-BPF算法为截断投影数据的3D图像快速重建提供了方法.

关键词:

X射线光学 /
CT /
图像重建 /
GPU

Efficient reconstruction from truncated data in circular cone-beam CT

1.
National Digital Switching System Engineering and Technological Research Center, Zhengzhou 450002, China;

2.
Faculty of Infrastructure Engineering, Dalian University of Technology, Dalian 116000, China

Fund Project:

Project supported by the National High Technology Research and Development Program of China (Grant No. 2012AA011603).

Received Date: 15 October 2012

Accepted Date: 04 January 2013

Published Online: 05 May 2013

Abstract: In circular cone-beam computed tomography (CT), to solve the 3D image reconstruction from truncated projection data which has no truncation along PI-line, backprojection-filtration (BPF) algorithm is a preferred choice. However, in its performance the integral interval of backprojection is variable for different PI-line, rendering the parallelism performance of backprojection low. So it cannot satisfy the requirement of fast image reconstruction in practical CT system. In this paper, a tent BPF (T-BPF) algorithm is developed based on the data rebinning method, which was performed by first rearranging the cone-beam data to tent-like parallel-beam data, and then applying the proposed BPF-type algorithm to reconstruct images from the rearranged data. T-BPF turns the variable view-angle integral interval of backprojection into a fixed integral interval, and there are no relations in the loops of backprojection calculation, which means the parallelism performance of T-BPF is an improvement over that of the original BPF algorithm. The results of experiments show that compared with the conventional CPU implementation, the GPU accelerated method provides images of the same quality with a speedup factor 1036 for the reconstruction of 2563 Shepp-Logan model. The speedup factor is an improvement in the original BPF algorithm. T-BPF provides a solution for the 3D fast reconstruction from truncated data.

Keywords:

全文HTML

参考文献 (27)

搜索

留言板