This study focuses on the Badui region in eastern Tibet, an area with complex topography featuring numerous valleys, ravines, and frequent geological hazards. Given the economic expansion in this region, advanced techniques are essential for analyzing the distribution of geological hazards and developing early warning of geological hazards. The research employs ESBAS-InSAR (Enhanced Small Baseline Subset Interferometric Synthetic Aperture Radar) technology, which provides more ascending and descending data than traditional SBAS-InSAR, allowing reprojection into vertical and horizontal components. Following dimensionality reduction through PCA (Principal Component Analysis) and k-means clustering, the horizontal displacements were categorized into four clusters, and the vertical displacements were categorized into five clusters. Time series data of vertical and horizontal displacements, rainfall, and NDVI (Normalized Difference Vegetation Index) were then used to assess 16 displacement prediction models. The Temporal Fusion Transformer (TFT) model demonstrated the best predictive performance. To further improve accuracy, 11 static variables such as clusters, elevation, slope, aspect, distance from faults, and time varying known categorical variable, earthquake times, were added as the TFT input variables. Results indicate that optimized TFT model reduces the RMSE (Root Mean Square Error) from 3.4842 to 2.1707, the MAPE (Mean Absolute Percentage Error) from 2625.6399 to 2154.5505, and the MAE (Mean Absolute Error) from 2.4392 to 2.3731. Overall, this study provides a framework for multivariate, multistep forecasting of diverse deformation modes across large areas and identifies distinct landslide deformation patterns through clustering, thereby enhancing the prediction of landslide deformation.