from sklearn.cluster import DBSCAN
from sklearn import metrics
from collections import namedtuple
import numpy as np
import matplotlib.pyplot as plt
from itertools import groupby

import csv_parser
import recog
import metric

if __name__ == '__main__':
    data = csv_parser.parse_data_from_csv('test0714.csv')
    groups = groupby(data, key = lambda x: x.hwid)
    entries = []
    for k, grp in groups:
        data1 = list(grp)
        entries1 = recog.recognize_entries(data)
        for i in entries1:
            entries.append(i)
    for ent in entries:
        print(ent)
    x = []
    for e in entries:
        x.append([e.lon, e.lat])
    x = np.array(x)
    db = DBSCAN(eps = 10/6400000, min_samples = 10,
            metric = metric.spherical_distance).fit(x)
    labels = db.labels_
    core_samples_mask = np.zeros_like(db.labels_, dtype = bool)
    core_samples_mask[db.core_sample_indices_] = True
    n_clusters_ = len(set(labels)) - (1 if -1 in labels else 0)
    n_noise_ = list(labels).count(-1)

    print('Estimated number of clusters: %d' % n_clusters_)
    print('Estimated number of noise points: %d' % n_noise_)
    if n_clusters_ == 0:
        print('can not get any clusters')
        plt.plot(x[:,0], x[:,1], 'o')
        plt.show()
        exit(0)
    print("Silhouette Coefficient: %0.3f"
              % metrics.silhouette_score(x, labels))

    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)

        if k != -1:
            xy = x[class_member_mask]
            print("cluster %d:" % k)
            center = np.array([np.average(xy[:, 0]), np.average(xy[:, 1])])
            print("center = %f, %f" % (center[0], center[1]))
            center_p = namedtuple("point", "lon lat")(center[0], center[1])
            rad = 0
            for loc in xy:
                p = namedtuple("point", "lon lat")(loc[0], loc[1])
                rad = max(rad, metric.spherical_distance(center_p, p))
            rad *= 6400 * 1000
            print("radius = %f m" % rad)
    plt.show()