The goafs caused by coal mining cause great harm to the surface farmland, buildings, and personal safety. The existing monitoring methods cost a lot of workforce and material resources. Therefore, this paper proposes an inversion approach for establishing the locations of underground goafs and the parameters of the probability integral method (PIM), thus integrating distributed scatter interferometric synthetic aperture radar (DS-InSAR) data and the PIM. Firstly, a large amount of surface deformation observation data above the goaf are obtained by DS-InSAR, and the line-of-sight deformation is regarded as the true value. Secondly, according to the obtained surface deformations, the ranges of eight goaf location parameters and three PIM parameters are set. Thirdly, a correlation function between the surface deformation and the underground goaf location is constructed. Finally, a particle swarm optimization algorithm is used to search for the optimal parameters in the range of the set parameters to meet the requirement for minimum error between the surface deformation calculated by PIM and the line-of-sight deformation obtained by DS-InSAR. These optimal parameters are thus regarded as the real values of the position of the underground goaf and the PIM parameters. The simulation results show that the maximum relative error between the position of the goaf and the PIM parameters is 2.11%. Taking the 93,604 working face of the Zhangshuanglou coal mine in the Peibei mining area as the research object and 12 Sentinel-1A images as the data source, the goaf location and PIM parameters of the working face were successfully inverted. The inversion results show that the maximum relative error in the goaf location parameters was 16.61%, and the maximum relative error in the PIM parameters was 26.67%.