In view of the limitation of inversion algorithm and the influence of thick cloud cover, MODIS aerosol products miss a lot of effective data. An ordinary Kriging-Natural proximity interpolation method based on reflectance statistical model is proposed. The response of MODIS blue band apparent reflectance to surface reflectivity and aerosol optical depth(AOD) is simulated by using 6S radiation transfer model. A statistical model of reflectivity and AOD is established. The MOD04_L2 and reflectivity data of regional aerosol products in China in November 2017 were used to carry out the experiment, and the time-space matching AERONET data was used for cross-validation. The results show that the MOD04_L2 missing data of China regions is better with this method, more than 60% of the interpolation results are within the expected error limits. The method not only effectively solves the problem of missing effective data of AOD in low reflectivity regions, but also is not affected by the improper assumption of aerosol type.