[发明专利]一种多通道SAR-GMTI地形干涉相位补偿方法及系统

权利要求书

1.一种多通道SAR-GMTI地形干涉相位补偿方法，其特征在于，包括：

步骤S1：获取经配准后的多通道SAR-GMTI图像，选取任意通道为参考通道，依次处理其他待补偿通道数据；

步骤S2：计算参考通道与待补偿通道数据的干涉相位，并通过二维带通滤波器处理干涉相位频谱提取干涉相位估计值，进而通过相位补偿更新待补偿通道数据；

步骤S3：通过重复上述操作多次实现待补偿通道数据的迭代更新，从而获取干涉相位直方图分布，实现多通道SAR-GMTI系统地形慢变干涉相位补偿；

判断所述直方图分布函数是否保持稳定或当前是否已进行预设次数迭代：

若不满足要求，则增加所述带通滤波器窗长，令重新进入步骤S3，继续进行迭代计算；

其中s_m,proc(x,y)表示第m个待处理通道图像数据；m＝1,2,… N-1，其中N表示多通道SAR系统的通道个数；表示对第m通道SAR图像数据s_m,proc(x,y)进行相位补偿得到更新后的通道数据；

若满足所述要求，则令m＝m+1，进入步骤S3，直至处理完成所有待处理通道的地形干涉相位补偿，输出补偿后的多通道SAR-GMTI距离-多普勒数据结果。

2.根据权利要求1所述的多通道SAR-GMTI地形干涉相位补偿方法，其特征在于，所述干涉相位计算表达式为：

式中，为二维干涉相位函数，表示对参考通道的SAR图像数据取共轭后的数据，s_comp,m(t,f_a)，m＝2,3,…,N-1，表示其他待补偿通道SAR图像数据，其中N表示多通道SAR系统的通道个数，t表示距离快时间变量，f_a为多普勒频域分量。

3.根据权利要求2所述的多通道SAR-GMTI地形干涉相位补偿方法，其特征在于，对所述干涉相位函数进行二维傅里叶变换获取二维频谱，设置带通滤波器提取干涉相位频谱峰值区域，获取干涉相位估计值。

4.根据权利要求3所述的多通道SAR-GMTI地形干涉相位补偿方法，其特征在于，根据所述干涉相位估计值，通过相位补偿更新待补偿通道SAR图像数据。

5.根据权利要求1所述的多通道SAR-GMTI地形干涉相位补偿方法，其特征在于，通过重复带通滤波提取干涉相位的方式迭代更新待补偿通道SAR图像数据，在每次迭代中，依次更新干涉相位函数及其频谱、带通滤波器大小、干涉相位估计值以及待补偿SAR图像数据。

6.一种多通道SAR-GMTI地形干涉相位补偿系统，其特征在于，包括：

模块M1：获取经配准后的多通道SAR-GMTI图像，选取任意通道为参考通道，依次处理其他待补偿通道数据；

模块M2：计算参考通道与待补偿通道数据的干涉相位，并通过二维带通滤波器处理干涉相位频谱提取干涉相位估计值，进而通过相位补偿更新待补偿通道数据；

模块M3：通过重复上述操作多次实现待补偿通道数据的迭代更新，从而获取干涉相位直方图分布，实现多通道SAR-GMTI系统地形慢变干涉相位补偿；

判断所述直方图分布函数是否保持稳定或当前是否已进行预设次数迭代：

若不满足要求，则增加所述带通滤波器窗长，令重新进入步骤S3，继续进行迭代计算；

其中s_m,proc(x,y)表示第m个待处理通道图像数据；m＝1,2,…N-1，其中N表示多通道SAR系统的通道个数；表示对第m通道SAR图像数据s_m,proc(x,y)进行相位补偿得到更新后的通道数据；

若满足所述要求，则令m＝m+1，进入步骤S3，直至处理完成所有待处理通道的地形干涉相位补偿，输出补偿后的多通道SAR-GMTI距离-多普勒数据结果。