code: def classify_clusters(points, n=10):
      
      arr = [[p.x, p.y] for p in points.values]
      clf = KMeans(n_clusters=n)
      clf.fit(arr)
      classes = clf.predict(arr)
      return classes
docs: Return an array of K-Means cluster classes for an array of `shapely.geometry.Point` objects.