[发明专利]一种基于最大周期图的PS-InSAR精确搜索方法

权利要求书

1.一种基于最大周期图的PS-InSAR精确搜索方法，其特征在于：包括以下步骤：

步骤一：对长时间序列差分干涉图像进行相位滤波，利用离差振幅准则筛选出永久散射体像素；

步骤二：基于最大周期图准则粗估计永久散射体的高程与形变速率值；具体为：PS-InSAR解决的是以下目标函数的估计问题，

其中，M为长时间序列SAR单视复图像的总幅数，λ为电磁波的波长，g_m为第m幅差分干涉SAR图像，b_⊥m为垂直基线，r和θ为PS点与主卫星天线对应的斜矩和侧视角，t_m为辅图像与主图像的成像时间间隔，h为PS点的高程值，v为PS点沿雷达视线方向上的形变速率值，为了精确估计以上时序相关系数γ_p(h,v)最大时的高程h及形变速率v，首先，利用最大化周期图法，解出每一幅长时间序列差分干涉图像的周期图谱估计：

经过傅里叶变换将复信号从t-b空间变换到h-v空间，然后将频谱峰值处对应的高程和形变速率作为最大化时序相关系数这一目标函数的初步估计值；

步骤三：以最大周期图粗估计的结果作为初始值，基于最大时序相关系数准则无约束精确搜索高程与形变速率值，具体为：

重新建立目标函数，

f_p(h,v)＝-|γ_p(h,v)|²(1.5)

γ_p(h,v)为时序相关系数，h为PS点的高程值，v为PS点沿雷达视线方向上的形变速率值，将目标函数γ_p的最大化问题转化为f_p的最小化问题，采用拟牛顿法，

其中，以d为搜索方向，H_i+1为f_p的Hesse矩阵的逆矩阵的一个近似矩阵，将初始值作为搜索的起始点，调用Matlab里的fminunc求解最优问题，以获得精确的高程和形变速率测量值。

2.根据权利要求1所述的一种基于最大周期图的PS-InSAR精确搜索方法，其特征在于：步骤一具体为：对M幅长时间序列SAR单视复图像进行配准、差分干涉操作以及差分干涉相位滤波，得到处理后的差分干涉SAR图像，接着，计算出每一幅观测图像的归一化幅度方差μ为：

其中，M为长时间序列SAR单视复图像的总幅数，g_m为第m幅差分干涉SAR图像，而振幅离差指数D_A与归一化振幅方差μ存在以下关系：

D_A＝1/(1-μ)-1 (1.2)

设定一定的阈值k，某个像素的振幅离差指数小于阈值，则将该像素选为永久散射体像素。