Fitting a manifold of large reach to noisy data