An Improved Phase Filter for Differential SAR Interferometry Based on an Iterative Method

Phase quality is a key element in the analysis of the deformation of the Earth's surface carried out with differential synthetic aperture radar interferometry. Various decorrelation sources may degrade the surface deformation estimates, and thus, phase filters are needed for this kind of application. The well-known Goldstein filter is the most widely used due to its simple implementation and computational efficiency. In the past years, improved filters have been proposed, which are based on this filter but introduce variations in the data processing. The effectiveness of these filters mostly depends on the size of the filtering window, the weight of the smoothed spectrum, and the kernel used to filter the spectrum. In this paper, we evaluate the performance of four of these filters and present a new method that outperforms all of them. The proposed filter is based on an iterative method in which the original phase is denoised progressively with adaptive filtering windows of different sizes. The effectiveness of the filter is controlled by the interferometric coherence, a direct indicator of the phase quality. Moreover, we introduce some modifications regarding the processing of the power spectrum. Specifically, we propose to smooth the original phase using a new filter which is based on a Chebyshev interpolation scheme. The performance of the new filter has been tested on both simulated and real interferograms, acquired by RADARSAT-2 and the Uninhabited Aerial Vehicle Synthetic Aperture Radar, which mapped two different geological events that caused surface deformation.This work was supported in part by the Spanish Ministry of Economy, Industry and Competitiveness, in part by the State Agency of Research (AEI), in part by the European Funds for Regional Development under Project TIN2014-55413-C2-2-P and Project TEC2017-85244-C2-1-P, in part by the U.K. Natural Environmental Research Council through the Looking Inside the Continents under Grant NE/K011006/1, in part by the Rapid deployment of a seismic array in Ecuador following the April 16th 2016 M7.8 Pedernales earthquake under Grant NE/P008828/1, and in part by the Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics under Grant COMET, GA/13/M/031