A method for Interferometric Synthetic Aperture Radar, InSAR, interferometric phase filtering in a wavelet domain in conjunction with local frequency estimation, the method comprising: Step A: transforming an InSAR interferometric phase φ into e jφ in a complex field, and determining a real part and an imaginary part of the interferometric phase e jφ in the complex field, respectively; Step B: performing the local frequency estimation on the interferometric phase φ to determine a frequency range in which the interferometric phase is located; Step C: performing wavelet decomposition with a scale of s on the real part and the imaginary part of the interferometric phase in the complex field, respectively, to determine frequency ranges of wavelet coefficients c m, n i of different sub-bands, wherein m, n are positions of the wavelet coefficients, and i is the scale of the decomposition within a range of 1 to s; Step D: determining, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, sub-bands in which useful information and noise are located, respectively, according to the frequency range in which the interferometric phase is located and the frequency ranges of the wavelet coefficients of different sub-bands; Step E: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, a general threshold shrinkage processing on the wavelet coefficients of the sub-bands in which the noise is located, respectively; Step F: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, a neighborhood threshold shrinkage processing on the wavelet coefficients of the sub-bands in which the useful information is located, respectively; Step G: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, wavelet reconstruction on the wavelet coefficients of the sub-bands in which the noise is located after the general threshold shrinkage processing and the wavelet coefficients of the sub-bands in which the useful information is located after the neighborhood threshold shrinkage processing collectively, to determine a real part and an imaginary part of a filtered interferometric phase in the complex field; and Step H: determining the filtered InSAR interferometric phase from the real part and the imaginary part of the filtered interferometric phase in the complex field. Step A: transforming an InSAR interferometric phase φ into e jφ in a complex field, and determining a real part and an imaginary part of the interferometric phase e jφ in the complex field, respectively; Step B: performing the local frequency estimation on the interferometric phase φ to determine a frequency range in which the interferometric phase is located; Step C: performing wavelet decomposition with a scale of s on the real part and the imaginary part of the interferometric phase in the complex field, respectively, to determine frequency ranges of wavelet coefficients c m, n i of different sub-bands, wherein m, n are positions of the wavelet coefficients, and i is the scale of the decomposition within a range of 1 to s; Step D: determining, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, sub-bands in which useful information and noise are located, respectively, according to the frequency range in which the interferometric phase is located and the frequency ranges of the wavelet coefficients of different sub-bands; Step E: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, a general threshold shrinkage processing on the wavelet coefficients of the sub-bands in which the noise is located, respectively; Step F: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, a neighborhood threshold shrinkage processing on the wavelet coefficients of the sub-bands in which the useful information is located, respectively; Step G: performing, for the real part and the imaginary part of the interferometric phase in the complex field, respectively, wavelet reconstruction on the wavelet coefficients of the sub-bands in which the noise is located after the general threshold shrinkage processing and the wavelet coefficients of the sub-bands in which the useful information is located after the neighborhood threshold shrinkage processing collectively, to determine a real part and an imaginary part of a filtered interferometric phase in the complex field; and Step H: determining the filtered InSAR interferometric phase from the real part and the imaginary part of the filtered interferometric phase in the complex field.