obtaining a time-series SAR data set, and performing registration and L-looks processing on the time-series SAR data set to obtain an L-looks intensity data set and an interferometric data set respectively by a processor;

taking the L-looks intensity data set as a reference, obtaining a preset digital elevation model (DEM) and a preset land cover image, performing registration and geocoding on the preset DEM to obtain a digital elevation in a SAR image coordinate system, and performing registration and geocoding on the preset land cover image to obtain a land cover image in the SAR image coordinate system by the processor;

performing a differential operation on the interferometric data set based on the digital elevation in the SAR image coordinate system to obtain a differential interferometric data set by the processor; and

estimating a covariance matrix of the time-series SAR data set at each spatial pixel position based on the land cover image in the SAR image coordinate system and the differential interferometric data set, estimating and obtaining an optimized time-series phase according to a preset maximum likelihood estimator for a time-series phase, and obtaining a ground deformation monitoring result based on the optimized time-series phase by the processor.