Python scikit-DBSCAN：错误的坐标或聚类

Question

我正在写 DBSCAN，但遇到了一些奇怪的问题。（2 个问题）

这是我的代码：

第一部分有问题，如果我添加X = StandardScaler().fit_transform(X) 结果坐标是错误的！但是如果我没有添加这段代码，它总是会是一个聚类（但结果的坐标是正确的！）。我尝试调整 esp 或 min_samples，但没有改变。

dataSet = []
fileIn = open('data')
for line in fileIn.readlines():
    lineArr = line.strip().split('\t')
    dataSet.append([float(lineArr[0]), float(lineArr[1])]) 
numSamples = len(dataSet)
X = np.array(dataSet) 
X = StandardScaler().fit_transform(X)

(添加) (无)

db = DBSCAN(eps=0.5, min_samples=10).fit(X)
core_samples_mask = np.zeros_like(db.labels_, dtype=bool)
core_samples_mask[db.core_sample_indices_] = True
labels = db.labels_
print(labels)

counters = {}
for item in labels:
if item in counters:
    counters[item] += 1
else:
    counters[item] = 1
print ("Count of different cluster:(#r,g,b,a)")#r,g,b,a
print (counters)   
n_clusters_ = len(set(labels)) - (1 if -1 in labels else 0)

print('Estimated number of clusters: %d' % n_clusters_)

第二个问题是我尝试绘制我计算的坐标，但我不知道为什么它显示结果如此错误！

clusters = [np.mean(X[labels == i],axis=0) for i in range(n_clusters_)]

outliers = X[labels == 0]
print(clusters)

for i in range(n_clusters_):
    plt.plot(clusters[i],'*',markersize=20)

unique_labels = set(labels)
colors = [plt.cm.Spectral(each)
          for each in np.linspace(0, 1, len(unique_labels))]
for k, col in zip(unique_labels, colors):
    if k == -1:
        # Black used for noise.
        col = [0, 0, 0, 1]

    class_member_mask = (labels == k)

    xy = X[class_member_mask & core_samples_mask]
    plt.plot(xy[:, 0], xy[:, 1], 'o',             markerfacecolor=tuple(col),markeredgecolor='k', markersize=14)
    xy = X[class_member_mask & ~core_samples_mask]
    plt.plot(xy[:, 0], xy[:, 1], 'o', markerfacecolor=tuple(col),markeredgecolor='k', markersize=6)
plt.title('Estimated number of clusters: %d' % n_clusters_)
plt.show()

请帮我谢谢你！

Answer 1

您显然需要为您的坐标系适当地选择epsilon。如果您缩放数据，epsilon 将不再相同。您认为最简单的一个技巧（但这不是“正确”的解决方案）是使用 unscaled 数据简单地计算平均值。但是DBSCAN集群的方式无论如何都不可靠。

从您的坐标轴来看，您可能需要将 epsilon 减少 100 倍。

由于您的数据显然是坐标，您应该使用半正弦距离，因为地球并不平坦，并根据对您的问题有意义的距离选择 epsilon。确切的缩放比例可能有点棘手。它可能是弧度，因此您需要将英里转换为弧度以转换距离。