发明内容
本发明的目的就是要解决上述现有技术存在的不足,提供一种用于胶囊内窥镜的消化道出血图像检测方法,该检测方法可有效识别出图像出血,提高了医生的工作效率。
本发明为实现该目的,提供了一种用于胶囊内窥镜的消化道出血图像检测方法,其特征在于,它包括如下步骤:
步骤1:患者吞服胶囊内窥镜后,胶囊内窥镜在消化道内采集图像,并将图像通过无线通信送入计算机,然后计算机对采集到的消化道图像进行数据预处理;
步骤2:对预处理后的消化道图像分别进行点状出血检测和面状出血检测,确定消化道内是否存在出血,并对点状出血区域和面状出血区域进行标识;
所述点状出血检测的具体步骤为:
步骤201:采用基于Haar特征或LBP特征或Hog特征的Adaboost级联算法训练一个级联分类器,用该级联分类器来对预处理后的消化道图像进行扫描,找出预处理后的消化道图像中的疑似点状出血区域;
步骤202:对得到的这些疑似点状出血区域采用疑似点状出血区域过滤算法进行排除,排除掉非点状出血区域;
步骤203:对疑似点状出血区域中没有被步骤202排除的区域采用疑似点状出血区域分类算法进行分类,分类为点状出血块和非点状出血块,并将点状出血块进行标记;
所述面状出血检测的具体步骤为:
步骤2001:对预处理后的消化道图像进行超像素分割得到多个超像素图像块,超像素分割采用简单线性迭代聚类算法(SLIC,simple linear iterative clustering),得到的超像素图像块的个数与初始设置的种子数有关;
步骤2002:计算每个超像素图像块在RGB颜色空间内的R分量均值meanR、G分量均值meanG和B分量均值meanB,另外,还计算出比值F1、比值F2和比值F3,其中,比值F1=meanR/meanG,比值F2=meanR/meanB,比值F3=meanR/(meanR+meanG+meanB);
步骤2003:计算每个超像素图像块在HSI颜色空间内H分量均值meanH、S分量均值meanS、I分量均值meanI,另外,还计算出比值F4、比值F5和比值F6,其中,比值F4=meanS/meanH,比值F5=meanS/meanI,比值F6=meanS/(meanH+meanS+meanI);
步骤2004:如果上述均值meanR小于阈值TmeanR或比值F1小于阈值TF1或比值F2小于阈值TF2,那么直接判定该超像素图像块不是出血超像素块,否则,该超像素图像块为疑似出血超像素块;
步骤2005:对于疑似出血超像素图像块利用已训练好的支持向量机模型(SVM,Support Vector Machine),输入比值F1、比值F2、比值F3、比值F4、比值F5、比值F6以及均值meanR、均值meanG、均值meanB进行出血超像素块与非出血超像素块的分类,并将出血超像素块进行标记。
本发明具有如下有益效果:
1、本发明可以对胶囊内窥镜采集到的消化道图像进行智能出血的判断,能够辅助医生进行疾病的诊断,减少了医生的工作量。
2、目前还没有针对胃部图像的出血检测方法。现有的出血图像检测都是针对小肠图像的,小肠空间较小,胶囊内窥镜贴着小肠壁拍摄到的病灶部位比较大,较明显。然而胃部空腔体积大,胶囊内窥镜采集到的图像因为拍摄距离的原因,病灶部位有时会比较小,现有的方法难以识别。本发明结合点状出血检测和面状出血检测,可以同时检测小面积的出血和大面积的出血,提高了出血检测的识别率,增强了图像检测方法的适应性。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明:
本发明的一种用于胶囊内窥镜的消化道出血图像检测方法,它包括如下步骤:
步骤1:患者吞服胶囊内窥镜后,胶囊内窥镜在消化道内采集图像,并将图像通过无线通信送入计算机,然后计算机对采集到的消化道图像进行数据预处理;
步骤2:对预处理后的消化道图像分别进行点状出血检测和面状出血检测,确定消化道内是否存在出血,并对点状出血区域和面状出血区域进行标识;
所述点状出血检测的具体步骤为:
步骤201:采用基于Haar特征或LBP特征或Hog特征的Adaboost级联算法训练一个级联分类器,用该级联分类器来对预处理后的消化道图像进行扫描,找出预处理后的消化道图像中的疑似点状出血区域;
步骤202:对得到的这些疑似点状出血区域采用疑似点状出血区域过滤算法进行排除,排除掉非点状出血区域;
步骤203:对疑似点状出血区域中没有被步骤202排除的区域采用疑似点状出血区域分类算法进行分类,分类为点状出血块和非点状出血块,并将点状出血块进行标记;
所述面状出血检测的具体步骤为:
步骤2001:对预处理后的消化道图像进行超像素分割得到多个超像素图像块,超像素分割采用简单线性迭代聚类算法,得到的超像素图像块的个数与初始设置的种子数有关;
步骤2002:计算每个超像素图像块在RGB颜色空间内的R分量均值meanR、G分量均值meanG和B分量均值meanB,另外,还计算出比值F1、比值F2和比值F3,其中,比值F1=meanR/meanG,比值F2=meanR/meanB,比值F3=meanR/(meanR+meanG+meanB);
步骤2003:计算每个超像素图像块在HSI颜色空间内H分量均值meanH、S分量均值meanS、I分量均值meanI,另外,还计算出比值F4、比值F5和比值F6,其中,比值F4=meanS/meanH,比值F5=meanS/meanI,比值F6=meanS/(meanH+meanS+meanI);
步骤2004:如果上述均值meanR小于阈值TmeanR或比值F1小于阈值TF1或比值F2小于阈值TF2,那么直接判定该超像素图像块不是出血超像素块,否则,该超像素图像块为疑似出血超像素块;上述阈值TmeanR的取值范围为40~80,阈值TF1的取值范围为1.0~1.6,阈值TF2的取值范围为1.0~1.6;
步骤2005:对于疑似出血超像素图像块利用已训练好的支持向量机模型(SVM,Support Vector Machine),输入比值F1、比值F2、比值F3、比值F4、比值F5、比值F6以及均值meanR、均值meanG、均值meanB进行出血超像素块与非出血超像素块的分类,并将出血超像素块进行标记。
上述技术方案的步骤202中的疑似点状出血区域过滤算法为:将疑似点状出血区域的图像转换到灰度空间,计算灰度图像的灰度均值,如果灰度均值小于阈值Tgray(阈值Tgray取值范围为40~60),判定该矩形区域为非点状出血区域,并对该非点状出血区域进行排除;
如果灰度均值不小于阈值Tgray,那么对疑似点状出血区域进行进一步的过滤排除,方法为:在RGB颜色空间下,分别获得消化道图像的单个颜色通道分量,即R、G和B分量,计算图像R分量与G分量的比值得到rg,R分量与B分量的比值得到rb,并计算rg的均值rgMean和rb的均值rbMean;
分别根据rg与rgMean和rb与rbMean得到两个掩膜图rgMask和rbMask,计算方法如下:
将掩膜图rgMask和rbMask作与运算得到两个掩膜相交的区域,如果相交区域为空,那么判定该矩形区域为非点状出血区域,实现了对疑似出血矩形区域中非出血区域的进一步排除。
上述技术方案的步骤203:对疑似点状出血区域中没有被步骤202排除的区域采用疑似点状出血区域分类算法进行分类的具体方法为:利用人类视觉系统理论算法(HumanVisual System,HVS)求取人类视觉系统掩膜Mask(HVS)(Mask(HVS)由五个子掩膜共同组成),并利用人类视觉系统掩膜对疑似点状出血区域进行分割,产生前景数据和背景数据,图像中被人类视觉系统掩膜Mask(HVS)覆盖到的部分称为前景数据,未覆盖到的部分称为背景数据,然后对分割后的疑似点状出血区域进行特征提取,提取的特征包括前景区所占的面积比特征f_arearatio、对比度特征f_contrast、偏心率特征f_eccentricity、颜色特征f_rgb和形态学特征f_morpho,然后将上述前景区所占的面积比特征f_arearatio、对比度特征f_contrast、偏心率特征f_eccentricity、颜色特征f_rgb和形态学特征f_morpho送入已训练好的SVM模型中进行分类,将疑似点状出血区域中没有被步骤202排除的区域分类为点状出血区域和非点状出血区域。
上述技术方案中,所述人类视觉系统掩膜Mask(HVS)包括R/G掩膜Mask(R/G)、R/B掩膜Mask(R/B)、边缘掩膜Mask(edge)、轮廓掩膜Mask(outline)和强度掩膜Mask(Intensity),其中,R/G掩膜Mask(R/G)中R和G分别为彩色RGB图像的R分量和G分量,二者相除得到Map(R/G),然后对Map(R/G)进行二值化,二值化计算的阈值设为Map(R/G)的中值,得到R/G掩膜,具体计算公式为:
其中median(Map(R/G))是指求取Map(R/G)的中值;
所述R/B掩膜Mask(R/B)中R和B分别为彩色RGB图像的R分量和B分量,二者相除得到Map(R/B),然后对Map(R/B)进行二值化,二值化计算的阈值设为Map(R/B)的中值,得到R/B掩膜,具体计算公式为:
其中median(Map(R/B))是指求取Map(R/B)的中值;
所述边缘掩膜Mask(edge)用四个滤波核与疑似点状出血区域图像进行卷积操作得到四个边缘图,对每个像素点取四个边缘图中的最大值作为疑似点状出血区域图像的边缘edge_I,然后对疑似点状出血区域图像的边缘edge_I,利用视觉显著度理论求取显著的边缘得到边缘掩膜Mask(edge),具体计算公式为:
其中median是求中值的函数,std是求标准差的函数;
所述轮廓掩膜Mask(outline),通过添加一个半径为图像宽度的一半的圆形掩膜,去掉四周的边界信息得到;
所述强度掩膜Mask(Intensity)的具体计算公式为:
其中I=(R+G+B)/3,R、G、B分别为彩色RGB图像的R分量、G分量和B分量,mean为求均值的函数,std为求标准差的函数;
所述的人类视觉系统掩膜Mask(HVS)由五个子掩膜Mask(R/G)、Mask(R/B)、Mask(edge)、Mask(outline)和Mask(Intensity)相交得到,具体计算公式为:
Mask(HVS)=Mask(R/G)∩Mask(R/B)∩Mask(edge)∩Mask(outline)∩Mask(Intensity).。
上述技术方案中对分割后的疑似点状出血区域进行特征提取,前景区所占的面积比特征f_arearatio、对比度特征f_contrast、偏心率特征f_eccentricity、颜色特征f_rgb和形态学特征f_morpho的具体计算方法为:
所述前景区所占的面积比特征f_arearatio的计算公式为:
其中area函数为求取区域面积的函数;area(foreground)为前景区的面积,即人类视觉系统掩膜Mask(HVS)的面积;width和height分别为疑似点状出血区域图像的宽和高。
对比度特征f_contrast的计算公式为:
其中I_gray为疑似点状出血区域图像的灰度值,max和min分别为求取最大值和最小值的函数,mean表示求均值函数,mean(foreground)表示前景区域的均值,mean(background)表示背景区域的均值,abs表示求绝对值;
偏心率特征f_eccentricity的获取方法为先得到人类视觉系统掩膜Mask(HVS)中面积最大的连通区域,然后求取这个最大连通区域的偏心率;
颜色特征f_rgb由r_mean、rg_mean、rb_mean三个特征组成,其中r_mean,rg_mean,rb_mean的计算公式分别为:
r_mean=mean(R)
其中,mean(R)、mean(G)和mean(B)分别表示R、G、B三个分量的均值;
所述形态学特征f_morpho的获取方式为,先求取人类视觉系统掩膜Mask(HVS)内所有点到区域中心(centerX,centerY)的归一化距离dist,其计算公式为:
Mask(HVS)(x)表示人类视觉系统掩膜的横坐标,Mask(HVS)(y)表示人类视觉系统掩膜的纵坐标,width表示疑似点状出血区域图像的宽,height表示疑似点状出血区域图像的高,sqrt表示表示开方运算;
求归一化距离dist的最小值dist_min,最大值dist_max,均值dist_mean、标准差dist_std以及归一化后直方图的熵dist_h,其计算公式如下:
dist_min=min(dist)
dist_max=max(dist)
dist_mean=mean(dist)
mean为求均值的函数
dist_std=std(dist)
std表示求标准差的函数
其中hist为求直方图的函数,sum为求和的函数;形态学特征f_morpho是将上述特征进行组合:
f_morpho=[dist_min,1-dist_max,1-dist_mean,dist_std,dist_h].
上述技术方案的步骤201中,采用基于Haar特征或LBP特征或Hog特征的Adaboost级联算法训练一个级联分类器,用该级联分类器来对预处理后的消化道图像进行扫描,找出预处理后的消化道图像中的疑似点状出血区域的具体方法为:先选取Haar特征或LBP特征或Hog特征,然后利用带有标签的样本(带有标签的意思是已知样本的类别,如出血样本、非出血样本)训练Adaboost级联分类器,得到一个训练好的分类模型,然后利用这个训练好的分类模型在消化道图像中找疑似点状出血矩形区域。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。