[发明专利]多基线干涉合成孔径雷达干涉相位展开方法

摘要

本发明公开了一种多基线干涉合成孔径雷达干涉相位展开方法。其过程为对卫星接收M幅回波数据分别进行SAR成像处理；选择SAR成像处理后的任意一幅图像作为主图像，以主图像为参考用相关法，分别对其它的图像进行粗配准，得到M－1幅粗配准SAR图像；利用粗配准后的SAR图像构造最优加权联合数据矢量jx(i,w_opt)；根据最优加权联合数据矢量估计协方差矩阵Cov_jx(i,w_opt)并对协方差矩阵进行特征分解，得出代价函数；将代价函数的输出功率最小值所对应的干涉相位作为干涉相位展开结果；对每一像素重复执行上述相关过程得到整个图像的干涉相位展开结果。本发明具有运算量小，在SAR图像配准精度很差时准确地展开干涉相位的优点，可用于对地面高分辨率、高精度三维地形重构和地面动目标检测。

主权项

1.一种多基线干涉合成孔径雷达干涉相位展开方法，包括如下过程：A.对卫星接收的M幅回波数据分别进行SAR成像处理，M≥3；B.选择SAR成像处理后的任意一幅图像作为主图像，以该主图像为参考用相关法，分别对其它的图像进行粗配准，得到M-1幅粗配准SAR图像；C.利用粗配准后的SAR图像构造最优加权联合数据矢量jx(i，w_opt)；jx(i,wopt)=[x1(i),x^2(i,wopt),x^3(i,wopt),···,x^M(i,wopt)]T]]>式中，x₁(i)为主图像，分别优化后的第2，第3，…，第M幅SAR复图像数据，w_opt为最优加权向量，i表示当前要计算的主图像和粗配准图像像素标号，上标T表示转置操作；D.根据最优加权联合数据矢量jx(i，w_opt)估计协方差矩阵并对协方差矩阵进行特征分解，得Covjx(i,wopt)=E[jx(i,wopt)jx(i,wopt)H]]]>式中，E[·]表示统计平均，上标H表示共轭转置；E.对协方差矩阵Cov_jx(i，w_opt)进行特征分解，得出代价函数：式中，张成协方差矩阵的信号子空间，{β_n^(l)}_l＝2^M张成噪声子空间；F.通过代价函数进行干涉相位展开，即将代价函数的输出功率最小值所对应的干涉相位作为干涉相位展开结果；G.对主图像和粗配准SAR图像的每一像素分别执行过程C～F操作，得到整个图像的干涉相位展开结果。