In this paper, we deal with two different problems. First, we provide the convergence rates of multiresolution approximations, with respect to the supremum norm, for the class of elliptic splines defined in Ref. 10, and in particular for polyharmonic splines. Secondly, we consider the problem of recovering a function from a sample of noisy data. To this end, we define a linear and smooth estimator obtained from a multiresolution process based on polyharmonic splines. We discuss its asymptotic properties and we prove that it converges to the unknown function almost surely.