The solution of linear inverse problems in satellite geodesy by means of spherical spline approximation