[发明专利]一种自适应的高精度干涉SAR相位估计方法

权利要求书

1.一种自适应的高精度干涉SAR相位估计方法，其特征是它包含以下几个步骤：

步骤1、采用传统的最大相干系数法进行InSAR的图像粗配准：

采用传统的最大相干系数法对InSAR副复图像进行图像粗配准，得到粗配准后副图像S₂；具体方法是：在主图像S₁中选定任一像素点作为参考点，以参考点为中心取一个P×P大小的窗口；在配准前的副图像S′₂中对应位置处选取一个Q×Q大小的搜索窗口，这里Q？P保证在副图像中搜索到与主图像相匹配的P×P大小窗口，P表示主图像S₁中参考窗口的大小，Q表示配准前的副图像S′₂中搜索窗口的大小；

在配准前副图像S′₂中选定任一像素点作为参考点，以参考点为中心取一个P×P大小的窗口，在整个Q×Q大小的搜索窗口进行滑窗移动，计算每次滑窗后主副图像中P×P各像素点的复相干系数γ：

γ=Σi=1PΣj=1PS1′(i,j)S2′*(i+u,j+v)Σi=1PΣj=1P|S1′(i,j)|2Σi=1PΣj=1P|S2′(i+u,j+v)|2]]>

其中S′₁(i，j)表示主图像S₁中P×P窗口中第i行、第j列的像素值，S′₂(i+u，j+v)表示副图像S′₂中滑窗移动u行、v列的搜索窗口的像素值，γ每次滑窗后主副图像中P×P各像素点的复相干系数，*表示共轭，||表示求模值；最后通过每次滑窗后主副图像中P×P各像素点的复相干系数γ表达式计算参考窗口与搜索得到的各窗口的相干系数，每次滑窗后主副图像中P×P各像素点的复相干系数γ取最大值所对应的位置即为粗配准的误差偏移量(u，v)；根据得到的误差偏移位置(u，v)，对配准前的副图像S′₂中所有像素点移动u行、v列，得到配准前的副图像S′₂的粗配准图像S₂；

步骤2、建立维纳滤波器信号模型：

由步骤1得到的粗配准图像S₂，建立维纳滤波器的信号模型；具体过程如下：取粗配准图像S₂中对应像素i为中心的一个3×3的矩阵作为维纳滤波器的输入信号，建立如下滤波器模型：

构造回波信号矢量

S‾2(i)=[S2(i-4),S2(i-3),...,S2(i+4)]T]]>

其中S₂(i-4)，S₂(i-3)，…，S₂(i+4)表示粗配准图像S₂中3×3大小矩阵中的像素点，T表示矩阵的转置；

构造维纳滤波器参量ω(i)

ω(i)＝[ω₁，ω₂，…，ω₉]^H

其中ω₁，ω₂，…，ω₉分别表示维纳滤波器9个抽头的权系数，H表示矩阵的共轭转置；

将维纳滤波器参量ω(i)与回波信号矢量相乘得到主图像S₁中像素i的估计值s₂(i)，即s₂(i)＝ω(i)H；将得到的主图像S₁中像素i的估计值s₂(i)与主图像S₁中对应像素i的值s₁(i)进行相减，得到主副图像对应像素点i的估计误差值e(i)

e(i)=s2(i)-s1(i)=ω(i)HS‾2(i)-s1(i)]]>

其中H表示矩阵的共轭转置；e(i)的均方值最小时对应着最佳滤波器参量ω_opt(i)；

步骤3、计算维纳滤波器的最佳滤波器参量ω_opt(i)：

求取粗配准后副图像回波信号矢量的自相关矩阵求取粗配准后副图像回波信号矢量与主图像像素点s₁(i)的互相关矢量

上式中，E[]表示统计平均，H表示矩阵的共轭转置，*表示矩阵的共轭；

对粗配准后副图像回波信号矢量的自相关矩阵求逆，并与互相关矢量相乘，得到最佳滤波器参量ω_opt(i)，即其中-1表示求矩阵的逆；

步骤4、实现InSAR图像精配准

将步骤3得到的最佳滤波器参量ω_opt(i)，与由步骤2构造的回波信号矢量相乘，得到副图像精配准后像素点值s₂(i)

s2(i)=(ωopt(i))HS‾2(i)]]>

其中ω_opt(i)表示像素点i的最佳滤波器参量，是步骤2中副图像中对应像素i的3×3的矩阵构成的回波信号矢量，H表示矩阵的共轭转置；

步骤5、构造接收信号的协方差矩阵

利用步骤4精配准后副图像的相应像素点s₂(i)，构造联合最佳加权矢量s(i)，s(i)＝[s₁(i)，s₂(i)]^T，这里T表示矩阵的转置；

利用最佳加权矢量s(i)，构造接收信号协方差矩阵C_s(i)为：

其中E[]表示统计平均，表示干涉相位导向矢量，表示主副图像中对应像素点i之间的干涉相位，x(i)表示加权矢量s(i)中的信号矢量，v(i)表示加权矢量s(i)中加性高斯噪声矢量，R(i)＝E[x(i)x(i)^H]表示信号矢量x(i)的自相关矩阵，σ²表示噪声矢量v(i)的方差，I₂表示二阶的单位矩阵，e表示Hadamard积，H表示矩阵的共轭转置；

步骤6、利用MUSIC算法分离信号和噪声子空间

将步骤5得到的接收信号的协方差矩阵C_s(i)，利用MUSIC算法分解得到信号和噪声子空间；具体步骤如下：对协方差矩阵C_s(i)进行特征值分解，求取它的最小特征值λ_min，与最小特征值λ_min对应的特征向量为噪声子空间记为u₂；

步骤7、构造谱函数估计干涉相位

将噪声子空间记为E_N＝u₂，利用干涉相位导向矢量与噪声子空间相乘，构造谱函数这里表示干涉相位估计值，T表示矩阵的转置，在[0，2π)内改变相位的值，P_MUSIC峰值位置对应的相位即为主副图像中对应像素点i的干涉相位值