[发明专利]基于相关加权联合子空间投影的InSAR干涉相位估计方法

权利要求书

1.一种基于相关加权联合子空间投影的InSAR干涉相位估计方法，其特征在于：所述的InSAR干涉相位估计方法包括按顺序进行的下列步骤：

1)用传统的相关法以主图像作为参考对副图像进行粗配准，得到粗配准图像；

2)利用主图像和粗配准图像构造相关加权观测矢量；

3)计算相关加权联合协方差矩阵和加权联合相关函数矩阵；

4)计算噪声子空间和信号子空间；

5)用噪声子空间投影来构造代价函数，并通过代价函数估计干涉相位。

2.根据权利要求1所述的基于相关加权联合子空间投影的InSAR干涉相位估计方法，其特征在于：在步骤2)中，所述的构造相关加权观测矢量的方法是：

a.分别计算主图像中像素a与粗配准图像中像素m的相关系数：

r31(m,a)=|Σk=-KKs3(m+k)s1*(a+k)|Σk=-KK|s3(m+k)|2Σk=-KK|s1(i+k)|2m=a-5,a-4,a-3,a-1,a,a+1,a+3,a+4,a+5a=i-1,i,i+3,i+4]]>

b.用所得到的相关系数r₃₁(m，a)分别对粗配准图像中的像素m进行加权，得到元素s_3W(a，w)，公式如下：

s3W(a,w)=waTVa]]>

其中：

V_a＝[s₂(a-5)，s₂(a-4)，s₂(a-3)，s₂(a-1)，s₂(a)，

        s₂(a+1)，s₂(a+3)，s₂(a+4)，s₂(a+5)]^T

                          (a＝i-1，i，i+3，i+4)

w_a＝[r₃₁(a-5，a)，r₃₁(a-4，a)，r₃₁(a-3，a)，r₃₁(a-1，a)，

       r₃₁(a，a)，r₃₁(a+1，a)，r₃₁(a+3，a)，r₃₁(a+4，a)，r₃₁(a+5，a)]^T

                          (a＝i-1，i，i+3，i+4)

c.用主图像中的像素a(a＝i-1，i，i+3，i+4)和得到的元素s_3W(a，w)构造相关加权观测矢量：

js(i,w)=[s1(i-1),s2W(i-1),s1(i),s2W(i),s1(i+3),s2W(i+3),s1(i+4),s2W(i+4)]T]]>

=[s1(i-1),wi-1TVi-1,s1(i),wiTVi,s1(i+3),wi+3TVi+3,s1(i+4),wi+4TVi+4]T]]>

其中：s₁和s₃分别为主图像和粗配准图像所对应的复图像数据，

i表示当前要计算的主图像和粗配准图像的像素标号，

2K+1为从相邻像素中获得的独立同分布样本数，

上标T表示转置操作，上标*表示共轭操作。