Python KMeans聚类问题分析

作者:落叶_小唱 时间:2023-11-24 14:15:52 

今天用python实现了一下简单的聚类分析,顺便熟悉了numpy数组操作和绘图的一些技巧,在这里做个记录。


from pylab import *
from sklearn.cluster import KMeans

## 利用numpy.append()函数实现matlab多维数组合并的效果,axis 参数值为 0 时是 y 轴方向合并,参数值为 1 时是 x 轴方向合并,分别对应matlab [A ; B] 和 [A , B]的效果

#创建5个随机的数据集
x1=append(randn(500,1)+5,randn(500,1)+5,axis=1)
x2=append(randn(500,1)+5,randn(500,1)-5,axis=1)
x3=append(randn(500,1)-5,randn(500,1)+5,axis=1)
x4=append(randn(500,1)-5,randn(500,1)-5,axis=1)
x5=append(randn(500,1),randn(500,1),axis=1)

# 下面用较笨的方法把5个数据集合并成 (2500,2)大小的数组data
data=append(x1,x2,axis=0)
data=append(data,x3,axis=0)
data=append(data,x4,axis=0)
data=append(data,x5,axis=0)

plot(x1[:,0],x1[:,1],'oc',markersize=0.8)
plot(x2[:,0],x2[:,1],'og',markersize=0.8)
plot(x3[:,0],x3[:,1],'ob',markersize=0.8)
plot(x4[:,0],x4[:,1],'om',markersize=0.8)
plot(x5[:,0],x5[:,1],'oy',markersize=0.8)

k=KMeans(n_clusters=5,random_state=0).fit(data)
t=k.cluster_centers_ # 获取数据中心点

plot(t[:,0],t[:,1],'r*',markersize=16) # 显示这5个中心点,五角星标记~

title('KMeans Clustering')
box(False)

xticks([])  # 去掉坐标轴的标记
yticks([])

show()

结果如下:

Python KMeans聚类问题分析

2017/01/11更新

今天重新试运行程序的出现报错了,提示导入NUMPY_MKL失败,因为之前用命令pip install -U numpy手动更新了numpy,最初的是在http://www.lfd.uci.edu/~gohlke/pythonlibs/#numpy 里下载的numpy-1.11.2+mkl-cp27-cp27m-win_amd64.whl 文件安装的,只要重新安装回去就可以了

2017/1/18更新

python中还有一个叫plotly 的package,可以通过pip install plotly 或 pip3 install plotly(Python3.X) ,使用这个package可以绘制精美的图像,官网中有很多例子介绍,同时plotly 还支持matlab,R等,但是个人觉得plotly 的绘图语法相比matplotlib 的繁琐,需要照着例程来修改才比较方便,不过如果只是要想数据可视化更好看的话参考官网例程并做修改也无妨,下面是来自官网的一段示例代码:


import plotly.plotly as py
import plotly.graph_objs as go
import plotly
import numpy as np

#生成三组高斯分布(Gaussian Distribution)点集

x0 = np.random.normal(2, 0.45, 300)
y0 = np.random.normal(2, 0.45, 300)

x1 = np.random.normal(6, 0.8, 200)
y1 = np.random.normal(6, 0.8, 200)

x2 = np.random.normal(4, 0.3, 200)
y2 = np.random.normal(4, 0.3, 200)

#创建图形对象 graph object
trace0 = go.Scatter(
x=x0,
y=y0,
mode='markers',
)
trace1 = go.Scatter(
x=x1,
y=y1,
mode='markers'
)
trace2 = go.Scatter(
x=x2,
y=y2,
mode='markers'
)
trace3 = go.Scatter(
x=x1,
y=y0,
mode='markers'
)
#布局是一个字典,字典关键字keys包括:'shapes', 'showlegend'
layout = {
'shapes': [
 {
  'type': 'circle',
  'xref': 'x',
  'yref': 'y',
  'x0': min(x0),
  'y0': min(y0),
  'x1': max(x0),
  'y1': max(y0),
  'opacity': 0.2,
  'fillcolor': 'blue',
  'line': {
   'color': 'blue',
  },
 },
 {
  'type': 'circle',
  'xref': 'x',
  'yref': 'y',
  'x0': min(x1),
  'y0': min(y1),
  'x1': max(x1),
  'y1': max(y1),
  'opacity': 0.2,
  'fillcolor': 'orange',
  'line': {
   'color': 'orange',
  },
 },
 {
  'type': 'circle',
  'xref': 'x',
  'yref': 'y',
  'x0': min(x2),
  'y0': min(y2),
  'x1': max(x2),
  'y1': max(y2),
  'opacity': 0.2,
  'fillcolor': 'green',
  'line': {
   'color': 'green',
  },
 },
 {
  'type': 'circle',
  'xref': 'x',
  'yref': 'y',
  'x0': min(x1),
  'y0': min(y0),
  'x1': max(x1),
  'y1': max(y0),
  'opacity': 0.2,
  'fillcolor': 'red',
  'line': {
   'color': 'red',
  },
 },
],
'showlegend': False,
}
data = [trace0, trace1, trace2, trace3]
#图像包括数据部分和布局部分
fig = {
'data': data,
'layout': layout,
}
#使用离线的方式绘制图像,因为没有注册官方的网站,而且那个网站不容易进去,所以用离线绘制
plotly.offline.plot(fig, filename='clusters')

结果是通过浏览器打开图片的,可以保存到本地,如下图:

Python KMeans聚类问题分析

总结:plotly 这个库虽然语法比较繁琐,但是对数据显示要求较高的情况下可以充分利用,一般绘图的话使用matplotlib比较方便,特别是ipython模式下先执行from pylab import * 可以获得和MATLAB 类似的工作环境。

来源:http://blog.csdn.net/ouening/article/details/53940289

标签:Python,KMeans,聚类
0
投稿

猜你喜欢

  • pandas Dataframe行列读取的实例

    2021-06-11 01:00:45
  • 在MySQL中删除表的操作教程

    2024-01-15 09:37:46
  • Access的特点及其概念问答

    2009-09-10 19:00:00
  • 全面了解python中的类,对象,方法,属性

    2021-10-07 10:54:50
  • 19个ASP编程基础典型代码

    2008-10-23 15:46:00
  • Python参数解析模块sys、getopt、argparse使用与对比分析

    2021-12-19 17:39:58
  • Linux oracle数据库自动备份自动压缩脚本代码

    2024-01-19 01:44:40
  • XHTML 和 DOCTYPE 切换

    2007-05-31 09:30:00
  • 解决在keras中使用model.save()函数保存模型失败的问题

    2021-09-12 17:32:41
  • 基于python实现FTP文件上传与下载操作(ftp&sftp协议)

    2021-11-03 00:11:26
  • 可以在线创建文件夹吗?

    2009-11-01 18:07:00
  • 利用Go语言实现Raft日志同步

    2024-03-15 02:55:20
  • Python处理CSV与List的转换方法

    2021-11-22 11:53:48
  • Pytorch教程内置模型源码实现

    2022-09-04 12:58:50
  • 如何用Idea或者webstorm跑一个Vue项目(步骤详解)

    2022-08-20 17:06:23
  • 通过T_sql语句向其中一次填入一条数据或一次填入多条数据的方式填充数据

    2012-11-30 19:55:34
  • Python基于动态规划算法计算单词距离

    2021-09-18 11:18:16
  • Python数据可视化Pyecharts库实现桑葚图效果

    2022-01-05 23:39:44
  • 深入理解Vue 的条件渲染和列表渲染

    2024-04-09 10:46:37
  • SQL Server2019数据库备份与还原脚本(批量备份)

    2024-01-27 06:46:47
  • asp之家 网络编程 m.aspxhome.com