Absolute phase is equivalent to absolute range【1】.所以形变测量的方法也可以按照这两类来分。
一类是基于相位的:InSAR的方法和频谱分割的方法(MAI, Split-bandwidth, Spectrum diversity, Multifrequency processing在我看来都一样,MAI是Zebker提出Jung深挖,DLR的delta-K相关论文大多是频谱分割的,如果非要说区别,可能就是MAI是方位向的自孔径,频谱分割时方位向距离向都可以),它的误差来源和相位有关系,影响相位的都会引入噪声
另一类是直接测量距离的,就是OT(简称OT)。测量像素的偏移量,本质上还是直接去测量距离差,而不是通过相位来转化。方法的原理上与配准一样,计算主辅影像的互相关系数,最大值的点就是同名的,这又分成基于复数据的和幅度数据的incoherent cross-correlation.
【1】中有各个方法详细的理论精度定义。重点来讨论OT的发现。常规OT来源于配准,问题在于:1. 计算量巨大(祁连山一个TSX-SM像对跑了4天);2. 结果中会出现patch-like artifact【2】的现象(受区域内强散射点的影响), band-like artifact【2】( azimuth offset result腾哥的论文中是说这一现象是Squint-angle的变化导致的,很难建模去除。我感觉,它和MAI论文中的电离层误差很像,不知道是怎么回事。)和cross-hair artifact【3】。3.会存在很多异常值(原因可能是)
腾哥的PT OT可以在一定程度上解决计算速度的问题和patch-like artifact【2】,Zebker他们【3】是用高斯滤波器消除cross-hair artifact以及异常值。赵chaoying是用窗口的标准差阈值来滤除异常值。此外剔除异常值的方法还有中值滤波,以及简单粗暴的直接设置一个形变量阈值(腾哥)
为什么会有这么多异常值呢?一方面是配准错误,另一方面是地物的变化。地物的变化对基于幅度的OT来说是致命的。因为它确实比较简单,玩不出什么花来,所以方法上的创新不多,一是SBAS,二是PT,最近中南的汪师兄发了一篇RS,将它拓展到全极化数据中用在了美国那个著名的滑坡上,在滑坡边缘上会比较好一点。个人感觉还是有一些东西可以做点文章的。
好吧~先就这些了~
【1】Bamler R, Eineder M. Accuracy of differential shift estimation by correlation and split-bandwidth interferometry for wideband and delta-k SAR systems[J]. IEEE Geoscience & Remote Sensing Letters, 2005, 2(2):151-155.
【2】Wang T, Jonsson S. Improved SAR Amplitude Image Offset Measurements for Deriving Three-Dimensional Coseismic Displacements[J]. IEEE Journal of Selected Topics in Applied Earth Observations & Remote Sensing, 2015, 8(7):1-8.
【3】Yun S H, Zebker H, Segall P, et al. Interferogram formation in the presence of complex and large deformation[J]. Geophysical Research Letters, 2007, 34(12):237-254.