With global climate change and increased human activities, landslides increasingly threaten human safety and property. Precisely extracting large-scale spatiotemporal information on landslides is crucial for risk management. However, existing methods are either locally based or have coarse temporal resolution, which is insufficient for regional analysis. In this study, spatiotemporal information on landslides was extracted using multiple remote sensing data from Emilia, Italy. An automated algorithm for extracting spatial information of landslides was developed with NDVI datasets. Then, we established a landslide prediction model based on a hydrometeorological threshold of three-day soil moisture and three-day accumulated rainfall. Based on this model, the locations and dates of rainfall-induced landslides were identified. Then, we further matched these identified locations with the extracted landslides from remote sensing data and finally determined the occurrence time. This approach was validated with recorded landslides events in Emilia. Despite some temporal clustering, the overall trend matched historical records, accurately reflecting the dynamic impacts of rainfall and soil moisture on landslides. The temporal bias for 87.3% of identified landslides was within seven days. Furthermore, higher rainfall magnitude was associated with better temporal accuracy, validating the effectiveness of the model and the reliability of rainfall as a landslide predictor.