Compared with common dual-temporal satellite images, satellite time series images contain richer surface information and can alleviate the impact of foreign objects with the same spectrum and the same object with different spectra. Therefore, they play an important role in change detection. However, the change detection methods for satellite time series images are mostly based on pixels and ignore the spatial relationship between pixels and their surroundings. This causes noise in the change detection result. Accordingly, this study proposed a method of change detection based on spatial-temporal-spectral features(CDSTS) for satellite time series images. First, the temporal, spatial (textural and statistical), and spectral features of each pixel were extracted from Landsat time series images using a gray-level co-occurrence matrix and local statistical calculation methods. Then, anomalies of time series features were automatically screened according to the time series performance regularity of each pixel in different bands. These anomalies were then fused with the detection results of the continuous change detection and classification method (CCDC) to obtain high-precision changed/unchanged training sample points. Finally, the SVM classifier was trained using the training sample points and their corresponding spatial-temporal-spectral features for full graph classification. The results show that the CDSTS algorithm significantly outperforms the commonly used time series change detection algorithms CCDC and COLD (continuous monitoring of land disturbance) in terms of change detection precision, with the overall precision improved by 4.8 to 11.7 percentage points.