发明名称 一种利用多方向一阶偏导数的红外弱小目标检测方法
摘要 本发明涉及一种利用多方向一阶偏导数的红外弱小目标检测方法,它有三大步骤:首先,利用小面模型构造小范围内图像灰度的二元三次函数,在中心区域求解各方向的一阶偏导数,并设计所需系数的卷积模板,各个运算系数可以直接从卷积后的图像中得到;然后,根据小目标在各方向上的一阶偏导数特性,通过内积最大化原则设计增强卷积模板,在各个方向上对目标进行增强;最后,对个各方向上的结果进行相乘融合,在进一步增强目标的同时尽可能地抑制背景,得到最后结果。本发明可广泛应用于红外图像的弱小目标检测,具有广阔的市场前景与应用价值。
申请公布号 CN106548457A 申请公布日期 2017.03.29
申请号 CN201610898902.0 申请日期 2016.10.14
申请人 北京航空航天大学 发明人 白相志;毕研广
分类号 G06T5/00(2006.01)I;G06T5/20(2006.01)I;G06T7/00(2017.01)I;G06T7/194(2017.01)I 主分类号 G06T5/00(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种利用多方向一阶偏导数的红外弱小目标检测方法,其特征在于:该方法具体步骤如下:步骤一:定义两个集合R={‑2 ‑1 0 1 2}与C={‑2 ‑1 0 1 2},由小面模型可知,在集合R×C范围内的三次曲面函数可以由以下离散正交基构成:<img file="FDA0001130729890000011.GIF" wi="627" he="63" /><img file="FDA0001130729890000012.GIF" wi="758" he="87" />将S作为R×C范围内的对称邻域,则曲面函数表示如下:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>f</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>k</mi><mn>1</mn></msub><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><mi>r</mi><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><mi>c</mi><mo>+</mo><msub><mi>k</mi><mn>4</mn></msub><mrow><mo>(</mo><msup><mi>r</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><mi>r</mi><mi>c</mi><mo>+</mo><msub><mi>k</mi><mn>6</mn></msub><mrow><mo>(</mo><msup><mi>c</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>7</mn></msub><mrow><mo>(</mo><msup><mi>r</mi><mn>3</mn></msup><mo>-</mo><mfrac><mn>17</mn><mn>5</mn></mfrac><mi>r</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>8</mn></msub><mrow><mo>(</mo><msup><mi>r</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>c</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msub><mi>k</mi><mn>9</mn></msub><mi>r</mi><mrow><mo>(</mo><msup><mi>c</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>10</mn></msub><mo>(</mo><mrow><msup><mi>c</mi><mn>3</mn></msup><mo>-</mo><mfrac><mn>17</mn><mn>5</mn></mfrac><mi>c</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001130729890000013.GIF" wi="1934" he="279" /></maths>当前中心点(0,0)的水平、垂直及各方向α上的一阶导数为:<maths num="0002"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><mi>f</mi></mrow><mrow><mo>&part;</mo><mi>r</mi></mrow></mfrac><mo>=</mo><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mfrac><mn>17</mn><mn>5</mn></mfrac><msub><mi>k</mi><mn>7</mn></msub><mo>-</mo><mn>2</mn><msub><mi>k</mi><mn>9</mn></msub></mrow>]]></math><img file="FDA0001130729890000014.GIF" wi="486" he="119" /></maths><maths num="0003"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><mi>f</mi></mrow><mrow><mo>&part;</mo><mi>c</mi></mrow></mfrac><mo>=</mo><msub><mi>k</mi><mn>3</mn></msub><mo>-</mo><mfrac><mn>17</mn><mn>5</mn></mfrac><msub><mi>k</mi><mn>10</mn></msub><mo>-</mo><mn>2</mn><msub><mi>k</mi><mn>8</mn></msub></mrow>]]></math><img file="FDA0001130729890000015.GIF" wi="510" he="119" /></maths><maths num="0004"><math><![CDATA[<mrow><msubsup><mi>f</mi><mi>&alpha;</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>f</mi></mrow><mrow><mo>&part;</mo><mi>r</mi></mrow></mfrac><mi>s</mi><mi>i</mi><mi>n</mi><mi>&alpha;</mi><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>f</mi></mrow><mrow><mo>&part;</mo><mi>c</mi></mrow></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mi>&alpha;</mi></mrow>]]></math><img file="FDA0001130729890000016.GIF" wi="510" he="119" /></maths>对于k<sub>2</sub>,k<sub>3</sub>,k<sub>7</sub>,k<sub>8</sub>,k<sub>9</sub>以及k<sub>10</sub>,由最小二乘法设计它们对应的卷积模板为w<sub>2</sub>,w<sub>3</sub>,w<sub>7</sub>,w<sub>8</sub>,w<sub>9</sub>与w<sub>10</sub>,各多项式系数可直接由卷积操作得到;步骤二:根据小目标在各方向上的一阶偏导数特性,通过内积最大化原则设计增强卷积模板,在各个方向上对目标进行增强:将原图像在各个方向上求解一阶导数之后,分别设计4个方向上的滤波器;由于目标一般所占面积较小,实际分布的有效体积受到限制,因此考虑在7×7范围内设计增强卷积模板;卷积的过程可以看作是两个高维向量的内积,将各个方向的目标区域看作是一维列向量V<sub>t</sub>,卷积增强模板看作是一维列向量V<sub>c</sub>,由其方向上的单位向量<img file="FDA0001130729890000017.GIF" wi="27" he="46" />与向量模k相乘,需要设计的量为<img file="FDA0001130729890000018.GIF" wi="43" he="51" />由夹角公式:<maths num="0005"><math><![CDATA[<mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><mo>&lt;</mo><msub><mi>V</mi><mi>t</mi></msub><mo>,</mo><msub><mi>V</mi><mi>c</mi></msub><mo>&gt;</mo></mrow><mrow><mo>|</mo><msub><mi>V</mi><mi>t</mi></msub><mo>|</mo><mo>|</mo><msub><mi>V</mi><mi>c</mi></msub><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA0001130729890000021.GIF" wi="382" he="127" /></maths>得到:<maths num="0006"><math><![CDATA[<mrow><mo>&lt;</mo><msub><mi>V</mi><mi>t</mi></msub><mo>,</mo><mi>k</mi><mover><mi>e</mi><mo>&RightArrow;</mo></mover><mo>&gt;</mo><mo>=</mo><mo>|</mo><msub><mi>V</mi><mi>t</mi></msub><mo>|</mo><mi>k</mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow>]]></math><img file="FDA0001130729890000022.GIF" wi="478" he="63" /></maths>令k=1,当V<sub>t</sub>固定时,为了使目标内积最大,等价于求cosθ的最大值,所以应将增强卷积模板<img file="FDA0001130729890000023.GIF" wi="29" he="46" />设计为与目标的各方向一阶导数相似的形状,实际中的k可以根据实际情况选取;步骤三:对个各方向上的结果进行相乘融合,在进一步增强目标的同时尽可能地抑制背景,得到最后结果:设在前两个步骤中已得到各方向上的增强图,最后的结果图f<sub>result</sub>可表示为:f<sub>result</sub>=f<sub>0</sub>·f<sub>45</sub>·f<sub>90</sub>·f<sub>135</sub>。
地址 100191 北京市海淀区学院路37号