Abstract: ‍ ‍Atmospheric disturbance is one of the main error sources in ground-based interferometric synthetic aperture radar （GB-InSAR） measurements. Under complex atmospheric conditions， conventional atmospheric phase compensation methods are no longer applicable. This paper proposes an improved method based on PS （Permanent Scatterer） technology. Firstly， the method establishes adjacent PS pairs through the Delaunay triangulation network， and analyzes the standard deviation of the differential phase sequence of each adjacent PS pair. Those noisy points or partial deformation points in the PS set can be filtered out. Then， the reference atmospheric phase sequence is estimated based on the meteorological data. By calculating the MIC （Maximal Information Coefficient） of the interference phase sequence of the remaining PS and the reference atmospheric phase sequence， together with the spatial coherence， dual thresholds are used to achieve the stable PS selection. Finally， the Kriging interpolation fitting method is used to achieve atmospheric phase estimation under complex atmospheric conditions. The method proposed in this paper is used to process the measured data in the Jiudaoguai area. Compared with the conventional methods， the atmospheric phase compensation error is effectively reduced， and the deformation measurement accuracy of GB-InSAR can be better guaranteed.