Fitting and filling of 3D datasets with volume constraints using radial basis functions under tension