搜索

x

留言板

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

姓名
邮箱
手机号码
标题
留言内容
验证码

半覆盖锥束CT中扁平物体的高效反投影滤波重建

张峰 闫镔 汪先超 江桦 魏星

引用本文:
Citation:

半覆盖锥束CT中扁平物体的高效反投影滤波重建

张峰, 闫镔, 汪先超, 江桦, 魏星

An efficient reconstruction algorithm for flat object based on backprojection filtration method in circular half-cover computed tomography

Zhang Feng, Yan Bin, Wang Xian-Chao, Jiang Hua, Wei Xing
PDF
导出引用
  • 全覆盖圆轨迹扫描的成像视野受探测器宽度限制, 对于大物体的成像效率较低.半覆盖扫描可以将成像视野扩展近1倍, 图像重建首推使用反投影滤波型算法. 反投影滤波型算法按PI线重建, 各PI线积分区间的不一致性导致通信和计算消耗大, 影响重建效率. 针对半覆盖成像中扁平形状物体的重建问题, 提出了一种改进的反投影滤波型算法, 且证明了当扁平物体的厚度小于2Rsin(2π/Np) (R为扫描半径, Np 为圆扫描一周均匀采集的投影数量)时, PI线积分区间的不一致性在数值计算过程中的误差是可以忽略的. 改进后的算法相比原半覆盖反投影滤波算法具有两个明显的优势: 一是数值计算过程中角度循环移至PI线循环之外, 算法的通信需求显著降低; 二是投影数据求导、反投影和沿PI线滤波三个步骤均能够并行计算, 算法的并行性得到增强.数值仿真与实际数据的实验结果表明, 本文算法与原半覆盖反投影滤波算法的重建精度相当, 但计算效率提高了4.6倍.
    In circular full-cover cone-beam computed tomography (CT), the field-of-view (FOV) is limited by the width of planar detector, resulting in low imaging efficiency for large object. The FOV can be doubled by half-cover scanning, in which the back-projection filtration (BPF) algorithm based on the concept of PI-line is the best choice for image reconstruction. However, the integral intervals of different PI-lines are unequal in the BPF algorithm, leading to heavy communication consuming and calculation. As a result, the reconstruction efficiency by use of the BPF algorithm is low. In this paper, an efficient image reconstruction strategy based on the BPF algorithm for flat object is proposed. With the method, we demonstrate that the inequality of integral interval of PI-line can be ignored in the discrete implementation of the BPF algorithm when the thickness of flat object is less than 2Rsin(2π/Np) (R is the scanning radius and Np is the number of uniform sampled projections in a full circle). Compared with the original BPF algorithm for half-cover scanning, our method has two major advantages: the first one is that the outer loop is the sample angle while the inner loop is the PI-line, which reduces the communication consuming for computer significantly; the second one is that the derivative of projection, back-projection and inverse Hilbert transform along the PI-line can be computed using parallel computing techniques readily. The results of numerical simulation and real data experiment indicate that the computational efficiency of the proposed method is 5.6 times that for original BPF algorithm and the reconstruction errors of the two methods are comparable.
    • 基金项目: 国家高技术研究发展计划(批准号: 2012AA011603)资助的课题.
    • Funds: Project supported by the National High Technology Research and Development Program of China (Grant No. 2012AA011603).
    [1]

    Katsevich A 2003 Int. J. Math. 21 1305

    [2]

    Zou Y, Pan X, Sidky E Y 2005 J. Opt. Soc. Am. A 22 2732

    [3]

    Yu H Y, Ye Y B, Wang G 2008 J. X-ray Sci. Technol. 16 243

    [4]

    Zeng L, Liu B D, Liu L H, Xiang C B 2010 J. X-ray Sci. Technol. 18 266

    [5]

    Wang X C, Yan B, Li L, Hu G E 2012 Chin. Phys. B 21 118702

    [6]

    Wang X C, Hu G E, Yan B, Han Y, Li L, Bao S L 2013 IEEE Trans. Nucl. Sci. 60 174

    [7]

    Xu Q, Yu H Y, Mou X Q, Zhang L 2012 IEEE Trans. Med. Imaging 31 1682

    [8]

    Zou X B, Zeng L, Li Z J 2009 J. X-ray Sci. Technol.17 233

    [9]

    Han Y, Yan B, Li L, Yu C Q, Li J X, Bao S L 2012 Chin. Phys. B 21 068701

    [10]

    Jia P X, Zhang F, Yan B, Bao S L 2010 Chin. Phys. B 19 087802

    [11]

    Wang G 2002 Med. Phys. 29 1634

    [12]

    Zou Y, Pan X C, Sidky E Y 2005 Phys. Med. Biol. 50 13

    [13]

    Li L, Chen Z Q, Zhang L, Xing Y X, Kang K J 2007 Appl. Radiat. Isot. 65 1041

    [14]

    Guo J Q, Zeng L, Zou X B 2011 J. X-ray Sci. Technol. 19 293

    [15]

    Danielsson P E, Edholm P, Seger M 1997 International Meeting On Fully Three-dimensional Image Reconstruction in Radiology and Nuclear Medicine Pennsylvania, USA, June 25-28, 1997 p141

    [16]

    Wang X C, Li L, Yu C Q, Yan B, Bao S L 2012 J. X-ray Sci. Technol. 20 69

    [17]

    Pan X C, Xia D, Zou Y Yu L F 2004 Phys. Med. Biol. 49 4349

    [18]

    Yu L F, Zou Y, Sidky E Y, Pelizzari C A, Munro P, Pan X C 2006 IEEE Trans. Med. Imaging 25 869

    [19]

    Tuy H K 1983 SIAM J. Appl. Math. 43 546

    [20]

    Kuen H T, Wayne L 2010 Proceedings of the 18th Annual ACM/SIGDA International Symposium on Field Programmable Gate Arrays New York, USA, February 20-27, 2010 p115

    [21]

    Zhang F, Lu L Z, Li Q L, Yan B, Li L 2012 2nd World Congress on Computer Science and Information Engineering Changchun, China, June 17-19, 2012 p687

    [22]

    Peter T 1996 Ph. D. Dissertation (Lyngby: Technical University of Denmark)

  • [1]

    Katsevich A 2003 Int. J. Math. 21 1305

    [2]

    Zou Y, Pan X, Sidky E Y 2005 J. Opt. Soc. Am. A 22 2732

    [3]

    Yu H Y, Ye Y B, Wang G 2008 J. X-ray Sci. Technol. 16 243

    [4]

    Zeng L, Liu B D, Liu L H, Xiang C B 2010 J. X-ray Sci. Technol. 18 266

    [5]

    Wang X C, Yan B, Li L, Hu G E 2012 Chin. Phys. B 21 118702

    [6]

    Wang X C, Hu G E, Yan B, Han Y, Li L, Bao S L 2013 IEEE Trans. Nucl. Sci. 60 174

    [7]

    Xu Q, Yu H Y, Mou X Q, Zhang L 2012 IEEE Trans. Med. Imaging 31 1682

    [8]

    Zou X B, Zeng L, Li Z J 2009 J. X-ray Sci. Technol.17 233

    [9]

    Han Y, Yan B, Li L, Yu C Q, Li J X, Bao S L 2012 Chin. Phys. B 21 068701

    [10]

    Jia P X, Zhang F, Yan B, Bao S L 2010 Chin. Phys. B 19 087802

    [11]

    Wang G 2002 Med. Phys. 29 1634

    [12]

    Zou Y, Pan X C, Sidky E Y 2005 Phys. Med. Biol. 50 13

    [13]

    Li L, Chen Z Q, Zhang L, Xing Y X, Kang K J 2007 Appl. Radiat. Isot. 65 1041

    [14]

    Guo J Q, Zeng L, Zou X B 2011 J. X-ray Sci. Technol. 19 293

    [15]

    Danielsson P E, Edholm P, Seger M 1997 International Meeting On Fully Three-dimensional Image Reconstruction in Radiology and Nuclear Medicine Pennsylvania, USA, June 25-28, 1997 p141

    [16]

    Wang X C, Li L, Yu C Q, Yan B, Bao S L 2012 J. X-ray Sci. Technol. 20 69

    [17]

    Pan X C, Xia D, Zou Y Yu L F 2004 Phys. Med. Biol. 49 4349

    [18]

    Yu L F, Zou Y, Sidky E Y, Pelizzari C A, Munro P, Pan X C 2006 IEEE Trans. Med. Imaging 25 869

    [19]

    Tuy H K 1983 SIAM J. Appl. Math. 43 546

    [20]

    Kuen H T, Wayne L 2010 Proceedings of the 18th Annual ACM/SIGDA International Symposium on Field Programmable Gate Arrays New York, USA, February 20-27, 2010 p115

    [21]

    Zhang F, Lu L Z, Li Q L, Yan B, Li L 2012 2nd World Congress on Computer Science and Information Engineering Changchun, China, June 17-19, 2012 p687

    [22]

    Peter T 1996 Ph. D. Dissertation (Lyngby: Technical University of Denmark)

  • [1] 朱军高, 卢海洋, 赵媛, 赖美福, 古永力, 徐世祥, 周沧涛. 面向激光驱动质子束应用的弱聚焦磁场束线设计研究. 物理学报, 2022, 71(19): 194102. doi: 10.7498/aps.71.20220599
    [2] 周腊珍, 夏文静, 许倩倩, 陈赞, 李坊佐, 刘志国, 孙天希. 一种基于毛细管X光透镜的微型锥束CT扫描仪. 物理学报, 2022, 71(9): 090701. doi: 10.7498/aps.71.20212195
    [3] 罗晓飞, 王波, 彭宽, 肖嘉莹. 基于聚焦声场模型的光声层析成像时间延迟快速校正反投影方法. 物理学报, 2022, 71(7): 078102. doi: 10.7498/aps.71.20212019
    [4] 席晓琦, 韩玉, 李磊, 闫镔. 螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法. 物理学报, 2019, 68(8): 088701. doi: 10.7498/aps.68.20190055
    [5] 戚俊成, 陈荣昌, 刘宾, 陈平, 杜国浩, 肖体乔. 基于迭代重建算法的X射线光栅相位CT成像. 物理学报, 2017, 66(5): 054202. doi: 10.7498/aps.66.054202
    [6] 刘鑫, 易明皓, 郭金川. 线焦斑X射线源成像. 物理学报, 2016, 65(21): 219501. doi: 10.7498/aps.65.219501
    [7] 孙亚秀, 卓庆坤, 姜庆辉, 李千. 基于多导体传输线理论的差模激励新型线束串扰模型研究. 物理学报, 2015, 64(4): 044102. doi: 10.7498/aps.64.044102
    [8] 陈平, 韩焱, 潘晋孝. 基于对数解调的递变能量CT成像方法. 物理学报, 2015, 64(13): 138701. doi: 10.7498/aps.64.138701
    [9] 韩玉, 李磊, 闫镔, 席晓琦, 胡国恩. 一种基于Radon逆变换的半覆盖螺旋锥束CT重建算法. 物理学报, 2015, 64(5): 058704. doi: 10.7498/aps.64.058704
    [10] 杨富强, 张定华, 黄魁东, 王鹍, 徐哲. CT不完全投影数据重建算法综述. 物理学报, 2014, 63(5): 058701. doi: 10.7498/aps.63.058701
    [11] 尹彬, 柏云龙, 齐艳辉, 冯素春, 简水生. 拉锥型啁啾光纤光栅滤波器的研究. 物理学报, 2013, 62(21): 214213. doi: 10.7498/aps.62.214213
    [12] 汪先超, 闫镔, 刘宏奎, 李磊, 魏星, 胡国恩. 一种圆轨迹锥束CT中截断投影数据的高效重建算法. 物理学报, 2013, 62(9): 098702. doi: 10.7498/aps.62.098702
    [13] 黄魁东, 张定华, 李明君, 张华. 锥束CT平板探测器成像的余晖建模与校正方法. 物理学报, 2013, 62(21): 210702. doi: 10.7498/aps.62.210702
    [14] 王国光, 王丹, 何丽桥. 混沌中信号的投影滤波. 物理学报, 2010, 59(5): 3049-3056. doi: 10.7498/aps.59.3049
    [15] 向良忠, 邢达, 郭华, 杨思华. 高分辨率快速数字化光声CT乳腺肿瘤成像. 物理学报, 2009, 58(7): 4610-4617. doi: 10.7498/aps.58.4610
    [16] 李永青, 李希国, 刘紫玉, 罗培燕, 张鹏鸣. Jackiw-Pi模型的新涡旋解. 物理学报, 2007, 56(11): 6178-6182. doi: 10.7498/aps.56.6178
    [17] 肖 沛, 张增明, 孙 霞, 丁泽军. 投影电子束光刻中电子穿透掩膜的Monte Carlo模拟. 物理学报, 2006, 55(11): 5803-5809. doi: 10.7498/aps.55.5803
    [18] 董小伟, 裴 丽, 简水生. 非对称熔锥法制作光纤光栅辅助耦合器型上下话路滤波器. 物理学报, 2006, 55(9): 4739-4743. doi: 10.7498/aps.55.4739
    [19] 陈雁萍, 王传兵, 周国成. 损失锥-束流分布电子驱动的回旋激射不稳定性. 物理学报, 2005, 54(7): 3221-3227. doi: 10.7498/aps.54.3221
    [20] 李鲠颖, 邬学文. 粉末样品中半整数四极核中心线的自旋锁定现象. 物理学报, 1989, 38(2): 338-342. doi: 10.7498/aps.38.338
计量
  • 文章访问数:  4546
  • PDF下载量:  401
  • 被引次数: 0
出版历程
  • 收稿日期:  2013-03-17
  • 修回日期:  2013-05-02
  • 刊出日期:  2013-08-05

/

返回文章
返回