The invention discloses a method for obtaining a deformation sequence of a high temporal-spatial resolution by fusing GNSS and InSAR data. The method comprises the following steps: according to the InSAR data, selecting an appropriate trend field; fitting the detrended residual volume by utilizing a variation function model, and thus obtaining a spatial variation function; according to PGS point locations, the variation function and trend terms, constructing a space field matrix in Kriging Kalman filtering; estimating model parameters in a Kriging Kalman filtering model by using EM; carrying out Kriging Kalman filtering; carrying out interpolation in time and space, and thus obtaining the deformation sequence of the high temporal-spatial resolution. Compared with the current method, the method adopts the InSAR data with the high temporal-spatial resolution instead of adopting the GPS data with sparse point locations during the selection of the variation function and the trend terms. Thus, the method can effectively improve the precision of the space field and obtain the deformation sequence of the high temporal-spatial resolution.