[发明专利]高程层叠情况下稳健的绝对干涉相位估计方法

权利要求书

1.一种高程层叠情况下稳健的绝对干涉相位估计方法，其特征在于：所述的绝对干涉相位估计方法包括按顺序进行的下列步骤： 

1)对合成孔径雷达图像进行粗配准，以构造观测矢量； 

2)确定存在配准误差时的协方差矩阵和真实导向矢量； 

3)用估计的协方差矩阵和真实导向矢量进行波束扫描，估计绝对干涉相位。 

2.根据权利要求1所述的高程层叠情况下稳健的绝对干涉相位估计方法，其特征在于：在步骤1)中，所述的构造观测矢量方法是： 

用传统的相关法以SAR主图像s₁作为参考，对副图像1进行粗配准得到SAR粗配准图像s₂，对副图像2进行粗配准得到SAR粗配准图像s₃，利用SAR主图像s₁、SAR粗配准图像s₂和SAR粗配准图像s₃构造观测矢量；对于SAR主图像s₁中的一个像素i，则构造出的观测矢量表示为： 

js(i)＝[s₁(i)，s₂(i-3)，s₂(i-2)，s₂(i)，s₂(i+1)，s₃(i-1)，s₃(i)，s₃(i+2)，s₃(i+3)]^T。 

3.根据权利要求1所述的高程层叠情况下稳健的绝对干涉相位估计方法，其特征在于：在步骤2)中，所述的确定存在配准误差时的协方差矩阵和真实导向矢量的方法是： 

a.利用步骤1)中得到的观测矢量，估计其协方差矩阵C_js(i)，表示为： 

C_js(i)＝E{js(i)js^H(i)} 

＝E{[s₁(i)，s₂(i-3)，s₂(i-2)，s₂(i)，s₂(i+1)，s₃(i-1)，s₃(i)，s₃(i+2)，s₃(i+3)]^T[s₁(i)，s₂(i-3)，s₂(i-2)，s₂(i)，s₂(i+1)，s₃(i-1)，s₃(i)，s₃(i+2)，s₃(i+3)]^*} 

b.对于像素i存在配准误差时的真实导向矢量β(φ_ni)，其中n＝1，2，...，N_s，N_s表示layover的个数，假设N_s＝2，表示为： 

(n＝1，2，...，N_s) 

其中， 

为SAR主图像s₁中像素i与SAR粗配准图像s₂中像素p的相关系数； 

为SAR主图像s₁中像素i与SAR粗配准图像s₃中像素l的相关系数； 

其中，2K+1为从相邻像素中获得的独立同分布样本数；B_1m(m＝2，3)为第m颗卫星与作为参考卫星的第1颗卫星之间的垂直航向基线，φ_ni(n＝1，2，...，N_s)为第二颗卫星和参考卫星之间的绝对干涉相位，T表示转置操作，E{}表示统计平均，上标H表示共轭转置，*表示共轭操作。 

4.根据权利要求1所述的高程层叠情况下稳健的绝对干涉相位估计方法，其特征在于：在步骤3)中，用估计的协方差矩阵和真实导向矢量进行波束扫描的过程如下： 

a.以0.01弧度为步长，将φ_ni(n＝1，2，...，N_s)在[-π，π]范围内的每一个取值代入代价函数J＝β^H(φ_ni)C_js(i)β(φ_ni)(n＝1，2，...，N_s)中，分别计算其结果； 

其中， 

(n＝1，2，...，N_s) 

b.对所计算的结果进行排列，找出其中对应的N_s个最大值时所对应的干涉相位即为绝对干涉相位估计结果。