发明名称 面向小目标检测的快速核密度估计图像阈值分割方法
摘要 本发明涉及一种面向小目标检测的快速核密度估计图像阈值分割方法,该方法包括:读入图像,在计算机中得到该图像的灰度矩阵,并设置参数<i>gate</i>值;将图像中相同灰度级的像素点作为一个集合,如果该集合中的像素点数目<i>N<sub>i</sub></i>&gt;<i>gate</i>,则使用FRSDE进行压缩;否则使用RSDE进行压缩;构建关系矩阵<i>M</i>来表示图像上不同灰度级之间的相互关系。将目标函数的求极值问题转化为基于矩阵区域的元素最小化求和问题,从而得到最佳阈值。本发明与已有技术相比,具有过程简单、易于实现、具有鲁棒性好与求解效率高等优点。因此,本发明为高分辨率图像的小目标检测问题提供了一种可行的方案;同时,为复杂背景下小目标图像的检测提供了一种行之有效的技术基础。
申请公布号 CN102903112B 申请公布日期 2015.01.07
申请号 CN201210376033.7 申请日期 2012.09.29
申请人 江南大学 发明人 王骏;王士同;邓赵红;钱鹏江;应文豪;蒋亦樟;倪彤光
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 无锡市大为专利商标事务所(普通合伙) 32104 代理人 殷红梅
主权项 一种面向小目标检测的快速核密度估计图像阈值分割方法,其特征是,该方法包括如下步骤:步骤一:读入图像,在计算机中得到该图像的灰度矩阵,并设置参数gate值;步骤二:对于每个灰度级中的像素点,如果该灰度集像素点数目N<sub>i</sub>&gt;gate,则使用FRSDE进行压缩,并得到压缩后各像素点对应的权值α<sub>i</sub>;对于每个灰度级中的像素点,如果该灰度集像素点数目N<sub>i</sub>≤gate,则使用RSDE进行压缩,并得到压缩后各像素点对应的权值α<sub>i</sub>;步骤三:构建关系矩阵M;该关系矩阵M的每个元素M(i,l)如式(1)进行计算:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msubsup><mi>N</mi><mi>i</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></msubsup></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msubsup><mi>N</mi><mi>l</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></msubsup></munderover><msubsup><mi>N</mi><mi>i</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></msubsup><msubsup><mi>N</mi><mi>l</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></msubsup><msub><mi>&alpha;</mi><mi>j</mi></msub><msub><mi>&alpha;</mi><mi>k</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msubsup><mi>&sigma;</mi><mi>j</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581180980000011.GIF" wi="1562" he="164" /></maths>其中N<sub>i</sub><sup>(r)</sup>为集合V<sub>i</sub>压缩后的像素点数目,V<sub>i</sub>表示灰度级为i的像素点集合,α<sub>j</sub>和α<sub>k</sub>是系数,是步骤二得到的权值,满足<img file="FDA0000581180980000012.GIF" wi="235" he="164" />N<sub>l</sub><sup>(r)</sup>和α<sub>k</sub>的关系同理;K(x<sub>j</sub>‑x<sub>k</sub>,σ<sub>j</sub><sup>2</sup>+σ<sub>k</sub><sup>2</sup>)是高斯核函数;步骤四:将目标函数的求极值问题转化为基于矩阵区域的元素最小化求和问题,从而得到最佳阈值;具体为:对于每个T值,将关系矩阵M的上三角区域划分出区域A和区域B两个顶点相接的直角三角形,随T值依次变化区域A和区域B的边界,并计算区域A和区域B中的所有数据元素总和,如式(2):f(T)=sum(A)+sum(B)                 (2)优化的阈值t<sub>opt</sub>通过式(3)得到:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mtext></mtext><msub><mi>t</mi><mi>opt</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>T</mi></munder><mi>f</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581180980000013.GIF" wi="1195" he="96" /></maths>该式表示t<sub>opt</sub>取函数f(T)为最小时的T值。
地址 214122 江苏省无锡市滨湖区蠡湖大道1800号