亚星游戏官网-www.yaxin868.com

专利名称：基于复数后向投影的自旋目标三维成像方法
技术领域：
本发明属于信号处理技术领域，更进一步涉及雷达成像领域中的基于复数后向投影的自旋目标三维成像方法。本发明可以有效地对空间目标进行三维成像，并对目标形状、 尺寸等特征进行精确描述，为后续的目标识别提供有力保障。
背景技术：
当目标高速自旋时，传统的三维成像方法已经不再适用。在观测时间内，散射点相对于成像中心的距离和多普勒为时变量，较少的有效积累数据将导致图像质量下降，甚至使距离-瞬时多普勒算法失效。因此，有必要根据自旋目标的运动特征研究新的成像方法。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)采用傅里叶变换对回波进行距离脉压，获得距离-慢时间域回波；
4
3b)采用某单次距离-慢时间域回波，分别计算其与其它各次距离-慢时间域回波的滑动相关系数，具有最大滑动相关系数的回波与该回波之间的时间间隔为目标的自旋周期，对自旋周期取倒数并乘以2 π，得到自旋角频率；3c)对应最大滑动相关系数的回波与该次回波之间移动的距离数为一个周期内目标的平动量，补偿平动量实现包络粗对齐；3d)从第二次距离-慢时间域回波开始，对于单次距离-慢时间域回波，采用广义 Radon变换方法获取不同偏移量下该次回波与之前所有回波对应的二维图像并求图像熵， 使图像熵最小的距离偏移量为最优距离偏移量，根据最优偏移量对回波包络进行微调，将每次回波进行包络精对齐后，最终得到平动补偿后的目标回波；⑷算法选择4a)采用距离-多普勒算法对平动补偿后的回波成像，得到目标的距离-多普勒图像；4b)若距离-多普勒图像的方位维没有混叠，则选择CIRT算法，执行步骤(5)；4c)若距离-多普勒图像出现方位维混叠，则选择RIRT算法；(5) 二维切片成像5a)建立自旋目标的柱面坐标系，沿目标自转轴方向将自旋目标分成多个二维切片；5b)将二维切片按照步骤(4)选择的算法进行成像，得到二维切片图像；(6)重构回波并更新6a)对于所有的二维切片图像，提取图像域中最大峰值对应的三维坐标，记录为目标散射点的位置坐标，并在距离-慢时间域构造散射点回波；6b)运用最小范数准则对步骤6a)中所构造的散射点回波的幅度进行估计；6c)将记录的散射点的位置坐标与估计的回波幅度组合构成散射点回波；6d)从步骤3d)得到的目标回波中减去散射点回波，更新目标回波；(7)判断目标回波能量是否低于能量门限7a)若目标回波能量低于门限，则执行步骤⑶；7b)若目标回波能量高于门限，则执行步骤(5)；(8)操作结束，根据6a)中记录的目标散射点位置重构目标三维图像。与现有技术相比，本发明具有以下优点第一，本发明通过采用回波相干化处理，克服了现有技术中广义Radon变换方法受目标高速运动影响，存在距离像偏移和展宽，以及方位像散焦的不足，具有消除高速运动目标对图像的影响，获得聚焦良好图像的优点。第二，本发明通过采用以图像熵最小为准则，建立有效的搜索算法对目标回波的平动量进行精细补偿。克服了现有技术中广义Radon变换方法无法对目标平动进行准确补偿的不足，具有获得平动准确补偿后的自旋目标回波的优点。第三，本发明通过采用针对不同脉冲重复频率条件下的算法选择，克服了现有技术中广义Radon变换方法成像旁瓣较高、分辨率低的不足，具有提高不同脉冲重复频率条件下目标的成像率，并使图像最高分辨率达到波长量级，成像质量高的优点。


图1为本发明的流程图；图2为本发明的仿真图。
具体实施例方式下面结合附图1，对本发明具体实施方式
作进一步的详细描述步骤1，获取目标的ISAR回波，雷达以脉冲重复频率发射并接收脉冲，得到以距离为行向量以方位为列向量的ISAR回波； 步骤2，对目标的ISAR回波进行相干化处理2a)为消除速度所引起的一维距离像平移和展宽，在距离频域对回波乘以相位补偿因子
权利要求
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)得到的目标回波中减去散射点回波，更新目标回波；(7)判断目标回波能量是否低于能量门限7a)若目标回波能量低于门限，则执行步骤(8)； 7b)若目标回波能量高于门限，则执行步骤(5)；(8)操作结束，根据6a)中记录的目标散射点位置重构目标三维图像。
2.根据权利要求1所述的基于复数后向投影的自旋目标三维成像方法，其特征在于， 步骤4b)中所述的CIRT算法公式为
3.根据权利要求1所述的基于复数后向投影的自旋目标三维成像方法，其特征在于， 步骤如)中所述的RIRT算法公式为第一步，按照下式对平动补偿后的回波实包络沿距离域进行傅里叶变换S ( ξ , tm) =S ι S (r, tm) I exp (_j ξ r) dr其中，s(〖，tm)为平动补偿后的回波实包络沿距离域进行傅里叶变换后的结果，ξ为距离频域，支撑区为[-η，n]，tm为慢时间，|s(r，tm) I为自旋目标回波的实包络，r为距离，j为虚数单位；第二步，将s( ξ，tm)按照下式进行一维滤波与逆傅里叶变换
4.根据权利要求1所述的基于复数后向投影的自旋目标三维成像方法，其特征在于， 步骤(7)中所述的能量门限为回波总能量的5%。
全文摘要
本发明公开了一种基于复数后向投影的自旋目标三维成像方法，其步骤包括(1)雷达录取ISAR回波；(2)回波相干化；(3)平动补偿；(4)算法选择；(5)二维切片成像；(6)重构回波并更新；(7)判断目标回波能量是否低于能量门限；(8)操作结束，重构目标三维图像。本发明通过回波相干化处理、平动精细补偿以及针对不同脉冲重复频率进行算法选择，克服了广义Radon变换方法无法对目标回波的平动精确补偿，受目标高速运动影响，对应图像旁瓣较高、分辨率低的不足，具有平动补偿精确，消除目标高速运动对图像的影响，提高不同脉冲重复频率条件下的成像率，对应图像最高分辨率达到波长量级，成像质量高的优点。
文档编号G01S13/90GK102540189SQ201210001679
公开日2012年7月4日 申请日期2012年1月4日 优先权日2012年1月4日
发明者保铮, 刘妍, 周峰, 白雪茹 申请人:西安电子科技大学