CN103744962B - 一种实现K‑means聚类的方法及装置 - Google Patents
一种实现K‑means聚类的方法及装置 Download PDFInfo
- Publication number
- CN103744962B CN103744962B CN201410005652.4A CN201410005652A CN103744962B CN 103744962 B CN103744962 B CN 103744962B CN 201410005652 A CN201410005652 A CN 201410005652A CN 103744962 B CN103744962 B CN 103744962B
- Authority
- CN
- China
- Prior art keywords
- clustering
- clusters
- sample
- cluster
- sample set
- 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.)
- Active
Links
- 238000003064 k means clustering Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000007418 data mining Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本申请公开了一种实现K‑means聚类的方法及装置,包括:对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;利用获得的样本集形成K个集群的层次聚类,将层次聚类的K个集群的质心作为K‑means聚类的初始质心;根据获取的K‑means初始质心进行K‑means聚类。本发明通过选取需要数量的聚类集样本进行层次聚类,以获取层次聚类的质心作为K‑means聚类方法的初始质心,使K‑means聚类方法不受初始质心影响。另一方面,通过聚类算法过程中不满足聚类第一数值的聚类进行删除,从最大的聚类开始拆分相应个数的聚类,以保证聚类要求的聚类集群的个数要求。
Description
技术领域
本发明涉及数据挖掘技术,尤指一种实现K-means聚类的方法及装置。
背景技术
K-means聚类的方法作为数据挖掘领域的一种常见方法,在一些数据处理的过程中存在以下问题:
一方面,初始质心的选择对聚类结果影响很大。也就是指作为聚类方法的第一步,如果对质心的选取不是数据集中的正常数值,那么对于聚类方法获得的结果,由于初值选择的随机性,很可能得出的聚类效果是不理想的。图1为现有技术对同一聚类数据采用随机选取的初始质心进行聚类的四次聚类结果示意图,如图1所示,聚类结果a1、聚类结果a2和聚类结果a3是采用第一组随机选取初始质心的聚类结果;聚类结果b1、聚类结果b2和聚类结果b3是采用第二组随机选取初始质心的聚类结果;聚类结果c1、聚类结果c2和聚类结果c3是采用第三组随机选取初始质心的聚类结果;聚类结果d1、聚类结果d2和聚类结果d3是采用第四组随机选取初始质心的聚类结果。图2是现有技术采用随机选取初始质心的聚类结果示意图,如图2所示,聚类结果k1和聚类结果k2为聚类形成的两个结果,但是从图中可以看出,这两个聚类集群的理想效果应当属于一个聚类集群;而聚类结果k3和聚类结果K4从聚类的结果示意图上可以看出,其从聚类要达到的理想效果应当是聚类为3个集群。通过分析可见,在现有的K-means聚类方法中,如果一个数据集有K个“真实”的集群,那么随机选择的K个点作为聚类的初始质心,则每个初始质心恰好落在一个“真实”的集群上的概率很小。例如,一个数据集上有10个类,按照目前的聚类方法随机选择10个点作为初始质心。那么10个初始点正好落在10个类的集群上的概率为10!/1010=0.00036。
另一方面,K-means聚类方法对异常点的敏感性高。当数据集含存在一定数量的异常点时,则异常点对SSE的贡献必然很大(采用的是距离的平方)。由于异常点的存在导致在算法迭代的过程中,质心非常容易被异常点“吸附”过去,从而严重影响了聚类方法的分析效果。
再者,出现空集群问题。空集群是指在一次聚类方法实施的过程中,对于选取的一个或几个质心,任何一个数据点都没有被分配到该质心上,目前的K-means聚类方法无法避免聚类时出现空集群的现象。
综上所述,目前的K-means聚类方法,采用随机选取初始质心的选取进行聚类,在聚类过程中,异常点对聚类结果影响较大,且无法避免出现空集群的现象,严重影响了K-means聚类方法在数据处理过程中的使用。
发明内容
为了解决上述技术问题,本发明公开了一种实现K-means聚类的方法及装置。能够选取出适合K-means聚类的初始质心,降低K-means聚类初始质心对聚类结果的影响。
为了达到本申请的目的,本发明提供一种实现K-means聚类的方法,包括:
对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;
利用获得的样本集形成K个集群的层次聚类,将层次聚类的K个集群的质心作为K-means聚类的初始质心;
根据获取的K-means初始质心进行K-means聚类。
进一步地,满足预先设置的相似度数值为:预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
进一步地,在进行K-means聚类时,设置聚类集群的聚类第一数值,该方法还包括:
删除不满足聚类第一数值的聚类集群,在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类的集群个数不变。
另一方面,本申请还提供一种实现K-means聚类的装置,包括:样本单元、初始质心获取单元及K-means聚类单元;其中,
样本单元,用于对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;
初始质心获取单元,用于利用样本单元获得的样本集进行形成K个集群的层次聚类,将获取的层次聚类的K个集群的质心设置为K-means聚类的初始质心;
K-means聚类单元,用于根据初始质心获取单元设置的K-means聚类的初始质心进行K-means聚类。
进一步地,样本单元具体用于:
预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
进一步地,该装置还包括聚类处理单元,用于在进行K-means聚类时,设置聚类集群的第一数值,删除不满足第一数值的聚类集群;在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类集群个数不变。
本申请提出一种技术方案,包括:对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;利用获得的样本集形成K个集群的层次聚类,将层次聚类的K个集群的质心作为K-means聚类的初始质心;根据获取的K-means初始质心进行K-means聚类。本申请还包括与方法对应的装置。本发明通过选取需要数量的聚类集样本进行层次聚类,以获取层次聚类的质心作为K-means聚类方法的初始质心,使K-means聚类方法不受初始质心影响。另一方面,通过聚类算法过程中不满足聚类第一数值的聚类进行删除,从最大的聚类开始拆分相应个数的聚类,以保证聚类要求的聚类集群的个数要求。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为现有技术对同一聚类数据采用随机选取的初始质心进行聚类的四次聚类结果示意图;
图2是现有技术采用随机选取初始质心的聚类结果示意图;
图3为本发明实现K-means聚类的方法的流程图;
图4为本发明实现K-means聚类的聚类效果示意图;
图5为本发明实现K-means聚类的装置的结构框图。
具体实施方式
图3为本发明实现K-means聚类的方法的流程图,如图3所示,包括:
步骤300、对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值。
本步骤中,满足预先设置的相似度数值为:预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
需要说明的是,采用上述公式是为了获得进行层次聚类的样本集中的样本数值,样本数值的取值大小根据实际的需求进行相应的设定。假设以相似度为95%和聚类个数为20进行样本集的样本数值计算,可以获得样本数值为117。当设置的概率越大,样本数值也会呈现增大的趋势,同时,由于样本数值的增大,获得的聚类效果应当会相应的变为更为理想;需要强调的是采用其他概率的计算方法也可以得到本发明样本集的样本数值,当获得的样本集的样本数值与本发明计算样本数值相同时,采用本发明层次聚类对获取初始质心,其技术效果是相同的,应当属于本发明方法保护的范围。
另外,从本发明的计算公式可以得出,对于样本集的大小和聚类数据的数据个数无关。也就是说,如果在大数据或其他数据较多的数据挖掘领域使用本发明方法,本发明方法是具有恒优的聚类效果的,不会因为数据集群数据量大小的改变而影响本发明方法的技术效果。
步骤301、利用获得的样本集形成K个集群的层次聚类,将层次聚类的K个集群的质心作为K-means聚类的初始质心。
步骤302、根据获取的K-means初始质心进行K-means聚类。
在进行K-means聚类时,设置聚类集群的聚类第一数值,本发明方法还包括:
删除不满足聚类第一数值的聚类集群,在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类集群个数不变。
需要说明的是,聚类的第一数值,可以是一个动态的数值,针对集群数据及集群个数,和对聚类的效果要求,设置在聚类过程的一个随着聚类的进展,不断增加变化的取值;或者,还可以是一个固定的取值,在对聚类数据进行聚类到一定程度时,设置一个固定的取值,用来判断聚类是否满足一定的聚类效果,如果无法满足,则根据本发明方法,进行删除处理。
图4为本发明实现K-means聚类的聚类效果示意图,如图4所示,通过计算机模拟仿真方式,对本发明实现K-means聚类效果进行比较分析,效果示意图为维度为2维,聚类个数为15,数据集大小为5000的一次K-means聚类的聚类结果。
通过对K-means聚类的初始质心采用本发明方法进行获取,解决了初始质心的选择对结果影响很大,从聚类结果K1、聚类结果K2、聚类结果K3、聚类结果K4与图2的聚类结果k1、聚类结果k2、聚类结果k3、聚类结果k4进行比较可以很清楚的看出本发明的技术效果。
另外,在聚类过程中,采用对在进行K-means聚类时,设置聚类集群的聚类第一数值,删除不满足聚类第一数值的聚类集群,在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类集群个数不变。解决了K-means聚类方法对异常点的敏感性问题。同时,通过上述方法的处理,聚类也避免了空集群的问题,使聚类效果更加理想。
图5为本发明实现K-means聚类的装置的结构框图,如图5所示,包括:样本单元、初始质心获取单元及K-means聚类单元;其中,
样本单元,用于对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值。
样本单元具体用于:预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
初始质心获取单元,用于利用样本单元获得的样本集进行形成K个集群的层次聚类,将获取的层次聚类的K个集群的质心设置为K-means聚类的初始质心。
K-means聚类单元,用于根据初始质心获取单元设置的K-means聚类的初始质心进行K-means聚类。
本发明装置还包括聚类处理单元,用于在进行K-means聚类时,设置聚类集群的第一数值,删除不满足第一数值的聚类集群;在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类集群个数不变。
需要说明的是,聚类处理单元可以是与K-means聚类单元连接,进行交互通信的一个单元,也可以是内嵌在K-means聚类单元内部的一个单元。另外,聚类单元可以是软件实现的功能单元,也可以是由相关芯片形成的功能单元。
虽然本申请所揭露的实施方式如上,但所述的内容仅为便于理解本申请而采用的实施方式,并非用以限定本申请。任何本申请所属领域内的技术人员,在不脱离本申请所揭露的精神和范围的前提下,可以在实施的形式及细节上进行任何的修改与变化,但本申请的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (4)
1.一种实现K-means聚类的方法,其特征在于,包括:
对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;
利用获得的样本集形成K个集群的层次聚类,将层次聚类的K个集群的质心作为K-means聚类的初始质心;
根据获取的K-means初始质心进行K-means聚类;
所述满足预先设置的相似度数值为:预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
2.根据权利要求1所述的方法,其特征在于,在进行K-means聚类时,设置聚类集群的聚类第一数值,该方法还包括:
删除不满足聚类第一数值的聚类集群,在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类的集群个数不变。
3.一种实现K-means聚类的装置,其特征在于,包括:样本单元、初始质心获取单元及K-means聚类单元;其中,
样本单元,用于对聚类数据进行随机抽样得到样本集,样本集与聚类数据满足预先设置的相似度数值;
初始质心获取单元,用于利用样本单元获得的样本集进行形成K个集群的层次聚类,将获取的层次聚类的K个集群的质心设置为K-means聚类的初始质心;
K-means聚类单元,用于根据初始质心获取单元设置的K-means聚类的初始质心进行K-means聚类;
所述样本单元具体用于:
预先设置相似度数值,通过计算公式(1-(1-1/K)s)K的数值大于或等于预先设置的相似度数值,获得样本集的样本数值;
其中,S为需要计算的样本集的样本数值。
4.根据权利要求3所述的装置,其特征在于,该装置还包括聚类处理单元,用于在进行K-means聚类时,设置聚类集群的第一数值,删除不满足第一数值的聚类集群;在删除聚类集群时,将最大集群开始的相应个数的集群进行拆分,以满足K-means聚类集群个数不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410005652.4A CN103744962B (zh) | 2014-01-06 | 2014-01-06 | 一种实现K‑means聚类的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410005652.4A CN103744962B (zh) | 2014-01-06 | 2014-01-06 | 一种实现K‑means聚类的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103744962A CN103744962A (zh) | 2014-04-23 |
CN103744962B true CN103744962B (zh) | 2017-01-25 |
Family
ID=50501980
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410005652.4A Active CN103744962B (zh) | 2014-01-06 | 2014-01-06 | 一种实现K‑means聚类的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103744962B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104063518B (zh) * | 2014-07-14 | 2017-06-09 | 南京弘数信息科技有限公司 | 一种针对大数据的分解组合聚类方法 |
CN107430824B (zh) * | 2015-02-06 | 2021-05-04 | 意识教育以色列公司 | 用于评价响应的半自动系统和方法 |
CN105183765A (zh) * | 2015-07-30 | 2015-12-23 | 成都鼎智汇科技有限公司 | 一种基于大数据的话题抽取方法 |
CN105068991A (zh) * | 2015-07-30 | 2015-11-18 | 成都鼎智汇科技有限公司 | 一种基于大数据的舆情发现方法 |
CN105118332B (zh) * | 2015-09-30 | 2018-01-23 | 中国民用航空总局第二研究所 | 一种基于聚类分析法的空中交通管制模拟仿真异常检测方法和装置 |
CN106484838A (zh) * | 2016-09-30 | 2017-03-08 | 中国南方电网有限责任公司 | 基于数据挖掘的安全检查标准库动态更新方法 |
CN108062555A (zh) * | 2016-11-08 | 2018-05-22 | 南京理工大学 | 基于Spark流式聚类的监测数据预警系统 |
CN110880005B (zh) * | 2018-09-05 | 2023-06-23 | 阿里巴巴集团控股有限公司 | 向量索引建立方法及装置和向量检索方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254020A (zh) * | 2011-07-22 | 2011-11-23 | 西安电子科技大学 | 基于特征权重的全局k-均值聚类方法 |
-
2014
- 2014-01-06 CN CN201410005652.4A patent/CN103744962B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254020A (zh) * | 2011-07-22 | 2011-11-23 | 西安电子科技大学 | 基于特征权重的全局k-均值聚类方法 |
Non-Patent Citations (1)
Title |
---|
聚类技术在车险业务分析中的应用研究;彭崇;《中国优秀硕士学位论文全文数据库 信息科技辑》;20070915(第3期);第1、33-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103744962A (zh) | 2014-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103744962B (zh) | 一种实现K‑means聚类的方法及装置 | |
CN108090567B (zh) | 电力通信系统故障诊断方法及装置 | |
US10380117B2 (en) | Event occurrence place estimation method, computer-readable recording medium storing event occurrence place estimation program, and event occurrence place estimation apparatus | |
US9177020B2 (en) | Gathering index statistics using sampling | |
KR101850993B1 (ko) | 클러스터 기반 키워드 산출 방법 및 장치 | |
JP2019511773A (ja) | サービス・パラメータ選択方法及び関連するデバイス | |
CN108805174A (zh) | 聚类方法及装置 | |
CN107346410B (zh) | 图片处理方法及装置 | |
CN104598632A (zh) | 热点事件检测方法和装置 | |
CN106301979B (zh) | 检测异常渠道的方法和系统 | |
CN109670153A (zh) | 一种相似帖子的确定方法、装置、存储介质及终端 | |
CN113094415B (zh) | 数据抽取方法、装置、计算机可读介质及电子设备 | |
CN107085603B (zh) | 一种数据处理方法及装置 | |
CN106651408B (zh) | 一种数据分析方法及装置 | |
CN109669875B (zh) | 应用崩溃测试方法、装置、电子设备及存储介质 | |
CN107154971A (zh) | 一种应用下载方法及装置 | |
CN111080362A (zh) | 广告监测系统及方法 | |
CN104869531B (zh) | 一种群组成员更新方法和装置 | |
CN108763400B (zh) | 基于对象行为和主题偏好的对象划分方法及装置 | |
CN105653533B (zh) | 一种更新分类关联词集合的方法和装置 | |
CN116089985A (zh) | 一种分布式日志的加密存储方法、装置、设备及介质 | |
CN112887426B (zh) | 信息流的推送方法、装置、电子设备以及存储介质 | |
CN115344627A (zh) | 数据筛选方法、装置、电子设备及存储介质 | |
US20150052238A1 (en) | Device Compatibility Management | |
CN103678545A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180817 Address after: 200436 Room 411, No. three, JIANGCHANG Road, Jingan District, Shanghai, 411 Patentee after: Shanghai wave Cloud Computing Service Co., Ltd. Address before: 100085 floor 1, C 2-1, No. 2, Shang Di Road, Haidian District, Beijing. Patentee before: Electronic information industry Co.,Ltd of the tide (Beijing) |