Fit regions and functions of bounded variation