By applying the proposed algorithms and associated processing code,MODIS radiance of channel 1 to 7 are corrected for atmospheric effects to generate the surface reflectance product. A simple statistical algorithm for the estimation of MODIS spectral reflectance is developed based on the statistics of MODIS surface reflectance and the radiance at the top of atmosphere. The algorithm is validated by comparing with ground-based measurements. The suitability and limitation of the proposed method are discussed.