A regularization method is proposed for the problem of reconstructing a univariate function from approximate average values over subintervals of its domain. The paper establishes the unique solvability of the method, the characterization of the solution as a quadratic spline, as well as error estimates. Numerical examples are provided to illustrate and compare the computational performance of the method.