发明名称 基于几何部件及径向梯度角的遥感检测飞机的方法
摘要 本发明公布了一种基于几何部件及径向梯度角的遥感检测飞机的方法。该方法首先根据飞机特有的几何轮廓建立几何原子库,结合稀疏表示原理,将已知飞机轮廓稀疏表示成几何原子的组合,并以这些几何原子为部件,构建飞机星型部件模型;然后采用滑动窗口法对测试图像进行扫描,通过径向‑梯度角投票算法估算出每个图像块的主方向,用于调整部件模型方向,并用调整后的模型计算出每个位置的总相似度分布图,最后通过最大后验概率推理,判断检测区域是否存在飞机,并进一步得出飞机的位置。该方法利用稀疏表示原理自适应选取几何部件,并采用轮廓为引导的径向梯度角投票算法估计主方向,实现了复杂背景下不同方向的飞机遥感检测。其算法简单,检测速度快,鲁棒性更好,检测结果更准确。
申请公布号 CN103714344B 申请公布日期 2017.02.08
申请号 CN201310703176.9 申请日期 2013.12.19
申请人 西南交通大学 发明人 林煜东;和红杰;尹忠科;陈帆
分类号 G06K9/62(2006.01)I;G06K9/32(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 成都博通专利事务所 51208 代理人 陈树明
主权项 一种基于几何部件及径向梯度角的遥感检测飞机的方法,其步骤如下:A、飞机星型部件模型的建立:A1、将已知的遥感图像中的已知飞机图像校正为机头竖直向上,得到校正图像;A2、将校正图像中飞机部分的像素值置为1,背景部分像素值置为0,得到已知飞机的二值轮廓P<sub>B</sub>;A3、对已知飞机的二值轮廓P<sub>B</sub>进行稀疏分解,得到最优原子图像集<img file="FDA0001121497710000011.GIF" wi="606" he="70" />N<sup>Opt</sup>为最优原子图像个数,i为最优原子图像<img file="FDA0001121497710000012.GIF" wi="83" he="63" />的索引编号,<img file="FDA0001121497710000013.GIF" wi="85" he="63" />为第i个最优原子图像;A4、取第i个最优原子图像<img file="FDA0001121497710000014.GIF" wi="83" he="70" />里的非零区域得到第i个几何图形<img file="FDA0001121497710000015.GIF" wi="107" he="70" />以第i个几何图形<img file="FDA0001121497710000016.GIF" wi="81" he="70" />作为第i个部件Part<sub>i</sub>,以已知飞机的二值轮廓P<sub>B</sub>为根部件Root,进行星型模型关联;关联时,第i个部件Part<sub>i</sub>的形状为第i个几何图形<img file="FDA0001121497710000017.GIF" wi="83" he="62" />的形状<img file="FDA0001121497710000018.GIF" wi="179" he="63" />第i个部件Part<sub>i</sub>到根部件Root中心(x<sup>R</sup>,y<sup>R</sup>)的相对坐标位置为第i个几何图形<img file="FDA0001121497710000019.GIF" wi="83" he="62" />中心到原子图像<img file="FDA00011214977100000110.GIF" wi="89" he="63" />中心的相对坐标<img file="FDA00011214977100000111.GIF" wi="332" he="63" />从而第i个部件Part<sub>i</sub>的主方向为第i个几何图形<img file="FDA00011214977100000112.GIF" wi="84" he="70" />的主方向<img file="FDA00011214977100000113.GIF" wi="126" he="61" />其中,x<sup>R</sup>,y<sup>R</sup>分别为根部件Root中心在校正图像中的横坐标和纵坐标;A5、重复A4,将所有的几何图形<img file="FDA00011214977100000114.GIF" wi="83" he="62" />均作为部件Part<sub>i</sub>与根部件Root进行关联,建立得到飞机星型部件模型;B、遥感图像中的飞机检测B1、通过卫星遥感得到检测区域的遥感图像P,取大小与已知飞机的二值轮廓P<sub>B</sub>相同的检测窗口,在遥感图像P中进行逐点扫描,获取子图像块P<sub>sub</sub>;B2、采用径向‑梯度角投票算法计算出子图像块P<sub>sub</sub>的主方向ξ<sub>sub</sub>;B3、将飞机星型模型旋转至机头方向与子图像块P<sub>sub</sub>的主方向ξ<sub>sub</sub>一致,计算旋转后的飞机星型模型中各部件Part<sub>i</sub>的主方向<img file="FDA00011214977100000115.GIF" wi="131" he="63" />及与根部件Root中心(x<sup>R</sup>,y<sup>R</sup>)的相对坐标<img file="FDA00011214977100000116.GIF" wi="334" he="63" /><img file="FDA00011214977100000117.GIF" wi="373" he="63" /><maths num="0001"><math><![CDATA[<mrow><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><msubsup><mi>Px</mi><mi>R</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>Py</mi><mi>R</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><msub><mi>cos&xi;</mi><mrow><mi>s</mi><mi>u</mi><mi>b</mi></mrow></msub></mrow></mtd><mtd><mrow><msub><mi>sin&xi;</mi><mrow><mi>s</mi><mi>u</mi><mi>b</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>sin&xi;</mi><mrow><mi>s</mi><mi>u</mi><mi>b</mi></mrow></msub></mrow></mtd><mtd><mrow><msub><mi>cos&xi;</mi><mrow><mi>s</mi><mi>u</mi><mi>b</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><msup><mi>Px</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>Py</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msup></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00011214977100000118.GIF" wi="822" he="151" /></maths>B4、在飞机星型模型中根据第i个部件Part<sub>i</sub>旋转后的主方向<img file="FDA00011214977100000119.GIF" wi="137" he="71" />第i个部件Part<sub>i</sub>旋转后与根部件Root中心(x<sup>R</sup>,y<sup>R</sup>)的相对坐标<img file="FDA0001121497710000021.GIF" wi="332" he="63" />及第i个部件的形状<img file="FDA0001121497710000022.GIF" wi="179" he="63" />采用几何方法计算得到飞机星型模型旋转后第i个部件所覆盖的第i个区域<img file="FDA0001121497710000023.GIF" wi="179" he="67" />B5、将飞机星型模型中第i个区域<img file="FDA0001121497710000024.GIF" wi="147" he="64" />内的像素置为1,其它像素置为0,生成第i个二值部件图<img file="FDA0001121497710000025.GIF" wi="277" he="63" />即:<maths num="0002"><math><![CDATA[<mrow><msup><mi>Pic</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>&Element;</mo><msubsup><mi>Reg</mi><mi>R</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001121497710000026.GIF" wi="691" he="156" /></maths>式中,x,y分别为二值部件图<img file="FDA0001121497710000027.GIF" wi="254" he="70" />在校正图像中的横坐标和纵坐标;B6、计算第i个二值部件图<img file="FDA0001121497710000028.GIF" wi="256" he="64" />的方向梯度直方图,生成第i个部件滤波器F<sub>i</sub>:<maths num="0003"><math><![CDATA[<mrow><msub><mi>F</mi><mi>i</mi></msub><mo>=</mo><mi>H</mi><mi>O</mi><mi>G</mi><mrow><mo>(</mo><msup><mi>Pic</mi><mrow><msub><mi>Part</mi><mi>i</mi></msub></mrow></msup><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001121497710000029.GIF" wi="501" he="83" /></maths>式中,<img file="FDA00011214977100000210.GIF" wi="181" he="54" />为方向梯度直方图特征提取算子;计算遥感图像P的子图像块P<sub>sub</sub>的方向梯度直方图H<sub>sub</sub>:H<sub>sub</sub>=HOG(P<sub>sub</sub>)B7、计算子图像块P<sub>sub</sub>与旋转后的飞机星型模型中的所有部件的总相似度Score:<maths num="0004"><math><![CDATA[<mrow><mi>S</mi><mi>c</mi><mi>o</mi><mi>e</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mrow><mi>O</mi><mi>p</mi><mi>t</mi></mrow></msup></munderover><mo>&lt;</mo><msub><mi>F</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>H</mi><mrow><mi>s</mi><mi>u</mi><mi>b</mi></mrow></msub><mo>&gt;</mo></mrow>]]></math><img file="FDA00011214977100000211.GIF" wi="422" he="142" /></maths>式中,<img file="FDA00011214977100000212.GIF" wi="70" he="71" />为内积算子;B8、重复B1~B7,计算出遥感图像P中所有检测窗口的总相似度Score,得到总相似度Score分布图S;C、当总相似度Score分布图S中有大于设定阈值的总相似度Score,判定检测区域中存在飞机,其对应检测窗口在遥感图像中的位置即为飞机所在位置。
地址 610031 四川省成都市二环路北一段111号