Our model can describe possible nonlinear relationships besides spatial dependence. The two-step estimation procedures is proposed to estimate the unknowns. We establish the asymptotic normality of the estimators. We establish the simultaneous confidence bands of the nonparametric function. Monte Carlo simulations support the proposed method.