Although the modified Goldstein filter based on the local signal-to-noise （SNR） has been proved to be superior to the classical Goldstein and Baran filters with more comprehensive filter parameter, its adaptation is not always sufficient in the reduction of phase noise. In this work, the local SNR-based Goldstein filter is further developed with the improvements in the definition of the local SNR and the adaption of the filtering patch size. What＇s more, for preventing the loss of the phase signal caused by the excessive filtering, an iteration filtering operation is also introduced in this new algorithm. To evaluate the performance of the proposed algorithm, both a simulated digital elevation model （DEM） interferogram and real SAR deformation interferogram spanning the L＇ Aquila earthquake are carried out. The quantitative results from the simulated and real data reveal that up to 79.5% noises can be reduced by the new filter, indicating 9%-32% improvements over the previous local SNR-based Goldstein filter. This demonstrates that the new filter is not only equipped with sufficient adaption, but also can suppress the phase noise without the sacrifice of the phase signal.
SUN QianLI Zhi-weiZHU Jian-junDING Xiao-liHU JunXU Bing
The subsidence of the mining area was monitored by analyzing the phase of permanent scatters （PS） which maintained high coherence in magnitude of SAR images.A hew method of spatial unwrapping was presented which used the subsidence rates calculated on comer reflector （CR） points as constraints for PS network to perform the spatial unwrapping using the parametric adjustment method.The algorithm achieved the integration of CR data and PSInSAR algorithm.The colliery dense distributed area around Baisha reservoir was chosen as the study area in the experiment.The time series of subsidence from February in 2007 to February in 2010 is successfully inversed by using the periodic function to simulate the linear and nonlinear components of the deformation.The simulation results show that the accuracy can be ± 2.1 mm with the leveling data being used as the external validation data.
To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to derive the coseismic slip distribution of the earthquake fault. Firstly, based on the homogeneous elastic half-space model, the fault geometric parameters were solved by the genetic algorithm. The best fitting model shows that the fault is a 13.7 km×14.1 km rectangular fault, in 139.3° strike direction and 50.2° southwest-dipping. Secondly, fixing the optimal fault geometric parameters, the fault plane was extended and discretized into 16× 16 patches, each with a size of 1 kmx 1 krn, and the non-uniform slip distribution of the fault was inverted by the steepest descent method with an appropriate smoothing ratio based on the layered crustal structure model. The preferred solution shows that the fault is mainly a normal fault with slight right-lateral strike slip, the maximum slip of 1.01 m is located in the depth of 8.28 km, the average rake is -100.9°, and the total geodetic moment is about 3.34× 1018 N.m （Mw 6.28）. The results are much closer than previous studies in comparison with the seismological estimation. These demonstrate that the coseismic fault slip distribution of the L＇Aauila earthauake inverted by the crustal model considering layered characters is reliable.
WANG Yong-zheZHU Jian-junOU Zi-qiangLI Zhi-weiXING Xue-min
