It seems you are working with a 1000-dimensional feature space, which would be hard to plot.
You could start by plotting only two dimensions and use plt.scatter and check the results first.
To reduce your feature space you could use e.g. PCA, but I’m not sure how this would affect the kmeans centroids and data. Might still be worth a try.