Land surface temperature is one of the important parameters of scientific research such as resource environment, climate change and terrestrial ecosystem. MODIS LST(Land Surface Temperature, LST) products are important data sources for land surface temperature related research. The land surface temperature information of MODIS LST products is lost in the cloud coverage area. Therefore, the land surface temperature estimation of cloud coverage areas has become a frontier research problem of thermal infrared remote sensing. In order to solve the problem of missing land surface temperature information in the cloud occlusion area of MODIS LST products. In this paper, the Qinling area is used as the research area and the experimental data of MOD11A2 from 2001 to 2017 is selected. In the traditional Inverse Distance Weighting(IDW), Regular Spline(SPLINE), Ordinary Kriging(OK) and Trend Surface(TREND) spatial interpolation method, the important influence factor of elevation is introduced. Through a large number of spatial interpolation experiments, the traditional spatial interpolation method is improved, and a MODIS LST spatial interpolation method based on DEM correction is formed. Analysis of spatial interpolation results indicates:(1) The spatial interpolation accuracy is from high to low: OK> SPLINE > IDW>TREND, and the accuracy of the OK, SPLINE, IDW, and TREND methods based on DEM correction is increased by about 0.38°C, 0.31°C, 0.32°C, and 0.78°C, respectively;(2) The accuracy of spatial interpolation results shows seasonal differences. The interpolation accuracy is higher in summer, July, and August, and the interpolation accuracy is the lowest in January.(3)The interpolation accuracy has a certain relationship with the cloud area. When the cloud coverage area is less than 1.1km~2, the interpolation error of the DEM+OK interpolation method is less than 0.55°C, and when the cloud coverage area is less than 3.1km~2, the spatial interpolation error is less than 1°C. When the cloud coverage area is less than 2.7 km~2, the interpolation error of the DEM+SPLINE method is less than 0.55°C, and the interpolation error of the DEM+SPLINE method is less than 1°C when the cloud coverage area is less than 10.4 km~2. When the cloud coverage is 1.1~2.7 km~2, the interpolation accuracy of DEM+SPLINE interpolation method is higher than of the DEM+OK interpolation method.