In order to achieve more accurate spatial reconstruction of ionospheric total electron content (TEC) and promote improved satellite positioning and ranging applications, a high accuracy spatial reconstruction (HASR) method for TEC is proposed based on the surface theory. The core theory of this method is as follows: (a) Any surface can be uniquely determined by its first and second fundamental quantities; (b) By direct difference approximation, differential equations are transformed into algebraic equations to solve Gauss equations faster. At the same time, taking parts of Europe as an example, the proposed HASR method is used to determine the correlation coefficients and the number of iterations of the model by using the relative root mean square error (RRMSE) as the evaluation criterion. The statistical results show that the TEC predicted by the HASR method is highly consistent with the actual observed values of ionospheric observation stations, and the prediction RRMSE is 9.75%. Compared with the Kriging interpolation with scale factor, the prediction accuracy of the HASR method is improved by 8.5%. We hope this method can provide ideas for the spatial reconstruction of other ionospheric parameters and further promote the realization of complete and accurate space weather forecast.