CN102903104A - 一种基于减法聚类的快速图像分割方法 - Google Patents
一种基于减法聚类的快速图像分割方法 Download PDFInfo
- Publication number
- CN102903104A CN102903104A CN2012103378380A CN201210337838A CN102903104A CN 102903104 A CN102903104 A CN 102903104A CN 2012103378380 A CN2012103378380 A CN 2012103378380A CN 201210337838 A CN201210337838 A CN 201210337838A CN 102903104 A CN102903104 A CN 102903104A
- Authority
- CN
- China
- Prior art keywords
- pixel
- pixels
- sampled pixel
- weight matrix
- density
- 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
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于减法聚类的快速图像分割方法。本发明首先将所有像素点归一化到一个超立方体中,对等待聚类的所有像素进行等间隔均匀采样后重组;在重组的像素中,计算采样像素点两两之间的密度权值矩阵及其逆阵以及采样像素与剩余未采样像素之间的密度权值矩阵。然后计算逼近的未采样像素两两之间的密度权值矩阵和计算所有像素点的密度值。最后计算所有像素的最大密度值并获得聚类中心,为找出新的聚类中心,需对每个像素点的密度值进行衰减,该过程不断迭代,根据终止条件停止迭代。本发明与经典的减法聚类方法相比,本发明在不影响聚类结果的情况下,对于较大规模数据集,大大提高减法聚类方法的实时性。
Description
技术领域
本发明涉及数字图像处理技术领域,具体来说是一种基于减法聚类的快速图像分割方法。
背景技术
数字图像分割是数字图像处理领域的一块重要技术。图像分割就是将图像细分为构成它的子区域或对象。常常被当成一个重要环节被嵌入机器视觉系统的处理过程中。k均值聚类、模糊c均值聚类经常被用于数字图像的分割,是大家常用的方法。然而由于k均值聚类和模糊c均值聚类这两种图像分割方法对初始的聚类中心很敏感,不同的初始值很容易得到不同的结果。由于减法聚类具有获取初始聚类中心的特点,所以可以和k均值聚类或模糊c均值聚类方法结合。因此,在对图像进行分割前,可首先用减法聚类方法获得初始图像的像素聚类中心,然后再结合k均值聚类、模糊c均值聚类对图像进行分割。近几年,减法聚类比较常见的应用是与模糊c均值聚类结合,结合具体应用对数据进行聚类并对减法聚类做相应的改进。
然而,现在的数字图像很多都是高清,数据量大,分割的时间复杂度高,所以在用减法聚类对数字图像进行分割前,需考虑一种能显著降低时间复杂度的数字图像分割方法。
发明内容
本发明针对数字图像分割时间复杂度高的问题,提出了一种基于减法聚类的快速图像分割方法。
一种基于减法聚类的快速图像分割方法,包括:
步骤(1)将所有像素点进行归一化,对等待分割的所有像素进行等间隔均匀采样后重组。
步骤(2)在重组的像素中,计算采样像素点两两之间的密度权值矩阵及其逆矩阵、采样像素与剩余未采样像素之间的密度权值矩阵。
步骤(3)根据奈斯特龙逼近原理,获得未采样像素两两之间密度权值矩阵的逼近矩阵。
步骤(4)用步骤(3)生成的未采样像素两两之间的密度权值矩阵、采样像素点两两之间的密度权值矩阵以及采样像素与剩余像素之间的密度权值矩阵计算所有像素点的密度值。
步骤(5)利用步骤(4)得到的所有像素的密度值,计算所有像素的最大密度值并获得聚类中心,为找出新的聚类中心,根据减法聚类,对每个像素点的密度值进行衰减,该过程不断迭代,根据终止条件停止迭代,从而完成图像分割。
附图说明
图1为本发明的方法流程图;
图2为本发明图像分割实验结果。其中:(a) Hestain测试图 (b) Fabric测试图 (c) Hestain经典减法聚类分割图 (d) Fabric经典减法聚类分割图 (e) Hestain本发明分割图 (f) Fabric本发明分割图。
具体实施方式
以下结合附图对本发明作进一步说明。
如图1所示,一种基于减法聚类的快速图像分割方法,包括:
(1)将所有像素点归一化到一个超立方体中,对等待聚类的所有像素进行等间隔均匀采样后重组。
(2)在重组的像素中,计算采样像素点两两之间的密度权值矩阵及其逆阵以及采样像素与剩余未采样像素之间的密度权值矩阵。
(3)根据奈斯特龙逼近原理,将计算逼近的未采样像素两两之间的密度权值矩阵,所述的奈斯特龙逼近原理由Fowlkes C.等人提出。
(4)用(3)生成的未采样像素两两之间的密度权值矩阵,采样像素点两两之间的密度权值矩阵以及采样像素与剩余像素之间的密度权值矩阵计算所有像素点的密度值。
(5)利用(4)得到的所有像素的密度值,计算所有像素的最大密度值并获得聚类中心,为找出新的聚类中心,需对每个像素点的密度值进行衰减,该过程不断迭代,根据终止条件停止迭代。
将所有像素点归一化到一个超立方体中,对等待聚类的所有像素进行等间隔均匀采样后重组的具体过程包括以下步骤:
1.计算所有像素各个维度的最大最小值;
2.利用1)计算的最大最小值归一化所有像素;
利用获得的不同维度的最大最小值对所有像素进行归一化,使得所有像素的不同维度值归一化到[0,1]区间。
3.对归一化后的所有像素进行均匀采样;
4.归一化后的像素经采样后分为采样像素和未采样像素两部分;
5.将所有采样像素放在新像素集的前面,未采样像素放在新像素集的后面。
在重组的像素中,计算采样像素点两两之间的密度权值矩阵及其逆阵以及采样像素与剩余未采样像素之间的密度权值矩阵的具体过程包括以下步骤:
1.计算采样像素两两之间的密度权值矩阵;
(1)
半径定义了该像素点在不同维度一个邻域,即影响半径。
2.求取采样像素两两之间的密度权值矩阵的逆矩阵;
3.计算采样像素与剩余未采样像素之间的密度权值矩阵;
根据奈斯特龙逼近原理,将计算逼近的未采样像素两两之间的密度权值矩阵的具体过程包括以下步骤:
1.求采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵;
2.将采样像素与剩余未采样像素之间的密度权值矩阵的转置矩阵乘于采样像素两两之间的密度权值矩阵的逆矩阵;
3.将计算结果乘于采样像素与剩余未采样像素之间的密度权值矩阵,从而获得未采样像素两两之间密度权值矩阵的逼近矩阵,计算公式如式(2):
(2)
用生成的未采样像素两两之间的密度权值矩阵,采样像素点两两之间的密度权值矩阵以及采样像素与剩余像素之间的密度权值矩阵计算所有像素点的密度值的具体过程包括以下步骤:
1.将采样像素两两之间的密度权值矩阵、采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵、采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵及未采样像素两两之间密度权值矩阵的逼近矩阵进行组合生成所有采样像素两两之间的密度权值矩阵;
即利用公式(3)生成所有采样像素两两之间的密度权值矩阵:
2.将所有采样像素两两之间的密度权值矩阵的每行元素进行求和;
3.将2)的求和结果映射为所有像素的密度值。
利用得到的所有像素的密度值,计算所有像素的最大密度值并获得聚类中心,为找出新的聚类中心,需对每个像素点的密度值进行衰减,该过程不断迭代,根据终止条件停止迭代。并在k个初始聚类中心的基础上,进一步用k均值聚类对图像进行分割的具体过程包括以下步骤:
1.计算所有像素的最大密度值,从而获取第一个聚类中心;
对式(5)进行排序,计算所有像素的最大密度值,选择最大密度值点为第1个聚类中心:
2.对所有像素的密度值进行衰减,并计算所有像素的最大密度值,从而获取第二个聚类中心;
(6)
3.上述的过程不断迭代,根据终止条件停止迭代。
4.并在k个初始聚类中心的基础上,进一步用公开使用的k均值聚类技术对图像做进一步分割。
为了说明本发明针对图像的分割效果及时间复杂度性能,针对标准图像进行分割并做了时间复杂度性能对比以验证本发明的实用性。
图2(a)名为Hestain,图像格式1521143,类别数为3。图2(b)名为Fabric,图像格式3202403,类别数为6。图像原为R-G-B空间的图像,本发明聚类时将其转至L-a-b空间,并对图像像素L-a-b空间中的a、b色度分量进行聚类。因此,调整后Hestain图像像素总个数为17328,维度为2。同理,Fabric图像像素总个数为76800,维度为2。
对于Hestain图像,像素总数为17328,采样间隔为50,采样像素数为347,剩余像素数为16981,和设置为0.3,取1.25,取0.5。对于Fabric图像,像素总数为76800个,采样间隔为500,像素数为154,剩余像素数为76646,和设置为0.15,取1.25,取0.5。图2给出了本发明方法与经典减法聚类方法的实验对比。从图2(c)与图2(e)及图2(d)与图2(f)的聚类结果来看,本发明的减法聚类方法与经典的减法聚类法对Hestain与Fabric两幅图像具有几乎一致的聚类效果。
表1和表2分别给出了对应图Hestain和图Fabric的分割参数。从聚类中心及其对应的密度值上看,个别聚类中心的密度值有细微差异,但本发明方法与经典减法聚类均具有同样的聚类中心。在分割效果一样的情况下,再看看两种分割方法的平均耗时情况(重复10次的平均值),可以发现本发明方法的加速效果非常明显,特别是对像素数达76800个的情况,性能加速比高达252.56。
Claims (6)
1.一种基于减法聚类的快速图像分割方法,其特征在于该方法包括如下步骤:
步骤(1)将所有像素点进行归一化,对等待分割的所有像素进行等间隔均匀采样后重组;
步骤(2)在重组的像素中,计算采样像素点两两之间的密度权值矩阵及其逆矩阵、采样像素与剩余未采样像素之间的密度权值矩阵;
步骤(3)根据奈斯特龙逼近原理,获得未采样像素两两之间密度权值矩阵的逼近矩阵;
步骤(4)用步骤(3)生成的未采样像素两两之间的密度权值矩阵、采样像素点两两之间的密度权值矩阵以及采样像素与剩余像素之间的密度权值矩阵计算所有像素点的密度值;
步骤(5)利用步骤(4)得到的所有像素的密度值,计算所有像素的最大密度值并获得聚类中心,为找出新的聚类中心,根据减法聚类,对每个像素点的密度值进行衰减,该过程不断迭代,根据终止条件停止迭代,从而完成图像分割。
2.根据权利要求1所述的一种基于减法聚类的快速图像分割方法,其特征在于:步骤(1)包括以下步骤:
1-1.计算所有像素各个维度的最大值和最小值;
1-2.利用1-1计算的最大值和最小值归一化所有像素;
1-3.对归一化后的所有像素进行均匀采样;
1-4.将均匀采样像素后分为采样像素和未采样像素两部分;
1-5.将所有采样像素放在新像素集的前面,未采样像素放在新像素集的后面。
4.根据权利要求1所述的一种基于减法聚类的快速图像分割方法,其特征在于:步骤(3)包括以下步骤:
3-1.求采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵;
3-2.将采样像素与剩余未采样像素之间的密度权值矩阵的转置矩阵乘于采样像素两两之间的密度权值矩阵的逆矩阵;
3-3.将3-2的计算结果乘于采样像素与剩余未采样像素之间的密度权值矩阵,从而获得未采样像素两两之间密度权值矩阵的逼近矩阵。
5.根据权利要求1所述的一种基于减法聚类的快速图像分割方法,其特征在于:步骤(4)包括以下步骤:
4-1.将采样像素两两之间的密度权值矩阵、采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵、采样像素与剩余未采样像素之间的密度权值矩阵转置矩阵及未采样像素两两之间密度权值矩阵的逼近矩阵进行组合生成所有采样像素两两之间的密度权值矩阵;
4-2.将所有采样像素两两之间的密度权值矩阵的每行元素进行求和;
4-3.将4-2的求和结果映射为所有像素的密度值。
6.根据权利要求1所述的一种基于减法聚类的快速图像分割方法,其特征在于:步骤(5)包括以下步骤:
5-1.计算所有像素的最大密度值,从而获取第一个聚类中心;
5-2.对所有像素的密度值进行衰减,并计算所有像素的最大密度值,从而获取第二个聚类中心;
5-3.5-2的过程不断迭代,根据终止条件停止迭代;
5-4.根据k均值方法对所有像素做进一步分割。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210337838.0A CN102903104B (zh) | 2012-09-13 | 2012-09-13 | 一种基于减法聚类的快速图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210337838.0A CN102903104B (zh) | 2012-09-13 | 2012-09-13 | 一种基于减法聚类的快速图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102903104A true CN102903104A (zh) | 2013-01-30 |
CN102903104B CN102903104B (zh) | 2015-04-29 |
Family
ID=47575317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210337838.0A Expired - Fee Related CN102903104B (zh) | 2012-09-13 | 2012-09-13 | 一种基于减法聚类的快速图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102903104B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196497A (zh) * | 2013-04-09 | 2013-07-10 | 厦门大学 | 一种基于核减法聚类的多功能传感器样本选择方法 |
CN103559716A (zh) * | 2013-11-12 | 2014-02-05 | 广州太普信息技术有限公司 | 一种缺陷图像的自动分割方法 |
CN106570015A (zh) * | 2015-10-09 | 2017-04-19 | 杭州海康威视数字技术股份有限公司 | 图像搜索方法和装置 |
CN112578248A (zh) * | 2020-12-26 | 2021-03-30 | 广东电网有限责任公司电力科学研究院 | 一种局部放电定位误差修正方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763514A (zh) * | 2010-01-15 | 2010-06-30 | 西安电子科技大学 | 基于特征重要度排序谱聚类的图像分割方法 |
CN102622756A (zh) * | 2012-03-14 | 2012-08-01 | 西安电子科技大学 | 基于全变分谱聚类的sar图像分割方法 |
-
2012
- 2012-09-13 CN CN201210337838.0A patent/CN102903104B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763514A (zh) * | 2010-01-15 | 2010-06-30 | 西安电子科技大学 | 基于特征重要度排序谱聚类的图像分割方法 |
CN102622756A (zh) * | 2012-03-14 | 2012-08-01 | 西安电子科技大学 | 基于全变分谱聚类的sar图像分割方法 |
Non-Patent Citations (3)
Title |
---|
孙志海 等: "基于减法聚类算法的视频运动目标定位", 《光电工程》 * |
尤全增: "基于Nystrom扩展的大规模谱聚类算法", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李俊英: "谱聚类方法在图像分割中的应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196497A (zh) * | 2013-04-09 | 2013-07-10 | 厦门大学 | 一种基于核减法聚类的多功能传感器样本选择方法 |
CN103559716A (zh) * | 2013-11-12 | 2014-02-05 | 广州太普信息技术有限公司 | 一种缺陷图像的自动分割方法 |
CN106570015A (zh) * | 2015-10-09 | 2017-04-19 | 杭州海康威视数字技术股份有限公司 | 图像搜索方法和装置 |
CN106570015B (zh) * | 2015-10-09 | 2020-02-21 | 杭州海康威视数字技术股份有限公司 | 图像搜索方法和装置 |
CN112578248A (zh) * | 2020-12-26 | 2021-03-30 | 广东电网有限责任公司电力科学研究院 | 一种局部放电定位误差修正方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102903104B (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111199214B (zh) | 一种残差网络多光谱图像地物分类方法 | |
CN109993748B (zh) | 一种基于点云处理网络的三维网格物体分割方法 | |
CN110322453B (zh) | 基于位置注意力和辅助网络的3d点云语义分割方法 | |
CN104023230B (zh) | 一种基于梯度关联性的无参考图像质量评价方法 | |
CN109934826A (zh) | 一种基于图卷积网络的图像特征分割方法 | |
CN109886937B (zh) | 基于超像素分割图像识别的绝缘子缺陷检测方法 | |
CN111339924B (zh) | 一种基于超像素和全卷积网络的极化sar图像分类方法 | |
CN111091147B (zh) | 一种图像分类方法、装置及设备 | |
CN102903104A (zh) | 一种基于减法聚类的快速图像分割方法 | |
CN114881861B (zh) | 基于双采样纹理感知蒸馏学习的不均衡图像超分方法 | |
CN105959175A (zh) | 基于GPU加速的kNN算法的网络流量分类方法 | |
CN111695624A (zh) | 数据增强策略的更新方法、装置、设备及存储介质 | |
CN110458189A (zh) | 压缩感知和深度卷积神经网络电能质量扰动分类方法 | |
CN105913451B (zh) | 一种基于图模型的自然图像超像素分割方法 | |
CN112215199A (zh) | 基于多感受野与密集特征聚合网络的sar图像舰船检测方法 | |
CN113963009A (zh) | 基于可形变划块的局部自注意力的图像处理方法和模型 | |
CN103927765B (zh) | 用于定位图像中条码区域的方法和装置 | |
CN113673465A (zh) | 图像检测方法、装置、设备及可读存储介质 | |
CN105868778A (zh) | 一种基于频域波动特征的中值滤波检测方法 | |
CN116703744B (zh) | 一种基于卷积神经网络的遥感影像匀光匀色方法和装置 | |
CN105354803B (zh) | 截断直方图均衡实现装置及方法 | |
CN104331711A (zh) | 基于多尺度模糊测度与半监督学习的sar图像识别方法 | |
CN114283431B (zh) | 一种基于可微分二值化的文本检测方法 | |
CN105872317A (zh) | 一种基于时域的视频图像滤波方法及装置 | |
CN112862718B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150429 Termination date: 20170913 |
|
CF01 | Termination of patent right due to non-payment of annual fee |