[发明专利]基于复数后向投影的自旋目标三维成像方法

说明书

技术领域

本发明属于信号处理技术领域，更进一步涉及雷达成像领域中的基于复数后向投影的自旋目标三维成像方法。本发明可以有效地对空间目标进行三维成像，并对目标形状、尺寸等特征进行精确描述，为后续的目标识别提供有力保障。

背景技术

当目标高速自旋时，传统的三维成像方法已经不再适用。在观测时间内，散射点相对于成像中心的距离和多普勒为时变量，较少的有效积累数据将导致图像质量下降，甚至使距离-瞬时多普勒算法失效。因此，有必要根据自旋目标的运动特征研究新的成像方法。

Qi Wang等人在文献“High-Resolution Three-Dimensional Radar Imaging for Rapidly Spinning Targets”(IEEE trans.on GRS，VOL.46，no.1，pp.22-30，2008)中提出将广义Radon变换(GRT)与CLEAN技术相结合进行目标三维成像的方法。该方法利用旋转点目标回波在距离慢时间域的周期性变化规律，通过检测正弦曲线估计散射点空间位置，并结合修正的CLEAN技术对目标进行三维成像。但该方法存在的不足是，广义Radon变换方法所得图像的主瓣较宽，分辨率低，影响对空间目标散射点分布特征的准确描述。

发明内容

本发明的目的在于克服现有技术的不足，提出一种基于复数后向投影的自旋目标三维成像方法。该方法弥补了广义Radon变换方法具有存在距离像偏移和展宽，方位像散焦，无法对目标平动进行准确补偿，成像旁瓣较高的不足，充分利用自旋散射点回波在距离-慢时间域具有正弦包络和相位这一特征，采用宽带复数后向投影算法沿不同的正弦曲线进行相干积分，并结合修正的CLEAN技术对目标进行三维成像。

实现本发明的基本思路是：首先将自旋目标的ISAR回波进行回波相干化，以补偿目标高速运动对回波产生的影响，接着以图像熵最小为准则，通过建立有效的搜索算法对目标回波的平动量进行精细补偿以得到转台模型下的目标回波然后根据不同尺寸、不同自旋频率目标对脉冲重复频率的要求，进行二维成像算法的选择，最后对目标沿转轴方向采用相应的二维成像方法对所有二维切片进行成像，最终结合修正的CLEAN技术对目标进行三维成像。

本发明的具体步骤如下：

1.基于复数后向投影的自旋目标三维成像方法，包括如下步骤：

(1)雷达录取ISAR回波；

(2)回波相干化

2a)在距离频域对回波距离向的平移和展宽分别补偿相位因子；

2b)在慢时间域对回波补偿相位因子，使得回波完全相干化；

(3)平动补偿

3a)采用傅里叶变换对回波进行距离脉压，获得距离-慢时间域回波；

3b)采用某单次距离-慢时间域回波，分别计算其与其它各次距离-慢时间域回波的滑动相关系数，具有最大滑动相关系数的回波与该回波之间的时间间隔为目标的自旋周期，对自旋周期取倒数并乘以2π，得到自旋角频率；

3c)对应最大滑动相关系数的回波与该次回波之间移动的距离数为一个周期内目标的平动量，补偿平动量实现包络粗对齐；

3d)从第二次距离-慢时间域回波开始，对于单次距离-慢时间域回波，采用广义Radon变换方法获取不同偏移量下该次回波与之前所有回波对应的二维图像并求图像熵，使图像熵最小的距离偏移量为最优距离偏移量，根据最优偏移量对回波包络进行微调，将每次回波进行包络精对齐后，最终得到平动补偿后的目标回波；

(4)算法选择

4a)采用距离-多普勒算法对平动补偿后的回波成像，得到目标的距离-多普勒图像；

4b)若距离-多普勒图像的方位维没有混叠，则选择CIRT算法，执行步骤(5)；

4c)若距离-多普勒图像出现方位维混叠，则选择RIRT算法；

(5)二维切片成像

5a)建立自旋目标的柱面坐标系，沿目标自转轴方向将自旋目标分成多个二维切片；

5b)将二维切片按照步骤(4)选择的算法进行成像，得到二维切片图像；

(6)重构回波并更新

6a)对于所有的二维切片图像，提取图像域中最大峰值对应的三维坐标，记录为目标散射点的位置坐标，并在距离-慢时间域构造散射点回波；

6b)运用最小范数准则对步骤6a)中所构造的散射点回波的幅度进行估计；

6c)将记录的散射点的位置坐标与估计的回波幅度组合构成散射点回波；

6d)从步骤3d)得到的目标回波中减去散射点回波，更新目标回波；