[发明专利]一种基于非线性优化策略的SBAS‑DInSAR方法

权利要求书

1.一种基于非线性优化策略的SBAS-DInSAR方法，其特征在于，包括以下几个步骤：

步骤一：输入包含地表形变信息的长时间SAR图像序列，设共有K+1幅SAR图像，获取时刻分别为t₀,t₁,…,t_K，其中K是正整数；对SAR图像进行预处理，包括配准、干涉、去除参考面相位、去除地形相位，获得差分干涉相位图像序列；设所述差分干涉相位图像序列包括M幅差分干涉相位图像，其中M是正整数；第j幅差分干涉相位图像中主图像和辅图像的获取时刻分别为和j是小于或等于M的正整数，IE_j和IS_j是小于或等于K的非负整数，且IE_j＞IS_j；

步骤二：对所述差分干涉相位图像序列，基于振幅信息提取高相干点，设高相干点的数量为H，其中H是正整数；第h个高相干点记为x_h，第j幅差分干涉相位图像中第h个高相干点的非缠绕差分干涉相位记为其中h是小于或等于H的正整数；

步骤三：对第j幅差分干涉相位图像的第h个高相干点x_h有：

ΔΦjuw(xh)=Φuw(xh,tIEj)-Φuw(xh,tISj)]]>

其中，和分别为高相干点x_h在时刻和时刻由于地表形变引起的相位；

步骤四：建立目标函数：

f(ζh)=1MlT·exp[i(4πλD·ζh-ΔΦhw)]]]>

式中，

其中，i是虚数单位；(D)_pq为矩阵D第p行第q列的元素，p为大于零且小于或等于M的正整数，q为大于零且小于或等于K的正整数,t_q和t_q-1分别是第q幅和第q-1幅SAR图像的获取时刻，IS_p和IE_p分别是第p幅差分干涉相位图像中主图像和辅图像的序号，第p幅差分干涉相位图像由获取时刻为和的SAR图像干涉生成；λ为雷达波长；ζ_h为目标函数f(ζ_h)的自变量；为由1构成的列向量；的第j个分量表示第j幅差分干涉相位图像中第h个高相干点的缠绕差分干涉相位；

步骤五：通过拟牛顿法对所述目标函数进行第一次优化，提取高相干点的线性形变速率；

步骤六：以线性形变速率作为迭代初值，对所述目标函数进行第二次优化，计算总形变速率。