CN101999885B - 一种自动分离动静脉血管的内源光学成像方法 - Google Patents
一种自动分离动静脉血管的内源光学成像方法 Download PDFInfo
- Publication number
- CN101999885B CN101999885B CN201010598859A CN201010598859A CN101999885B CN 101999885 B CN101999885 B CN 101999885B CN 201010598859 A CN201010598859 A CN 201010598859A CN 201010598859 A CN201010598859 A CN 201010598859A CN 101999885 B CN101999885 B CN 101999885B
- Authority
- CN
- China
- Prior art keywords
- frequency
- image
- pixel
- optical imaging
- arteries
- 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
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种自动分离动静脉血管的内源光学成像方法,目的是提供一种有效的自动分离生物组织动静脉血管的内源光学成像方法,技术方案是采用窄带准单色光照明,对采集到的图像序列进行频谱分析,利用图像序列的低频振荡信号的幅值分布特征,通过频谱变换,局部自适应对比度增强和自动阈值分割等操作实现了动脉血管和静脉血管的识别与自动分割。采用本发明可实现动静脉血管的自动分离,能凸显并分离很多小动脉和小静脉的血管结构,避免了使用多种中心波长的光照明时成像装置设计的复杂性。
Description
技术领域:
本发明属于生物医学成像方法,尤其指采用窄带准单色光照明,利用内源生理信号谱值特征的动静脉血管自动分离方法。
背景技术:
脑皮层动静脉血管的自动分离对于生物医学基础研究和血管相关疾病的临床诊断治疗具有重要的现实意义。动脉肿瘤,动脉硬化等血管疾病日益威胁着人们的生命安全,通过血管分析方法能够及时可靠的诊断和预防此类疾病。目前,对动静脉血管进行自动分离的方法的报道较为有限。根据动脉血和静脉血含氧饱和度的不同而具有的光谱吸收特性差异,可通过不同中心波长的准单色光检测目标血管的血氧饱和度的值,实现动脉血管和静脉血管的自动分离(Narasimha-Iyer,H等,2007.Automatic identification of retinal arteries and veins fromdual-wavelength images using structural and functional features(双光谱图像中的结构特征和功能特征用于视网膜动脉和静脉的自动识别),IEEE Transactions on Biomedical Engineering,54(8),1427-1435),附图1给出了上述方法的流程图。但是上述方法对细小血管的分离效果很差,而且受限于预先提取血管网络结构,在现有的血管网络提取方法下,无法分离出小动脉和小静脉。以MRA成像(磁共振造影成像)为技术手段的动静脉分离方法,主要依赖于动脉和静脉血管中血液的流向相反特性(Svensson,J等,2002.Separation of arteries and veinsusing flow-induced phase effects in contrast-enhanced MRA of the lower extremities(下肢MRA成像中基于血流相位特征的动静脉分离),Magnetic Resonance Imaging,20,49-57),但该特性仅存在于人体的四肢血管和部分心血管系统中,在复杂的大脑血管网络中是不存在的,因此无法应用于脑皮层中动静脉的分离。生物组织中内源信号包含的频率值很丰富,在生理信号的频谱分布中,低于1Hz的信号称为低频振荡信号,这些振荡信号以0.1Hz频率为中心,包括超低频振荡信号(VLFO),低频振荡信号(LFO),其中0.1Hz低频振荡信号的存在具有普遍性,是生物体内源振荡信号的主要成分之一;通过光学成像系统采集的光学信号是根据不同的血液成分对不同波长的光的吸收率和反射率不同来反映动脉血和静脉血的参数变化情况,如血流,血容,血氧浓度等,从而反映血管的动态变化过程。但目前还没有利用生物组织中低频振荡信号的谱值特征进行动静脉血管的自动分离的公开报导。
发明内容:
本发明要解决的技术问题是提供一种有效的自动分离生物组织动静脉血管的内源光学成像方法,采用窄带准单色光照明,对采集到的图像序列进行频谱分析,利用图像序列的低频振荡信号的幅值分布特征,实现动脉血管和静脉血管的识别与自动分割。
为解决上述技术问题,本发明动静脉自动分离的内源光学成像方法包括以下步骤:
第一步,将窄带准单色光照射到被测对象上,用CCD或CMOS相机通过光学成像系统,以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的光学图像,每帧图像的曝光时间不高于100ms,帧间隔时间ΔT不高于500ms,采集的帧数N≥20;
第二步,设单幅图像的尺寸为IR×C,其中R为图像矩阵的行数,C为图像矩阵的列数,则该图像共有R×C个像素,相应的图像时间序列中包含R×C个像素时间序列;平均像素时间序列中第n个数据点的平均灰度值按公式一有:
第三步,利用快速傅里叶变换算法(A.V.奥本海姆,R.W.谢弗(美).离散时间信号处理.黄建国,刘树樘 译.北京:科学出版社,1998)将第二步所得的平均像素时间序列变换到频谱域,确定频率区间(0.05Hz,1Hz)中具有最大幅值的频率值fm;fm通常位于0.1Hz附近,称为0.1-Hz低频振荡。
第四步,利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,获得fm频率处的幅值Am(r,c),称为特征频率幅值,其中r为像素所在位置的行,c为像素所在位置的列;以每个像素时间序列对应的特征频率幅值Am(r,c)为灰度,组成二维的特征频率幅值图Am;
第五步,利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,计算频率区间[1Hz,1/(2·ΔT)]中所有离散频率点的谱功率总和Wg(r,c),其中各频率点的谱功率为相应频率幅值的平方;计算特征频率比值Ra(r,c),Ra(r,c)为各像素时间序列对应的特征频率fm的谱值(Am(r,c))2与Wg(r,c)的比值,Ra(r,c)=(Am(r,c))2/Wg(r,c);以各像素时间序列对应的特征频率比值Ra(r,c)为灰度,组成二维的特征频率比值图Ra;
第六步,为提高动脉血管和静脉血管的灰度分布差异,利用局部对比度增强方法(C.Sinthanayothin,F.J.Boyce,H.L.Cook,1999.Automated localisation of the optic disk,fovea,andretinal blood vessels from digital colour fundus images(眼底数字彩色图像中视盘、中央凹和视网膜血管的自动定位).British Journal of Ophthalmology 83,902-910)对由第四步和第五步得到的特征频率幅值图Am和特征频率比值图Ra做局部对比度增强,具体操作如下:
步骤6.1以图像中第r行,第c列处的像素为中心,选取像素大小为49×49的空间滑动窗W,按公式二得到中心像素的新灰度值p(r,c),
其中,pmin为待操作图像所有像素灰度值的最小值,pmax为待操作图像所有像素灰度值的最大值,Ψw(p)按公式三计算,
其中和σw分别为W滑动窗所覆盖的所有像素的灰度平均值和方差;
步骤6.2遍历图像Am和Ra中所有的像素,将每个像素的灰度值按步骤6.1的方法进行更新,得到局部对比度增强后的特征频率幅值图IAm和特征频率比值图IRa。
第七步,使用阈值分割方法(章毓晋著.图像分割.北京:科学出版社,2001)将第六步中所得局部对比度增强后的特征频率幅值图IAm作二值化处理,得到二维动脉血管结构图像Iartery,其中动脉血管区域像素值为1,其他区域像素值为0;
第八步,使用阈值分割方法将第六步中所得图像IRa作二值化处理,得到二维静脉血管结构图像Ivein,其中静脉血管区域像素值为1,其他区域像素值为0;
第九步,为体现动静脉分离后的效果,可以将图像Iartery和图像Ivein按如下方式构建新的RGB图像CRGB,动脉区域和静脉区域分别用R(红色)、G(绿色)、B(蓝色)三基色中的一种来表示。设置动脉区域的颜色值为红色:(1,0,0),静脉区域的颜色值为蓝色(0,0,1),其他区域的颜色值为白色(1,1,1);从而绘制出一幅以不同颜色标记动静脉血管的RGB图像。
本发明以窄带准单色光照明,用CCD或CMOS相机通过内源光学成像系统连续采集若干帧光学图像,对图像序列做频谱分析,获取0.1-Hz振荡的特征频率幅值图和比值图,以及局部对比度增强等操作,实现生物组织中动静脉血管的自动分离。与其他现有的方法相比,采用本发明可以达到以下技术效果:
1.采用一种准单色光照明,避免了使用多种中心波长的光照明时成像装置设计的复杂性。
2.利用低频振荡信号的幅值分布特性,实现了动静脉血管的自动分离。
3.由于低频振荡信号能够反应动脉和静脉的功能性差异,因此采用本发明能凸显并分离很多小动脉和小静脉的血管结构,可避免造影剂注射的MRA成像方法区分动静脉血管时产生的副作用,而这些小动脉和小静脉的血管结构由于相对于背景软组织的对比度很低,仅依靠光学图像的灰度特征是很难提取和分离的。
4.本发明有助于对生物组织中二维血管的形态和动静脉血管中血流分布,及血管中血液动力学变化进行实时、动态和高时空分辨率的监测。
附图说明
图1为已报导的视网膜动静脉分离方法流程图(Narasimha-Iyer,H等,2007.Automaticidentification of retinal arteries and veins from dual-wavelength images using structural andfunctional features(双光谱图像中的结构特征和功能特征用于视网膜动脉和静脉的自动识别),IEEE Transactions on Biomedical Engineering,54(8),1427-1435)。
图2为本发明总体流程图。
图3为采用已报导的内源光学成像系统(Optical Imaging Inc.,Germantown,NY,USA)采集去除颅骨的大鼠顶叶脑皮层图像并用本发明所提出的自动分离动静脉血管的方法进行处理的结果。(a)脑皮层灰度图像;(b)平均时间序列的频谱分布;(c)各像素时间序列的fm幅值分布图Am;(d)各像素时间序列的fm功率比值图Ra。(e)用不同颜色表示的动静脉血管分离结果图。
具体实施方式
图2为本发明总体流程图,本发明包括以下步骤:
1、将窄带准单色光照射到被测对象上,用CCD或CMOS相机通过光学成像系统,以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的光学图像;
4、利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,计算每个像素时间序列在特征频率fm处的幅值,组成特征频率幅值图Am;
5、利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,计算频率区间[1Hz,1/(2·ΔT)]中所有离散频率点的谱功率总和Wg(r,c),其中各频率点的谱功率为相应频率幅值的平方;计算特征频率比值Ra(r,c),Ra(r,c)为各像素序列对应的特征频率fm的谱值(Am(r,c))2与Wg(r,c)的比值,Ra(r,c)=(Am(r,c))2/Wg(r,c);以各像素序列对应的特征频率比值Ra(r,c)为灰度,组成二维的特征频率比值图Ra;
6、利用局部对比度增强方法对由第4步和第5步得到的特征频率幅值图Am和特征频率比值图Ra做局部对比度增强(公式二和公式三),分别得到增强图像IAm和IRa;
7、使用阈值分割方法将第6步中所得图像IAm作二值化处理,得到二维动脉血管结构图像Iartery;
8、使用阈值分割方法将第6步中所得图像IRa作二值化处理,得到二维静脉血管结构图像Ivein;
9、将得到的动脉和静脉以不同的颜色区别标记,以显示动静脉分离的效果。
动物实验:
如图3所示,图3-(a)是(546±10)nm准单色光照射下的单帧光学灰度图像;图中所示标尺长度为0.5mm。我们采用购买自Optical Imaging公司(Germantown,NY,USA)的光学成像与采集装置,光学采集装置置于物理隔震台上。实验对象采用Sprague-Dawley大鼠,固定于立体定位架上,以卤素灯为光源,经滤光片(546±10nm)滤波后斜入射到去除头骨后的大鼠顶叶皮层上。经大鼠顶叶皮层反射后的光线由光学成像系统成像,并通过图像采集卡与计算机通信,由计算机存储采集到的图像序列。曝光时间30ms,帧时间间隔66.67ms,连续采集1500帧光学灰度图像。冷光源经滤光片滤光得到的绿色光束(波长为546±10nm)通过光纤照射到待采集区域。经过血管区域反射的光线通过CCD镜头成像并由采集系统完成血管图像数据的采集和存储。图像序列采集频率为15Hz,采集过程是在静息(无刺激)状态下完成的。将采集到的图像时间序列变换到频域,构建特征频率幅值图和比值图,局部自适应对比度增强等操作,构建一幅由不同颜色标记动静脉血管的伪彩色图像。图3-(b)是按本发明第二步和第三步对1500帧光学灰度图像序列中的每一帧图像做灰度平均,得到有1500个采样点的平均灰度时间序列,再用快速傅里叶变换算法将平均灰度时间序列变换到频域,搜索0.1Hz频率附近的峰值幅度,由此确定具有最大幅值的特征频率的位置fm;图3-(c)是按本发明第四步对每个像素位置处对应的灰度时间序列用FFT方法做频谱变换,计算fm处的幅值作为该像素位置处的灰度值,组成的特征频率幅值图并做局部对比度增强(第六步)后的结果;图3-(d)是按本发明第五步计算每个像素位置处对应的特征频率fm的幅值与≥1Hz频率区间内谱功率之和的比值图,并做局部对比度增强(第六步)后的结果。按本发明第八步和第九步对图3-(c)和图3-(d)分别做二值化处理,前者得到二维动脉血管结构图,后者得到二维静脉血管结构图。由图3-(c)和图3-(d)可见,0.1-Hz振荡信号在动脉区域的幅值要高于其在静脉和皮层其他区域的幅值,而静脉区域的0.1-Hz振荡信号的幅值是最低的,因此,对图3-(b)和图3-(c)通过自动阈值分割方法作二值化处理,可分别得到动脉和静脉血管结构。按发明内容所述第九步,将分割出来的动脉区域用红色表示,静脉区域用蓝色表示,其他皮层区域用白色表示,构建一幅RGB伪彩色图像,如附图3-(e)所示。由附图3-(e)可知,所分割出的动脉血管和静脉血管的走向与其解剖学特征相吻合,说明动脉血管和静脉血管得到了很好的区分,证明了本发明所提出的自动分离动静脉血管的内源光学成像方法的有效性。
Claims (2)
1.一种自动分离动静脉血管的内源光学成像方法,其特征在于包括以下步骤:
第一步,将窄带准单色光照射到被测对象上,用CCD或CMOS相机通过光学成像系统,
以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的光学图像,每帧图像的曝光时间不高于100ms,帧间隔时间ΔT不高于500ms,采集的帧数N≥20;
其中pn,i为第n幅图像的第i个像素的灰度值,为N帧图像所有像素的平均灰度值,R为图像矩阵的行数,C为图像矩阵的列数,单幅图像的尺寸为IR×C,该图像共有R×C个像素,相应的图像时间序列中包含R×C个像素时间序列;
第四步,利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,获得fm频率处的特征频率幅值Am(r,c),其中r为像素所在位置的行,c为像素所在位置的列;以每个像素时间序列对应的特征频率幅值Am(r,c)为灰度,组成二维的特征频率幅值图Am;
第五步,利用快速傅里叶变换算法将每个像素时间序列变换到频谱域,计算频率区间[1Hz,1/(2·ΔT)]中所有离散频率点的谱功率总和Wg(r,c),其中各频率点的谱功率为相应频率幅值的平方;计算特征频率比值Ra(r,c),Ra(r,c)为各像素时间序列对应的特征频率fm的谱值(Am(r,c))2与Wg(r,c)的比值,Ra(r,c)=(Am(r,c))2/Wg(r,c);以各像素时间序列对应的特征频率比值Ra(r,c)为灰度,组成二维的特征频率比值图Ra;
第六步,利用局部对比度增强方法对由第四步和第五步得到的特征频率幅值图Am和特征频率比值图Ra做局部对比度增强,具体操作如下:
步骤6.1以图像中第r行,第c列处的像素为中心,选取像素大小为49×49的空间滑动窗W,按公式二得到中心像素的新灰度值p(r,c),
其中,pmin为待操作图像所有像素灰度值的最小值,pmax为待操作图像所有像素灰度值的最大值,Ψw(p)按公式三计算,
步骤6.2遍历图像Am和Ra中所有的像素,将每个像素的灰度值按步骤6.1的方法进行更新,得到局部对比度增强后的特征频率幅值图IAm和特征频率比值图IRa;
第七步,使用阈值分割方法将第六步中所得局部对比度增强后的特征频率幅值图IAm作二值化处理,得到二维动脉血管结构图像Iartery,其中动脉血管区域像素值为1,其他区域像素值为0;
第八步,使用阈值分割方法将第六步中所得图像IRa作二值化处理,得到二维静脉血管结构图像Ivein,其中静脉血管区域像素值为1,其他区域像素值为0。
2.如权利要求1所述的一种自动分离动静脉血管的内源光学成像方法,其特征在于将图像Iartery和图像Ivein按如下方式构建新的RGB图像CRGB,动脉区域和静脉区域分别用红色R、绿色G、蓝色B三基色中的一种来表示,设置动脉区域的颜色值为红色:(1,0,0),静脉区域的颜色值为蓝色(0,0,1),其他区域的颜色值为白色(1,1,1);从而绘制出一幅以不同颜色标记动静脉血管的RGB图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010598859A CN101999885B (zh) | 2010-12-21 | 2010-12-21 | 一种自动分离动静脉血管的内源光学成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010598859A CN101999885B (zh) | 2010-12-21 | 2010-12-21 | 一种自动分离动静脉血管的内源光学成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101999885A CN101999885A (zh) | 2011-04-06 |
CN101999885B true CN101999885B (zh) | 2011-12-21 |
Family
ID=43807763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010598859A Active CN101999885B (zh) | 2010-12-21 | 2010-12-21 | 一种自动分离动静脉血管的内源光学成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101999885B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104739370A (zh) * | 2013-12-31 | 2015-07-01 | 环达电脑(上海)有限公司 | 具有血管识别定位功能的电子装置及其方法 |
CN105184747B (zh) * | 2015-09-09 | 2019-05-10 | 天津光电高斯通信工程技术股份有限公司 | 低照度图像对比度的提升方法 |
CN106725263B (zh) * | 2016-12-15 | 2018-07-06 | 深圳开立生物医疗科技股份有限公司 | 应用于内窥镜系统的成像方法 |
CN107788950B (zh) * | 2017-09-30 | 2023-08-18 | 浙江大学 | 基于自适应阈值分割的血流成像方法与系统 |
CN108852385B (zh) * | 2018-03-13 | 2022-03-04 | 中国科学院上海应用物理研究所 | 一种x射线造影方法及基于x射线造影的动态阅片方法 |
CN109394268B (zh) * | 2018-12-07 | 2021-05-11 | 刘志红 | 息肉危害程度映射平台 |
CN109949322B (zh) * | 2019-03-27 | 2023-06-23 | 中山大学 | 一种基于磁共振t1增强图像的脑血管图像分割方法 |
CN110477852B (zh) * | 2019-08-02 | 2024-06-18 | 佛山科学技术学院 | 一种虹膜血管成像系统 |
CN111932554B (zh) * | 2020-07-31 | 2024-03-22 | 青岛海信医疗设备股份有限公司 | 一种肺部血管分割方法、设备及存储介质 |
CN111681242B (zh) * | 2020-08-14 | 2020-11-17 | 北京至真互联网技术有限公司 | 视网膜血管动静脉区分方法和装置、设备 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6023968A (en) * | 1995-08-23 | 2000-02-15 | Diasonics Ultrasound, Inc. | Real-time two-dimensional resistance and pulsatility mapping |
JP3548473B2 (ja) * | 1999-11-19 | 2004-07-28 | 日本電信電話株式会社 | 眼底画像の動静脈識別方法、記録媒体及び装置 |
EP1851720B1 (en) * | 2005-02-11 | 2009-07-29 | Koninklijke Philips Electronics N.V. | Method of automatic extraction of the pulmonary artery tree from 3d medical images |
JP5196323B2 (ja) * | 2009-02-23 | 2013-05-15 | 日本光電工業株式会社 | 血中酸素飽和度測定装置 |
-
2010
- 2010-12-21 CN CN201010598859A patent/CN101999885B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN101999885A (zh) | 2011-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101999885B (zh) | 一种自动分离动静脉血管的内源光学成像方法 | |
JP4615865B2 (ja) | 静止したバックグラウンド中の移動物質の特徴づけ | |
Salem et al. | Novel and adaptive contribution of the red channel in pre-processing of colour fundus images | |
CN101697871B (zh) | 一种自动分割动静脉血管的激光成像方法和装置 | |
CN104768446A (zh) | 用光学相干断层扫描血管造影对局部循环进行量化 | |
Blöcher et al. | An online PPGI approach for camera based heart rate monitoring using beat-to-beat detection | |
Mansoor et al. | Real-time low cost infrared vein imaging system | |
Ortíz et al. | Support system for the preventive diagnosis of hypertensive retinopathy | |
Bousefsaf et al. | Peripheral vasomotor activity assessment using a continuous wavelet analysis on webcam photoplethysmographic signals | |
CN103544686A (zh) | 基于相位一致性的眼底图像微动脉瘤检测方法 | |
Tripathi et al. | Automatic detection of exudates in retinal fundus images using differential morphological profile | |
Noronha et al. | Enhancement of retinal fundus image to highlight the features for detection of abnormal eyes | |
Karthika et al. | Retinal image analysis using contourlet transform and multistructure elements morphology by reconstruction | |
CN103006175A (zh) | 一种基于相位一致性的眼底图像视盘定位方法 | |
Tong et al. | Detail-preserving arterial pulse wave measurement based biorthogonal wavelet decomposition from remote rgb observations | |
He et al. | Spatiotemporal monitoring of changes in oxy/deoxy-hemoglobin concentration and blood pulsation on human skin using smartphone-enabled remote multispectral photoplethysmography | |
Zardadi et al. | Unsupervised segmentation of retinal blood vessels using the human visual system line detection model | |
Burton et al. | Feasibility study of remote contactless perfusion imaging with consumer-grade mobile camera | |
Zhao et al. | Separation of cortical arteries and veins in optical neurovascular imaging | |
Deng et al. | Fatigue detection based on infrared video pupillography | |
Wang et al. | Cerebral artery–vein separation using 0.1-Hz oscillation in dual-wavelength optical imaging | |
US8403862B2 (en) | Time-based imaging | |
Bansal et al. | An efficient automatic intensity based method for detection of macula in retinal images | |
Hu et al. | Separation of arteries and veins in the cerebral cortex using physiological oscillations by optical imaging of intrinsic signal | |
CN114092855A (zh) | 基于深度学习的非接触智能人体心率预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |