基于OMP的SAR目标图像几何尺寸估计方法与流程

本发明属于通信技术领域，更进一步涉及雷达目标特征提取技术领域中的一种基于正交匹配追踪算法omp的合成孔径雷达sar(syntheticapertureradar)目标图像几何尺寸估计方法。本发明可用于提取雷达接收到的单幅图像中某一静止目标的外层轮廓，利用提取的轮廓估计该静止目标的几何尺寸。

背景技术：

目标的几何尺寸特征是雷达目标直观的物理特征，对于高分辨雷达图像目标切片，估计其中静止目标的几何尺寸在sar图像目标鉴别和分类过程中具有重要的应用意义。对于高分辨雷达图像目标的尺寸估计问题，实际得到的高分辨雷达图像数据一般都受到目标姿态敏感性和噪声的影响，如何能从实际雷达目标图像中获得准确的目标几何尺寸信息是个很重要的研究课题。目前已有的高分辨雷达目标物理尺寸估计方法有基于目标运动模型对空间锥体目标估计物理几何尺寸的。但这些估计方法需要知道目标的运动模型，不适用于静止目标。基于正交匹配追踪算法omp的雷达目标图像几何尺寸估计方法不需要知道目标的运动模型就能获得目标的物理尺寸信息，并且基于正交匹配追踪算法omp的雷达目标尺寸估计方法能够实现目标的去噪处理，消除噪声对目标尺寸估计的影响。

西安电子科技大学在其申请的专利文献“基于激光一维距离像反演动态锥形目标姿态角和几何尺寸”(专利申请号cn201410188024.4，公开号cn103954969a)中提出了一种反演目标几何尺寸的方法。该方法的实现过程是：首先建立目标的激光一维距离像模型，通过目标一维距离像的成像特征，提取其峰值坐标，作为特征矢量；设定遗传参数和遗传计算变量范围，针对某一时刻的目标进行距离像采样；根据一维距离像峰值大小、峰值点在入射方向的投影长度、以及雷达阵列的角度关系进行遗传计算然后取均值，根据目标的运动模型反演得到目标的几何尺寸。该方法通过函数反演，能够精确得到目标的几何形态特征。但是，该方法存在的不足之处是需要知道目标的运动模型，不适用于静止目标，所以该方法不能普遍适用。

中国科学院计算技术研究所在其申请的专利文献“一种sar图像舰船目标几何特征并行提取方法与装置”(专利申请号cn201710182152.1，申请公开号cn107103283a)中提出一种sar图像目标几何特征提取方法。该方法的实现过程是：生成像素在roi内的目标原始坐标，根据设定的阈值进行目标筛选；根据目标坐标偏移量和目标原始坐标，计算生成目标偏移后坐标；根据预先设定好的旋转步进值和目标偏移后坐标，计算相应旋转角度下目标在x轴上的投影长度；根据该投影长度，计算相应旋转角度下y轴上的投影长度，得到目标尺寸特征。该方法存在的不足之处是，由于该方法在进行目标筛选时，直接根据设定的阈值将目标和背景分开，无法消除噪声对目标几何尺寸估计的影响，不具有噪声稳健性。

技术实现要素：

本发明的目的在于针对上述现有技术存在的不足，提出一种基于正交匹配追踪算法omp的合成孔径雷达sar目标图像几何尺寸估计方法。

实现本发明目的的思路是，利用正交匹配追踪算法omp在合成孔径雷达sar目标图像上获取观测目标的重构目标二值图像，根据外接矩形方法，在重构目标二值图像上估计出观测目标的几何尺寸，消除了噪声在尺寸估计中的负面影响，提高了雷达目标尺寸估计方法的噪声稳健性。同时本发明在图像域对目标进行尺寸估计，不需要知道目标的运动模型，具有普遍适用的优点。

本发明的具体步骤如下：

(1)选取一个只包含雷达观测目标的合成孔径雷达sar图像中的一个切片矩阵；

(2)利用正交匹配追踪算法omp，生成目标的原子矩阵和散射中心幅度向量：

(2a)将所选切片矩阵中的每列依次首尾相连，组成一个切片列向量；

(2b)利用正交傅里叶基公式，计算超分辨基矩阵中的每个元素值；

(2c)将切片列向量作为初始残差信号向量，用超分辨基矩阵中每列向量分别与初始残差信号向量做内积操作，提取所有内积结果中最大值对应的超分辨基矩阵的一列元素，加入一个空矩阵中，组成当前原子矩阵；

(2d)利用残差信号更新公式，计算更新后的残差信号向量，用超分辨基矩阵中每列向量分别与更新后的残差信号向量做内积操作，将所有内积结果中最大值对应的超分辨基矩阵的一列元素，加入当前原子矩阵中；

(2e)将更新后的残差信号向量做转置操作后与自身相乘，得到更新后的残差信号的功率值；

(2f)判断更新后的残差信号的功率值是否小于终止门限值，若是，则得最终的原子矩阵后执行步骤(2g)，否则，执行步骤(2d)；

(2g)将最终原子矩阵取逆操作后与切片列向量的快速傅立叶变换相乘，得到散射中心幅度向量；

(3)获得重构目标二值图像：

将散射中心幅度向量与最终原子矩阵相乘，得到一个重构目标矩阵，对重构目标矩阵做二值化操作，得到重构目标二值图像；

(4)生成目标区域：

将重构目标二值图像投影到一个二维平面坐标系中，将重构目标二值图像中所有值为1的像素点，组成目标区域；

(5)获取坐标向量：

取出目标区域最外层轮廓上所有像素点在x轴上和y轴上对应的坐标值，分别沿逆时针方向依次排列，组成横轴和纵轴的坐标向量；

(6)获得外接矩形：

(6a)利用外接矩形方法，获得二维平面坐标系中的一个矩形；

(6b)利用旋转公式，对横轴和纵轴坐标向量依次从1°到90°范围内每隔1°进行一次旋转操作；

(6c)对每次旋转后的横轴和纵轴坐标向量通过外接矩形方法，获得外接矩形；

(7)外接矩形的周长：

(7a)将旋转后的横轴坐标向量最大值与最小值之差作为外接矩形的长度值；

(7b)将旋转后的纵轴坐标向量最大值与最小值之差作为外接矩形的宽度值；

(7c)将长度值与宽度值相加后乘以2作为外接矩形的周长值；

(8)获得目标几何尺寸的估计结果：

将所有外接矩形中周长值最小的外接矩形对应的长度值与宽度值，作为合成孔径雷达sar图像目标几何尺寸的估计结果。

本发明与现有的技术相比具有以下优点：

第一，由于本发明采用正交匹配追踪算法omp得到雷达观测目标的重构目标图像，去除了合成孔径雷达sar图像中的噪声信号，克服了现有技术sar图像目标几何特征并行提取方法噪声影响目标几何尺寸估计结果的问题，使得本发明具噪声稳健性强的优点。

第二，由于本发明采用外接矩形方法，获得雷达观测目标的外接矩形，利用雷达目标的外接矩形估计雷达目标的几何尺寸，克服了现有技术采用一维距离像反演几何尺寸估计方法需要已知目标运动模型的问题，使得本发明具有对静止雷达目标普遍适用的优点。

附图说明

图1是本发明的流程图；

图2是本发明的仿真图。

具体实施方式

下面结合附图对本发明作进一步的描述。

参照附图1，对本发明的实现步骤作进一步的描述。

步骤1，选取一个只包含雷达观测目标的合成孔径雷达sar图像中的一个切片矩阵。

步骤2，利用正交匹配追踪算法omp，生成目标的原子矩阵和散射中心幅度向量。

所述正交匹配追踪算法omp具体步骤如下：

第1步，将所选切片矩阵中的每列依次首尾相连，组成一个切片列向量。

第2步，利用正交傅里叶基公式，计算超分辨基矩阵中的每个元素值。

所述正交傅里叶基公式如下：

其中，d表示超分辨基矩阵的第l行第m列的元素，e表示以自然常数e为底的指数操作，j表示虚数单位符号，π表示圆周率，l表示超分辨基矩阵的行序号，l＝1,…,n，m表示超分辨基矩阵的列序号，m＝1,…,m，m表示超分辨基矩阵的总列数。

第3步，将切片列向量作为初始残差信号向量，用超分辨基矩阵中每列向量分别与初始残差信号向量做内积操作，提取所有内积结果中最大值对应的超分辨基矩阵的一列元素，加入一个空矩阵中，组成当前原子矩阵。

第4步，利用残差信号更新公式，计算更新后的残差信号向量，用超分辨基矩阵中每列向量分别与更新后的残差信号向量做内积操作，将所有内积结果中最大值对应的超分辨基矩阵的一列元素，加入当前原子矩阵中。

所述残差更新公式如下：

r＝fft(x)-ef⁺fft(x)

其中，r表示更新后的残差信号向量，fft(·)表示做快速傅立叶变换操作，x表示切片列向量，e表示超分辨基矩阵，f表示当前原子矩阵，+表示穆尔-彭罗斯moore-penrose逆操作。

第5步，将更新后的残差信号向量做转置操作后与自身相乘，得到更新后的残差信号的功率值。

第6步，判断更新后的残差信号的功率值是否小于终止门限值，若是，则得最终的原子矩阵后执行第7步，否则，执行第4步。

所述终止门限的计算公式如下：

其中，th表示终止门限值，n表示切片列向量的长度，sn表示噪声向量，h表示共轭转置操作，l表示对合成孔径雷达sar图像切片矩阵的第一列做快速傅里叶变换后得到噪声向量的长度。

第7步，将最终原子矩阵取逆操作后与切片列向量的快速傅立叶变换相乘，得到散射中心幅度向量。

步骤3，获得重构目标二值图像。

将散射中心幅度向量与最终原子矩阵相乘，得到一个重构目标矩阵，对重构目标矩阵做二值化操作，得到重构目标二值图像。

步骤4，生成目标区域。

将重构目标二值图像投影到一个二维平面坐标系中，将重构目标二值图像中所有值为1的像素点，组成目标区域。

步骤5，获取坐标向量。

取出目标区域最外层轮廓上所有像素点在x轴上和y轴上对应的坐标值，分别沿逆时针方向依次排列，组成横轴和纵轴的坐标向量。

步骤6，获得外接矩形。

利用外接矩形方法，获得二维平面坐标系中的一个矩形。

利用下式，对横轴和纵轴坐标向量依次从1°到90°范围内每隔1°进行一次旋转操作：

其中，表示旋转后的坐标向量，x表示原始坐标向量，cos表示余弦操作，sin表示正弦操作，β表示旋转角度。

对每次旋转后的横轴和纵轴坐标向量通过外接矩形方法，获得外接矩形。

所述外接矩形方法的步骤如下：

第1步，分别在x轴和y轴上找到横轴和纵轴坐标向量中的最小值，构成二维平面坐标系中一个起始点；

第2步，在x轴上找到横轴坐标向量的最大值，在y轴上找到纵轴坐标向量的最小值，构成二维平面坐标系中一个中间点；

第3步，分别在x轴和y轴上找到横轴和纵轴坐标向量的最大值，构成二维平面坐标系中一个递进点；

第4步，在x轴上找到横轴坐标向量的最小值，在y轴上找到纵轴坐标向量的最大值，构成二维平面坐标系中一个终止点；

第5步，将起始点，中间点，递进点，终止点依次连接成一个矩形。

步骤7，计算外接矩形的周长。

将旋转后的横轴坐标向量最大值与最小值之差作为外接矩形的长度值。

将旋转后的纵轴坐标向量最大值与最小值之差作为外接矩形的宽度值。

将长度值与宽度值相加后乘以2作为外接矩形的周长值。

步骤8，获得目标几何尺寸的估计结果。

将所有外接矩形中周长值最小的外接矩形对应的长度值与宽度值，作为合成孔径雷达sar图像目标几何尺寸的估计结果。

下面结合仿真实验对本发明的效果做进一步的说明。

1.仿真实验条件：

本发明的仿真实验的硬件测试平台是：处理器为intelcorei7cpu，主频为3.40ghz，内存8gb，软件平台为：windows7操作系统和matlabr2017a。

2.仿真实验内容：

本发明的仿真实验是利用宽带雷达采集的坦克、装甲运输车和步兵战车等10类不同地面目标的回波信号，采用本发明的方法对每类目标的雷达图像进行尺寸估计。10类不同地面目标的回波信号依次是步兵战车bmp2、装甲运输车btr70、装甲运输车btr60、自行榴弹炮2s1、装甲侦察车brdm2、坦克t62、自行高炮zsu234、推土机d7、货运卡车zil131、坦克t72。

图2是本发明的仿真图，图2(a)是利用本发明方法对10类地面目标的长度估计的相对误差结果图，图2(a)中的横坐标表示10类地面目标，纵坐标表示相对误差结果，图2(a)中以实线表示的曲线代表每个目标的长度相对误差结果，以虚线表示的直线代表所有目标的平均长度相对误差。图2(b)是利用本发明方法对上述10类地面目标中除推土机d7之外的9类地面目标的宽度估计相对误差结果图，图2(b)的横坐标表示9类地面目标，纵坐标表示相对误差结果，图中以实线表示的曲线代表每个目标的宽度相对误差结果，以虚线表示的直线代表所有目标的平均宽度相对误差。

3.仿真结果分析：

为了定量评价本发明方法的效果，按照下式，计算本发明仿真实验对不同目标的长度和宽度尺寸的相对误差结果：

其中，error表示相对误差结果，m表示按照本发明估计出来的尺寸值，n表示目标实际的尺寸值，error值越小，说明对雷达目标的尺寸估计结果越精确。

由图2(a)和图2(b)可见，每幅图中的实线对应的值大部分小于5％，虚线对应的值小于5％，表明本发明方法得到的雷达观测目标的尺寸估计结果较为准确，说明本发明方法确实能够有效的对雷达观测目标的几何尺寸进行估计。