Land surface temperature (LST) is a vital parameter associated with the land\u2013atmosphere interface. The Moderate Resolution Imaging Spectroradiometer (MODIS) LST product can provide precise LST with high time resolution, and is widely applied in various remote sensing temperature research. However, due to its inability to penetrate the cloud and fog, its quality is not able to meet the requirements of actual research. Hence, obtaining continuous and cloudless MODIS LST datasets remains challenging for researchers. The critical point is to reconstruct missing pixels. To compare the performance of different methods, first, three kinds of methods were used to reconstruct the missing pixels, namely, temporal, spatial, and spatiotemporal methods. The predicted values using these methods were validated by the automatic weather system data (AWS) in the Heihe river basin of China. The results demonstrated that, compared with other methods, linear temporal interpolation using Aqua data had the best performance in MODIS LST reconstruction in the Heihe river basin, with an RMSE of 7.13 K and an R2 of 0.82, and the NSE and PBias were 0.78 and \u22120.76%, respectively. Furthermore, the interpolation method was improved using adaptive windows and robust regression. First, the international Geosphere\u2013Biosphere Program (IGBP) classification was employed to distinguish the different land surface types. Then, the invalid LST values were reconstructed using adjacent days\u2019 effective LST values combined with a robust regression. Finally, a mean filter was applied to eliminate outliers. The overall results combined with ERA5 data were validated by AWS, with an RMSE of 6.96 K and an R2 of 0.79 and the NSE and PBias were 0.77 and \u22120.20%, respectively. The validation demonstrated that the scheme proposed in this paper is able to accurately reconstruct the missing values and improve the accuracy of the interpolation method to a certain extent when reconstructing MODIS LST.