РЕДАКТИРОВАТЬ: известно, что не является полной реализацией ОПТИКИ.
Я сделал быстрый поиск и обнаружил следующее ( Оптика ). Я не могу ручаться за его качество, однако алгоритм кажется довольно простым, поэтому вы сможете быстро его проверить / адаптировать.
Вот краткий пример того, как построить кластеры на выходе алгоритма оптики:
def cluster(order, distance, points, threshold):
''' Given the output of the options algorithm,
compute the clusters:
@param order The order of the points
@param distance The relative distances of the points
@param points The actual points
@param threshold The threshold value to cluster on
@returns A list of cluster groups
'''
clusters = [[]]
points = sorted(zip(order, distance, points))
splits = ((v > threshold, p) for i,v,p in points)
for iscluster, point in splits:
if iscluster: clusters[-1].append(point)
elif len(clusters[-1]) > 0: clusters.append([])
return clusters
rd, cd, order = optics(points, 4)
print cluster(order, rd, points, 38.0)