In this paper, we study the k-means clustering scheme based on the observations of a phenomenon modelled by a sequence of random fields X 1 , · · · , X n taking values in a Hilbert space. In the k-means algorithm, clustering is performed by computing a Voronoi partition associated with centers that minimize an empirical criterion, called distorsion. The performance of the method is evaluated by comparing a theoretical distorsion of empirically optimal centers to the theoretical optimal distorsion. Our first result states that, provided the underlying distribution satisfies an exponential moment condition, an upper bound for the above performance criterion is O(1/ √ n). Then, motivated by a broad range of applications and computational matters, we use a Hölder property shared by classical random fields in stochastic modelling to construct a numerically simple algorithm that computes empirical centers based on a discretized version of the data. With a judicious choice of the discretization, we are abble to recover the same performance than in the non-discretized case.