Предположим, у нас есть следующие примерные данные и модель kmeans:
from pyspark.ml.linalg import Vectors
from pyspark.ml.clustering import KMeans
import pyspark.sql.functions as F
data = [(Vectors.dense([0.0, 0.0]),), (Vectors.dense([1.0, 1.0]),),
(Vectors.dense([9.0, 8.0]),), (Vectors.dense([8.0, 9.0]),),
(Vectors.dense([10.0, 1.5]),), (Vectors.dense([11, 0.0]),) ]
df = spark.createDataFrame(data, ["features"])
n_centres = 2
kmeans = KMeans().setK(n_centres).setSeed(1)
kmModel = kmeans.fit(df)
df_pred = kmModel.transform(df)
df_pred.show()
+----------+----------+
| features|prediction|
+----------+----------+
| [0.0,0.0]| 1|
| [1.0,1.0]| 1|
| [9.0,8.0]| 0|
| [8.0,9.0]| 0|
|[10.0,1.5]| 0|
|[11.0,0.0]| 0|
+----------+----------+
Теперь давайте добавим столбец, содержащий координаты центров:
l_clusters = kmModel.clusterCenters()
# Let's convert the list of centers to a dict, each center is a list of float
d_clusters = {int(i):[float(l_clusters[i][j]) for j in range(len(l_clusters[i]))]
for i in range(len(l_clusters))}
# Let's create a dataframe containing the centers and their coordinates
df_centers = spark.sparkContext.parallelize([(k,)+(v,) for k,v in
d_clusters.items()]).toDF(['prediction','center'])
df_pred = df_pred.withColumn('prediction',F.col('prediction').cast(IntegerType()))
df_pred = df_pred.join(df_centers,on='prediction',how='left')
df_pred.show()
+----------+----------+------------+
|prediction| features| center|
+----------+----------+------------+
| 0| [8.0,9.0]|[9.5, 4.625]|
| 0|[10.0,1.5]|[9.5, 4.625]|
| 0| [9.0,8.0]|[9.5, 4.625]|
| 0|[11.0,0.0]|[9.5, 4.625]|
| 1| [1.0,1.0]| [0.5, 0.5]|
| 1| [0.0,0.0]| [0.5, 0.5]|
+----------+----------+------------+
Наконец, мы можем использовать udf для вычисления расстояния между элементами столбца и координатами центра:
get_dist = F.udf(lambda features, center :
float(features.squared_distance(center)),FloatType())
df_pred = df_pred.withColumn('dist',get_dist(F.col('features'),F.col('center')))
df_pred.show()
+----------+----------+------------+---------+
|prediction| features| center| dist|
+----------+----------+------------+---------+
| 0|[11.0,0.0]|[9.5, 4.625]|23.640625|
| 0| [9.0,8.0]|[9.5, 4.625]|11.640625|
| 0| [8.0,9.0]|[9.5, 4.625]|21.390625|
| 0|[10.0,1.5]|[9.5, 4.625]|10.015625|
| 1| [1.0,1.0]| [0.5, 0.5]| 0.5|
| 1| [0.0,0.0]| [0.5, 0.5]| 0.5|
+----------+----------+------------+---------+