CN1367468A - 基于k-均值聚类及信息融合的角膜病灶图像分割方法 - Google Patents
基于k-均值聚类及信息融合的角膜病灶图像分割方法 Download PDFInfo
- Publication number
- CN1367468A CN1367468A CN 02103794 CN02103794A CN1367468A CN 1367468 A CN1367468 A CN 1367468A CN 02103794 CN02103794 CN 02103794 CN 02103794 A CN02103794 A CN 02103794A CN 1367468 A CN1367468 A CN 1367468A
- Authority
- CN
- China
- Prior art keywords
- image
- cluster
- cornea
- brightness
- focus
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 210000004087 cornea Anatomy 0.000 title claims description 60
- 238000005267 amalgamation Methods 0.000 title 1
- 238000003709 image segmentation Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000003708 edge detection Methods 0.000 claims abstract description 4
- 230000003902 lesion Effects 0.000 claims description 37
- 230000004927 fusion Effects 0.000 claims description 13
- 238000005286 illumination Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 7
- 208000025865 Ulcer Diseases 0.000 claims description 5
- 231100000397 ulcer Toxicity 0.000 claims description 5
- CRHLEZORXKQUEI-UHFFFAOYSA-N dialuminum;cobalt(2+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Al+3].[Co+2].[Co+2] CRHLEZORXKQUEI-UHFFFAOYSA-N 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000011664 signaling Effects 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 2
- 238000007429 general method Methods 0.000 claims description 2
- 238000009499 grossing Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 abstract description 29
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000003745 diagnosis Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 8
- 238000000638 solvent extraction Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- 206010011026 Corneal lesion Diseases 0.000 description 3
- 239000000872 buffer Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 229910017052 cobalt Inorganic materials 0.000 description 2
- 239000010941 cobalt Substances 0.000 description 2
- 238000004040 coloring Methods 0.000 description 2
- GNBHRKFJIUUOQI-UHFFFAOYSA-N fluorescein Chemical compound O1C(=O)C2=CC=CC=C2C21C1=CC=C(O)C=C1OC1=CC(O)=CC=C21 GNBHRKFJIUUOQI-UHFFFAOYSA-N 0.000 description 2
- 238000003064 k means clustering Methods 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 150000001868 cobalt Chemical class 0.000 description 1
- GUTLYIVDDKVIGB-UHFFFAOYSA-N cobalt atom Chemical compound [Co] GUTLYIVDDKVIGB-UHFFFAOYSA-N 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
一种基于改进的k-均值聚类及信息融合的角膜病灶图像分割方法,涉及数字图像处理领域及眼病的诊断治疗,特征在于利用外眼角膜病灶图像特征的色度信息,通过改进的k-均值聚类完成图像的初次分割,即自学习确定初始色度聚类中心,再通过全图搜索确定最佳聚类中心,进行聚类分割;其次利用图像的亮度信息,通过改进的k-均值聚类完成图像的二次分割,即自学习确定初始亮度聚类中心,再通过全图搜索确定最佳聚类中心,进行聚类分割;利用一种信息融合技术,根据一种新的判断准则完成对色度和亮度信息分割结果的融合;利用数学形态学,设计结构算子及图像边缘检测算子,来完成对图像分割结果的优化。该方法复杂度低,分割效果理想,具有较好的鲁棒性。
Description
技术领域
本发明涉及医学图像领域及外眼病的临床诊断治疗,设计了一种结合外眼病图像特征,基于改进的k-均值聚类及信息融合的角膜病灶图像分割方法,将病灶从角膜病灶图像中分割出来以便于随后的特征分析。
背景技术
与角膜病灶图像分割相关的背景技术如下:
角膜病灶图像分割即角膜病灶图象中的角膜区域分割,是图象分割方法的一种具体应用。根据使用知识的特点与层次,通常的图象分割一般分为数据驱动与模型驱动两大类型。诸如边缘检测、区域聚类等算法。
分割算法中的区域聚类通常采用动态聚类法:先选择若干样品作为聚类中心,再按某种聚类准则。其中k-均值算法是一种典型的动态聚类法,此种算法能使聚类域中所有样品到聚类中心的距离平方和最小,其步骤如下所示:(1)选取k个聚类中心Z1 1,Z2 1,Λ,Zk 1:(上角标记为聚类中的迭代次数);(2)对于样本X(设进行到第k次迭代)
如果|X-Zj k|<|X-Zj k| 则X∈Sj k (i,j=1,2,Λ,k且i≠j) (1)
式中nj为Sj所包含的样本数;(4)如果Zj k+1≠Zj k,j=1,2,Λ,k则回到第二步,将全部样品重新分类,重新迭代计算;
如果Zj k+1=Zj k,j=1,2,Λ,k 则结束。
这种k-均值聚类算法的优点是:它能够动态聚类,具有一定的自适应性。但是,k-均值聚类的结果易受聚类中心的个数k及初始聚类中心的影响,同时也受样本的几何形状及排列次序的影响。因此,算法能否收敛取决于样本的特性和其能够形成不同区域的个数。
通常图像分割还需利用数学形态学进行结果优化。数学形态学的基本运算有膨胀、腐蚀、开、闭等,其中膨胀和腐蚀是两种最基本的运算,有直观的几何背景,根据它们可定义其它运算。
设A为原始图象,B是对原图像进行运算的结构元素,且Ac是A被b的平移,Ab表示A的补集,
表示B的反射。A被B膨胀的结果: A被B腐蚀的结果: 二值图象的开运算: foB=(fΘB)B (5)开运算的处理效果常表现为使目标轮廓平滑,消除毛刺和孤立点。二值图象的闭运算: f·B=(fB)ΘB (6)闭运算的处理效果常表现为使目标轮廓平滑,填平小沟,弥合空洞和裂缝。
彩色图像的分割通常需要进行彩色空间变换,彩色空间是以数值方式描述色彩的模型,如RGB空间,YIQ空间,HIS空间,CY空间等。由于人眼不能直接感知红绿蓝三色的比例,所以在RGB空间进行图象分割不符合人眼视觉。通常的图像分割算法都利用HIS空间或CY空间进行彩色变换。C-Y色空间与RGB色空间的转换关系如下(Y是亮度分量,R-Y和B-Y是色差分量): 在C-Y色空间内定义为:
饱和度分量S:
亮度分量Y: Y=0.299*R+0.587*G+0.114*B (10)其中θ表示彩色彼此相互区分的特性,不同波长的单色光具有不同的色调。而且,传统的彩色图像的分割算法很少利用信息融合的技术,只是简单地对图像的分割结果进行叠加。
发明内容
为了能将角膜病灶图像中的病灶自动分割出来,本发明充分考虑了角膜病灶图像的特点,利用图象处理中的现有技术加以创新,设计了一种新的分割方法,本发明基于上述技术的研究,利用改进的k-均值聚类、数学形态学及信息融合等方法,实现了角膜病灶图像的自动分割。
本发明的技术特征思路为:1、利用外眼角膜病灶图像特征的色度信息,通过改进的k-均值聚类完成图像的初次分割,即自学习确定初始色度聚类中心,然后再通过全图搜索确定最佳聚类中心,进行聚类分割;2、其次利用图像的亮度信息,通过改进的k-均值聚类完成图像的二次分割,即自学习确定初始亮度聚类中心,然后再通过全图搜索确定最佳聚类中心,进行聚类分割;3、利用一种信息融合技术,根据一种新的判断准则完成对色度和亮度信息分割结果的融合;4、利用数学形态学,设计结构算子及图像边缘检测算子,来完成对图像分割结果的优化;
本发明的技术方案参见图1、图2。这种基于改进的k-均值聚类及信息融合的角膜病灶图像分割方法,是由数码相机完成采集角膜病灶图像(一般指钴兰光病灶图像),并将角膜病灶图像的光学信号转换为电信号图象输入到计算机进行处理、传输等操作,计算机处理主要是通过USB接口软件、在k-均值聚类及数学形态学的基础上对角膜病灶图像进行读/写处理,处理后的角膜病灶图像或对角膜病灶图象进行分割处理后输出到缓存器,经显示器显示结果,本发明的特征在于它还依次包括下述步骤:(1)计算机从USB接口读入角膜病灶图象信号,并保存在内存中;(2)采用通用的直方图均衡的方法对图像进行预处理,由于该系统中的裂隙灯光源条件的限制,拍摄出的角膜病灶图像受光照不均的影响较大,因此在对角膜病灶图像进行病灶分割之前,必须对图像进行预处理;(3)进行色度空间转化,将图像从RGB空间转换到CY空间来得到其相应的色度图像和亮度图像;由于角膜病灶图像是彩色图象,而彩色图象与灰度图象的分割有着极大的差别:单纯利用颜色信息对彩色图象进行分割,在光照均匀的前提下,分割效果较好,但在光照不均或光照条件差的条件下,效果较差;而仅利用亮度信息对其分割,会丢失相关彩色信息,造成误分。所以本发明充分利用颜色和亮度信息完成图像分割。同时对于彩色图像而言,由于人眼不能直接感知红绿蓝三色的比例,所以在RGB空间进行图象分割不符合人眼视觉。CY空间和HIS空间都定义了颜色属性的色度、亮度及饱和度,与人类的视觉感知紧密相关;但是由于从RGB空间到HIS空间的转变较为复杂,而从RGB空间到CY空间的转变较为简单,而且速度很快。因此图像在CY空间当中进行色彩分割,效果会更好;(4)对色度图像进行角膜病灶的分割,其步骤为:设角膜病灶图像的大小为M*N(其中M为图像的宽度,N为图像的高度),同时设图像的像素为所需要聚类的样本,根据其特点,可将样本分为两类:A、通过自学习确定两个初始聚类中心Z1 1,Z2 1,设角膜病灶图像的像素为X,其三基色为R、G、B,其色度值为Xu,其计算参见背景技术中的公式9,如果G>R&&G>B则X∈S1 1,且S1类中心的矢量
如果B>G&&B>R则X∈S2 1,且S2类中心的矢量
接下来再利用传统的k-均值聚类方法进行聚类:B、如果|Xu-Z1 k|<|Xu-Z2 k| 则X∈S1 k,C、计算两个聚类中心的新向量值, D、如果|Zj k+1-Zj k|>ε j=1,2, 则回到B,重新迭代计算,
如果|Zj k+1-Zj k|≤ε j=1,2, 则结束;
分割的主要思路见图3。因角膜病灶图象(一般指钴兰光病灶图像)有其独特的特征:①图象的背景大都偏蓝;②所分割的病灶颜色以黄绿色为主;③图象的颜色信息受亮度影响较大;④病灶的形状大都是封闭的、不规则的。并且角膜病灶图像可分为病灶和非病灶两部分,即目标和背景两部分。所以k-均值算法的聚类中心应为:k=2。同时本发明利用这些已有的先验知识,对k-均值算法加以改进:
由于眼科病灶图象是卤素灯全光谱光源通过钴蓝滤色片拍摄而成,这种钴蓝滤色片对红光和绿光有一定的抑制作用,并且对于红光的抑制比绿光更为明显,同时由于滴染的荧光素作用,角膜中的病灶着染为绿色。所以眼科病灶图象以蓝绿分量为主,而红色分量低于蓝绿分量;(5)对亮度图像进行角膜病灶的分割;由于裂隙灯硬件拍摄条件的限制,拍摄出的图像受光照的影响很大,在对亮度进行聚类之前,应对其先进行亮度直方图均衡处理,通过亮度直方图均衡处理过的图像对比度和亮度都得到了明显的改善,这样便于以后的图像分割;经过直方图均衡处理过的图象,再对其亮度进行k-均值聚类,聚类时只有第一步自学习确定两个初始聚类中心与色度的聚类过程不同,其他都相似,
如果Y<aver(Y)+ε (ε=10)则
自“搜索全角膜的每个像素的亮度值”后的其他步骤与色度的聚类过程相似;(6)将色彩分割和亮度分割这两种方法得到的结果进行信息融合:
设外眼图像的角膜区域为全集W,色度的k均值聚类结果为A(AW),亮度的k均值聚类结果为B(BW),最终的聚类结果为S,在光照条件均匀的情况下,色度和亮度的聚类结果都比较理想,S应为二者的交集,但是在光照条件不理想的情况下,色度的聚类结果失真较大,亮度的聚类结果较为理想,故S应为亮度的聚类结果,而不应取二者的交集,(7)利用数学形态学对得到的图像结果进行优化处理;
根据眼科病灶图象特点选择十字形的4-连通结构算子B1,3×3方形的8-连通结构算子B2,设原始图像为f,设计了剔除噪声点的数学形态学的边缘检测算子:f·B-(fΘB1)ΘB1即(fB1)ΘB2-(fΘB1)ΘB1,
利用该公式得到分割图象的边缘,由于溃疡病灶基本都是闭合的,所以边缘如果是线状、枝状等非闭合形状,应全部剔除,接下来填充图象,并对图象进行开运算,使分割图象轮廓平滑,消除毛刺和孤立点,
再对图象进行闭运算,填补分割图象中的空洞和裂缝;(8)输出并显示图像分割结果。
图6利用图像的分割结果描述了5-8的分割过程。
本发明提出的方法复杂度低,能够满足分割角膜病灶图象的要求,分割效果较理想,而且在光照不均匀的条件下也能得到理想的结果,具有较好的鲁棒性。
附图说明
图1是外眼病灶图像分割系统框图
1、数码相机,2、USB接口,3、计算机处理器,4、输出缓存,5、角膜病灶分割,6、显示器,7、分割结果;
图2是本发明方法主程序流程图;
图3是本发明方法中色度的k均值聚类子程序流程图;
图4是本发明方法中亮度的k均值聚类子程序流程图;
图5是本发明中信息融合得到分割结果S子程序流程图;
图6是角膜病灶图像的分割过程示例;
图中(a)角膜病灶原始图像,(b)色度的k-均值聚类分割结果A,(c)亮度的k-均值聚类分割结果B,(d)A和B进行信息融合的聚类结果S,(e)进行数学形态学运算后的最终分割结果;
图6(a)为输入的角膜病灶原始图象,白色圆圈所包括的范围为角膜,目的是将病灶从角膜当中分割出来;其中溃疡(如B区),也就是感兴趣区病灶,须将其分割;而周围的过渡区(如A区),并不是溃疡,不需要分割;
图7是本发明中一幅角膜病灶图像及其分割结果;
图中(a)角膜病灶原始图像,(b)分割的病灶与角膜对比图,(c)病灶分割结果;
其中图b中的三个白色圆圈是角膜当中描述病灶位置的三个常用区域,两条红色的直线分别为角膜的横轴和纵轴,绿线则是描述病灶的形状轮廓。
具体实施方式
本发明的实施方案参见图1-5,原始角膜病灶图象可以是通过数码相机实时采集到的图象,也可以是实现通过数码相机采集到保存在计算机硬盘里的图象。
图1中数码相机和USB接口都是市售的,主要完成采集角膜病灶图象,将角膜病灶图象的光学信号转换为电信号图象输入到计算机,便于计算机处理、传输等操作。计算机处理主要是通过USB接口软件对角膜病灶图象进行读/写处理,处理后的角膜病灶图象输出到缓存器,便于显示。显示器是图象的输出设备,人眼可以通过显示器观看原始角膜病灶图象和分割结果图象,角膜病灶分割是对计算机读入的角膜病灶图象进行分割处理,输出分割结果。
角膜病灶分割主要通过软件来实现。下面结合实例详细描述角膜病灶分割的过程。我们已通过新的外眼图像采集系统采集了150幅角膜病灶图像,并存在计算机硬盘当中,其程序流程图可参见方法流程图,步骤如下:
第一步:首先向病眼滴染荧光素,这样可以着染溃疡病灶,其次在钴兰光的条件下,利用裂隙灯聚焦,然后通过光电接口用数码相机对病人的眼睛进行拍摄,然后通过USB接口传到计算机硬盘当中,这样完成了角膜病灶图像采集过程;
第二步:对角膜病灶图像进行图像预处理,采用传统的直方图均衡,提高图像的对比度和亮度,便于接下来的色度和亮度图像的病灶分割;
第三步:利用已有的色度空间变换公式(见上),将角膜病灶图像从RGB空间转换到CY空间得到色度图像和亮度图像;
第四步:用改进的k-均值聚类对亮度图像进行病灶分割;具体思路见上述方法简介中的步骤(4),根据色度信息,确定初始聚类中心。再通过传统的k-均值聚类算法进行聚类得到初次分割图像1;
第五步:利用改进的k-均值聚类对亮度图像进行病灶分割,首先对图像进行亮度直方图均衡预处理。然后对预处理过的图像进行自学习聚类,首先计算图像的平均亮度,将大于平均亮度的像素值作为集合1进行聚类病计算聚类中心1,将小于平均亮度的像素值作为集合2进行聚类病计算聚类中心2。然后再通过全图搜索进行聚类,确定最佳聚类中心,通过上述方法得到分割的病灶图像2;
第六步:将上述方法得到的两幅图像进行信息融合,得到最佳结果。如果在光照均匀的情况下,融合图像取两者的交集;如果在光照不均的情况下,由于利用色度信息分割的图像1效果较差,融合图像应为亮度分割结果,即病灶图像2;
第七步:利用数学形态学对融合图像进行优化——使融合图象轮廓平滑,消除毛刺和孤立点并填补分割图象中的空洞和裂缝;
第八步:得到分割图像结果,见图7;
第九步:得到最佳病灶图像之后,计算病灶中所包含的像素数,同时计算角膜中所包含的像素数,计算病灶占角膜中的百分比。
Claims (1)
1、一种基于k-均值聚类及信息融合的角膜病灶图像分割方法,是由数码相机完成采集角膜病灶图像(一般指钴兰光病灶图像),并将角膜病灶图像的光学信号转换为电信号图象输入到计算机进行处理,本发明的特征在于它还依次包括下述步骤:(1)计算机从USB接口读入角膜病灶图象信号,并保存在内存中;(2)采用通用的直方图均衡的方法对图像进行预处理;(3)进行色度空间转化,将图像从RGB空间转换到CY空间来得到其相应的色度图像和亮度图像;(4)对色度图像进行角膜病灶的分割,其步骤为:设角膜病灶图像的大小为M*N(其中M为图像的宽度,N为图像的高度),同时设图像的像素为所需要聚类的样本,根据其特点,可将样本分为两类:A、通过自学习确定两个初始聚类中心Z1 1,Z2 1:设角膜病灶图像的像素为X,其三基色为R、G、B,其色度值为Xu:如果G>R&&G>B则X∈S1 1,且S1类中心的矢量
如果B>G&&B>R则X∈S2 1,且S2类中心的矢量
接下来再利用传统的k-均值聚类方法进行聚类:B、如果|Xu-Z1 k|<|Xu-Z2 k| 则X∈S1 k,C、计算两个聚类中心的新向量值 D、如果|Zj k+1-Zj k|>ε j=1,2, 则回到B,重新迭代计算,如果|Zj k+1-Zj k|≤ε j=1,2, 则结束;(5)对亮度图像进行角膜病灶的分割,分割过程如下:
在对亮度进行聚类之前,应对其先进行亮度直方图均衡处理,经过直方图均衡处理过的图象,再对其亮度进行k-均值聚类,聚类时只有第一步自学习确定两个初始聚类中心与色度的聚类过程不同,其他都相似,
通过自学习确定两个初始聚类中心Z1 1,Z2 1,
自“搜索全角膜的每个像素的亮度值”后的其他步骤与色度的聚类过程相似;(6)将色彩分割和亮度分割这两种方法得到的结果进行信息融合:
设外眼图像的角膜区域为全集W,色度的k均值聚类结果为A(AW),亮度的k均值聚类结果为B(BW),最终的聚类结果为S,在光照条件均匀的情况下,色度和亮度的聚类结果都比较理想,S应为二者的交集,但是在光照条件不理想的情况下,色度的聚类结果失真较大,亮度的聚类结果较为理想,故S应为亮度的聚类结果,而不应取二者的交集,(7)利用数学形态学对得到的图像结果进行优化处理,根据眼科病灶图象特点选择十字形的4-连通结构算子B1,3×3方形的8-连通结构算子B2,设原始图像为f,设计了剔除噪声点的数学形态学的边缘检测算子:f·B-(fΘB1)ΘB1即(fB1)ΘB2-(fΘB1)ΘB1,
利用该公式得到分割图象的边缘,由于溃疡病灶基本都是闭合的,所以边缘如果是线状、枝状等非闭合形状,应全部剔除,接下来填充图象,并对图象进行开运算,使分割图象轮廓平滑,消除毛刺和孤立点,再对图象进行闭运算,填补分割图象中的空洞和裂缝;(8)输出并显示图像分割结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021037949A CN1139898C (zh) | 2002-03-25 | 2002-03-25 | 基于k-均值聚类及信息融合的角膜病灶图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021037949A CN1139898C (zh) | 2002-03-25 | 2002-03-25 | 基于k-均值聚类及信息融合的角膜病灶图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1367468A true CN1367468A (zh) | 2002-09-04 |
CN1139898C CN1139898C (zh) | 2004-02-25 |
Family
ID=4739900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB021037949A Expired - Fee Related CN1139898C (zh) | 2002-03-25 | 2002-03-25 | 基于k-均值聚类及信息融合的角膜病灶图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1139898C (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1313964C (zh) * | 2004-07-05 | 2007-05-02 | 南京大学 | 基于聚类学习器集成的数字图像分割方法 |
CN100359532C (zh) * | 2003-05-29 | 2008-01-02 | 上海交通大学 | 多目标图像的分层聚类方法 |
CN1875378B (zh) * | 2003-11-12 | 2010-05-05 | 英国电讯有限公司 | 图像中的对象检测 |
CN101216886B (zh) * | 2008-01-11 | 2010-06-09 | 北京航空航天大学 | 一种基于谱分割理论的镜头聚类方法 |
CN101923712A (zh) * | 2010-08-03 | 2010-12-22 | 苏州大学 | 基于粒子群优化的K-means聚类算法的基因芯片图像分割法 |
CN101719277B (zh) * | 2009-12-31 | 2011-11-30 | 华中科技大学 | 一种遗传模糊聚类图像分割方法 |
CN101694720B (zh) * | 2009-10-13 | 2012-02-08 | 西安电子科技大学 | 基于空间关联条件概率融合的多时相sar图像变化检测方法 |
CN101203185B (zh) * | 2005-04-15 | 2012-05-23 | R·戈多伊斯 | 角膜内皮细胞分析样本的统计学检验方法 |
CN102496146A (zh) * | 2011-11-28 | 2012-06-13 | 南京大学 | 一种基于视觉共生的图像分割方法 |
CN102608116A (zh) * | 2012-01-30 | 2012-07-25 | 辽宁中医药大学 | 一种基于灰度图像分析的胃溃疡模型评价方法 |
CN101665210B (zh) * | 2008-09-03 | 2013-05-15 | 株式会社日立制作所 | 图像识别装置 |
CN103292741A (zh) * | 2013-05-29 | 2013-09-11 | 哈尔滨工程大学 | 一种基于k均值颜色聚类的物体表面三维轮廓的结构光视觉测量方法 |
CN104063599A (zh) * | 2014-06-26 | 2014-09-24 | 江苏中烟工业有限责任公司 | 一种用于评价烟叶质量的指标筛选及处理方法 |
CN104166400A (zh) * | 2014-07-11 | 2014-11-26 | 杭州精久科技有限公司 | 基于多传感器融合的视觉导引agv系统 |
CN105096322A (zh) * | 2015-07-26 | 2015-11-25 | 郭新 | 基于谱聚类的边缘检测方法 |
CN107122597A (zh) * | 2017-04-12 | 2017-09-01 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种角膜受损智能诊断系统 |
CN109064473A (zh) * | 2018-07-26 | 2018-12-21 | 华南理工大学 | 一种2.5d超声宽景图像分割方法 |
CN109949906A (zh) * | 2019-03-22 | 2019-06-28 | 上海鹰瞳医疗科技有限公司 | 病理切片图像处理和模型训练方法及设备 |
CN111652841A (zh) * | 2020-04-24 | 2020-09-11 | 温州医科大学附属眼视光医院 | 一种基于图像处理的角膜损伤检测方法 |
CN112132153A (zh) * | 2020-09-22 | 2020-12-25 | 山东大学 | 基于聚类和形态学处理的番茄果实识别方法及系统 |
CN117422717A (zh) * | 2023-12-19 | 2024-01-19 | 长沙韶光芯材科技有限公司 | 一种掩膜板污渍智能定位方法和系统 |
-
2002
- 2002-03-25 CN CNB021037949A patent/CN1139898C/zh not_active Expired - Fee Related
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100359532C (zh) * | 2003-05-29 | 2008-01-02 | 上海交通大学 | 多目标图像的分层聚类方法 |
CN1875378B (zh) * | 2003-11-12 | 2010-05-05 | 英国电讯有限公司 | 图像中的对象检测 |
CN1313964C (zh) * | 2004-07-05 | 2007-05-02 | 南京大学 | 基于聚类学习器集成的数字图像分割方法 |
CN101203185B (zh) * | 2005-04-15 | 2012-05-23 | R·戈多伊斯 | 角膜内皮细胞分析样本的统计学检验方法 |
CN101216886B (zh) * | 2008-01-11 | 2010-06-09 | 北京航空航天大学 | 一种基于谱分割理论的镜头聚类方法 |
CN101665210B (zh) * | 2008-09-03 | 2013-05-15 | 株式会社日立制作所 | 图像识别装置 |
CN101694720B (zh) * | 2009-10-13 | 2012-02-08 | 西安电子科技大学 | 基于空间关联条件概率融合的多时相sar图像变化检测方法 |
CN101719277B (zh) * | 2009-12-31 | 2011-11-30 | 华中科技大学 | 一种遗传模糊聚类图像分割方法 |
CN101923712A (zh) * | 2010-08-03 | 2010-12-22 | 苏州大学 | 基于粒子群优化的K-means聚类算法的基因芯片图像分割法 |
CN102496146A (zh) * | 2011-11-28 | 2012-06-13 | 南京大学 | 一种基于视觉共生的图像分割方法 |
CN102496146B (zh) * | 2011-11-28 | 2014-03-05 | 南京大学 | 一种基于视觉共生的图像分割方法 |
CN102608116A (zh) * | 2012-01-30 | 2012-07-25 | 辽宁中医药大学 | 一种基于灰度图像分析的胃溃疡模型评价方法 |
CN103292741A (zh) * | 2013-05-29 | 2013-09-11 | 哈尔滨工程大学 | 一种基于k均值颜色聚类的物体表面三维轮廓的结构光视觉测量方法 |
CN103292741B (zh) * | 2013-05-29 | 2015-12-23 | 哈尔滨工程大学 | 一种基于k均值颜色聚类的物体表面三维轮廓的结构光视觉测量方法 |
CN104063599A (zh) * | 2014-06-26 | 2014-09-24 | 江苏中烟工业有限责任公司 | 一种用于评价烟叶质量的指标筛选及处理方法 |
CN104063599B (zh) * | 2014-06-26 | 2017-02-15 | 江苏中烟工业有限责任公司 | 一种用于评价烟叶质量的指标筛选及处理方法 |
CN104166400A (zh) * | 2014-07-11 | 2014-11-26 | 杭州精久科技有限公司 | 基于多传感器融合的视觉导引agv系统 |
CN104166400B (zh) * | 2014-07-11 | 2017-02-22 | 杭州精久科技有限公司 | 基于多传感器融合的视觉导引agv系统 |
CN105096322A (zh) * | 2015-07-26 | 2015-11-25 | 郭新 | 基于谱聚类的边缘检测方法 |
CN105096322B (zh) * | 2015-07-26 | 2018-11-23 | 郭新 | 基于谱聚类的边缘检测方法 |
CN107122597A (zh) * | 2017-04-12 | 2017-09-01 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种角膜受损智能诊断系统 |
CN107122597B (zh) * | 2017-04-12 | 2020-04-24 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种角膜受损智能诊断系统 |
CN109064473A (zh) * | 2018-07-26 | 2018-12-21 | 华南理工大学 | 一种2.5d超声宽景图像分割方法 |
CN109064473B (zh) * | 2018-07-26 | 2022-03-29 | 华南理工大学 | 一种2.5d超声宽景图像分割方法 |
CN109949906A (zh) * | 2019-03-22 | 2019-06-28 | 上海鹰瞳医疗科技有限公司 | 病理切片图像处理和模型训练方法及设备 |
CN111652841A (zh) * | 2020-04-24 | 2020-09-11 | 温州医科大学附属眼视光医院 | 一种基于图像处理的角膜损伤检测方法 |
CN112132153A (zh) * | 2020-09-22 | 2020-12-25 | 山东大学 | 基于聚类和形态学处理的番茄果实识别方法及系统 |
CN117422717A (zh) * | 2023-12-19 | 2024-01-19 | 长沙韶光芯材科技有限公司 | 一种掩膜板污渍智能定位方法和系统 |
CN117422717B (zh) * | 2023-12-19 | 2024-02-23 | 长沙韶光芯材科技有限公司 | 一种掩膜板污渍智能定位方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN1139898C (zh) | 2004-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1139898C (zh) | 基于k-均值聚类及信息融合的角膜病灶图像分割方法 | |
Hsu et al. | The role of domain knowledge in the detection of retinal hard exudates | |
JP4649532B2 (ja) | 認知光源とハイライトとを利用した照明色度変換装置及びそのための方法並びにそのためのプログラムを記録したコンピュータが読める記録媒体 | |
US7068840B2 (en) | Determination of an illuminant of digital color image by segmentation and filtering | |
CN104537373B (zh) | 舌下络脉诊断用多光谱舌下图像特征提取方法 | |
CN106934418B (zh) | 一种基于卷积递归网络的绝缘子红外诊断方法 | |
CN106388781A (zh) | 一种皮肤肤色及其色素沉淀情况的检测方法 | |
CN105030200A (zh) | 一种基于绿色单色光源的中医舌质舌苔分离方法 | |
CN109584210B (zh) | 多光谱三维静脉显像系统 | |
CN112598761B (zh) | 一种纺织服装色彩搭配协同设计与智能推荐方法 | |
WO2023142615A1 (zh) | 图像处理方法、装置、设备、可读存储介质及程序产品 | |
US20130342694A1 (en) | Method and system for use of intrinsic images in an automotive driver-vehicle-assistance device | |
CN1162798C (zh) | 基于多类支持向量机的中医舌色、苔色、舌苔厚度分析方法 | |
Bai et al. | Automatic segmentation of cervical region in colposcopic images using K-means | |
CN110246116A (zh) | 数字病理切片由he染色到ihc染色的计算机自动生成方法 | |
JP3200705B2 (ja) | 色座標空間構造、色量子化及可変グレー領域の設定方法 | |
CN105701806A (zh) | 基于深度图像的帕金森震颤运动特征检测方法及系统 | |
CN115082451A (zh) | 一种基于图像处理的不锈钢汤勺缺陷检测方法 | |
CN111415352A (zh) | 一种基于深度级联网络的癌转移全景病理切片分析方法 | |
CN107463913B (zh) | 基于Lab的图像胶原纤维识别计算系统及方法 | |
CN100405402C (zh) | 舌头彩色数码照片的舌苔的分割提取方法 | |
CN110223253B (zh) | 一种基于图像增强的去雾方法 | |
CN106204476B (zh) | 一种低照度彩色图像增强方法 | |
CN116542888A (zh) | 一种多聚焦彩色图像获取方法及产品 | |
CN114511567B (zh) | 舌体与舌苔图像识别分离方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
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 | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |