发明名称 SAR图像目标识别方法
摘要 本发明公开一种SAR图像目标识别方法,主要解决现有技术的样本特征维数过大,计算复杂,时间消耗长的问题。其实现步骤为:(1)对SAR图像初始化,得到训练样本矩阵和测试样本矩阵;(2)利用Johnson–Lindenstrauss推论和有限等距RIP条件,构造稀疏测量矩阵;(3)根据稀疏测量矩阵,对训练样本矩阵和测试样本矩阵进行降维并归一化,获得归一化后的样本矩阵;(4)根据归一化后的样本矩阵,利用最近邻分类器,得到测试样本的类别标号。本发明与现有技术相比降低了样本的特征维数和计算复杂度,提高了SAR目标识别的精度和运算速度,可用于图像处理。
申请公布号 CN103426001B 申请公布日期 2016.08.10
申请号 CN201310391126.1 申请日期 2013.08.30
申请人 西安电子科技大学 发明人 于昕;焦李成;韩文婷;马文萍;马晶晶;侯彪
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种SAR图像自动目标识别方法,包括如下步骤:(1)输入SAR图像,并对其进行Haar小波变换,得到训练样本矩阵X∈R<sup>m×N</sup>和测试样本矩阵X′∈R<sup>m×N′</sup>,其中,R表示实数集合,m表示原样本的特征维数,N表示训练样本的个数,N′表示测试样本的个数;(2)根据Johnson–Lindenstrauss推论和RIP条件得出一个稀疏测量矩阵W∈R<sup>n×m</sup>,其矩阵元素为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>a</mi><mi>b</mi></mrow></msub><mo>=</mo><msqrt><mi>s</mi></msqrt><mo>&times;</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>p</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>s</mi></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>p</mi><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mi>s</mi></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>p</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>s</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000978790330000011.GIF" wi="678" he="402" /></maths>其中,w<sub>ab</sub>表示稀疏测量矩阵W的第a行第b列的元素,a∈{1,2,…,n},n为稀疏测量矩阵的维数,b∈{1,2,…,m},s为稀疏测量矩阵的量化参数,p为稀疏测量矩阵中元素出现的概率;(3)将步骤(2)中得到的稀疏测量矩阵W,分别与训练样本矩阵X和测试样本矩阵X′相乘进行降维,得到降维后的训练样本矩阵Z<sub>train</sub>∈R<sup>n×N</sup>和测试样本矩阵Z<sub>test</sub>∈R<sup>n×N′</sup>,并对该降维后的训练样本矩阵和测试样本矩阵进行归一化,得到归一化的训练样本矩阵<img file="FDA0000978790330000012.GIF" wi="235" he="70" />和测试样本矩阵<img file="FDA0000978790330000013.GIF" wi="251" he="71" />(4)对N个训练样本设置类别标号Y∈R<sup>1×N</sup>;(5)计算归一化后的测试样本矩阵<img file="FDA0000978790330000014.GIF" wi="84" he="71" />和N个已知类别的训练样本之间的距离:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>d</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>x</mi><mi>l</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi></mrow></munder><mo>|</mo><mo>|</mo><msubsup><mi>x</mi><mi>l</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>j</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>|</mo><mo>|</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>c</mi><mo>,</mo></mrow></mtd><mtd><mrow><mi>l</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>N</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000978790330000015.GIF" wi="1230" he="103" /></maths>其中,x′<sub>l</sub>是归一化后测试样本矩阵<img file="FDA0000978790330000016.GIF" wi="78" he="71" />的第l个样本,<img file="FDA0000978790330000017.GIF" wi="62" he="78" />是训练样本矩阵<img file="FDA0000978790330000018.GIF" wi="91" he="70" />中的第j个样本,属于类别标号Y中的第i类,c是训练样本的类别数;(6)将步骤(5)中得到的距离,带入下式中,得到测试样本的最终类别标号:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><mi>arg</mi><mo>&lsqb;</mo><munder><mrow><msub><mi>mind</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>x</mi><mi>l</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>c</mi></mrow></munder><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000978790330000021.GIF" wi="413" he="135" /></maths>其中,k是测试样本x′<sub>l</sub>的最终类别标号。
地址 710071 陕西省西安市太白南路2号