The proposed method has three steps. Firstly, smooth the MODIS-NDVI time-series data with the Savitzky-Golay (S-G) filter of TIMESAT (Jönsson and Eklundh, 2004). Derive phenological indices from the fitted MODIS-NDVI time-series data using a dynamic-threshold method via MATLAB. Resample time series MODIS-NDVI and phonological indices to 30 m resolution. Then, determine the critical stages of rice growth based on the MODIS-NDVI time-series data, phenological indices, and the observed