OPTICS聚类Python代码

# -*- coding: gbk -*-
import numpy as np
import matplotlib.pyplot as plt
import operator
from scipy.spatial.distance import pdist
from scipy.spatial.distance import squareform
import sys
from sklearn import metrics
import pandas as pd


def fit_plot_Kmean_model(datalist):
    data = pd.read_csv(datalist, header=None, error_bad_lines=False)
    # 必须添加header=None，否则默认把第一行数据处理成列名导致缺失
    data_list = data.values.tolist()
    data_list = np.array(data_list)
    return data_list

def compute_squared_EDM(X):
  return squareform(pdist(X,metric='euclidean'))

# # 显示决策图
# def plotReachability(data,eps):
#     plt.figure()
#     plt.plot(range(0,len(data)), data)
#     plt.plot([0, len(data)], [eps, eps])
#     plt.show()

# 显示分类的类别
def plotFeature(eps,minPts,data,labels,savePath):
    eps=float(eps)
    minPts=int(minPts)
    plt.figure(figsize=(12, 12))
    plt.subplot(221)
    plt.scatter(data[:, 0], data[:, 1])
    plt.title(&#