The micro-deformation monitoring radar is usually based on Permanent Scatterer (PS) technology to realize deformation inversion. When the region is continuously monitored for a long time, the radar image amplitude and pixel variance will change significantly with time. Therefore, it is difficult to select phase-stable scatterers by conventional amplitude deviation methods, as they can seriously affect the accuracy of deformation inversion. For different regions studied within the same scenario, using a PS selection method based on the same threshold often increases the size of the deformation error. Therefore, this paper proposes a new PS selection method based on the Gaussian Mixture Model (GMM). Firstly, PS candidates (PSCs) are selected based on the pixels’ amplitude information. Then, the amplitude deviation index of each PSC is calculated, and each pixel’s probability values in different Gaussian distributions are acquired through iterations. Subsequently, the cluster types of pixels with larger probability values are designated as low-amplitude deviation pixels. Finally, the coherence coefficient and phase stability of low-amplitude deviation pixels are calculated. By comparing the probability values of each of the pixels in different Gaussian distributions, the cluster type with the larger probability, such as high-coherence pixels and high-phase stability pixels, is selected and designated as the final PS. Our analysis of the measured data revealed that the proposed method not only increased the number of PSs in the group, but also improved the stability of the number of PSs between groups.