CN108596237A - 一种基于颜色和血管的lci激光内镜下的结肠息肉分类方法 - Google Patents

一种基于颜色和血管的lci激光内镜下的结肠息肉分类方法 Download PDF

Info

Publication number
CN108596237A
CN108596237A CN201810351707.5A CN201810351707A CN108596237A CN 108596237 A CN108596237 A CN 108596237A CN 201810351707 A CN201810351707 A CN 201810351707A CN 108596237 A CN108596237 A CN 108596237A
Authority
CN
China
Prior art keywords
polyp
color
colon
feature vector
value
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.)
Granted
Application number
CN201810351707.5A
Other languages
English (en)
Other versions
CN108596237B (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.)
Pla Third O Seven Hospital
Beijing University of Posts and Telecommunications
Original Assignee
Pla Third O Seven Hospital
Beijing University of Posts and Telecommunications
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 Pla Third O Seven Hospital, Beijing University of Posts and Telecommunications filed Critical Pla Third O Seven Hospital
Priority to CN201810351707.5A priority Critical patent/CN108596237B/zh
Publication of CN108596237A publication Critical patent/CN108596237A/zh
Application granted granted Critical
Publication of CN108596237B publication Critical patent/CN108596237B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30028Colon; Small intestine
    • G06T2207/30032Colon polyp

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供了一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,通过提取LCI激光内镜下的结肠息肉图片的颜色特征向量;将腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型;将非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型;提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型和第二高斯混合模型得到结果值,得出分类结果的技术方案,本发明基于LCI激光内镜,提取颜色和血管密度作为特征向量为结肠息肉分类,该方法所需样本的数量小,训练任务简单,计算量小,在实时监测中无需为此搭设服务器,降低成本。

Description

一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法
技术领域
本发明涉及互联网技术领域,尤其涉及一种基于颜色和血管的LCI 激光内镜下的结肠息肉分类方法。
背景技术
结直肠癌(CRC)是全球癌症中致死的重要成因。尽管早期进行息肉切除手术就可以大大降低其发病率,但只有腺瘤性息肉应该切除,而其余90%的增生息肉在临床上被认为是对健康没有影响的,切除的风险可能超过其好处。故正确地检测和预测结肠息肉类型显得尤为重要。
现有技术中,大多均针对白光(white-light,WL)或蓝激光成像(blue laserimaging,BLI)的激光内镜下的结肠息肉分类方法:
a.提取HOG特征,使用支持向量机(Support vector machine,SVM) 分类;
b.使用神经网络(neural network,NN)构造分类器。
现有的最新实验结果大多将医学上的先验知识与计算机视觉技术剥离开来,使用时下流行的卷积神经网络直接构建分类器。但由于结肠镜图片本身具有较高相似性,这一方法需要大量的训练样本(105数量级及以上)才能达到比较理想的效果。而在临床上收集特定病例图片具有相当的难度,指定类别的息肉高度依赖于病患的实际情况。事实上,现有的公开数据集的数据量不仅大小都在百张左右浮动,其数据图片更是多由同一视频分帧截取得到,训练模型对不同形态的息肉具有低容错性。
发明人在研究的过程中发现:彩色联动成像(linked color imaging, LCI)激光内镜是一种新型激光光源,相较于白光(white-light,WL)、蓝激光成像(blue laserimaging,BLI)等目前应用较广的光源来说,对红色分量更为敏感,拍摄出的下消化道图像具有色彩鲜明、轮廓清晰等优点。
光电容积描记术(photoplethysmography,PPG)是一种简单的光学技术,用于检测外周循环血液的体积变化。由于血液比周围组织更强烈地吸收光纤,所以通过PPG传感器检测为光强度的变化。
高斯混合模型是一种应用广泛的模型,它是用高斯概率密度函数精确地量化事物,以将事物分解为若干的基于高斯概率密度函数形成的模型。
本发明提出了一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,为一种区别增生和腺瘤息肉的分类方法,旨在小样本、小计算量的情况下,基于LCI激光内镜图片的现有特点,综合资深内镜医师的先验临床知识,提取颜色和血管密度特征作为特征向量,完成分类任务。
发明内容
为了达到上述目的,本发明提供了一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,充分利用了LCI激光内镜色彩鲜明、轮廓清晰的特点,提取颜色特征向量和血管密度特征值,相较于现有的支持向量机、神经网络等方法:1.所需样本的数量小,仅为其百分之一甚至千分之一;2.训练任务简单,计算量小,在实时监测中医院不需为此搭设服务器,大大降低成本。
本发明提供了一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,包括:
提取LCI激光内镜下的结肠息肉图片的颜色特征向量;
将结肠息肉图片中的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型A;将结肠息肉图片中的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型B;
提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型A和第二高斯混合模型B得到结果值,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果。
其中,所述颜色特征向量,包括但不限于:LCI激光内镜下的结肠息肉图片中的R、G、B、H、S和/或V分量。
进一步的,所述提取LCI激光内镜下的结肠息肉图片的颜色特征向量的步骤包括:
截取LCI激光内镜下的结肠息肉图片中的息肉区域,息肉区域不包含非息肉的肠道组织;
识别并存储所述息肉区域中所有像素点的R、G、B颜色分量值;将所述息肉区域中RGB颜色空间变换为HSV颜色空间,获取息肉区域中所有像素点的H、S、V颜色分量值;
将获取的R、G、B、H、S、V六个颜色分量值按顺序组成六维的颜色特征向量,使得息肉区域内每一个像素点对应一个六维颜色特征向量。
进一步的,将所述息肉区域中RGB颜色空间变换为HSV颜色空间的步骤,具体包括:
基于Matlab的rgb2hsv()函数,对息肉区域做RGB颜色空间到HSV 颜色空间的颜色空间变换:
v=max
其中,r是颜色的红坐标,g是颜色的绿坐标,b是颜色的蓝坐标, r,g,b的坐标范围值为0到1之间的实数;max等价于r,g和b中的最大者;min等价于r,g和b中的最小者。
进一步的,提取血管密度特征值的步骤包括:
基于PPG技术对LCI激光内镜下具有结肠息肉的视频提取反映微血管中血容量变化的BVP信号视频样本,并在BVP信号视频样本中标记血管密度特征值。
进一步的,提取血管密度特征值的具体步骤包括:
对LCI激光内镜下的一段具有结肠息肉的内镜视频,截取不包含息肉边缘的息肉区域视频样本;
将视频样本去除噪声和抖动,获取BVP信号视频样本;
将BVP信号视频样本分解为R(red),G(green),B(blue)三个颜色通道,分别计算三个通道分量的空间信号平均强度序列;
使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱;
基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值。
进一步的,计算三个通道分量的空间信号平均强度序列的方法包括:
x(n)=[xR(n),xG(n),xB(n)]T,
其中,x(n)为强度序列,n为视频样本的帧数序号,xR(n),xG(n),xB(n)分别对应R,G,B三个颜色通道的空间信号平均强度序列;T为视频样本的时间。
进一步的,使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱,包括:
S(f)=[SR(f),SG(f),SB(f)]T
其中SR(f),SG(f),SB(f)分别对应R,G,B三个颜色通道的能量谱,各个通道的计算遵循如下公式:
S(f)=|X(f)|2=|FFT(x(n))|2
其中,FFT代表快速傅里叶变换,X(f)表示空间信号平均强度序列 x(n)经过变换后的值。
进一步的,基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值,包括:
其中,fHR为能量谱S(f)幅度中极值处的频率值,T为视频样本的时间长度,i表示的视频样本的序号;R、G、B分别表示三个彩色通道。
进一步的,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果,具体包括:
将提取的颜色特征向量分别输入第一高斯混合模型A得到第一结果值;将提取的颜色特征向量分别输入第二高斯混合模型B得到第二结果值;
将第一结果值与第二结果值相减后得到结果值;所述结果值与血管密度特征值的带权值相加后,与存储的分类阈值进行对比,若高于分类阈值则认为该图片中的结肠息肉为腺瘤性结肠息肉;否则,为非腺瘤性结肠息肉。
本发明提供的一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,通过提取LCI激光内镜下的结肠息肉图片的颜色特征向量;将结肠息肉图片中的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型A;将结肠息肉图片中的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型B;提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型A和第二高斯混合模型B得到结果值,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果的技术方案,充分利用了LCI激光内镜色彩鲜明、轮廓清晰的特点,提取颜色和血管密度作为特征向量,相较于现有的支持向量机、神经网络等方法: 1.所需样本的数量小,仅为其百分之一甚至千分之一;2.训练任务简单,计算量小,在实时监测中医院不需为此搭设服务器,大大降低成本。
附图说明
图1为根据本发明的一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法的流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例一
参照图1,图1示出了本发明提供了一种基于颜色和血管的LCI 激光内镜下的结肠息肉分类方法,包括:步骤S110到步骤S130。
在步骤S110中,提取LCI激光内镜下的结肠息肉图片的颜色特征向量。
在步骤S120中,将结肠息肉图片中的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型A;将结肠息肉图片中的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型B。
具体的,预备好的训练样本分别为腺瘤性结肠息肉内镜图片集与非腺瘤性结肠息肉内镜图片集。使用python语句初始化两个独立同分布的高斯混合模型,以训练集的六维颜色特征向量作为输入进行分别训练。其中一个模型以全部的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到模型A。另一个模型以全部的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到模型B。
例如:
from sklearn.mixture import GaussianMixture as GMM
#建立两个独立同分布的高斯混合模型
gmm_1=GMM(n_components=512,tol=1e-4,max_iter=20000)
gmm_2=GMM(n_components=512,tol=1e-4,max_iter=20000)
#分别训练两个模型
高斯混合模型可以用如下公式表示:
其中,模型共有K个高斯分布分量,N(x|μkk)称为模型的第k 个高斯分布分量,μkk为该分量对应的高斯分布参数,πk为该分量在模型中的权重。高斯混合模型的训练过程即在已知x(样本)和p(x)(标签)的条件下,采用EM算法(Expectation-Maximizationalgorithm) 对πkkk进行参数估计。
gmm_1=gmm_1.fit(adenomatous)
gmm_2=gmm_2.fit(inflammatory)
其中,adenomatous代表所有腺瘤性息肉图片的六维特征向量, inflammatory代表所有非腺瘤性息肉的六维特征向量。
在步骤S130中,提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型A 和第二高斯混合模型B得到结果值,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果。
其中,所述颜色特征向量,包括但不限于:LCI激光内镜下的结肠息肉图片中的R、G、B、H、S和/或V分量。
进一步的,所述提取LCI激光内镜下的结肠息肉图片的颜色特征向量的步骤包括:
截取LCI激光内镜下的结肠息肉图片中的息肉区域,息肉区域不包含非息肉的肠道组织;
识别并存储所述息肉区域中所有像素点的R、G、B颜色分量值;将所述息肉区域中RGB颜色空间变换为HSV颜色空间,获取息肉区域中所有像素点的H、S、V颜色分量值;
其中,RGB:将截取后的图片读入到Matlab软件中,图片被识别为一个行数为三的矩阵,通过Matlab截取区域中所有像素点的R、G、 B颜色分量值即分别为矩阵的第1、2、3行值。
HSV空间变换:通过Matlab的rgb2hsv()函数,对截取区域做RGB 颜色空间到HSV颜色空间的颜色空间变换,变换后可直接获取区域中所有像素点的H、S、V颜色分量值。
将获取的R、G、B、H、S、V六个颜色分量值按顺序组成六维的颜色特征向量,使得息肉区域内每一个像素点对应一个六维颜色特征向量。
若截取区域含有四个点,设R,G,B,H,S,V分别为(r1,r2,r3,r4), (g1,g2,g3,g4)等以此类推。则六维特征向量为 A=[r1,g1,b1,h1,s1,v1],B=[r2,g2,b2,h2,s2,v2]等依次类推。
进一步的,将所述息肉区域中RGB颜色空间变换为HSV颜色空间的步骤,具体包括:
基于Matlab的rgb2hsv()函数,对息肉区域做RGB颜色空间到HSV 颜色空间的颜色空间变换:
v=max
其中,r是颜色的红坐标,g是颜色的绿坐标,b是颜色的蓝坐标, r,g,b的坐标范围值为0到1之间的实数;max等价于r,g和b中的最大者;min等价于r,g和b中的最小者。
进一步的,提取血管密度特征值的步骤包括:
基于PPG技术对LCI激光内镜下具有结肠息肉的视频提取反映微血管中血容量变化的BVP(blood volume pulse,BVP)信号视频样本,并在BVP信号视频样本中标记血管密度特征值。
进一步的,提取血管密度特征值的具体步骤包括:
对LCI激光内镜下的一段具有结肠息肉的内镜视频,截取不包含息肉边缘的息肉区域视频样本;
将视频样本去除噪声和抖动,获取BVP信号视频样本;
具体的a.对一段具有结肠息肉的内镜视频,人工截取其息肉区域的一个长度为T秒的视频样本,区域形状应为正方形且不包含息肉边缘。b.为了去除图像光照不均和噪声造成的影响,将样本通过一个0.5Hz-4Hz的滤波器进行时域滤波处理,去除噪声和抖动,以获取BVP 信号。也可以用ICA、小波变换等更复杂方法滤波得到BVP信号。
将BVP信号视频样本分解为R(red),G(green),B(blue)三个颜色通道,分别计算三个通道分量的空间信号平均强度序列;
使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱;
基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值。
进一步的,计算三个通道分量的空间信号平均强度序列的方法包括:
x(n)=[xR(n),xG(n),xB(n)]T,
其中,x(n)为强度序列,n为视频样本的帧数序号,xR(n),xG(n),xB(n)分别对应R,G,B三个颜色通道的空间信号平均强度序列;T为视频样本的时间。
进一步的,使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱,包括:
S(f)=[SR(f),SG(f),SB(f)]T
其中SR(f),SG(f),SB(f)分别对应R,G,B三个颜色通道的能量谱,各个通道的计算遵循如下公式:
S(f)=|X(f)|2=|FFT(x(n))|2
其中,FFT代表快速傅里叶变换,X(f)表示空间信号平均强度序列 x(n)经过变换后的值。
进一步的,基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值,包括:
其中,fHR为能量谱S(f)幅度中极值处的频率值,T为视频样本的时间长度,i表示的视频样本的序号;R、G、B分别表示三个彩色通道,使用BVP信号的功率值Pi表征序号为i的息肉的样本区域血管密度的相对大小,即血管密度特征值。
进一步的,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果,具体包括:
将提取的颜色特征向量分别输入第一高斯混合模型A得到第一结果值;将提取的颜色特征向量分别输入第二高斯混合模型B得到第二结果值;
将第一结果值与第二结果值相减后得到结果值;所述结果值与血管密度特征值的带权值相加后,与存储的分类阈值进行对比,若高于分类阈值则认为该图片中的结肠息肉为腺瘤性结肠息肉;否则,为非腺瘤性结肠息肉。
例如:设X为未知图片中提取出的颜色特征向量,Px为未知图片提取出的血管密度值,threshold为选定的比较阈值。
高斯混合模型可以用如下公式表示:
得到分数的过程即在已知x(颜色特征分量)πkkθ,k(模型对应参数)的条件下计算score=log2(p(x))。
Threshold的选取:将训练样本送入训练好的模型进行打分并与 Px相加,根据分数绘制ROC曲线,选取真阳性率等于真阴性率的分数作为threshold。
score_A=gmm_1.score(X)
socre_B=gmm_2.socre(X)
score_C=score_A–score_B
score_D=score_C+Px
if score_D>=threshold:Picture represents adenomatous polyps
if score_D<threshold:Picture represents inflammatory polyps。
本发明提供的一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,通过提取LCI激光内镜下的结肠息肉图片的颜色特征向量;将结肠息肉图片中的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型A;将结肠息肉图片中的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型B;提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型A和第二高斯混合模型B得到结果值,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果的技术方案,充分利用了LCI激光内镜色彩鲜明、轮廓清晰的特点,提取颜色和血管密度作为特征向量,相较于现有的支持向量机、神经网络等方法: 1.所需样本的数量小,仅为其百分之一甚至千分之一;2.训练任务简单,计算量小,在实时监测中医院不需为此搭设服务器,大大降低成本。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (10)

1.一种基于颜色和血管的LCI激光内镜下的结肠息肉分类方法,其特征在于,包括:
提取LCI激光内镜下的结肠息肉图片的颜色特征向量;
将结肠息肉图片中的腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第一高斯混合模型A;将结肠息肉图片中的非腺瘤性结肠息肉内镜图片集的颜色特征向量作为训练样本进行训练,得到第二高斯混合模型B;
提取未知类型的结肠息肉图片的颜色特征向量和血管密度特征值;将提取的颜色特征向量分别输入第一高斯混合模型A和第二高斯混合模型B得到结果值,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果。
2.如权利要求1所述的方法,其特征在于,所述颜色特征向量,包括但不限于:LCI激光内镜下的结肠息肉图片中的R、G、B、H、S和/或V分量。
3.如权利要求1或2所述的方法,其特征在于,所述提取LCI激光内镜下的结肠息肉图片的颜色特征向量的步骤包括:
截取LCI激光内镜下的结肠息肉图片中的息肉区域,息肉区域不包含非息肉的肠道组织;
识别并存储所述息肉区域中所有像素点的R、G、B颜色分量值;将所述息肉区域中RGB颜色空间变换为HSV颜色空间,获取息肉区域中所有像素点的H、S、V颜色分量值;
将获取的R、G、B、H、S、V六个颜色分量值按顺序组成六维的颜色特征向量,使得息肉区域内每一个像素点对应一个六维颜色特征向量。
4.如权利要求3所述的方法,其特征在于,将所述息肉区域中RGB颜色空间变换为HSV颜色空间的步骤,具体包括:
基于Matlab的rgb2hsv()函数,对息肉区域做RGB颜色空间到HSV颜色空间的颜色空间变换:
v=max
其中,r是颜色的红坐标,g是颜色的绿坐标,b是颜色的蓝坐标,r,g,b的坐标范围值为0到1之间的实数;max等价于r,g和b中的最大者;min等价于r,g和b中的最小者。
5.如权利要求1所述的方法,其特征在于,提取血管密度特征值的步骤包括:
基于PPG技术对LCI激光内镜下具有结肠息肉的视频提取反映微血管中血容量变化的BVP信号视频样本,并在BVP信号视频样本中标记血管密度特征值。
6.如权利要求1或5所述的方法,其特征在于,提取血管密度特征值的具体步骤包括:
对LCI激光内镜下的一段具有结肠息肉的内镜视频,截取不包含息肉边缘的息肉区域视频样本;
将视频样本去除噪声和抖动,获取BVP信号视频样本;
将BVP信号视频样本分解为R(red),G(green),B(blue)三个颜色通道,分别计算三个通道分量的空间信号平均强度序列;
使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱;
基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值。
7.如权利要求6所述的方法,其特征在于,计算三个通道分量的空间信号平均强度序列的方法包括:
x(n)=[xR(n),xG(n),xB(n)]T,
其中,x(n)为强度序列,n为视频样本的帧数序号,xR(n),xG(n),xB(n)分别对应R,G,B三个颜色通道的空间信号平均强度序列;T为视频样本的时间。
8.如权利要求6所述的方法,其特征在于,使用快速傅里叶变换(FFT)的方法得到三个通道分量的空间信号平均强度序列的能量谱,包括:
S(f)=[SR(f),SG(f),SB(f)]T
其中SR(f),SG(f),SB(f)分别对应R,G,B三个颜色通道的能量谱,各个通道的计算遵循如下公式:
S(f)=|X(f)|2=|FFT(x(n))|2
其中,FFT代表快速傅里叶变换,X(f)表示空间信号平均强度序列x(n)经过变换后的值。
9.如权利要求6所述的方法,其特征在于,基于能量谱,计算息肉区域中各彩色通道信号的归一化和功率,得到血管密度特征值,包括:
其中,fHR为能量谱S(f)幅度中极值处的频率值,T为视频样本的时间长度,i表示的视频样本的序号;R、G、B分别表示三个彩色通道。
10.如权利要求1-9之一所述的方法,其特征在于,将所述结果值和所述血管密度特征值相加后与存储的分类阈值做对比得出分类结果,具体包括:
将提取的颜色特征向量分别输入第一高斯混合模型A得到第一结果值;将提取的颜色特征向量分别输入第二高斯混合模型B得到第二结果值;
将第一结果值与第二结果值相减后得到结果值;所述结果值与血管密度特征值的带权值相加后,与存储的分类阈值进行对比,若高于分类阈值则认为该图片中的结肠息肉为腺瘤性结肠息肉;否则,为非腺瘤性结肠息肉。
CN201810351707.5A 2018-04-19 2018-04-19 一种基于颜色和血管的lci激光内镜下的结肠息肉分类装置 Active CN108596237B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810351707.5A CN108596237B (zh) 2018-04-19 2018-04-19 一种基于颜色和血管的lci激光内镜下的结肠息肉分类装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810351707.5A CN108596237B (zh) 2018-04-19 2018-04-19 一种基于颜色和血管的lci激光内镜下的结肠息肉分类装置

Publications (2)

Publication Number Publication Date
CN108596237A true CN108596237A (zh) 2018-09-28
CN108596237B CN108596237B (zh) 2019-11-15

Family

ID=63613570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810351707.5A Active CN108596237B (zh) 2018-04-19 2018-04-19 一种基于颜色和血管的lci激光内镜下的结肠息肉分类装置

Country Status (1)

Country Link
CN (1) CN108596237B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109447973A (zh) * 2018-10-31 2019-03-08 腾讯科技(深圳)有限公司 一种结肠息肉图像的处理方法和装置及系统
CN109636856A (zh) * 2019-01-17 2019-04-16 天津大学 基于hog特征融合算子的物体六维度位姿信息联合测量方法
CN111028219A (zh) * 2019-12-10 2020-04-17 浙江同花顺智能科技有限公司 一种结肠图像识别方法、装置以及相关设备
CN111839445A (zh) * 2019-04-25 2020-10-30 天津御锦人工智能医疗科技有限公司 一种基于图像识别的结肠镜手术中窄带成像检测方法
WO2021073279A1 (zh) * 2019-10-18 2021-04-22 平安科技(深圳)有限公司 数字病理图像染色归一化方法、系统、电子装置及存储介质
CN115496748A (zh) * 2022-11-08 2022-12-20 武汉楚精灵医疗科技有限公司 小肠图像的肠段识别方法、装置及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101048797A (zh) * 2004-08-24 2007-10-03 美国西门子医疗解决公司 在存在标记排泄物或折叠结肠区域的情况下用于结肠壁提取的系统和方法
CN101655912A (zh) * 2009-09-17 2010-02-24 上海交通大学 小波变换的计算机生成图像与自然图像的检测方法
CN102184404A (zh) * 2011-04-29 2011-09-14 汉王科技股份有限公司 掌纹图像中手掌区域的获取方法及装置
CN105893925A (zh) * 2015-12-01 2016-08-24 乐视致新电子科技(天津)有限公司 基于肤色的人手检测方法及装置
CN106097335A (zh) * 2016-06-08 2016-11-09 安翰光电技术(武汉)有限公司 消化道病灶图像识别系统及识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101048797A (zh) * 2004-08-24 2007-10-03 美国西门子医疗解决公司 在存在标记排泄物或折叠结肠区域的情况下用于结肠壁提取的系统和方法
CN101655912A (zh) * 2009-09-17 2010-02-24 上海交通大学 小波变换的计算机生成图像与自然图像的检测方法
CN102184404A (zh) * 2011-04-29 2011-09-14 汉王科技股份有限公司 掌纹图像中手掌区域的获取方法及装置
CN105893925A (zh) * 2015-12-01 2016-08-24 乐视致新电子科技(天津)有限公司 基于肤色的人手检测方法及装置
CN106097335A (zh) * 2016-06-08 2016-11-09 安翰光电技术(武汉)有限公司 消化道病灶图像识别系统及识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TAKUTO SUZUKI等: "Linked-color imaging improves endoscopic visibility ofcolorectal nongranular flat lesions", 《ELSEVIER》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109447973A (zh) * 2018-10-31 2019-03-08 腾讯科技(深圳)有限公司 一种结肠息肉图像的处理方法和装置及系统
US11468563B2 (en) 2018-10-31 2022-10-11 Tencent Technology (Shenzhen) Company Limited Colon polyp image processing method and apparatus, and system
US11748883B2 (en) 2018-10-31 2023-09-05 Tencent Technology (Shenzhen) Company Limited Colon polyp image processing method and apparatus, and system
CN109636856A (zh) * 2019-01-17 2019-04-16 天津大学 基于hog特征融合算子的物体六维度位姿信息联合测量方法
CN109636856B (zh) * 2019-01-17 2022-03-04 天津大学 基于hog特征融合算子的物体六维度位姿信息联合测量方法
CN111839445A (zh) * 2019-04-25 2020-10-30 天津御锦人工智能医疗科技有限公司 一种基于图像识别的结肠镜手术中窄带成像检测方法
WO2021073279A1 (zh) * 2019-10-18 2021-04-22 平安科技(深圳)有限公司 数字病理图像染色归一化方法、系统、电子装置及存储介质
CN111028219A (zh) * 2019-12-10 2020-04-17 浙江同花顺智能科技有限公司 一种结肠图像识别方法、装置以及相关设备
CN115496748A (zh) * 2022-11-08 2022-12-20 武汉楚精灵医疗科技有限公司 小肠图像的肠段识别方法、装置及存储介质
CN115496748B (zh) * 2022-11-08 2023-03-14 武汉楚精灵医疗科技有限公司 小肠图像的肠段识别方法、装置及存储介质

Also Published As

Publication number Publication date
CN108596237B (zh) 2019-11-15

Similar Documents

Publication Publication Date Title
CN108596237B (zh) 一种基于颜色和血管的lci激光内镜下的结肠息肉分类装置
Tosta et al. Segmentation methods of H&E-stained histological images of lymphoma: A review
Sánchez et al. A novel automatic image processing algorithm for detection of hard exudates based on retinal image analysis
US7260248B2 (en) Image processing using measures of similarity
Qureshi et al. Detection of glaucoma based on cup-to-disc ratio using fundus images
Abbas et al. Plasmodium species aware based quantification of malaria parasitemia in light microscopy thin blood smear
CN110619301A (zh) 一种基于双模态信号的情绪自动识别方法
CN110120056A (zh) 基于自适应直方图阈值和轮廓检测的血液白细胞分割方法
CN111798425B (zh) 基于深度学习的胃肠道间质瘤中核分裂象智能检测方法
CN109241898B (zh) 腔镜视像的目标定位方法和系统、存储介质
Kavitha et al. Hierarchical classifier for soft and hard exudates detection of retinal fundus images
Maghsoudi et al. A computer aided method to detect bleeding, tumor, and disease regions in Wireless Capsule Endoscopy
CN110148126A (zh) 基于颜色分量组合和轮廓拟合的血液白细胞分割方法
Alaguselvi et al. Performance analysis of automated lesion detection of diabetic retinopathy using morphological operation
JP2013108926A (ja) 病理診断支援装置、病理診断支援方法、及び病理診断支援プログラム
CN112990015A (zh) 一种病变细胞自动识别方法、装置和电子设备
Ghosh et al. Block based histogram feature extraction method for bleeding detection in wireless capsule endoscopy
Rosebrock et al. Quantitative analysis of TDLUs using adaptive morphological shape techniques
CN116912260B (zh) 基于人工智能的肉鸡养殖健康状态检测方法
Anandgaonkar et al. Brain tumor detection and identification from T1 post contrast MR images using cluster based segmentation
CN104850861B (zh) 基于rx异常探测和纹理分析的真菌性角膜炎图像识别方法
Kayal et al. An approach to detect hard exudates using normalized cut image segmentation technique in digital retinal fundus image
Abdul-Nasir et al. The cascaded moving k-means and fuzzy c-means clustering algorithms for unsupervised segmentation of malaria images
Resita et al. Color RGB and structure GLCM method to feature extraction system in endoscope image for the diagnosis support of otitis media disease
Al-Sharfaa et al. Localization of optic disk and exudates detection in retinal fundus images

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