基于OMP的SAR目标图像几何尺寸估计方法技术

本发明专利技术公开了一种基于正交匹配追踪算法OMP的合成孔径雷达SAR目标图像几何尺寸估计方法，其实现步骤为：(1)选取图像切片矩阵；(2)获得目标的原子矩阵和散射中心幅度向量；(3)获得重构目标图像；(4)将重构目标图像投影到坐标系生成目标区域；(5)获取坐标向量；(6)对坐标向量进行旋转；(7)获得最小外接矩形；(8)将最小外接矩形的尺寸值作为目标几何尺寸的估计值。本发明专利技术对静止雷达观测目标普遍适用，同时本发明专利技术消除了由于噪声而引起的负面影响，具体良好的噪声稳健性，提高了合成孔径雷达SAR目标几何尺寸的估计精度。

【技术实现步骤摘要】
基于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表示超分辨基本文档来自技高网...

【技术保护点】
1.一种基于正交匹配追踪算法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图像目标几何尺寸的估计结果。...

【技术特征摘要】
1.一种基于正交匹配追踪算法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)获得目标几何...

【专利技术属性】
技术研发人员：杜兰，崔曼，杨栋文，王家东，
申请(专利权)人：西安电子科技大学，西安中电科西电科大雷达技术协同创新研究院有限公司，
类型：发明
国别省市：陕西,61