Knot selection for regression splines is crucial to the approximation power of splines. For univariate splines, many knot selection algorithms are available. However, extending those univariate algorithms to the multivariate case can be challenging. In this work, an estimation procedure with knot selection for multivariate regression splines is proposed, where the knot selection part is adapted from an existing univariate knot selection algorithm. Simulation results are included to demonstrate the performance of the proposed method.