B-spline function-based approach for GPS tropospheric tomography