基于clean算法散射中心提取的isar目标轮廓提取方法
【技术领域】
[0001] 本发明属于雷达成像技术领域，涉及基于CLEAN算法散射中心提取的ISAR目标轮 廓提取方法。
【背景技术】
[0002] 逆合成孔径雷达成像是一种获取运动目标高分辨图像的有效手段，特别是应用在 在雷达目标检测识别与分类中。逆合成孔径雷达（inversesyntheticapertureradar， ISAR)成像在目标分类识别和检测中都有重要应用价值。在实际雷达成像中，目标既有平动 又有转动，平移运动使目标上散射点回波产生相同的多普勒频率偏移，对成像没有帮助必 须补偿掉。而转动则是实现成像的基础，它使得目标上不同方位的散射点产生不同的多普 勒分量，因而可以得到方位的多普勒变化信息。补偿掉平移运动之后，目标的运动近似为围 绕雷达转动，目标可以看作是转台目标。即可获得其聚焦的二维ISAR像，用于后续的轮廓 特征提取。
[0003] 作为一种重要的自动目标识别途径，ISAR像的轮廓提取是ISAR像目标特性反演 和识别的基础。很多现有文献对ISAR目标轮廓提取进行了深入的研究。然而这些轮廓提 取方法大都是基于形态学运算，需要将图像转换成二值图，很少有文献从ISAR回波数据本 身入手。雷达目标的光学图像能显示出比较完整的目标轮廓形状，然而雷达图像则通常表 现为目标上稀疏的散射中心分布，这些散射中心对应了目标的棱边、尖顶等的相对位置，反 映了目标的形状结构。
【发明内容】
[0004] 本发明的目的在于提供一种基于CLEAN算法散射中心提取的ISAR目标轮廓提取 方法，解决了现有的轮廓提取方法大都基于形态学运算，其准确性受到限制的问题。
[0005] 本发明所采用的技术方案是，基于CLEAN算法散射中心提取的ISAR目标轮廓提取 方法，按照以下步骤进行：
[0006] 步骤1:利用二维CLEAN算法提取ISAR目标的散射中心位置信息，根据提取的N 个散射强点对ISAR图像进行重构得到最终的分辨图像，
[0007]
[0008] An代表位于坐标（xn，yn)处的散射强度最高的点的幅度，h(x，y)为点扩散函数，
[0009] 步骤2 :利用a-shape算法得到落在目标边缘处的散射中心位置；
[0010] 步骤3 :通过橡皮筋算法将这些散射中心连成闭合曲线，并对结果进行平滑处理 获得平滑的目标轮廓线。
[0011] 进一步的，步骤2的具体过程是：基于a-Shape算法得到的凸包获取落在目标边 缘处的所有散点，点集S内，过任意两点ΡρΡ2绘制半径为α的圆，如果这个圆内没有其他 点，则认为点Pi，匕是边界点，其连线是边界线段，若已知Pi，Ρ2的坐标（Xdyi)，（x2,y2)，即 可得到圆心P3的坐标（x3,y3)，
[0012]
[0013]
[0014] Η是中间量。
[0015] 进一步的，步骤3具体包括：通过a-shape获取位于边界处的所有点，此处将边 界点所在点集定义为Ps，找到所有凸包，定义为Cs，并利用线段，定义为Ls，按照顺时针顺 序连接所有凸包，剩余的未包含的点集定义为Rs，根据以上的定义，橡皮筋算法包括以下步 骤：
[0016] 步骤a.将Rs中所有的点按照X坐标递增的方式排列；
[0017] 步骤b.搜索Rs中第一个点，见图2(b)中的点3,并从Ls中找出与该点最接近的 边，此处为c边；
[0018] 步骤c.如果三角形内角A、B、C满足以下条件：A<90°，B<90°andC>Θ， 用三角形另外两条边a和b替换第三条边c，之后将点3从Rs移到Cs中，并且将其按照顺 时针方向插入点1和2中间，其中Θ为角度门限以确保轮廓的光滑度；
[0019] 步骤d.重复步骤a到c直至所有Rs中的点都处理完毕；如果Rs= {0}，转到步 骤f，否则转到步骤e;
[0020] 步骤e.利用给定的角度门限0，RS中可能还有剩余的点，接下来以步长d0减小 Θ并重复步骤a到d直至Rs= {0};
[0021] 步骤f.最后利用平滑的曲线代替折线获得平滑的闭合轮廓模型。
[0022] 本发明的有益效果是将二维CLEAN高分辨技术应用于ISAR像轮廓提取当中，通过 本文提出的橡皮筋算法将位于目标边缘处的散射中心用闭合的平滑曲线连接以获得ISAR 目标轮廓提取当中，便于后续的ISAR目标分类识别。
【附图说明】
[0023] 图1是a-shape算法流程图。
[0024] 其中，图1 (a)是原始散列点集，图1 (b)是a-shape结果。
[0025] 图2是目标轮廓提取示意图。
[0026] 其中，图2(a)α端点提取结果，图2(b)用三角形的两条边替代第三条边，图2(c) 线段连接结果，图2(d)最终平滑轮廓。
[0027] 图3是实测数据ISAR成像结果图。
[0028] 其中，图3 (a)雅克-42几何模型，图3 (b)ISAR成像结果。
[0029] 图4是散射中心位置提取结果图。
[0030] 其中，图4(a)实际散射中心位置，（b) 10点圆形插值后的结果。
[0031] 图5是实测数据ISAR目标轮廓提取示意图。
[0032] 其中，图5(a)a-shape结果，图5(b)橡皮筋算法结果，图5(c)线段连接结果，图 5(d)最终平滑轮廓。
[0033] 图6是本文方法与传统snake方法对比结果图。
[0034] 其中，图6 (a)snake结果，（b)本文结果。
【具体实施方式】
[0035] 下面结合附图和【具体实施方式】对本发明进行详细说明。
[0036] 本文提出一种基于CLEAN和橡皮筋的目标轮廓提取的新方法，通过平动补偿获得 聚焦的ISAR像。提取目标的二维散射中心位置。散射中心提取方法在很多文献中均有所提 及。本文利用CLEAN算法从二维ISAR中获取目标的散射中心位置。利用a-shape方法， 获取位于边界处的散射中心的位置，并连成闭合曲线。此处利用本文提出的一种名为橡皮 筋的新方法，从散点的凸包入手，通过拉动橡皮筋的方式获得闭合轮廓线。最后，利用实测 数据验证本文的可行性及准确性。
[0037] 具体来说，本发明按照以下3个步骤进行：
[0038] 步骤1 :在光学区，雷达目标的频率响应特性可以表示为各个散射中心的频率响 应之和。故此本文通过提取散射中心位置代替原始ISAR图像。散射中心提取方法很多，本 文采用的二维CLEAN算法就是为了解决旁瓣抑制与分辨率之间的矛盾而提出的。CLEAN算 法稳定性高，算法简单，被广泛应用于雷达成像中。
[0039] 令In代表第η步的二维残余图像，An代表位于坐标（Xn，yn)处的散射强度最高的 点的幅度，h(x，y)为点扩散函数。则残余的图像可以表示为
[0040] Ιη+ι=In-An ·h(x-xn,y-yn) (1)
[0041] CLEAN技术的基本原理为：依次提取散射强度最大的散射点并将其从图像域中抠 除，从而生成新的图像，重复上述过程直到终止条件满足。实际数据处理中，算法是对二维 数据进行的搜索迭代过程，需要多次将强点对应的响应提出，需要预先设定某一特定阈值 或循环步数作为退出循环的条件，当残余复图像幅度小于阈值时终止计算。令In代表第η 步的二维残余图像，总代表位于（xn，yn)处的散射强度最高的点的幅度，则残余的图像可以 表示为
[0042] In=Iηι~Αη ·h(x-xn,y-yn) (2)
[0043] 即第n步的二维残余图像等于第n-1步残余图像减去最强的散射点。根据提取的 N个散射强点对ISAR图像进行重构得到最终的高分辨图像。
[0044]
[0045] 步骤2 :基于2D(二维）U-ESPRIT的散射中心提取；
[0046] 本苄基于a-shape算法得到的凸包获取落在目标边缘处的所有散点。点集S内， 过任意两点Pi，P2绘制半径为α的圆，如果这个圆内没有其他