Near-surface PM2.5 estimates remain a global scientific research challenge due to their effect on human fitness and atmospheric environmental quality. However, practical near-surface PM2.5 estimates are impeded by the incomplete monitoring data. In this study, we propose the trigonometric polynomial fitting (TPF) method to estimate near-surface PM2.5 concentrations in south-central China during 2015. We employ 10-fold cross-validation (CV) to assess the reliability of TPF in estimating practical PM2.5 values. When compared to alternative methods such as the orthogonal polynomial fitting (OBF) method based on Chebyshev basis functions, Kriging interpolation, and radial basis function (RBF) interpolation, our results show that utilizing TPF31, with a maximum order of 3 in the x direction and a maximum order of 1 in the y direction, leads to superior efficiency through error minimization. TPF31 reduces MAE and RMSE by 1.93%, 24%, 6.96% and 3.6%, 23.07%, 10.43%, respectively, compared to the other three methods. In addition, the TPF31 method effectively reconstructs the spatial distribution of PM2.5 concentrations in the unevenly distributed observation stations of Inner Mongolia and the marginal regions of the study area. The reconstructed spatial distribution is remarkably smooth. Despite the non-uniform distribution of observation stations and the presence of missing data, the TPF31 method demonstrates exceptional effectiveness in accurately capturing the inherent physical attributes of spatial distribution. The theoretical and experimental results emphasize that the TPF method holds significant potential for accurately reconstructing the spatial distribution of PM2.5 in China.