留言板

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

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

基于加权全变差最小化的中子外部CT重建算法

林强 杨民 张晓敏 唐彬 刘斌

林强, 杨民, 张晓敏, 等 . 基于加权全变差最小化的中子外部CT重建算法[J]. 北京航空航天大学学报, 2021, 47(6): 1199-1209. doi: 10.13700/j.bh.1001-5965.2020.0126
引用本文: 林强, 杨民, 张晓敏, 等 . 基于加权全变差最小化的中子外部CT重建算法[J]. 北京航空航天大学学报, 2021, 47(6): 1199-1209. doi: 10.13700/j.bh.1001-5965.2020.0126
LIN Qiang, YANG Min, ZHANG Xiaomin, et al. Neutron exterior CT reconstruction algorithm based on weighted total variation minimization[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(6): 1199-1209. doi: 10.13700/j.bh.1001-5965.2020.0126(in Chinese)
Citation: LIN Qiang, YANG Min, ZHANG Xiaomin, et al. Neutron exterior CT reconstruction algorithm based on weighted total variation minimization[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(6): 1199-1209. doi: 10.13700/j.bh.1001-5965.2020.0126(in Chinese)

基于加权全变差最小化的中子外部CT重建算法

doi: 10.13700/j.bh.1001-5965.2020.0126
基金项目: 

国家自然科学基金 11675012

国家自然科学基金 U1730116

核安全项目 41426020401

详细信息
    通讯作者:

    杨民, E-mail: minyang.ndt@buaa.edu.cn

  • 中图分类号: TP391.4

Neutron exterior CT reconstruction algorithm based on weighted total variation minimization

Funds: 

National Natural Science Foundation of China 11675012

National Natural Science Foundation of China U1730116

Nuclear Safety Project 41426020401

More Information
  • 摘要:

    针对中子外部CT检测需求,提出了基于加权方向全变差(WDTV)最小化的中子外部CT重建算法。首先,利用平行束数据对称性原理将正弦图中旋转中心另一侧缺失的投影数据进行补充;其次,采用传统FBP算法进行了外部CT扫描模式下的CT重建。为了抑制传统FBP、SART算法重建图像中径向边缘伪影,采用WDTV算法计算沿径向和周向方向的局部方向差分并作加权和。另外,在重建模型的WDTV项中引入2个权重参数以控制径向边缘和切向边缘的不同响应强度。最后,采用所提算法研究了重建图像质量与投影数据量的关系。计算机仿真和真实冷中子实验表明:基于WDTV最小化的中子外部CT重建算法能够有效抑制径向边缘伪影,获得高质量重建图像。

     

  • 图 1  中子外部CT扫描原理

    Figure 1.  Principle of neutron exterior CT scanning

    图 2  外部CT扫描模体

    Figure 2.  Sample of exterior CT scanning

    图 3  模体外部CT扫描示意

    Figure 3.  Schematic diagram of exterior CT scanning of sample

    图 4  模体FBP算法重建图像

    Figure 4.  Reconstruction image of sample by FBP algorithm

    图 5  局部方向差分计算方法[23]

    Figure 5.  Calculation method of local directional difference[23]

    图 6  探测器宽度为128时不同算法重建结果

    Figure 6.  Reconstruction results from different algorithms when detector width is 128

    图 7  不同探测器宽度外部CT扫描示意图

    Figure 7.  Schematic diagram of exterior CT scanning with different detector width

    图 8  不同探测器宽度时SART-WDTV算法重建图像

    Figure 8.  Images reconstructed by SART-WDTV algorithm with different detector width

    图 9  不同探测器宽度重建图像的MSE和PSNR

    Figure 9.  MSE and PSNR of reconstruction images with different detector width

    图 10  齿轮样品投影正弦图与SART-TV算法重建图像

    Figure 10.  Complete projection sinogram and SART-TV algorithm reconstruction image of gear sample

    图 11  正弦图宽度为200时FBP算法与SART算法重建图像

    Figure 11.  Images reconstructed by FBP and SART algorithms when sinogram width is 200

    图 12  正弦图宽度为200时SART-TV算法与SART-WDTV算法重建图像

    Figure 12.  Images reconstructed by SART-TV and SART-WDTV algorithms when sinogram width is 200

    图 13  正弦图宽度为122时FBP算法与SART算法重建图像

    Figure 13.  Images reconstructed by FBP and SART algorithms when sinogram width is 122

    图 14  正弦图宽度为122时SART-TV算法与SART-WDTV算法重建图像

    Figure 14.  Images reconstructed by SART-TV and SART-WDTV algorithms when sinogram width is 122

    表  1  外部CT扫描参数

    Table  1.   Parameters of exterior CT scanning

    CT模式 探元尺寸 像素数量 角度范围 投影数量 模体大小
    平行束 1像素 128像素 0°~360° 720 512像素×512像素
    下载: 导出CSV

    表  2  模体外部CT重建图像的MSE和PSNR

    Table  2.   MSE and PSNR of exterior CT reconstruction images of sample

    重建算法 MSE PSNR
    FBP 3 050.040 13.287 7
    SART 1 319.476 16.926 8
    SART-TV 1 034.756 17.982 4
    SART-WDTV 874.629 8 18.712 6
    下载: 导出CSV

    表  3  样品外部CT重建图像的MSE和PSNR

    Table  3.   MSE and PSNR of exterior CT reconstruction images of sample

    重建算法 正弦图宽度200 正弦图宽度122
    MSE PSNR MSE PSNR
    FBP 1 048.667 17.924 4 1 438.552 16.551 5
    SART 762.909 19.306 1 1 436.950 16.556 4
    SART-TV 499.984 21.141 2 925.194 18.468 5
    SART-WDTV 431.721 21.778 8 923.005 18.478 8
    下载: 导出CSV
  • [1] YASUDA R, MATSUBAYASHI M, NAKATA M, et al. Application of neutron imaging plate and neutron CT methods on nuclear fuels and materials[J]. IEEE Transactions on Nuclear Science, 2005, 52(1): 313-316. doi: 10.1109/TNS.2005.844297
    [2] SINHA V, SRIVASTAVA A, LEE H K. A novel method for NDT applications using NXCT system at the Missouri University of Science & Technology[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2014, 750(1): 43-55. http://www.sciencedirect.com/science/article/pii/S0168900214002721
    [3] TAMAKI M. Conceptual monochromatic digital neutron radiography using continuous cold neutron beam[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2005, 542(1-3): 32-37. doi: 10.1016/j.nima.2005.01.008
    [4] 魏国海, 韩松柏, 陈东风, 等. 中子照相技术在核燃料元件无损检测中的应用[J]. 核技术, 2012, 35(11): 821-826. https://www.cnki.com.cn/Article/CJFDTOTAL-HJSU201211005.htm

    WEI G H, HAN S B, CHEN D F, et al. Application of neutron radiography for non-destructive testing nuclear fuel elements[J]. Nuclear Techniques, 2012, 35(11): 821-826(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HJSU201211005.htm
    [5] STROBL M, MANKE I, KARDJILOV N, et al. Advances in neutron radiography and tomography[J]. Journal of Physics D: Applied Physics, 2009, 42(24): 243001. doi: 10.1088/0022-3727/42/24/243001
    [6] 郭雨濛, 曾理. CT外问题的原理与算法综述[J]. 中国体视学与图像分析, 2017, 22(3): 320-326. https://www.cnki.com.cn/Article/CJFDTOTAL-ZTSX201703011.htm

    GUO Y M, ZENG L. Review on the principle and algorithms of the exterior problem of CT[J]. Chinese Journal of Stereology and Image Analysis, 2017, 22(3): 320-326(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZTSX201703011.htm
    [7] 倪风岳. CUDA加速CV图像分割和外部CT图像重建算法研究[D]. 重庆: 重庆大学, 2011: 2-30.

    NI F Y. Algorithm study on accelerate CV image segmentation and exterior industrial image reconstruction by CUDA[D]. Chongqing: Chongqing University, 2011: 2-30(in Chinese).
    [8] GUO J Q, ZENG L, LIU B D. High-quality image reconstruction from exterior helical cone-beam CT data for NDE of industrial pipelines[J]. Insight-Non-Destructive Testing and Condition Monitoring, 2011, 53(10): 534-541. doi: 10.1784/insi.2011.53.10.534
    [9] 陈庆贵, 卢洪义, 齐强, 等. 固体火箭发动机界面脱粘切向CT检测[J]. 固体火箭技术, 2016, 39(3): 347-352. https://www.cnki.com.cn/Article/CJFDTOTAL-GTHJ201603010.htm

    CHEN Q G, LU H Y, QI Q, et al. Tangential CT inspection of interface debonding of SRM[J]. Journal of Solid Rocket Technology, 2016, 39(3): 347-352(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GTHJ201603010.htm
    [10] CHEN B X, YANG M, ZHANG Z, et al. Constrained TV-minimization reconstruction from exterior CT data[C]//2013 IEEE Nuclear Science Symposium and Medical Imaging Conference. Piscataway: IEEE Press, 2013: 1-3.
    [11] QUINTO E T. Local algorithms in exterior tomography[J]. Journal of Computational & Applied Mathematics, 2007, 199(1): 141-148. http://www.ams.org/mathscinet-getitem?mr=2267538
    [12] GUO J Q, ZENG L, LIU B D. The improved TVM image reconstruction algorithm of exterior circular cone-beam CT[C]//Proceedings of the 13th International Congress for Stereology. Beijing: Chinese Society for Stereology, 2011: 458-462.
    [13] CORMACK A M. Representation of a function by its line integrals, with some radiological applications[J]. Journal of Applied Physics, 1963, 34(9): 2722-2727. doi: 10.1063/1.1729798
    [14] QUINTO E T. Tomographic reconstructions from incomplete data-numerical inversion of the exterior Radon transform[J]. Inverse Problems, 1988, 4(3): 867-876. doi: 10.1088/0266-5611/4/3/019
    [15] QUINTO E T. Singularities of the X-ray transform and limited data tomography in R2 and R3[J]. SIAM Journal on Mathematical Analysis, 1993, 24(5): 1215-1225. doi: 10.1137/0524069
    [16] QUINTO E T. Exterior and limited-angle tomography in non-destructive evaluation[J]. Inverse Problems, 1998, 14(2): 339-353. doi: 10.1088/0266-5611/14/2/009
    [17] QUINTO E T. Singular value decompositions and inversion methods for the exterior Radon transform and a spherical transform[J]. Journal of Mathematical Analysis and Applications, 1983, 95(2): 437-448. doi: 10.1016/0022-247X(83)90118-X
    [18] VAINBERG E, KAZAK I, KUROZAEV V. Reconstruction of the internal three-dimensional structure of objects based on real-time integral projections[J]. Soviet Journal of Nondestructive Testing, 1981, 17(6): 415-423. http://www.researchgate.net/publication/280136857_RECONSTRUCTION_OF_THE_INTERNAL_THREE-DIMENSIONAL_STRUCTURE_OF_OBJECTS_BASED_ON_REAL-TIME_INTEGRAL_PROJECTIONS
    [19] 蔺鲁萍, 王永革. 不完全角度CT图像重建的模型与算法[J]. 北京航空航天大学学报, 2017, 43(4): 823-830. doi: 10.13700/j.bh.1001-5965.2016.0232

    LIN L P, WANG Y G. CT image reconstruction model and algorithm from few views[J]. Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(4): 823-830(in Chinese). doi: 10.13700/j.bh.1001-5965.2016.0232
    [20] SIDKY E Y, KAO C M, PAN X. Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT[J]. Journal of X-Ray Science and Technology, 2006, 14(2): 119-139. http://arxiv.org/abs/0904.4495v1
    [21] ZENG L, LIU B, LIU L, et al. A new iterative reconstruction algorithm for 2D exterior fan-beam CT[J]. Journal of X-Ray Science and Technology, 2010, 18(3): 267-277. http://europepmc.org/abstract/med/20714085
    [22] GUO Y, ZENG L. Improved iterative image reconstruction algorithm for the exterior problem of computed tomography[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2017, 842: 96-108. http://www.sciencedirect.com/science/article/pii/S0168900216310993
    [23] GUO Y, ZENG L, WANG C, et al. Image reconstruction model for the exterior problem of computed tomography based on weighted directional total variation[J]. Applied Mathematical Modelling, 2017, 52: 358-377. doi: 10.1016/j.apm.2017.07.057
    [24] 闫镔, 李磊. CT图像重建算法[M]. 北京: 科学出版社, 2014: 5-14.

    YAN B, LI L. CT image reconstruction algorithm[M]. Beijing: Science Press, 2014: 5-14(in Chinese).
    [25] HUYNH-THU Q, GHANBARI M. The accuracy of PSNR in predicting video quality for different video scenes and frame rates[J]. Telecommunication Systems, 2012, 49(1): 35-48. doi: 10.1007/s11235-010-9351-x
    [26] HORE A, ZIOU D, Image quality metrics: PSNR vs. SSIM[C]//2010 20th International Conference on Pattern Recognition. Piscataway: IEEE Press, 2010: 2366-2369.
  • 加载中
图(14) / 表(3)
计量
  • 文章访问数:  512
  • HTML全文浏览量:  111
  • PDF下载量:  41
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-07
  • 录用日期:  2020-07-03
  • 网络出版日期:  2021-06-20

目录

    /

    返回文章
    返回
    常见问答