CN109758119B - 一种肿瘤微血管图像处理装置 - Google Patents

一种肿瘤微血管图像处理装置 Download PDF

Info

Publication number
CN109758119B
CN109758119B CN201910080990.7A CN201910080990A CN109758119B CN 109758119 B CN109758119 B CN 109758119B CN 201910080990 A CN201910080990 A CN 201910080990A CN 109758119 B CN109758119 B CN 109758119B
Authority
CN
China
Prior art keywords
image processing
imaging
fiber coupler
optical fiber
tumor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910080990.7A
Other languages
English (en)
Other versions
CN109758119A (zh
Inventor
秦嘉
蔡佳龙
易遥
安林
叶新荣
邓卓健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Weiren Medical Technology Co ltd
Original Assignee
Guangdong Weiren Medical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Weiren Medical Technology Co ltd filed Critical Guangdong Weiren Medical Technology Co ltd
Priority to CN201910080990.7A priority Critical patent/CN109758119B/zh
Publication of CN109758119A publication Critical patent/CN109758119A/zh
Application granted granted Critical
Publication of CN109758119B publication Critical patent/CN109758119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种肿瘤微血管图像处理装置,包括肿瘤微血管成像***和图像处理模块,所述图像处理模块采用如下方法对获取的血管图像进行处理,包括步骤:S1,对结构信号进行大血管***影补偿;S2,利用Hessian滤波进行三维管状形态物体增强。本发明的图像处理模块对成像观察到的局部血管几何特征对血流信息进行针对性地增强和抑制,使得真实血流部分得到保留,从而可以实现活体提取与量化肿瘤表层微血管的参数,从而可以指导早期肿瘤的筛查及诊断、***光学活检和术中指导切除范围等。

Description

一种肿瘤微血管图像处理装置
技术领域
本发明属于图像处理装置,具体涉及一种肿瘤微血管的图像处理装置。
背景技术
血管生成在肿瘤生长和扩散过程中至关重要,肿瘤细胞诱导血管生成往往发生在肿瘤形成初期到中期阶段。因此,肿瘤血管形态学的可视化有助于提高诊断的准确性。
当前的肿瘤成像技术,如X线计算机断层摄影(CT),正电子发射断层摄(PET)和磁共振成像(MRI)等,往往分辨率低,达不到可靠检测<2mm的肿瘤沉积物的作用;同时神经外科医生或许也有这样的困扰:在清除肿瘤时,该如何精确地界定肿瘤的切除区域?清除不彻底,肿瘤会复发,切除了正常的生理组织,可能给病人造成终身伤害。因此需要一种先进的技术,即时、非侵入性、精准评估肿瘤的发生与分布情况。
扫频源光学相干层析成像***(Swept-source optical coherence tomography,SS-OCT),集光学、超灵敏探测及计算机图像处理技术于一体,是一种可以快速提供拥有微米级分辨率的活体组织形态的断面图像的辅助诊疗工具,对人体无创无损。通过与内镜、细针、导管、腹腔镜等装置结合,在肿瘤外科中拥有巨大的成像潜力。
发明内容
本发明的目的是,提供一种肿瘤微血管的图像处理装置,从而便于获得准确的肿瘤微血管图像信息,为临床诊断和治疗提供更好的支持。
为了达到上述技术目的,本发明的技术方案如下:
一种肿瘤微血管图像处理装置,包括肿瘤微血管成像***和图像处理模块,所述图像处理模块用于接收和处理所述肿瘤微血管成像***产生的图像信号并进行处理,所述图像处理模块采用如下方法对获取的血管图像进行处理,包括步骤:
S1,对结构信号进行大血管***影补偿;
S2,利用Hessian滤波进行三维管状形态物体增强。
本发明提供的处理装置具备如下有益效果:
本发明的图像处理模块对成像观察到的局部血管几何特征对血流信息进行针对性地增强和抑制,通过该处理,血管下方非血流伪影信号得到抑制而对应的真实血流部分得到保留,从而可以实现活体提取与量化肿瘤表层微血管的参数,如肿瘤微血管半径、密度,微血管血流量与血氧含量等,揭示在健康人和癌症患者肿瘤形态学和血管疾病学的差异,从而可以指导早期肿瘤的筛查及诊断、***光学活检和术中指导切除范围等。
优选地,所述图像处理模块的图像处理速率为40-215fps。
优选地,所述处理装置的灵敏度>120dB,成像***的成像范围2mm×2mm,成像深度3mm。
优选地,所述成像***的横向分辨率为16.0μm,轴向分辨率为6μm。
作为优选方式,所述肿瘤微血管成像***为SS-OTC成像***。
进一步地,所述SS-OCT成像***包括扫频光源,第一光纤耦合器,主成像干涉***,标定干涉***,探测臂;其中主成像干涉***包括第二光纤耦合器、参考臂、样品臂、探测臂、第三光纤耦合器及第一平衡探测器;标定干涉***包括马赫-曾德尔干涉仪、第四光纤耦合器及第二平衡探测器;
所述扫频光源提供的光束,经光纤传输入第一光纤耦合器后,分为两部分,分别进入主成像干涉***和标定干涉***;
进入主成像干涉***的光束首先被第二光纤耦合器分为两部分,分别进入参考臂和样品臂,进入参考臂的光经反射后原路返回;进入样品臂的光,最后到达探测臂对样品进行扫描,然后原路返回;参考臂和样品臂向后返回的光进入第三光纤耦合器,最后经第一平衡探测器进行光电转化形成成像干涉信号;
进入标定干涉***的光束,先进入马赫-曾德尔干涉仪,然后进入第四光纤耦合器,最后经第二平衡探测器进行光电转化形成标定干涉信号;
所述干涉信号和标定干涉信号被所述图像处理模块接收成像。
优选地,所述扫频光源,波长为1310nm,带宽长度为110nm,输出功率15mW,扫描速度为40-220kHz(40000-220000A-lines/s)。
优选地,所述第一光纤耦合器为95/5光纤耦合器,95%的光进入主成像干涉***,5%的光进入标定干涉***。
优选地,所述第二光纤耦合器为75/25光纤耦合器,75%的光进入样品臂,25%的光进入参考臂。
优选地,所述探测臂内设置有x-y扫描式振镜,直径12mm。
附图说明
图1为本发明实施例1提供的SS-OTC成像***的肿瘤微血管图像处理装置结构示意图。
具体实施方式
下面结合附图和具体实施例详细说明本发明的技术方案。
血管生成在肿瘤生长和扩散过程中至关重要,肿瘤细胞诱导血管生成往往发生在肿瘤形成初期到中期阶段。因此,肿瘤血管形态学的可视化有助于提高诊断的准确性。现有的肿瘤血管成像技术,由于各种各样的限制,采集到的图像总存在一定的噪音,因此需要进行处理,才能使信息更准确。
本发明旨在提供一种肿瘤微血管图像处理装置,包括肿瘤微血管成像***和图像处理模块,所述图像处理模块用于接收和处理所述肿瘤微血管成像***产生的图像信号并进行处理,所述图像处理模块采用如下方法对获取的血管图像进行处理,包括步骤:
S1,对结构信号进行大血管***影补偿;
通过对结构信号进行大血管***影补偿可以恢复大血管下结构信号跟周围组织的一致性,具体可以通过提高横断面的成像分辨率和改进***的成像深度进行补偿。
S2,利用Hessian滤波进行三维管状形态物体增强。
Hessian矩阵的特征值能够很好地描述眼底图像的血管信息,利用Hessian滤波进行三维管状形态物体增强,可以得到每个体素点对应的血管指数V(r)(0~1)。V(r)接近于1表示体素点处属于血管的可能性较大,接近于0表示血管的可能性小。对归一化的血流信号进行对应的加权:
Fnew(r)=Fnorm(r)Vl(r)
其中l是一个非线性指数加权因子。通过该方法,血管下方非血流伪影信号得到抑制而对应的真实血流部分得到保留,从而使图像更准确。
本发明的图像处理模块对成像观察到的局部血管几何特征对血流信息进行针对性地增强和抑制,通过该处理,血管下方非血流伪影信号得到抑制而对应的真实血流部分得到保留,从而可以实现活体提取与量化肿瘤表层微血管的参数,如肿瘤微血管半径、密度,微血管血流量与血氧含量等,揭示在健康人和癌症患者肿瘤形态学和血管疾病学的差异,从而可以指导早期肿瘤的筛查及诊断、***光学活检手术中指导切除范围等。
优选地,所述图像处理模块的图像处理速率为40-215fps。
优选地,所述处理装置的灵敏度>120dB,成像***的成像范围2mm×2mm,成像深度3mm。
优选地,所述成像***的横向分辨率为16.0μm,轴向分辨率为6μm。
本发明的微肿瘤成像***可以采用X线计算机断层摄影(CT),正电子发射断层摄(PET)和磁共振成像(MRI),或者光线相干层成像技术(OCT)等,优选地,采用OCT成像***,成像过程友好,无创无损,耗时短,成像分辨率高,而且费用便宜,OCT成像又包括光谱域OCT(SD-OCT)、频域OCT(FD-OCT)或扫频源OCT(SS-OCT),鉴于SS-OCT具有更长波长、更快速度的扫频光源以及能够精确逐点探测的平衡探测器,因此进一步优选地,所述肿瘤微血管成像***为SS-OTC成像***。
下面结合实施例说明本发明优选的实施方式。
实施例1
参照图1所示,一种肿瘤微血管图像处理装置,包括肿瘤微血管成像***和图像处理模块12。
其中肿瘤微血管成像***为SS-OTC成像***;包括:
扫频光源1,其波长为1310nm,带宽长度为110nm,输出功率15mW,扫描速度为40-220kHz(40000-220000A-lines/s);
第一光纤耦合器2,第一光纤耦合器2为95/5光纤耦合器;
主成像干涉***,主成像干涉***包括第二光纤耦合器3、参考臂、样品臂、探测臂4、第三光纤耦合器7及第一平衡探测器8;其中第二光纤耦合器3为75/25光纤耦合器,参考臂内设有一环形器13,聚焦透镜5和平面镜6;样品臂内设有一环形器13,探测臂4内设置有x-y扫描式振镜,直径12mm。
标定干涉***包括马赫-曾德尔干涉仪9、第四光纤耦合器10及第二平衡探测器11。
扫频光源1提供的光束,经光纤传输入第一光纤耦合器2后,分为两部分,95%的光进入主成像干涉***,5%的光进入标定干涉***;进入主成像干涉***的光束首先被第二光纤耦合器3分为两部分,75%的光进入样品臂,25%的光进入参考臂,进入参考臂的光经环形器13后通过透镜5聚焦到平面镜6反射后原路返回;进入样品臂的光,经过环形器13后最后到达探测臂4对样品进行扫描,然后原路返回;参考臂和样品臂向后返回的光分别经过对应的环形器13导向,然后汇入第三光纤耦合器7,分为两束等强度的干涉光谱,最后经第一平衡探测器8进行光电转化形成成像干涉信号;
进入标定干涉***的光束,先进入马赫-曾德尔干涉仪9,然后进入第四光纤耦合器10,分为两束等强度的干涉光谱,最后经第二平衡探测器11进行光电转化形成标定干涉信号;
成像干涉信号和标定干涉信号被图像处理模块12接收成像。图像处理模块对获取的血管图像进行处理,包括步骤:
S1,对结构信号进行大血管***影补偿;
通过对结构信号进行大血管***影补偿可以恢复大血管下结构信号跟周围组织的一致性,具体可以通过提高横断面的成像分辨率和改进***的成像深度进行补偿。
S2,利用Hessian滤波进行三维管状形态物体增强。
利用Hessian滤波进行三维管状形态物体增强,可以得到每个体素点对应的血管指数V(r)(0~1)。V(r)接近于1表示体素点处属于血管的可能性较大,接近于0表示血管的可能性小。对归一化的血流信号进行对应的加权:
Fnew(r)=Fnorm(r)Vl(r)
其中l是一个非线性指数加权因子。通过该方法,血管下方非血流伪影信号得到抑制而对应的真实血流部分得到保留,从而使图像更准确。
图像处理模块12可以包括图像采集卡和图像处理卡,可以选用NI PCIe-1433型号的产品。
实施例1提供的图像处理装置具备如下有益效果:
1、与光谱域OCT(SD-OCT)、频域OCT(FD-OCT),SS-OCT具有更长波长、更快速度的扫频光源以及能够精确逐点探测的平衡探测器,因此它具有更快的成像速度、更深的成像深度、更高的信噪比、更强的灵敏度。
2、SS-OCT***成像过程友好,无创无损,耗时短,成像分辨率高,且费用便宜。
3、对肿瘤微血管实现成像和量化,可以指导早期肿瘤的筛查及诊断、***光学活检和术中指导切除范围。
本发明图像处置装置与其他不同的成像技术结合起来,还可以实现检测更为细小的微血管网络,从而将发明应用到其他人体如胃癌、乳腺癌、***癌等各类肿瘤的检测。
根据上述说明书的揭示和教导,本发明所属领域的技术人员还可以对上述实施方式进行变更和修改。因此,本发明并不局限于上面揭示和描述的具体实施方式,对本发明的一些修改和变更也应当落入本发明的权利要求的保护范围内。此外,尽管本说明书中使用了一些特定的术语,但这些术语只是为了方便说明,并不对本发明构成任何限制。

Claims (7)

1.一种肿瘤微血管图像处理装置,包括肿瘤微血管成像***和图像处理模块,所述图像处理模块用于接收和处理所述肿瘤微血管成像***产生的图像信号并进行处理,其特征在于,所述图像处理模块采用如下方法对获取的血管图像进行处理,包括步骤:
S1,对结构信号进行大血管***影补偿;
通过对结构信号进行大血管***影补偿可以恢复大血管下结构信号跟周围组织的一致性,具体可以通过提高横断面的成像分辨率和改进***的成像深度进行补偿;
S2,利用Hessian滤波进行三维管状形态物体增强;
利用Hessian滤波进行三维管状形态物体增强,可以得到每个体素点对应的血管指数V(r);V(r)接近于1表示体素点处属于血管的可能性较大,接近于0表示血管的可能性小;
所述肿瘤微血管成像***为SS-OTC成像***;
所述SS-OCT成像***包括扫频光源,第一光纤耦合器,主成像干涉***,标定干涉***,探测臂;
其中,主成像干涉***包括第二光纤耦合器、参考臂、样品臂、探测臂、第三光纤耦合器及第一平衡探测器;
标定干涉***包括马赫-曾德尔干涉仪、第四光纤耦合器及第二平衡探测器;
所述扫频光源提供的光束,经光纤传输入第一光纤耦合器后,分为两部分,分别进入主成像干涉***和标定干涉***;
进入主成像干涉***的光束首先被第二光纤耦合器分为两部分,分别进入参考臂和样品臂,进入参考臂的光经反射后原路返回;进入样品臂的光,最后到达探测臂对样品进行扫描,然后原路返回;参考臂和样品臂向后返回的光进入第三光纤耦合器,最后经第一平衡探测器进行光电转化形成成像干涉信号;
进入标定干涉***的光束,先进入马赫-曾德尔干涉仪,然后进入第四光纤耦合器,最后经第二平衡探测器进行光电转化形成标定干涉信号;
所述干涉信号和标定干涉信号被所述图像处理模块接收成像;
所述第一光纤耦合器为95/5光纤耦合器,95%的光进入主成像干涉***,5%的光进入标定干涉***。
2.根据权利要求 1 所述的肿瘤微血管图像处理装置,其特征在于, 所述第二光纤耦合器为75/25光纤耦合器,75%的光进入样品臂,25%的光进入参考臂。
3.根据权利要求1所述的肿瘤微血管图像处理装置,其特征在于,所述图像处理模块的图像处理速率为40-215fps。
4.根据权利要求1所述的肿瘤微血管图像处理装置,其特征在于,所述处理装置的灵敏度 >120 dB,成像***的成像范围2 mm×2 mm,成像深度3 mm。
5.根据权利要求1所述的肿瘤微血管图像处理装置,其特征在于,所述成像***的横向分辨率为16.0μm,轴向分辨率为6 μm。
6.根据权利要求1所述的肿瘤微血管图像处理装置,其特征在于,所述扫频光源波长为1310 nm,带宽长度为110 nm,输出功率15 mW,扫描速度为40-220kHz。
7.根据权利要求1所述的肿瘤微血管图像处理装置,其特征在于,所述探测臂内设置有x-y扫描式振镜,直径12 mm。
CN201910080990.7A 2019-01-28 2019-01-28 一种肿瘤微血管图像处理装置 Active CN109758119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910080990.7A CN109758119B (zh) 2019-01-28 2019-01-28 一种肿瘤微血管图像处理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910080990.7A CN109758119B (zh) 2019-01-28 2019-01-28 一种肿瘤微血管图像处理装置

Publications (2)

Publication Number Publication Date
CN109758119A CN109758119A (zh) 2019-05-17
CN109758119B true CN109758119B (zh) 2022-03-08

Family

ID=66454448

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910080990.7A Active CN109758119B (zh) 2019-01-28 2019-01-28 一种肿瘤微血管图像处理装置

Country Status (1)

Country Link
CN (1) CN109758119B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110403576B (zh) * 2019-08-01 2024-01-30 中国医学科学院北京协和医院 三维光声成像在乳腺肿瘤评分***中的应用
CN110793941A (zh) * 2019-10-10 2020-02-14 成都贝瑞光电科技股份有限公司 一种智能分辨激光光学层析成像方法的***

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN200970231Y (zh) * 2006-11-27 2007-11-07 浙江大学 一种扩展光学相干层析成像动态范围的***
CN101803908A (zh) * 2010-03-01 2010-08-18 浙江大学 基于色散调制的无镜像光学频域成像***及方法
CN102151121B (zh) * 2011-01-28 2012-08-15 浙江大学 基于干涉光谱相位信息的光谱标定方法及***
CN102657519B (zh) * 2012-05-11 2013-11-20 浙江大学 基于oct的大动态范围流速的高灵敏度测量***及方法
CN104239874B (zh) * 2014-09-29 2017-11-03 青岛海信医疗设备股份有限公司 一种器官血管识别方法及装置
CN105989598A (zh) * 2015-02-13 2016-10-05 中国科学院沈阳自动化研究所 基于局部强化主动轮廓模型的眼底图像血管分割方法
CN105030201B (zh) * 2015-07-27 2017-10-13 浙江大学 基于静态区域信息的扫频oct数字相位矫正方法与***
CN108245130B (zh) * 2016-12-28 2021-06-04 南京理工大学 一种光学相干断层血管造影装置及方法

Also Published As

Publication number Publication date
CN109758119A (zh) 2019-05-17

Similar Documents

Publication Publication Date Title
Hu et al. Photoacoustic imaging and characterization of the microvasculature
Pogue et al. Quantitative hemoglobin tomography with diffuse near-infrared spectroscopy: pilot results in the breast
Jesser et al. High resolution imaging of transitional cell carcinoma with optical coherence tomography: feasibility for the evaluation of bladder pathology.
JP7058264B2 (ja) 血管内光干渉断層撮影法を用いた冠動脈プラークの自動特徴分析およびリスク評価のためのシステムおよび方法
Choi et al. Improved microcirculation imaging of human skin in vivo using optical microangiography with a correlation mapping mask
CN100477955C (zh) 用于光电检测肿瘤的装置
US6795195B1 (en) System and method for tomographic imaging of dynamic properties of a scattering medium
US9858668B2 (en) Guidewire artifact removal in images
JP2016507272A (ja) 機能的利得の測定技術及び表示
CN105748040B (zh) 立体结构功能成像***
Jain et al. Full-field optical coherence tomography for the analysis of fresh unstained human lobectomy specimens
CN205215191U (zh) 一种血管内成像的光学相干层析装置
Turchin Methods of biomedical optical imaging: from subcellular structures to tissues and organs
Jung et al. Optical coherence tomography for rapid tissue screening and directed histological sectioning
CN109758119B (zh) 一种肿瘤微血管图像处理装置
Lin et al. Photoacoustic imaging
CN105411530A (zh) 新型光学相干层析装置
Patterson et al. Optoacoustic characterization of prostate cancer in an in vivo transgenic murine model
CN114468960A (zh) 集光学、光声、oct、荧光多模态前向成像一体化的内窥镜装置及其成像方法
WO2018161384A1 (zh) 一种呼吸道oct***
Fakhoury et al. Photoacoustic imaging for cutaneous melanoma assessment: a comprehensive review
Ma et al. Multiscale confocal photoacoustic dermoscopy to evaluate skin health
Kiseleva et al. Cross-polarization optical coherent tomography in comparative in vivo and ex vivo studies of optical properties of normal and tumorous brain tissues
Pantanowitz et al. High‐resolution imaging of the thyroid gland using optical coherence tomography
CN210130817U (zh) 一种基于ss-oct肿瘤微血管成像装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A tumor microvascular image processing device

Effective date of registration: 20230227

Granted publication date: 20220308

Pledgee: Foshan rural commercial bank Limited by Share Ltd.

Pledgor: GUANGDONG WEIREN MEDICAL TECHNOLOGY Co.,Ltd.

Registration number: Y2023980033570

PE01 Entry into force of the registration of the contract for pledge of patent right