搜索

x

留言板

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

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

抑制傅里叶变换法恢复的X射线相衬像中的伪影

杨君 吴浩 罗琨皓 郭金川 宗方轲

引用本文:
Citation:

抑制傅里叶变换法恢复的X射线相衬像中的伪影

杨君, 吴浩, 罗琨皓, 郭金川, 宗方轲

Suppression of artifacts in X-ray phase-contrast images retrieved by Fourier transform

Yang Jun, Wu Hao, Luo Kun-Hao, Guo Jin-Chuan, Zong Fang-Ke
PDF
HTML
导出引用
  • 在基于光栅的X射线相衬信号的恢复方法中, 主要有相移法和傅里叶变换法两种方法. 相移法具有精度高、噪声小的优点, 但由于至少需要三幅图像才能恢复出相衬信号, 样品所受的辐射剂量大. 而傅里叶变换法只需一幅图像即可恢复出物体的相衬信号, 具有快速、实时的优点, 但恢复出的信号精度低, 易受伪影影响. 因此, 本文利用两幅图像傅里叶变换法恢复X射线相衬信号, 该方法能够有效地抑制相衬信号中由于频谱混叠所产生的伪影. 另外, 通过增加载波条纹的频率, 能够拉大频域中的频谱间隔, 从而进一步抑制伪影的产生.
    Over the last two decades, the grating-based phase-contrast imaging has aroused the interest of a number of researchers. It could provide an access to three complementary signals simultaneously: the conventional absorption contrast, the differential phase contrast related to refraction of incident wave, and the dark-field contrast that relates to ultra small angle scattering in a sample. The grating-based phase-contrast signals have higher contrast sensitivity for some types of soft samples than the absorption signals. Dark-field signals have better diagnostic effects in the detection of lung tumors, pneumothorax and the identification of microcalcifications in breast. There are two main phase retrieval methods in grating-based X-ray phase-contrast imaging, i.e. phase stepping method and Fourier transform method. The phase signals retrieved by phase stepping is high precise and has low noise. But the sample suffers high dose due to at least three exposures. The phase signals retrieved by Fourier transform is low-dose due to the fact that only one image with sample is needed, but it is easily affected by artifacts when the size of the filtering window is too large. However, when the size of the filtering window is too small, the high-frequency information of the phase-contrast image will be lost, and the image will become blurred. A trade-off between definitions of the image and artifacts should be made. Since the phase-contrast signal and the dark-field signal of the sample are modulated by carrier fringes, the frequency spectrum of the detected image consists of many different harmonics. The artifacts in the retrieved signals originate from the spectrum aliasing between primary peak around zero spatial frequency and first-order harmonic peaks. Therefore, the subtraction between two images with phase difference can remove the primary peak, and the artifacts in the phase-contrast signals and dark-field signals will be suppressed. In order to further suppress the artifacts, we increase the frequency of carrier fringes, which results in a larger distance between first-order harmonic peaks in frequency domain. We finally attain artifact-free phase-contrast images and dark-field images while maintaining high definition of the images. The method proposed here is not only applicable to incoherent imaging system, but also to Talbot-Lau interferometer, and it would be useful in fast and low-dose X-ray phase-contrast and dark-field imaging.
      通信作者: 郭金川, jcguo@szu.edu.cn ; 宗方轲, zongfk168@163.com
    • 基金项目: 国家自然科学基金(批准号: 11674232)和广东省基础与应用基础研究基金(批准号: 2019A1515011785)资助的课题
      Corresponding author: Guo Jin-Chuan, jcguo@szu.edu.cn ; Zong Fang-Ke, zongfk168@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11674232) and Guangdong Provincial Basic and Applied Basic Research Foundation, China (Grant No. 2019A1515011785)
    [1]

    David C, Nohammer E, Solak H H, Ziegler E 2002 Appl. Phys. Lett. 81 3287Google Scholar

    [2]

    Momose A, Kawamoto S, Koyama I, Hamaishi Y, Takai K, Suzuki Y 2003 Jpn. J. Appl. Phys. 42 L866Google Scholar

    [3]

    Pfeiffer F, Weitkamp T, Bunk O, David C 2006 Nat. Phys. 2 258Google Scholar

    [4]

    Pfeiffer F, Bech M, Bunk O, Kraft P, Eikenberry E F, Bronnimann C, Grunzweig C, David C 2008 Nat. Mater. 7 134Google Scholar

    [5]

    Bech M, Tapfer A, Pauwels B, Bruyndonckx P, Sasov A, Pfeiffer F 2013 Sci. Rep. 3 3209Google Scholar

    [6]

    Anton G, Michel T, Pelzer G, Radicke M, Rieger J, Weber T 2013 Z. Med. Phys. 23 228Google Scholar

    [7]

    Yang J, Guo J C, Lei Y H, Yi M H, Chen L 2017 Chin. Phys. B 26 028701Google Scholar

    [8]

    Weitkamp T, Diaz A, David C, Pfeiffer F, Stampanoni M, Cloetens P, Ziegler E 2005 Opt. Express 13 6296Google Scholar

    [9]

    Takeda M, Ina H, Kobayashi S 1982 J. Opt. Soc. Am. 72 156Google Scholar

    [10]

    Wen H, Bennett E E, Hegedus M M, Carroll S C 2008 IEEE Trans. Med. Imaging 27 997Google Scholar

    [11]

    Wen H, Bennett E E, Hegedus M M, Rapacchi S 2009 Radiology 251 910Google Scholar

    [12]

    Lim H, Park Y, Cho H, Je U, Hong D, Park C, Woo T, Lee M, Kim J, Chung N, Kim J, Kim J 2015 Opt. Commun. 348 85Google Scholar

    [13]

    Lim H W, Lee H W, Cho H S, Je U K, Park C K, Kim K S, Kim G A, Park S Y, Lee D Y, Park Y O, Woo T H, Lee S H, Chung W H, Kim J W, Kim J G 2017 Nucl. Instrum. Methods Phys. Res., Sect. A 850 89Google Scholar

    [14]

    Lim H, Lee H, Cho H, Seo C, Je U, Park C, Kim K, Kim G, Park S, Lee D, Kang S, Lee M 2017 J. Korean Phys. Soc. 71 722Google Scholar

    [15]

    Seifert M, Gallersdörfer M, Ludwig V, Schuster M, Horn F, Pelzer G, Rieger J, Michel T, Anton G 2018 J. Imaging 4 62Google Scholar

    [16]

    Seifert M, Ludwig V, Gallersdorfer M, Hauke C, Hellbach K, Horn F, Pelzer G, Radicke M, Rieger J, Sutter S M, Michel T, Anton G 2018 Phys. Med. Biol. 63 185010Google Scholar

    [17]

    Li J, Su X Y, Guo L R 1990 Opt. Eng. 29 1439Google Scholar

    [18]

    陈文静, 苏显渝, 曹益平, 向立群 2004 中国激光 31 740Google Scholar

    Chen W J, Su X Y, Cao Y P, Xiang L Q 2004 Chin. J. Las. 31 740Google Scholar

    [19]

    Zhu P, Zhang K, Wang Z, Liu Y, Liu X, Wu Z, McDonald S A, Marone F, Stampanoni M 2010 Proc. Natl. Acad. Sci. U. S. A. 107 13576Google Scholar

    [20]

    Wang Z, Gao K, Ge X, Wu Z, Chen H, Wang S, Zhu P, Yuan Q, Huang W, Zhang K, Wu Z 2013 J. Phys. D: Appl. Phys. 46 494003Google Scholar

    [21]

    杜杨, 雷耀虎, 刘鑫, 郭金川, 牛憨笨 2013 物理学报 62 06872Google Scholar

    Yang D, Lei Y H, Liu X, Guo J C, Niu H B 2013 Acta Phys. Sin. 62 06872Google Scholar

    [22]

    Momose A, Yashiro W, Takeda Y, Suzuki Y, Hattori T 2006 Jpn. J. Appl. Phys. 45 5254Google Scholar

  • 图 1  非相干X射线成像系统示意图

    Fig. 1.  Schematic diagram of incoherent X-ray imaging.

    图 2  一般情形下载波条纹的频谱图

    Fig. 2.  The Fourier spectrum of carrier fringe patterns in general case.

    图 3  实际情况下发生的频谱混叠

    Fig. 3.  Spectrum aliasing between different harmonic peaks in practice.

    图 4  频谱混叠对恢复出相位的影响 (a), (b)和(c)分别代表载波条纹周期与探测器像素尺寸比值r为3, 4和5时的情形

    Fig. 4.  The impact of spectrum aliasing on phase retrieval. (a), (b) and (c) denote the cases, in which the ratios of the carrier fringe period to size of detector pixel are 3, 4 and 5, respectively.

    图 5  无频谱混叠时恢复出来的相位分布

    Fig. 5.  The phase retrieval when no spectrum aliasing.

    图 6  (a) 两幅图像傅里叶变换法所恢复出的PMMA相衬像; (b)单幅图像傅里叶变换法所恢复出的PMMA相衬像

    Fig. 6.  (a) The phase-contrast image of PMMA retrieved by Fourier transform with two images; (b) the phase-contrast image of PMMA retrieved by Fourier transform with one image.

    图 7  (a) 图6中白色方框区域按行取平均后绘制的曲线; (b) 图6中黑色方框区域按行取平均后绘制的曲线

    Fig. 7.  (a) Curves from averaging the area of white rectangle in figure 6 by row; (b) curves from averaging the area of black rectangle in Figure 6 by row.

    图 8  载波条纹周期与探测器像素尺寸比值r = 3时, 鸡翅的频谱 (a)、相衬像(c)和暗场像(e). r = 5时, 鸡翅的频谱(b)、相衬像(d)和暗场像(f)

    Fig. 8.  Fourier spectrum (a), phase-contrast image (c) and dark-field image (e) of a chicken wing when r = 3. Fourier spectrum (b), phase-contrast image (d) and dark-field image (f) of a chicken wing when r = 5.

    表 1  两种不同傅里叶变换方法的定量比较

    Table 1.  Quantitative comparison between two kinds of Fourier transform algorithms.

    背景相位
    均值/rad
    背景相位标
    准差/rad
    横截面峰
    峰值/rad
    单幅图像
    傅里叶变换
    0.35020.00590.2412
    两幅图像
    傅里叶变换
    0.25260.00170.1112
    下载: 导出CSV
  • [1]

    David C, Nohammer E, Solak H H, Ziegler E 2002 Appl. Phys. Lett. 81 3287Google Scholar

    [2]

    Momose A, Kawamoto S, Koyama I, Hamaishi Y, Takai K, Suzuki Y 2003 Jpn. J. Appl. Phys. 42 L866Google Scholar

    [3]

    Pfeiffer F, Weitkamp T, Bunk O, David C 2006 Nat. Phys. 2 258Google Scholar

    [4]

    Pfeiffer F, Bech M, Bunk O, Kraft P, Eikenberry E F, Bronnimann C, Grunzweig C, David C 2008 Nat. Mater. 7 134Google Scholar

    [5]

    Bech M, Tapfer A, Pauwels B, Bruyndonckx P, Sasov A, Pfeiffer F 2013 Sci. Rep. 3 3209Google Scholar

    [6]

    Anton G, Michel T, Pelzer G, Radicke M, Rieger J, Weber T 2013 Z. Med. Phys. 23 228Google Scholar

    [7]

    Yang J, Guo J C, Lei Y H, Yi M H, Chen L 2017 Chin. Phys. B 26 028701Google Scholar

    [8]

    Weitkamp T, Diaz A, David C, Pfeiffer F, Stampanoni M, Cloetens P, Ziegler E 2005 Opt. Express 13 6296Google Scholar

    [9]

    Takeda M, Ina H, Kobayashi S 1982 J. Opt. Soc. Am. 72 156Google Scholar

    [10]

    Wen H, Bennett E E, Hegedus M M, Carroll S C 2008 IEEE Trans. Med. Imaging 27 997Google Scholar

    [11]

    Wen H, Bennett E E, Hegedus M M, Rapacchi S 2009 Radiology 251 910Google Scholar

    [12]

    Lim H, Park Y, Cho H, Je U, Hong D, Park C, Woo T, Lee M, Kim J, Chung N, Kim J, Kim J 2015 Opt. Commun. 348 85Google Scholar

    [13]

    Lim H W, Lee H W, Cho H S, Je U K, Park C K, Kim K S, Kim G A, Park S Y, Lee D Y, Park Y O, Woo T H, Lee S H, Chung W H, Kim J W, Kim J G 2017 Nucl. Instrum. Methods Phys. Res., Sect. A 850 89Google Scholar

    [14]

    Lim H, Lee H, Cho H, Seo C, Je U, Park C, Kim K, Kim G, Park S, Lee D, Kang S, Lee M 2017 J. Korean Phys. Soc. 71 722Google Scholar

    [15]

    Seifert M, Gallersdörfer M, Ludwig V, Schuster M, Horn F, Pelzer G, Rieger J, Michel T, Anton G 2018 J. Imaging 4 62Google Scholar

    [16]

    Seifert M, Ludwig V, Gallersdorfer M, Hauke C, Hellbach K, Horn F, Pelzer G, Radicke M, Rieger J, Sutter S M, Michel T, Anton G 2018 Phys. Med. Biol. 63 185010Google Scholar

    [17]

    Li J, Su X Y, Guo L R 1990 Opt. Eng. 29 1439Google Scholar

    [18]

    陈文静, 苏显渝, 曹益平, 向立群 2004 中国激光 31 740Google Scholar

    Chen W J, Su X Y, Cao Y P, Xiang L Q 2004 Chin. J. Las. 31 740Google Scholar

    [19]

    Zhu P, Zhang K, Wang Z, Liu Y, Liu X, Wu Z, McDonald S A, Marone F, Stampanoni M 2010 Proc. Natl. Acad. Sci. U. S. A. 107 13576Google Scholar

    [20]

    Wang Z, Gao K, Ge X, Wu Z, Chen H, Wang S, Zhu P, Yuan Q, Huang W, Zhang K, Wu Z 2013 J. Phys. D: Appl. Phys. 46 494003Google Scholar

    [21]

    杜杨, 雷耀虎, 刘鑫, 郭金川, 牛憨笨 2013 物理学报 62 06872Google Scholar

    Yang D, Lei Y H, Liu X, Guo J C, Niu H B 2013 Acta Phys. Sin. 62 06872Google Scholar

    [22]

    Momose A, Yashiro W, Takeda Y, Suzuki Y, Hattori T 2006 Jpn. J. Appl. Phys. 45 5254Google Scholar

  • [1] 李文文, 惠宁菊, 李存霞, 刘洋河, 方妍, 李凌青, 王彦龙, 唐远河. 多普勒非对称空间外差仪探测高层大气风速的三种方法比较研究. 物理学报, 2023, 72(24): 240601. doi: 10.7498/aps.72.20231292
    [2] 王子硕, 刘磊, 刘晨博, 刘珂, 钟志, 单明广. 数字差分-积分快速相位解包裹算法研究. 物理学报, 2023, 72(18): 184201. doi: 10.7498/aps.72.20230473
    [3] 陈星宇, 周昕, 白星, 余展, 王玉杰, 李欣家, 刘洋, 孙铭泽. 傅里叶鬼成像与正弦鬼成像的等价性分析. 物理学报, 2023, 72(14): 144202. doi: 10.7498/aps.72.20222317
    [4] 姚春霞, 何其利, 张锦, 付天宇, 吴朝, 王山峰, 黄万霞, 袁清习, 刘鹏, 王研, 张凯. 免分析光栅一次曝光相位衬度成像方法. 物理学报, 2021, 70(2): 028701. doi: 10.7498/aps.70.20201170
    [5] 王振, 杜艳君, 丁艳军, 彭志敏. 基于傅里叶变换的波长扫描腔衰荡光谱. 物理学报, 2019, 68(20): 204204. doi: 10.7498/aps.68.20191062
    [6] 于文婷, 张娟, 唐军. 动态突触、神经耦合与时间延迟对神经元发放的影响. 物理学报, 2017, 66(20): 200201. doi: 10.7498/aps.66.200201
    [7] 荣锋, 谢艳娜, 邰雪凤, 耿磊. 双能X射线光栅相衬成像的研究. 物理学报, 2017, 66(1): 018701. doi: 10.7498/aps.66.018701
    [8] 于慧, 张瑞, 李克武, 薛锐, 王志斌. 双强度调制静态傅里叶变换偏振成像光谱系统测量原理及仿真. 物理学报, 2017, 66(5): 054201. doi: 10.7498/aps.66.054201
    [9] 王华英, 于梦杰, 刘飞飞, 江亚男, 宋修法, 高亚飞. 基于同态信号处理的数字全息广义线性重建算法研究. 物理学报, 2013, 62(23): 234207. doi: 10.7498/aps.62.234207
    [10] 尹增谦, 武臣, 宫琬钰, 龚之珂, 王永杰. Voigt线型函数及其最大值的研究. 物理学报, 2013, 62(12): 123301. doi: 10.7498/aps.62.123301
    [11] 杜杨, 雷耀虎, 刘鑫, 郭金川, 牛憨笨. 硬X射线光栅微分干涉相衬成像两步相移算法的理论与实验研究. 物理学报, 2013, 62(6): 068702. doi: 10.7498/aps.62.068702
    [12] 杨强, 刘鑫, 郭金川, 雷耀虎, 黄建衡, 牛憨笨. 无吸收光栅的X射线相位衬度成像实验研究. 物理学报, 2012, 61(16): 160702. doi: 10.7498/aps.61.160702
    [13] 张淑娜, 罗震岳, 沈伟东, 刘旭, 章岳光. 利用白光干涉技术测量块状材料的群折射率. 物理学报, 2011, 60(1): 014221. doi: 10.7498/aps.60.014221
    [14] 张记全, 张延惠, 周慧, 贾正茂, 林圣路. 圆环弹子球量子谱的衍射效应. 物理学报, 2009, 58(9): 5965-5969. doi: 10.7498/aps.58.5965
    [15] 相里斌, 袁艳, 吕群波. 傅里叶变换光谱成像仪光谱传递函数研究. 物理学报, 2009, 58(8): 5399-5405. doi: 10.7498/aps.58.5399
    [16] 赵保银, 吕百达. 使用离焦望远镜系统合成轴上平顶光束的一种新方法. 物理学报, 2008, 57(5): 2919-2924. doi: 10.7498/aps.57.2919
    [17] 刘 斌, 金伟其, 董立泉. 热成像系统前置栅网结构的衍射效应分析. 物理学报, 2008, 57(9): 5578-5583. doi: 10.7498/aps.57.5578
    [18] 谢红兰, 高鸿奕, 陈建文, 王寯越, 朱佩平, 熊诗圣, 洗鼎昌, 徐至展. 具有原子分辨率的x射线荧光全息术的数值模拟研究. 物理学报, 2003, 52(9): 2223-2228. doi: 10.7498/aps.52.2223
    [19] 金 鹏, 潘士宏, 梁基本. SIN+ GaAs结构中的Franz-Keldysh振荡的傅里叶变换研究. 物理学报, 2000, 49(9): 1821-1828. doi: 10.7498/aps.49.1821
    [20] 刘 刚, 卢 柯, 梁志德. X射线线形的解卷分析与Bragg角的精确计算. 物理学报, 2000, 49(8): 1520-1523. doi: 10.7498/aps.49.1520
计量
  • 文章访问数:  5425
  • PDF下载量:  71
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-27
  • 修回日期:  2020-12-11
  • 上网日期:  2021-05-09
  • 刊出日期:  2021-05-20

/

返回文章
返回