CN109829931A - 一种基于区域增长pcnn的视网膜血管分割方法 - Google Patents
一种基于区域增长pcnn的视网膜血管分割方法 Download PDFInfo
- Publication number
- CN109829931A CN109829931A CN201910013381.XA CN201910013381A CN109829931A CN 109829931 A CN109829931 A CN 109829931A CN 201910013381 A CN201910013381 A CN 201910013381A CN 109829931 A CN109829931 A CN 109829931A
- Authority
- CN
- China
- Prior art keywords
- image
- vessel
- gray level
- blood
- pixel
- 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
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Image Processing (AREA)
Abstract
一种基于区域增长PCNN的视网膜血管分割方法,包括从目标视网膜血管图像的未标记像素中选择种子点;增大PCNN模型的连接强度,使用连接强度增大后的PCNN模型以种子点为起始点提取目标视网膜血管图像中的血管特征,直到增大后的连接强度大于第一预设阈值;若本次迭代提取的血管特征不同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征对应的像素标记为同一标签,直到所有像素都标记有标签;其中,第一预设条件为血管边缘像素的数量占血管像素总数量的比例小于或等于第二预设阈值;第二预设条件为血管面积占整张图像面积的比例小于或等于第三预设阈值。本发明实现血管区域的自动生长,提高了视网膜血管分割的精度。
Description
技术领域
本发明属于图像处理技术领域,更具体地,涉及一种基于区域增长PCNN的视网膜血管分割方法
背景技术
研究表明,各类眼科疾病以及心脑血管疾病会使得眼底视网膜血管产生形变、出血等不同程度的影响。临床上,医疗人员能够从检眼镜采集的彩色眼底图像中提取视网膜血管,通过对血管形态状况的分析达到诊断这类疾病的目的。
由于受眼底图像采集技术的限制,图像中往往存在大量噪声,再加之视网膜血管自身结构复杂多变,使得视网膜血管的分割变得困难重重。传统方法中依靠人工手动分割视网膜血管,不仅工作量巨大,受主观因素影响严重,而且极为耗时。因此,利用计算机技术,找到一种能够快速、准确分割视网膜的算法,实现对眼底图像血管特征的实时提取,对辅助医疗人员诊断眼科疾病、心脑血管疾病等具有重要作用。
视网膜图像中的血管分割存在以下难点:
(1)血管自身结构特征复杂。视网膜血管的弯曲程度、形状不一,呈树状分布,使得分割起来有一定的困难,传统方法对于血管分割精度不够;
(2)已有的基于PCNN和区域增长的方法[1]分割出的微小血管数量很少,血管分割出现断裂的地方无法继续生长。
发明内容
本发明为解决上述技术问题而提供一种基于区域增长PCNN的视网膜血管分割方法。
为了实现上述发明目的,本发明采取如下技术方案:
一种基于区域增长PCNN的视网膜血管分割方法,包括:
从目标视网膜血管图像的未标记像素中选择种子点;
增大PCNN模型的连接强度,使用连接强度增大后的所述PCNN模型以所述种子点为起始点提取所述目标视网膜血管图像中的血管特征,直到增大后的连接强度大于第一预设阈值;
若本次迭代提取的血管特征不同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征对应的像素标记为同一标签,直到所述目标视网膜血管图像的所有像素都标记有标签;
其中,所述第一预设条件为本次迭代提取的血管特征中血管边缘像素的数量占血管像素总数量的比例小于或等于第二预设阈值;
所述第二预设条件为本次迭代提取的血管特征中血管面积占所述目标视网膜血管图像面积的比例小于或等于第三预设阈值。
若本次迭代提取的血管特征不同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征对应的像素标记为同一标签,直到所述目标视网膜血管图像的所有像素都标记有标签的步骤还包括:
若本次迭代提取的血管特征同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征删除。
从目标视网膜血管图像的未标记像素中选择种子点的步骤之前还包括:
将待处理视网膜血管图像的各通道乘以相应预设权重后相加,获取所述待处理视网膜血管图像的灰度图像;
根据所述灰度图像,获取所述目标视网膜血管图像。
根据所述灰度图像,获取所述目标视网膜血管图像的步骤具体包括:
分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波;
将所述二维高斯滤波方法的滤波结果与二维Gabor滤波方法的滤波结果进行融合,获取融合图像;
将所述灰度图像与所述融合图像进行相减,将相减结果进行取反操作后,乘以预设掩膜图像,获取所述目标视网膜血管图像。
分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波的步骤之前还包括:
若基于四邻域法获知所述灰度图像中各像素为边界,则将所述灰度图像中为边界的各像素进行边缘膨胀;
基于对比度受限制的直方图均衡化方法对边缘膨胀后的所述灰度图像进行对比度增强,获取预处理后的所述灰度图像;
相应地,分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波的步骤具体包括:
分别基于二维高斯滤波方法和二维Gabor滤波方法对预处理后的所述灰度图像进行滤波。
若基于四邻域法获知所述灰度图像中各像素为边界,则将所述灰度图像中为边界的各像素进行边缘膨胀的步骤具体包括:
若所述灰度图像中各像素的四邻域中至少有一个值为0且不全为0,则获知所述灰度图像中各像素为边界;
将所述灰度图像中为边界的各像素进行边缘膨胀。
基于对比度受限制的直方图均衡化方法对边缘膨胀后的所述灰度图像进行对比度增强,获取预处理后的所述灰度图像的步骤具体包括:
将边缘膨胀后的所述灰度图像划分为多个块;
计算出各所述块的灰度直方图,并对各所述块的对灰度直方图进行裁剪;
对各所述块裁剪后的灰度直方图进行均衡化;
基于线性插值法将均衡化后灰度直方图对应的块进行连接,获取预处理后的所述灰度图像。
所述二维高斯滤波方法的第i个高斯卷积核K′i(x,y)的公式为:
K′i(x,y)=Ki(x,y)-mi;
Ki(x,y)=-exp(-u2/(2σ2));
其中,Ki(x,y)为第i个高斯核矩阵中的系数,mi为第i个高斯核的均值,为第i个高斯核的离散点,u、v分别为第i个高斯核在x轴和y轴上的值,σ为高斯核距离x轴坐标中心的偏移量,N={[u,v]||u|≤3σ,|v|≤L/2},L为在y轴方向上被截断的血管长度,A为N中的点数。
将所述二维高斯滤波方法的滤波结果与二维Gabor滤波方法的滤波结果进行融合,获取融合图像的步骤具体包括:
将所述二维高斯滤波方法的滤波结果和二维Gabor滤波方法的滤波结果分别乘以相应预设权重后相加,获取融合图像。
将相减结果进行取反操作后,乘以预设掩膜图像,获取所述目标视网膜血管图像的步骤具体包括:
对预设掩模图像相继进行开运算、闭运算和腐蚀操作;
将处理后的所述预设掩模图像与取反结果进行相乘,获取所述目标视网膜血管图像。
本发明一种视网膜血管分割方法,有益效果如下:
(1)从DRIVE数据库中随机抽出10幅图,通过准确度(Acc),灵敏度(Sen),特异性(Spe)的计算,反映出预处理中使用高斯滤波与Gabor滤波相结合的办法,比仅使用Gabor滤波或高斯滤波处理结果好。
(2)本发明在对视网膜血管预处理的基础上首先选择出一定的血管区域作为初始种子区域;接着,将带有快速连接机制的PCNN与种子区域增长思想相结合,通过设置PCNN中的连接强度系数最大值和停止条件,如果生长无效则会删除掉那部分伪血管,解决了以前方法中断点不能继续生长的问题,实现眼底图像中血管区域的自动生长,从而有效地提取出视网膜图像中的血管。
附图说明
图1是视网膜血管分割流程图。
图2是彩色眼底视网膜图像原图。
图3是图2按比例提取三个通道的图。
图4是图3边缘膨胀后的图。
图5是图4经CLAHE处理后的图。
图6是图5二维高斯匹配后的结果图。
图7是图6二维Gabor匹配滤波处理后结果图。
图8是图7二维Gabor二维高斯融合后的图。
图9是图8相减的图。
图10是图9掩膜图。
图11是图10最终预处理图像。
图12是图11最终分割的图像。
图13是PCNN神经元模型图。
图14是文献[1]提取出视网膜图像中的血管图。
图15是本发明提取出视网膜图像中的血管图。
具体实施方式
一种基于区域增长PCNN的视网膜血管分割方法,括以下步骤,流程图如图1所示:
步骤1:从标准图像库DRIVE中选择一幅彩色眼底视网膜图像作为待处理视网膜血管图像,如图2所示;
步骤2:按比例提取红绿蓝通道图像Y=0.299R+0.587G+0.114B,如图3所示;
步骤3:用四邻域法判断是否为边界,若判断为边界,则将边缘膨胀,如图4所示;
步骤4:采用对比度受限制的直方图均衡化(CLAHE)的方法增强图像的对比度,如图5所示;
步骤5:采用二维高斯滤波的方法增强视网膜血管与背景之间的对比度,如图6所示
步骤6:采用二维Gabor滤波的方法进一步增强视网膜血管与背景之间的对比度,如图7所示;
步骤7:将二维Gabor与二维高斯匹配滤波结果相融合,如图8所示
步骤8:将按红绿蓝成相应比例提取的图像与经过匹配滤波处理融合的图像相减,如图9所示
步骤9:进行取反操作,乘以掩膜图像,掩膜图像如图10所示,即可得到预处理结果,即目标视网膜血管图像,如图11所示;
步骤10:将带有快速连接机制的PCNN与种子区域增长思想相结合,通过阈值操作首先选择出可靠的血管区域作为种子,然后通过PCNN中的连接强度系数和停止条件,实现血管区域的自动生长,从而有效提取视网膜血管,如图12所示。
进一步的,所述步骤2具体为:考虑到虽然绿色通道内含有大量有用信息,但是红色通道和蓝色通道仍然含有有用信息,因此从三个不同通道内按照红色提取29.9%,绿色提取58.7%,蓝色提取11.4%
进一步的,所述步骤3具体为:
步骤3.1:用四邻域的方法判断,如若四个值里有一个值为0但不全为0即为图像的边界
步骤3.2:若为边界则使用膨胀腐蚀的方法使边界处理更加平滑
进一步的,所述步骤4具体为:
步骤4.1:将原始图像分为M×N个大小相等的连续但不重叠的块;
步骤4.2:计算出每个块的直方图;
步骤4.3:对每个子区域进行灰度直方图“剪切”;
步骤4.4:对每个块在进行直方图均衡化;
步骤4.5:把各个块用线性插值法连接成一幅图像,即为变换后的图像。
进一步的,所述步骤5具体为:
步骤5.1:将血管每15度旋转一次,得到从0度到180度的12个匹配滤波器;
所述的二维匹配滤波器的高斯核是一个二维矩阵,假设q=[x,y]是这个高斯核矩阵中的离散点,φi表示n个匹配滤波器中的第i个高斯核的旋转角度。在计算高斯核旋转矩阵中的系数时,假设旋转中心在(0,0)点处,旋转矩阵为:
其中,
旋转后的坐标系中的点为:
步骤5.2:计算匹配滤波器的高斯卷积核;
由于高斯曲线两端沿x轴正负方向无限扩展,需要对其进行截断:u=±3σ,|v|≤L/2。因此,第i个高斯核矩阵中的系数为:
其中,N={[u,v]||u|≤3σ,|v|≤L/2}。
式中,K(x,y)为核函数,u、v分别为在x轴和y轴上的值,σ是高斯函数距离x轴坐标中心的偏移量。σ越大,函数值越向x轴延伸,反之,σ越小,函数值就越向x轴收缩。在y轴方向上被截断的血管长度由L表示。
匹配滤波时,为了使图像原来的背景灰度特性不发生改变,要求卷积核系数平均值为零,因此,匹配滤波器的高斯核的卷积核为:
Ki'(x,y)=Ki(x,y)-mi
其中:
为高斯核的均值,A为N中的点数。
步骤5.3:匹配滤波器设计好之后,利用该滤波器对图像进行卷积运算,然后选择其中最大的卷积值作为增强图像的像素值。假设f(x,y)为原始图像,f'(x,y)为进行匹配滤波增强之后的图像,用n个高斯匹配滤波器对图像中的血管增强的计算公式为:
其中u、v分别为x轴和y轴上的值。
进一步的步骤6所示:
步骤6.1所示:将血管每5度旋转一次,得到从-90度到90度的36个匹配滤波器
所述的二维gabor滤波器的高斯核也是一个二维矩阵,在计算高斯核旋转矩阵中的系数时,假设旋转中心在(0,0)点处,旋转矩阵为:
其中,
旋转后的坐标系中的点为:
步骤6.2所示:计算匹配滤波器的高斯卷积核;
步骤5.3匹配滤波器设计好之后,利用该滤波器对图像进行卷积运算,选其中最大的卷积核值作为增强图像的像素值。假设f(x,y)为原始图像,f'(x,y)为进行匹配滤波增强之后的图像,用n个高斯匹配滤波器对图像中的血管增强的计算公式为:
仅使用高斯滤波在微弱血管区域响应上有一定的优势,但是响应产生的背景噪声也会较多,仅使用Gabor滤波器组,微弱血管区域响应不足,且容易产生虚假响应,因此将两者融合,各自发挥自身优势。如图14所示。
进一步的步骤7所示:
将上述分别用高斯滤波之后的结果与经过Gabor滤波之后的结果融合即图像相加,如图15所示,通过对比选用高斯滤波之后的结果与Gabor滤波之后的结果6:4进行融合。
进一步的,所述步骤8具体为:
步骤8.1:将经过三个通道后的图像中的每个像素点减去经过CLAHE与匹配滤波融合处理后的图像对应位置的像素点;
步骤8.2:进行取反运算,即1-F(x,y);
进一步的,所述步骤9具体为:
步骤9.1:使用3×3大小的掩膜进行开运算。
步骤9.2:使用3×3大小的掩膜进行闭运算。
步骤9.3:使用3×3大小的掩膜腐蚀操作。
步骤9.4:将处理后的掩膜乘以处理后的视网膜血管图像得到最终预处理图像11
步骤10:将带有快速连接机制的PCNN与区域生长思想相结合
进一步的,步骤10中所述的PCNN简化模型是20世纪90年代形成和发展的与传统人工神经网络有着根本不同的新型神经网络,它是根据猫的大脑皮层中神经元同步发放脉冲这一现象提出的,在图像分割、边缘检测、降噪、特征提取等领域中有着广泛的应用。PCNN应用于眼底图像处理时,其神经元数目与图像像素数目一致,各神经元与各像素一一对应。PCNN模型可以分为三个部分:输入域、调制域和脉冲产生域。输入域又可以分为:反馈输入域和连接输入域。如图12所示。
进一步的,PCNN用迭代公式表示如下:
Fij[n]=Sij (1)
Uij[n]=Fij[n]{1+βLij[n]} (3)
式中,Sij为外部刺激,即点(i,j)对应像素的灰度值,Fij为神经元的反馈输入项,Lij为神经元的连接输入项,表示一个神经元的周围神经元的输出进行加权求和,Uij为内部活动项,Yij为神经元的脉冲输出,θij为动态阈值。VL为连接输入域的放大系数,β为神经元之间的连接强度系数,Vθ和αθ分别为动态阈值的放大系数和衰减常数因子,Wijkl为连接加权系数矩阵。
通过将内部活动项与动态阈值进行比较,来决定神经元是否输出脉冲,当Uij>θij时,神经元被激发输出脉冲,同时会使得阈值增大,以阻止该神经元再次被激发,不断地迭代,直到某一时刻Uij再次大于θij时,神经元被再次激发,输出脉冲。
基于区域增长的PCNN算法(RG-PCNN)对简化PCNN模型进行了一些改进,在耦合部分如公式(9)连接系数β被集合βn取代:
Uij[n]=Fij[n](1+βn(t)Lij[n]) (6)
在每次迭代中,β值是变化的,先给β设定一个初值,使得第一次发出脉冲的神经元能够就近捕获下一个神经元,然后反复的用设定的较小的δβ来增加它的值。将β变化前后激发的神经元对应的像素区域的边缘比例以及面积比例分别进行比较,如果小于某个阈值,就继续增加β,否则结束迭代过程。
本发明使用单程条件来解决算法的终止,即每个神经元允许发放脉冲,状态从0变1只有一次,算法过程中,将一个区域对应的神经元同时激发。发放脉冲的迭代次数作为区域标记,迭代次数存储在一个矩阵P中,所有元素的初始值设置为0。
为了实现时间独立性,本发明采用发出脉冲的神经元确定的阈值Wn,n表示迭代次数。其他神经元的阈值设置一个很大的阈值Ω,以防止同一位置重复激发。
将Wn设置为没有激发状态神经元对应的像素的最大值,即将一个强度最大的没有激发过的神经元的初始状态设置为激发状态。如果没激发的神经元具有相同的最大强度值,则只选择一个激发。
由于在PCNN模型中连接域的脉冲传输具有延时性,这样会导致图像中出现断点或者产生不连续的区域。为了减少这些影响,采用快速连接模型,它可以使得已经激发的神经元所引起的自动波在一次迭代过程中充分传播,进而使所有被捕获的神经元可以在同一次迭代过程中被激发从而产生脉冲。
为了得到更好的血管分割结果,本文选取以下几种情况作为局部迭代及全局迭代终止条件:
(1)所有的神经元都被激发;
(2)β大于βmax;
(3)生成血管边缘像素点数量占血管像素点总数的比例e小于设定的阈值;
(4)血管占整幅图的面积比m大于设定阈值;
(5)快速连接终止。
其中(1)为全局迭代条件,(2)(3)(4)(5)为局部迭代条件。
条件(1)即整幅图片中每个像素都被划分到不同区域,不存在没有划分的像素,该条件可以通过判断P中有没有为0的元素来确定。因为P中的初始值都为0,每经过一次迭代,迭代次数赋给激发神经元对应的像素;条件(2)即每次迭代的连接系数强度β大于给定值βmax则终止。相比于原始RG-PCNN,本文根据实际应用场景增加了(3)(4)两个终止条件。其中(3)是为了判断当前种子生长合适的生长区域,如果种子生长的血管边沿像素点总个数与血管总像素点个数的比值小于设定的阈值e则删除掉这部分伪血管,终止本次迭代,转而从剩下未激发的像素点中选择最大的种子点继续上述步骤。(4)是限制血管生长总面积占整幅图的比例小于设定的阈值m,防止血管过生长。条件(5)即没有神经元再被捕获。比较迭代前后激发的神经元个数,如果相等,则停止,如果迭代后激发的神经元数量仍在增加,则迭代继续。
进一步的,计算等式(1)-(5)式进行第一次迭代,重复计算:
Uij[n]=Fij[n]{1+β(n)Lij[n]} (10)
直到Y不再发生变化为止,然后按照公式(12)更新阈值,开始新一轮的迭代:
所述步骤10具体为:
步骤10.1:设置连接强度系数初始值β、连接系数的增加项δβ、连接强度系数最大值的初始值βmax,生成血管边缘像素点数量占血管像素点总数的比例e,血管占整幅图的面积比m;停止阈值T=255;
步骤10.2:将预处理结果进行归一化Sy,作为PCNN的外部输入Fij[n]=Sy;
步骤10.3:选取的种子点,将Sij中最大的像素点作为种子点将它的值赋为1,Y(ijseed)=1;
步骤10.4:设定PCNN阈值Wn取没有激发状态神经元对应像素的最大值,即Tij[n]=Wn,其他神经元阈值被设置为一个很大的常数Ω(默认为500);
步骤10.5:引入快速连接机制进行迭代分割,即定义一个与PCNN_Y同阶的矩阵Y0,使得Y0=1,当β≤βmax时,如果Y0与PCNN_Y中存在不相等的值,则Y0=PCNN_Y,然后引入快速连接机制进行迭代,重新计算PCNN_Y;
步骤10.6:反复循环,不断地比较迭代后的PCNN_Y和迭代前的值Y0,直到两者相等为止;
步骤10.7:用较小的常数δβ来增加它的值,并将Y0的值再赋为1,开始新的循环;
步骤10.8:重复上述过程,直到β>βmax为止,血管完成第一轮生长;
步骤10.9:最后通过计算生成血管边缘像素点数量占血管像素点总数的比例e小于等于0.2时并且血管占整幅图的面积比m小于等于0.11时,生长无效,删除掉一部分伪血管,开始新一轮地生长,并设置了T=255作为停止条件。
整个算法描述如下:
设置β0、δβ、βmax、T、e、m;
将预处理结果归一化为Sy,作为PCNN的外部输入:Fij[n]=Sy;
选取的种子点,将Sij像素值中最大的点对应像素值赋为1,即Y(ijseed)=1;
找到没有激发神经元对应像素值中最大值置为阈值Wn,即Tij=Wn,其他神经元的阈值被设置为一个很大的常数Ω(默认为500)
whileβ≤βmax
while PCNN_Y!=Y0
Y0=PCNN_Y
Uij[n]=Fij[n]{1+β(n)Lij[n]}
end。
β=β+δβ
Y0=1
end
统计生成血管边缘像素点比例e,血管占整幅图的面积比m:
if sum(e>0.2&&m>0.11)
n=n+1
统计P==0索引值
end
本发明的有益效果是:
(1):从DRIVE数据库中随机抽出10幅图,通过准确度(Acc),灵敏度(Sen),特异性(Spe)的计算,反映出预处理中使用高斯滤波与Gabor滤波相结合的办法,比仅使用Gabor滤波或高斯滤波处理结果好。
表1 Gabor滤波的结果
表2高斯滤波的结果
表3 Gabor滤波与高斯滤波融合的结果
(2):提供了一种基于区域增长PCNN的视网膜血管分割方法。本发明在对视网膜血管预处理的基础上首先选择出一定的血管区域作为初始种子区域;接着,将带有快速连接机制的PCNN与种子区域增长思想相结合,通过设置PCNN中的连接强度系数最大值和停止条件,如果生长无效则会删除掉那部分伪血管,解决了以前方法中断点不能继续生长的问题,实现眼底图像中血管区域的自动生长,从而有效地提取出视网膜图像中的血管。
Claims (10)
1.一种基于区域增长PCNN的视网膜血管分割方法,特征在于,包括:
从目标视网膜血管图像的未标记像素中选择种子点;
增大PCNN模型的连接强度,使用连接强度增大后的所述PCNN模型以所述种子点为起始点提取所述目标视网膜血管图像中的血管特征,直到增大后的连接强度大于第一预设阈值;
若本次迭代提取的血管特征不同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征对应的像素标记为同一标签,直到所述目标视网膜血管图像的所有像素都标记有标签;
其中,所述第一预设条件为本次迭代提取的血管特征中血管边缘像素的数量占血管像素总数量的比例小于或等于第二预设阈值;
所述第二预设条件为本次迭代提取的血管特征中血管面积占所述目标视网膜血管图像面积的比例小于或等于第三预设阈值。
2.根据权利要求1所述的方法,其特征在于,若本次迭代提取的血管特征不同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征对应的像素标记为同一标签,直到所述目标视网膜血管图像的所有像素都标记有标签的步骤还包括:
若本次迭代提取的血管特征同时满足第一预设条件和第二预设条件,则将本次迭代提取的血管特征删除。
3.根据权利要求1所述的方法,其特征在于,从目标视网膜血管图像的未标记像素中选择种子点的步骤之前还包括:
将待处理视网膜血管图像的各通道乘以相应预设权重后相加,获取所述待处理视网膜血管图像的灰度图像;
根据所述灰度图像,获取所述目标视网膜血管图像。
4.根据权利要求3所述的方法,其特征在于,根据所述灰度图像,获取所述目标视网膜血管图像的步骤具体包括:
分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波;
将所述二维高斯滤波方法的滤波结果与二维Gabor滤波方法的滤波结果进行融合,获取融合图像;
将所述灰度图像与所述融合图像进行相减,将相减结果进行取反操作后,乘以预设掩膜图像,获取所述目标视网膜血管图像。
5.根据权利要求4所述的方法,其特征在于,分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波的步骤之前还包括:
若基于四邻域法获知所述灰度图像中各像素为边界,则将所述灰度图像中为边界的各像素进行边缘膨胀;
基于对比度受限制的直方图均衡化方法对边缘膨胀后的所述灰度图像进行对比度增强,获取预处理后的所述灰度图像;
相应地,分别基于二维高斯滤波方法和二维Gabor滤波方法对所述灰度图像进行滤波的步骤具体包括:
分别基于二维高斯滤波方法和二维Gabor滤波方法对预处理后的所述灰度图像进行滤波。
6.根据权利要求5所述的方法,其特征在于,若基于四邻域法获知所述灰度图像中各像素为边界,则将所述灰度图像中为边界的各像素进行边缘膨胀的步骤具体包括:
若所述灰度图像中各像素的四邻域中至少有一个值为0且不全为0,则获知所述灰度图像中各像素为边界;
将所述灰度图像中为边界的各像素进行边缘膨胀。
7.根据权利要求5所述的方法,其特征在于,基于对比度受限制的直方图均衡化方法对边缘膨胀后的所述灰度图像进行对比度增强,获取预处理后的所述灰度图像的步骤具体包括:
将边缘膨胀后的所述灰度图像划分为多个块;
计算出各所述块的灰度直方图,并对各所述块的对灰度直方图进行裁剪;
对各所述块裁剪后的灰度直方图进行均衡化;
基于线性插值法将均衡化后灰度直方图对应的块进行连接,获取预处理后的所述灰度图像。
8.根据权利要求4所述的方法,其特征在于,所述二维高斯滤波方法的第i个高斯卷积核K'i(x,y)的公式为:
K'i(x,y)=Ki(x,y)-mi;
Ki(x,y)=-exp(-u2/(2σ2));
其中,Ki(x,y)为第i个高斯核矩阵中的系数,mi为第i个高斯核的均值,为第i个高斯核的离散点,u、v分别为第i个高斯核在x轴和y轴上的值,σ为高斯核距离x轴坐标中心的偏移量,N={[u,v]||u|≤3σ,|v|≤L/2},L为在y轴方向上被截断的血管长度,A为N中的点数。
9.根据权利要求4所述的方法,其特征在于,将所述二维高斯滤波方法的滤波结果与二维Gabor滤波方法的滤波结果进行融合,获取融合图像的步骤具体包括:
将所述二维高斯滤波方法的滤波结果和二维Gabor滤波方法的滤波结果分别乘以相应预设权重后相加,获取融合图像。
10.根据权利要求4所述的方法,其特征在于,将相减结果进行取反操作后,乘以预设掩膜图像,获取所述目标视网膜血管图像的步骤具体包括:
对预设掩模图像相继进行开运算、闭运算和腐蚀操作;
将处理后的所述预设掩模图像与取反结果进行相乘,获取所述目标视网膜血管图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910013381.XA CN109829931B (zh) | 2019-01-07 | 2019-01-07 | 一种基于区域增长pcnn的视网膜血管分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910013381.XA CN109829931B (zh) | 2019-01-07 | 2019-01-07 | 一种基于区域增长pcnn的视网膜血管分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109829931A true CN109829931A (zh) | 2019-05-31 |
CN109829931B CN109829931B (zh) | 2023-07-11 |
Family
ID=66860796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910013381.XA Active CN109829931B (zh) | 2019-01-07 | 2019-01-07 | 一种基于区域增长pcnn的视网膜血管分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109829931B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110288588A (zh) * | 2019-07-01 | 2019-09-27 | 齐鲁工业大学 | 基于灰度方差和标准差的视网膜图像血管分割方法及*** |
CN114155193A (zh) * | 2021-10-27 | 2022-03-08 | 北京医准智能科技有限公司 | 一种基于特征强化的血管分割方法及装置 |
CN117974692A (zh) * | 2024-03-29 | 2024-05-03 | 贵州毅丹恒瑞医药科技有限公司 | 一种基于区域生长的眼科医学影像处理方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1483727A2 (en) * | 2002-03-04 | 2004-12-08 | Isis Innovation Limited | Unsupervised data segmentation |
WO2005048190A1 (en) * | 2003-11-13 | 2005-05-26 | Centre Hospitalier De L'universite De Montreal (Chum) | Automatic multi-dimensional intravascular ultrasound image segmentation method |
US20120195480A1 (en) * | 2011-02-01 | 2012-08-02 | Universidade Da Coruna | Method, apparatus, and system for retinal image analysis |
US20120249959A1 (en) * | 2011-03-31 | 2012-10-04 | The Hong Kong Polytechnic University | Apparatus and method for non-invasive diabetic retinopathy detection and monitoring |
WO2015003225A1 (en) * | 2013-07-10 | 2015-01-15 | Commonwealth Scientific And Industrial Research Organisation | Quantifying a blood vessel reflection parameter of the retina |
WO2015063054A1 (en) * | 2013-10-30 | 2015-05-07 | Agfa Healthcare | Vessel segmentation method |
JP2015202236A (ja) * | 2014-04-15 | 2015-11-16 | キヤノン株式会社 | 眼底画像解析装置及び解析方法 |
CN106815853A (zh) * | 2016-12-14 | 2017-06-09 | 海纳医信(北京)软件科技有限责任公司 | 对眼底图像中视网膜血管的分割方法和装置 |
CN107016676A (zh) * | 2017-03-13 | 2017-08-04 | 三峡大学 | 一种基于pcnn的视网膜血管图像分割方法和*** |
CN108090899A (zh) * | 2017-12-27 | 2018-05-29 | 重庆大学 | 一种血管提取及去噪方法 |
CN108109159A (zh) * | 2017-12-21 | 2018-06-01 | 东北大学 | 一种基于hessian矩阵和区域增长相结合的视网膜血管分割*** |
CN109003279A (zh) * | 2018-07-06 | 2018-12-14 | 东北大学 | 一种基于K-Means聚类标注和朴素贝叶斯模型的眼底视网膜血管分割方法及*** |
-
2019
- 2019-01-07 CN CN201910013381.XA patent/CN109829931B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1483727A2 (en) * | 2002-03-04 | 2004-12-08 | Isis Innovation Limited | Unsupervised data segmentation |
WO2005048190A1 (en) * | 2003-11-13 | 2005-05-26 | Centre Hospitalier De L'universite De Montreal (Chum) | Automatic multi-dimensional intravascular ultrasound image segmentation method |
US20120195480A1 (en) * | 2011-02-01 | 2012-08-02 | Universidade Da Coruna | Method, apparatus, and system for retinal image analysis |
US20120249959A1 (en) * | 2011-03-31 | 2012-10-04 | The Hong Kong Polytechnic University | Apparatus and method for non-invasive diabetic retinopathy detection and monitoring |
WO2015003225A1 (en) * | 2013-07-10 | 2015-01-15 | Commonwealth Scientific And Industrial Research Organisation | Quantifying a blood vessel reflection parameter of the retina |
US20160260212A1 (en) * | 2013-10-30 | 2016-09-08 | Agfa Healthcare | Vessel segmentation method |
WO2015063054A1 (en) * | 2013-10-30 | 2015-05-07 | Agfa Healthcare | Vessel segmentation method |
JP2015202236A (ja) * | 2014-04-15 | 2015-11-16 | キヤノン株式会社 | 眼底画像解析装置及び解析方法 |
CN106815853A (zh) * | 2016-12-14 | 2017-06-09 | 海纳医信(北京)软件科技有限责任公司 | 对眼底图像中视网膜血管的分割方法和装置 |
CN107016676A (zh) * | 2017-03-13 | 2017-08-04 | 三峡大学 | 一种基于pcnn的视网膜血管图像分割方法和*** |
CN108109159A (zh) * | 2017-12-21 | 2018-06-01 | 东北大学 | 一种基于hessian矩阵和区域增长相结合的视网膜血管分割*** |
CN108090899A (zh) * | 2017-12-27 | 2018-05-29 | 重庆大学 | 一种血管提取及去噪方法 |
CN109003279A (zh) * | 2018-07-06 | 2018-12-14 | 东北大学 | 一种基于K-Means聚类标注和朴素贝叶斯模型的眼底视网膜血管分割方法及*** |
Non-Patent Citations (6)
Title |
---|
STEWART R D ETAL.: "Region growing with pulse-coupled neural networks: an alternative to seeded region growing", 《IEEE TRANSACTIONS ON NEURAL NETWORKS》 * |
STEWART R D ETAL.: "Region growing with pulse-coupled neural networks: an alternative to seeded region growing", 《IEEE TRANSACTIONS ON NEURAL NETWORKS》, vol. 13, no. 6, 31 December 2002 (2002-12-31), pages 1557 - 1562 * |
姜平 等: "基于局部自适应点特异度阈值的眼底图像血管分割方法研究", 《南京大学学报(自然科学)》 * |
姜平 等: "基于局部自适应点特异度阈值的眼底图像血管分割方法研究", 《南京大学学报(自然科学)》, vol. 51, no. 01, 30 January 2015 (2015-01-30), pages 171 - 179 * |
解立志;周明全;田沄;武仲科;王醒策;: "基于区域增长与局部自适应C-V模型的脑血管分割", 软件学报, no. 08, pages 237 - 246 * |
黄琳;沈建新;罗煦;: "视网膜图像中的血管自适应提取", 中国制造业信息化, no. 01, pages 68 - 71 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110288588A (zh) * | 2019-07-01 | 2019-09-27 | 齐鲁工业大学 | 基于灰度方差和标准差的视网膜图像血管分割方法及*** |
CN114155193A (zh) * | 2021-10-27 | 2022-03-08 | 北京医准智能科技有限公司 | 一种基于特征强化的血管分割方法及装置 |
CN117974692A (zh) * | 2024-03-29 | 2024-05-03 | 贵州毅丹恒瑞医药科技有限公司 | 一种基于区域生长的眼科医学影像处理方法 |
CN117974692B (zh) * | 2024-03-29 | 2024-06-07 | 贵州毅丹恒瑞医药科技有限公司 | 一种基于区域生长的眼科医学影像处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109829931B (zh) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106920227B (zh) | 基于深度学习与传统方法相结合的视网膜血管分割方法 | |
CN111815574B (zh) | 一种基于粗糙集神经网络的眼底视网膜血管图像分割方法 | |
WO2018082084A1 (zh) | 融合全卷积神经网络和条件随机场的脑肿瘤自动分割方法 | |
CN111292338B (zh) | 一种从眼底oct图像中分割脉络膜新生血管的方法及*** | |
CN107016676B (zh) | 一种基于pcnn的视网膜血管图像分割方法和*** | |
CN110276356A (zh) | 基于r-cnn的眼底图像微动脉瘤识别方法 | |
CN109087302A (zh) | 一种眼底图像血管分割方法及设备 | |
CN112016626B (zh) | 基于不确定度的糖尿病视网膜病变分类*** | |
CN106683080B (zh) | 一种视网膜眼底图像预处理方法 | |
CN108095683A (zh) | 基于深度学习的处理眼底图像的方法和装置 | |
CN109544540A (zh) | 一种基于图像分析技术的糖尿病视网膜图像质量检测方法 | |
CN109829931A (zh) | 一种基于区域增长pcnn的视网膜血管分割方法 | |
CN104732499A (zh) | 基于多尺度多方向的视网膜图像增强算法 | |
Subramani et al. | Fuzzy contextual inference system for medical image enhancement | |
Kumar et al. | Image processing in diabetic related causes | |
CN109087310A (zh) | 睑板腺纹理区域的分割方法、***、存储介质及智能终端 | |
Argade et al. | Automatic detection of diabetic retinopathy using image processing and data mining techniques | |
CN115775210A (zh) | 一种广域眼底相机图像融合方法、***及存储介质 | |
Acharya et al. | Swarm intelligence based adaptive gamma corrected (SIAGC) retinal image enhancement technique for early detection of diabetic retinopathy | |
Liu et al. | Retinal vessel segmentation using densely connected convolution neural network with colorful fundus images | |
Yue et al. | Deep pyramid network for low-light endoscopic image enhancement | |
Shaikha et al. | Optic Disc Detection and Segmentation in Retinal Fundus Image | |
Alimanov et al. | A hybrid approach for retinal image super-resolution | |
Kuruba et al. | Robust blood vessel detection with image enhancement using relative intensity order transformation and deep learning | |
Ramlan et al. | Image Enhancement Methods for Anterior Segment Photographed 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 |