发明名称 基于加权信息熵和马尔可夫随机场的红外目标分割方法
摘要 本发明涉及一种基于加权信息熵和马尔可夫随机场的红外目标分割方法,技术特征在于:首先利用加权信息熵得到感兴趣区域,继而利用高斯马尔可夫随机场在感兴趣区域内建模来捕捉图像的上下文信息,并在迭代中采用ICM优化算法,使得解空间大大减小、噪声干扰大大降低,从而迭代步数减少、分割精度有明显的改善、分割速度大大提高。许多红外目标分割实验,都验证了本发明的有效性。
申请公布号 CN101847259B 申请公布日期 2012.05.16
申请号 CN201010013647.X 申请日期 2010.01.21
申请人 西北工业大学 发明人 李映;毛星锦;张艳宁
分类号 G06T7/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于加权信息熵和马尔可夫随机场的红外目标分割方法,其特征在于步骤如下:步骤1基于加权信息熵的感兴趣区域提取:步骤a:以大小为m×n的图像块来平分M×N图像,m=2<sup>k</sup>,n=2<sup>l</sup>,k,l=1,2,3,...,得到平分后的若干个子图像;步骤b:对每幅子图像F(u,v),其中0≤u≤M/m-1,0≤v≤N/n-1,利用<maths num="0001"><![CDATA[<math><mrow><mi>H</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>S</mi><mo>=</mo><mn>0</mn></mrow><mn>255</mn></munderover><msup><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mover><mi>s</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>*</mo><msub><mi>p</mi><mi>s</mi></msub><mo>*</mo><mi>log</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>when p<sub>s</sub>=0,let p<sub>s</sub>*log(p<sub>s</sub>)=0计算每幅子图像的加权信息熵WIE值H(u,v),得到M/(m-1)×N/(n-1)大小的熵图像,并计算出熵图像的均值μ和方差δ;其中s表示灰度值,s∈[0,255],p<sub>s</sub>表示每种灰度级对应的概率,<img file="FSB00000586946400012.GIF" wi="24" he="53" />表示红外图像的灰度平均值;步骤c:当最大的熵值的熵图像H(u,v)≥H<sub>T</sub>作为种子,其中H<sub>T</sub>=μ+6×δ;如果H(u,v)<H<sub>T</sub>,将m=m/2,n=n/2,重复步骤a~b,直到熵图像的最大值满足H(u,v)≥H<sub>T</sub>或m=2或n=2;步骤d:在熵图像H(u,v)上,结合种子利用四邻域进行区域增长,得到包含目标的矩形区域;步骤2:在矩形区域内利用马尔可夫随机场建模,然后做随机初始分割,利用ICM进行迭代优化得到目标区域与背景区域的分割图像;所述迭代优化中的权值修正方案为W=800*0.97<sup>t</sup>+1,其中:t表示迭代步数,最优分割问题转化为求解能量最小问题,能量表示为两部分加权和的形式:E=W*E<sub>f</sub>+E<sub>r</sub>,其中<img file="FSB00000586946400013.GIF" wi="642" he="124" />表示能量中的特征成分,<maths num="0002"><![CDATA[<math><mrow><msub><mi>E</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>s</mi></munder><msub><mi>U</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>s</mi></munder><mo>[</mo><mi>&beta;</mi><munder><mi>&Sigma;</mi><mrow><mi>t</mi><mo>&Element;</mo><msub><mi>N</mi><mi>s</mi></msub></mrow></munder><mi>&delta;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>s</mi></msub><mo>,</mo><msub><mi>x</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>表示能量中的区域部分,其中μ<sub>m</sub>,σ<sub>m</sub>分别代表第m类的均值和标准差;如果x<sub>s</sub>=x<sub>t</sub>,则δ(x<sub>s</sub>,x<sub>t</sub>)=-1;如果x<sub>s</sub>≠x<sub>t</sub>,则δ(x<sub>s</sub>,x<sub>t</sub>)=1;β是事先定义好的常量。
地址 710072 陕西省西安市友谊西路127号