留言板

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

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

扫频光学相干层析视网膜图像配准去噪算法

蔡怀宇,韩晓艳,娄世良,汪毅,陈文光,陈晓冬

downloadPDF
蔡怀宇, 韩晓艳, 娄世良, 汪毅, 陈文光, 陈晓冬. 扫频光学相干层析视网膜图像配准去噪算法[J]. , 2021, 14(2): 289-297. doi: 10.37188/CO.2020-0130
引用本文: 蔡怀宇, 韩晓艳, 娄世良, 汪毅, 陈文光, 陈晓冬. 扫频光学相干层析视网膜图像配准去噪算法[J]. , 2021, 14(2): 289-297.doi:10.37188/CO.2020-0130
CAI Huai-yu, HAN Xiao-yan, LOU Shi-liang, WANG Yi, CHEN Wen-guang, CHEN Xiao-dong. Speckle noise reduction in swept-source optical coherence tomography by retinal image registration[J]. Chinese Optics, 2021, 14(2): 289-297. doi: 10.37188/CO.2020-0130
Citation: CAI Huai-yu, HAN Xiao-yan, LOU Shi-liang, WANG Yi, CHEN Wen-guang, CHEN Xiao-dong. Speckle noise reduction in swept-source optical coherence tomography by retinal image registration[J].Chinese Optics, 2021, 14(2): 289-297.doi:10.37188/CO.2020-0130

扫频光学相干层析视网膜图像配准去噪算法

doi:10.37188/CO.2020-0130
基金项目:国家重点研发计划(No. 2017YFC0109901);天津市自然科学基金项目(No. 15JCQNJC14200)
详细信息
    作者简介:

    蔡怀宇(1965—),女,湖南涟源人,博士,教授,硕士生导师,1991年、2000年于天津大学分别获得硕士、博士学位,主要从事信息光学、光电技术及仪器和图像处理等方面的研究。E-mail:hycai@tju.edu.cn

    韩晓艳(1996—),女,内蒙古乌兰察布人,天津大学精密仪器与光电工程技术学院硕士研究生,2018年于南京理工大学获得学士学位,主要从事光学相干层析成像方面的研究。E-mail:xyhan@tju.edu.cn

  • 中图分类号:TN247

Speckle noise reduction in swept-source optical coherence tomography by retinal image registration

Funds:Supported by National Key R&D Program of China (No. 2017YFC0109901); Natural Science Foundation Project of Tianjin (No. 15JCQNJC14200)
More Information
  • 摘要:多帧叠加平均处理是去除扫频光学相干层析系统散斑噪声、获得较为清晰结构信息的有效方法,但眼睛的震颤、漂移、微眼跳等生理特性和系统光路特性会使图像之间存在错位,导致叠加效果不佳、结构稳定性差,为此本文提出一种基于灰度分布信息和目标几何信息相结合的配准算法。该方法根据图像平均灰度分布提取包含目标信息的感兴趣区域,通过相位相关算法和基于分段拟合的灰度投影算法的双重作用校正图像的平移变换;通过拟合视网膜上边界作为特征点迭代确定最佳旋转参数,并再次重新估计平移参数,实现图像的刚性配准;最后通过轴向扫描一对一映射法以能量函数为约束条件实现图像的非刚性配准。对活体兔眼进行实验,结果表明,本文算法配准后的叠加图像边界清晰,结构信息增强,信噪比和对比度平均有效提高一倍多。本算法适用于强噪声视网膜B-Scans图像的配准,能满足多种类型OCT系统的叠加成像需要,具有较高的鲁棒性和图像配准精度。

  • 图 1OCT序列图像错位原因。(a)眼睛轴向运动;(b)眼睛垂轴运动;(c)眼睛旋转运动;(d)多角度光线入射

    Figure 1.Analysis of reasons for geometric transformation of OCT sequence images. (a) Axial eye movement; (b) vertical eye movement; (c) rotational eye movement; (d) multi-angle light incidence

    图 2配准算法流程图

    Figure 2.Flow chart of the registration algorithm

    图 3预处理过程。(a)原始视网膜B-Scan图像;(b)列平均灰度分布图;(c)ROI图像

    Figure 3.Pretreatment process. (a) Original B-Scan image of retinal image; (b) column average gray distribution map; (c) ROI image

    图 4分段拟合处理过程。(a)曲线分割结果(黑色虚线为分割线);(b)原始曲线段(最佳连接点A、黑色直线为包络线);(c)处理后的曲线段

    Figure 4.Fitting process of curve segments. (a) Results of curve segmentation (the dotted black line is the segmentation line); (b) original curve segment results (optimal connection point A, the straight black line is the envelope); (c) curve segment after processing

    图 5刚性配准结果。(a)上边界提取过程;(b)平移配准后的叠加图;(c)旋转配准后的叠加图(绿色为参考图像,紫色为目标图像)

    Figure 5.Results of rigid registration. (a) Upper boundary extraction process; (b) averaging image after translation registration; (c) averaging image after rotation registration. (Green is the reference image, purple is the target image)

    图 6非刚性配准结果。(a)刚性配准后的叠加图及局部放大图;(b)非刚性配准后的叠加图及局部放大图(绿色为参考图像,紫色为目标图像)

    Figure 6.Results of non-rigid registration. (a) Average image obtained after rigid registration; (b) average image obtained after non-rigid registration (The red dashed box represents an enlarged partial view; Green is the reference image, and purple is the target image)

    图 7本算法实验结果图。(a)原始单帧图;(b)配准后的叠加图

    Figure 7.The experimental results obtained by the algorithm proposed in this paper. (a) Original single frame image; (b) averaging image after registration

    图 8(a)区域选择及SNR(b)和CNR(c)随叠加帧数递增的曲线变化图

    Figure 8.(a) Area selection, SNR (b) and CNR (c) curve change with an increasing number of multiple frames

    表 1本文配准算法较对比算法的提升效果

    Table 1.Improvement performance of proposed algorithm comparing with other algorithms

    比较对象 SNR提高倍数 CNR提高倍数
    本文方法/文献[18]方法 1.12 1.19
    本文方法/文献[21]方法 1.50 1.48
    下载: 导出CSV
  • [1] SAHYOUN C C, SUBHASH H M, PERU D,et al. An experimental review of optical coherence tomography systems for noninvasive assessment of hard dental tissues[J].Caries Research, 2019, 54(1): 43-54.
    [2] PASAOGLU I, SATANA B, ALTAN C,et al. Lamina cribrosa surface position in idiopathic intracranial hypertension with swept-source optical coherence tomography[J].Indian Journal of Ophthalmology, 2019, 67(7): 1085-1088.doi:10.4103/ijo.IJO_1736_18
    [3] PARK J H, YOO C, JUNG J H,et al. The association between prelaminar tissue thickness and peripapillary choroidal thickness in untreated normal-tension glaucoma patients[J].Medicine, 2019, 98(1): e14044.doi:10.1097/MD.0000000000014044
    [4] 汪毅, 刘珊珊, 张玮茜, 等. 扫频光学相干层析角膜图像轮廓自动提取算法[J]. 物理学报,2019,68(20):204201.doi:10.7498/aps.68.20190731

    WANG Y, LIU SH SH, ZHANG W Q,et al. Automatic contour extraction algorithm for swept-source optical coherence tomography cornea image[J].Acta Physica Sinica, 2019, 68(20): 204201. (in Chinese)doi:10.7498/aps.68.20190731
    [5] VIRA J, MARCHESE A, SINGH R B,et al. Swept-source optical coherence tomography imaging of the retinochoroid and beyond[J].Expert Review of Medical Devices, 2020, 17(5): 413-426.doi:10.1080/17434440.2020.1755256
    [6] LOU SH L, CHEN X D, HAN X Y,et al. Fast retinal segmentation based on the wave algorithm[J].IEEE Access, 2020, 8: 53678-53686.doi:10.1109/ACCESS.2020.2981206
    [7] SCHMITT J M, XIANG S H, YUNG K M,et al. Speckle in optical coherence tomography[J].Journal of Biomedical Optics, 1999, 4(1): 95-105.doi:10.1117/1.429925
    [8] 蔡怀宇, 张玮茜, 陈晓冬, 等. 眼科光学相干层析成像的图像处理方法[J]. 中国光学,2019,12(4):731-740.doi:10.3788/co.20191204.0731

    CAI H Y, ZHANG W Q, CHEN X D,et al. Image processing method for ophthalmic optical coherence tomography[J].Chinese Optics, 2019, 12(4): 731-740. (in Chinese)doi:10.3788/co.20191204.0731
    [9] JORGENSEN T M, THOMADSEN J, CHRISTENSEN U,et al. Enhancing the signal-to-noise ratio in ophthalmic optical coherence tomography by image registration—method and clinical examples[J].Journal of Biomedical Optics, 2007, 12(4): 041208.doi:10.1117/1.2772879
    [10] SZKULMOWSKI M, WOJTKOWSKI M. Averaging techniques for OCT imaging[J].Optics Express, 2013, 21(8): 9757-9773.doi:10.1364/OE.21.009757
    [11] MARTINEZ-CONDE S, MACKNIK S L, HUBEL D H. The role of Fixational eye movements in visual perception[J].Nature Reviews Neuroscience, 2004, 5(3): 229-240.doi:10.1038/nrn1348
    [12] KLEIN T, HUBER R. High-speed OCT light sources and systems [Invited][J].Biomedical Optics Express, 2017, 8(2): 828-859.doi:10.1364/BOE.8.000828
    [13] 余轮, 魏丽芳. 眼底图像配准技术研究进展[J]. 生物医学工程学杂志,2011,28(5):1043-1047.

    YU L, WEI L F. Progress of research in retinal image registration[J].Journal of Biomedical Engineering, 2011, 28(5): 1043-1047. (in Chinese)
    [14] NIEMEIJER M, GARVIN M K, LEE K,et al. Registration of 3D spectral OCT volumes using 3D SIFT feature point matching[J].Proceedings of SPIE, 2009, 7259: 72591I.
    [15] NIEMEIJER M, LEE K, GARVIN M K,et al. Registration of 3D spectral OCT volumes combining ICP with a graph-based approach[J].Proceedings of SPIE, 2012, 8314: 83141A.
    [16] CHEN M, LANG A, YING H S,et al. Analysis of macular OCT images using deformable registration[J].Biomedical Optics Express, 2014, 5(7): 2196-2214.doi:10.1364/BOE.5.002196
    [17] KHANSARI M M, ZHANG J, QIAO Y CH,et al. Automated deformation-based analysis of 3D optical coherence tomography in diabetic retinopathy[J].IEEE Transactions on Medical Imaging, 2020, 39(1): 236-245.doi:10.1109/TMI.2019.2924452
    [18] CHENG J, DUAN L X, WONG D W K,et al.. Speckle reduction in optical coherence tomography by image registration and matrix completion[C].Proceedings of the 17thInternational Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, 2014: 162-169.
    [19] LEE S, LEBED E, SARUNIC M V,et al. Exact surface registration of retinal surfaces from 3-D optical coherence tomography images[J].IEEE Transactions on Biomedical Engineering, 2015, 62(2): 609-617.doi:10.1109/TBME.2014.2361778
    [20] LEZAMA J, MUKHERJEE D, MCNABB R P,et al. Segmentation guided registration of wide field-of-view retinal optical coherence tomography volumes[J].Biomedical Optics Express, 2016, 7(12): 4827-4846.doi:10.1364/BOE.7.004827
    [21] DU X Y, GONG L, SHI F,et al..Non-rigid Registration of Retinal OCT Images Using Conditional Correlation Ratio[M]. CARDOSO M J, ARBEL T, MELBOURNE A,et al.. Fetal, Infant and Ophthalmic Medical Image Analysis. Cham: Springer, 2017: 159-167.
    [22] OJANSIVU V, HEIKKILA J. Image registration using blur-invariant phase correlation[J].IEEE Signal Processing Letters, 2007, 14(7): 449-452.doi:10.1109/LSP.2006.891338
    [23] 孙辉, 马天玮. 基于相位相关的目标图像亚像元运动参数估计[J]. 液晶与显示,2011,26(6):858-862.doi:10.3788/YJYXS20112606.0858

    SUN H, MA T W. Sub-pixel motion estimation based on phase-only correlation[J].Chinese Journal of Liquid Crystals and Displays, 2011, 26(6): 858-862. (in Chinese)doi:10.3788/YJYXS20112606.0858
    [24] 万钇良, 王建立, 张楠, 等. 一种基于相位相关与子图像的偏振图像配准方法[J]. 液晶与显示,2019,34(5):530-536.doi:10.3788/YJYXS20193405.0530

    WAN Y L, WANG J L, ZHANG N,et al. Polarized image registration method based on phase correlation and sub-graph[J].Chinese Journal of Liquid Crystals and Displays, 2019, 34(5): 530-536. (in Chinese)doi:10.3788/YJYXS20193405.0530
    [25] PERONA P, MALIK J. Scale-space and edge detection using anisotropic diffusion[J].IEEE Transactions on Pattern Analysis and Machine Intelligence, 1990, 12(7): 629-639.doi:10.1109/34.56205
    [26] 王亚强, 陈波. 一种改进的各向异性扩散超声图像去噪算法[J]. 液晶与显示,2015,30(2):310-316.doi:10.3788/YJYXS20153002.0310

    WANG Y Q, CHEN B. Improved anisotropic diffusion ultrasound image denoising algorithm[J].Chinese Journal of Liquid Crystals and Displays, 2015, 30(2): 310-316. (in Chinese)doi:10.3788/YJYXS20153002.0310
    [27] PUVANATHASAN P, BIZHEVA K. Interval type-II fuzzy anisotropic diffusion algorithm for speckle noise reduction in optical coherence tomography images[J].Optics Express, 2009, 17(2): 733-746.doi:10.1364/OE.17.000733
  • 加载中
图(8)/ 表(1)
计量
  • 文章访问数:1205
  • HTML全文浏览量:318
  • PDF下载量:108
  • 被引次数:0
出版历程
  • 收稿日期:2020-07-23
  • 修回日期:2020-09-03
  • 网络出版日期:2021-03-05
  • 刊出日期:2021-03-23

目录

    /

      返回文章
      返回
        Baidu
        map