A two stage $k$-monotone B-spline regression estimator: Uniform Lipschitz property and optimal convergence rate