基于二分法的山区PS-InSAR形变参数估计

发布时间：2021-08-01 00:10

　　针对山区PS-InSAR第一次模型反演时需使估形变速率及高程残差进一步逼近最优取值的问题,提出了一种基于二分法的参数估计方法。其分为四步:第一步,固定其他参数,对待定参数取过大值a和过小值b以及二者的平均值c分别模型反演;第二步,依据空间连续性、分块性、椒盐噪声和均匀噪声等质量标准评判三组结果图的优劣次序;第三步,根据质量递减原则,判断最优值在前半段或后半段,之后将a、c或c、b作为新的过大、过小值,重复前面的步骤,直至结果图质量稳定;第四步,固定该参数,并按步骤一至三依次找出其他参数的最优取值。通过丽江市某滑坡区的PS-InSAR反演试验并与SBAS-InSAR得到的估计初值比较,证明了本方法是一种可行的使估计参数收敛的方法。 

【文章来源】：软件. 2020,41(04)

【文章页数】：7 页

【部分图文】：

新建地区影像图

严格来讲，引言中提出了三个问题：（1）调参幅度有限导致结果变化不明显；（2）即使变化明显，能归纳出什么样的标准；（3）即使标准明确，怎样使参数快速收敛。前文提到，当参数变化很小时，中间结果的变化往往也是细微的。但反过来想，当参数变化很大时，是否就意味着中间结果的变化会很大，进而可以确定一个标准，作为评判结果图像质量变好（差）的依据。因而下面的实验展示的是一个如何从零向最优取值逼近的过程，而非从初值（即SBAS-InSAR的结果）开始估计，就是为了可以观察到更明显的中间结果变化过程，以更好地理解质量判定标准的由来。但使用这些准则时，就需要从中间结果的细微变化中仔细区分了。

对于城市地区，当不确定有明显沉降时，形变速率和残余高度一般设置为±25 mm/year和±70 m，当有已知沉降或高层建筑物较多时，可以适当增加取值至±50 mm/year和±100 m。而新建地区地处峡谷深处，鲜有人烟，只有散落的低矮民房，因而将残余高度暂设为±20 m。由于两景影像可探测到的最大形变为波长的一半（28.15 mm），换算成年为867.6 mm/year，故将其设为过大值。这里将形变速率参数设为±400 mm/year、±500 mm/year、±867 mm/year，得到的形变速率图和高程残差图如下。图3 b不同参数的高程残差图(b1～b3:?20 m)


本文编号：3314450