As a type of earth observation technology, interferometric synthetic aperture radar (InSAR) is increasingly widely used in the field of geological disaster detection. However, the application of InSAR in low-coherence areas, such as alpine canyon areas and vegetation coverage areas, is subject to considerable limitations. How to accurately identify landslides from InSAR measurement data in these areas remains the subject of several challenges and shortcomings. Based on statistical analysis and spatial cluster analysis, in this paper, we propose an automatic landslide identification and gradation method suitable for low-coherence areas. The proposed method combines the small baseline subset InSAR (SBAS-InSAR) method and the interferogram stacking (stacking-InSAR) method to obtain a deformation map in the study area, using statistical analysis and spatial cluster analysis to extract deformation regions and landslide polygons to propose a landslide screening model (LSM) based on multivariate features to screen landslides and reduce the interference of noise in landslide identification, in addition to proposing a landslide gradation model (LGM) based on signum function to grade the identified landslides and provide support to distinguish landslides with different deformation degrees. The method was applied to landslide identification in the upper section of the Jinsha River basin, and 47 potential landslides were identified, including 15 high-risk landslides and 13 landslides endangering villages. The experimental results show that the proposed method can identify landslides accurately and hierarchically in low-coherence areas, providing support for geological hazard investigation agencies and local departments.