[发明专利]基于相关加权的干涉合成孔径雷达干涉相位估计方法

权利要求书

1.一种基于相关加权的干涉合成孔径雷达干涉相位估计方法，包括如下过程：

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

B.利用主图像(1)和粗配准图像(3)，构造相关加权观测矢量si(i，w(i))，

si(i，w(i))＝[s₁(i)，s_3W(i，w(i))]^T＝[s₁(i)，w^T(i)s_3v(i)]^T

式中s_3W(i，w(i))＝w^T(i)s_3v(i)

    w(i)＝[r₃₁(i-4，i)，r₃₁(i-3，i)，r₃₁(i-2，i)，r₃₁(i-1，i)，

                r₃₁(i，i)，r₃₁(i+1，i)，r₃₁(i+2，i)，r₃₁(i+3，i)，r₃₁(i+4，i)]^T

s_3v(i)＝[s₃(i-4)，s₃(i-3)，s₃(i-2)，s₃(i-1)，

                s_3v(i)，s₃(i+1)，s₃(i+2)，s₃(i+3)，s₃(i+4)]^T

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

m=i-4，i-3，i-2，i-1，i，i+1，i+2，i+3，i+4

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

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

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

上标T表示转置操作，上标·表示共轭操作；

C.根据相关加权观测矢量估计协方差矩阵为

C_si(i，w(i))＝E{si(i，w(i))si^H(i，w(i))}

式中E表示统计平均，上标H表示共轭转置操作；

D.由协方差矩阵得到代价函数：

J_BF＝a^H(φ_i)C_si(i，w(i))a(φ_i)

式中a(φi)=[1,ejφi]T]]>称为导向矢量，φ_i的取值范围是[-π，π]；

E.通过代价函数估计干涉相位，即用导向矢量和协方差矩阵进行波束扫描，并将波束形成的输出功率最大值所对应的干涉相位作为干涉相位估计结果。