Estimation of BVk functions from scattered data